]>
Commit | Line | Data |
---|---|---|
1 | #ifndef ALIANAPARTICLEISOLATION_H | |
2 | #define ALIANAPARTICLEISOLATION_H | |
3 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | //_________________________________________________________________________ | |
7 | ||
8 | // Class for the analysis of particle isolation | |
9 | // Input is selected particles put in AOD branch (AliAODPWG4ParticleCorrelation) | |
10 | // | |
11 | // Class created from old AliPHOSGammaJet | |
12 | // (see AliRoot versions previous Release 4-09) | |
13 | ||
14 | //-- Author: Gustavo Conesa (INFN-LNF) | |
15 | ||
16 | // --- ROOT system --- | |
17 | class TH2F; | |
18 | class TList ; | |
19 | class TObjString; | |
20 | ||
21 | // --- ANALYSIS system --- | |
22 | #include "AliAnaCaloTrackCorrBaseClass.h" | |
23 | class AliAODPWG4Particle; | |
24 | class AliAODPWG4ParticleCorrelation ; | |
25 | ||
26 | ||
27 | class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass { | |
28 | ||
29 | public: | |
30 | AliAnaParticleIsolation() ; // default ctor | |
31 | virtual ~AliAnaParticleIsolation() { ; } //virtual dtor | |
32 | ||
33 | // Main general methods | |
34 | ||
35 | void CalculateCaloUEBand (AliAODPWG4ParticleCorrelation * pCandidate, | |
36 | Float_t & etaBand, Float_t & phiBand) ; | |
37 | void CalculateCaloCellUEBand(AliAODPWG4ParticleCorrelation * pCandidate, | |
38 | Float_t & etaBand, Float_t & phiBand) ; | |
39 | void CalculateTrackUEBand (AliAODPWG4ParticleCorrelation * pCandidate, | |
40 | Float_t & etaBand, Float_t & phiBand) ; | |
41 | ||
42 | void CalculateCaloSignalInCone (AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumCluster) ; | |
43 | void CalculateCaloCellSignalInCone(AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumCell) ; | |
44 | void CalculateTrackSignalInCone (AliAODPWG4ParticleCorrelation * aodParticle, Float_t & coneptsumTrack ) ; | |
45 | ||
46 | Float_t CalculateExcessAreaFraction(const Float_t excess, const Float_t conesize); | |
47 | ||
48 | void CalculateNormalizeUEBandPerUnitArea(AliAODPWG4ParticleCorrelation * pCandidate, | |
49 | const Float_t coneptsumCluster,const Float_t coneptsumCell, const Float_t coneptsumTrack ) ; | |
50 | ||
51 | TObjString * GetAnalysisCuts() ; | |
52 | ||
53 | TList * GetCreateOutputObjects() ; | |
54 | ||
55 | void Init() ; | |
56 | ||
57 | void InitParameters() ; | |
58 | ||
59 | void MakeAnalysisFillAOD() ; | |
60 | ||
61 | void MakeAnalysisFillHistograms() ; | |
62 | ||
63 | ||
64 | void Print( const Option_t * opt ) const ; | |
65 | ||
66 | //Analysis specific methods | |
67 | ||
68 | void FillPileUpHistograms(Int_t clusterID) ; | |
69 | ||
70 | void FillTrackMatchingShowerShapeControlHistograms(AliAODPWG4ParticleCorrelation * pCandidate, | |
71 | const AliCaloTrackReader * reader, | |
72 | const AliCaloPID * pid) ; | |
73 | ||
74 | void MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ; | |
75 | ||
76 | // Analysis Setters and Getters | |
77 | ||
78 | TString GetCalorimeter() const { return fCalorimeter ; } | |
79 | Int_t GetNCones() const { return fNCones ; } | |
80 | Int_t GetNPtThresFrac() const { return fNPtThresFrac ; } | |
81 | Float_t GetConeSizes(Int_t i) const { return fConeSizes[i] ; } | |
82 | Float_t GetPtThresholds(Int_t i) const { return fPtThresholds[i] ; } | |
83 | Float_t GetSumPtThresholds(Int_t i) const { return fSumPtThresholds[i]; } | |
84 | Float_t GetPtFractions(Int_t i) const { return fPtFractions[i] ; } | |
85 | ||
86 | void SetCalorimeter(TString & det) { fCalorimeter = det ; } | |
87 | void SetNCones(Int_t ncs) { fNCones = ncs ; } | |
88 | void SetNPtThresFrac(Int_t npt) { fNPtThresFrac = npt ; } | |
89 | void SetConeSizes(Int_t i, Float_t r) { fConeSizes[i] = r ; } | |
90 | void SetPtThresholds(Int_t i, Float_t pt) { fPtThresholds[i] = pt ; } | |
91 | void SetPtFractions(Int_t i, Float_t pt) { fPtFractions[i] = pt ; } | |
92 | void SetSumPtThresholds(Int_t i, Float_t pt){ fSumPtThresholds[i] = pt ; } | |
93 | ||
94 | ||
95 | Bool_t IsReIsolationOn() const { return fReMakeIC ; } | |
96 | void SwitchOnReIsolation() { fReMakeIC = kTRUE ; } | |
97 | void SwitchOffReIsolation() { fReMakeIC = kFALSE ; } | |
98 | ||
99 | Bool_t IsSeveralIsolationOn() const { return fMakeSeveralIC ; } | |
100 | void SwitchOnSeveralIsolation() { fMakeSeveralIC = kTRUE ; } | |
101 | void SwitchOffSeveralIsolation() { fMakeSeveralIC = kFALSE ; } | |
102 | ||
103 | void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; } | |
104 | void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; } | |
105 | ||
106 | void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; } | |
107 | void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; } | |
108 | ||
109 | void SwitchOnSSHistoFill() { fFillSSHisto = kTRUE ; } | |
110 | void SwitchOffSSHistoFill() { fFillSSHisto = kFALSE ; } | |
111 | ||
112 | //Histogrammes setters and getters | |
113 | ||
114 | void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n) { | |
115 | fHistoNPtSumBins = n ; fHistoPtSumMax = max ; fHistoPtSumMin = min ; } | |
116 | ||
117 | Int_t GetHistoNPtSumBins() const { return fHistoNPtSumBins ; } | |
118 | Float_t GetHistoPtSumMin() const { return fHistoPtSumMin ; } | |
119 | Float_t GetHistoPtSumMax() const { return fHistoPtSumMax ; } | |
120 | ||
121 | void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n) { | |
122 | fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min ; } | |
123 | ||
124 | Int_t GetHistoNPtInConeBins() const { return fHistoNPtInConeBins; } | |
125 | Float_t GetHistoPtInConeMin() const { return fHistoPtInConeMin ; } | |
126 | Float_t GetHistoPtInConeMax() const { return fHistoPtInConeMax ; } | |
127 | ||
128 | private: | |
129 | ||
130 | TString fCalorimeter ; // Calorimeter where neutral particles in cone for isolation are; | |
131 | Bool_t fReMakeIC ; // Do isolation analysis | |
132 | Bool_t fMakeSeveralIC ; // Do analysis for different IC | |
133 | Bool_t fFillPileUpHistograms; // Fill pile-up related histograms | |
134 | Bool_t fFillTMHisto; // Fill track matching plots | |
135 | Bool_t fFillSSHisto; // Fill Shower shape plots | |
136 | ||
137 | // Analysis data members for multiple cones and pt thresholds | |
138 | Int_t fNCones ; //! Number of cone sizes to test | |
139 | Int_t fNPtThresFrac ; //! Number of ptThres and ptFrac to test | |
140 | ||
141 | Float_t fConeSizes[5] ; //! Array with cones to test | |
142 | Float_t fPtThresholds[5] ; //! Array with pt thresholds to test | |
143 | Float_t fPtFractions[5] ; //! Array with pt thresholds to test frac | |
144 | Float_t fSumPtThresholds[5] ; //! Array with pt thresholds to test frac | |
145 | ||
146 | //Histograms | |
147 | ||
148 | TH1F * fhEIso ; //! Number of isolated particles vs energy | |
149 | TH1F * fhPtIso ; //! Number of isolated particles vs pT | |
150 | TH2F * fhPtCentralityIso ; //! centrality vs pT | |
151 | TH2F * fhPtEventPlaneIso ; //! event plane angle vs pT | |
152 | TH2F * fhPtNLocMaxIso ; //! Number of isolated particles vs NLM in cluster | |
153 | TH2F * fhPhiIso ; //! Phi of isolated particles | |
154 | TH2F * fhEtaIso ; //! eta of isolated particles | |
155 | TH2F * fhEtaPhiIso ; //! eta vs phi of isolated particles | |
156 | TH2F * fhEtaPhiNoIso ; //! eta vs phi of not isolated leading particles | |
157 | TH1F * fhENoIso ; //! Number of not isolated leading particles vs Energy | |
158 | TH1F * fhPtNoIso ; //! Number of not isolated leading particles vs pT | |
159 | TH2F * fhPtNLocMaxNoIso ; //! Number of not isolated particles vs NLM in cluster | |
160 | TH1F * fhPtDecayIso ; //! Number of isolated Pi0 decay particles (invariant mass tag) | |
161 | TH1F * fhPtDecayNoIso ; //! Number of not isolated Pi0 decay leading particles (invariant mass tag) | |
162 | TH2F * fhEtaPhiDecayIso ; //! eta vs phi of isolated Pi0 decay particles | |
163 | TH2F * fhEtaPhiDecayNoIso ; //! eta vs phi of not isolated leading Pi0 decay particles | |
164 | ||
165 | TH2F * fhPtInCone ; //! Cluster/track Pt in the cone | |
166 | TH2F * fhPtClusterInCone ; //! Cluster Pt in the cone | |
167 | TH2F * fhPtCellInCone ; //! Cell amplitude in the cone | |
168 | TH2F * fhPtTrackInCone ; //! Track Pt in the cone | |
169 | TH2F * fhPtTrackInConeOtherBC ; //! Track Pt in the cone, tracks out of main BC Time window | |
170 | TH2F * fhPtTrackInConeOtherBCPileUpSPD ; //! Track Pt in the cone, tracks out of main BC Time window | |
171 | TH2F * fhPtTrackInConeBC0 ; //! Track Pt in the cone, tracks in BC=0 | |
172 | TH2F * fhPtTrackInConeVtxBC0 ; //! Track Pt in the cone, tracks in BC=0 | |
173 | TH2F * fhPtTrackInConeBC0PileUpSPD ; //! Track Pt in the cone, tracks in BC=0 | |
174 | TH2F * fhPtInConePileUp[7] ; //! Particle Pt in the cone, if event is from pile-up (SPD method) | |
175 | TH2F * fhPtInConeCent ; //! Particle Pt in the cone versus centrality | |
176 | TH2F * fhPerpConeSumPt ; //! Sum Pt in cone at the perpendicular phi region to trigger axis (phi +90) | |
177 | TH2F * fhPtInPerpCone ; //! Particle Pt in cone at the perpendicular phi region to trigger axis (phi +90) | |
178 | ||
179 | TH2F * fhEtaPhiInConeCluster ; //! Eta vs. phi of clusters in cone | |
180 | TH2F * fhEtaPhiCluster ; //! Eta vs. phi of all clusters | |
181 | TH2F * fhEtaPhiInConeTrack ; //! Eta vs. phi of tracks in cone | |
182 | TH2F * fhEtaPhiTrack ; //! Eta vs. phi of all tracks | |
183 | ||
184 | TH2F * fhEtaBandCluster ; //! Accumulated pT in Eta band to estimate UE in cone, only clusters | |
185 | TH2F * fhPhiBandCluster ; //! Accumulated pT in Phi band to estimate UE in cone, only clusters | |
186 | TH2F * fhEtaBandTrack ; //! Accumulated pT in Eta band to estimate UE in cone, only tracks | |
187 | TH2F * fhPhiBandTrack ; //! Accumulated pT in Phi band to estimate UE in cone, only tracks | |
188 | TH2F * fhEtaBandCell ; //! Accumulated pT in Eta band to estimate UE in cone, only cells | |
189 | TH2F * fhPhiBandCell ; //! Accumulated pT in Phi band to estimate UE in cone, only cells | |
190 | ||
191 | TH2F * fhConeSumPt ; //! Cluster and tracks Sum Pt Sum Pt in the cone | |
192 | TH2F * fhConeSumPtCellTrack ; //! Cells and tracks Sum Pt Sum Pt in the cone | |
193 | TH2F * fhConeSumPtCell ; //! Cells Sum Pt Sum Pt in the cone | |
194 | TH2F * fhConeSumPtCluster ; //! Clusters Sum Pt Sum Pt in the cone | |
195 | TH2F * fhConeSumPtTrack ; //! Tracks Sum Pt Sum Pt in the cone | |
196 | TH2F * fhConeSumPtEtaBandUECluster; //! Cluster Sum Pt in the eta band for clusters, before normalization | |
197 | TH2F * fhConeSumPtPhiBandUECluster; //! Cluster Sum Pt in the phi band for clusters, before normalization | |
198 | TH2F * fhConeSumPtEtaBandUETrack; //! Track Sum Pt in the eta band for tracks , before normalization | |
199 | TH2F * fhConeSumPtPhiBandUETrack; //! Track Sum Pt in the phi badn for tracks , before normalization | |
200 | TH2F * fhConeSumPtEtaBandUECell; //! Cell Sum amplitude in the eta band for cells, before normalization | |
201 | TH2F * fhConeSumPtPhiBandUECell; //! Cell Sum amplitude in the phi band for cells, before normalization | |
202 | ||
203 | TH2F * fhConeSumPtTrigEtaPhi ; //! Cluster and tracks Sum Pt Sum Pt in the cone, per eta-phi bin of trigger, | |
204 | TH2F * fhConeSumPtCellTrackTrigEtaPhi ; //! Cell and tracks Sum Pt Sum Pt in the cone, per eta-phi bin of trigger, | |
205 | TH2F * fhConeSumPtEtaBandUEClusterTrigEtaPhi; //! Cluster Sum Pt in the eta band for clusters, per eta-phi bin of trigger,before normalization | |
206 | TH2F * fhConeSumPtPhiBandUEClusterTrigEtaPhi; //! Cluster Sum Pt in the phi band for clusters, per eta-phi bin of trigger, before normalization | |
207 | TH2F * fhConeSumPtEtaBandUETrackTrigEtaPhi; //! Track Sum Pt in the eta band for tracks , per eta-phi bin of trigger, before normalization | |
208 | TH2F * fhConeSumPtPhiBandUETrackTrigEtaPhi; //! Track Sum Pt in the phi badn for tracks , per eta-phi bin of trigger, before normalization | |
209 | TH2F * fhConeSumPtEtaBandUECellTrigEtaPhi; //! Cluster Sum amplitude in the eta band for cells, per eta-phi bin of trigger, before normalization | |
210 | TH2F * fhConeSumPtPhiBandUECellTrigEtaPhi; //! Cluster Sum amplitude in the phi band for cells, per eta-phi bin of trigger, before normalization | |
211 | ||
212 | TH2F * fhConeSumPtEtaUESub; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger | |
213 | TH2F * fhConeSumPtPhiUESub; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger | |
214 | TH2F * fhConeSumPtEtaUESubTrigEtaPhi; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
215 | TH2F * fhConeSumPtPhiUESubTrigEtaPhi; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
216 | ||
217 | TH2F * fhConeSumPtEtaUESubTrackCell; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger | |
218 | TH2F * fhConeSumPtPhiUESubTrackCell; //! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs pT trigger | |
219 | TH2F * fhConeSumPtEtaUESubTrackCellTrigEtaPhi;//! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
220 | TH2F * fhConeSumPtPhiUESubTrackCellTrigEtaPhi;//! Cluster and tracks Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
221 | ||
222 | TH2F * fhConeSumPtEtaUESubCluster; //! Cluster Sum Pt in the cone after bkg subtraction, vs pT trigger | |
223 | TH2F * fhConeSumPtPhiUESubCluster; //! Cluster Sum Pt in the cone after bkg subtraction, vs pT trigger | |
224 | TH2F * fhConeSumPtEtaUESubClusterTrigEtaPhi; //! Cluster Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
225 | TH2F * fhConeSumPtPhiUESubClusterTrigEtaPhi; //! Cluster Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
226 | ||
227 | TH2F * fhConeSumPtEtaUESubCell; //! Cell Sum amplitude in the cone after bkg subtraction, vs pT trigger | |
228 | TH2F * fhConeSumPtPhiUESubCell; //! Cell Sum amplitude in the cone after bkg subtraction, vs pT trigger | |
229 | TH2F * fhConeSumPtEtaUESubCellTrigEtaPhi; //! Cell Sum amplitude in the cone after bkg subtraction, vs eta-phi trigger | |
230 | TH2F * fhConeSumPtPhiUESubCellTrigEtaPhi; //! Cell Sum amplitude in the cone after bkg subtraction, vs eta-phi trigger | |
231 | ||
232 | TH2F * fhConeSumPtEtaUESubTrack; //! Track Sum Pt in the cone after bkg subtraction, vs pT trigger | |
233 | TH2F * fhConeSumPtPhiUESubTrack; //! Track Sum Pt in the cone after bkg subtraction, vs pT trigger | |
234 | TH2F * fhConeSumPtEtaUESubTrackTrigEtaPhi; //! Track Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
235 | TH2F * fhConeSumPtPhiUESubTrackTrigEtaPhi; //! Track Sum Pt in the cone after bkg subtraction, vs eta-phi trigger | |
236 | ||
237 | TH2F * fhFractionTrackOutConeEta; //! Fraction of cone out of tracks acceptance in eta | |
238 | TH2F * fhFractionTrackOutConeEtaTrigEtaPhi; //! Fraction of cone out of tracks acceptance in eta, vs trigger eta-phi | |
239 | TH2F * fhFractionClusterOutConeEta; //! Fraction of cone out of clusters acceptance in eta | |
240 | TH2F * fhFractionClusterOutConeEtaTrigEtaPhi; //! Fraction of cone out of clusters acceptance in eta, vs trigger eta-phi | |
241 | TH2F * fhFractionClusterOutConePhi; //! Fraction of cone out of clusters acceptance in phi | |
242 | TH2F * fhFractionClusterOutConePhiTrigEtaPhi; //! Fraction of cone out of clusters acceptance in phi, vs trigger eta-phi | |
243 | ||
244 | TH2F * fhFractionCellOutConeEta; //! Fraction of cone out of cells acceptance in eta | |
245 | TH2F * fhFractionCellOutConeEtaTrigEtaPhi; //! Fraction of cone out of cells acceptance in eta, vs trigger eta-phi | |
246 | TH2F * fhFractionCellOutConePhi; //! Fraction of cone out of cells acceptance in phi | |
247 | TH2F * fhFractionCellOutConePhiTrigEtaPhi; //! Fraction of cone out of cells acceptance in phi, vs trigger eta-phi | |
248 | ||
249 | TH2F * fhConeSumPtClustervsTrack ; //! Cluster vs tracks Sum Pt Sum Pt in the cone | |
250 | TH2F * fhConeSumPtEtaUESubClustervsTrack ; //! Cluster vs tracks Sum Pt Sum Pt in the cone, after subtraction in eta band | |
251 | TH2F * fhConeSumPtPhiUESubClustervsTrack ; //! Cluster vs tracks Sum Pt Sum Pt in the cone, after subtraction in phi band | |
252 | TH2F * fhConeSumPtCellvsTrack; //! Cell vs tracks Sum Pt Sum Pt in the cone | |
253 | TH2F * fhConeSumPtEtaUESubCellvsTrack ; //! Cell vs tracks Sum Pt Sum Pt in the cone, after subtraction in eta band | |
254 | TH2F * fhConeSumPtPhiUESubCellvsTrack ; //! Cell vs tracks Sum Pt Sum Pt in the cone, after subtraction in phi band | |
255 | ||
256 | TH2F * fhEtaBandClustervsTrack ; //! Accumulated pT in Eta band to estimate UE in cone, clusters vs tracks | |
257 | TH2F * fhPhiBandClustervsTrack ; //! Accumulated pT in Phi band to estimate UE in cone, clusters vs tracks | |
258 | TH2F * fhEtaBandNormClustervsTrack ; //! Accumulated pT in Eta band to estimate UE in cone, normalized to cone size, clusters vs tracks | |
259 | TH2F * fhPhiBandNormClustervsTrack ; //! Accumulated pT in Phi band to estimate UE in cone, normalized to cone size, clusters vs tracks | |
260 | TH2F * fhEtaBandCellvsTrack ; //! Accumulated pT in Eta band to estimate UE in cone, cells vs tracks | |
261 | TH2F * fhPhiBandCellvsTrack ; //! Accumulated pT in Phi band to estimate UE in cone, cells vs tracks | |
262 | TH2F * fhEtaBandNormCellvsTrack ; //! Accumulated pT cell in Eta band to estimate UE in cone, normalized to cone size, clusters vs tracks | |
263 | TH2F * fhPhiBandNormCellvsTrack ; //! Accumulated pT cell in Phi band to estimate UE in cone, normalized to cone | |
264 | ||
265 | TH2F * fhConeSumPtSubvsConeSumPtTotPhiTrack; //! Tracks, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
266 | TH2F * fhConeSumPtSubNormvsConeSumPtTotPhiTrack; //! Tracks, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
267 | TH2F * fhConeSumPtSubvsConeSumPtTotEtaTrack; //! Tracks, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
268 | TH2F * fhConeSumPtSubNormvsConeSumPtTotEtaTrack; //! Tracks, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
269 | TH2F * fhConeSumPtSubvsConeSumPtTotPhiCluster; //! Clusters, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
270 | TH2F * fhConeSumPtSubNormvsConeSumPtTotPhiCluster; //! Clusters, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
271 | TH2F * fhConeSumPtSubvsConeSumPtTotEtaCluster; //! Clusters, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
272 | TH2F * fhConeSumPtSubNormvsConeSumPtTotEtaCluster; //! Clusters, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
273 | TH2F * fhConeSumPtSubvsConeSumPtTotPhiCell; //! Cells, phi band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
274 | TH2F * fhConeSumPtSubNormvsConeSumPtTotPhiCell; //! Cells, phi band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
275 | TH2F * fhConeSumPtSubvsConeSumPtTotEtaCell; //! Cells, eta band: sum pT in cone after bkg sub vs sum pT in cone before bkg sub | |
276 | TH2F * fhConeSumPtSubNormvsConeSumPtTotEtaCell; //! Cells, eta band: sum pT in cone after bkg sub normalized by sum pT in cone before bkg sub vs sum pT in cone before bkg sub | |
277 | ||
278 | ||
279 | //MC | |
280 | TH1F * fhPtIsoPrompt; //! Number of isolated prompt gamma | |
281 | TH2F * fhPhiIsoPrompt; //! Phi of isolated prompt gamma | |
282 | TH2F * fhEtaIsoPrompt; //! eta of isolated prompt gamma | |
283 | TH1F * fhPtThresIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt threshold | |
284 | TH1F * fhPtFracIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt frac | |
285 | TH2F * fhPtSumIsolatedPrompt[5]; //! Isolated prompt gamma with threshold on cone pt sume | |
286 | TH1F * fhPtIsoFragmentation; //! Number of isolated fragmentation gamma | |
287 | TH2F * fhPhiIsoFragmentation; //! Phi of isolated fragmentation gamma | |
288 | TH2F * fhEtaIsoFragmentation; //! eta of isolated fragmentation gamma | |
289 | TH1F * fhPtThresIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt threshold | |
290 | TH1F * fhPtFracIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt frac | |
291 | TH2F * fhPtSumIsolatedFragmentation[5]; //! Isolated fragmentation gamma with threshold on cone pt sume | |
292 | TH1F * fhPtIsoPi0; //! Number of isolated pi0 (2 gamma) | |
293 | TH2F * fhPhiIsoPi0; //! Phi of isolated pi0 (2 gamma) | |
294 | TH2F * fhEtaIsoPi0; //! eta of isolated pi0 (2 gamma) | |
295 | TH1F * fhPtThresIsolatedPi0[5][5]; //! Isolated pi0 (2 gamma) with pt threshold | |
296 | TH1F * fhPtFracIsolatedPi0[5][5]; //! Isolated pi0 (2 gamma) with pt frac | |
297 | TH2F * fhPtSumIsolatedPi0[5]; //! Isolated pi0 (2 gamma) with threshold on cone pt sum | |
298 | TH1F * fhPtIsoPi0Decay; //! Number of isolated pi0 decay gamma | |
299 | TH2F * fhPhiIsoPi0Decay; //! Phi of isolated pi0 decay gamma | |
300 | TH2F * fhEtaIsoPi0Decay; //! eta of isolated pi0 decay gamma | |
301 | TH1F * fhPtThresIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt threshold | |
302 | TH1F * fhPtFracIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt frac | |
303 | TH2F * fhPtSumIsolatedPi0Decay[5]; //! Isolated pi0 decay gamma with threshold on cone pt sume | |
304 | TH1F * fhPtIsoEtaDecay; //! Number of isolated eta decay gamma | |
305 | TH2F * fhPhiIsoEtaDecay; //! Phi of isolated eta decay gamma | |
306 | TH2F * fhEtaIsoEtaDecay; //! eta of isolated eta decay gamma | |
307 | TH1F * fhPtThresIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt threshold | |
308 | TH1F * fhPtFracIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt frac | |
309 | TH2F * fhPtSumIsolatedEtaDecay[5]; //! Isolated eta fecay gamma with threshold on cone pt sume | |
310 | TH1F * fhPtIsoOtherDecay; //! Number of isolated other decay gamma | |
311 | TH2F * fhPhiIsoOtherDecay; //! Phi of isolated other decay gamma | |
312 | TH2F * fhEtaIsoOtherDecay; //! eta of isolated other decay gamma | |
313 | TH1F * fhPtThresIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt threshold | |
314 | TH1F * fhPtFracIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt frac | |
315 | TH2F * fhPtSumIsolatedOtherDecay[5]; //! Isolated OtherDecay gamma with threshold on cone pt sume | |
316 | // TH1F * fhPtIsoConversion; //! Number of isolated Conversion gamma | |
317 | // TH2F * fhPhiIsoConversion; //! Phi of isolated Conversion gamma | |
318 | // TH2F * fhEtaIsoConversion; //! eta of isolated Conversion gamma | |
319 | // TH1F * fhPtThresIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt threshold | |
320 | // TH1F * fhPtFracIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt frac | |
321 | // TH2F * fhPtSumIsolatedConversion[5]; //! Isolated Conversion gamma with threshold on cone pt sume | |
322 | TH1F * fhPtIsoHadron; //! Number of isolated Hadron | |
323 | TH2F * fhPhiIsoHadron; //! Phi of isolated Hadron | |
324 | TH2F * fhEtaIsoHadron; //! eta of isolated Hadron | |
325 | TH1F * fhPtThresIsolatedHadron[5][5]; //! Isolated Hadron gamma with pt threshold | |
326 | TH1F * fhPtFracIsolatedHadron[5][5]; //! Isolated Hadron gamma with pt frac | |
327 | TH2F * fhPtSumIsolatedHadron[5]; //! Isolated Hadron gamma with threshold on cone pt sume | |
328 | ||
329 | // Multi Cut analysis Several IC | |
330 | TH1F * fhPtNoIsoPi0; //! Number of not isolated leading pi0 (2 gamma) | |
331 | TH1F * fhPtNoIsoPi0Decay; //! Number of not isolated leading pi0 decay gamma | |
332 | TH1F * fhPtNoIsoEtaDecay; //! Number of not isolated leading eta decay gamma | |
333 | TH1F * fhPtNoIsoOtherDecay; //! Number of not isolated leading other decay gamma | |
334 | TH1F * fhPtNoIsoPrompt; //! Number of not isolated leading prompt gamma | |
335 | TH1F * fhPtIsoMCPhoton; //! Number of isolated leading gamma | |
336 | TH1F * fhPtNoIsoMCPhoton; //! Number of not isolated leading gamma | |
337 | // TH1F * fhPtNoIsoConversion; //! Number of not isolated leading conversion gamma | |
338 | TH1F * fhPtNoIsoFragmentation; //! Number of not isolated leading fragmentation gamma | |
339 | TH1F * fhPtNoIsoHadron; //! Number of not isolated leading hadrons | |
340 | ||
341 | TH2F * fhSumPtLeadingPt[5] ; //! Sum Pt in the cone | |
342 | TH2F * fhPtLeadingPt[5] ; //! Particle Pt in the cone | |
343 | TH2F * fhPerpSumPtLeadingPt[5] ; //! Sum Pt in the cone at the perpendicular phi region to trigger axis (phi +90) | |
344 | TH2F * fhPerpPtLeadingPt[5]; //! Sum Pt in the cone at the perpendicular phi region to trigger axis (phi +90) | |
345 | ||
346 | TH1F * fhPtThresIsolated[5][5] ; //! Isolated particle with pt threshold | |
347 | TH1F * fhPtFracIsolated[5][5] ; //! Isolated particle with pt threshold frac | |
348 | TH1F * fhPtSumIsolated[5][5] ; //! Isolated particle with threshold on cone pt sum | |
349 | ||
350 | TH2F * fhEtaPhiPtThresIso[5][5] ; //! eta vs phi of isolated particles with pt threshold | |
351 | TH2F * fhEtaPhiPtThresDecayIso[5][5] ; //! eta vs phi of isolated particles with pt threshold | |
352 | TH1F * fhPtPtThresDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt threshold | |
353 | ||
354 | TH2F * fhEtaPhiPtFracIso[5][5] ; //! eta vs phi of isolated particles with pt frac | |
355 | TH2F * fhEtaPhiPtFracDecayIso[5][5] ; //! eta vs phi of isolated particles with pt frac | |
356 | TH1F * fhPtPtFracDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt fra | |
357 | ||
358 | TH2F * fhEtaPhiPtSumIso[5][5] ; //! eta vs phi of isolated particles with pt sum | |
359 | TH2F * fhEtaPhiPtSumDecayIso[5][5] ; //! eta vs phi of isolated particles with pt sum | |
360 | TH1F * fhPtPtSumDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum | |
361 | ||
362 | TH2F * fhEtaPhiSumDensityIso[5][5]; //! Isolated particle with threshold on cone sum density | |
363 | TH2F * fhEtaPhiSumDensityDecayIso[5][5]; //! Isolated particle with threshold on cone sum density | |
364 | TH1F * fhPtSumDensityIso[5][5]; //! Isolated particle with threshold on cone sum density | |
365 | TH1F * fhPtSumDensityDecayIso[5][5]; //! Isolated decay particle with threshold on cone sum density | |
366 | ||
367 | TH1F * fhPtFracPtSumIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum | |
368 | TH1F * fhPtFracPtSumDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum | |
369 | TH2F * fhEtaPhiFracPtSumIso[5][5]; //! Isolated particle with threshold on cone sum density | |
370 | TH2F * fhEtaPhiFracPtSumDecayIso[5][5]; //! Isolated particle with threshold on cone sum density | |
371 | ||
372 | // Track matching studies | |
373 | TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E | |
374 | TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E | |
375 | TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV | |
376 | TH2F * fhdEdx[2] ; //! matched track dEdx vs cluster E | |
377 | TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut | |
378 | TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle | |
379 | ||
380 | // Shower Shape histograms | |
381 | TH2F * fhELambda0[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously) | |
382 | TH2F * fhELambda1[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously) | |
383 | TH2F * fhELambda0SSBkg; //! Shower shape of non isolated photons close to isolation threshold (do not apply SS cut previously) | |
384 | TH2F * fhELambda0TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously) | |
385 | TH2F * fhELambda1TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously) | |
386 | TH2F * fhELambda0MCPhoton[2]; //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously) | |
387 | TH2F * fhELambda0MCPi0[2]; //! Shower shape of (non) isolated photon candidates originated by particle 2 merged pi0 photons (do not apply SS cut previously) | |
388 | TH2F * fhELambda0MCPi0Decay[2]; //! Shower shape of (non) isolated photon candidates originated by particle pi0 decay photon (do not apply SS cut previously) | |
389 | TH2F * fhELambda0MCEtaDecay[2]; //! Shower shape of (non) isolated photon candidates originated by particle eta decay photon (do not apply SS cut previously) | |
390 | TH2F * fhELambda0MCOtherDecay[2]; //! Shower shape of (non) isolated photon candidates originated by particle other decay photon (do not apply SS cut previously) | |
391 | TH2F * fhELambda0MCHadron[2]; //! Shower shape of (non) isolated photon candidates originated by particle other hadrons (do not apply SS cut previously) | |
392 | ||
393 | // Local maxima | |
394 | TH2F * fhNLocMax[2]; //! number of maxima in selected clusters | |
395 | TH2F * fhELambda0LocMax1[2] ; //! E vs lambda0 of selected cluster, 1 local maxima in cluster | |
396 | TH2F * fhELambda1LocMax1[2] ; //! E vs lambda1 of selected cluster, 1 local maxima in cluster | |
397 | TH2F * fhELambda0LocMax2[2] ; //! E vs lambda0 of selected cluster, 2 local maxima in cluster | |
398 | TH2F * fhELambda1LocMax2[2] ; //! E vs lambda1 of selected cluster, 2 local maxima in cluster | |
399 | TH2F * fhELambda0LocMaxN[2] ; //! E vs lambda0 of selected cluster, N>2 local maxima in cluster | |
400 | TH2F * fhELambda1LocMaxN[2] ; //! E vs lambda1 of selected cluster, N>2 local maxima in cluster | |
401 | ||
402 | // Pile-up | |
403 | TH1F * fhEIsoPileUp[7] ; //! Number of isolated particles | |
404 | TH1F * fhPtIsoPileUp[7] ; //! Number of isolated particles | |
405 | TH1F * fhENoIsoPileUp[7] ; //! Number of not isolated particles | |
406 | TH1F * fhPtNoIsoPileUp[7] ; //! Number of not isolated particles | |
407 | TH2F * fhTimeENoCut; //! time of cluster vs E, no cut | |
408 | TH2F * fhTimeESPD; //! time of cluster vs E, IsSPDPileUp | |
409 | TH2F * fhTimeESPDMulti; //! time of cluster vs E, IsSPDPileUpMulti | |
410 | TH2F * fhTimeNPileUpVertSPD; //! time of cluster vs n pile-up vertices from SPD | |
411 | TH2F * fhTimeNPileUpVertTrack; //! time of cluster vs n pile-up vertices from Tracks | |
412 | TH2F * fhTimeNPileUpVertContributors; //! time of cluster vs n pile-up vertex from SPD contributors | |
413 | TH2F * fhTimePileUpMainVertexZDistance; //! time of cluster vs difference of z main vertex and pile-up vertex | |
414 | TH2F * fhTimePileUpMainVertexZDiamond; //! time of cluster vs difference of z diamond and pile-up vertex | |
415 | ||
416 | //Histograms settings | |
417 | Int_t fHistoNPtSumBins; // Number of bins in PtSum histograms | |
418 | Float_t fHistoPtSumMax; // PtSum maximum in histogram | |
419 | Float_t fHistoPtSumMin; // PtSum minimum in histogram | |
420 | Int_t fHistoNPtInConeBins; // Number of bins in PtInCone histogram | |
421 | Float_t fHistoPtInConeMax; // PtInCone maximum in histogram | |
422 | Float_t fHistoPtInConeMin; // PtInCone maximum in histogram | |
423 | ||
424 | AliAnaParticleIsolation( const AliAnaParticleIsolation & iso) ; // cpy ctor | |
425 | AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment | |
426 | ||
427 | ClassDef(AliAnaParticleIsolation,23) | |
428 | } ; | |
429 | ||
430 | ||
431 | #endif //ALIANAPARTICLEISOLATION_H | |
432 | ||
433 | ||
434 |