Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
CaloAna.h
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CaloAna.h
1 #ifndef CALOANA_H__
2 #define CALOANA_H__
3 
4 #include <fun4all/SubsysReco.h>
5 #include <vector>
6 // Forward declarations
7 class Fun4AllHistoManager;
8 class PHCompositeNode;
9 class TFile;
10 class TNtuple;
11 class TTree;
12 class TH2F;
13 class TH1F;
14 class TH1;
15 class TProfile2D;
16 class TLorentzVector;
17 
18 class CaloAna : public SubsysReco
19 {
20  public:
22  CaloAna(const std::string &name = "CaloAna", const std::string &fname = "MyNtuple.root");
23 
25  virtual ~CaloAna();
26 
28  int Init(PHCompositeNode *);
29 
31  int process_event(PHCompositeNode *);
32 
34  int End(PHCompositeNode *);
35 
36  int process_g4hits(PHCompositeNode *);
37  int process_g4cells(PHCompositeNode *);
38  int process_towers(PHCompositeNode *);
39  int process_clusters(PHCompositeNode *);
40 
41  void Detector(const std::string &name) { detector = name; }
42  void set_timing_cut_width(const int &t) { _range = t;}
43  void set_vertex_cut(const float &v) { _vz = v;}
44  void apply_vertex_cut(bool Vtx_cut) { m_vtxCut = Vtx_cut; }
45 
46  std::pair<double, double> fitHistogram(TH1F* h) ;
47  void fitEtaSlices(std::string infile, std::string outfile,std::string cdbFile);
48 
49  protected:
50  std::string detector;
51  std::string outfilename;
52  int Getpeaktime(TH1 *h);
53  Fun4AllHistoManager *hm = nullptr;
54  TFile *outfile = nullptr;
55  TH2F* h_emcal_mbd_correlation = nullptr;
56  TH2F* h_ohcal_mbd_correlation = nullptr;
57  TH2F* h_ihcal_mbd_correlation = nullptr;
58  TH2F* h_emcal_hcal_correlation = nullptr;
59  TH2F* h_emcal_zdc_correlation = nullptr;
60 
61  TH1F* h_InvMass = nullptr;
62  TH1F* h_InvMassMix = nullptr;
63 
64  TH2F* h_cemc_etaphi = nullptr;
65  TH2F* h_hcalin_etaphi = nullptr;
66  TH2F* h_hcalout_etaphi = nullptr;
67  TH2F* h_cemc_etaphi_wQA = nullptr;
68  TH2F* h_hcalin_etaphi_wQA = nullptr;
69  TH2F* h_hcalout_etaphi_wQA = nullptr;
70  TH1* h_totalzdc_e;
71 
72  TProfile2D* h_cemc_etaphi_time = nullptr;
73  TProfile2D* h_hcalin_etaphi_time = nullptr;
74  TProfile2D* h_hcalout_etaphi_time = nullptr;
75 
76  TProfile2D* h_cemc_etaphi_badChi2 = nullptr;
77  TProfile2D* h_hcalin_etaphi_badChi2 = nullptr;
78  TProfile2D* h_hcalout_etaphi_badChi2 = nullptr;
79 
80  TH1* hzdctime;
81  TH1* hmbdtime;
82  TH1* hemcaltime;
83  TH1* hihcaltime;
84  TH1* hohcaltime;
85 
86  TH1* hzdctime_cut;
87  TH1* hmbdtime_cut;
88  TH1* hemcaltime_cut;
89  TH1* hihcaltime_cut;
90  TH1* hohcaltime_cut;
91 
92  TH1* hvtx_z_raw;
93  TH1* hvtx_z_cut;
94 
95  TH1* hzdcSouthraw;
96  TH1* hzdcNorthraw;
97  TH1* hzdcSouthcalib;
98  TH1* hzdcNorthcalib;
99 
100  TH1F* h_clusE;
101  TH2F* h_etaphi_clus;
102 
103  TNtuple *g4hitntuple = nullptr;
104  TNtuple *g4cellntuple = nullptr;
105  TTree *towerntuple = nullptr;
106  TNtuple *clusterntuple = nullptr;
107  std::vector<float> m_energy;
108  std::vector<int> m_etabin;
109  std::vector<int> m_phibin;
110  std::vector<int> m_time;
111 
112  std::vector<float> m_hcalin_energy;
113  std::vector<int> m_hcalin_etabin;
114  std::vector<int> m_hcalin_phibin;
115  std::vector<int> m_hcalin_time;
116 
117  std::vector<float> m_hcalout_energy;
118  std::vector<int> m_hcalout_etabin;
119  std::vector<int> m_hcalout_phibin;
120  std::vector<int> m_hcalout_time;
121 
122 
123  std::vector<float> m_zdc_energy;
124  std::vector<int> m_zdc_index;
125  std::vector<int> m_zdc_side;
126 
127  std::vector<float> m_bbc_energy;
128  std::vector<int> m_bbc_type;
129  std::vector<int> m_bbc_side;
130  int _eventcounter;
131  int _range = 1;
132  float _vz = 0.;
133  bool m_vtxCut = false;
134  bool dynMaskClus = false;
135 
136  TH1F* h_pt1;
137  TH1F* h_pt2;
138  TH1F* h_nclusters;
139  TH1F* h_mass_eta_lt[96];
140  TH1F* h_emcal_e_eta;
141 
142  float target_pi0_mass = 0.145;
143 
144 };
145 
146 #endif