SiCalo/INTTCalo/INTT_Calo_Ana.h

0001 #ifndef INTT_CALO_ANA_H
0002 #define INTT_CALO_ANA_H
0003
0004 #include <string>
0005
0006 string idatafile; // input data file (INTT-Calo nDST)
0007 int NmaxEvent;    // # of events in the data file
0008 string rootfile;  // output ROOT file
0009
0010 // imode==1 : run AnalyzeEvent()
0011 // imode==2 : Event Display
0012 int imode = 2;
0013
0014 //pointer to the TTree for the data.0
0015 TTree *evtTree;
0016 TTree *caloTree;  //EMCal clusters
0017 TTree *topoTree;  //topoClusters (EMC+iHCal+oHCal)
0018 TTree *SiClusTree;
0019 TTree *SiClusAllTree;
0020 TTree *trackTree;
0021 TTree *truthTree; // Simulation only. Truth info
0022
0023 // data in evtTree
0024 Int_t evt;
0025 float xgvtx;  //Global vertex X
0026 float ygvtx;  //Global vertex Y
0027 float zgvtx;  //Global vertex Z
0028
0029 // data in caloTree
0030 vector<float> *energy=nullptr;
0031 vector<float> *emc_phi=nullptr;
0032 vector<float> *emc_eta=nullptr;
0033 vector<float> *emc_r=nullptr;
0034 vector<float> *emc_z  =nullptr;
0035 vector<vector<int>>   *emc_tower_ieta = nullptr;
0036 vector<vector<int>>   *emc_tower_iphi = nullptr;
0037 vector<vector<float>> *emc_tower_e = nullptr;
0038 vector<vector<float>> *emc_tower_x = nullptr;
0039 vector<vector<float>> *emc_tower_y = nullptr;
0040 vector<vector<float>> *emc_tower_z = nullptr;
0041 vector<vector<float>> *emc_tower_r = nullptr;
0042 vector<vector<float>> *emc_tower_eta = nullptr;
0043 vector<vector<float>> *emc_tower_phi = nullptr;
0044 vector<vector<float>> *emc_tower_time = nullptr;
0045
0046
0047 // data in topoTree
0048 vector<float> *top_e  =nullptr;
0049 vector<float> *top_x  =nullptr;
0050 vector<float> *top_y  =nullptr;
0051 vector<float> *top_z  =nullptr;
0052 vector<float> *top_r  =nullptr;
0053 vector<float> *top_phi=nullptr;
0054 vector<float> *top_eta=nullptr;
0055 vector<float> *emc_e  =nullptr;
0056 vector<float> *ihc_e  =nullptr;
0057 vector<float> *ohc_e  =nullptr;
0058
0059 // data in SiClustTree
0060 vector<int> *layer=nullptr;
0061 vector<int> *trkid=nullptr;
0062 vector<float> *si_x=nullptr;
0063 vector<float> *si_y=nullptr;
0064 vector<float> *si_z=nullptr;
0065
0066 // data in SiClustAllTree
0067 vector<int> *a_layer=nullptr;
0068 //vector<int> *a_trkid=nullptr;
0069 vector<float> *a_si_x=nullptr;
0070 vector<float> *a_si_y=nullptr;
0071 vector<float> *a_si_z=nullptr;
0072 vector<float> *a_si_t=nullptr;
0073
0074 // data in trackTree
0075 vector<float> *eta0=nullptr;
0076 vector<float> *phi0=nullptr;
0077 vector<float> *pt0=nullptr;
0078 vector<float> *vx0=nullptr;
0079 vector<float> *vy0=nullptr;
0080 vector<float> *vz0=nullptr;
0081
0082 // SIM truth info
0083 vector<int> *truth_pid = nullptr;
0084 vector<int> *truth_id  = nullptr;
0085 vector<float> *truth_px = nullptr;
0086 vector<float> *truth_py = nullptr;
0087 vector<float> *truth_pz = nullptr;
0088 vector<float> *truth_x = nullptr;
0089 vector<float> *truth_y = nullptr;
0090 vector<float> *truth_z = nullptr;
0091 vector<float> *truth_e = nullptr;
0092 vector<float> *truth_pt = nullptr;
0093 vector<float> *truth_eta = nullptr;
0094 vector<float> *truth_phi = nullptr;
0095 vector<float> *truth_primary_id = nullptr;
0096 vector<float> *truth_parent_id = nullptr;
0097 vector<float> *truth_vtxid = nullptr;
0098 vector<float> *truth_vtx_x = nullptr;
0099 vector<float> *truth_vtx_y = nullptr;
0100 vector<float> *truth_vtx_z = nullptr;
0101
0102
0103
0104 int i_event;   // current event that is analyzed/displayed
0105 // functions
0106
0107 void INTT_Calo_Ana(void);
0108
0109 #endif