DQMHistAnalysisPXDER.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Description : An example module for DQM histogram analysis             *
 * Author: The Belle II Collaboration                                     *
 * Author : Bjoern Spruck, Peter Kodys                                    *
 *                                                                        *
 * Prepared for Belle II geometry                                         *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <dqm/analysis/modules/DQMHistAnalysis.h>
 
#include "TFile.h"
#include "TH1I.h"
#include "TH1F.h"
 
#include <vector>
 
namespace Belle2 {
  class DQMHistAnalysisPXDERModule : public DQMHistAnalysisModule {
 
  public:
 
    DQMHistAnalysisPXDERModule();
    /* Destructor */
    ~DQMHistAnalysisPXDERModule();
 
  private:
 
    void initialize() override final;
    void beginRun() override final;
    void event() override final;
    void endRun() override final;
    void terminate() override final;
 
  private:
 
    //TDirectory* m_oldDir;
 
//     TH1I* m_fHitMapCountsFlag;
//     TH1I* m_fHitMapClCountsFlag;
    TH1I* m_fFiredFlag = nullptr;
    TH1I* m_fClustersFlag = nullptr;
    TH1I* m_fStartRowFlag = nullptr;
    TH1I* m_fChargStartRowFlag = nullptr;
    TH1I* m_fStartRowCountFlag = nullptr;
    TH1I* m_fClusterChargeFlag = nullptr;
    TH1I* m_fPixelSignalFlag = nullptr;
    TH1I* m_fClusterSizeUFlag = nullptr;
    TH1I* m_fClusterSizeVFlag = nullptr;
    TH1I* m_fClusterSizeUVFlag = nullptr;
 
    // Name the histograms, we have to find them anyway every event
//     std::string m_hitMapCounts;
//     std::string m_hitMapClCounts;
 
    std::vector <std::string> m_fired;
    std::vector <std::string> m_clusters;
    std::vector <std::string> m_startRow;
    std::vector <std::string> m_chargStartRow;
    std::vector <std::string> m_startRowCount;
    std::vector <std::string> m_clusterCharge;
    std::vector <std::string> m_pixelSignal;
    std::vector <std::string> m_clusterSizeU;
    std::vector <std::string> m_clusterSizeV;
    std::vector <std::string> m_clusterSizeUV;
 
 
    std::vector <std::string> m_ref_fired;
    std::vector <std::string> m_ref_clusters;
    std::vector <std::string> m_ref_startRow;
    std::vector <std::string> m_ref_chargStartRow;
    std::vector <std::string> m_ref_startRowCount;
    std::vector <std::string> m_ref_clusterCharge;
    std::vector <std::string> m_ref_pixelSignal;
    std::vector <std::string> m_ref_clusterSizeU;
    std::vector <std::string> m_ref_clusterSizeV;
    std::vector <std::string> m_ref_clusterSizeUV;
 
    //int m_nPixels;
    int c_nVXDLayers = 0;
    int c_nPXDLayers = 0;
    int c_nSVDLayers = 0;
    int c_firstVXDLayer = 0;
    int c_lastVXDLayer = 0;
    int c_firstPXDLayer = 0;
    int c_lastPXDLayer = 0;
    int c_firstSVDLayer = 0;
    int c_lastSVDLayer = 0;
    int c_nPXDSensors = 0;
 
    void getIDsFromIndex(const int Index, int& Layer, int& Ladder, int& Sensor)  const;
    int SetFlag(int Type, int bin, double* pars, double ratio, const std::string& name_hist, const std::string& name_refhist,
                TH1I* flaghist);
//     int SetFlag(int Type, int bin, double* pars, double ratio, TH1I* hist, TH1I* refhist, TH1I* flaghist);
 
    std::string m_refFileName;
    TFile* m_refFile = nullptr;
 
    TH1* GetHisto(TString histoname);
 
  };
 
}