OverlapResidualsModule.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 <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 svdDeltaRes_U = 0,          
          svdClCharge_U_int = 0,      
          svdClSNR_U_int = 0,         
          svdClTime_U_int = 0,        
          svdRes_U_int = 0,           
          svdClIntStrPos_U_int = 0,   
          svdClPos_U_int = 0,         
          svdClPosErr_U_int = 0,      
          svdTruePos_U_int = -99,     
          svdClPhi_U_int = 0,         
          svdClZ_U_int = 0,           
          svdTrkd0 = 0,         
          svdTrkz0 = 0,         
          svdTrkpT = 0,         
          svdTrkpCM = 0,        
          svdTrkTraversedLength_U_int = 0,   
          svdTrkPos_U_int = 0,        
          svdTrkPosOS_U_int = 0,      
          svdTrkPosErr_U_int = 0,     
          svdTrkPosErrOS_U_int = 0,   
          svdTrkQoP_U_int = 0,        
          svdTrkPrime_U_int = 0,      
          svdTrkPrimeOS_U_int = 0,    
          svdTrkPosUnbiased_U_int = 0,   
          svdTrkPosErrUnbiased_U_int = 0,   
          svdTrkQoPUnbiased_U_int = 0,   
          svdTrkPrimeUnbiased_U_int = 0,   
          svdClCharge_U_ext = 0,      
          svdClSNR_U_ext = 0,         
          svdClTime_U_ext = 0,        
          svdRes_U_ext = 0,           
          svdClIntStrPos_U_ext = 0,   
          svdClPos_U_ext = 0,         
          svdClPosErr_U_ext = 0,      
          svdTruePos_U_ext = -99,     
          svdClPhi_U_ext = 0,         
          svdClZ_U_ext = 0,           
          svdTrkTraversedLength_U_ext = 0,   
          svdTrkPos_U_ext = 0,        
          svdTrkPosOS_U_ext = 0,      
          svdTrkPosErr_U_ext = 0,     
          svdTrkPosErrOS_U_ext = 0,   
          svdTrkQoP_U_ext = 0,        
          svdTrkPrime_U_ext = 0,      
          svdTrkPrimeOS_U_ext = 0,    
          svdTrkPosUnbiased_U_ext = 0,   
          svdTrkPosErrUnbiased_U_ext = 0,   
          svdTrkQoPUnbiased_U_ext = 0,   
          svdTrkPrimeUnbiased_U_ext = 0;   
    int svdTrkPXDHits = 0,   
        svdTrkSVDHits = 0,   
        svdTrkCDCHits = 0;   
    unsigned int svdLayer_U_int = 0,    
                 svdLadder_U_int = 0,   
                 svdSensor_U_int = 0,   
                 svdSize_U_int = 0,     
                 svdLayer_U_ext = 0,    
                 svdLadder_U_ext = 0,   
                 svdSensor_U_ext = 0,   
                 svdSize_U_ext = 0;     
    std::vector<float> svdStripCharge_U_int,     
        svdStrip6Samples_U_int,   
        svdStripTime_U_int,       
        svdStripPosition_U_int,   
        svdStripCharge_U_ext,     
        svdStrip6Samples_U_ext,   
        svdStripTime_U_ext,       
        svdStripPosition_U_ext;   
    /* Branches of SVD v-clusters tree */
    float svdDeltaRes_V = 0,          
          svdClCharge_V_int = 0,      
          svdClSNR_V_int = 0,         
          svdClTime_V_int = 0,        
          svdRes_V_int = 0,           
          svdClIntStrPos_V_int = 0,   
          svdClPos_V_int = 0,         
          svdClPosErr_V_int = 0,      
          svdTruePos_V_int = -99,     
          svdClPhi_V_int = 0,         
          svdClZ_V_int = 0,           
          svdTrkTraversedLength_V_int = 0,   
          svdTrkPos_V_int = 0,        
          svdTrkPosOS_V_int = 0,      
          svdTrkPosErr_V_int = 0,     
          svdTrkPosErrOS_V_int = 0,   
          svdTrkQoP_V_int = 0,        
          svdTrkPrime_V_int = 0,      
          svdTrkPrimeOS_V_int = 0,    
          svdTrkPosUnbiased_V_int = 0,   
          svdTrkPosErrUnbiased_V_int = 0,   
          svdTrkQoPUnbiased_V_int = 0,   
          svdTrkPrimeUnbiased_V_int = 0,   
          svdClCharge_V_ext = 0,      
          svdClSNR_V_ext = 0,         
          svdClTime_V_ext = 0,        
          svdRes_V_ext = 0,           
          svdClIntStrPos_V_ext = 0,   
          svdClPos_V_ext = 0,         
          svdClPosErr_V_ext = 0,      
          svdTruePos_V_ext = -99,     
          svdClPhi_V_ext = 0,         
          svdClZ_V_ext = 0,           
          svdTrkTraversedLength_V_ext = 0,   
          svdTrkPos_V_ext = 0,        
          svdTrkPosOS_V_ext = 0,      
          svdTrkPosErr_V_ext = 0,     
          svdTrkPosErrOS_V_ext = 0,   
          svdTrkQoP_V_ext = 0,        
          svdTrkPrime_V_ext = 0,      
          svdTrkPrimeOS_V_ext = 0,    
          svdTrkPosUnbiased_V_ext = 0,   
          svdTrkPosErrUnbiased_V_ext = 0,   
          svdTrkQoPUnbiased_V_ext = 0,   
          svdTrkPrimeUnbiased_V_ext = 0;   
    unsigned int svdLayer_V_int = 0,    
                 svdLadder_V_int = 0,   
                 svdSensor_V_int = 0,   
                 svdSize_V_int = 0,     
                 svdLayer_V_ext = 0,    
                 svdLadder_V_ext = 0,   
                 svdSensor_V_ext = 0,   
                 svdSize_V_ext = 0;     
    std::vector<float> svdStripCharge_V_int,     
        svdStrip6Samples_V_int,   
        svdStripTime_V_int,       
        svdStripPosition_V_int,   
        svdStripCharge_V_ext,     
        svdStrip6Samples_V_ext,   
        svdStripTime_V_ext,       
        svdStripPosition_V_ext;   
  };
}