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