]>
Commit | Line | Data |
---|---|---|
4998becf | 1 | #ifndef ALIANALYSISETCUTS_H |
2 | #define ALIANALYSISETCUTS_H | |
3 | //_________________________________________________________________________ | |
4 | // Utility Class for transverse energy studies | |
5 | // Base class for ESD & MC analysis | |
6 | // - cuts for reconstruction and MonteCarlo | |
7 | // | |
8 | //*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL) | |
9 | //_________________________________________________________________________ | |
2fbf38ac | 10 | |
4998becf | 11 | #include "TNamed.h" |
2fbf38ac | 12 | |
4998becf | 13 | class AliAnalysisEtCuts : public TNamed |
2fbf38ac | 14 | { |
4998becf | 15 | public: |
2fbf38ac | 16 | |
4998becf | 17 | AliAnalysisEtCuts(); |
18 | virtual ~AliAnalysisEtCuts(); | |
2fbf38ac | 19 | |
0fa8c632 | 20 | virtual void SetPbPbDefaults(); |
21 | ||
4998becf | 22 | // Getters |
23 | // Common | |
24 | Double_t GetCommonEtaCut() const { return fCommonEtaCut; } | |
25 | Double_t GetCommonClusterEnergyCut() const { return fCommonClusterEnergyCut; } | |
26 | Double_t GetCommonTrackPtCut() const { return fCommonTrackPtCut; } | |
27 | Int_t GetCommonSingleCell() const { return fCommonSingleCell; } | |
e16aac77 | 28 | |
4998becf | 29 | // GeometryPhos |
30 | Double_t GetGeometryPhosEtaAccCut() const { return fGeometryPhosEtaAccCut; } | |
31 | Double_t GetGeometryPhosPhiAccMinCut() const { return fGeometryPhosPhiAccMinCut; } | |
32 | Double_t GetGeometryPhosPhiAccMaxCut() const { return fGeometryPhosPhiAccMaxCut; } | |
33 | Double_t GetGeometryPhosDetectorRadius() const { return fGeometryPhosDetectorRadius; } | |
34 | // GeometryEmcal | |
35 | Double_t GetGeometryEmcalEtaAccCut() const { return fGeometryEmcalEtaAccCut; } | |
36 | Double_t GetGeometryEmcalPhiAccMinCut() const { return fGeometryEmcalPhiAccMinCut; } | |
37 | Double_t GetGeometryEmcalPhiAccMaxCut() const { return fGeometryEmcalPhiAccMaxCut; } | |
38 | Double_t GetGeometryEmcalDetectorRadius() const { return fGeometryEmcalDetectorRadius; } | |
39 | // Reconstructed | |
40 | Double_t GetReconstructedVertexXCut() const { return fReconstructedVertexXCut; } | |
41 | Double_t GetReconstructedVertexYCut() const { return fReconstructedVertexYCut; } | |
42 | Double_t GetReconstructedVertexZCut() const { return fReconstructedVertexZCut; } | |
43 | Double_t GetReconstructedIPxyCut() const { return fReconstructedIPxyCut; } | |
44 | Double_t GetReconstructedIPzCut() const { return fReconstructedIPzCut; } | |
45 | Int_t GetReconstructedNTpcClustersCut() const { return fReconstructedNTpcClustersCut; } | |
46 | Int_t GetReconstructedNItsClustersCut() const { return fReconstructedNItsClustersCut; } | |
e16aac77 | 47 | Double_t GetReconstructedPidCut() const { return fReconstructedPidCut; } |
4998becf | 48 | // ReconstructedPhos |
ef647350 | 49 | Char_t GetPhosClusterType() const { return fReconstructedPhosClusterType; } |
4998becf | 50 | Double_t GetReconstructedPhosClusterEnergyCut() const { return fReconstructedPhosClusterEnergyCut; } |
51 | Double_t GetReconstructedPhosSingleCellEnergyCut() const { return fReconstructedPhosSingleCellEnergyCut; } | |
393b8149 | 52 | Double_t GetPhosTrackDistanceCut() const { return fPhosTrackDistanceCut; } |
53 | Double_t GetPhosTrackDxCut() const { return fPhosTrackDxCut; } | |
54 | Double_t GetPhosTrackDzCut() const { return fPhosTrackDzCut; } | |
ef647350 | 55 | Double_t GetPhosTrackRCut() const { return fPhosTrackRCut; } |
56 | ||
57 | Double_t GetPhosBadDistanceCut() const { return fPhosBadDistanceCut; } | |
58 | ||
4998becf | 59 | // ReconstructedEmcal |
ef647350 | 60 | Char_t GetEmcalClusterType() const { return fReconstructedEmcalClusterType; } |
4998becf | 61 | Double_t GetReconstructedEmcalClusterEnergyCut() const { return fReconstructedEmcalClusterEnergyCut; } |
62 | Double_t GetReconstructedEmcalSingleCellEnergyCut() const { return fReconstructedEmcalSingleCellEnergyCut; } | |
393b8149 | 63 | Double_t GetEmcalTrackDistanceCut() const { return fEmcalTrackDistanceCut; } |
64 | Double_t GetEmcalTrackDxCut() const { return fEmcalTrackDxCut; } | |
65 | Double_t GetEmcalTrackDzCut() const { return fEmcalTrackDzCut; } | |
66 | ||
4998becf | 67 | // MonteCarlo |
68 | Double_t GetMonteCarloSingleChargedParticle() const { return fMonteCarloSingleChargedParticle; } | |
69 | Double_t GetMonteCarloNeutralParticle() const { return fMonteCarloNeutralParticle; } | |
ce546038 | 70 | // Hist: TTree and histogram info |
71 | Bool_t GetHistMakeTree() const { return fHistMakeTree; } | |
8ea99ab0 | 72 | Bool_t GetHistMakeTreeDeposit() const { return fHistMakeTreeDeposit; } |
0fa8c632 | 73 | // |
74 | Int_t GetHistNbinsMult() const { return fHistNbinsMult; } | |
75 | Double_t GetHistMinMult() const { return fHistMinMult; } | |
d2b935c4 | 76 | Double_t GetHistMaxMult() const { return fHistMaxMult; } |
0fa8c632 | 77 | // |
78 | Int_t GetHistNbinsTotEt() const { return fHistNbinsTotEt; } | |
79 | Double_t GetHistMinTotEt() const { return fHistMinTotEt; } | |
80 | Double_t GetHistMaxTotEt() const { return fHistMaxTotEt; } | |
81 | // | |
82 | Int_t GetHistNbinsParticleEt() const { return fHistNbinsParticleEt; } | |
83 | Double_t GetHistMinParticleEt() const { return fHistMinParticleEt; } | |
84 | Double_t GetHistMaxParticleEt() const { return fHistMaxParticleEt; } | |
85 | // | |
86 | Int_t GetHistNbinsParticlePt() const { return fHistNbinsParticlePt; } | |
d2b935c4 | 87 | Double_t GetHistMinParticlePt() const { return fHistMinParticlePt; } |
0fa8c632 | 88 | Double_t GetHistMaxParticlePt() const { return fHistMaxParticlePt; } |
393b8149 | 89 | |
ef647350 | 90 | |
91 | ||
537e541d | 92 | Short_t GetDetectorPhos() const { return fgkDetectorPhos; } |
93 | Short_t GetDetectorEmcal() const { return fgkDetectorEmcal; } | |
2fbf38ac | 94 | |
4998becf | 95 | // Setters |
96 | // Common | |
97 | void SetCommonEtaCut(const Double_t val) { fCommonEtaCut = val; } | |
98 | void SetCommonClusterEnergyCut(const Double_t val) { fCommonClusterEnergyCut = val; } | |
99 | void SetCommonTrackPtCut(const Double_t val) { fCommonTrackPtCut = val; } | |
100 | void SetCommonSingleCell(const Int_t val) { fCommonSingleCell = val;} | |
101 | // GeometryPhos | |
102 | void SetGeometryPhosEtaAccCut(const Double_t val) { fGeometryPhosEtaAccCut = val; } | |
103 | void SetGeometryPhosPhiAccMinCut(const Double_t val) { fGeometryPhosPhiAccMinCut = val; } | |
104 | void SetGeometryPhosPhiAccMaxCut(const Double_t val) { fGeometryPhosPhiAccMaxCut = val; } | |
105 | void SetGeometryPhosDetectorRadius(const Double_t val) { fGeometryPhosDetectorRadius = val; } | |
106 | // GeometryEmcal | |
107 | void SetGeometryEmcalEtaAccCut(const Double_t val) { fGeometryEmcalEtaAccCut = val; } | |
108 | void SetGeometryEmcalPhiAccMinCut(const Double_t val) { fGeometryEmcalPhiAccMinCut = val; } | |
109 | void SetGeometryEmcalPhiAccMaxCut(const Double_t val) { fGeometryEmcalPhiAccMaxCut = val; } | |
110 | void SetGeometryEmcalDetectorRadius(const Double_t val) { fGeometryEmcalDetectorRadius = val; } | |
111 | // Reconstructed | |
112 | void SetReconstructedVertexXCut(const Double_t val) { fReconstructedVertexXCut = val; } | |
113 | void SetReconstructedVertexYCut(const Double_t val) { fReconstructedVertexYCut = val; } | |
114 | void SetReconstructedVertexZCut(const Double_t val) { fReconstructedVertexZCut = val; } | |
115 | void SetReconstructedIPxyCut(const Double_t val) { fReconstructedIPxyCut = val; } | |
116 | void SetReconstructedIPzCut(const Double_t val) { fReconstructedIPzCut = val; } | |
117 | void SetReconstructedNTpcClustersCut(const Int_t val) { fReconstructedNTpcClustersCut = val; } | |
118 | void SetReconstructedNItsClustersCut(const Int_t val) { fReconstructedNItsClustersCut = val; } | |
e16aac77 | 119 | void SetReconstrucedPidCut(const Double_t val) { fReconstructedPidCut = val; } |
4998becf | 120 | // ReconstructedPhos |
121 | void SetReconstructedPhosClusterType(const Char_t val) { fReconstructedPhosClusterType = val; } | |
122 | void SetReconstructedPhosClusterEnergyCut(const Double_t val) { fReconstructedPhosClusterEnergyCut = val; } | |
123 | void SetReconstructedPhosSingleCellEnergyCut(const Double_t val) { fReconstructedPhosSingleCellEnergyCut = val; } | |
393b8149 | 124 | void SetPhosTrackDistanceCut(const Double_t val) { fPhosTrackDistanceCut = val; } |
125 | void SetPhosTrackDxCut(const Double_t val) { fPhosTrackDxCut = val; } | |
126 | void SetPhosTrackDzCut(const Double_t val) { fPhosTrackDzCut = val; } | |
ef647350 | 127 | void SetPhosTrackRCut(const Double_t val) { fPhosTrackRCut = val; } |
128 | ||
129 | void SetPhosBadDistanceCut(const Double_t val) { fPhosBadDistanceCut = val; } | |
393b8149 | 130 | |
4998becf | 131 | // ReconstructedEmcal |
132 | void SetReconstructedEmcalClusterType(const Char_t val) { fReconstructedEmcalClusterType = val; } | |
133 | void SetReconstructedEmcalClusterEnergyCut(const Double_t val) { fReconstructedEmcalClusterEnergyCut = val; } | |
134 | void SetReconstructedEmcalSingleCellEnergyCut(const Double_t val) { fReconstructedEmcalSingleCellEnergyCut = val; } | |
393b8149 | 135 | void SetEmcalTrackDistanceCut(const Double_t val) { fEmcalTrackDistanceCut = val; } |
4998becf | 136 | // MonteCarlo |
137 | void SetMonteCarloSingleChargedParticle(const Double_t val) { fMonteCarloSingleChargedParticle = val; } | |
138 | void SetMonteCarloNeutralParticle(const Double_t val) { fMonteCarloNeutralParticle = val; } | |
ce546038 | 139 | // Hist: TTree and histogram info |
140 | void SetHistMakeTree(const Bool_t val) { fHistMakeTree = val; } | |
8ea99ab0 | 141 | void SetHistMakeTreeDeposit(const Bool_t val) { fHistMakeTreeDeposit = val; } |
0fa8c632 | 142 | // |
143 | void SetHistNbinsMult(const Int_t val) { fHistNbinsMult = val; } | |
144 | void SetHistMinMult(const Double_t val) { fHistMinMult = val; } | |
145 | void SetHistMaxMult(const Double_t val) { fHistMaxMult = val; } | |
146 | // | |
147 | void SetHistNbinsTotEt(const Int_t val) { fHistNbinsTotEt = val; } | |
148 | void SetHistMinTotEt(const Double_t val) { fHistMinTotEt = val; } | |
149 | void SetHistMaxTotEt(const Double_t val) { fHistMaxTotEt = val; } | |
150 | // | |
151 | void SetHistNbinsParticleEt(const Int_t val) { fHistNbinsParticleEt = val; } | |
152 | void SetHistMinParticleEt(const Double_t val) { fHistMinParticleEt = val; } | |
153 | void SetHistMaxParticleEt(const Double_t val) { fHistMaxParticleEt = val; } | |
154 | // | |
155 | void SetHistNbinsParticlePt(const Int_t val) { fHistNbinsParticlePt = val; } | |
156 | void SetHistMinParticlePt(const Double_t val) { fHistMinParticlePt = val; } | |
157 | void SetHistMaxParticlePt(const Double_t val) { fHistMaxParticlePt = val; } | |
158 | ||
2fbf38ac | 159 | |
4998becf | 160 | protected: |
2fbf38ac | 161 | |
4998becf | 162 | // Common |
163 | Double_t fCommonEtaCut; // Eta cut | |
164 | Double_t fCommonClusterEnergyCut; // Cluster Energy cut | |
165 | Double_t fCommonTrackPtCut; // Track Pt | |
166 | Int_t fCommonSingleCell; // Single Cell (1) | |
393b8149 | 167 | Double_t fEmcalTrackDistanceCut; // EMCal track distance |
168 | Double_t fEmcalTrackDxCut; // EMCal track distance in x | |
169 | Double_t fEmcalTrackDzCut; // EMCal track distance in z | |
170 | ||
171 | Double_t fPhosTrackDistanceCut; // PHOS track distance | |
172 | Double_t fPhosTrackDxCut; // PHOS track distance in x | |
173 | Double_t fPhosTrackDzCut; // PHOS track distance in z | |
ef647350 | 174 | Double_t fPhosTrackRCut; // PHOS track distance in r (using the parametrized track distance) |
175 | ||
176 | Double_t fPhosBadDistanceCut; // PHOS distance to bad channel | |
177 | ||
4998becf | 178 | // GeometryPhos |
179 | Double_t fGeometryPhosEtaAccCut; // PHOS Eta Acc cut | |
180 | Double_t fGeometryPhosPhiAccMinCut; // PHOS Phi Acc Min cut | |
181 | Double_t fGeometryPhosPhiAccMaxCut; // PHOS Phi Acc Max cut | |
182 | Double_t fGeometryPhosDetectorRadius; // PHOS Detector Radius | |
2fbf38ac | 183 | |
4998becf | 184 | // GeometryEmcal |
185 | Double_t fGeometryEmcalEtaAccCut; // EMCal Eta Acc cut | |
186 | Double_t fGeometryEmcalPhiAccMinCut; // EMCal Phi Acc Min cut | |
187 | Double_t fGeometryEmcalPhiAccMaxCut; // EMCal Phi Acc Max cut | |
188 | Double_t fGeometryEmcalDetectorRadius; // EMCal Detector Radius | |
189 | ||
190 | // Reconstructed | |
191 | Double_t fReconstructedVertexXCut; // vertex X cut | |
192 | Double_t fReconstructedVertexYCut; // vertex Y cut | |
193 | Double_t fReconstructedVertexZCut; // vertex Z cut | |
194 | Double_t fReconstructedIPxyCut; // IP xy cut | |
195 | Double_t fReconstructedIPzCut; // IP z cut | |
196 | Int_t fReconstructedNTpcClustersCut; // # of TPC clusters cut | |
197 | Int_t fReconstructedNItsClustersCut; // # of ITS clusters cut | |
e16aac77 | 198 | Double_t fReconstructedPidCut; // cut on pid prob |
4998becf | 199 | |
200 | // ReconstructedPhos | |
201 | Char_t fReconstructedPhosClusterType; // PHOS cluster type | |
202 | Double_t fReconstructedPhosClusterEnergyCut; // PHOS cluster energy | |
203 | Double_t fReconstructedPhosSingleCellEnergyCut; // PHOS single cell energy | |
393b8149 | 204 | Double_t fReconstructedPhosTrackDistanceTightCut; // PHOS track distance |
205 | Double_t fReconstructedPhosTrackDistanceMediumCut; // PHOS track distance | |
206 | Double_t fReconstructedPhosTrackDistanceLooseCut; // PHOS track distance | |
4998becf | 207 | |
208 | // ReconstructedEmcal | |
209 | Char_t fReconstructedEmcalClusterType; // EMCal cluster type | |
210 | Double_t fReconstructedEmcalClusterEnergyCut; // EMCal cluster energy | |
211 | Double_t fReconstructedEmcalSingleCellEnergyCut; // EMCal single cell energy | |
393b8149 | 212 | Double_t fReconstructedEmcalTrackDistanceTightCut; // EMCAL track distance |
213 | Double_t fReconstructedEmcalTrackDistanceMediumCut; // EMCAL track distance | |
214 | Double_t fReconstructedEmcalTrackDistanceLooseCut; // EMCAL track distance | |
4998becf | 215 | |
216 | // MonteCarlo | |
217 | Double_t fMonteCarloSingleChargedParticle; // MC charged | |
218 | Double_t fMonteCarloNeutralParticle; // MC neutral | |
219 | ||
ce546038 | 220 | // Hist: TTree and histogram info |
221 | Bool_t fHistMakeTree; // whether to make a summary tree or not | |
8ea99ab0 | 222 | Bool_t fHistMakeTreeDeposit; // whether to make a summary tree of energy deposit or not |
223 | ||
0fa8c632 | 224 | Int_t fHistNbinsMult; // number of bins in multiplicity histograms |
225 | Double_t fHistMinMult; // minimum value in multiplicity histograms | |
226 | Double_t fHistMaxMult; // maximum value in multiplicity histograms | |
227 | ||
228 | Int_t fHistNbinsTotEt; // number of bins in event Et histograms | |
229 | Double_t fHistMinTotEt; // minimum value in event Et histograms | |
230 | Double_t fHistMaxTotEt; // maximum value in event Et histograms | |
231 | ||
232 | Int_t fHistNbinsParticleEt; // number of bins in particle Et histograms | |
233 | Double_t fHistMinParticleEt; // minimum value in particle Et histograms | |
234 | Double_t fHistMaxParticleEt; // maximum value in particle Et histograms | |
235 | ||
236 | Int_t fHistNbinsParticlePt; // number of bins in particle Pt histograms | |
237 | Double_t fHistMinParticlePt; // minimum value in particle Pt histograms | |
238 | Double_t fHistMaxParticlePt; // maximum value in particle Pt histograms | |
ce546038 | 239 | |
393b8149 | 240 | // Detector definition |
537e541d | 241 | static const Short_t fgkDetectorPhos = -1; // PHOS |
242 | static const Short_t fgkDetectorEmcal = 1; // EMCAL | |
393b8149 | 243 | |
4998becf | 244 | private: |
245 | //Declare private to avoid compilation warning | |
246 | AliAnalysisEtCuts & operator = (const AliAnalysisEtCuts & g) ;//copy assignment | |
247 | AliAnalysisEtCuts(const AliAnalysisEtCuts & g) ; // copy ctor | |
248 | ||
e573e46d | 249 | ClassDef(AliAnalysisEtCuts, 1); |
2fbf38ac | 250 | }; |
251 | ||
4998becf | 252 | #endif // ALIANALYSISETCUTS_H |