DQMHistAnalysisCDCEpics.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/database/DBArray.h>
#include <framework/database/DBObjPtr.h>
#include <dqm/core/DQMHistAnalysis.h>
 
#include <TROOT.h>
#include <TLine.h>
#include <TH2Poly.h>
#include <TH2F.h>
#include <TH1F.h>
#include <TStyle.h>
#include <TLine.h>
#include <numeric>
#include <iostream>
 
namespace Belle2 {
  class DQMHistAnalysisCDCEpicsModule final : public DQMHistAnalysisModule {
 
  public:
 
    DQMHistAnalysisCDCEpicsModule();
 
    ~DQMHistAnalysisCDCEpicsModule();
 
    void initialize() override final;
 
    void beginRun() override final;
 
    void event() override final;
 
    void endRun() override final;
 
    void terminate() override final;
 
    void getHistStyle(TH1F*& htemp, std::string label, double max) const
    {
      gStyle->SetOptStat("ne");
      if (strcmp(label.data(), "adc") == 0)htemp->GetYaxis()->SetRangeUser(max * 0.25, max * 2.25);
      else if (strcmp(label.data(), "tdc") == 0)htemp->GetYaxis()->SetRangeUser(max * 0.90, max * 1.10);
      htemp->Sumw2(0);
    };
 
    float getHistMedian(TH1D* h) const;
 
    TH2Poly* createEffiTH2Poly(const TString& name, const TString& title) ;
 
    void fillEffiTH2Poly(TH2F* hist, TH2Poly* attached, TH2Poly* expected, TH2Poly* efficiency) ;
 
    void fillEffiTH2(TH2F* hist, TH2F* attached, TH2F* expected, TH2F* efficiency) ;
 
 
  protected:
 
    //Canvas for DQM analysis IR plots
    TCanvas* c_hist_adc = nullptr; 
    TH1F* m_hist_adc = nullptr; 
    TCanvas* c_hist_tdc = nullptr; 
    TH1F* m_hist_tdc = nullptr; 
    TCanvas* c_hist_crphi = nullptr; 
    TH1D* m_hist_crphi = nullptr; 
    TCanvas* c_hist_hitsphi = nullptr; 
    TCanvas* c_hist_effphi = nullptr; 
    TH1D* m_hist_effphi = nullptr; 
    TCanvas* c_hist_skimphi[8] = {nullptr}; 
    TH1D* m_hist_skimphi[8] = {nullptr}; 
    TCanvas* c_hist_attach_eff[4] = {nullptr}; 
    TH2F* m_hist_attach_eff[3] = {nullptr}; 
    TH2Poly* m_hist_attach_eff_Poly[3] = {nullptr}; 
    TH1F* m_hist_wire_attach_eff_1d = nullptr; 
    TLine* m_line_ladc  = nullptr; 
    TLine* m_line_hadc  = nullptr; 
    TLine* m_line_ltdc  = nullptr; 
    TLine* m_line_htdc  = nullptr; 
    std::string m_histoDir = ""; 
    std::string m_histoADC = ""; 
    std::string m_histoTDC = ""; 
    std::string m_histoPhiIndex = ""; 
    std::string m_histoPhiEff = ""; 
    std::string m_histoHitsPhi = ""; 
    std::string m_histoTrackingWireEff = ""; 
    bool        m_doTH2PolyTrackingWireEff = false; 
    std::string m_pvPrefix = ""; 
    std::string m_refDir = ""; 
    std::string m_refNamePhi = ""; 
    double m_firstEffBoundary = 0.08; 
    double m_secondEffBoundary = 0.72; 
    MonitoringObject* m_monObj = nullptr; 
    TFile* m_fileRefPhi = nullptr; 
    TH2F* m_histref_phiindex = nullptr; 
    TH1D* m_hist_refphi = nullptr; 
    int m_minevt;
    double m_minadc;
    double m_maxadc;
    double m_mintdc;
    double m_maxtdc;
    double m_phistop;
    double m_phialarm;
    double m_phiwarn;
    TH1D* m_hADCs[300]; 
    TH1D* m_hTDCs[300]; 
  };
} // Belle2 namespace