SVDHotStripFinderModule.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#ifndef SVDHOTSTRIPFINDERMODULE_H
#define SVDHOTSTRIPFINDERMODULE_H
 
#include <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/dataobjects/EventMetaData.h>
 
#include <svd/dataobjects/SVDShaperDigit.h>
#include <vxd/dataobjects/VxdID.h>
 
#include <svd/dataobjects/SVDHistograms.h>
#include <svd/dataobjects/SVDSummaryPlots.h>
 
#include <string>
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TList.h>
 
namespace Belle2 {
  class SVDHotStripFinderModule : public Module {
  public:
    SVDHotStripFinderModule();
 
    virtual ~SVDHotStripFinderModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    virtual void endRun() override;
 
    virtual void terminate() override;
 
  private:
 
    /* user-defined parameters */
    StoreArray<SVDShaperDigit> m_storeDigits; 
    StoreObjPtr<EventMetaData> m_eventMetaData; 
    std::string m_rootFileName;   
    std::string m_ShaperDigitName; 
    float m_thr; 
    int m_base; 
    bool m_useHSFinderV1 = true; 
    float m_absThr; 
    float m_relOccPrec; 
    bool m_verbose; 
    /* ROOT file related parameters */
    TFile* m_rootFilePtr = nullptr; 
    //dbobject related parameters
    float m_zs; 
    int m_firstExp; 
    int m_firstRun; 
    int m_lastExp; 
    int m_lastRun; 
    static const int m_nLayers = 4; 
    static const int m_nLadders = 16; 
    static const int m_nSensors = 5; 
    static const int m_nSides = 2; 
    TList* m_histoList_occu = nullptr;
    SVDHistograms<TH1F>* hm_occupancy = nullptr; 
    SVDHistograms<TH1F>* hm_hot_strips = nullptr; 
    SVDHistograms<TH1F>* hm_occupancy_after = nullptr; 
    SVDHistograms<TH1F>* hm_occAll = nullptr; 
    SVDHistograms<TH1F>* hm_occHot = nullptr; 
    SVDHistograms<TH1F>* hm_occAfter = nullptr; 
    SVDHistograms<TH1F>* hm_dist = nullptr; 
    SVDHistograms<TH1F>* hm_dist1 = nullptr; 
    SVDHistograms<TH2F>* hm_dist12 = nullptr; 
    SVDSummaryPlots* m_hHotStripsSummary = nullptr; 
    TH1F*  h_tot_dqm = nullptr;   
    TH1F*  h_tot_dqm1 = nullptr;  
    TH1F*  h_tot_dist = nullptr;  
    TH1F*  h_tot_dist1 = nullptr; 
    TH2F*  h_tot_dist12 = nullptr; 
    TH1F*  h_nevents = nullptr;  
    TH1F* createHistogram1D(const char* name, const char* title,
                            Int_t nbins, Double_t min, Double_t max,
                            const char* xtitle, TList* histoList = nullptr);
 
    TH2F* createHistogram2D(const char* name, const char* title,
                            Int_t nbinsX, Double_t minX, Double_t maxX, const char* titleX,
                            Int_t nbinsY, Double_t minY, Double_t maxY, const char* titleY,
                            TList* histoList = nullptr);  
    bool theHSFinder(double* stripOccAfterAbsCut, int* hsflag, int nstrips);
 
 
 
  };
}
#endif