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 <TEllipse.h>
#include <TLatex.h>
#include <TPad.h>
#include <numeric>
#include <iostream>
#include <cdc/geometry/CDCGeometryPar.h>
#include <cdc/geometry/CDCGeometryParConstants.h>
 
namespace Belle2 {
 

  constexpr unsigned kNumLayers = c_maxNFieldLayers; 
  constexpr unsigned kNumBoards = c_nBoards;
  constexpr std::array<int, 9> slindex = { 8, 14, 20, 26, 32, 38, 44, 50, 56 }; 

  class DQMHistAnalysisCDCEpicsModule final : public DQMHistAnalysisModule {
 
  public:

    DQMHistAnalysisCDCEpicsModule();

    void initialize() override final;

    void beginRun() override final;

    void event() override final;

    void endRun() override final;

    void terminate() override final;

    static void getHistStyle(TH1F*& htemp, const std::string& label, double max)
    {
      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);
    };

    static float getHistMedian(TH1* h);

    static TH2Poly* createEffiTH2Poly(const TString& name, const TString& title) ;

    static 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* m_canv_md_ladc = nullptr; 
    TH1F* m_histmd_ladc = nullptr; 
 
    TCanvas* m_canv_adc = nullptr; 
    TH1F* m_hist_adc = nullptr; 
 
    TCanvas* m_canv_tdc = nullptr; 
    TH1F* m_hist_tdc = nullptr; 
 
    TCanvas* m_canv_crphi = nullptr; 
    TH1D* m_hist_crphi = nullptr; 
 
    TCanvas* m_canv_hitsphi = nullptr; 
 
    TCanvas* m_canv_effphi = nullptr; 
    TH1D* m_hist_effphi = nullptr; 
 
    TCanvas* m_canv_skimphi[8] = {nullptr}; 
    TH1D* m_hist_skimphi[8] = {nullptr}; 
 
    TCanvas* m_canv_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; 
    double m_lbinEdges[kNumLayers + 1] = {0.0}; 
 
    TLine* m_line_ladc_sl01  = nullptr; 
    TLine* m_line_hadc_sl01  = nullptr; 
    TLine* m_line_ltdc_sl01  = nullptr; 
    TLine* m_line_htdc_sl01  = nullptr; 
 
    TLine* m_line_ladc_sl28  = nullptr; 
    TLine* m_line_hadc_sl28  = nullptr; 
    TLine* m_line_ltdc_sl28  = nullptr; 
    TLine* m_line_htdc_sl28  = nullptr; 
 
    std::string m_name_dir = ""; 
    std::string m_name_refdir = ""; 
    std::string m_name_pvpfx = ""; 
    std::string m_fname_refphi = ""; 
    std::string m_hname_ladc = ""; 
    std::string m_hname_badc = ""; 
    std::string m_hname_btdc = ""; 
    std::string m_hname_idxphi = ""; 
    std::string m_hname_effphi = ""; 
    std::string m_hname_hitsphi = ""; 
    std::string m_histoTrackingWireEff = ""; 
    bool  m_doTH2PolyTrackingWireEff = false; 
    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_minphibinsfrac;
 
    double m_minadc_sl01 = 20.0;  
    double m_maxadc_sl01 = 40.0;  
    double m_mintdc_sl01 = 4700.0;
    double m_maxtdc_sl01 = 5200.0;
 
    double m_minadc_sl28 = 40.0;  
    double m_maxadc_sl28 = 180.0; 
    double m_mintdc_sl28 = 4600.0;
    double m_maxtdc_sl28 = 5000.0;
 
    double m_phiwarn = 0.05;
    double m_phialarm = 0.15;
    std::vector<TLine*> m_lines;
 
    TH1D* m_hists_lADC[kNumLayers]; 
    TH1D* m_hists_bADC[kNumBoards]; 
    TH1D* m_hists_bTDC[kNumBoards]; 
  };

} // Belle2 namespace