OverlapResidualsModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Filippo Dattola                                          *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef OVERLAPRESIDUALSMODULE_H
#define OVERLAPRESIDUALSMODULE_H
 
#include <pxd/dataobjects/PXDCluster.h>
#include <svd/dataobjects/SVDCluster.h>
 
#include <framework/core/HistoModule.h>
#include <framework/datastore/StoreArray.h>
 
#include <TH1F.h>
#include <TH2F.h>
#include <TTree.h>
 
namespace Belle2 {
  class OverlapResidualsModule : public HistoModule {  // <- derived from HistoModule
 
  public:
    OverlapResidualsModule();
    void initialize() override;
    void event() override;
    void defineHisto() override;
 
  private:
    bool m_ExpertLevel;
    std::string m_recoTracksStoreArrayName{"RecoTracks"};
    StoreArray<PXDCluster> m_pxdcluster;
    StoreArray<SVDCluster> m_svdcluster;
    TH1F* h_U_DeltaRes = nullptr;
    TH1F* h_V_DeltaRes = nullptr;
    TH1F* h_U_DeltaRes_PXD = nullptr;
    TH1F* h_U_DeltaRes_PXD_Lyr1 = nullptr;
    TH1F* h_U_DeltaRes_PXD_Lyr2 = nullptr;
    TH1F* h_V_DeltaRes_PXD = nullptr;
    TH1F* h_V_DeltaRes_PXD_Lyr1 = nullptr;
    TH1F* h_V_DeltaRes_PXD_Lyr2 = nullptr;
    TH1F* h_U_DeltaRes_SVD = nullptr;
    TH1F* h_U_DeltaRes_SVD_Lyr3 = nullptr;
    TH1F* h_U_DeltaRes_SVD_Lyr4 = nullptr;
    TH1F* h_U_DeltaRes_SVD_Lyr5 = nullptr;
    TH1F* h_U_DeltaRes_SVD_Lyr6 = nullptr;
    TH1F* h_V_DeltaRes_SVD = nullptr;
    TH1F* h_V_DeltaRes_SVD_Lyr3 = nullptr;
    TH1F* h_V_DeltaRes_SVD_Lyr4 = nullptr;
    TH1F* h_V_DeltaRes_SVD_Lyr5 = nullptr;
    TH1F* h_V_DeltaRes_SVD_Lyr6 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr6 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr5 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr4 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr3 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr2 = nullptr;
    TH2F* h_DeltaResUPhi_Lyr1 = nullptr;
    TH2F* h_DeltaResUz_Lyr6 = nullptr;
    TH2F* h_DeltaResUz_Lyr5 = nullptr;
    TH2F* h_DeltaResUz_Lyr4 = nullptr;
    TH2F* h_DeltaResUz_Lyr3 = nullptr;
    TH2F* h_DeltaResUz_Lyr2 = nullptr;
    TH2F* h_DeltaResUz_Lyr1 = nullptr;
    TH2F* h_DeltaResVz_Lyr6 = nullptr;
    TH2F* h_DeltaResVz_Lyr5 = nullptr;
    TH2F* h_DeltaResVz_Lyr4 = nullptr;
    TH2F* h_DeltaResVz_Lyr3 = nullptr;
    TH2F* h_DeltaResVz_Lyr2 = nullptr;
    TH2F* h_DeltaResVz_Lyr1 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr6 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr5 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr4 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr3 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr2 = nullptr;
    TH2F* h_DeltaResVPhi_Lyr1 = nullptr;
    TH1F* h_SVDstrips_Mult = nullptr;
    TH1F* h_U_Cl1Cl2_DeltaRes[5] = {nullptr};
    TH1F* h_V_Cl1Cl2_DeltaRes[5] = {nullptr};
    TH2F* h_Lyr6[17][6] = {nullptr}; //[Ladder][sensor number]
    TH2F* h_Lyr5[13][5] = {nullptr}; //[Ladder][sensor number]
    TH2F* h_Lyr4[11][4] = {nullptr}; //[Ladder][sensor number]
    TH2F* h_Lyr3[8][3] = {nullptr}; //[Ladder][sensor number]
    TH2F* h_Lyr2[13][3] = {nullptr}; //[Ladder][sensor number]
    TH2F* h_Lyr1[9][3] = {nullptr}; //[Ladder][sensor number]
    TTree* t_PXD = nullptr;
    TTree* t_SVD_U = nullptr;
    TTree* t_SVD_V = nullptr;
    /* Branches of PXD tree */
    float deltaResU_PXD = 0, 
          intResU_PXD = 0,   
          intResV_PXD = 0,   
          intU_PXD = 0,      
          intV_PXD = 0,      
          intPhi_PXD = 0,    
          intZ_PXD = 0,      
          extResU_PXD = 0,   
          extResV_PXD = 0,   
          extU_PXD = 0,      
          extV_PXD = 0,      
          extPhi_PXD = 0,    
          extZ_PXD = 0;      
    unsigned int intLayer_PXD = 0,  
                 intLadder_PXD = 0, 
                 intSensor_PXD = 0, 
                 extLayer_PXD = 0,  
                 extLadder_PXD = 0, 
                 extSensor_PXD = 0; 
    /* Branches of SVD u-clusters tree */
    float deltaRes_SVD_U = 0, 
          intRes_SVD_U = 0,   
          intClPos_SVD_U = 0,      
          intClPosErr_SVD_U = 0,      
          intTruePos_SVD_U = -99,      
          intClPhi_SVD_U = 0,   
          intClZ_SVD_U = 0,     
          intTrkPos_SVD_U = 0, 
          intTrkPosErr_SVD_U = 0, 
          intTrkQoP_SVD_U = 0, 
          intTrkPrime_SVD_U = 0,     
          extRes_SVD_U = 0,   
          extClPos_SVD_U = 0,      
          extClPosErr_SVD_U = 0,      
          extTruePos_SVD_U = -99,      
          extClPhi_SVD_U = 0,   
          extClZ_SVD_U = 0,     
          extTrkPos_SVD_U = 0, 
          extTrkPosErr_SVD_U = 0, 
          extTrkQoP_SVD_U = 0, 
          extTrkPrime_SVD_U = 0;     
    unsigned int intLayer_SVD_U = 0,  
                 intLadder_SVD_U = 0, 
                 intSensor_SVD_U = 0, 
                 intSize_SVD_U = 0, 
                 extLayer_SVD_U = 0,  
                 extLadder_SVD_U = 0, 
                 extSensor_SVD_U = 0, 
                 extSize_SVD_U = 0; 
    /* Branches of SVD v-clusters tree */
    float deltaRes_SVD_V = 0, 
          intRes_SVD_V = 0,   
          intClPos_SVD_V = 0,      
          intClPosErr_SVD_V = 0,      
          intTruePos_SVD_V = -99,      
          intClPhi_SVD_V = 0,   
          intClZ_SVD_V = 0,     
          intTrkPos_SVD_V = 0, 
          intTrkPosErr_SVD_V = 0, 
          intTrkQoP_SVD_V = 0, 
          intTrkPrime_SVD_V = 0,     
          extRes_SVD_V = 0,   
          extClPos_SVD_V = 0,      
          extClPosErr_SVD_V = 0,      
          extTruePos_SVD_V = -99,      
          extClPhi_SVD_V = 0,   
          extClZ_SVD_V = 0,     
          extTrkPos_SVD_V = 0, 
          extTrkPosErr_SVD_V = 0, 
          extTrkQoP_SVD_V = 0, 
          extTrkPrime_SVD_V = 0;     
    unsigned int intLayer_SVD_V = 0,  
                 intLadder_SVD_V = 0, 
                 intSensor_SVD_V = 0, 
                 intSize_SVD_V = 0, 
                 extLayer_SVD_V = 0,  
                 extLadder_SVD_V = 0, 
                 extSensor_SVD_V = 0, 
                 extSize_SVD_V = 0; 
  };
}
#endif /* OVERLAPRESIDUALS_H */