change name of variable, conflicts with enum with detector types
[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaParticleHadronCorrelation.h
CommitLineData
045396c8 1#ifndef ALIANAPARTICLEHADRONCORRELATION_H
2#define ALIANAPARTICLEHADRONCORRELATION_H
3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
045396c8 5
6//_________________________________________________________________________
55c20a99 7// Class that contains the algorithm for the analysis of
8// particle - hadron correlations
9// Particle (for example direct gamma) must be found in a previous analysis
10//
b6afb6b1 11//-- Author: Gustavo Conesa (LNF-INFN) (LPSC-IN2P3-CNRS)
12// Yaxian Mao (LPSC-IN2P3-CNRS) and (CNWU) first usable implementation.
13// Xiangrong Zhu (CNWU), implementtion of own mixing.
55c20a99 14//
b6afb6b1 15
045396c8 16// --- Analysis system ---
55c20a99 17
745913ae 18#include "AliAnaCaloTrackCorrBaseClass.h"
045396c8 19class AliAODPWG4ParticleCorrelation ;
20
745913ae 21class AliAnaParticleHadronCorrelation : public AliAnaCaloTrackCorrBaseClass {
045396c8 22
23 public:
c5693f62 24
029dea5a 25 AliAnaParticleHadronCorrelation() ; // default ctor
26 virtual ~AliAnaParticleHadronCorrelation() ; // virtual dtor
045396c8 27
28 // General methods
907b38cd 29
045396c8 30 TObjString * GetAnalysisCuts();
31
32 TList * GetCreateOutputObjects();
33
065234f0 34 void Init();
35
045396c8 36 void InitParameters();
31864468 37
38 void FillEventMixPool() ;
065234f0 39
045396c8 40 void MakeAnalysisFillHistograms() ;
41
42 void Print(const Option_t * opt) const;
43
44 // Main analysis methods
45
6b013448 46 Bool_t FindLeadingOppositeHadronInWindow(AliAODPWG4ParticleCorrelation * particle);
907b38cd 47
55bf1615 48 Bool_t GetDecayPhotonMomentum (Int_t indexPhoton1, Int_t indexPhoton2, Int_t idetector);
045396c8 49
5267de40 50 void MakeChargedCorrelation (AliAODPWG4ParticleCorrelation * particle) ;
045396c8 51
bf01409a 52 void MakeNeutralCorrelation (AliAODPWG4ParticleCorrelation * particle) ;
045396c8 53
4a296f4a 54 void MakeMCChargedCorrelation (Int_t triggerMCLable, Int_t histoIndex, Bool_t lostDecayPair) ;
6b013448 55
56 void MakeChargedMixCorrelation(AliAODPWG4ParticleCorrelation * particle) ;
029dea5a 57
907b38cd 58 // Filling histogram methods
59
22ad7981 60 void FillChargedAngularCorrelationHistograms (Float_t ptAssoc, Float_t ptTrig, Int_t assocBin,
6ed1ad3e 61 Float_t phiAssoc, Float_t phiTrig, Float_t deltaPhi,
22ad7981 62 Float_t etaAssoc, Float_t etaTrig,
e7d2373d 63 Int_t decayTag, Float_t hmpidSignal, Int_t outTOF,
d4a8ff9a 64 Int_t cenbin, Int_t mcTag);
907b38cd 65
029dea5a 66 void FillChargedEventMixPool();
67
6ed1ad3e 68 Bool_t FillChargedMCCorrelationHistograms (Float_t mcAssocPt, Float_t mcAssocPhi, Float_t mcAssocEta,
4a296f4a 69 Float_t mcTrigPt, Float_t mcTrigPhi, Float_t mcTrigEta,
70 Int_t histoIndex, Bool_t lostDecayPair);
907b38cd 71
22ad7981 72 void FillChargedMomentumImbalanceHistograms (Float_t ptTrig, Float_t ptAssoc,
e7d2373d 73 Float_t deltaPhi, Int_t cenbin, Int_t charge,
74 Int_t assocBin, Int_t decayTag,
75 Int_t outTOF, Int_t mcTag );
907b38cd 76
22ad7981 77 void FillChargedUnderlyingEventHistograms (Float_t ptTrig, Float_t ptAssoc,
d4a8ff9a 78 Float_t deltaPhi, Int_t cenbin, Int_t outTOF);
907b38cd 79
22ad7981 80 void FillChargedUnderlyingEventSidesHistograms(Float_t ptTrig, Float_t ptAssoc,
81 Float_t deltaPhi);
907b38cd 82
349eb23d 83 void FillDecayPhotonCorrelationHistograms (Float_t ptAssoc, Float_t phiAssoc, Bool_t bChargedOrNeutral);
907b38cd 84
d07ffd54 85 void FillNeutralEventMixPool();
86
2bb7ac98 87
22ad7981 88 void FillNeutralUnderlyingEventSidesHistograms(Float_t ptTrig, Float_t ptAssoc,
22ad7981 89 Float_t zT, Float_t hbpZT,
90 Float_t deltaPhi);
b1f720a7 91
22ad7981 92 Int_t GetMCTagHistogramIndex(Int_t tag);
4a296f4a 93 static const Int_t fgkNmcTypes = 8;
0ff0cd12 94
f3c47a7f 95 Bool_t IsTriggerTheEventLeadingParticle();
065234f0 96
045396c8 97 // Parameter setter and getter
98
1f8591af 99 Float_t GetMinimumTriggerPt() const { return GetMinPt() ; }
100 Float_t GetMaximumTriggerPt() const { return GetMaxPt() ; }
101 void SetTriggerPtRange(Float_t min, Float_t max)
102 { SetMinPt(min), SetMaxPt(max) ; }
103
f7408d50 104
907b38cd 105 Float_t GetMaximumAssociatedPt() const { return fMaxAssocPt ; }
106 Float_t GetMinimumAssociatedPt() const { return fMinAssocPt ; }
f7408d50 107 void SetAssociatedPtRange(Float_t min, Float_t max)
108 { fMaxAssocPt = max ; fMinAssocPt = min ; }
109
907b38cd 110 Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; }
111 Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; }
f7408d50 112 void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
113 { fDeltaPhiMaxCut = phimax ; fDeltaPhiMinCut = phimin ; }
66e64043 114
f7408d50 115 // Leading Hadron
116 Double_t GetLeadHadronPhiMaxCut() const { return fMaxLeadHadPhi ; }
117 Double_t GetLeadHadronPhiMinCut() const { return fMinLeadHadPhi ; }
118 void SetLeadHadronPhiCut(Float_t min, Float_t max)
119 { fMaxLeadHadPhi = max ; fMinLeadHadPhi = min ; }
120
121 Double_t GetLeadHadronPtMinCut() const { return fMinLeadHadPt ; }
122 Double_t GetLeadHadronPtMaxCut() const { return fMaxLeadHadPt ; }
123 void SetLeadHadronPtCut(Float_t min, Float_t max)
124 { fMaxLeadHadPt = max ; fMinLeadHadPt = min ; }
045396c8 125
98604e3a 126 Bool_t IsLeadHadronCutOn() const { return fSelectLeadingHadronAngle ; }
f7408d50 127 void SwitchOnLeadHadronSelection() { fSelectLeadingHadronAngle = kTRUE ; }
128 void SwitchOffLeadHadronSelection() { fSelectLeadingHadronAngle = kFALSE ; }
66e64043 129
98604e3a 130 void SwitchOnFillLeadHadronHistograms() { fFillLeadHadOppositeHisto = kTRUE ; }
131 void SwitchOffFillLeadHadronHistograms(){ fFillLeadHadOppositeHisto = kFALSE ; }
132
f7408d50 133 // UE
66e64043 134
f7408d50 135 Double_t GetUeDeltaPhiMaxCut() const { return fUeDeltaPhiMaxCut ; }
136 Double_t GetUeDeltaPhiMinCut() const { return fUeDeltaPhiMinCut ; }
66e64043 137
045396c8 138 void SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax)
907b38cd 139 { fUeDeltaPhiMaxCut = uephimax ; fUeDeltaPhiMinCut = uephimin ; }
045396c8 140
907b38cd 141 Bool_t IsSeveralUEOn() const { return fMakeSeveralUE ; }
142 void SwitchOnSeveralUECalculation() { fMakeSeveralUE = kTRUE ; }
143 void SwitchOffSeveralUECalculation() { fMakeSeveralUE = kFALSE ; }
045396c8 144
145 // Do trigger-neutral correlation
907b38cd 146 Bool_t DoNeutralCorr() const { return fNeutralCorr ; }
147 void SwitchOnNeutralCorr() { fNeutralCorr = kTRUE ; }
148 void SwitchOffNeutralCorr() { fNeutralCorr = kFALSE ; }
045396c8 149
150 // Taking the absolute leading as the trigger or not
907b38cd 151 Bool_t DoAbsoluteLeading() const { return fMakeAbsoluteLeading ; }
152 void SwitchOnAbsoluteLeading() { fMakeAbsoluteLeading = kTRUE ; }
153 void SwitchOffAbsoluteLeading() { fMakeAbsoluteLeading = kFALSE ; }
045396c8 154
3f150b4b 155 // Taking the near side leading as the trigger or not
907b38cd 156 Bool_t DoNearSideLeading() const { return fMakeNearSideLeading ; }
157 void SwitchOnNearSideLeading() { fMakeNearSideLeading = kTRUE ; }
158 void SwitchOffNearSideLeading() { fMakeNearSideLeading = kFALSE ; }
3f150b4b 159
045396c8 160 // Do decay-hadron correlation if it is pi0 trigger
907b38cd 161 Bool_t IsPi0Trigger() const { return fPi0Trigger ; }
162 void SwitchOnPi0TriggerDecayCorr() { fPi0Trigger = kTRUE ; }
e7d2373d 163 void SwitchOffPi0TriggerDecayCorr() { fPi0Trigger = kFALSE ; }
164
907b38cd 165 Bool_t IsDecayTrigger() const { return fDecayTrigger ; }
166 void SwitchOnDecayTriggerDecayCorr() { fDecayTrigger = kTRUE ; }
167 void SwitchOffDecayTriggerDecayCorr() { fDecayTrigger = kFALSE ; }
e7d2373d 168 void SetNDecayBits(Int_t n) { fNDecayBits = n ; }
169 void SetDecayBits(Int_t i, UInt_t bit) { if(i < 4) fDecayBits[i] = bit ; }
170
907b38cd 171 Bool_t IsHMPIDCorrelation() const { return fHMPIDCorrelation ; }
172 void SwitchOnHMPIDCorrelation() { fHMPIDCorrelation = kTRUE ; }
173 void SwitchOffHMPIDCorrelation() { fHMPIDCorrelation = kFALSE ; }
045396c8 174
907b38cd 175 void SwitchOnFillBradHistograms() { fFillBradHisto = kTRUE ; }
176 void SwitchOffFillBradHistograms() { fFillBradHisto = kFALSE ; }
177
178 Bool_t OnlyIsolated() const { return fSelectIsolated ; }
179 void SelectIsolated(Bool_t s) { fSelectIsolated = s ; }
045396c8 180
907b38cd 181 void SetPi0AODBranchName(TString n) { fPi0AODBranchName = n ; }
045396c8 182
05d0d05d 183 void SetNAssocPtBins(Int_t n) ;
184 void SetAssocPtBinLimit(Int_t ibin, Float_t pt) ;
283f989c 185
186 Bool_t IsMixStoredInReaderOn() const { return fUseMixStoredInReader ; }
187 void SwitchOnUseMixStoredInReader() { fUseMixStoredInReader = kTRUE ; }
188 void SwitchOffUseMixStoredInReader() { fUseMixStoredInReader = kFALSE; }
189
2bb7ac98 190 void SwitchOnFillNeutralInMixedEvent() { fFillNeutralEventMixPool = kTRUE ; }
191 void SwitchOffFillNeutralInMixedEvent(){ fFillNeutralEventMixPool = kFALSE ; }
192
e4c0a0eb 193 void SetM02Cut(Float_t min=0, Float_t max=10) { fM02MinCut = min ; fM02MaxCut = max ; }
194
2bb7ac98 195 void SwitchOnCorrelationVzBin() { fCorrelVzBin = kTRUE ; }
196 void SwitchOffCorrelationVzBin() { fCorrelVzBin = kFALSE ; }
e4c0a0eb 197
b6afb6b1 198 void SwitchOnFillTriggerAODWithReferences() { fFillAODWithReferences = kTRUE ; }
199 void SwitchOffFillTriggerAODWithReferences() { fFillAODWithReferences = kFALSE ; }
2b65bd0e 200
f3c47a7f 201 void SwitchOnCheckNeutralClustersForLeading() { fCheckLeadingWithNeutralClusters = kTRUE ; }
202 void SwitchOffCheckNeutralClustersForLeading(){ fCheckLeadingWithNeutralClusters = kFALSE ; }
b1f720a7 203
a22f3f34 204 void SwitchOnFillEtaGapHistograms() { fFillEtaGapsHisto = kTRUE ; }
205 void SwitchOffFillEtaGapHistograms() { fFillEtaGapsHisto = kFALSE ; }
206
207 void SwitchOnFillPtImbalancePerPtABinHistograms() { fFillMomImbalancePtAssocBinsHisto = kTRUE ; }
208 void SwitchOffFillPtImbalancePerPtABinHistograms() { fFillMomImbalancePtAssocBinsHisto = kFALSE ; }
209
0ff0cd12 210 void SetMCGenType(Int_t min = 0, Int_t max = 6) { if(min >= 0 && min < fgkNmcTypes) fMCGenTypeMin = min ;
211 if(max >= 0 && max < fgkNmcTypes) fMCGenTypeMax = max ; }
be2df971 212
045396c8 213 private:
2b65bd0e 214
215 Bool_t fFillAODWithReferences; // Add to the trigger particle AOD the reference to the tracks or neutrals in correlation.
f3c47a7f 216 Bool_t fCheckLeadingWithNeutralClusters;// Compare the trigger candidate to Leading pT with the clusters pT, by default only charged
66e64043 217 Float_t fMaxAssocPt ; // Maximum associated hadron pt
218 Float_t fMinAssocPt ; // Minimum associated hadron pt
045396c8 219 Double_t fDeltaPhiMaxCut ; // Minimum Delta Phi Gamma-Hadron
220 Double_t fDeltaPhiMinCut ; // Maximum Delta Phi Gamma-Hadron
221 Bool_t fSelectIsolated ; // Select only trigger particles isolated
222 Bool_t fMakeSeveralUE ; // Do analysis for several underlying events contribution
223 Double_t fUeDeltaPhiMaxCut ; // Minimum Delta Phi Gamma-Underlying Hadron
224 Double_t fUeDeltaPhiMinCut ; // Maximum Delta Phi Gamma-Underlying Hadron
225 TString fPi0AODBranchName; // Name of AOD branch with pi0, not trigger
226 Bool_t fNeutralCorr ; // switch the analysis with neutral particles
227 Bool_t fPi0Trigger ; // switch the analysis with decay photon from pi0 trigger
907b38cd 228 Bool_t fDecayTrigger ; // switch the analysis with decay photon from photon trigger
e7d2373d 229 Int_t fNDecayBits ; // in case of study of decay triggers, select the decay bit
230 UInt_t fDecayBits[4] ; // in case of study of decay triggers, select the decay bit
3f150b4b 231 Bool_t fMakeAbsoluteLeading ; // requesting absolute leading triggers
232 Bool_t fMakeNearSideLeading ; // requesting near side leading (+-90ยบ from trigger particle) triggers
045396c8 233 Int_t fLeadingTriggerIndex ; // Store here per event the trigger index, to avoid too many loops
907b38cd 234 Bool_t fHMPIDCorrelation ; // Correlate with particles on HMPID or its acceptance
235 Bool_t fFillBradHisto ; // DPhi histograms calculated differently
05d0d05d 236 Int_t fNAssocPtBins ; // Number of associated pT bins under study
f85f4afb 237 Float_t fAssocPtBinLimit[20] ; // Associated pT under study
7d306daf 238 Bool_t fCorrelVzBin ; // Fill one histogram per vz bin
045396c8 239
d07ffd54 240 TList ** fListMixTrackEvents ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for tracks in stored events for mixing
241 TList ** fListMixCaloEvents ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for calo clusters in stored events for mixing
242
283f989c 243 Bool_t fUseMixStoredInReader; // Signal if in the current event the pool was filled
2bb7ac98 244 Bool_t fFillNeutralEventMixPool; // Add clusters to pool if requested
283f989c 245
e4c0a0eb 246 Float_t fM02MaxCut ; // Study photon clusters with l0 smaller than cut
247 Float_t fM02MinCut ; // Study photon clusters with l0 larger than cut
248
f7408d50 249 Bool_t fSelectLeadingHadronAngle; // Select events with leading particle within a range
98604e3a 250 Bool_t fFillLeadHadOppositeHisto; // Fill histograms for leading hadrons in opposite side of trigger
251
f7408d50 252 Float_t fMinLeadHadPhi; // Minimum angle between the trigger and leading hadron
253 Float_t fMaxLeadHadPhi; // Maximum ange between the trigger and leading hadron
254 Float_t fMinLeadHadPt; // Minimum pT of leading hadron
255 Float_t fMaxLeadHadPt; // Maximum pT of leading hadron
256
a22f3f34 257 Bool_t fFillEtaGapsHisto; // Fill azimuthal correlation histograms in 2 eta gaps, |eta|>0.8 and |eta|<0.01
258 Bool_t fFillMomImbalancePtAssocBinsHisto; // momentum imbalance histograms in bins of pT associated
259
0ff0cd12 260 Int_t fMCGenTypeMin; // Of the fgkNmcTypes possible types, select those between fMCGenTypeMin and fMCGenTypeMax
261 Int_t fMCGenTypeMax; // Of the fgkNmcTypes possible types, select those between fMCGenTypeMin and fMCGenTypeMax
be2df971 262
349eb23d 263 TVector3 fTrackVector; //! track momentum vector
264 TLorentzVector fMomentum; //! trigger momentum
265 TLorentzVector fDecayMom1; //! decay particle momentum
266 TLorentzVector fDecayMom2; //! decay particle momentum
267
045396c8 268 //Histograms
269
55c20a99 270 //trigger particles
26118504 271 TH1F * fhPtTriggerInput; //! pT distribution of trigger particles before selection
272 TH1F * fhPtTriggerSSCut; //! pT distribution of trigger particles after shower shape selection
273 TH1F * fhPtTriggerIsoCut; //! pT distribution of trigger particles after isolation cut selection
274 TH1F * fhPtTriggerFidCut; //! pT distribution of trigger particles after fiducial selection
55c20a99 275 TH1F * fhPtTrigger; //! pT distribution of trigger particles
276 TH1F * fhPtTriggerVtxBC0; //! pT distribution of trigger particles
277 TH1F * fhPtTriggerPileUp[7]; //! pT distribution of trigger particles
278 TH2F * fhPtTriggerVzBin; //! pT distribution of trigger particles vs vz bin
279 TH2F * fhPtTriggerBin; //! pT distribution of trigger particles, vs mixing bin
280 TH2F * fhPhiTrigger; //! phi distribution vs pT of trigger particles
281 TH2F * fhEtaTrigger; //! eta distribution vs pT of trigger particles
282
0ff0cd12 283 TH1F * fhPtTriggerMC[fgkNmcTypes]; //! pT distribution of trigger particles, check the origin of the cluster : "Photon","Pi0","Pi0Decay","EtaDecay","OtherDecay","Electron","Hadron"
55c20a99 284
48c5da66 285 TH1F * fhPtDecayTrigger[4]; //! pT distribution of trigger particles, tagged as decay
0ff0cd12 286 TH1F * fhPtDecayTriggerMC[4][fgkNmcTypes];//! pT distribution of trigger particles, tagged as decay, check the origin of the cluster
48c5da66 287
55c20a99 288 TH2F * fhPtTriggerCentrality; //! pT distribution of trigger particles vs centrality
289 TH2F * fhPtTriggerEventPlane; //! pT distribution of trigger particles vs centrality
290 TH2F * fhTriggerEventPlaneCentrality; //! event plane vs centrality for trigger particles
291
292 TH1F * fhPtTriggerMixed; //! pT distribution of trigger particles, used in mixing
293 TH2F * fhPtTriggerMixedVzBin; //! pT distribution of trigger particles, used in mixing, vs vz bin
294 TH2F * fhPtTriggerMixedBin; //! pT distribution of trigger particles vs mixing bin
295 TH2F * fhPhiTriggerMixed; //! phi distribution vs pT of trigger particles, used in mixing
296 TH2F * fhEtaTriggerMixed; //! eta distribution vs pT of trigger particles, used in mixing
764ab1f4 297
55c20a99 298 // Leading hadron in the opposite side of the trigger
a22f3f34 299 TH2F * fhPtLeadingOppositeHadron; //! pT trigger : pT distribution of leading hadron oposite to trigger
300 TH2F * fhPtDiffPhiLeadingOppositeHadron; //! pT trigger : difference phi distribution of leading hadron oposite and trigger
301 TH2F * fhPtDiffEtaLeadingOppositeHadron; //! pT trigger: difference eta distribution of leading hadron oposite and trigger
98604e3a 302 TH1F * fhPtNoLeadingOppositeHadron; //! pT trigger for events without opposite hadrons
303 TH2F * fhEtaPhiNoLeadingOppositeHadron; //! location of trigger when no hadron is found on the opposite side
d0b625bc 304
045396c8 305 //trigger-charged histograms
306 TH2F * fhDeltaPhiDeltaEtaCharged ; //! differences of eta and phi between trigger and charged hadrons
307 TH2F * fhPhiCharged ; //! Phi distribution of charged particles
308 TH2F * fhEtaCharged ; //! Eta distribution of charged particles
309 TH2F * fhDeltaPhiCharged ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
310 TH2F * fhDeltaEtaCharged ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT
311 TH2F * fhDeltaPhiChargedPt ; //! Difference of charged particle phi and trigger particle phi as function of charged particle pT
312 TH2F * fhDeltaPhiUeChargedPt ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of charged particle pT
9623bf34 313 TH1F * fhUePart; //! UE particles distribution vs pt trig
3f150b4b 314 TH2F * fhXECharged ; //! Trigger particle -charged hadron momentum imbalance histogram
727a309a 315 TH2F * fhXECharged_Cone2 ; //! Trigger particle -charged hadron momentum imbalance histogram in cone2 (5pi/6-7pi/6)
3f150b4b 316 TH2F * fhXEUeCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
317 TH2F * fhXEPosCharged ; //! Trigger particle -positive charged hadron momentum imbalance histogram
318 TH2F * fhXENegCharged ; //! Trigger particle -negative charged hadron momentum imbalance histogram
319 TH2F * fhPtHbpXECharged ; //! Trigger particle -charged hadron momentum HBP histogram
727a309a 320 TH2F * fhPtHbpXECharged_Cone2 ; //! Trigger particle -charged hadron momentum HBP histogram in cone2 (5pi/6-7pi/6)
3f150b4b 321 TH2F * fhPtHbpXEUeCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram
322 TH2F * fhZTCharged ; //! Trigger particle -charged hadron momentum imbalance histogram
323 TH2F * fhZTUeCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
324 TH2F * fhZTPosCharged ; //! Trigger particle -positive charged hadron momentum imbalance histogram
325 TH2F * fhZTNegCharged ; //! Trigger particle -negative charged hadron momentum imbalance histogram
326 TH2F * fhPtHbpZTCharged ; //! Trigger particle -charged hadron momentum HBP histogram
327 TH2F * fhPtHbpZTUeCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram
2bb7ac98 328
0ff0cd12 329 TH2F * fhXEChargedMC[fgkNmcTypes] ; //! Trigger particle -charged hadron momentum imbalance histogram, check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0), hadron, rest of photons (prompt, FSR, ISR)
330 TH2F * fhDeltaPhiChargedMC[fgkNmcTypes];//! Trigger particle -charged hadron delta phi histogram, check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0), hadron, rest of photons (prompt, FSR, ISR)
b1f720a7 331
332 TH2F * fhDeltaPhiDeltaEtaChargedPtA3GeV;//! differences of eta and phi between trigger and charged hadrons, pTa > 3 GeV
333 TH2F * fhDeltaPhiChargedPtA3GeV ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV
334 TH2F * fhDeltaEtaChargedPtA3GeV ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT, pTa > 3 GeV
335
17af6e24 336 // Events tagged as pileup by SDD,EMCal, or combination
337 TH2F * fhDeltaPhiChargedPileUp[7] ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
338 TH2F * fhDeltaEtaChargedPileUp[7] ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT
fedea415 339 TH2F * fhDeltaPhiChargedPtA3GeVPileUp[7] ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV
340 TH2F * fhDeltaEtaChargedPtA3GeVPileUp[7] ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT, pTa > 3 GeV
17af6e24 341 TH2F * fhXEChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram
342 TH2F * fhXEUeChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram
343 TH2F * fhZTChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram
344 TH2F * fhZTUeChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram
345 TH2F * fhPtTrigChargedPileUp[7] ; //! trigger and correlated particl pt, to be used for mean value for kt
fedea415 346
347 TH2F * fhDeltaPhiChargedOtherBC ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
348 TH2F * fhDeltaPhiChargedPtA3GeVOtherBC ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV
349 TH2F * fhXEChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram
350 TH2F * fhXEUeChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram
351 TH2F * fhZTChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram
352 TH2F * fhZTUeChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram
353 TH2F * fhPtTrigChargedOtherBC ; //! trigger and correlated particl pt, to be used for mean value for kt
354
2a9171b5 355 TH2F * fhDeltaPhiChargedBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
356 TH2F * fhDeltaPhiChargedPtA3GeVBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV
357 TH2F * fhXEChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
358 TH2F * fhXEUeChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
359 TH2F * fhZTChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
360 TH2F * fhZTUeChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
361 TH2F * fhPtTrigChargedBC0 ; //! trigger and correlated particl pt, to be used for mean value for kt
362
cc944149 363 TH2F * fhDeltaPhiChargedVtxBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
364 TH2F * fhDeltaPhiChargedPtA3GeVVtxBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV
365 TH2F * fhXEChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
366 TH2F * fhXEUeChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
367 TH2F * fhZTChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
368 TH2F * fhZTUeChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram
369 TH2F * fhPtTrigChargedVtxBC0 ; //! trigger and correlated particl pt, to be used for mean value for kt
370
045396c8 371 //if several UE calculation is on, most useful for jet-jet events contribution
372 TH2F * fhDeltaPhiUeLeftCharged ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of charged particle pT
6ed1ad3e 373 TH2F * fhDeltaPhiUeLeftUpCharged; //! Difference of charged particle from underlying events phi and trigger particle phi
8cc41381 374 TH2F * fhDeltaPhiUeRightUpCharged; //! Difference of charged particle from underlying events phi and trigger particle phi
375 TH2F * fhDeltaPhiUeLeftDownCharged; //! Difference of charged particle from underlying events phi and trigger particle phi
376 TH2F * fhDeltaPhiUeRightDownCharged; //! Difference of charged particle from underlying events phi and trigger particle phi
3f150b4b 377 TH2F * fhXEUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
6ed1ad3e 378 TH2F * fhXEUeLeftUpCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
b1f720a7 379 TH2F * fhXEUeRightUpCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
380 TH2F * fhXEUeLeftDownCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
381 TH2F * fhXEUeRightDownCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
3f150b4b 382 TH2F * fhPtHbpXEUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram
6ed1ad3e 383 TH2F * fhZTUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram
384 TH2F * fhPtHbpZTUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram
3f150b4b 385
045396c8 386 //for pout and kt extraction
05d0d05d 387 TH2F * fhPtTrigPout ; //! Pout =associated pt*sin(delta phi) distribution vs trigger pt
045396c8 388 TH2F * fhPtTrigCharged ; //! trigger and correlated particl pt, to be used for mean value for kt
389
390 //if different multiplicity analysis asked
d4a8ff9a 391 TH2F ** fhDeltaPhiChargedMult ; //![GetNCentrBin()] differences of phi between trigger and charged hadrons: multiplicity bin
392 TH2F ** fhDeltaEtaChargedMult ; //![GetNCentrBin()] differences of eta between trigger and charged hadrons: multiplicity bin
393 TH2F ** fhXEMult ; //![GetNCentrBin()] Trigger particle -charged hadron momentum imbalance histogram: multiplicity bin
394 TH2F ** fhXEUeMult ; //![GetNCentrBin()] Trigger particle -UE charged hadron momentum imbalance histogram: multiplicity bin
395 TH2F ** fhZTMult ; //![GetNCentrBin()] Trigger particle -charged hadron momentum imbalance histogram: multiplicity bin
396 TH2F ** fhZTUeMult ; //![GetNCentrBin()] Trigger particle -UE charged hadron momentum imbalance histogram: multiplicity bin
3f150b4b 397
05d0d05d 398 TH2F * fhAssocPtBkg; //! Trigger pT vs associated pT for background
8cc41381 399 TH2F ** fhDeltaPhiDeltaEtaAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function eta difference, for different associated bins
a22f3f34 400 TH2F ** fhDeltaPhiAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins
401 TH2F ** fhDeltaPhiAssocPtBinDEta08; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins for Delta eta > 0.8
402 TH2F ** fhDeltaPhiAssocPtBinDEta0 ; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins for Delta eta = 0
403 TH2F ** fhDeltaPhiAssocPtBinHMPID; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins, track with HMPID
404 TH2F ** fhDeltaPhiAssocPtBinHMPIDAcc; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins, track with HMPIDAcc
8cc41381 405 TH2F ** fhDeltaPhiBradAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi Brad (?) for different associated pt bins
66e64043 406 TH2F * fhDeltaPhiBrad; //! Trigger pT vs dPhi Brad (?) for different associated pt bins
a22f3f34 407 TH2F ** fhXEAssocPtBin ; //![fNAssocPtBins] Trigger pT vs xE for different associated pt bins
408 TH2F ** fhZTAssocPtBin ; //![fNAssocPtBins] Trigger pT vs zT for different associated pt bins
409 TH2F ** fhXEVZ ; //![GetNZvertBin()] Trigger pT vs xE for different vz bins
410 TH2F ** fhZTVZ ; //![GetNZvertBin()] Trigger pT vs zT for different vz bins
3f150b4b 411
045396c8 412 //trigger-neutral histograms
413 TH2F * fhDeltaPhiDeltaEtaNeutral ; //! differences of eta and phi between trigger and neutral hadrons (pi0)
414 TH2F * fhPhiNeutral ; //! Phi distribution of neutral particles
415 TH2F * fhEtaNeutral ; //! Eta distribution of neutral particles
416 TH2F * fhDeltaPhiNeutral ; //! Difference of neutral particle phi and trigger particle phi as function of trigger particle pT
417 TH2F * fhDeltaEtaNeutral ; //! Difference of neutral particle eta and trigger particle eta as function of trigger particle pT
418 TH2F * fhDeltaPhiNeutralPt ; //! Difference of neutral particle phi and trigger particle phi as function of neutral particle particle pT
419 TH2F * fhDeltaPhiUeNeutralPt ; //! Difference of neutral particle phi and trigger particle phi as function of neutral particle particle pT
3f150b4b 420 TH2F * fhXENeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram
421 TH2F * fhXEUeNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram
907b38cd 422 TH2F * fhPtHbpXENeutral ; //! Trigger particle - neutral particle momentum HBP histogram
423 TH2F * fhPtHbpXEUeNeutral ; //! Trigger particle - underlying neutral hadron momentum HBP histogram
3f150b4b 424 TH2F * fhZTNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram
425 TH2F * fhZTUeNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram
907b38cd 426 TH2F * fhPtHbpZTNeutral ; //! Trigger particle - neutral particle momentum HBP histogram
427 TH2F * fhPtHbpZTUeNeutral ; //! Trigger particle - underlying neutral hadron momentum HBP histogram
045396c8 428
6ed1ad3e 429 // If several UE calculation is on,
045396c8 430 TH2F * fhDeltaPhiUeLeftNeutral ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of neutral particle pT
6ed1ad3e 431 TH2F * fhXEUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram
432 TH2F * fhPtHbpXEUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram
433 TH2F * fhZTUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram
434 TH2F * fhPtHbpZTUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram
045396c8 435
e7d2373d 436 // Pi0/Eta trigger correlation, recover input photons
437 TH2F * fhPtPi0DecayRatio ; //! for pi0 trigger pt and ratio of decay photon pt
438 TH2F * fhDeltaPhiPi0DecayCharged ; //! Difference of charged particle phi and decay photon from pi0/eta trigger
439 TH2F * fhXEPi0DecayCharged ; //! Trigger particle (decay from pi0/eta trigger)-charged hadron momentum imbalance histogram
440 TH2F * fhZTPi0DecayCharged ; //! Trigger particle (decay from pi0/eta trigger)-charged hadron momentum imbalance histogram
907b38cd 441
e7d2373d 442 TH2F * fhDeltaPhiPi0DecayNeutral ; //! Difference of neutral particle phi and decay photon from pi0/eta trigger
443 TH2F * fhXEPi0DecayNeutral ; //! Trigger particle (decay from pi0/eta trigger)-neutral hadron momentum imbalance histogram
444 TH2F * fhZTPi0DecayNeutral ; //! Trigger particle (decay from pi0/eta trigger)-neutral hadron momentum imbalance histogram
3f150b4b 445
e7d2373d 446 // Decay photon trigger correlation
447 TH2F * fhDeltaPhiDecayCharged[4] ; //! Difference of charged particle phi and photon decay trigger
448 TH2F * fhXEDecayCharged[4] ; //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram
449 TH2F * fhZTDecayCharged[4] ; //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram
450 TH2F ** fhDeltaPhiDecayChargedAssocPtBin;//![fNAssocPtBins*GetNZvertBin()] Tagged as decay (fDecayBits[0]) Trigger pT vs dPhi for different associated pt bins
907b38cd 451
be2df971 452 // If the data is MC, correlation with generated particles
453 // check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0),
454 // hadron, rest of photons (prompt, FSR, ISR)
0ff0cd12 455 TH1F * fhMCPtTrigger[fgkNmcTypes]; //! MC pure pT distribution of trigger particles
456 TH2F * fhMCPhiTrigger[fgkNmcTypes]; //! MC pure Phi distribution of trigger particles
457 TH2F * fhMCEtaTrigger[fgkNmcTypes]; //! MC pure Eta distribution of trigger particles
458 TH1F * fhMCPtTriggerNotLeading[fgkNmcTypes]; //! MC pure pT distribution of trigger not leading particles
459 TH2F * fhMCPhiTriggerNotLeading[fgkNmcTypes]; //! MC pure Phi distribution of trigger not leading particles
460 TH2F * fhMCEtaTriggerNotLeading[fgkNmcTypes]; //! MC pure Eta distribution of trigger not leading particles
461 TH2F * fhMCEtaCharged[fgkNmcTypes]; //! MC pure particles charged primary pt vs eta (both associated)
462 TH2F * fhMCPhiCharged[fgkNmcTypes]; //! MC pure particles charged primary pt vs phi (both associated)
463 TH2F * fhMCDeltaEtaCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs delta eta (associated-trigger)
464 TH2F * fhMCDeltaPhiCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs delta phi (associated-trigger)
465 TH2F * fhMCDeltaPhiDeltaEtaCharged[fgkNmcTypes]; //! MC pure particles charged associated primary pt vs delta phi (associated-trigger), in away side
466 TH2F * fhMCDeltaPhiChargedPt[fgkNmcTypes]; //! MC pure particles charged delta phi vs delta eta (associated-trigger)
467 TH2F * fhMCPtXECharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs xE
468 TH2F * fhMCPtXEUeCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs xE (underlying event)
469 TH2F * fhMCPtXEUeLeftCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs xE (underlying event,left cone)
470 TH2F * fhMCPtHbpXECharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/xE)
471 TH2F * fhMCPtHbpXEUeCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event)
472 TH2F * fhMCPtHbpXEUeLeftCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event, left cone)
473 TH1F * fhMCUePart[fgkNmcTypes]; //! MC pure UE particles distribution vs pt trig
474 TH2F * fhMCPtZTCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs zT
475 TH2F * fhMCPtZTUeCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs zT (underlying event)
476 TH2F * fhMCPtZTUeLeftCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs zT (underlying event, left cone)
477 TH2F * fhMCPtHbpZTCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/zT)
478 TH2F * fhMCPtHbpZTUeCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event)
479 TH2F * fhMCPtHbpZTUeLeftCharged[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event, left cone)
480 TH2F * fhMCPtTrigPout[fgkNmcTypes]; //! MC pure particles charged trigger primary pt vs pOut
481 TH2F * fhMCPtAssocDeltaPhi[fgkNmcTypes]; //! MC pure particles charged associated primary pt vs delta phi (associated-trigger)
045396c8 482
029dea5a 483 // Mixing
484 TH1I * fhNEventsTrigger; //! number of analyzed triggered events
7aaa6c97 485 TH2F * fhNtracksMB; //! total number of tracks in MB events
486 TH2F * fhNclustersMB; //! total number of clusters in MB events
029dea5a 487 TH2F * fhMixDeltaPhiCharged ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT
488 TH2F * fhMixDeltaPhiDeltaEtaCharged ; //! Difference of charged particle phi and trigger particle phi as function eta difference
08f220c8 489 TH2F * fhMixXECharged; //! xE for mixed event
18838011 490 TH2F * fhMixXEUeCharged; //! xE for mixed event in Ue region
08f220c8 491 TH2F * fhMixHbpXECharged; //! ln(1/xE) for mixed event
8cc41381 492 TH2F ** fhMixDeltaPhiChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function of trigger particle pT, for different associated bins
493 TH2F ** fhMixDeltaPhiChargedAssocPtBinDEta08; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function of trigger particle pT, for different associated bins, delta eta > 0.8
494 TH2F ** fhMixDeltaPhiChargedAssocPtBinDEta0; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function of trigger particle pT, for different associated bins, delta eta = 0
495 TH2F ** fhMixDeltaPhiDeltaEtaChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function eta difference, for different associated bins
029dea5a 496
efb540af 497 TH1I * fhEventBin; //! Number of triggers in a particular event bin (cen,vz,rp)
498 TH1I * fhEventMixBin; //! Number of triggers mixed in a particular bin (cen,vz,rp)
499 TH1I * fhEventMBBin; //! Number of MB events in a particular bin (cen,vz,rp)
2e876d85 500
3f150b4b 501 AliAnaParticleHadronCorrelation( const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor
502 AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ; // cpy assignment
045396c8 503
e7d2373d 504 ClassDef(AliAnaParticleHadronCorrelation,35)
045396c8 505} ;
506
507
508#endif //ALIANAPARTICLEHADRONCORRELATION_H
509
510
511