SVDLocalCalibrationsCheckModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Laura Zani, Giulia Casarosa                              *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <framework/core/Module.h>
#include <vxd/dataobjects/VxdID.h>
 
#include <svd/dataobjects/SVDHistograms.h>
#include <svd/dataobjects/SVDAPVHistograms.h>
#include <svd/dataobjects/SVDSummaryPlots.h>
 
#include <string>
#include <TFile.h>
#include <TTree.h>
#include <TBranch.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TList.h>
#include <TLegend.h>
 
namespace Belle2 {
  class SVDLocalCalibrationsCheckModule : public Module {
 
  public:
 
    SVDLocalCalibrationsCheckModule();
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    /* ROOT file related parameters */
    TFile* m_rootFilePtrREF = nullptr; 
    TFile* m_rootFilePtrCHECK = nullptr; 
    TTree* m_treeREF = nullptr; 
    TTree* m_treeCHECK = nullptr; 
    //branches REF
    TBranch* b_runREF = nullptr; 
    TBranch* b_ladderREF = nullptr; 
    TBranch* b_layerREF = nullptr; 
    TBranch* b_sensorREF = nullptr; 
    TBranch* b_sideREF = nullptr; 
    TBranch* b_stripREF = nullptr; 
    TBranch* b_maskREF = nullptr; 
    TBranch* b_pedestalREF = nullptr; 
    TBranch* b_gainREF = nullptr; 
    TBranch* b_noiseREF = nullptr; 
    TBranch* b_noiseElREF = nullptr; 
    TBranch* b_calPeakADCREF = nullptr; 
    TBranch* b_calPeakTimeREF = nullptr; 
    TBranch* b_pulseWidthREF = nullptr; 
    //branches CHECK
    TBranch* b_runCHECK = nullptr; 
    TBranch* b_ladderCHECK = nullptr; 
    TBranch* b_layerCHECK = nullptr; 
    TBranch* b_sensorCHECK = nullptr; 
    TBranch* b_sideCHECK = nullptr; 
    TBranch* b_stripCHECK = nullptr; 
    TBranch* b_maskCHECK = nullptr; 
    TBranch* b_pedestalCHECK = nullptr; 
    TBranch* b_gainCHECK = nullptr; 
    TBranch* b_noiseCHECK = nullptr; 
    TBranch* b_noiseElCHECK = nullptr; 
    TBranch* b_calPeakTimeCHECK = nullptr; 
    TBranch* b_calPeakADCCHECK = nullptr; 
    TBranch* b_pulseWidthCHECK = nullptr; 
    //branch variables
    UInt_t m_runREF = -1; 
    UInt_t m_layerREF = -1; 
    UInt_t m_ladderREF = -1; 
    UInt_t m_sensorREF = -1; 
    UInt_t m_sideREF = -1; 
    UInt_t m_stripREF = -1; 
    float m_maskREF = -1; 
    float m_noiseREF = -1; 
    float m_noiseElREF = -1; 
    float m_pedestalREF = -1; 
    float m_gainREF = -1; 
    float m_calPeakTimeREF = -1; 
    float m_calPeakADCREF = -1; 
    float m_pulseWidthREF = -1; 
    //branch variables
    UInt_t m_runCHECK = -1; 
    UInt_t m_layerCHECK = -1; 
    UInt_t m_ladderCHECK = -1; 
    UInt_t m_sensorCHECK = -1; 
    UInt_t m_sideCHECK = -1; 
    UInt_t m_stripCHECK = -1; 
    float m_maskCHECK = -1; 
    float m_noiseCHECK = -1; 
    float m_noiseElCHECK = -1; 
    float m_pedestalCHECK = -1; 
    float m_gainCHECK = -1; 
    float m_calPeakTimeCHECK = -1; 
    float m_calPeakADCCHECK = -1; 
    float m_pulseWidthCHECK = -1; 
    std::string m_rootFileNameREF = "SVDLocalCalibrationMonitor_experiment5_run92.root";   
    std::string m_rootFileNameCHECK = "SVDLocalCalibrationMonitor_experiment5_run408.root";   
    std::string m_idFileNameREF = "refID";   
    std::string m_idFileNameCHECK = "checkID";   
    std::string m_outputPdfName = "SVDLocalCalibrationCheck.pdf"; 
    bool m_plotGoodAPVs = false; 
    //analsyis parameters
    int m_cutN_out = -1; 
    float m_cutNoise_ave = -1; 
    float m_cutNoise_out = -1; 
    float m_cutCalpeakADC_ave = -1; 
    float m_cutCalpeakADC_out = -1; 
    float m_cutPedestal_ave = -1; 
    float m_cutPedestal_out = -1; 
  private:
 
    void printConfiguration(); 
    const int m_apvColors[6] = { 1, 2, 8 , kBlue, 6, 28}; 
    void   setAPVHistoStyles(SVDAPVHistograms<TH1F>* m_APVhistos); 
    void createLegends(); 
    TLegend* m_leg2D = nullptr; 
    TLegend* m_legU = nullptr;  
    TLegend* m_legV = nullptr;  
    void printFirstPage();  
    void printLayerPage(int layer);  
    void printPage(VxdID theVxdID, TList* listUBAD, TList* listVBAD, TList* listUGOOD, TList* listVGOOD, TString variable,
                   bool isL3);  
    void printSummaryPages();  
    void printLastPage();  
    int hasAnyProblem(TH1F* h, float cutAve, float cutCOUNT);  
    SVDHistograms<TH2F>* m_h2NoiseREF = nullptr; 
    SVDHistograms<TH2F>* m_h2NoiseCHECK = nullptr; 
    SVDAPVHistograms<TH1F>* m_hNoiseDIFF = nullptr; 
    SVDSummaryPlots* m_hNoiseSummary = nullptr; 
    SVDHistograms<TH2F>* m_h2CalpeakADCREF = nullptr; 
    SVDHistograms<TH2F>* m_h2CalpeakADCCHECK = nullptr; 
    SVDAPVHistograms<TH1F>* m_hCalpeakADCDIFF = nullptr; 
    SVDSummaryPlots* m_hCalpeakADCSummary = nullptr; 
    SVDHistograms<TH2F>* m_h2CalpeakTimeREF = nullptr; 
    SVDHistograms<TH2F>* m_h2CalpeakTimeCHECK = nullptr; 
    SVDAPVHistograms<TH1F>* m_hCalpeakTimeDIFF = nullptr; 
    SVDSummaryPlots* m_hCalpeakTimeSummary = nullptr; 
    SVDHistograms<TH2F>* m_h2PedestalREF = nullptr; 
    SVDHistograms<TH2F>* m_h2PedestalCHECK = nullptr; 
    SVDAPVHistograms<TH1F>* m_hPedestalDIFF = nullptr; 
    SVDSummaryPlots* m_hPedestalSummary = nullptr; 
  };
}