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
36 TObjString * GetAnalysisCuts() ;
38 TList * GetCreateOutputObjects() ;
42 void InitParameters() ;
44 void MakeAnalysisFillAOD() ;
46 void MakeAnalysisFillHistograms() ;
48 void Print( const Option_t * opt ) const ;
50 //Analysis specific methods
52 void FillPileUpHistograms(Int_t clusterID) ;
54 void FillTrackMatchingShowerShapeControlHistograms(const Bool_t isolated,
55 const Int_t clusterID,
58 const TObjArray * plCTS,
59 const TObjArray * plNe,
60 AliAODPWG4ParticleCorrelation * pCandidate,
61 const AliCaloTrackReader * reader,
62 const AliCaloPID * pid) ;
66 void MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ;
68 // Analysis Setters and Getters
70 TString GetCalorimeter() const { return fCalorimeter ; }
71 Int_t GetNCones() const { return fNCones ; }
72 Int_t GetNPtThresFrac() const { return fNPtThresFrac ; }
73 Float_t GetConeSizes(Int_t i) const { return fConeSizes[i] ; }
74 Float_t GetPtThresholds(Int_t i) const { return fPtThresholds[i] ; }
75 Float_t GetSumPtThresholds(Int_t i) const { return fSumPtThresholds[i]; }
76 Float_t GetPtFractions(Int_t i) const { return fPtFractions[i] ; }
78 void SetCalorimeter(TString & det) { fCalorimeter = det ; }
79 void SetNCones(Int_t ncs) { fNCones = ncs ; }
80 void SetNPtThresFrac(Int_t npt) { fNPtThresFrac = npt ; }
81 void SetConeSizes(Int_t i, Float_t r) { fConeSizes[i] = r ; }
82 void SetPtThresholds(Int_t i, Float_t pt) { fPtThresholds[i] = pt ; }
83 void SetPtFractions(Int_t i, Float_t pt) { fPtFractions[i] = pt ; }
84 void SetSumPtThresholds(Int_t i, Float_t pt){ fSumPtThresholds[i] = pt ; }
87 Bool_t IsReIsolationOn() const { return fReMakeIC ; }
88 void SwitchOnReIsolation() { fReMakeIC = kTRUE ; }
89 void SwitchOffReIsolation() { fReMakeIC = kFALSE ; }
91 Bool_t IsSeveralIsolationOn() const { return fMakeSeveralIC ; }
92 void SwitchOnSeveralIsolation() { fMakeSeveralIC = kTRUE ; }
93 void SwitchOffSeveralIsolation() { fMakeSeveralIC = kFALSE ; }
95 void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; }
96 void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; }
98 void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
99 void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
101 void SwitchOnSSHistoFill() { fFillSSHisto = kTRUE ; }
102 void SwitchOffSSHistoFill() { fFillSSHisto = kFALSE ; }
104 //Histogrammes setters and getters
106 void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n) {
107 fHistoNPtSumBins = n ; fHistoPtSumMax = max ; fHistoPtSumMin = min ; }
109 Int_t GetHistoNPtSumBins() const { return fHistoNPtSumBins ; }
110 Float_t GetHistoPtSumMin() const { return fHistoPtSumMin ; }
111 Float_t GetHistoPtSumMax() const { return fHistoPtSumMax ; }
113 void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
114 fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min ; }
116 Int_t GetHistoNPtInConeBins() const { return fHistoNPtInConeBins; }
117 Float_t GetHistoPtInConeMin() const { return fHistoPtInConeMin ; }
118 Float_t GetHistoPtInConeMax() const { return fHistoPtInConeMax ; }
122 TString fCalorimeter ; // Calorimeter where neutral particles in cone for isolation are;
123 Bool_t fReMakeIC ; // Do isolation analysis
124 Bool_t fMakeSeveralIC ; // Do analysis for different IC
125 Bool_t fFillPileUpHistograms; // Fill pile-up related histograms
126 Bool_t fFillTMHisto; // Fill track matching plots
127 Bool_t fFillSSHisto; // Fill Shower shape plots
129 // Analysis data members for multiple cones and pt thresholds
130 Int_t fNCones ; //! Number of cone sizes to test
131 Int_t fNPtThresFrac ; //! Number of ptThres and ptFrac to test
133 Float_t fConeSizes[5] ; //! Array with cones to test
134 Float_t fPtThresholds[5] ; //! Array with pt thresholds to test
135 Float_t fPtFractions[5] ; //! Array with pt thresholds to test frac
136 Float_t fSumPtThresholds[5] ; //! Array with pt thresholds to test frac
140 TH1F * fhEIso ; //! Number of isolated particles
141 TH1F * fhPtIso ; //! Number of isolated particles
142 TH2F * fhPhiIso ; //! Phi of isolated particles
143 TH2F * fhEtaIso ; //! eta of isolated particles
144 TH2F * fhEtaPhiIso ; //! eta vs phi of isolated particles
145 TH2F * fhEtaPhiNoIso ; //! eta vs phi of not isolated leading particles
146 TH1F * fhPtNoIso ; //! Number of not isolated leading particles
147 TH1F * fhPtDecayIso ; //! Number of isolated Pi0 decay particles (invariant mass tag)
148 TH1F * fhPtDecayNoIso ; //! Number of not isolated Pi0 decay leading particles (invariant mass tag)
149 TH2F * fhEtaPhiDecayIso ; //! eta vs phi of isolated Pi0 decay particles
150 TH2F * fhEtaPhiDecayNoIso ; //! eta vs phi of not isolated leading Pi0 decay particles
151 TH2F * fhConeSumPt ; //! Sum Pt in the cone
152 TH2F * fhPtInCone ; //! Particle Pt in the cone
153 TH2F * fhPtInConeCent ; //! Particle Pt in the cone versus centrality
154 TH2F * fhFRConeSumPt ; //! Sum Pt in the forward region cone (phi +90)
155 TH2F * fhPtInFRCone ; //! Particle Pt in the forward region cone (phi +90 )
156 TH2F * fhPhiUEConeSumPt ; //! UE Pt sum in phi band around cone
157 TH2F * fhEtaUEConeSumPt ; //! UE Pt sum in eta band around cone
158 TH2F * fhEtaBand ; //! Eta band to estimate UE in cone
159 TH2F * fhPhiBand ; //! Phi band to estimate UE in cone
160 TH2F * fhConeSumPtEtaUESub; //! Sum Pt in the cone after bkg subtraction
161 TH2F * fhConeSumPtPhiUESub; //! Sum Pt in the cone after bkg subtraction
164 TH1F * fhPtIsoPrompt; //! Number of isolated prompt gamma
165 TH2F * fhPhiIsoPrompt; //! Phi of isolated prompt gamma
166 TH2F * fhEtaIsoPrompt; //! eta of isolated prompt gamma
167 TH1F * fhPtThresIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt threshold
168 TH1F * fhPtFracIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt frac
169 TH2F * fhPtSumIsolatedPrompt[5]; //! Isolated prompt gamma with threshold on cone pt sume
170 TH1F * fhPtIsoFragmentation; //! Number of isolated fragmentation gamma
171 TH2F * fhPhiIsoFragmentation; //! Phi of isolated fragmentation gamma
172 TH2F * fhEtaIsoFragmentation; //! eta of isolated fragmentation gamma
173 TH1F * fhPtThresIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt threshold
174 TH1F * fhPtFracIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt frac
175 TH2F * fhPtSumIsolatedFragmentation[5]; //! Isolated fragmentation gamma with threshold on cone pt sume
176 TH1F * fhPtIsoPi0; //! Number of isolated pi0 (2 gamma)
177 TH2F * fhPhiIsoPi0; //! Phi of isolated pi0 (2 gamma)
178 TH2F * fhEtaIsoPi0; //! eta of isolated pi0 (2 gamma)
179 TH1F * fhPtThresIsolatedPi0[5][5]; //! Isolated pi0 (2 gamma) with pt threshold
180 TH1F * fhPtFracIsolatedPi0[5][5]; //! Isolated pi0 (2 gamma) with pt frac
181 TH2F * fhPtSumIsolatedPi0[5]; //! Isolated pi0 (2 gamma) with threshold on cone pt sum
182 TH1F * fhPtIsoPi0Decay; //! Number of isolated pi0 decay gamma
183 TH2F * fhPhiIsoPi0Decay; //! Phi of isolated pi0 decay gamma
184 TH2F * fhEtaIsoPi0Decay; //! eta of isolated pi0 decay gamma
185 TH1F * fhPtThresIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt threshold
186 TH1F * fhPtFracIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt frac
187 TH2F * fhPtSumIsolatedPi0Decay[5]; //! Isolated pi0 decay gamma with threshold on cone pt sume
188 TH1F * fhPtIsoEtaDecay; //! Number of isolated eta decay gamma
189 TH2F * fhPhiIsoEtaDecay; //! Phi of isolated eta decay gamma
190 TH2F * fhEtaIsoEtaDecay; //! eta of isolated eta decay gamma
191 TH1F * fhPtThresIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt threshold
192 TH1F * fhPtFracIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt frac
193 TH2F * fhPtSumIsolatedEtaDecay[5]; //! Isolated eta fecay gamma with threshold on cone pt sume
194 TH1F * fhPtIsoOtherDecay; //! Number of isolated other decay gamma
195 TH2F * fhPhiIsoOtherDecay; //! Phi of isolated other decay gamma
196 TH2F * fhEtaIsoOtherDecay; //! eta of isolated other decay gamma
197 TH1F * fhPtThresIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt threshold
198 TH1F * fhPtFracIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt frac
199 TH2F * fhPtSumIsolatedOtherDecay[5]; //! Isolated OtherDecay gamma with threshold on cone pt sume
200 // TH1F * fhPtIsoConversion; //! Number of isolated Conversion gamma
201 // TH2F * fhPhiIsoConversion; //! Phi of isolated Conversion gamma
202 // TH2F * fhEtaIsoConversion; //! eta of isolated Conversion gamma
203 // TH1F * fhPtThresIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt threshold
204 // TH1F * fhPtFracIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt frac
205 // TH2F * fhPtSumIsolatedConversion[5]; //! Isolated Conversion gamma with threshold on cone pt sume
206 TH1F * fhPtIsoHadron; //! Number of isolated Hadron
207 TH2F * fhPhiIsoHadron; //! Phi of isolated Hadron
208 TH2F * fhEtaIsoHadron; //! eta of isolated Hadron
209 TH1F * fhPtThresIsolatedHadron[5][5]; //! Isolated Hadron gamma with pt threshold
210 TH1F * fhPtFracIsolatedHadron[5][5]; //! Isolated Hadron gamma with pt frac
211 TH2F * fhPtSumIsolatedHadron[5]; //! Isolated Hadron gamma with threshold on cone pt sume
213 // Multi Cut analysis Several IC
214 TH1F * fhPtNoIsoPi0; //! Number of not isolated leading pi0 (2 gamma)
215 TH1F * fhPtNoIsoPi0Decay; //! Number of not isolated leading pi0 decay gamma
216 TH1F * fhPtNoIsoEtaDecay; //! Number of not isolated leading eta decay gamma
217 TH1F * fhPtNoIsoOtherDecay; //! Number of not isolated leading other decay gamma
218 TH1F * fhPtNoIsoPrompt; //! Number of not isolated leading prompt gamma
219 TH1F * fhPtIsoMCPhoton; //! Number of isolated leading gamma
220 TH1F * fhPtNoIsoMCPhoton; //! Number of not isolated leading gamma
221 // TH1F * fhPtNoIsoConversion; //! Number of not isolated leading conversion gamma
222 TH1F * fhPtNoIsoFragmentation; //! Number of not isolated leading fragmentation gamma
223 TH1F * fhPtNoIsoHadron; //! Number of not isolated leading hadrons
225 TH2F * fhSumPtLeadingPt[5] ; //! Sum Pt in the cone
226 TH2F * fhPtLeadingPt[5] ; //! Particle Pt in the cone
227 TH2F * fhFRSumPtLeadingPt[5] ; //! Sum Pt in the forward region cone (phi +90)
228 TH2F * fhFRPtLeadingPt[5]; //! Particle Pt in the forward region cone (phi +90 )
230 TH1F * fhPtThresIsolated[5][5] ; //! Isolated particle with pt threshold
231 TH1F * fhPtFracIsolated[5][5] ; //! Isolated particle with pt threshold frac
232 TH1F * fhPtSumIsolated[5][5] ; //! Isolated particle with threshold on cone pt sum
234 TH2F * fhEtaPhiPtThresIso[5][5] ; //! eta vs phi of isolated particles with pt threshold
235 TH2F * fhEtaPhiPtThresDecayIso[5][5] ; //! eta vs phi of isolated particles with pt threshold
236 TH1F * fhPtPtThresDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt threshold
238 TH2F * fhEtaPhiPtFracIso[5][5] ; //! eta vs phi of isolated particles with pt frac
239 TH2F * fhEtaPhiPtFracDecayIso[5][5] ; //! eta vs phi of isolated particles with pt frac
240 TH1F * fhPtPtFracDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt fra
242 TH2F * fhEtaPhiPtSumIso[5][5] ; //! eta vs phi of isolated particles with pt sum
243 TH2F * fhEtaPhiPtSumDecayIso[5][5] ; //! eta vs phi of isolated particles with pt sum
244 TH1F * fhPtPtSumDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
246 TH2F * fhEtaPhiSumDensityIso[5][5]; //! Isolated particle with threshold on cone sum density
247 TH2F * fhEtaPhiSumDensityDecayIso[5][5]; //! Isolated particle with threshold on cone sum density
248 TH1F * fhPtSumDensityIso[5][5]; //! Isolated particle with threshold on cone sum density
249 TH1F * fhPtSumDensityDecayIso[5][5]; //! Isolated decay particle with threshold on cone sum density
251 TH1F * fhPtFracPtSumIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
252 TH1F * fhPtFracPtSumDecayIso[5][5] ; //! Number of isolated Pi0 decay particles (invariant mass tag) with pt sum
253 TH2F * fhEtaPhiFracPtSumIso[5][5]; //! Isolated particle with threshold on cone sum density
254 TH2F * fhEtaPhiFracPtSumDecayIso[5][5]; //! Isolated particle with threshold on cone sum density
256 // Track matching studies
257 TH2F * fhTrackMatchedDEta[2] ; //! Eta distance between track and cluster vs cluster E
258 TH2F * fhTrackMatchedDPhi[2] ; //! Phi distance between track and cluster vs cluster E
259 TH2F * fhTrackMatchedDEtaDPhi[2] ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
260 TH2F * fhdEdx[2] ; //! matched track dEdx vs cluster E
261 TH2F * fhEOverP[2]; //! matched track E cluster over P track vs cluster E, after dEdx cut
262 TH2F * fhTrackMatchedMCParticle[2]; //! Trace origin of matched particle
264 // Shower Shape histograms
265 TH2F * fhELambda0[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously)
266 TH2F * fhELambda1[2]; //! Shower shape of (non) isolated photons (do not apply SS cut previously)
267 TH2F * fhELambda0SSBkg; //! Shower shape of non isolated photons close to isolation threshold (do not apply SS cut previously)
268 TH2F * fhELambda0TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
269 TH2F * fhELambda1TRD[2]; //! Shower shape of (non) isolated photons, SM behind TRD (do not apply SS cut previously)
270 TH2F * fhELambda0MCPhoton[2]; //! Shower shape of (non) isolated photon candidates originated by particle photon not decay (do not apply SS cut previously)
271 TH2F * fhELambda0MCPi0[2]; //! Shower shape of (non) isolated photon candidates originated by particle 2 merged pi0 photons (do not apply SS cut previously)
272 TH2F * fhELambda0MCPi0Decay[2]; //! Shower shape of (non) isolated photon candidates originated by particle pi0 decay photon (do not apply SS cut previously)
273 TH2F * fhELambda0MCEtaDecay[2]; //! Shower shape of (non) isolated photon candidates originated by particle eta decay photon (do not apply SS cut previously)
274 TH2F * fhELambda0MCOtherDecay[2]; //! Shower shape of (non) isolated photon candidates originated by particle other decay photon (do not apply SS cut previously)
275 TH2F * fhELambda0MCHadron[2]; //! Shower shape of (non) isolated photon candidates originated by particle other hadrons (do not apply SS cut previously)
278 TH2F * fhNLocMax[2]; //! number of maxima in selected clusters
279 TH2F * fhELambda0LocMax1[2] ; //! E vs lambda0 of selected cluster, 1 local maxima in cluster
280 TH2F * fhELambda1LocMax1[2] ; //! E vs lambda1 of selected cluster, 1 local maxima in cluster
281 TH2F * fhELambda0LocMax2[2] ; //! E vs lambda0 of selected cluster, 2 local maxima in cluster
282 TH2F * fhELambda1LocMax2[2] ; //! E vs lambda1 of selected cluster, 2 local maxima in cluster
283 TH2F * fhELambda0LocMaxN[2] ; //! E vs lambda0 of selected cluster, N>2 local maxima in cluster
284 TH2F * fhELambda1LocMaxN[2] ; //! E vs lambda1 of selected cluster, N>2 local maxima in cluster
287 TH2F * fhTimeENoCut; //! time of cluster vs E, no cut
288 TH2F * fhTimeESPD; //! time of cluster vs E, IsSPDPileUp
289 TH2F * fhTimeESPDMulti; //! time of cluster vs E, IsSPDPileUpMulti
290 TH2F * fhTimeNPileUpVertSPD; //! time of cluster vs n pile-up vertices from SPD
291 TH2F * fhTimeNPileUpVertTrack; //! time of cluster vs n pile-up vertices from Tracks
292 TH2F * fhTimeNPileUpVertContributors; //! time of cluster vs n pile-up vertex from SPD contributors
293 TH2F * fhTimePileUpMainVertexZDistance; //! time of cluster vs difference of z main vertex and pile-up vertex
294 TH2F * fhTimePileUpMainVertexZDiamond; //! time of cluster vs difference of z diamond and pile-up vertex
296 //Histograms settings
297 Int_t fHistoNPtSumBins; // Number of bins in PtSum histograms
298 Float_t fHistoPtSumMax; // PtSum maximum in histogram
299 Float_t fHistoPtSumMin; // PtSum minimum in histogram
300 Int_t fHistoNPtInConeBins; // Number of bins in PtInCone histogram
301 Float_t fHistoPtInConeMax; // PtInCone maximum in histogram
302 Float_t fHistoPtInConeMin; // PtInCone maximum in histogram
304 AliAnaParticleIsolation( const AliAnaParticleIsolation & iso) ; // cpy ctor
305 AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
307 ClassDef(AliAnaParticleIsolation,17)
311 #endif //ALIANAPARTICLEISOLATION_H