eclDQM.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
 
/* ECL headers. */
#include <ecl/dbobjects/ECLCrystalCalib.h>
#include <ecl/geometry/ECLGeometryPar.h>
#include <ecl/mapper/ECLChannelMapper.h>
 
/* Basf2 headers. */
#include <mdst/dataobjects/SoftwareTriggerResult.h>
#include <framework/core/HistoModule.h>
#include <framework/core/ModuleParam.templateDetails.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/database/DBObjPtr.h>
 
class TH1F;
class TH2F;
class TProfile;
 
namespace Belle2 {
 
  class EventMetaData;
  class ECLDigit;
  class ECLDsp;
  class ECLTrig;
  class ECLCalDigit;
  class TRGSummary;

  class ECLDQMModule : public HistoModule {  
 
  public:

    ECLDQMModule();

    virtual ~ECLDQMModule();

    virtual void initialize() override;
    virtual void beginRun() override;
    virtual void event() override;
    virtual void endRun() override;
    virtual void terminate() override;

    virtual void defineHisto() override;
 
  private:
    bool isRandomTrigger();

    bool fillInvMassHistogram();

    ECL::ECLGeometryPar* m_geom{nullptr};
    StoreObjPtr<EventMetaData> m_eventmetadata;
    StoreObjPtr<TRGSummary> m_l1Trigger;
    ECL::ECLChannelMapper mapper;
    StoreArray<ECLDigit> m_ECLDigits;
    StoreArray<ECLDsp> m_ECLDsps;
    StoreArray<ECLTrig> m_ECLTrigs;
    StoreArray<ECLCalDigit> m_ECLCalDigits;
    DBObjPtr<ECLCrystalCalib> m_calibrationThrApsd;

    int m_iEvent{ -1};
    std::string m_histogramDirectoryName;
    double m_EnergyUpperThr;
    bool m_DPHYTTYP{0};
    std::string m_pi0PListName;

    std::vector<double> m_HitThresholds = {};
    std::vector<double> m_TotalEnergyThresholds = {};
    std::vector<double> m_TimingThresholds = {};
    std::vector<double> m_HitNumberUpperLimits = {};
    std::vector<std::string> m_WaveformOption;
    std::vector<double> ecltot = {};
    std::vector<double> nhits = {};
    std::vector<int> v_totalthrApsd = {};
 

    int m_DspArray[ECLElementNumbers::c_NCrystals][31] = {};
    double m_PedestalMean[ECLElementNumbers::c_NCrystals] = {};
    double m_PedestalRms[ECLElementNumbers::c_NCrystals] = {};

    TH1F* h_evtot{nullptr};
    TH1F* h_evtot_logic{nullptr};
    TH1F* h_evtot_rand{nullptr};
    TH1F* h_evtot_dphy{nullptr};
    TH1F* h_quality{nullptr};
    TH1F* h_quality_other{nullptr};
    TH1F* h_bad_quality{nullptr};
    TH1F* h_trigtag1{nullptr};
    TH1F* h_adc_hits{nullptr};
    TH1F* h_time_crate_Thr1GeV_large{nullptr};
    TH1F* h_pi0_mass{nullptr};

    std::vector<TH1F*> h_cids = {};
    std::vector<TH1F*> h_edeps = {};
    std::vector<TH1F*> h_time_barrels = {};
    std::vector<TH1F*> h_time_endcaps = {};
    std::vector<TH1F*> h_ncevs = {};
    std::vector<TH1F*> h_cells = {};
    TH1F* h_cell_psd_norm{nullptr};
    std::vector<TH1F*> h_time_crate_Thr1GeV = {};
    TH2F* h_trigtime_trigid{nullptr};
    TH2F* h_trigtag2_trigid{nullptr};
    TProfile* h_pedmean_cellid{nullptr};
    TProfile* h_pedrms_cellid{nullptr};
    TProfile* h_pedrms_thetaid{nullptr};
  };

}; // end Belle2 namespace