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