]>
Commit | Line | Data |
---|---|---|
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 | //_________________________________________________________________________ | |
7 | // Class that contains the algorithm for the analysis of particle - hadron correlations | |
8 | // Particle (for example direct gamma) must be found in a previous analysis | |
9 | //-- Author: Gustavo Conesa (INFN-LNF) | |
10 | ||
11 | // Modified by Yaxian Mao: | |
12 | // 1. add the UE subtraction for corrlation study | |
13 | // 2. change the correlation variable | |
14 | // 3. Only use leading particle(cluster/track) as trigger for correlation (2010/07/02) | |
15 | // 4. Make decay photon-hadron correlations where decay contribute pi0 mass (2010/09/09) | |
16 | // 5. fill the pout to extract kt at the end, also to study charge asymmetry(2010/10/06) | |
17 | // 6. Add the possibility for event selection analysis based on vertex and multiplicity bins (10/10/2010) | |
18 | // 7. change the way of delta phi cut for UE study due to memory issue (reduce histograms) | |
19 | // 8. Add the possibility to request the absolute leading particle at the near side or not, set trigger bins, general clean-up (08/2011) | |
20 | ||
045396c8 | 21 | // --- Analysis system --- |
745913ae | 22 | #include "AliAnaCaloTrackCorrBaseClass.h" |
045396c8 | 23 | class AliAODPWG4ParticleCorrelation ; |
24 | ||
745913ae | 25 | class AliAnaParticleHadronCorrelation : public AliAnaCaloTrackCorrBaseClass { |
045396c8 | 26 | |
27 | public: | |
c5693f62 | 28 | |
029dea5a | 29 | AliAnaParticleHadronCorrelation() ; // default ctor |
30 | virtual ~AliAnaParticleHadronCorrelation() ; // virtual dtor | |
045396c8 | 31 | |
32 | // General methods | |
907b38cd | 33 | |
045396c8 | 34 | TObjString * GetAnalysisCuts(); |
35 | ||
36 | TList * GetCreateOutputObjects(); | |
37 | ||
065234f0 | 38 | void Init(); |
39 | ||
045396c8 | 40 | void InitParameters(); |
31864468 | 41 | |
42 | void FillEventMixPool() ; | |
065234f0 | 43 | |
045396c8 | 44 | void MakeAnalysisFillHistograms() ; |
45 | ||
46 | void Print(const Option_t * opt) const; | |
47 | ||
48 | // Main analysis methods | |
49 | ||
907b38cd | 50 | Bool_t GetDecayPhotonMomentum(const AliAODPWG4Particle* trigger, TLorentzVector & mom1,TLorentzVector & mom2); |
51 | ||
2b65bd0e | 52 | Bool_t MakeChargedCorrelation (AliAODPWG4ParticleCorrelation * aodParticle, const TObjArray* pl) ; |
045396c8 | 53 | |
2b65bd0e | 54 | Bool_t MakeNeutralCorrelation (AliAODPWG4ParticleCorrelation * aodParticle, const TObjArray* pl) ; |
045396c8 | 55 | |
7ed3a248 | 56 | void MakeMCChargedCorrelation(AliAODPWG4ParticleCorrelation * aodParticle) ; |
045396c8 | 57 | |
7ed3a248 | 58 | void MakeChargedMixCorrelation(AliAODPWG4ParticleCorrelation *aodParticle) ; |
029dea5a | 59 | |
907b38cd | 60 | // Filling histogram methods |
61 | ||
22ad7981 | 62 | void FillChargedAngularCorrelationHistograms (Float_t ptAssoc, Float_t ptTrig, Int_t assocBin, |
63 | Float_t phiAssoc, Float_t phiTrig, Float_t & deltaPhi, | |
64 | Float_t etaAssoc, Float_t etaTrig, | |
65 | Bool_t decay, Float_t hmpidSignal, Int_t outTOF, | |
66 | Int_t nTracks, Int_t mcTag); | |
907b38cd | 67 | |
029dea5a | 68 | void FillChargedEventMixPool(); |
69 | ||
22ad7981 | 70 | Bool_t FillChargedMCCorrelationHistograms (Float_t mcAssocPt, Float_t mcAssocPhi, Float_t mcAssocEta, |
71 | Float_t mcTrigPt, Float_t mcTrigPhi, Float_t mcTrigEta ); | |
907b38cd | 72 | |
73 | ||
22ad7981 | 74 | void FillChargedMomentumImbalanceHistograms (Float_t ptTrig, Float_t ptAssoc, |
75 | Float_t xE, Float_t hbpXE, | |
76 | Float_t zT, Float_t hbpZT, | |
727a309a | 77 | Float_t pout, Float_t deltaPhi, Int_t nTracks, Int_t charge, |
22ad7981 | 78 | Int_t assocBin, Bool_t decay, |
79 | Int_t outTOF, Int_t mcTag ); | |
907b38cd | 80 | |
22ad7981 | 81 | void FillChargedUnderlyingEventHistograms (Float_t ptTrig, Float_t ptAssoc, |
82 | Float_t deltaPhi, Int_t nTracks, Int_t outTOF); | |
907b38cd | 83 | |
22ad7981 | 84 | void FillChargedUnderlyingEventSidesHistograms(Float_t ptTrig, Float_t ptAssoc, |
85 | Float_t deltaPhi); | |
907b38cd | 86 | |
22ad7981 | 87 | void FillDecayPhotonCorrelationHistograms (Float_t ptAssoc, Float_t phiAssoc, |
88 | TLorentzVector mom1, TLorentzVector mom2, | |
89 | Bool_t bChargedOrNeutral); | |
907b38cd | 90 | |
91 | ||
22ad7981 | 92 | void FillNeutralAngularCorrelationHistograms (Float_t ptAssoc, Float_t ptTrig, |
93 | Float_t phiAssoc, Float_t phiTrig, Float_t & deltaPhi, | |
94 | Float_t etaAssoc, Float_t etaTrig); | |
907b38cd | 95 | |
d07ffd54 | 96 | void FillNeutralEventMixPool(); |
97 | ||
2bb7ac98 | 98 | |
22ad7981 | 99 | void FillNeutralUnderlyingEventSidesHistograms(Float_t ptTrig, Float_t ptAssoc, |
100 | Float_t xE, Float_t hbpXE, | |
101 | Float_t zT, Float_t hbpZT, | |
102 | Float_t deltaPhi); | |
b1f720a7 | 103 | |
22ad7981 | 104 | Int_t GetMCTagHistogramIndex(Int_t tag); |
045396c8 | 105 | |
f3c47a7f | 106 | Bool_t IsTriggerTheEventLeadingParticle(); |
065234f0 | 107 | |
045396c8 | 108 | // Parameter setter and getter |
109 | ||
1f8591af | 110 | Float_t GetMinimumTriggerPt() const { return GetMinPt() ; } |
111 | Float_t GetMaximumTriggerPt() const { return GetMaxPt() ; } | |
112 | void SetTriggerPtRange(Float_t min, Float_t max) | |
113 | { SetMinPt(min), SetMaxPt(max) ; } | |
114 | ||
f7408d50 | 115 | |
907b38cd | 116 | Float_t GetMaximumAssociatedPt() const { return fMaxAssocPt ; } |
117 | Float_t GetMinimumAssociatedPt() const { return fMinAssocPt ; } | |
f7408d50 | 118 | void SetAssociatedPtRange(Float_t min, Float_t max) |
119 | { fMaxAssocPt = max ; fMinAssocPt = min ; } | |
120 | ||
907b38cd | 121 | Double_t GetDeltaPhiMaxCut() const { return fDeltaPhiMaxCut ; } |
122 | Double_t GetDeltaPhiMinCut() const { return fDeltaPhiMinCut ; } | |
f7408d50 | 123 | void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax) |
124 | { fDeltaPhiMaxCut = phimax ; fDeltaPhiMinCut = phimin ; } | |
66e64043 | 125 | |
f7408d50 | 126 | // Leading Hadron |
127 | Double_t GetLeadHadronPhiMaxCut() const { return fMaxLeadHadPhi ; } | |
128 | Double_t GetLeadHadronPhiMinCut() const { return fMinLeadHadPhi ; } | |
129 | void SetLeadHadronPhiCut(Float_t min, Float_t max) | |
130 | { fMaxLeadHadPhi = max ; fMinLeadHadPhi = min ; } | |
131 | ||
132 | Double_t GetLeadHadronPtMinCut() const { return fMinLeadHadPt ; } | |
133 | Double_t GetLeadHadronPtMaxCut() const { return fMaxLeadHadPt ; } | |
134 | void SetLeadHadronPtCut(Float_t min, Float_t max) | |
135 | { fMaxLeadHadPt = max ; fMinLeadHadPt = min ; } | |
045396c8 | 136 | |
f7408d50 | 137 | Bool_t IsLeadHadronCutOn() const { return fSelectLeadingHadronAngle ; } |
138 | void SwitchOnLeadHadronSelection() { fSelectLeadingHadronAngle = kTRUE ; } | |
139 | void SwitchOffLeadHadronSelection() { fSelectLeadingHadronAngle = kFALSE ; } | |
66e64043 | 140 | |
f7408d50 | 141 | // UE |
66e64043 | 142 | |
f7408d50 | 143 | Double_t GetUeDeltaPhiMaxCut() const { return fUeDeltaPhiMaxCut ; } |
144 | Double_t GetUeDeltaPhiMinCut() const { return fUeDeltaPhiMinCut ; } | |
66e64043 | 145 | |
045396c8 | 146 | void SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax) |
907b38cd | 147 | { fUeDeltaPhiMaxCut = uephimax ; fUeDeltaPhiMinCut = uephimin ; } |
045396c8 | 148 | |
907b38cd | 149 | Bool_t IsSeveralUEOn() const { return fMakeSeveralUE ; } |
150 | void SwitchOnSeveralUECalculation() { fMakeSeveralUE = kTRUE ; } | |
151 | void SwitchOffSeveralUECalculation() { fMakeSeveralUE = kFALSE ; } | |
045396c8 | 152 | |
153 | // Do trigger-neutral correlation | |
907b38cd | 154 | Bool_t DoNeutralCorr() const { return fNeutralCorr ; } |
155 | void SwitchOnNeutralCorr() { fNeutralCorr = kTRUE ; } | |
156 | void SwitchOffNeutralCorr() { fNeutralCorr = kFALSE ; } | |
045396c8 | 157 | |
158 | // Taking the absolute leading as the trigger or not | |
907b38cd | 159 | Bool_t DoAbsoluteLeading() const { return fMakeAbsoluteLeading ; } |
160 | void SwitchOnAbsoluteLeading() { fMakeAbsoluteLeading = kTRUE ; } | |
161 | void SwitchOffAbsoluteLeading() { fMakeAbsoluteLeading = kFALSE ; } | |
045396c8 | 162 | |
3f150b4b | 163 | // Taking the near side leading as the trigger or not |
907b38cd | 164 | Bool_t DoNearSideLeading() const { return fMakeNearSideLeading ; } |
165 | void SwitchOnNearSideLeading() { fMakeNearSideLeading = kTRUE ; } | |
166 | void SwitchOffNearSideLeading() { fMakeNearSideLeading = kFALSE ; } | |
3f150b4b | 167 | |
045396c8 | 168 | // Do decay-hadron correlation if it is pi0 trigger |
907b38cd | 169 | Bool_t IsPi0Trigger() const { return fPi0Trigger ; } |
170 | void SwitchOnPi0TriggerDecayCorr() { fPi0Trigger = kTRUE ; } | |
171 | void SwitchOffPi0TriggerDecayCorr() { fPi0Trigger = kFALSE ; } | |
172 | ||
173 | Bool_t IsDecayTrigger() const { return fDecayTrigger ; } | |
174 | void SwitchOnDecayTriggerDecayCorr() { fDecayTrigger = kTRUE ; } | |
175 | void SwitchOffDecayTriggerDecayCorr() { fDecayTrigger = kFALSE ; } | |
176 | ||
177 | Bool_t IsHMPIDCorrelation() const { return fHMPIDCorrelation ; } | |
178 | void SwitchOnHMPIDCorrelation() { fHMPIDCorrelation = kTRUE ; } | |
179 | void SwitchOffHMPIDCorrelation() { fHMPIDCorrelation = kFALSE ; } | |
045396c8 | 180 | |
907b38cd | 181 | void SwitchOnFillBradHistograms() { fFillBradHisto = kTRUE ; } |
182 | void SwitchOffFillBradHistograms() { fFillBradHisto = kFALSE ; } | |
183 | ||
184 | Bool_t OnlyIsolated() const { return fSelectIsolated ; } | |
185 | void SelectIsolated(Bool_t s) { fSelectIsolated = s ; } | |
045396c8 | 186 | |
907b38cd | 187 | void SetPi0AODBranchName(TString n) { fPi0AODBranchName = n ; } |
045396c8 | 188 | |
05d0d05d | 189 | void SetNAssocPtBins(Int_t n) ; |
190 | void SetAssocPtBinLimit(Int_t ibin, Float_t pt) ; | |
283f989c | 191 | |
192 | Bool_t IsMixStoredInReaderOn() const { return fUseMixStoredInReader ; } | |
193 | void SwitchOnUseMixStoredInReader() { fUseMixStoredInReader = kTRUE ; } | |
194 | void SwitchOffUseMixStoredInReader() { fUseMixStoredInReader = kFALSE; } | |
195 | ||
2bb7ac98 | 196 | void SwitchOnFillNeutralInMixedEvent() { fFillNeutralEventMixPool = kTRUE ; } |
197 | void SwitchOffFillNeutralInMixedEvent(){ fFillNeutralEventMixPool = kFALSE ; } | |
198 | ||
e4c0a0eb | 199 | void SetM02Cut(Float_t min=0, Float_t max=10) { fM02MinCut = min ; fM02MaxCut = max ; } |
200 | ||
2bb7ac98 | 201 | void SwitchOnCorrelationVzBin() { fCorrelVzBin = kTRUE ; } |
202 | void SwitchOffCorrelationVzBin() { fCorrelVzBin = kFALSE ; } | |
e4c0a0eb | 203 | |
b1f720a7 | 204 | void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; } |
560498d1 | 205 | void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; } |
2b65bd0e | 206 | |
207 | void SwitchOnFillTriggerAODWithReferences() { fFillAODWithReferences = kTRUE ; } | |
208 | void SwitchOffFillTriggerAODWithReferences(){ fFillAODWithReferences = kFALSE ; } | |
209 | ||
f3c47a7f | 210 | void SwitchOnCheckNeutralClustersForLeading() { fCheckLeadingWithNeutralClusters = kTRUE ; } |
211 | void SwitchOffCheckNeutralClustersForLeading(){ fCheckLeadingWithNeutralClusters = kFALSE ; } | |
b1f720a7 | 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 fMinTriggerPt ; // Minimum trigger hadron pt |
218 | Float_t fMaxAssocPt ; // Maximum associated hadron pt | |
219 | Float_t fMinAssocPt ; // Minimum associated hadron pt | |
045396c8 | 220 | Double_t fDeltaPhiMaxCut ; // Minimum Delta Phi Gamma-Hadron |
221 | Double_t fDeltaPhiMinCut ; // Maximum Delta Phi Gamma-Hadron | |
222 | Bool_t fSelectIsolated ; // Select only trigger particles isolated | |
223 | Bool_t fMakeSeveralUE ; // Do analysis for several underlying events contribution | |
224 | Double_t fUeDeltaPhiMaxCut ; // Minimum Delta Phi Gamma-Underlying Hadron | |
225 | Double_t fUeDeltaPhiMinCut ; // Maximum Delta Phi Gamma-Underlying Hadron | |
226 | TString fPi0AODBranchName; // Name of AOD branch with pi0, not trigger | |
227 | Bool_t fNeutralCorr ; // switch the analysis with neutral particles | |
228 | Bool_t fPi0Trigger ; // switch the analysis with decay photon from pi0 trigger | |
907b38cd | 229 | Bool_t fDecayTrigger ; // switch the analysis with decay photon from photon trigger |
3f150b4b | 230 | Bool_t fMakeAbsoluteLeading ; // requesting absolute leading triggers |
231 | Bool_t fMakeNearSideLeading ; // requesting near side leading (+-90º from trigger particle) triggers | |
045396c8 | 232 | Int_t fLeadingTriggerIndex ; // Store here per event the trigger index, to avoid too many loops |
907b38cd | 233 | Bool_t fHMPIDCorrelation ; // Correlate with particles on HMPID or its acceptance |
234 | Bool_t fFillBradHisto ; // DPhi histograms calculated differently | |
05d0d05d | 235 | Int_t fNAssocPtBins ; // Number of associated pT bins under study |
f85f4afb | 236 | Float_t fAssocPtBinLimit[20] ; // Associated pT under study |
7d306daf | 237 | Bool_t fCorrelVzBin ; // Fill one histogram per vz bin |
045396c8 | 238 | |
d07ffd54 | 239 | TList ** fListMixTrackEvents ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for tracks in stored events for mixing |
240 | TList ** fListMixCaloEvents ; //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for calo clusters in stored events for mixing | |
241 | ||
283f989c | 242 | Bool_t fUseMixStoredInReader; // Signal if in the current event the pool was filled |
2bb7ac98 | 243 | Bool_t fFillNeutralEventMixPool; // Add clusters to pool if requested |
283f989c | 244 | |
e4c0a0eb | 245 | Float_t fM02MaxCut ; // Study photon clusters with l0 smaller than cut |
246 | Float_t fM02MinCut ; // Study photon clusters with l0 larger than cut | |
247 | ||
b1f720a7 | 248 | Bool_t fFillPileUpHistograms; // Fill pile-up related histograms |
249 | ||
f7408d50 | 250 | Bool_t fSelectLeadingHadronAngle; // Select events with leading particle within a range |
251 | Float_t fMinLeadHadPhi; // Minimum angle between the trigger and leading hadron | |
252 | Float_t fMaxLeadHadPhi; // Maximum ange between the trigger and leading hadron | |
253 | Float_t fMinLeadHadPt; // Minimum pT of leading hadron | |
254 | Float_t fMaxLeadHadPt; // Maximum pT of leading hadron | |
255 | ||
045396c8 | 256 | //Histograms |
257 | ||
258 | //leading particles | |
26118504 | 259 | TH1F * fhPtTriggerInput; //! pT distribution of trigger particles before selection |
260 | TH1F * fhPtTriggerSSCut; //! pT distribution of trigger particles after shower shape selection | |
261 | TH1F * fhPtTriggerIsoCut; //! pT distribution of trigger particles after isolation cut selection | |
262 | TH1F * fhPtTriggerFidCut; //! pT distribution of trigger particles after fiducial selection | |
7d306daf | 263 | TH1F * fhPtLeading; //! pT distribution of leading particles |
cc944149 | 264 | TH1F * fhPtLeadingVtxBC0; //! pT distribution of leading particles |
17af6e24 | 265 | TH1F * fhPtLeadingPileUp[7]; //! pT distribution of leading particles |
7d306daf | 266 | TH2F * fhPtLeadingVzBin; //! pT distribution of leading particles vs vz bin |
d0b625bc | 267 | TH2F * fhPtLeadingBin; //! pT distribution of leading particles, vs mixing bin |
045396c8 | 268 | TH2F * fhPhiLeading; //! phi distribution vs pT of leading particles |
269 | TH2F * fhEtaLeading; //! eta distribution vs pT of leading particles | |
270 | ||
5a951426 | 271 | TH1F * fhPtLeadingMC[7]; //! pT distribution of leading particles, check the origin of the cluster : "Photon","Pi0","Pi0Decay","EtaDecay","OtherDecay","Electron","Hadron" |
764ab1f4 | 272 | |
68cd2624 | 273 | TH2F * fhPtLeadingCentrality; //! pT distribution of leading particles vs centrality |
274 | TH2F * fhPtLeadingEventPlane; //! pT distribution of leading particles vs centrality | |
275 | TH2F * fhLeadingEventPlaneCentrality; //! event plane vs centrality for leading particles | |
276 | ||
7d306daf | 277 | TH1F * fhPtLeadingMixed; //! pT distribution of leading particles, used in mixing |
278 | TH2F * fhPtLeadingMixedVzBin; //! pT distribution of leading particles, used in mixing, vs vz bin | |
d0b625bc | 279 | TH2F * fhPtLeadingMixedBin; //! pT distribution of leading particles vs mixing bin |
280 | TH2F * fhPhiLeadingMixed; //! phi distribution vs pT of leading particles, used in mixing | |
281 | TH2F * fhEtaLeadingMixed; //! eta distribution vs pT of leading particles, used in mixing | |
282 | ||
045396c8 | 283 | //trigger-charged histograms |
284 | TH2F * fhDeltaPhiDeltaEtaCharged ; //! differences of eta and phi between trigger and charged hadrons | |
285 | TH2F * fhPhiCharged ; //! Phi distribution of charged particles | |
286 | TH2F * fhEtaCharged ; //! Eta distribution of charged particles | |
287 | TH2F * fhDeltaPhiCharged ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT | |
288 | TH2F * fhDeltaEtaCharged ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT | |
289 | TH2F * fhDeltaPhiChargedPt ; //! Difference of charged particle phi and trigger particle phi as function of charged particle pT | |
290 | TH2F * fhDeltaPhiUeChargedPt ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of charged particle pT | |
9623bf34 | 291 | TH1F * fhUePart; //! UE particles distribution vs pt trig |
3f150b4b | 292 | TH2F * fhXECharged ; //! Trigger particle -charged hadron momentum imbalance histogram |
727a309a | 293 | TH2F * fhXECharged_Cone2 ; //! Trigger particle -charged hadron momentum imbalance histogram in cone2 (5pi/6-7pi/6) |
3f150b4b | 294 | TH2F * fhXEUeCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram |
295 | TH2F * fhXEPosCharged ; //! Trigger particle -positive charged hadron momentum imbalance histogram | |
296 | TH2F * fhXENegCharged ; //! Trigger particle -negative charged hadron momentum imbalance histogram | |
297 | TH2F * fhPtHbpXECharged ; //! Trigger particle -charged hadron momentum HBP histogram | |
727a309a | 298 | TH2F * fhPtHbpXECharged_Cone2 ; //! Trigger particle -charged hadron momentum HBP histogram in cone2 (5pi/6-7pi/6) |
3f150b4b | 299 | TH2F * fhPtHbpXEUeCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram |
300 | TH2F * fhZTCharged ; //! Trigger particle -charged hadron momentum imbalance histogram | |
301 | TH2F * fhZTUeCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
302 | TH2F * fhZTPosCharged ; //! Trigger particle -positive charged hadron momentum imbalance histogram | |
303 | TH2F * fhZTNegCharged ; //! Trigger particle -negative charged hadron momentum imbalance histogram | |
304 | TH2F * fhPtHbpZTCharged ; //! Trigger particle -charged hadron momentum HBP histogram | |
305 | TH2F * fhPtHbpZTUeCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram | |
2bb7ac98 | 306 | |
b1f720a7 | 307 | TH2F * fhXEChargedMC[7] ; //! 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) |
308 | TH2F * fhDeltaPhiChargedMC[7] ; //! 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) | |
309 | ||
310 | TH2F * fhDeltaPhiDeltaEtaChargedPtA3GeV;//! differences of eta and phi between trigger and charged hadrons, pTa > 3 GeV | |
311 | TH2F * fhDeltaPhiChargedPtA3GeV ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV | |
312 | TH2F * fhDeltaEtaChargedPtA3GeV ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT, pTa > 3 GeV | |
313 | ||
17af6e24 | 314 | // Events tagged as pileup by SDD,EMCal, or combination |
315 | TH2F * fhDeltaPhiChargedPileUp[7] ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT | |
316 | TH2F * fhDeltaEtaChargedPileUp[7] ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT | |
fedea415 | 317 | TH2F * fhDeltaPhiChargedPtA3GeVPileUp[7] ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV |
318 | TH2F * fhDeltaEtaChargedPtA3GeVPileUp[7] ; //! Difference of charged particle eta and trigger particle eta as function of trigger particle pT, pTa > 3 GeV | |
17af6e24 | 319 | TH2F * fhXEChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram |
320 | TH2F * fhXEUeChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram | |
321 | TH2F * fhZTChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram | |
322 | TH2F * fhZTUeChargedPileUp[7] ; //! Trigger particle -charged hadron momentum imbalance histogram | |
323 | TH2F * fhPtTrigChargedPileUp[7] ; //! trigger and correlated particl pt, to be used for mean value for kt | |
fedea415 | 324 | |
325 | TH2F * fhDeltaPhiChargedOtherBC ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT | |
326 | TH2F * fhDeltaPhiChargedPtA3GeVOtherBC ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV | |
327 | TH2F * fhXEChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram | |
328 | TH2F * fhXEUeChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram | |
329 | TH2F * fhZTChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram | |
330 | TH2F * fhZTUeChargedOtherBC ; //! Trigger particle -charged hadron momentum imbalance histogram | |
331 | TH2F * fhPtTrigChargedOtherBC ; //! trigger and correlated particl pt, to be used for mean value for kt | |
332 | ||
2a9171b5 | 333 | TH2F * fhDeltaPhiChargedBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT |
334 | TH2F * fhDeltaPhiChargedPtA3GeVBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV | |
335 | TH2F * fhXEChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
336 | TH2F * fhXEUeChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
337 | TH2F * fhZTChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
338 | TH2F * fhZTUeChargedBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
339 | TH2F * fhPtTrigChargedBC0 ; //! trigger and correlated particl pt, to be used for mean value for kt | |
340 | ||
cc944149 | 341 | TH2F * fhDeltaPhiChargedVtxBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT |
342 | TH2F * fhDeltaPhiChargedPtA3GeVVtxBC0 ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT, pTa > 3 GeV | |
343 | TH2F * fhXEChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
344 | TH2F * fhXEUeChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
345 | TH2F * fhZTChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
346 | TH2F * fhZTUeChargedVtxBC0 ; //! Trigger particle -charged hadron momentum imbalance histogram | |
347 | TH2F * fhPtTrigChargedVtxBC0 ; //! trigger and correlated particl pt, to be used for mean value for kt | |
348 | ||
045396c8 | 349 | //if several UE calculation is on, most useful for jet-jet events contribution |
350 | TH2F * fhDeltaPhiUeLeftCharged ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of charged particle pT | |
351 | TH2F * fhDeltaPhiUeRightCharged ; //! Difference of charged particle from underlying events phi and trigger particle phi | |
8cc41381 | 352 | TH2F * fhDeltaPhiUeLeftUpCharged; //! Difference of charged particle from underlying events phi and trigger particle phi |
353 | TH2F * fhDeltaPhiUeRightUpCharged; //! Difference of charged particle from underlying events phi and trigger particle phi | |
354 | TH2F * fhDeltaPhiUeLeftDownCharged; //! Difference of charged particle from underlying events phi and trigger particle phi | |
355 | TH2F * fhDeltaPhiUeRightDownCharged; //! Difference of charged particle from underlying events phi and trigger particle phi | |
3f150b4b | 356 | TH2F * fhXEUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram |
357 | TH2F * fhXEUeRightCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
b1f720a7 | 358 | TH2F * fhXEUeLeftUpCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram |
359 | TH2F * fhXEUeRightUpCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
360 | TH2F * fhXEUeLeftDownCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
361 | TH2F * fhXEUeRightDownCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
3f150b4b | 362 | TH2F * fhPtHbpXEUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram |
363 | TH2F * fhPtHbpXEUeRightCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram | |
364 | TH2F * fhZTUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
365 | TH2F * fhZTUeRightCharged ; //! Trigger particle -underlying charged hadron momentum imbalance histogram | |
366 | TH2F * fhPtHbpZTUeLeftCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram | |
367 | TH2F * fhPtHbpZTUeRightCharged ; //! Trigger particle -underlying charged hadron momentum HBP histogram | |
368 | ||
045396c8 | 369 | //for pout and kt extraction |
05d0d05d | 370 | TH2F * fhPtTrigPout ; //! Pout =associated pt*sin(delta phi) distribution vs trigger pt |
045396c8 | 371 | TH2F * fhPtTrigCharged ; //! trigger and correlated particl pt, to be used for mean value for kt |
372 | ||
373 | //if different multiplicity analysis asked | |
05d0d05d | 374 | TH2F ** fhTrigDeltaPhiCharged ; //![GetMultiBin()] differences of phi between trigger and charged hadrons |
375 | TH2F ** fhTrigDeltaEtaCharged ; //![GetMultiBin()] differences of eta between trigger and charged hadrons | |
907b38cd | 376 | TH2F ** fhTrigXECorr ; //![GetMultiBin()] Trigger particle -charged hadron momentum imbalance histogram |
377 | TH2F ** fhTrigXEUeCorr ; //![GetMultiBin()] Trigger particle -UE charged hadron momentum imbalance histogram | |
378 | TH2F ** fhTrigZTCorr ; //![GetMultiBin()] Trigger particle -charged hadron momentum imbalance histogram | |
379 | TH2F ** fhTrigZTUeCorr ; //![GetMultiBin()] Trigger particle -UE charged hadron momentum imbalance histogram | |
3f150b4b | 380 | |
05d0d05d | 381 | TH2F * fhAssocPtBkg; //! Trigger pT vs associated pT for background |
8cc41381 | 382 | TH2F ** fhDeltaPhiDeltaEtaAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function eta difference, for different associated bins |
383 | TH2F ** fhDeltaPhiAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt bins | |
384 | TH2F ** fhDeltaPhiAssocPtBinDEta08; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt bins for Delta eta > 0.8 | |
385 | TH2F ** fhDeltaPhiAssocPtBinDEta0 ; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt bins for Delta eta = 0 | |
386 | TH2F ** fhDeltaPhiAssocPtBinHMPID; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt bins, track with HMPID | |
2bb7ac98 | 387 | TH2F ** fhDeltaPhiAssocPtBinHMPIDAcc; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt bins, track with HMPIDAcc |
8cc41381 | 388 | TH2F ** fhDeltaPhiBradAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi Brad (?) for different associated pt bins |
66e64043 | 389 | TH2F * fhDeltaPhiBrad; //! Trigger pT vs dPhi Brad (?) for different associated pt bins |
8cc41381 | 390 | TH2F ** fhXEAssocPtBin ; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs xE for different associated pt bins |
391 | TH2F ** fhZTAssocPtBin ; //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs zT for different associated pt bins | |
3f150b4b | 392 | |
045396c8 | 393 | //trigger-neutral histograms |
394 | TH2F * fhDeltaPhiDeltaEtaNeutral ; //! differences of eta and phi between trigger and neutral hadrons (pi0) | |
395 | TH2F * fhPhiNeutral ; //! Phi distribution of neutral particles | |
396 | TH2F * fhEtaNeutral ; //! Eta distribution of neutral particles | |
397 | TH2F * fhDeltaPhiNeutral ; //! Difference of neutral particle phi and trigger particle phi as function of trigger particle pT | |
398 | TH2F * fhDeltaEtaNeutral ; //! Difference of neutral particle eta and trigger particle eta as function of trigger particle pT | |
399 | TH2F * fhDeltaPhiNeutralPt ; //! Difference of neutral particle phi and trigger particle phi as function of neutral particle particle pT | |
400 | TH2F * fhDeltaPhiUeNeutralPt ; //! Difference of neutral particle phi and trigger particle phi as function of neutral particle particle pT | |
3f150b4b | 401 | TH2F * fhXENeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram |
402 | TH2F * fhXEUeNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram | |
907b38cd | 403 | TH2F * fhPtHbpXENeutral ; //! Trigger particle - neutral particle momentum HBP histogram |
404 | TH2F * fhPtHbpXEUeNeutral ; //! Trigger particle - underlying neutral hadron momentum HBP histogram | |
3f150b4b | 405 | TH2F * fhZTNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram |
406 | TH2F * fhZTUeNeutral ; //! Trigger particle - neutral hadron momentum imbalance histogram | |
907b38cd | 407 | TH2F * fhPtHbpZTNeutral ; //! Trigger particle - neutral particle momentum HBP histogram |
408 | TH2F * fhPtHbpZTUeNeutral ; //! Trigger particle - underlying neutral hadron momentum HBP histogram | |
045396c8 | 409 | |
045396c8 | 410 | //if several UE calculation is on, most useful for jet-jet events contribution |
411 | TH2F * fhDeltaPhiUeLeftNeutral ; //! Difference of charged particle from underlying events phi and trigger particle phi as function of neutral particle pT | |
412 | TH2F * fhDeltaPhiUeRightNeutral ; //! Difference of charged particle from underlying events phi and trigger particle phi | |
3f150b4b | 413 | TH2F * fhXEUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram |
414 | TH2F * fhXEUeRightNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram | |
415 | TH2F * fhPtHbpXEUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram | |
416 | TH2F * fhPtHbpXEUeRightNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram | |
417 | TH2F * fhZTUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram | |
418 | TH2F * fhZTUeRightNeutral ; //! Trigger particle -underlying neutral hadron momentum imbalance histogram | |
419 | TH2F * fhPtHbpZTUeLeftNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram | |
420 | TH2F * fhPtHbpZTUeRightNeutral ; //! Trigger particle -underlying neutral hadron momentum HBP histogram | |
045396c8 | 421 | |
422 | //for decay photon trigger correlation | |
423 | TH2F * fhPtPi0DecayRatio ; //! for pi0 pt and ratio of decay photon pt | |
424 | TH2F * fhDeltaPhiDecayCharged ; //! Difference of charged particle phi and decay trigger | |
3f150b4b | 425 | TH2F * fhXEDecayCharged ; //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram |
907b38cd | 426 | TH2F * fhZTDecayCharged ; //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram |
427 | ||
045396c8 | 428 | TH2F * fhDeltaPhiDecayNeutral ; //! Difference of neutral particle phi and decay trigger |
3f150b4b | 429 | TH2F * fhXEDecayNeutral ; //! Trigger particle (decay from pi0)-neutral hadron momentum imbalance histogram |
430 | TH2F * fhZTDecayNeutral ; //! Trigger particle (decay from pi0)-neutral hadron momentum imbalance histogram | |
431 | ||
8cc41381 | 432 | TH2F ** fhDeltaPhiDecayChargedAssocPtBin;//![fNAssocPtBins*GetNZvertBin()] Tagged as decay Trigger pT vs dPhi for different associated pt bins |
433 | TH2F ** fhXEDecayChargedAssocPtBin ; //![fNAssocPtBins*GetNZvertBin()] Tagged as decay Trigger pT vs xE for different associated pt bins | |
434 | TH2F ** fhZTDecayChargedAssocPtBin ; //![fNAssocPtBins*GetNZvertBin()] Tagged as decay Trigger pT vs xE for different associated pt bins | |
907b38cd | 435 | |
045396c8 | 436 | //if the data is MC, fill MC information |
437 | TH2F * fh2phiLeadingParticle; //! #phi resolution for triggers | |
9623bf34 | 438 | TH1F * fhMCPtLeading; //! MC pure pT distribution of leading particles |
08f220c8 | 439 | TH2F * fhMCPhiLeading; //! MC pure Phi distribution of leading particles |
440 | TH2F * fhMCEtaLeading; //! MC pure Eta distribution of leading particles | |
907b38cd | 441 | TH2F * fhMCEtaCharged; //! MC pure particles charged primary pt vs eta (both associated) |
442 | TH2F * fhMCPhiCharged; //! MC pure particles charged primary pt vs phi (both associated) | |
443 | TH2F * fhMCDeltaEtaCharged; //! MC pure particles charged trigger primary pt vs delta eta (associated-trigger) | |
444 | TH2F * fhMCDeltaPhiCharged; //! MC pure particles charged trigger primary pt vs delta phi (associated-trigger) | |
445 | TH2F * fhMCDeltaPhiDeltaEtaCharged; //! MC pure particles charged associated primary pt vs delta phi (associated-trigger), in away side | |
446 | TH2F * fhMCDeltaPhiChargedPt; //! MC pure particles charged delta phi vs delta eta (associated-trigger) | |
447 | TH2F * fhMCPtXECharged; //! MC pure particles charged trigger primary pt vs xE | |
9623bf34 | 448 | TH2F * fhMCPtXEUeCharged; //! MC pure particles charged trigger primary pt vs xE (underlying event) |
c78eef61 | 449 | TH2F * fhMCPtXEUeLeftCharged; //! MC pure particles charged trigger primary pt vs xE (underlying event,left cone) |
450 | TH2F * fhMCPtXEUeRightCharged; //! MC pure particles charged trigger primary pt vs xE (underlying event,right cone) | |
907b38cd | 451 | TH2F * fhMCPtHbpXECharged; //! MC pure particles charged trigger primary pt vs ln(1/xE) |
9623bf34 | 452 | TH2F * fhMCPtHbpXEUeCharged; //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event) |
c78eef61 | 453 | TH2F * fhMCPtHbpXEUeLeftCharged; //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event, left cone) |
454 | TH2F * fhMCPtHbpXEUeRightCharged; //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event, right cone) | |
9623bf34 | 455 | TH1F * fhMCUePart; //! MC pure UE particles distribution vs pt trig |
907b38cd | 456 | TH2F * fhMCPtZTCharged; //! MC pure particles charged trigger primary pt vs zT |
c78eef61 | 457 | TH2F * fhMCPtZTUeCharged; //! MC pure particles charged trigger primary pt vs zT (underlying event) |
458 | TH2F * fhMCPtZTUeLeftCharged; //! MC pure particles charged trigger primary pt vs zT (underlying event, left cone) | |
459 | TH2F * fhMCPtZTUeRightCharged; //! MC pure particles charged trigger primary pt vs zT (underlying event, right cone) | |
907b38cd | 460 | TH2F * fhMCPtHbpZTCharged; //! MC pure particles charged trigger primary pt vs ln(1/zT) |
c78eef61 | 461 | TH2F * fhMCPtHbpZTUeCharged; //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event) |
462 | TH2F * fhMCPtHbpZTUeLeftCharged; //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event, left cone) | |
463 | TH2F * fhMCPtHbpZTUeRightCharged; //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event, right cone) | |
907b38cd | 464 | TH2F * fhMCPtTrigPout ; //! MC pure particles charged trigger primary pt vs pOut |
465 | TH2F * fhMCPtAssocDeltaPhi ; //! MC pure particles charged associated primary pt vs delta phi (associated-trigger) | |
045396c8 | 466 | |
029dea5a | 467 | // Mixing |
468 | TH1I * fhNEventsTrigger; //! number of analyzed triggered events | |
2e876d85 | 469 | TH1F * fhNtracksMB; //! total number of tracks in MB events |
2bb7ac98 | 470 | TH1F * fhNclustersMB; //! total number of clusters in MB events |
029dea5a | 471 | TH2F * fhMixDeltaPhiCharged ; //! Difference of charged particle phi and trigger particle phi as function of trigger particle pT |
472 | TH2F * fhMixDeltaPhiDeltaEtaCharged ; //! Difference of charged particle phi and trigger particle phi as function eta difference | |
08f220c8 | 473 | TH2F * fhMixXECharged; //! xE for mixed event |
18838011 | 474 | TH2F * fhMixXEUeCharged; //! xE for mixed event in Ue region |
08f220c8 | 475 | TH2F * fhMixHbpXECharged; //! ln(1/xE) for mixed event |
8cc41381 | 476 | TH2F ** fhMixDeltaPhiChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function of trigger particle pT, for different associated bins |
477 | 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 | |
478 | 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 | |
479 | TH2F ** fhMixDeltaPhiDeltaEtaChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle phi as function eta difference, for different associated bins | |
029dea5a | 480 | |
2e876d85 | 481 | TH1I * fhEventBin; //! Number of real events in a particular bin (cen,vz,rp) |
482 | TH1I * fhEventMixBin; //! Number of mixed events in a particular bin (cen,vz,rp) | |
483 | ||
3f150b4b | 484 | AliAnaParticleHadronCorrelation( const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor |
485 | AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ; // cpy assignment | |
045396c8 | 486 | |
2b65bd0e | 487 | ClassDef(AliAnaParticleHadronCorrelation,32) |
045396c8 | 488 | } ; |
489 | ||
490 | ||
491 | #endif //ALIANAPARTICLEHADRONCORRELATION_H | |
492 | ||
493 | ||
494 |