SVDOccupancyAnalysisModule.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.                  *
 **************************************************************************/
 
#pragma once
 
#include <framework/core/Module.h>
#include <vxd/dataobjects/VxdID.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <svd/dataobjects/SVDShaperDigit.h>
#include <svd/dataobjects/SVDHistograms.h>
#include <svd/dataobjects/SVDSummaryPlots.h>
#include <framework/dataobjects/EventMetaData.h>
 
#include <mdst/dataobjects/SoftwareTriggerResult.h>
 
#include <svd/calibration/SVDNoiseCalibrations.h>
 
 
#include <string>
#include <TFile.h>
#include <TH1F.h>
#include <TH2F.h>
 
namespace Belle2 {
  class SVDOccupancyAnalysisModule : public Module {
 
  public:
 
    SVDOccupancyAnalysisModule();
 
    virtual ~SVDOccupancyAnalysisModule() {};
 
    virtual void initialize() override; 
    virtual void beginRun() override; 
    virtual void event() override; 
    virtual void endRun() override; 
    /* user-defined parameters */
    std::string m_rootFileName = "";   
    std::string m_ShaperDigitName = "SVDShaperDigits";   
    /* ROOT file related parameters */
    TFile* m_rootFilePtr = nullptr; 
    float m_group = 10000; 
    float m_minZS = 3; 
    float m_maxZS = 6; 
    int m_pointsZS = 7; 
    bool m_FADCmode = true; 
  private:
 
    StoreObjPtr<SoftwareTriggerResult> m_resultStoreObjectPointer;
    bool m_skipRejectedEvents = false; 
    int m_nEvents = 0; 
    StoreArray<SVDShaperDigit> m_svdShapers; 
    StoreObjPtr<EventMetaData> m_eventMetaData; 
    SVDNoiseCalibrations m_NoiseCal; 
    std::vector<VxdID> m_SVDModules;
 
    //layer summary
    float m_distr_Nbins = 10000; 
    float m_distr_min = 0; 
    float m_distr_max = 100; 
    TH1F* m_occ_L3U = nullptr; 
    TH1F* m_occ_L3V = nullptr; 
    TH1F* m_occ_L4U = nullptr; 
    TH1F* m_occ_L4V = nullptr; 
    TH1F* m_occ_L5U = nullptr; 
    TH1F* m_occ_L5V = nullptr; 
    TH1F* m_occ_L6U = nullptr; 
    TH1F* m_occ_L6V = nullptr; 
    //SHAPER
    SVDSummaryPlots* m_hit = nullptr; 
    SVDHistograms<TH1F>* m_histo_dist = nullptr; 
    SVDHistograms<TH1F>* m_histo_occ = nullptr; 
    SVDHistograms<TH1F>* m_histo_zsOcc = nullptr; 
    SVDHistograms<TH1F>* m_histo_zsOccSQ = nullptr; 
    SVDHistograms<TH2F>* m_histo_occtdep = nullptr; 
  };
}