]>
Commit | Line | Data |
---|---|---|
57189f04 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | ||
17 | //////////////////////////////////////////////////////////////////////// | |
18 | // // | |
19 | // Task for Heavy Flavour Electron-Hadron DeltaPhi Correlation // | |
20 | // Non-Photonic Electron identified with Invariant mass // | |
21 | // analysis methos in function SelectPhotonicElectron // | |
22 | // DeltaPhi calculated in function ElectronHadCorrel // | |
23 | // // | |
24 | // Author: Deepa Thomas (Utrecht University) // | |
25 | // // | |
26 | //////////////////////////////////////////////////////////////////////// | |
27 | ||
28 | #include "TChain.h" | |
29 | #include "TTree.h" | |
30 | #include "TH2F.h" | |
31 | #include "TMath.h" | |
32 | #include "TCanvas.h" | |
33 | #include "THnSparse.h" | |
34 | #include "TLorentzVector.h" | |
35 | #include "TString.h" | |
36 | #include "TFile.h" | |
37 | ||
38 | #include "AliAnalysisTask.h" | |
39 | #include "AliAnalysisManager.h" | |
40 | ||
41 | #include "AliESDEvent.h" | |
42 | #include "AliESDHandler.h" | |
43 | #include "AliAODEvent.h" | |
44 | #include "AliAODHandler.h" | |
45 | ||
46 | #include "AliAnalysisTaskElecHadronCorrel.h" | |
47 | #include "TGeoGlobalMagField.h" | |
48 | #include "AliLog.h" | |
49 | #include "AliAnalysisTaskSE.h" | |
50 | #include "TRefArray.h" | |
51 | #include "TVector.h" | |
987053ce | 52 | |
57189f04 | 53 | #include "AliESDInputHandler.h" |
987053ce | 54 | #include "AliAODInputHandler.h" |
57189f04 | 55 | #include "AliESDpid.h" |
987053ce | 56 | #include "AliAODPid.h" |
57189f04 | 57 | #include "AliESDtrackCuts.h" |
58 | #include "AliPhysicsSelection.h" | |
987053ce | 59 | #include "AliCentralitySelectionTask.h" |
57189f04 | 60 | #include "AliESDCaloCluster.h" |
61 | #include "AliAODCaloCluster.h" | |
520a29a8 | 62 | #include "AliESDCaloTrigger.h" |
57189f04 | 63 | #include "AliEMCALRecoUtils.h" |
64 | #include "AliEMCALGeometry.h" | |
65 | #include "AliGeomManager.h" | |
66 | #include "stdio.h" | |
67 | #include "TGeoManager.h" | |
68 | #include "iostream" | |
69 | #include "fstream" | |
3c56855b | 70 | |
3db57d71 | 71 | #include "AliEventPoolManager.h" |
3c56855b | 72 | |
520a29a8 | 73 | #include "AliCentrality.h" |
57189f04 | 74 | #include "AliMagF.h" |
75 | ||
76 | #include "AliKFParticle.h" | |
77 | #include "AliKFVertex.h" | |
20791315 | 78 | #include "AliSelectNonHFE.h" |
57189f04 | 79 | |
80 | #include "AliPID.h" | |
81 | #include "AliPIDResponse.h" | |
82 | #include "AliHFEcontainer.h" | |
83 | #include "AliHFEcuts.h" | |
84 | #include "AliHFEpid.h" | |
85 | #include "AliHFEpidBase.h" | |
86 | #include "AliHFEpidQAmanager.h" | |
87 | #include "AliHFEtools.h" | |
88 | #include "AliCFContainer.h" | |
89 | #include "AliCFManager.h" | |
987053ce | 90 | #include "AliVEvent.h" |
91 | #include "AliStack.h" | |
92 | #include "AliMCEvent.h" | |
93 | #include "TProfile.h" | |
94 | #include "AliESDVZERO.h" | |
95 | #include "AliAODVZERO.h" | |
96 | #include "TVector3.h" | |
97 | #include "TRandom2.h" | |
57189f04 | 98 | |
a9123929 | 99 | ClassImp(AliAnalysisTaskElecHadronCorrel) |
3db57d71 | 100 | ClassImp(AliehDPhiBasicParticle) |
a9123929 | 101 | //________________________________________________________________________ |
520a29a8 | 102 | AliAnalysisTaskElecHadronCorrel::AliAnalysisTaskElecHadronCorrel(const char *name) |
57189f04 | 103 | : AliAnalysisTaskSE(name) |
987053ce | 104 | ,fVevent(0) |
57189f04 | 105 | ,fESD(0) |
987053ce | 106 | ,fAOD(0) |
57189f04 | 107 | ,fGeom(0) |
20791315 | 108 | ,fpidResponse(0) |
57189f04 | 109 | ,fOutputList(0) |
520a29a8 | 110 | ,fTrackCuts1(new AliESDtrackCuts) |
111 | ,fTrackCuts2(new AliESDtrackCuts) | |
57189f04 | 112 | ,fCuts(0) |
113 | ,fIdentifiedAsOutInz(kFALSE) | |
114 | ,fPassTheEventCut(kFALSE) | |
115 | ,fRejectKinkMother(kFALSE) | |
116 | ,fVz(0.0) | |
117 | ,fCFM(0) | |
118 | ,fPID(0) | |
119 | ,fPIDqa(0) | |
57189f04 | 120 | ,fInvmassCut(0.01) |
3c56855b | 121 | ,fCentrality(0) |
122 | ,fCentralityMin(0) | |
123 | ,fCentralityMax(0) | |
124 | ,fkCentralityMethod(0) | |
a9123929 | 125 | ,fTPCnsigEleMin(-2) |
126 | ,fTPCnsigEleMax(2) | |
127 | ,fTPCnsigHadMin(-10) | |
128 | ,fTPCnsigHadMax(-3.5) | |
129 | ,fM02CutMin(0.03) | |
130 | ,fM02CutMax(0.5) | |
131 | ,fM20CutMin(0.03) | |
132 | ,fM20CutMax(0.3) | |
133 | ,fDispCutMin(0) | |
134 | ,fDispCutMax(1) | |
135 | ,fEovPMin(0.8) | |
136 | ,fEovPMax(1.2) | |
137 | ,fTriggerCentral(kTRUE) | |
6ec70f8c | 138 | ,fTriggerMB(kTRUE) |
28ea601a | 139 | ,fTPCNClsHad(80) |
140 | ,fAssoEleITSref(kTRUE) | |
ca49cfe0 | 141 | ,fAssoElecTPCNCls(80) |
20791315 | 142 | ,fNonHFE(0) |
3db57d71 | 143 | ,fPoolMgr(0x0) |
a9123929 | 144 | ,fNoEvents(0) |
145 | // ,fTrkpt(0) | |
146 | ,fTrkEovPAft(0) | |
147 | // ,fTrkEovPBefHad(0) | |
148 | // ,fdEdxBef(0) | |
987053ce | 149 | ,fSemiIncElecDphi(0) |
487ae960 | 150 | ,fSemiIncElecDphi1(0) |
151 | ,fSemiIncElecDphi2(0) | |
152 | ,fSemiIncElecDphi3(0) | |
153 | ,fSemiIncElecDphi4(0) | |
57189f04 | 154 | ,fPhotElecDphi(0) |
487ae960 | 155 | ,fPhotElecDphi1(0) |
156 | ,fPhotElecDphi2(0) | |
157 | ,fPhotElecDphi3(0) | |
158 | ,fPhotElecDphi4(0) | |
57189f04 | 159 | ,fInclusiveElecDphi(0) |
487ae960 | 160 | ,fInclusiveElecDphi1(0) |
161 | ,fInclusiveElecDphi2(0) | |
162 | ,fInclusiveElecDphi3(0) | |
163 | ,fInclusiveElecDphi4(0) | |
33369c14 | 164 | ,fInclusiveElecDphiEtaFS(0) |
165 | ,fInclusiveElecDphiEtaFS1(0) | |
166 | ,fInclusiveElecDphiEtaFS2(0) | |
167 | ,fInclusiveElecDphiEtaFS3(0) | |
168 | ,fInclusiveElecDphiEtaFS4(0) | |
57189f04 | 169 | ,fDphiULSMassLow(0) |
487ae960 | 170 | ,fDphiULSMassLow1(0) |
171 | ,fDphiULSMassLow2(0) | |
172 | ,fDphiULSMassLow3(0) | |
173 | ,fDphiULSMassLow4(0) | |
57189f04 | 174 | ,fDphiLSMassLow(0) |
487ae960 | 175 | ,fDphiLSMassLow1(0) |
176 | ,fDphiLSMassLow2(0) | |
177 | ,fDphiLSMassLow3(0) | |
178 | ,fDphiLSMassLow4(0) | |
520a29a8 | 179 | ,fDphiULSMassLowNoPartner(0) |
487ae960 | 180 | ,fDphiULSMassLowNoPartner1(0) |
181 | ,fDphiULSMassLowNoPartner2(0) | |
182 | ,fDphiULSMassLowNoPartner3(0) | |
183 | ,fDphiULSMassLowNoPartner4(0) | |
520a29a8 | 184 | ,fDphiLSMassLowNoPartner(0) |
487ae960 | 185 | ,fDphiLSMassLowNoPartner1(0) |
186 | ,fDphiLSMassLowNoPartner2(0) | |
187 | ,fDphiLSMassLowNoPartner3(0) | |
188 | ,fDphiLSMassLowNoPartner4(0) | |
57189f04 | 189 | ,fPhotoElecPt(0) |
190 | ,fSemiInclElecPt(0) | |
520a29a8 | 191 | ,fInclusiveElecPt(0) |
192 | ,fULSElecPt(0) | |
193 | ,fLSElecPt(0) | |
20791315 | 194 | ,fDCAMetPhotElecDphi(0) |
195 | ,fDCAMetPhotElecDphi1(0) | |
196 | ,fDCAMetPhotElecDphi2(0) | |
197 | ,fDCAMetPhotElecDphi3(0) | |
198 | ,fDCAMetPhotElecDphi4(0) | |
199 | ,fDCAMetDphiULSMassLow(0) | |
200 | ,fDCAMetDphiULSMassLow1(0) | |
201 | ,fDCAMetDphiULSMassLow2(0) | |
202 | ,fDCAMetDphiULSMassLow3(0) | |
203 | ,fDCAMetDphiULSMassLow4(0) | |
204 | ,fDCAMetDphiLSMassLow(0) | |
205 | ,fDCAMetDphiLSMassLow1(0) | |
206 | ,fDCAMetDphiLSMassLow2(0) | |
207 | ,fDCAMetDphiLSMassLow3(0) | |
208 | ,fDCAMetDphiLSMassLow4(0) | |
209 | ,fDCAMetDphiULSMassLowNoPartner(0) | |
210 | ,fDCAMetDphiULSMassLowNoPartner1(0) | |
211 | ,fDCAMetDphiULSMassLowNoPartner2(0) | |
212 | ,fDCAMetDphiULSMassLowNoPartner3(0) | |
213 | ,fDCAMetDphiULSMassLowNoPartner4(0) | |
214 | ,fDCAMetDphiLSMassLowNoPartner(0) | |
215 | ,fDCAMetDphiLSMassLowNoPartner1(0) | |
216 | ,fDCAMetDphiLSMassLowNoPartner2(0) | |
217 | ,fDCAMetDphiLSMassLowNoPartner3(0) | |
218 | ,fDCAMetDphiLSMassLowNoPartner4(0) | |
219 | ,fDCAMetPhotoElecPt(0) | |
220 | ,fDCAMetULSElecPt(0) | |
221 | ,fDCAMetLSElecPt(0) | |
bf697dda | 222 | ,fSemiIncElecDphiEta1(0) |
223 | ,fSemiIncElecDphiEta11(0) | |
224 | ,fSemiIncElecDphiEta12(0) | |
225 | ,fSemiIncElecDphiEta13(0) | |
226 | ,fSemiIncElecDphiEta14(0) | |
227 | ,fPhotElecDphiEta1(0) | |
228 | ,fPhotElecDphiEta11(0) | |
229 | ,fPhotElecDphiEta12(0) | |
230 | ,fPhotElecDphiEta13(0) | |
231 | ,fPhotElecDphiEta14(0) | |
232 | ,fInclusiveElecDphiEta1(0) | |
233 | ,fInclusiveElecDphiEta11(0) | |
234 | ,fInclusiveElecDphiEta12(0) | |
235 | ,fInclusiveElecDphiEta13(0) | |
236 | ,fInclusiveElecDphiEta14(0) | |
237 | ,fDphiULSMassLowEta1(0) | |
238 | ,fDphiULSMassLowEta11(0) | |
239 | ,fDphiULSMassLowEta12(0) | |
240 | ,fDphiULSMassLowEta13(0) | |
241 | ,fDphiULSMassLowEta14(0) | |
242 | ,fDphiLSMassLowEta1(0) | |
243 | ,fDphiLSMassLowEta11(0) | |
244 | ,fDphiLSMassLowEta12(0) | |
245 | ,fDphiLSMassLowEta13(0) | |
246 | ,fDphiLSMassLowEta14(0) | |
247 | ,fDphiULSMassLowNoPartnerEta1(0) | |
248 | ,fDphiULSMassLowNoPartnerEta11(0) | |
249 | ,fDphiULSMassLowNoPartnerEta12(0) | |
250 | ,fDphiULSMassLowNoPartnerEta13(0) | |
251 | ,fDphiULSMassLowNoPartnerEta14(0) | |
252 | ,fDphiLSMassLowNoPartnerEta1(0) | |
253 | ,fDphiLSMassLowNoPartnerEta11(0) | |
254 | ,fDphiLSMassLowNoPartnerEta12(0) | |
255 | ,fDphiLSMassLowNoPartnerEta13(0) | |
256 | ,fDphiLSMassLowNoPartnerEta14(0) | |
257 | ,fSemiIncElecDphiEta2(0) | |
258 | ,fSemiIncElecDphiEta21(0) | |
259 | ,fSemiIncElecDphiEta22(0) | |
260 | ,fSemiIncElecDphiEta23(0) | |
261 | ,fSemiIncElecDphiEta24(0) | |
262 | ,fPhotElecDphiEta2(0) | |
263 | ,fPhotElecDphiEta21(0) | |
264 | ,fPhotElecDphiEta22(0) | |
265 | ,fPhotElecDphiEta23(0) | |
266 | ,fPhotElecDphiEta24(0) | |
267 | ,fInclusiveElecDphiEta2(0) | |
268 | ,fInclusiveElecDphiEta21(0) | |
269 | ,fInclusiveElecDphiEta22(0) | |
270 | ,fInclusiveElecDphiEta23(0) | |
271 | ,fInclusiveElecDphiEta24(0) | |
272 | ,fDphiULSMassLowEta2(0) | |
273 | ,fDphiULSMassLowEta21(0) | |
274 | ,fDphiULSMassLowEta22(0) | |
275 | ,fDphiULSMassLowEta23(0) | |
276 | ,fDphiULSMassLowEta24(0) | |
277 | ,fDphiLSMassLowEta2(0) | |
278 | ,fDphiLSMassLowEta21(0) | |
279 | ,fDphiLSMassLowEta22(0) | |
280 | ,fDphiLSMassLowEta23(0) | |
281 | ,fDphiLSMassLowEta24(0) | |
282 | ,fDphiULSMassLowNoPartnerEta2(0) | |
283 | ,fDphiULSMassLowNoPartnerEta21(0) | |
284 | ,fDphiULSMassLowNoPartnerEta22(0) | |
285 | ,fDphiULSMassLowNoPartnerEta23(0) | |
286 | ,fDphiULSMassLowNoPartnerEta24(0) | |
287 | ,fDphiLSMassLowNoPartnerEta2(0) | |
288 | ,fDphiLSMassLowNoPartnerEta21(0) | |
289 | ,fDphiLSMassLowNoPartnerEta22(0) | |
290 | ,fDphiLSMassLowNoPartnerEta23(0) | |
291 | ,fDphiLSMassLowNoPartnerEta24(0) | |
292 | // ,fTrackPtBefTrkCuts(0) | |
293 | // ,fTrackPtAftTrkCuts(0) | |
7e02de85 | 294 | ,fTPCnsigma(0) |
bf697dda | 295 | // ,fNCellv1(0) |
296 | // ,fClsEv1(0) | |
297 | // ,fNClusv1(0) | |
298 | ,fInvmassLS1(0) | |
299 | // ,fInvmassLS2(0) | |
300 | // ,fInvmassLS3(0) | |
301 | // ,fInvmassLS4(0) | |
302 | // ,fInvmassLS5(0) | |
303 | ,fInvmassULS1(0) | |
20791315 | 304 | ,fDCAMetInvmassLS1(0) |
305 | ,fDCAMetInvmassULS1(0) | |
bf697dda | 306 | // ,fInvmassULS2(0) |
307 | // ,fInvmassULS3(0) | |
308 | // ,fInvmassULS4(0) | |
309 | // ,fInvmassULS5(0) | |
520a29a8 | 310 | ,fcentrality(0) |
311 | ,fElecPhi(0) | |
3c56855b | 312 | ,fElecPhiTPChalf(0) |
bf697dda | 313 | ,fElecPhiPt(0) |
314 | // ,fElecPhiTPC(0) | |
315 | // ,fElecPhiTPCEovP(0) | |
520a29a8 | 316 | ,fHadronPhi(0) |
3c56855b | 317 | ,fHadronPhiTPChalf(0) |
bf697dda | 318 | ,fHadronPhiPt(0) |
319 | /* ,fTrackHFEcuts(0) | |
320 | ,fTrakPhiSPD1(0) | |
321 | ,fTrakPhiSPD2(0) | |
322 | ,fTrakPhiSPDOr(0) | |
323 | ,fTrakPhiSPDAnd(0) | |
324 | ,fTrackHFEcutsITS(0) | |
325 | */ | |
3db57d71 | 326 | ,fNoMixedEvents(0) |
a9123929 | 327 | ,fMixStat(0) |
328 | ,fMixStat1(0) | |
329 | ,fMixedIncElecDphi(0) | |
330 | ,fMixedIncElecDphi1(0) | |
331 | ,fMixedIncElecDphi2(0) | |
332 | ,fMixedIncElecDphi3(0) | |
333 | ,fMixedIncElecDphi4(0) | |
334 | ,fMixedPhotElecDphi(0) | |
335 | ,fMixedPhotElecDphi1(0) | |
336 | ,fMixedPhotElecDphi2(0) | |
337 | ,fMixedPhotElecDphi3(0) | |
338 | ,fMixedPhotElecDphi4(0) | |
339 | ,fMixedSemiIncElecDphi(0) | |
340 | ,fMixedSemiIncElecDphi1(0) | |
341 | ,fMixedSemiIncElecDphi2(0) | |
342 | ,fMixedSemiIncElecDphi3(0) | |
343 | ,fMixedSemiIncElecDphi4(0) | |
344 | ,fMixedDphiULSMassLow(0) | |
345 | ,fMixedDphiULSMassLow1(0) | |
346 | ,fMixedDphiULSMassLow2(0) | |
347 | ,fMixedDphiULSMassLow3(0) | |
348 | ,fMixedDphiULSMassLow4(0) | |
349 | ,fMixedDphiLSMassLow(0) | |
350 | ,fMixedDphiLSMassLow1(0) | |
20791315 | 351 | ,fMixedDphiLSMassLow2(0) |
352 | ,fMixedDphiLSMassLow3(0) | |
353 | ,fMixedDphiLSMassLow4(0) | |
354 | ,fDCAMetMixedPhotElecDphi(0) | |
355 | ,fDCAMetMixedPhotElecDphi1(0) | |
356 | ,fDCAMetMixedPhotElecDphi2(0) | |
357 | ,fDCAMetMixedPhotElecDphi3(0) | |
358 | ,fDCAMetMixedPhotElecDphi4(0) | |
359 | ,fDCAMetMixedDphiULSMassLow(0) | |
360 | ,fDCAMetMixedDphiULSMassLow1(0) | |
361 | ,fDCAMetMixedDphiULSMassLow2(0) | |
362 | ,fDCAMetMixedDphiULSMassLow3(0) | |
363 | ,fDCAMetMixedDphiULSMassLow4(0) | |
364 | ,fDCAMetMixedDphiLSMassLow(0) | |
365 | ,fDCAMetMixedDphiLSMassLow1(0) | |
366 | ,fDCAMetMixedDphiLSMassLow2(0) | |
367 | ,fDCAMetMixedDphiLSMassLow3(0) | |
368 | ,fDCAMetMixedDphiLSMassLow4(0) | |
3c56855b | 369 | ,fHadronPt(0) |
370 | ,fCentralityPass(0) | |
371 | ,fCentralityNoPass(0) | |
0841d73f | 372 | ,fHadronDphi(0) |
373 | ,fHadronDphi1(0) | |
374 | ,fHadronDphi2(0) | |
375 | ,fHadronDphi3(0) | |
376 | ,fHadronDphi4(0) | |
377 | ,fPiPt(0) | |
33369c14 | 378 | ,fHadronDphiNoSS(0) |
379 | ,fHadronDphiNoSS1(0) | |
380 | ,fHadronDphiNoSS2(0) | |
381 | ,fHadronDphiNoSS3(0) | |
382 | ,fHadronDphiNoSS4(0) | |
383 | ,fPiPtNoSS(0) | |
bf697dda | 384 | ,fEovPWoSS(0) |
385 | ,fEovPWSS(0) | |
386 | ,fEovPHadWoSS(0) | |
387 | ,fEovPHadWSS(0) | |
0c884279 | 388 | ,fHadronDphiEta1(0) |
389 | ,fHadronDphiEta11(0) | |
390 | ,fHadronDphiEta12(0) | |
391 | ,fHadronDphiEta13(0) | |
392 | ,fHadronDphiEta14(0) | |
393 | ,fHadronDphiNoSSEta1(0) | |
394 | ,fHadronDphiNoSSEta11(0) | |
395 | ,fHadronDphiNoSSEta12(0) | |
396 | ,fHadronDphiNoSSEta13(0) | |
397 | ,fHadronDphiNoSSEta14(0) | |
398 | ,fHadronDphiEta2(0) | |
399 | ,fHadronDphiEta21(0) | |
400 | ,fHadronDphiEta22(0) | |
401 | ,fHadronDphiEta23(0) | |
402 | ,fHadronDphiEta24(0) | |
403 | ,fHadronDphiNoSSEta2(0) | |
404 | ,fHadronDphiNoSSEta21(0) | |
405 | ,fHadronDphiNoSSEta22(0) | |
406 | ,fHadronDphiNoSSEta23(0) | |
20791315 | 407 | ,fHadronDphiNoSSEta24(0) |
408 | //,fSparseElectron(0) | |
409 | // ,fvalueElectron(0) | |
57189f04 | 410 | { |
411 | //Named constructor | |
520a29a8 | 412 | |
57189f04 | 413 | fPID = new AliHFEpid("hfePid"); |
bf697dda | 414 | //fvalueElectron = new Double_t[6]; |
57189f04 | 415 | |
416 | // Define input and output slots here | |
417 | // Input slot #0 works with a TChain | |
418 | DefineInput(0, TChain::Class()); | |
419 | // Output slot #0 id reserved by the base class for AOD | |
420 | // Output slot #1 writes into a TH1 container | |
421 | // DefineOutput(1, TH1I::Class()); | |
422 | DefineOutput(1, TList::Class()); | |
423 | // DefineOutput(3, TTree::Class()); | |
424 | } | |
425 | ||
426 | //________________________________________________________________________ | |
427 | AliAnalysisTaskElecHadronCorrel::AliAnalysisTaskElecHadronCorrel() | |
428 | : AliAnalysisTaskSE("DefaultAnalysis_AliAnalysisElecHadCorrel") | |
987053ce | 429 | ,fVevent(0) |
57189f04 | 430 | ,fESD(0) |
987053ce | 431 | ,fAOD(0) |
57189f04 | 432 | ,fGeom(0) |
20791315 | 433 | ,fpidResponse(0) |
57189f04 | 434 | ,fOutputList(0) |
520a29a8 | 435 | ,fTrackCuts1(new AliESDtrackCuts) |
436 | ,fTrackCuts2(new AliESDtrackCuts) | |
57189f04 | 437 | ,fCuts(0) |
438 | ,fIdentifiedAsOutInz(kFALSE) | |
439 | ,fPassTheEventCut(kFALSE) | |
440 | ,fRejectKinkMother(kFALSE) | |
4e01b68c | 441 | ,fVz(0.0) |
442 | ,fCFM(0) | |
443 | ,fPID(0) | |
444 | ,fPIDqa(0) | |
57189f04 | 445 | ,fInvmassCut(0.01) |
3c56855b | 446 | ,fCentrality(0) |
447 | ,fCentralityMin(0) | |
448 | ,fCentralityMax(0) | |
a9123929 | 449 | ,fkCentralityMethod(0) |
450 | ,fTPCnsigEleMin(-2) | |
451 | ,fTPCnsigEleMax(2) | |
452 | ,fTPCnsigHadMin(-10) | |
453 | ,fTPCnsigHadMax(-3.5) | |
454 | ,fM02CutMin(0.03) | |
455 | ,fM02CutMax(0.5) | |
456 | ,fM20CutMin(0.03) | |
457 | ,fM20CutMax(0.3) | |
458 | ,fDispCutMin(0) | |
459 | ,fDispCutMax(1) | |
460 | ,fEovPMin(0.8) | |
461 | ,fEovPMax(1.2) | |
462 | ,fTriggerCentral(kTRUE) | |
6ec70f8c | 463 | ,fTriggerMB(kTRUE) |
28ea601a | 464 | ,fTPCNClsHad(80) |
465 | ,fAssoEleITSref(kTRUE) | |
ca49cfe0 | 466 | ,fAssoElecTPCNCls(80) |
20791315 | 467 | ,fNonHFE(0) |
3db57d71 | 468 | ,fPoolMgr(0x0) |
bf697dda | 469 | ,fNoEvents(0) |
470 | // ,fTrkpt(0) | |
471 | ,fTrkEovPAft(0) | |
472 | // ,fTrkEovPBefHad(0) | |
473 | // ,fdEdxBef(0) | |
987053ce | 474 | ,fSemiIncElecDphi(0) |
487ae960 | 475 | ,fSemiIncElecDphi1(0) |
476 | ,fSemiIncElecDphi2(0) | |
477 | ,fSemiIncElecDphi3(0) | |
478 | ,fSemiIncElecDphi4(0) | |
57189f04 | 479 | ,fPhotElecDphi(0) |
487ae960 | 480 | ,fPhotElecDphi1(0) |
481 | ,fPhotElecDphi2(0) | |
482 | ,fPhotElecDphi3(0) | |
483 | ,fPhotElecDphi4(0) | |
57189f04 | 484 | ,fInclusiveElecDphi(0) |
487ae960 | 485 | ,fInclusiveElecDphi1(0) |
486 | ,fInclusiveElecDphi2(0) | |
487 | ,fInclusiveElecDphi3(0) | |
488 | ,fInclusiveElecDphi4(0) | |
33369c14 | 489 | ,fInclusiveElecDphiEtaFS(0) |
490 | ,fInclusiveElecDphiEtaFS1(0) | |
491 | ,fInclusiveElecDphiEtaFS2(0) | |
492 | ,fInclusiveElecDphiEtaFS3(0) | |
493 | ,fInclusiveElecDphiEtaFS4(0) | |
57189f04 | 494 | ,fDphiULSMassLow(0) |
487ae960 | 495 | ,fDphiULSMassLow1(0) |
496 | ,fDphiULSMassLow2(0) | |
497 | ,fDphiULSMassLow3(0) | |
498 | ,fDphiULSMassLow4(0) | |
57189f04 | 499 | ,fDphiLSMassLow(0) |
487ae960 | 500 | ,fDphiLSMassLow1(0) |
501 | ,fDphiLSMassLow2(0) | |
502 | ,fDphiLSMassLow3(0) | |
503 | ,fDphiLSMassLow4(0) | |
520a29a8 | 504 | ,fDphiULSMassLowNoPartner(0) |
487ae960 | 505 | ,fDphiULSMassLowNoPartner1(0) |
506 | ,fDphiULSMassLowNoPartner2(0) | |
507 | ,fDphiULSMassLowNoPartner3(0) | |
508 | ,fDphiULSMassLowNoPartner4(0) | |
520a29a8 | 509 | ,fDphiLSMassLowNoPartner(0) |
487ae960 | 510 | ,fDphiLSMassLowNoPartner1(0) |
511 | ,fDphiLSMassLowNoPartner2(0) | |
512 | ,fDphiLSMassLowNoPartner3(0) | |
513 | ,fDphiLSMassLowNoPartner4(0) | |
57189f04 | 514 | ,fPhotoElecPt(0) |
515 | ,fSemiInclElecPt(0) | |
520a29a8 | 516 | ,fInclusiveElecPt(0) |
517 | ,fULSElecPt(0) | |
518 | ,fLSElecPt(0) | |
20791315 | 519 | ,fDCAMetPhotElecDphi(0) |
520 | ,fDCAMetPhotElecDphi1(0) | |
521 | ,fDCAMetPhotElecDphi2(0) | |
522 | ,fDCAMetPhotElecDphi3(0) | |
523 | ,fDCAMetPhotElecDphi4(0) | |
524 | ,fDCAMetDphiULSMassLow(0) | |
525 | ,fDCAMetDphiULSMassLow1(0) | |
526 | ,fDCAMetDphiULSMassLow2(0) | |
527 | ,fDCAMetDphiULSMassLow3(0) | |
528 | ,fDCAMetDphiULSMassLow4(0) | |
529 | ,fDCAMetDphiLSMassLow(0) | |
530 | ,fDCAMetDphiLSMassLow1(0) | |
531 | ,fDCAMetDphiLSMassLow2(0) | |
532 | ,fDCAMetDphiLSMassLow3(0) | |
533 | ,fDCAMetDphiLSMassLow4(0) | |
534 | ,fDCAMetDphiULSMassLowNoPartner(0) | |
535 | ,fDCAMetDphiULSMassLowNoPartner1(0) | |
536 | ,fDCAMetDphiULSMassLowNoPartner2(0) | |
537 | ,fDCAMetDphiULSMassLowNoPartner3(0) | |
538 | ,fDCAMetDphiULSMassLowNoPartner4(0) | |
539 | ,fDCAMetDphiLSMassLowNoPartner(0) | |
540 | ,fDCAMetDphiLSMassLowNoPartner1(0) | |
541 | ,fDCAMetDphiLSMassLowNoPartner2(0) | |
542 | ,fDCAMetDphiLSMassLowNoPartner3(0) | |
543 | ,fDCAMetDphiLSMassLowNoPartner4(0) | |
544 | ,fDCAMetPhotoElecPt(0) | |
545 | ,fDCAMetULSElecPt(0) | |
546 | ,fDCAMetLSElecPt(0) | |
bf697dda | 547 | ,fSemiIncElecDphiEta1(0) |
548 | ,fSemiIncElecDphiEta11(0) | |
549 | ,fSemiIncElecDphiEta12(0) | |
550 | ,fSemiIncElecDphiEta13(0) | |
551 | ,fSemiIncElecDphiEta14(0) | |
552 | ,fPhotElecDphiEta1(0) | |
553 | ,fPhotElecDphiEta11(0) | |
554 | ,fPhotElecDphiEta12(0) | |
555 | ,fPhotElecDphiEta13(0) | |
556 | ,fPhotElecDphiEta14(0) | |
557 | ,fInclusiveElecDphiEta1(0) | |
558 | ,fInclusiveElecDphiEta11(0) | |
559 | ,fInclusiveElecDphiEta12(0) | |
560 | ,fInclusiveElecDphiEta13(0) | |
561 | ,fInclusiveElecDphiEta14(0) | |
562 | ,fDphiULSMassLowEta1(0) | |
563 | ,fDphiULSMassLowEta11(0) | |
564 | ,fDphiULSMassLowEta12(0) | |
565 | ,fDphiULSMassLowEta13(0) | |
566 | ,fDphiULSMassLowEta14(0) | |
567 | ,fDphiLSMassLowEta1(0) | |
568 | ,fDphiLSMassLowEta11(0) | |
569 | ,fDphiLSMassLowEta12(0) | |
570 | ,fDphiLSMassLowEta13(0) | |
571 | ,fDphiLSMassLowEta14(0) | |
572 | ,fDphiULSMassLowNoPartnerEta1(0) | |
573 | ,fDphiULSMassLowNoPartnerEta11(0) | |
574 | ,fDphiULSMassLowNoPartnerEta12(0) | |
575 | ,fDphiULSMassLowNoPartnerEta13(0) | |
576 | ,fDphiULSMassLowNoPartnerEta14(0) | |
577 | ,fDphiLSMassLowNoPartnerEta1(0) | |
578 | ,fDphiLSMassLowNoPartnerEta11(0) | |
579 | ,fDphiLSMassLowNoPartnerEta12(0) | |
580 | ,fDphiLSMassLowNoPartnerEta13(0) | |
581 | ,fDphiLSMassLowNoPartnerEta14(0) | |
582 | ,fSemiIncElecDphiEta2(0) | |
583 | ,fSemiIncElecDphiEta21(0) | |
584 | ,fSemiIncElecDphiEta22(0) | |
585 | ,fSemiIncElecDphiEta23(0) | |
586 | ,fSemiIncElecDphiEta24(0) | |
587 | ,fPhotElecDphiEta2(0) | |
588 | ,fPhotElecDphiEta21(0) | |
589 | ,fPhotElecDphiEta22(0) | |
590 | ,fPhotElecDphiEta23(0) | |
591 | ,fPhotElecDphiEta24(0) | |
592 | ,fInclusiveElecDphiEta2(0) | |
593 | ,fInclusiveElecDphiEta21(0) | |
594 | ,fInclusiveElecDphiEta22(0) | |
595 | ,fInclusiveElecDphiEta23(0) | |
596 | ,fInclusiveElecDphiEta24(0) | |
597 | ,fDphiULSMassLowEta2(0) | |
598 | ,fDphiULSMassLowEta21(0) | |
599 | ,fDphiULSMassLowEta22(0) | |
600 | ,fDphiULSMassLowEta23(0) | |
601 | ,fDphiULSMassLowEta24(0) | |
602 | ,fDphiLSMassLowEta2(0) | |
603 | ,fDphiLSMassLowEta21(0) | |
604 | ,fDphiLSMassLowEta22(0) | |
605 | ,fDphiLSMassLowEta23(0) | |
606 | ,fDphiLSMassLowEta24(0) | |
607 | ,fDphiULSMassLowNoPartnerEta2(0) | |
608 | ,fDphiULSMassLowNoPartnerEta21(0) | |
609 | ,fDphiULSMassLowNoPartnerEta22(0) | |
610 | ,fDphiULSMassLowNoPartnerEta23(0) | |
611 | ,fDphiULSMassLowNoPartnerEta24(0) | |
612 | ,fDphiLSMassLowNoPartnerEta2(0) | |
613 | ,fDphiLSMassLowNoPartnerEta21(0) | |
614 | ,fDphiLSMassLowNoPartnerEta22(0) | |
615 | ,fDphiLSMassLowNoPartnerEta23(0) | |
616 | ,fDphiLSMassLowNoPartnerEta24(0) | |
617 | // ,fTrackPtBefTrkCuts(0) | |
618 | // ,fTrackPtAftTrkCuts(0) | |
20791315 | 619 | ,fTPCnsigma(0) |
bf697dda | 620 | // ,fNCellv1(0) |
621 | // ,fClsEv1(0) | |
622 | // ,fNClusv1(0) | |
623 | ,fInvmassLS1(0) | |
624 | // ,fInvmassLS2(0) | |
625 | // ,fInvmassLS3(0) | |
626 | // ,fInvmassLS4(0) | |
627 | // ,fInvmassLS5(0) | |
628 | ,fInvmassULS1(0) | |
629 | // ,fInvmassULS2(0) | |
630 | // ,fInvmassULS3(0) | |
631 | // ,fInvmassULS4(0) | |
632 | // ,fInvmassULS5(0) | |
20791315 | 633 | ,fDCAMetInvmassLS1(0) |
634 | ,fDCAMetInvmassULS1(0) | |
520a29a8 | 635 | ,fcentrality(0) |
636 | ,fElecPhi(0) | |
3c56855b | 637 | ,fElecPhiTPChalf(0) |
bf697dda | 638 | ,fElecPhiPt(0) |
639 | // ,fElecPhiTPC(0) | |
640 | // ,fElecPhiTPCEovP(0) | |
520a29a8 | 641 | ,fHadronPhi(0) |
3c56855b | 642 | ,fHadronPhiTPChalf(0) |
bf697dda | 643 | ,fHadronPhiPt(0) |
644 | /* ,fTrackHFEcuts(0) | |
645 | ,fTrakPhiSPD1(0) | |
646 | ,fTrakPhiSPD2(0) | |
647 | ,fTrakPhiSPDOr(0) | |
648 | ,fTrakPhiSPDAnd(0) | |
649 | ,fTrackHFEcutsITS(0) | |
650 | */ | |
3db57d71 | 651 | ,fNoMixedEvents(0) |
a9123929 | 652 | ,fMixStat(0) |
653 | ,fMixStat1(0) | |
654 | ,fMixedIncElecDphi(0) | |
655 | ,fMixedIncElecDphi1(0) | |
656 | ,fMixedIncElecDphi2(0) | |
657 | ,fMixedIncElecDphi3(0) | |
658 | ,fMixedIncElecDphi4(0) | |
659 | ,fMixedPhotElecDphi(0) | |
660 | ,fMixedPhotElecDphi1(0) | |
661 | ,fMixedPhotElecDphi2(0) | |
662 | ,fMixedPhotElecDphi3(0) | |
663 | ,fMixedPhotElecDphi4(0) | |
664 | ,fMixedSemiIncElecDphi(0) | |
665 | ,fMixedSemiIncElecDphi1(0) | |
666 | ,fMixedSemiIncElecDphi2(0) | |
667 | ,fMixedSemiIncElecDphi3(0) | |
668 | ,fMixedSemiIncElecDphi4(0) | |
669 | ,fMixedDphiULSMassLow(0) | |
670 | ,fMixedDphiULSMassLow1(0) | |
671 | ,fMixedDphiULSMassLow2(0) | |
672 | ,fMixedDphiULSMassLow3(0) | |
673 | ,fMixedDphiULSMassLow4(0) | |
674 | ,fMixedDphiLSMassLow(0) | |
675 | ,fMixedDphiLSMassLow1(0) | |
676 | ,fMixedDphiLSMassLow2(0) | |
677 | ,fMixedDphiLSMassLow3(0) | |
20791315 | 678 | ,fMixedDphiLSMassLow4(0) |
679 | ,fDCAMetMixedPhotElecDphi(0) | |
680 | ,fDCAMetMixedPhotElecDphi1(0) | |
681 | ,fDCAMetMixedPhotElecDphi2(0) | |
682 | ,fDCAMetMixedPhotElecDphi3(0) | |
683 | ,fDCAMetMixedPhotElecDphi4(0) | |
684 | ,fDCAMetMixedDphiULSMassLow(0) | |
685 | ,fDCAMetMixedDphiULSMassLow1(0) | |
686 | ,fDCAMetMixedDphiULSMassLow2(0) | |
687 | ,fDCAMetMixedDphiULSMassLow3(0) | |
688 | ,fDCAMetMixedDphiULSMassLow4(0) | |
689 | ,fDCAMetMixedDphiLSMassLow(0) | |
690 | ,fDCAMetMixedDphiLSMassLow1(0) | |
691 | ,fDCAMetMixedDphiLSMassLow2(0) | |
692 | ,fDCAMetMixedDphiLSMassLow3(0) | |
693 | ,fDCAMetMixedDphiLSMassLow4(0) | |
3c56855b | 694 | ,fHadronPt(0) |
695 | ,fCentralityPass(0) | |
696 | ,fCentralityNoPass(0) | |
0841d73f | 697 | ,fHadronDphi(0) |
698 | ,fHadronDphi1(0) | |
699 | ,fHadronDphi2(0) | |
700 | ,fHadronDphi3(0) | |
701 | ,fHadronDphi4(0) | |
33369c14 | 702 | ,fPiPt(0) |
703 | ,fHadronDphiNoSS(0) | |
704 | ,fHadronDphiNoSS1(0) | |
705 | ,fHadronDphiNoSS2(0) | |
706 | ,fHadronDphiNoSS3(0) | |
707 | ,fHadronDphiNoSS4(0) | |
708 | ,fPiPtNoSS(0) | |
bf697dda | 709 | ,fEovPWoSS(0) |
710 | ,fEovPWSS(0) | |
711 | ,fEovPHadWoSS(0) | |
0c884279 | 712 | ,fEovPHadWSS(0) |
713 | ,fHadronDphiEta1(0) | |
714 | ,fHadronDphiEta11(0) | |
715 | ,fHadronDphiEta12(0) | |
716 | ,fHadronDphiEta13(0) | |
717 | ,fHadronDphiEta14(0) | |
718 | ,fHadronDphiNoSSEta1(0) | |
719 | ,fHadronDphiNoSSEta11(0) | |
720 | ,fHadronDphiNoSSEta12(0) | |
721 | ,fHadronDphiNoSSEta13(0) | |
722 | ,fHadronDphiNoSSEta14(0) | |
723 | ,fHadronDphiEta2(0) | |
724 | ,fHadronDphiEta21(0) | |
725 | ,fHadronDphiEta22(0) | |
726 | ,fHadronDphiEta23(0) | |
727 | ,fHadronDphiEta24(0) | |
728 | ,fHadronDphiNoSSEta2(0) | |
729 | ,fHadronDphiNoSSEta21(0) | |
730 | ,fHadronDphiNoSSEta22(0) | |
731 | ,fHadronDphiNoSSEta23(0) | |
20791315 | 732 | ,fHadronDphiNoSSEta24(0) |
733 | //,fSparseElectron(0) | |
734 | // ,fvalueElectron(0) | |
57189f04 | 735 | { |
520a29a8 | 736 | //Default constructor |
737 | fPID = new AliHFEpid("hfePid"); | |
bf697dda | 738 | //fvalueElectron = new Double_t[6]; |
57189f04 | 739 | |
520a29a8 | 740 | // Constructor |
4e01b68c | 741 | // Define input and output slots here |
742 | // Input slot #0 works with a TChain | |
743 | DefineInput(0, TChain::Class()); | |
744 | // Output slot #0 id reserved by the base class for AOD | |
745 | // Output slot #1 writes into a TH1 container | |
746 | // DefineOutput(1, TH1I::Class()); | |
747 | DefineOutput(1, TList::Class()); | |
748 | //DefineOutput(3, TTree::Class()); | |
57189f04 | 749 | } |
750 | //_________________________________________ | |
751 | ||
752 | AliAnalysisTaskElecHadronCorrel::~AliAnalysisTaskElecHadronCorrel() | |
753 | { | |
4e01b68c | 754 | //Destructor |
755 | ||
756 | delete fOutputList; | |
757 | delete fGeom; | |
758 | delete fPID; | |
759 | delete fCFM; | |
760 | delete fPIDqa; | |
761 | delete fTrackCuts1; | |
762 | delete fTrackCuts2; | |
bf697dda | 763 | // delete fSparseElectron; |
764 | // delete []fvalueElectron; | |
57189f04 | 765 | } |
766 | //_________________________________________ | |
767 | ||
768 | void AliAnalysisTaskElecHadronCorrel::UserExec(Option_t*) | |
769 | { | |
4e01b68c | 770 | //Main loop |
771 | //Called for each event | |
772 | ||
773 | // create pointer to event | |
987053ce | 774 | fAOD = dynamic_cast<AliAODEvent*>(InputEvent()); |
4e01b68c | 775 | fESD = dynamic_cast<AliESDEvent*>(InputEvent()); |
bf697dda | 776 | |
987053ce | 777 | if(!(fESD || fAOD)){ |
778 | printf("ERROR: fESD & fAOD not available\n"); | |
779 | return; | |
780 | } | |
781 | fVevent = dynamic_cast<AliVEvent*>(InputEvent()); | |
782 | if (!fVevent) { | |
783 | printf("ERROR: fVEvent not available\n"); | |
4e01b68c | 784 | return; |
785 | } | |
520a29a8 | 786 | |
4e01b68c | 787 | if(!fCuts){ |
788 | AliError("HFE cuts not available"); | |
789 | return; | |
790 | } | |
791 | ||
792 | if(!fPID->IsInitialized()){ | |
793 | // Initialize PID with the given run number | |
794 | AliWarning("PID not initialised, get from Run no"); | |
a9123929 | 795 | |
987053ce | 796 | if(IsAODanalysis())fPID->InitializePID(fAOD->GetRunNumber()); |
797 | else fPID->InitializePID(fESD->GetRunNumber()); | |
4e01b68c | 798 | } |
799 | ||
520a29a8 | 800 | |
6ec70f8c | 801 | if(fTriggerMB && fTriggerCentral){ |
a9123929 | 802 | // trigger selection |
803 | if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & (AliVEvent::kCentral))) return; | |
804 | } | |
805 | ||
6ec70f8c | 806 | if(fTriggerMB && !fTriggerCentral){ |
a9123929 | 807 | // trigger selection |
808 | if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & (AliVEvent::kSemiCentral))) return; | |
809 | } | |
20791315 | 810 | |
6ec70f8c | 811 | if(!fTriggerMB){ |
20791315 | 812 | // trigger selection |
813 | if(!(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected() & (AliVEvent::kEMCEGA))) return; | |
6ec70f8c | 814 | } |
987053ce | 815 | // centrality selection |
b8c2fb6d | 816 | // SetCentralityParameters(0., 7., fkCentralityMethod); |
987053ce | 817 | Bool_t pass = kFALSE; |
818 | CheckCentrality(fVevent,pass); | |
e4b0faf2 | 819 | if(!pass)return; |
820 | ||
987053ce | 821 | Int_t fNOtrks = fVevent->GetNumberOfTracks(); |
822 | const AliVVertex *pVtx = fVevent->GetPrimaryVertex(); | |
520a29a8 | 823 | |
4e01b68c | 824 | Double_t pVtxZ = -999; |
825 | pVtxZ = pVtx->GetZ(); | |
520a29a8 | 826 | |
4e01b68c | 827 | if(TMath::Abs(pVtxZ)>10) return; |
828 | fNoEvents->Fill(0); | |
829 | ||
830 | if(fNOtrks<2) return; | |
831 | ||
20791315 | 832 | fpidResponse = fInputHandler->GetPIDResponse(); |
833 | if(!fpidResponse){ | |
4e01b68c | 834 | AliDebug(1, "Using default PID Response"); |
20791315 | 835 | fpidResponse = AliHFEtools::GetDefaultPID(kFALSE, fInputEvent->IsA() == AliAODEvent::Class()); |
a9123929 | 836 | } |
837 | ||
20791315 | 838 | fPID->SetPIDResponse(fpidResponse); |
a9123929 | 839 | |
840 | fCFM->SetRecEventInfo(fVevent); | |
841 | ||
842 | AliCentrality *fCentrality2 = (AliCentrality*)fAOD->GetCentrality(); | |
b8c2fb6d | 843 | Double_t centvalue1 = fCentrality2->GetCentralityPercentile(fkCentralityMethod); |
3db57d71 | 844 | |
a9123929 | 845 | //Event mixing |
7e02de85 | 846 | AliEventPool* pool; |
847 | pool = fPoolMgr->GetEventPool(centvalue1, pVtxZ); // Get the buffer associated with the current centrality and z-vtx | |
a9123929 | 848 | if (!pool) |
7e02de85 | 849 | { |
a9123929 | 850 | AliFatal(Form("No pool found for centrality = %f, zVtx = %f", centvalue1, pVtxZ)); |
7e02de85 | 851 | return; |
852 | } | |
a9123929 | 853 | // Look for kink mother for AOD |
28ea601a | 854 | Int_t numberofvertices = fAOD->GetNumberOfVertices(); |
855 | Double_t listofmotherkink[numberofvertices]; | |
856 | Int_t numberofmotherkink = 0; | |
a9123929 | 857 | if(IsAODanalysis()){ |
a9123929 | 858 | for(Int_t ivertex=0; ivertex < numberofvertices; ivertex++) { |
859 | AliAODVertex *aodvertex = fAOD->GetVertex(ivertex); | |
860 | if(!aodvertex) continue; | |
861 | if(aodvertex->GetType()==AliAODVertex::kKink) { | |
862 | AliAODTrack *mother = (AliAODTrack *) aodvertex->GetParent(); | |
863 | if(!mother) continue; | |
864 | Int_t idmother = mother->GetID(); | |
865 | listofmotherkink[numberofmotherkink] = idmother; | |
866 | numberofmotherkink++; | |
867 | } | |
868 | } | |
869 | } | |
870 | ||
871 | // Track loop | |
872 | for (Int_t iTracks = 0; iTracks < fVevent->GetNumberOfTracks(); iTracks++) { | |
873 | AliVParticle* Vtrack = fVevent->GetTrack(iTracks); | |
874 | if (!Vtrack) { | |
875 | printf("ERROR: Could not receive track %d\n", iTracks); | |
876 | continue; | |
877 | } | |
878 | AliVTrack *track = dynamic_cast<AliVTrack*>(Vtrack); | |
879 | AliESDtrack *etrack = dynamic_cast<AliESDtrack*>(Vtrack); | |
880 | AliAODTrack *atrack = dynamic_cast<AliAODTrack*>(Vtrack); | |
881 | ||
7e02de85 | 882 | if(!track && !etrack && !atrack) continue; |
883 | ||
a9123929 | 884 | if(IsAODanalysis()) |
28ea601a | 885 | if(atrack){ |
886 | if(!atrack->TestFilterMask(AliAODTrack::kTrkGlobalNoDCA)) continue; | |
887 | } | |
a9123929 | 888 | if(track->Pt()<1) continue; |
889 | ||
890 | // fTrackPtBefTrkCuts->Fill(track->Pt()); | |
891 | ||
892 | // RecKine: ITSTPC cuts | |
893 | if(!ProcessCutStep(AliHFEcuts::kStepRecKineITSTPC, track)) continue; | |
894 | ||
895 | // Reject kink mother | |
896 | if(fRejectKinkMother) { // Quick and dirty fix to reject both kink mothers and daughters | |
897 | if(IsAODanalysis()){ | |
898 | Bool_t kinkmotherpass = kTRUE; | |
899 | for(Int_t kinkmother = 0; kinkmother < numberofmotherkink; kinkmother++) { | |
900 | if(track->GetID() == listofmotherkink[kinkmother]) { | |
901 | kinkmotherpass = kFALSE; | |
902 | continue; | |
903 | } | |
bf697dda | 904 | } |
a9123929 | 905 | if(!kinkmotherpass) continue; |
906 | } | |
907 | else{ | |
908 | if(etrack->GetKinkIndex(0) != 0) continue; | |
909 | } | |
910 | } | |
911 | // RecPrim | |
912 | // if(!ProcessCutStep(AliHFEcuts::kStepRecPrim, track)) continue; //gives warning for AOD, so not using | |
913 | ||
914 | // HFE cuts: TPC PID cleanup | |
915 | if(!ProcessCutStep(AliHFEcuts::kStepHFEcutsTPC, track)) continue; | |
916 | ||
917 | // fTrackHFEcuts->Fill(track->Phi()); | |
918 | ||
919 | // HFEcuts: ITS layers cuts | |
920 | if(!ProcessCutStep(AliHFEcuts::kStepHFEcutsITS, track)) continue; | |
921 | ||
922 | // fTrackHFEcutsITS->Fill(track->Phi()); | |
923 | //fTrackPtAftTrkCuts->Fill(track->Pt()); | |
924 | ||
925 | Double_t fClsE = -999, p = -999, fEovP=-999, pt = -999, dEdx=-999, fTPCnSigma=0; | |
926 | pt = track->Pt(); | |
927 | p = track->P(); | |
928 | dEdx = track->GetTPCsignal(); | |
929 | fTPCnSigma = fPID->GetPIDResponse() ? fPID->GetPIDResponse()->NumberOfSigmasTPC(track, AliPID::kElectron) : 1000; | |
930 | ||
20791315 | 931 | if(pt<2) continue; |
a9123929 | 932 | //TPC electron phi |
933 | // if(fTPCnSigma >= -2 && fTPCnSigma <= 2){ | |
934 | // fElecPhiTPC->Fill(track->Phi()); | |
935 | // } | |
936 | ||
937 | //eta cut (-0.7,0.7) | |
938 | if(track->Eta() < -0.7 || track->Eta() > 0.7) continue; | |
939 | ||
940 | // Track extrapolation to EMCAL | |
941 | Int_t fClsId = track->GetEMCALcluster(); | |
942 | if(fClsId <0) continue; | |
943 | AliVCluster *cluster = fVevent->GetCaloCluster(fClsId); | |
944 | if(!cluster->IsEMCAL()) continue; | |
945 | if(TMath::Abs(cluster->GetTrackDx())>0.05 || TMath::Abs(cluster->GetTrackDz())>0.05) continue; | |
946 | // fdEdxBef->Fill(p,dEdx); | |
7e02de85 | 947 | fTPCnsigma->Fill(p,fTPCnSigma); |
a9123929 | 948 | |
949 | // fTrkpt->Fill(pt); | |
950 | fClsE = cluster->E(); | |
951 | fEovP = fClsE/p; | |
952 | ||
953 | //--------THnsparse--------- | |
954 | //fvalueElectron[0] = pt; | |
955 | //fvalueElectron[1] = fTPCnSigma; | |
956 | //fvalueElectron[2] = fEovP; | |
957 | //fvalueElectron[3] = cluster->GetM20(); | |
958 | //fvalueElectron[4] = cluster->GetM02(); | |
959 | //fvalueElectron[5] = cluster->GetDispersion(); | |
960 | ||
961 | //fSparseElectron->Fill(fvalueElectron); | |
962 | ||
963 | //---------------- | |
964 | //EovP distribution for Had | |
965 | if((fTPCnSigma > fTPCnsigHadMin) && (fTPCnSigma < fTPCnsigHadMax)) fEovPHadWoSS->Fill(pt,fEovP); | |
966 | if(((fTPCnSigma > fTPCnsigHadMin) && (fTPCnSigma < fTPCnsigHadMax)) && ((cluster->GetM20()> fM20CutMin) && (cluster->GetM20()< fM20CutMax)) && ((cluster->GetM02()> fM02CutMin) && (cluster->GetM02()< fM02CutMax)) && ((cluster->GetDispersion()< fDispCutMax))) fEovPHadWSS->Fill(pt,fEovP); | |
967 | ||
968 | //Dphi distribution hadrons with shower shape cuts | |
969 | if(((fTPCnSigma > fTPCnsigHadMin) && (fTPCnSigma < fTPCnsigHadMax)) && ((cluster->GetM20()> fM20CutMin) && (cluster->GetM20()< fM20CutMax)) && ((cluster->GetM02()> fM02CutMin) && (cluster->GetM02()< fM02CutMax)) && ((cluster->GetDispersion()< fDispCutMax))&&(fEovP >= fEovPMin && fEovP <= fEovPMax)){ | |
970 | ElectronHadCorrel(iTracks, track, fHadronDphi, fHadronDphi1,fHadronDphi2,fHadronDphi3,fHadronDphi4); | |
971 | ElectronHadCorrelEtaBins(iTracks, track, fHadronDphiEta1, fHadronDphiEta11,fHadronDphiEta12,fHadronDphiEta13,fHadronDphiEta14,fHadronDphiEta2, fHadronDphiEta21,fHadronDphiEta22,fHadronDphiEta23,fHadronDphiEta24); | |
972 | fPiPt->Fill(pt); | |
973 | } | |
974 | ||
975 | //Dphi distribution hadrons without shower shape cuts | |
976 | if((fTPCnSigma > fTPCnsigHadMin) && (fTPCnSigma < fTPCnsigHadMax)&&(fEovP >= fEovPMin && fEovP <= fEovPMax)){ | |
977 | ElectronHadCorrel(iTracks, track, fHadronDphiNoSS, fHadronDphiNoSS1,fHadronDphiNoSS2,fHadronDphiNoSS3,fHadronDphiNoSS4); | |
978 | ElectronHadCorrelEtaBins(iTracks, track, fHadronDphiNoSSEta1, fHadronDphiNoSSEta11,fHadronDphiNoSSEta12,fHadronDphiNoSSEta13,fHadronDphiNoSSEta14,fHadronDphiNoSSEta2, fHadronDphiNoSSEta21,fHadronDphiNoSSEta22,fHadronDphiNoSSEta23,fHadronDphiNoSSEta24); | |
979 | fPiPtNoSS->Fill(pt); | |
980 | } | |
981 | ||
982 | //Electron id with TPC | |
983 | if(fTPCnSigma < fTPCnsigEleMin || fTPCnSigma > fTPCnsigEleMax) continue; | |
984 | fEovPWoSS->Fill(pt,fEovP); | |
985 | // fElecPhiTPCEovP->Fill(track->Phi()); | |
986 | ||
987 | //Electron id with shower shape | |
988 | if(cluster->GetM20()< fM20CutMin || cluster->GetM20()> fM20CutMax || cluster->GetM02()< fM02CutMin || cluster->GetM02()> fM02CutMax || cluster->GetDispersion()> fDispCutMax) continue; | |
989 | fEovPWSS->Fill(pt,fEovP); | |
990 | ||
991 | //Electron id with E/p | |
992 | if(fEovP < fEovPMin || fEovP > fEovPMax) continue; | |
993 | ||
994 | fTrkEovPAft->Fill(pt,fEovP); | |
995 | fElecPhi->Fill(track->Phi()); | |
996 | fElecPhiPt->Fill(track->Phi(),track->Pt()); | |
997 | if (track->Eta() >0 && track->Eta() <0.7) fElecPhiTPChalf->Fill(track->Phi()); | |
998 | ||
999 | HadronInfo(iTracks); | |
1000 | ||
1001 | Bool_t fFlagPhotonicElec = kFALSE; | |
1002 | // select photonic electron | |
1003 | SelectPhotonicElectron(iTracks,track,fFlagPhotonicElec); | |
20791315 | 1004 | SelectPhotonicElectronDCAMet(iTracks,track,fpidResponse); |
a9123929 | 1005 | |
1006 | //Inclusive electron-hadron correlation | |
1007 | ElectronHadCorrel(iTracks, track, fInclusiveElecDphi, fInclusiveElecDphi1,fInclusiveElecDphi2,fInclusiveElecDphi3,fInclusiveElecDphi4); | |
1008 | fInclusiveElecPt->Fill(pt); | |
3db57d71 | 1009 | MixedEvent(track,fMixedIncElecDphi, fMixedIncElecDphi1,fMixedIncElecDphi2, fMixedIncElecDphi3, fMixedIncElecDphi4); |
a9123929 | 1010 | |
1011 | //Dphi in Eta bins | |
1012 | ElectronHadCorrelEtaBins(iTracks, track, fInclusiveElecDphiEta1, fInclusiveElecDphiEta11,fInclusiveElecDphiEta12,fInclusiveElecDphiEta13,fInclusiveElecDphiEta14,fInclusiveElecDphiEta2, fInclusiveElecDphiEta21,fInclusiveElecDphiEta22,fInclusiveElecDphiEta23,fInclusiveElecDphiEta24); | |
1013 | ||
1014 | //Inclusive electron-hadron correlation far eta side | |
1015 | ElectronHadCorrelEtaFarSide(iTracks, track, fInclusiveElecDphiEtaFS, fInclusiveElecDphiEtaFS1,fInclusiveElecDphiEtaFS2,fInclusiveElecDphiEtaFS3,fInclusiveElecDphiEtaFS4); | |
3db57d71 | 1016 | MixedEvent(track,fMixedIncElecDphi, fMixedIncElecDphi1,fMixedIncElecDphi2, fMixedIncElecDphi3, fMixedIncElecDphi4); |
a9123929 | 1017 | |
1018 | // photonic electron | |
1019 | if(fFlagPhotonicElec){ | |
1020 | //Electron hadron correlation | |
1021 | ElectronHadCorrel(iTracks, track, fPhotElecDphi,fPhotElecDphi1,fPhotElecDphi2,fPhotElecDphi3,fPhotElecDphi4); | |
1022 | fPhotoElecPt->Fill(pt); | |
3db57d71 | 1023 | MixedEvent(track,fMixedPhotElecDphi, fMixedPhotElecDphi1,fMixedPhotElecDphi2, fMixedPhotElecDphi3, fMixedPhotElecDphi4); |
a9123929 | 1024 | |
1025 | //Dphi in Eta bins | |
1026 | ElectronHadCorrelEtaBins(iTracks, track, fPhotElecDphiEta1,fPhotElecDphiEta11,fPhotElecDphiEta12,fPhotElecDphiEta13,fPhotElecDphiEta14,fPhotElecDphiEta2,fPhotElecDphiEta21,fPhotElecDphiEta22,fPhotElecDphiEta23,fPhotElecDphiEta24); | |
1027 | } | |
1028 | ||
1029 | // Semi inclusive electron | |
1030 | if(!fFlagPhotonicElec){ | |
1031 | //Electron hadron correlation | |
1032 | ElectronHadCorrel(iTracks, track, fSemiIncElecDphi, fSemiIncElecDphi1,fSemiIncElecDphi2,fSemiIncElecDphi3,fSemiIncElecDphi4); | |
1033 | fSemiInclElecPt->Fill(pt); | |
3db57d71 | 1034 | MixedEvent(track,fMixedSemiIncElecDphi,fMixedSemiIncElecDphi1,fMixedSemiIncElecDphi2, fMixedSemiIncElecDphi3, fMixedSemiIncElecDphi4); |
a9123929 | 1035 | |
1036 | //Dphi in Eta bins | |
1037 | ElectronHadCorrelEtaBins(iTracks, track, fSemiIncElecDphiEta1, fSemiIncElecDphiEta11,fSemiIncElecDphiEta12,fSemiIncElecDphiEta13,fSemiIncElecDphiEta14,fSemiIncElecDphiEta2, fSemiIncElecDphiEta21,fSemiIncElecDphiEta22,fSemiIncElecDphiEta23,fSemiIncElecDphiEta24); | |
1038 | } | |
1039 | } | |
1040 | /* //EMC clusters | |
1041 | Int_t clsNo = fVevent->GetNumberOfCaloClusters(); | |
1042 | fNClusv1->Fill(clsNo); | |
1043 | for(Int_t iclus=0; iclus<clsNo ; iclus++){ | |
1044 | AliVCluster* clus = fVevent->GetCaloCluster(iclus); | |
1045 | if(!clus->IsEMCAL()) continue; | |
1046 | fNCellv1->Fill(clus->GetNCells()); | |
1047 | fClsEv1->Fill(clus->E()); | |
1048 | } | |
1049 | */ | |
3db57d71 | 1050 | |
a9123929 | 1051 | TObjArray* tracksClone = CloneAndReduceTrackList(); |
1052 | tracksClone->SetOwner(); | |
1053 | pool->UpdatePool(tracksClone); | |
3db57d71 | 1054 | |
a9123929 | 1055 | PostData(1, fOutputList); |
57189f04 | 1056 | } |
1057 | //_________________________________________ | |
1058 | void AliAnalysisTaskElecHadronCorrel::UserCreateOutputObjects() | |
1059 | { | |
1060 | //Create histograms | |
987053ce | 1061 | |
1062 | AliDebug(3, "Creating Output Objects"); | |
1063 | // Automatic determination of the analysis mode | |
1064 | AliVEventHandler *inputHandler = dynamic_cast<AliVEventHandler *>(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()); | |
1065 | if(!TString(inputHandler->IsA()->GetName()).CompareTo("AliAODInputHandler")){ | |
1066 | SetAODAnalysis(); | |
1067 | } else { | |
1068 | SetESDAnalysis(); | |
1069 | } | |
1070 | printf("Analysis Mode: %s Analysis\n", IsAODanalysis() ? "AOD" : "ESD"); | |
520a29a8 | 1071 | |
57189f04 | 1072 | //--------Initialize PID |
1073 | fPID->SetHasMCData(kFALSE); | |
1074 | if(!fPID->GetNumberOfPIDdetectors()) | |
520a29a8 | 1075 | { |
1076 | fPID->AddDetector("TPC", 0); | |
1077 | fPID->AddDetector("EMCAL", 1); | |
1078 | } | |
1079 | ||
57189f04 | 1080 | fPID->SortDetectors(); |
1081 | fPIDqa = new AliHFEpidQAmanager(); | |
1082 | fPIDqa->Initialize(fPID); | |
520a29a8 | 1083 | |
57189f04 | 1084 | //--------Initialize correction Framework and Cuts |
1085 | fCFM = new AliCFManager; | |
1086 | const Int_t kNcutSteps = AliHFEcuts::kNcutStepsMCTrack + AliHFEcuts::kNcutStepsRecTrack + AliHFEcuts::kNcutStepsDETrack; | |
1087 | fCFM->SetNStepParticle(kNcutSteps); | |
1088 | for(Int_t istep = 0; istep < kNcutSteps; istep++) | |
1089 | fCFM->SetParticleCutsList(istep, NULL); | |
520a29a8 | 1090 | |
57189f04 | 1091 | if(!fCuts){ |
1092 | AliWarning("Cuts not available. Default cuts will be used"); | |
1093 | fCuts = new AliHFEcuts; | |
1094 | fCuts->CreateStandardCuts(); | |
1095 | } | |
987053ce | 1096 | |
1097 | if(IsAODanalysis()) fCuts->SetAOD(); | |
57189f04 | 1098 | fCuts->Initialize(fCFM); |
3db57d71 | 1099 | |
4e01b68c | 1100 | //Mixed event initialising |
1101 | Int_t trackDepth = 2000; | |
1102 | Int_t poolsize = 1000; | |
1103 | ||
1104 | Int_t nCentralityBins = 5; | |
3836f884 | 1105 | Int_t nCentralityBinsSC = 6; |
487ae960 | 1106 | Int_t nZvtxBins = 4; |
a9123929 | 1107 | Double_t CentralityBins[6]; |
1108 | Double_t CentralityBinsSC[7]; | |
1109 | Double_t vertexBins[5]; | |
1110 | if(fTriggerCentral) | |
1111 | { | |
1112 | CentralityBins[0] = 0; | |
1113 | CentralityBins[1] = 2; | |
1114 | CentralityBins[2] = 4; | |
1115 | CentralityBins[3] = 6; | |
1116 | CentralityBins[4] = 8; | |
1117 | CentralityBins[5] = 10; | |
1118 | ||
1119 | vertexBins[0] = -10; | |
1120 | vertexBins[1] = -5; | |
1121 | vertexBins[2] = 0; | |
1122 | vertexBins[3] = 5; | |
1123 | vertexBins[4] = 10; | |
1124 | fPoolMgr = new AliEventPoolManager(poolsize, trackDepth, nCentralityBins, (Double_t*) CentralityBins, nZvtxBins, (Double_t*) vertexBins); | |
1125 | } | |
1126 | ||
a9123929 | 1127 | if(!fTriggerCentral) |
1128 | { | |
1129 | CentralityBinsSC[0] = 20; | |
1130 | CentralityBinsSC[1] = 25; | |
1131 | CentralityBinsSC[2] = 30; | |
1132 | CentralityBinsSC[3] = 35; | |
1133 | CentralityBinsSC[4] = 40; | |
1134 | CentralityBinsSC[5] = 45; | |
1135 | CentralityBinsSC[6] = 50; | |
1136 | ||
1137 | vertexBins[0] = -10; | |
1138 | vertexBins[1] = -5; | |
1139 | vertexBins[2] = 0; | |
1140 | vertexBins[3] = 5; | |
1141 | vertexBins[4] = 10; | |
3836f884 | 1142 | fPoolMgr = new AliEventPoolManager(poolsize, trackDepth, nCentralityBinsSC, (Double_t*) CentralityBinsSC, nZvtxBins, (Double_t*) vertexBins); |
a9123929 | 1143 | } |
3db57d71 | 1144 | |
4e01b68c | 1145 | |
57189f04 | 1146 | //---------Output Tlist |
1147 | fOutputList = new TList(); | |
1148 | fOutputList->SetOwner(); | |
1149 | fOutputList->Add(fPIDqa->MakeList("PIDQA")); | |
520a29a8 | 1150 | |
57189f04 | 1151 | fNoEvents = new TH1F("fNoEvents","",1,0,1) ; |
1152 | fOutputList->Add(fNoEvents); | |
520a29a8 | 1153 | |
1154 | fcentrality = new TH1F("fcentrality","centrality", 100,0,100); | |
1155 | fOutputList->Add(fcentrality); | |
987053ce | 1156 | |
bf697dda | 1157 | // fTrkpt = new TH1F("fTrkpt","track pt",1000,0,50); |
1158 | // fOutputList->Add(fTrkpt); | |
1159 | ||
1160 | // fTrackPtBefTrkCuts = new TH1F("fTrackPtBefTrkCuts","track pt before track cuts",1000,0,50); | |
1161 | // fOutputList->Add(fTrackPtBefTrkCuts); | |
520a29a8 | 1162 | |
bf697dda | 1163 | // fTrackPtAftTrkCuts = new TH1F("fTrackPtAftTrkCuts","track pt after track cuts",1000,0,50); |
1164 | // fOutputList->Add(fTrackPtAftTrkCuts); | |
520a29a8 | 1165 | |
7e02de85 | 1166 | fTPCnsigma = new TH2F("fTPCnsigma", "TPC - n sigma",1000,0,50,200,-10,10); |
1167 | fOutputList->Add(fTPCnsigma); | |
987053ce | 1168 | |
94b12681 | 1169 | fTrkEovPAft = new TH2F("fTrkEovPAft","track E/p after HFE pid",1000,0,50,100,0,2); |
1170 | fOutputList->Add(fTrkEovPAft); | |
520a29a8 | 1171 | |
bf697dda | 1172 | // fTrkEovPBefHad = new TH2F("fTrkEovPBefHad","track E/p for TPCnsig < 3.5",1000,0,50,100,0,2); |
1173 | // fOutputList->Add(fTrkEovPBefHad); | |
520a29a8 | 1174 | |
bf697dda | 1175 | // fdEdxBef = new TH2F("fdEdxBef","track dEdx vs p before HFE pid",1000,0,50,150,0,150); |
1176 | // fOutputList->Add(fdEdxBef); | |
520a29a8 | 1177 | |
520a29a8 | 1178 | fElecPhi = new TH1F("fElecPhi", "Electron phi",1000,0,6.28); |
1179 | fOutputList->Add(fElecPhi); | |
1180 | ||
3c56855b | 1181 | fElecPhiPt = new TH2F("fElecPhiPt", "Electron phi vs pt; Electron phi; pt (GeV/c)",1000,0,6.28,1000,0,100); |
1182 | fOutputList->Add(fElecPhiPt); | |
1183 | ||
1184 | fElecPhiTPChalf = new TH1F("fElecPhiTPChalf", "Electron phi for 0<eta<0.7",1000,0,6.28); | |
1185 | fOutputList->Add(fElecPhiTPChalf); | |
1186 | ||
bf697dda | 1187 | /* fElecPhiTPC = new TH1F("fElecPhiTPC", "Electron phi after TPC cut",1000,0,6.28); |
1188 | fOutputList->Add(fElecPhiTPC); | |
4e01b68c | 1189 | |
bf697dda | 1190 | fElecPhiTPCEovP = new TH1F("fElecPhiTPCEovP", "Electron phi after TPC and E/p cut",1000,0,6.28); |
1191 | fOutputList->Add(fElecPhiTPCEovP); | |
1192 | */ | |
520a29a8 | 1193 | fHadronPhi = new TH1F("fHadronPhi", "Hadron phi",1000,0,6.28); |
1194 | fOutputList->Add(fHadronPhi); | |
1195 | ||
3c56855b | 1196 | fHadronPhiTPChalf = new TH1F("fHadronPhiTPChalf", "Hadron phi for 0<eta<0.9",1000,0,6.28); |
1197 | fOutputList->Add(fHadronPhiTPChalf); | |
1198 | ||
1199 | fHadronPhiPt = new TH2F("fHadronPhiPt", "Hadron phi vs pt; hadron phi; pt (GeV/c)",1000,0,6.28,1000,0,100); | |
1200 | fOutputList->Add(fHadronPhiPt); | |
1201 | ||
bf697dda | 1202 | /* |
1203 | fTrackHFEcuts = new TH1F("fTrackHFEcuts","Track phi for HFE cuts",1000,0,6.28); | |
1204 | fOutputList->Add(fTrackHFEcuts); | |
4e01b68c | 1205 | |
bf697dda | 1206 | fTrakPhiSPD1 = new TH1F("fTrakPhiSPD1","Track phi for hit in SPD layer 1",1000,0,6.28); |
1207 | fOutputList->Add(fTrakPhiSPD1); | |
4e01b68c | 1208 | |
bf697dda | 1209 | fTrakPhiSPD2 = new TH1F("fTrakPhiSPD2","Track phi for hit in SPD layer 2",1000,0,6.28); |
1210 | fOutputList->Add(fTrakPhiSPD2); | |
4e01b68c | 1211 | |
bf697dda | 1212 | fTrakPhiSPDOr = new TH1F("fTrakPhiSPDOr","Track phi for hit in any SPD layer",1000,0,6.28); |
1213 | fOutputList->Add(fTrakPhiSPDOr); | |
4e01b68c | 1214 | |
bf697dda | 1215 | fTrakPhiSPDAnd = new TH1F("fTrakPhiSPDAnd","Track phi for hit in both SPD layer",1000,0,6.28); |
1216 | fOutputList->Add(fTrakPhiSPDAnd); | |
4e01b68c | 1217 | |
bf697dda | 1218 | fTrackHFEcutsITS = new TH1F("fTrackHFEcutsITS","Track phi for HFE cuts + ITS HFE cuts",1000,0,6.28); |
1219 | fOutputList->Add(fTrackHFEcutsITS); | |
1220 | */ | |
b7039672 | 1221 | fSemiIncElecDphi = new TH2F("fSemiIncElecDphi", "Semi Inclusive elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
57189f04 | 1222 | fOutputList->Add(fSemiIncElecDphi); |
520a29a8 | 1223 | |
b7039672 | 1224 | fSemiIncElecDphi1 = new TH2F("fSemiIncElecDphi1", "Semi Inclusive elec-had Dphi correlation for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
487ae960 | 1225 | fOutputList->Add(fSemiIncElecDphi1); |
1226 | ||
b7039672 | 1227 | fSemiIncElecDphi2 = new TH2F("fSemiIncElecDphi2", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
487ae960 | 1228 | fOutputList->Add(fSemiIncElecDphi2); |
1229 | ||
b7039672 | 1230 | fSemiIncElecDphi3 = new TH2F("fSemiIncElecDphi3", "Semi Inclusive elec-had Dphi correlation for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
487ae960 | 1231 | fOutputList->Add(fSemiIncElecDphi3); |
1232 | ||
b7039672 | 1233 | fSemiIncElecDphi4 = new TH2F("fSemiIncElecDphi4", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
487ae960 | 1234 | fOutputList->Add(fSemiIncElecDphi4); |
1235 | ||
b7039672 | 1236 | fPhotElecDphi = new TH2F("fPhotElecDphi", "Photon elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
57189f04 | 1237 | fOutputList->Add(fPhotElecDphi); |
520a29a8 | 1238 | |
b7039672 | 1239 | fPhotElecDphi1 = new TH2F("fPhotElecDphi1", "Photon elec-had Dphi correlation for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
487ae960 | 1240 | fOutputList->Add(fPhotElecDphi1); |
1241 | ||
b7039672 | 1242 | fPhotElecDphi2 = new TH2F("fPhotElecDphi2", "Photon elec-had Dphi correlation for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
487ae960 | 1243 | fOutputList->Add(fPhotElecDphi2); |
1244 | ||
b7039672 | 1245 | fPhotElecDphi3 = new TH2F("fPhotElecDphi3", "Photon elec-had Dphi correlation for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
487ae960 | 1246 | fOutputList->Add(fPhotElecDphi3); |
1247 | ||
b7039672 | 1248 | fPhotElecDphi4 = new TH2F("fPhotElecDphi4", "Photon elec-had Dphi correlation for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
487ae960 | 1249 | fOutputList->Add(fPhotElecDphi4); |
1250 | ||
b7039672 | 1251 | fInclusiveElecDphi = new TH2F("fInclusiveElecDphi", "Inclusive elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
57189f04 | 1252 | fOutputList->Add(fInclusiveElecDphi); |
520a29a8 | 1253 | |
b7039672 | 1254 | fInclusiveElecDphi1 = new TH2F("fInclusiveElecDphi1", "Inclusive elec-had Dphi correlation for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
487ae960 | 1255 | fOutputList->Add(fInclusiveElecDphi1); |
1256 | ||
b7039672 | 1257 | fInclusiveElecDphi2 = new TH2F("fInclusiveElecDphi2", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
487ae960 | 1258 | fOutputList->Add(fInclusiveElecDphi2); |
1259 | ||
b7039672 | 1260 | fInclusiveElecDphi3 = new TH2F("fInclusiveElecDphi3", "Inclusive elec-had Dphi correlation for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
487ae960 | 1261 | fOutputList->Add(fInclusiveElecDphi3); |
1262 | ||
b7039672 | 1263 | fInclusiveElecDphi4 = new TH2F("fInclusiveElecDphi4", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
487ae960 | 1264 | fOutputList->Add(fInclusiveElecDphi4); |
1265 | ||
b7039672 | 1266 | fInclusiveElecDphiEtaFS = new TH2F("fInclusiveElecDphiEtaFS", "Inclusive elec-had Dphi correlation (hadron 1<eta<1.6)",200,0,20,100,-1.57,4.71); |
33369c14 | 1267 | fOutputList->Add(fInclusiveElecDphiEtaFS); |
1268 | ||
b7039672 | 1269 | fInclusiveElecDphiEtaFS1 = new TH2F("fInclusiveElecDphiEtaFS1", "Inclusive elec-had Dphi correlation for 2<pt^{asso}<4 (hadron 1<eta<1.6)",200,0,20,100,-1.57,4.71); |
33369c14 | 1270 | fOutputList->Add(fInclusiveElecDphiEtaFS1); |
1271 | ||
b7039672 | 1272 | fInclusiveElecDphiEtaFS2 = new TH2F("fInclusiveElecDphiEtaFS2", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<6 (hadron 1<eta<1.6)",200,0,20,100,-1.57,4.71); |
33369c14 | 1273 | fOutputList->Add(fInclusiveElecDphiEtaFS2); |
1274 | ||
b7039672 | 1275 | fInclusiveElecDphiEtaFS3 = new TH2F("fInclusiveElecDphiEtaFS3", "Inclusive elec-had Dphi correlation for 6<pt^{asso}<8 (hadron 1<eta<1.6)",200,0,20,100,-1.57,4.71); |
33369c14 | 1276 | fOutputList->Add(fInclusiveElecDphiEtaFS3); |
1277 | ||
b7039672 | 1278 | fInclusiveElecDphiEtaFS4 = new TH2F("fInclusiveElecDphiEtaFS4", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<10 (hadron 1<eta<1.6)",200,0,20,100,-1.57,4.71); |
33369c14 | 1279 | fOutputList->Add(fInclusiveElecDphiEtaFS4); |
1280 | ||
b7039672 | 1281 | fDphiULSMassLow = new TH2F("fDphiULSMassLow", "e-h Dphi ULS, mass<cut",200,0,20,100,-1.57,4.71); |
57189f04 | 1282 | fOutputList->Add(fDphiULSMassLow); |
b7039672 | 1283 | fDphiULSMassLow1 = new TH2F("fDphiULSMassLow1", "e-h Dphi ULS, mass<cut for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
487ae960 | 1284 | fOutputList->Add(fDphiULSMassLow1); |
b7039672 | 1285 | fDphiULSMassLow2 = new TH2F("fDphiULSMassLow2", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
487ae960 | 1286 | fOutputList->Add(fDphiULSMassLow2); |
b7039672 | 1287 | fDphiULSMassLow3 = new TH2F("fDphiULSMassLow3", "e-h Dphi ULS, mass<cut for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
487ae960 | 1288 | fOutputList->Add(fDphiULSMassLow3); |
b7039672 | 1289 | fDphiULSMassLow4 = new TH2F("fDphiULSMassLow4", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
487ae960 | 1290 | fOutputList->Add(fDphiULSMassLow4); |
b7039672 | 1291 | fDphiLSMassLow = new TH2F("fDphiLSMassLow", "e-h Dphi LS, mass<cut",200,0,20,100,-1.57,4.71); |
57189f04 | 1292 | fOutputList->Add(fDphiLSMassLow); |
b7039672 | 1293 | fDphiLSMassLow1 = new TH2F("fDphiLSMassLow1", "e-h Dphi LS, mass<cut for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
487ae960 | 1294 | fOutputList->Add(fDphiLSMassLow1); |
b7039672 | 1295 | fDphiLSMassLow2 = new TH2F("fDphiLSMassLow2", "e-h Dphi LS, mass<cut for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
487ae960 | 1296 | fOutputList->Add(fDphiLSMassLow2); |
b7039672 | 1297 | fDphiLSMassLow3 = new TH2F("fDphiLSMassLow3", "e-h Dphi LS, mass<cut for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
487ae960 | 1298 | fOutputList->Add(fDphiLSMassLow3); |
b7039672 | 1299 | fDphiLSMassLow4 = new TH2F("fDphiLSMassLow4", "e-h Dphi LS, mass<cut for 8<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
487ae960 | 1300 | fOutputList->Add(fDphiLSMassLow4); |
b7039672 | 1301 | fDphiULSMassLowNoPartner = new TH2F("fDphiULSMassLowNoPartner", "e-h Dphi ULS with no partner, mass<mass cut,",200,0,20,100,-1.57,4.71); |
520a29a8 | 1302 | fOutputList->Add(fDphiULSMassLowNoPartner); |
b7039672 | 1303 | fDphiULSMassLowNoPartner1 = new TH2F("fDphiULSMassLowNoPartner1", "e-h Dphi ULS with no partner, mass<mass cut for 2<pt^{asso}<4,",200,0,20,100,-1.57,4.71); |
487ae960 | 1304 | fOutputList->Add(fDphiULSMassLowNoPartner1); |
b7039672 | 1305 | fDphiULSMassLowNoPartner2 = new TH2F("fDphiULSMassLowNoPartner2", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<6,",200,0,20,100,-1.57,4.71); |
487ae960 | 1306 | fOutputList->Add(fDphiULSMassLowNoPartner2); |
b7039672 | 1307 | fDphiULSMassLowNoPartner3 = new TH2F("fDphiULSMassLowNoPartner3", "e-h Dphi ULS with no partner, mass<mass cut for 6<pt^{asso}<8,",200,0,20,100,-1.57,4.71); |
487ae960 | 1308 | fOutputList->Add(fDphiULSMassLowNoPartner3); |
b7039672 | 1309 | fDphiULSMassLowNoPartner4 = new TH2F("fDphiULSMassLowNoPartner4", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<10,",200,0,20,100,-1.57,4.71); |
487ae960 | 1310 | fOutputList->Add(fDphiULSMassLowNoPartner4); |
b7039672 | 1311 | fDphiLSMassLowNoPartner = new TH2F("fDphiLSMassLowNoPartner", "e-h Dphi LS with no partner, mass<mass cut",200,0,20,100,-1.57,4.71); |
520a29a8 | 1312 | fOutputList->Add(fDphiLSMassLowNoPartner); |
b7039672 | 1313 | fDphiLSMassLowNoPartner1 = new TH2F("fDphiLSMassLowNoPartner1", "e-h Dphi LS with no partner, mass<mass cut for 2<pt^{asso}<4,",200,0,20,100,-1.57,4.71); |
487ae960 | 1314 | fOutputList->Add(fDphiLSMassLowNoPartner1); |
b7039672 | 1315 | fDphiLSMassLowNoPartner2 = new TH2F("fDphiLSMassLowNoPartner2", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<6,",200,0,20,100,-1.57,4.71); |
487ae960 | 1316 | fOutputList->Add(fDphiLSMassLowNoPartner2); |
b7039672 | 1317 | fDphiLSMassLowNoPartner3 = new TH2F("fDphiLSMassLowNoPartner3", "e-h Dphi LS with no partner, mass<mass cut for 6<pt^{asso}<8,",200,0,20,100,-1.57,4.71); |
487ae960 | 1318 | fOutputList->Add(fDphiLSMassLowNoPartner3); |
b7039672 | 1319 | fDphiLSMassLowNoPartner4 = new TH2F("fDphiLSMassLowNoPartner4", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<10,",200,0,20,100,-1.57,4.71); |
487ae960 | 1320 | fOutputList->Add(fDphiLSMassLowNoPartner4); |
57189f04 | 1321 | fPhotoElecPt = new TH1F("fPhotoElecPt", "photonic electron pt",1000,0,100); |
1322 | fOutputList->Add(fPhotoElecPt); | |
57189f04 | 1323 | fSemiInclElecPt = new TH1F("fSemiInclElecPt", "Semi-inclusive electron pt",1000,0,100); |
1324 | fOutputList->Add(fSemiInclElecPt); | |
520a29a8 | 1325 | fInclusiveElecPt = new TH1F("fInclElecPt", "Inclusive electron pt",1000,0,100); |
1326 | fOutputList->Add(fInclusiveElecPt); | |
520a29a8 | 1327 | fULSElecPt = new TH1F("fULSElecPt", "ULS electron pt",1000,0,100); |
1328 | fOutputList->Add(fULSElecPt); | |
520a29a8 | 1329 | fLSElecPt = new TH1F("fLSElecPt", "LS electron pt",1000,0,100); |
1330 | fOutputList->Add(fLSElecPt); | |
987053ce | 1331 | |
b7039672 | 1332 | fDCAMetPhotElecDphi = new TH2F("fDCAMetPhotElecDphi", "Photon elec-had Dphi correlation DCA method",200,0,20,100,-1.57,4.71); |
20791315 | 1333 | fOutputList->Add(fDCAMetPhotElecDphi); |
b7039672 | 1334 | fDCAMetPhotElecDphi1 = new TH2F("fDCAMetPhotElecDphi1", "Photon elec-had Dphi correlation DCA method for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
20791315 | 1335 | fOutputList->Add(fDCAMetPhotElecDphi1); |
b7039672 | 1336 | fDCAMetPhotElecDphi2 = new TH2F("fDCAMetPhotElecDphi2", "Photon elec-had Dphi correlation DCA method for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
20791315 | 1337 | fOutputList->Add(fDCAMetPhotElecDphi2); |
b7039672 | 1338 | fDCAMetPhotElecDphi3 = new TH2F("fDCAMetPhotElecDphi3", "Photon elec-had Dphi correlation DCA method for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
20791315 | 1339 | fOutputList->Add(fDCAMetPhotElecDphi3); |
b7039672 | 1340 | fDCAMetPhotElecDphi4 = new TH2F("fDCAMetPhotElecDphi4", "Photon elec-had Dphi correlation DCA method for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
20791315 | 1341 | fOutputList->Add(fDCAMetPhotElecDphi4); |
b7039672 | 1342 | fDCAMetDphiULSMassLow = new TH2F("fDCAMetDphiULSMassLow", "e-h Dphi ULS, mass<cut",200,0,20,100,-1.57,4.71); |
20791315 | 1343 | fOutputList->Add(fDCAMetDphiULSMassLow); |
b7039672 | 1344 | fDCAMetDphiULSMassLow1 = new TH2F("fDCAMetDphiULSMassLow1", "e-h Dphi ULS, mass<cut for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
20791315 | 1345 | fOutputList->Add(fDCAMetDphiULSMassLow1); |
b7039672 | 1346 | fDCAMetDphiULSMassLow2 = new TH2F("fDCAMetDphiULSMassLow2", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
20791315 | 1347 | fOutputList->Add(fDCAMetDphiULSMassLow2); |
b7039672 | 1348 | fDCAMetDphiULSMassLow3 = new TH2F("fDCAMetDphiULSMassLow3", "e-h Dphi ULS, mass<cut for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
20791315 | 1349 | fOutputList->Add(fDCAMetDphiULSMassLow3); |
b7039672 | 1350 | fDCAMetDphiULSMassLow4 = new TH2F("fDCAMetDphiULSMassLow4", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
20791315 | 1351 | fOutputList->Add(fDCAMetDphiULSMassLow4); |
b7039672 | 1352 | fDCAMetDphiLSMassLow = new TH2F("fDCAMetDphiLSMassLow", "e-h Dphi LS, mass<cut",200,0,20,100,-1.57,4.71); |
20791315 | 1353 | fOutputList->Add(fDCAMetDphiLSMassLow); |
b7039672 | 1354 | fDCAMetDphiLSMassLow1 = new TH2F("fDCAMetDphiLSMassLow1", "e-h Dphi LS, mass<cut for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
20791315 | 1355 | fOutputList->Add(fDCAMetDphiLSMassLow1); |
b7039672 | 1356 | fDCAMetDphiLSMassLow2 = new TH2F("fDCAMetDphiLSMassLow2", "e-h Dphi LS, mass<cut for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
20791315 | 1357 | fOutputList->Add(fDCAMetDphiLSMassLow2); |
b7039672 | 1358 | fDCAMetDphiLSMassLow3 = new TH2F("fDCAMetDphiLSMassLow3", "e-h Dphi LS, mass<cut for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
20791315 | 1359 | fOutputList->Add(fDCAMetDphiLSMassLow3); |
b7039672 | 1360 | fDCAMetDphiLSMassLow4 = new TH2F("fDCAMetDphiLSMassLow4", "e-h Dphi LS, mass<cut for 8<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
20791315 | 1361 | fOutputList->Add(fDCAMetDphiLSMassLow4); |
b7039672 | 1362 | fDCAMetDphiULSMassLowNoPartner = new TH2F("fDCAMetDphiULSMassLowNoPartner", "e-h Dphi ULS with no partner, mass<mass cut,",200,0,20,100,-1.57,4.71); |
20791315 | 1363 | fOutputList->Add(fDCAMetDphiULSMassLowNoPartner); |
b7039672 | 1364 | fDCAMetDphiULSMassLowNoPartner1 = new TH2F("fDCAMetDphiULSMassLowNoPartner1", "e-h Dphi ULS with no partner, mass<mass cut for 2<pt^{asso}<4,",200,0,20,100,-1.57,4.71); |
20791315 | 1365 | fOutputList->Add(fDCAMetDphiULSMassLowNoPartner1); |
b7039672 | 1366 | fDCAMetDphiULSMassLowNoPartner2 = new TH2F("fDCAMetDphiULSMassLowNoPartner2", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<6,",200,0,20,100,-1.57,4.71); |
20791315 | 1367 | fOutputList->Add(fDCAMetDphiULSMassLowNoPartner2); |
b7039672 | 1368 | fDCAMetDphiULSMassLowNoPartner3 = new TH2F("fDCAMetDphiULSMassLowNoPartner3", "e-h Dphi ULS with no partner, mass<mass cut for 6<pt^{asso}<8,",200,0,20,100,-1.57,4.71); |
20791315 | 1369 | fOutputList->Add(fDCAMetDphiULSMassLowNoPartner3); |
b7039672 | 1370 | fDCAMetDphiULSMassLowNoPartner4 = new TH2F("fDCAMetDphiULSMassLowNoPartner4", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<10,",200,0,20,100,-1.57,4.71); |
20791315 | 1371 | fOutputList->Add(fDCAMetDphiULSMassLowNoPartner4); |
b7039672 | 1372 | fDCAMetDphiLSMassLowNoPartner = new TH2F("fDCAMetDphiLSMassLowNoPartner", "e-h Dphi LS with no partner, mass<mass cut",200,0,20,100,-1.57,4.71); |
20791315 | 1373 | fOutputList->Add(fDCAMetDphiLSMassLowNoPartner); |
b7039672 | 1374 | fDCAMetDphiLSMassLowNoPartner1 = new TH2F("fDCAMetDphiLSMassLowNoPartner1", "e-h Dphi LS with no partner, mass<mass cut for 2<pt^{asso}<4,",200,0,20,100,-1.57,4.71); |
20791315 | 1375 | fOutputList->Add(fDCAMetDphiLSMassLowNoPartner1); |
b7039672 | 1376 | fDCAMetDphiLSMassLowNoPartner2 = new TH2F("fDCAMetDphiLSMassLowNoPartner2", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<6,",200,0,20,100,-1.57,4.71); |
20791315 | 1377 | fOutputList->Add(fDCAMetDphiLSMassLowNoPartner2); |
b7039672 | 1378 | fDCAMetDphiLSMassLowNoPartner3 = new TH2F("fDCAMetDphiLSMassLowNoPartner3", "e-h Dphi LS with no partner, mass<mass cut for 6<pt^{asso}<8,",200,0,20,100,-1.57,4.71); |
20791315 | 1379 | fOutputList->Add(fDCAMetDphiLSMassLowNoPartner3); |
b7039672 | 1380 | fDCAMetDphiLSMassLowNoPartner4 = new TH2F("fDCAMetDphiLSMassLowNoPartner4", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<10,",200,0,20,100,-1.57,4.71); |
20791315 | 1381 | fOutputList->Add(fDCAMetDphiLSMassLowNoPartner4); |
1382 | fDCAMetPhotoElecPt = new TH1F("fDCAMetPhotoElecPt", "photonic electron pt",1000,0,100); | |
1383 | fOutputList->Add(fDCAMetPhotoElecPt); | |
1384 | fDCAMetULSElecPt = new TH1F("fDCAMetULSElecPt", "ULS electron pt",1000,0,100); | |
1385 | fOutputList->Add(fDCAMetULSElecPt); | |
1386 | fDCAMetLSElecPt = new TH1F("fDCAMetLSElecPt", "LS electron pt",1000,0,100); | |
1387 | fOutputList->Add(fDCAMetLSElecPt); | |
1388 | ||
b7039672 | 1389 | fSemiIncElecDphiEta1 = new TH2F("fSemiIncElecDphiEta1", "Semi Inclusive elec-had Dphi correlation for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1390 | fOutputList->Add(fSemiIncElecDphiEta1); |
b7039672 | 1391 | fSemiIncElecDphiEta11 = new TH2F("fSemiIncElecDphiEta11", "Semi Inclusive elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1392 | fOutputList->Add(fSemiIncElecDphiEta11); |
b7039672 | 1393 | fSemiIncElecDphiEta12 = new TH2F("fSemiIncElecDphiEta12", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1394 | fOutputList->Add(fSemiIncElecDphiEta12); |
b7039672 | 1395 | fSemiIncElecDphiEta13 = new TH2F("fSemiIncElecDphiEta13", "Semi Inclusive elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1396 | fOutputList->Add(fSemiIncElecDphiEta13); |
b7039672 | 1397 | fSemiIncElecDphiEta14 = new TH2F("fSemiIncElecDphiEta14", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1398 | fOutputList->Add(fSemiIncElecDphiEta14); |
b7039672 | 1399 | fPhotElecDphiEta1 = new TH2F("fPhotElecDphiEta1", "Photon elec-had Dphi correlation for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1400 | fOutputList->Add(fPhotElecDphiEta1); |
b7039672 | 1401 | fPhotElecDphiEta11 = new TH2F("fPhotElecDphiEta11", "Photon elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1402 | fOutputList->Add(fPhotElecDphiEta11); |
b7039672 | 1403 | fPhotElecDphiEta12 = new TH2F("fPhotElecDphiEta12", "Photon elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1404 | fOutputList->Add(fPhotElecDphiEta12); |
b7039672 | 1405 | fPhotElecDphiEta13 = new TH2F("fPhotElecDphiEta13", "Photon elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1406 | fOutputList->Add(fPhotElecDphiEta13); |
b7039672 | 1407 | fPhotElecDphiEta14 = new TH2F("fPhotElecDphiEta14", "Photon elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1408 | fOutputList->Add(fPhotElecDphiEta14); |
b7039672 | 1409 | fInclusiveElecDphiEta1 = new TH2F("fInclusiveElecDphiEta1", "Inclusive elec-had Dphi correlation for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1410 | fOutputList->Add(fInclusiveElecDphiEta1); |
b7039672 | 1411 | fInclusiveElecDphiEta11 = new TH2F("fInclusiveElecDphiEta11", "Inclusive elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1412 | fOutputList->Add(fInclusiveElecDphiEta11); |
b7039672 | 1413 | fInclusiveElecDphiEta12 = new TH2F("fInclusiveElecDphiEta12", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1414 | fOutputList->Add(fInclusiveElecDphiEta12); |
b7039672 | 1415 | fInclusiveElecDphiEta13 = new TH2F("fInclusiveElecDphiEta13", "Inclusive elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1416 | fOutputList->Add(fInclusiveElecDphiEta13); |
b7039672 | 1417 | fInclusiveElecDphiEta14 = new TH2F("fInclusiveElecDphiEta14", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1418 | fOutputList->Add(fInclusiveElecDphiEta14); |
b7039672 | 1419 | fDphiULSMassLowEta1 = new TH2F("fDphiULSMassLowEta1", "e-h Dphi ULS, mass<cut for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1420 | fOutputList->Add(fDphiULSMassLowEta1); |
b7039672 | 1421 | fDphiULSMassLowEta11 = new TH2F("fDphiULSMassLowEta11", "e-h Dphi ULS, mass<cut for 2<pt^{asso}<4 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1422 | fOutputList->Add(fDphiULSMassLowEta11); |
b7039672 | 1423 | fDphiULSMassLowEta12 = new TH2F("fDphiULSMassLowEta12", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<6 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1424 | fOutputList->Add(fDphiULSMassLowEta12); |
b7039672 | 1425 | fDphiULSMassLowEta13 = new TH2F("fDphiULSMassLowEta13", "e-h Dphi ULS, mass<cut for 6<pt^{asso}<8 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1426 | fOutputList->Add(fDphiULSMassLowEta13); |
b7039672 | 1427 | fDphiULSMassLowEta14 = new TH2F("fDphiULSMassLowEta14", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<10 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1428 | fOutputList->Add(fDphiULSMassLowEta14); |
b7039672 | 1429 | fDphiLSMassLowEta1 = new TH2F("fDphiLSMassLowEta1", "e-h Dphi LS, mass<cut for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1430 | fOutputList->Add(fDphiLSMassLowEta1); |
b7039672 | 1431 | fDphiLSMassLowEta11 = new TH2F("fDphiLSMassLowEta11", "e-h Dphi LS, mass<cut for 2<pt^{asso}<4 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1432 | fOutputList->Add(fDphiLSMassLowEta11); |
b7039672 | 1433 | fDphiLSMassLowEta12 = new TH2F("fDphiLSMassLowEta12", "e-h Dphi LS, mass<cut for 4<pt^{asso}<6 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1434 | fOutputList->Add(fDphiLSMassLowEta12); |
b7039672 | 1435 | fDphiLSMassLowEta13 = new TH2F("fDphiLSMassLowEta13", "e-h Dphi LS, mass<cut for 6<pt^{asso}<8 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1436 | fOutputList->Add(fDphiLSMassLowEta13); |
b7039672 | 1437 | fDphiLSMassLowEta14 = new TH2F("fDphiLSMassLowEta14", "e-h Dphi LS, mass<cut for 8<pt^{asso}<10 for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1438 | fOutputList->Add(fDphiLSMassLowEta14); |
b7039672 | 1439 | fDphiULSMassLowNoPartnerEta1 = new TH2F("fDphiULSMassLowNoPartnerEta1", "e-h Dphi ULS with no partner, mass<mass cut, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1440 | fOutputList->Add(fDphiULSMassLowNoPartnerEta1); |
b7039672 | 1441 | fDphiULSMassLowNoPartnerEta11 = new TH2F("fDphiULSMassLowNoPartnerEta11", "e-h Dphi ULS with no partner, mass<mass cut for 2<pt^{asso}<4, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1442 | fOutputList->Add(fDphiULSMassLowNoPartnerEta11); |
b7039672 | 1443 | fDphiULSMassLowNoPartnerEta12 = new TH2F("fDphiULSMassLowNoPartnerEta12", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<6, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1444 | fOutputList->Add(fDphiULSMassLowNoPartnerEta12); |
b7039672 | 1445 | fDphiULSMassLowNoPartnerEta13 = new TH2F("fDphiULSMassLowNoPartnerEta13", "e-h Dphi ULS with no partner, mass<mass cut for 6<pt^{asso}<8, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1446 | fOutputList->Add(fDphiULSMassLowNoPartnerEta13); |
b7039672 | 1447 | fDphiULSMassLowNoPartnerEta14 = new TH2F("fDphiULSMassLowNoPartnerEta14", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<10, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1448 | fOutputList->Add(fDphiULSMassLowNoPartnerEta14); |
b7039672 | 1449 | fDphiLSMassLowNoPartnerEta1 = new TH2F("fDphiLSMassLowNoPartnerEta1", "e-h Dphi LS with no partner, mass<mass cut, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1450 | fOutputList->Add(fDphiLSMassLowNoPartnerEta1); |
b7039672 | 1451 | fDphiLSMassLowNoPartnerEta11 = new TH2F("fDphiLSMassLowNoPartnerEta11", "e-h Dphi LS with no partner, mass<mass cut for 2<pt^{asso}<4, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1452 | fOutputList->Add(fDphiLSMassLowNoPartnerEta11); |
b7039672 | 1453 | fDphiLSMassLowNoPartnerEta12 = new TH2F("fDphiLSMassLowNoPartnerEta12", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<6, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1454 | fOutputList->Add(fDphiLSMassLowNoPartnerEta12); |
b7039672 | 1455 | fDphiLSMassLowNoPartnerEta13 = new TH2F("fDphiLSMassLowNoPartnerEta13", "e-h Dphi LS with no partner, mass<mass cut for 6<pt^{asso}<8, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1456 | fOutputList->Add(fDphiLSMassLowNoPartnerEta13); |
b7039672 | 1457 | fDphiLSMassLowNoPartnerEta14 = new TH2F("fDphiLSMassLowNoPartnerEta14", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<10, for |Eta < 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1458 | fOutputList->Add(fDphiLSMassLowNoPartnerEta14); |
1459 | ||
b7039672 | 1460 | fSemiIncElecDphiEta2 = new TH2F("fSemiIncElecDphiEta2", "Semi Inclusive elec-had Dphi correlation for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1461 | fOutputList->Add(fSemiIncElecDphiEta2); |
b7039672 | 1462 | fSemiIncElecDphiEta21 = new TH2F("fSemiIncElecDphiEta21", "Semi Inclusive elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1463 | fOutputList->Add(fSemiIncElecDphiEta21); |
b7039672 | 1464 | fSemiIncElecDphiEta22 = new TH2F("fSemiIncElecDphiEta22", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1465 | fOutputList->Add(fSemiIncElecDphiEta22); |
b7039672 | 1466 | fSemiIncElecDphiEta23 = new TH2F("fSemiIncElecDphiEta23", "Semi Inclusive elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1467 | fOutputList->Add(fSemiIncElecDphiEta23); |
b7039672 | 1468 | fSemiIncElecDphiEta24 = new TH2F("fSemiIncElecDphiEta24", "Semi Inclusive elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1469 | fOutputList->Add(fSemiIncElecDphiEta24); |
b7039672 | 1470 | fPhotElecDphiEta2 = new TH2F("fPhotElecDphiEta2", "Photon elec-had Dphi correlation for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1471 | fOutputList->Add(fPhotElecDphiEta2); |
b7039672 | 1472 | fPhotElecDphiEta21 = new TH2F("fPhotElecDphiEta21", "Photon elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1473 | fOutputList->Add(fPhotElecDphiEta21); |
b7039672 | 1474 | fPhotElecDphiEta22 = new TH2F("fPhotElecDphiEta22", "Photon elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1475 | fOutputList->Add(fPhotElecDphiEta22); |
b7039672 | 1476 | fPhotElecDphiEta23 = new TH2F("fPhotElecDphiEta23", "Photon elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1477 | fOutputList->Add(fPhotElecDphiEta23); |
b7039672 | 1478 | fPhotElecDphiEta24 = new TH2F("fPhotElecDphiEta24", "Photon elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1479 | fOutputList->Add(fPhotElecDphiEta24); |
b7039672 | 1480 | fInclusiveElecDphiEta2 = new TH2F("fInclusiveElecDphiEta2", "Inclusive elec-had Dphi correlation for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1481 | fOutputList->Add(fInclusiveElecDphiEta2); |
b7039672 | 1482 | fInclusiveElecDphiEta21 = new TH2F("fInclusiveElecDphiEta21", "Inclusive elec-had Dphi correlation for 2<pt^{asso}<4 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1483 | fOutputList->Add(fInclusiveElecDphiEta21); |
b7039672 | 1484 | fInclusiveElecDphiEta22 = new TH2F("fInclusiveElecDphiEta22", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<6 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1485 | fOutputList->Add(fInclusiveElecDphiEta22); |
b7039672 | 1486 | fInclusiveElecDphiEta23 = new TH2F("fInclusiveElecDphiEta23", "Inclusive elec-had Dphi correlation for 6<pt^{asso}<8 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1487 | fOutputList->Add(fInclusiveElecDphiEta23); |
b7039672 | 1488 | fInclusiveElecDphiEta24 = new TH2F("fInclusiveElecDphiEta24", "Inclusive elec-had Dphi correlation for 4<pt^{asso}<10 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1489 | fOutputList->Add(fInclusiveElecDphiEta24); |
b7039672 | 1490 | fDphiULSMassLowEta2 = new TH2F("fDphiULSMassLowEta2", "e-h Dphi ULS, mass<cut for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1491 | fOutputList->Add(fDphiULSMassLowEta2); |
b7039672 | 1492 | fDphiULSMassLowEta21 = new TH2F("fDphiULSMassLowEta21", "e-h Dphi ULS, mass<cut for 2<pt^{asso}<4 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1493 | fOutputList->Add(fDphiULSMassLowEta21); |
b7039672 | 1494 | fDphiULSMassLowEta22 = new TH2F("fDphiULSMassLowEta22", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<6 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1495 | fOutputList->Add(fDphiULSMassLowEta22); |
b7039672 | 1496 | fDphiULSMassLowEta23 = new TH2F("fDphiULSMassLowEta23", "e-h Dphi ULS, mass<cut for 6<pt^{asso}<8 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1497 | fOutputList->Add(fDphiULSMassLowEta23); |
b7039672 | 1498 | fDphiULSMassLowEta24 = new TH2F("fDphiULSMassLowEta24", "e-h Dphi ULS, mass<cut for 4<pt^{asso}<10 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1499 | fOutputList->Add(fDphiULSMassLowEta24); |
b7039672 | 1500 | fDphiLSMassLowEta2 = new TH2F("fDphiLSMassLowEta2", "e-h Dphi LS, mass<cut for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1501 | fOutputList->Add(fDphiLSMassLowEta2); |
b7039672 | 1502 | fDphiLSMassLowEta21 = new TH2F("fDphiLSMassLowEta21", "e-h Dphi LS, mass<cut for 2<pt^{asso}<4 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1503 | fOutputList->Add(fDphiLSMassLowEta21); |
b7039672 | 1504 | fDphiLSMassLowEta22 = new TH2F("fDphiLSMassLowEta22", "e-h Dphi LS, mass<cut for 4<pt^{asso}<6 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1505 | fOutputList->Add(fDphiLSMassLowEta22); |
b7039672 | 1506 | fDphiLSMassLowEta23 = new TH2F("fDphiLSMassLowEta23", "e-h Dphi LS, mass<cut for 6<pt^{asso}<8 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1507 | fOutputList->Add(fDphiLSMassLowEta23); |
b7039672 | 1508 | fDphiLSMassLowEta24 = new TH2F("fDphiLSMassLowEta24", "e-h Dphi LS, mass<cut for 4<pt^{asso}<10 for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1509 | fOutputList->Add(fDphiLSMassLowEta24); |
b7039672 | 1510 | fDphiULSMassLowNoPartnerEta2 = new TH2F("fDphiULSMassLowNoPartnerEta2", "e-h Dphi ULS with no partner, mass<mass cut, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1511 | fOutputList->Add(fDphiULSMassLowNoPartnerEta2); |
b7039672 | 1512 | fDphiULSMassLowNoPartnerEta21 = new TH2F("fDphiULSMassLowNoPartnerEta21", "e-h Dphi ULS with no partner, mass<mass cut for 2<pt^{asso}<4, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1513 | fOutputList->Add(fDphiULSMassLowNoPartnerEta21); |
b7039672 | 1514 | fDphiULSMassLowNoPartnerEta22 = new TH2F("fDphiULSMassLowNoPartnerEta22", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<6, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1515 | fOutputList->Add(fDphiULSMassLowNoPartnerEta22); |
b7039672 | 1516 | fDphiULSMassLowNoPartnerEta23 = new TH2F("fDphiULSMassLowNoPartnerEta23", "e-h Dphi ULS with no partner, mass<mass cut for 6<pt^{asso}<8, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1517 | fOutputList->Add(fDphiULSMassLowNoPartnerEta23); |
b7039672 | 1518 | fDphiULSMassLowNoPartnerEta24 = new TH2F("fDphiULSMassLowNoPartnerEta24", "e-h Dphi ULS with no partner, mass<mass cut for 4<pt^{asso}<10, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1519 | fOutputList->Add(fDphiULSMassLowNoPartnerEta24); |
b7039672 | 1520 | fDphiLSMassLowNoPartnerEta2 = new TH2F("fDphiLSMassLowNoPartnerEta2", "e-h Dphi LS with no partner, mass<mass cut, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1521 | fOutputList->Add(fDphiLSMassLowNoPartnerEta2); |
b7039672 | 1522 | fDphiLSMassLowNoPartnerEta21 = new TH2F("fDphiLSMassLowNoPartnerEta21", "e-h Dphi LS with no partner, mass<mass cut for 2<pt^{asso}<4, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1523 | fOutputList->Add(fDphiLSMassLowNoPartnerEta21); |
b7039672 | 1524 | fDphiLSMassLowNoPartnerEta22 = new TH2F("fDphiLSMassLowNoPartnerEta22", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<6, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1525 | fOutputList->Add(fDphiLSMassLowNoPartnerEta22); |
b7039672 | 1526 | fDphiLSMassLowNoPartnerEta23 = new TH2F("fDphiLSMassLowNoPartnerEta23", "e-h Dphi LS with no partner, mass<mass cut for 6<pt^{asso}<8, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1527 | fOutputList->Add(fDphiLSMassLowNoPartnerEta23); |
b7039672 | 1528 | fDphiLSMassLowNoPartnerEta24 = new TH2F("fDphiLSMassLowNoPartnerEta24", "e-h Dphi LS with no partner, mass<mass cut for 4<pt^{asso}<10, for |Eta > 0.8|",200,0,20,100,-1.57,4.71); |
bf697dda | 1529 | fOutputList->Add(fDphiLSMassLowNoPartnerEta24); |
1530 | ||
1531 | /* fNCellv1 = new TH1F("fNCellv1","Ncell in clus (v1); NCell; count",100,0,100) ; | |
1532 | fOutputList->Add(fNCellv1); | |
1533 | ||
1534 | fClsEv1 = new TH1F("fClsEv1", "Clus E(v1); Cls E; count",1000,0,100); | |
1535 | fOutputList->Add(fClsEv1); | |
1536 | ||
1537 | fNClusv1 = new TH1F("fNClusv1","Nclus in event (v1); NClus; count",500,0,500) ; | |
1538 | fOutputList->Add(fNClusv1); | |
1539 | */ | |
33369c14 | 1540 | fInvmassLS1 = new TH1F("fInvmassLS1", "Inv mass of LS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); |
520a29a8 | 1541 | fOutputList->Add(fInvmassLS1); |
1542 | ||
33369c14 | 1543 | fInvmassULS1 = new TH1F("fInvmassULS1", "Inv mass of ULS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); |
520a29a8 | 1544 | fOutputList->Add(fInvmassULS1); |
20791315 | 1545 | |
1546 | fDCAMetInvmassLS1 = new TH1F("fDCAMetInvmassLS1", "Inv mass of LS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); | |
1547 | fOutputList->Add(fDCAMetInvmassLS1); | |
1548 | ||
1549 | fDCAMetInvmassULS1 = new TH1F("fDCAMetInvmassULS1", "Inv mass of ULS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); | |
1550 | fOutputList->Add(fDCAMetInvmassULS1); | |
bf697dda | 1551 | /* |
1552 | fInvmassLS2 = new TH1F("fInvmassLS2", "Inv mass of LS (e,e) for pt^{e}>1; mass(GeV/c^2); counts;", 1000,0,1.0); | |
1553 | fOutputList->Add(fInvmassLS2); | |
487ae960 | 1554 | |
bf697dda | 1555 | fInvmassULS2 = new TH1F("fInvmassULS2", "Inv mass of ULS (e,e) for pt^{e}>1; mass(GeV/c^2); counts;", 1000,0,1.0); |
1556 | fOutputList->Add(fInvmassULS2); | |
487ae960 | 1557 | |
bf697dda | 1558 | fInvmassLS3 = new TH1F("fInvmassLS3", "Inv mass of LS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); |
1559 | fOutputList->Add(fInvmassLS3); | |
487ae960 | 1560 | |
bf697dda | 1561 | fInvmassULS3 = new TH1F("fInvmassULS3", "Inv mass of ULS (e,e) for pt^{e}>2; mass(GeV/c^2); counts;", 1000,0,1.0); |
1562 | fOutputList->Add(fInvmassULS3); | |
487ae960 | 1563 | |
bf697dda | 1564 | fInvmassLS4 = new TH1F("fInvmassLS4", "Inv mass of LS (e,e) for pt^{e}>3; mass(GeV/c^2); counts;", 1000,0,1.0); |
1565 | fOutputList->Add(fInvmassLS4); | |
487ae960 | 1566 | |
bf697dda | 1567 | fInvmassULS4 = new TH1F("fInvmassULS4", "Inv mass of ULS (e,e) for pt^{e}>3; mass(GeV/c^2); counts;", 1000,0,1.0); |
1568 | fOutputList->Add(fInvmassULS4); | |
487ae960 | 1569 | |
bf697dda | 1570 | fInvmassLS5 = new TH1F("fInvmassLS5", "Inv mass of LS (e,e) for pt^{e}>4; mass(GeV/c^2); counts;", 1000,0,1.0); |
1571 | fOutputList->Add(fInvmassLS5); | |
487ae960 | 1572 | |
bf697dda | 1573 | fInvmassULS5 = new TH1F("fInvmassULS5", "Inv mass of ULS (e,e) for pt^{e}>4; mass(GeV/c^2); counts;", 1000,0,1.0); |
1574 | fOutputList->Add(fInvmassULS5); | |
1575 | */ | |
3db57d71 | 1576 | fNoMixedEvents = new TH1F("fNoMixedEvents","",1,0,1) ; |
a9123929 | 1577 | fOutputList->Add(fNoMixedEvents); |
4e01b68c | 1578 | |
a9123929 | 1579 | fMixStat = new TH2F("fMixStat","no of events in pool vs Centrality;Nevent in pool;Centrality",200,0,200,5,0,10); |
1580 | fOutputList->Add(fMixStat); | |
4e01b68c | 1581 | |
a9123929 | 1582 | fMixStat1 = new TH2F("fMixStat1","no of events in pool vs zvtx;Nevents in pool;zvtx",200,0,200,4,-10,10); |
1583 | fOutputList->Add(fMixStat1); | |
4e01b68c | 1584 | |
b7039672 | 1585 | fMixedIncElecDphi = new TH2F("fMixedIncElecDphi", "Mixed event - Inclusive elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
a9123929 | 1586 | fOutputList->Add(fMixedIncElecDphi); |
4e01b68c | 1587 | |
b7039672 | 1588 | fMixedIncElecDphi1 = new TH2F("fMixedIncElecDphi1", "Mixed event - Inclusive elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
a9123929 | 1589 | fOutputList->Add(fMixedIncElecDphi1); |
3c56855b | 1590 | |
b7039672 | 1591 | fMixedIncElecDphi2 = new TH2F("fMixedIncElecDphi2", "Mixed event - Inclusive elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
a9123929 | 1592 | fOutputList->Add(fMixedIncElecDphi2); |
3c56855b | 1593 | |
b7039672 | 1594 | fMixedIncElecDphi3 = new TH2F("fMixedIncElecDphi3", "Mixed event - Inclusive elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
a9123929 | 1595 | fOutputList->Add(fMixedIncElecDphi3); |
4e01b68c | 1596 | |
b7039672 | 1597 | fMixedIncElecDphi4 = new TH2F("fMixedIncElecDphi4", "Mixed event - Inclusive elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
a9123929 | 1598 | fOutputList->Add(fMixedIncElecDphi4); |
3c56855b | 1599 | |
b7039672 | 1600 | fMixedSemiIncElecDphi = new TH2F("fMixedSemiIncElecDphi", "Mixed event - Semi Inclusive elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
a9123929 | 1601 | fOutputList->Add(fMixedSemiIncElecDphi); |
3c56855b | 1602 | |
b7039672 | 1603 | fMixedSemiIncElecDphi1 = new TH2F("fMixedSemiIncElecDphi1", "Mixed event - Semi Inclusive elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
a9123929 | 1604 | fOutputList->Add(fMixedSemiIncElecDphi1); |
4e01b68c | 1605 | |
b7039672 | 1606 | fMixedSemiIncElecDphi2 = new TH2F("fMixedSemiIncElecDphi2", "Mixed event - Semi Inclusive elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
a9123929 | 1607 | fOutputList->Add(fMixedSemiIncElecDphi2); |
3c56855b | 1608 | |
b7039672 | 1609 | fMixedSemiIncElecDphi3 = new TH2F("fMixedSemiIncElecDphi3", "Mixed event - Semi Inclusive elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
a9123929 | 1610 | fOutputList->Add(fMixedSemiIncElecDphi3); |
3c56855b | 1611 | |
b7039672 | 1612 | fMixedSemiIncElecDphi4 = new TH2F("fMixedSemiIncElecDphi4", "Mixed event - Semi Inclusive elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
a9123929 | 1613 | fOutputList->Add(fMixedSemiIncElecDphi4); |
4e01b68c | 1614 | |
b7039672 | 1615 | fMixedPhotElecDphi = new TH2F("fMixedPhotElecDphi", "Mixed event - Photo elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
a9123929 | 1616 | fOutputList->Add(fMixedPhotElecDphi); |
3c56855b | 1617 | |
b7039672 | 1618 | fMixedPhotElecDphi1 = new TH2F("fMixedPhotElecDphi1", "Mixed event - Photo elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
a9123929 | 1619 | fOutputList->Add(fMixedPhotElecDphi1); |
3c56855b | 1620 | |
b7039672 | 1621 | fMixedPhotElecDphi2 = new TH2F("fMixedPhotElecDphi2", "Mixed event - Photo elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
a9123929 | 1622 | fOutputList->Add(fMixedPhotElecDphi2); |
3c56855b | 1623 | |
b7039672 | 1624 | fMixedPhotElecDphi3 = new TH2F("fMixedPhotElecDphi3", "Mixed event - Photo elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
a9123929 | 1625 | fOutputList->Add(fMixedPhotElecDphi3); |
1626 | ||
b7039672 | 1627 | fMixedPhotElecDphi4 = new TH2F("fMixedPhotElecDphi4", "Mixed event - Photo elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
a9123929 | 1628 | fOutputList->Add(fMixedPhotElecDphi4); |
1629 | ||
b7039672 | 1630 | fMixedDphiULSMassLow = new TH2F("fMixedDphiULSMassLow", "Mixed event - ULS mass < cut elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
a9123929 | 1631 | fOutputList->Add(fMixedDphiULSMassLow); |
1632 | ||
b7039672 | 1633 | fMixedDphiULSMassLow1 = new TH2F("fMixedDphiULSMassLow1", "Mixed event - ULS mass < cut elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
a9123929 | 1634 | fOutputList->Add(fMixedDphiULSMassLow1); |
1635 | ||
b7039672 | 1636 | fMixedDphiULSMassLow2 = new TH2F("fMixedDphiULSMassLow2", "Mixed event - ULS mass < cut elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
a9123929 | 1637 | fOutputList->Add(fMixedDphiULSMassLow2); |
1638 | ||
b7039672 | 1639 | fMixedDphiULSMassLow3 = new TH2F("fMixedDphiULSMassLow3", "Mixed event - ULS mass < cut elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
a9123929 | 1640 | fOutputList->Add(fMixedDphiULSMassLow3); |
1641 | ||
b7039672 | 1642 | fMixedDphiULSMassLow4 = new TH2F("fMixedDphiULSMassLow4", "Mixed event - ULS mass < cut elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
a9123929 | 1643 | fOutputList->Add(fMixedDphiULSMassLow4); |
1644 | ||
b7039672 | 1645 | fMixedDphiLSMassLow = new TH2F("fMixedDphiLSMassLow", "Mixed event - LS mass < cut elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
a9123929 | 1646 | fOutputList->Add(fMixedDphiLSMassLow); |
1647 | ||
b7039672 | 1648 | fMixedDphiLSMassLow1 = new TH2F("fMixedDphiLSMassLow1", "Mixed event - LS mass < cut elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
a9123929 | 1649 | fOutputList->Add(fMixedDphiLSMassLow1); |
1650 | ||
b7039672 | 1651 | fMixedDphiLSMassLow2 = new TH2F("fMixedDphiLSMassLow2", "Mixed event - LS mass < cut elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
a9123929 | 1652 | fOutputList->Add(fMixedDphiLSMassLow2); |
1653 | ||
b7039672 | 1654 | fMixedDphiLSMassLow3 = new TH2F("fMixedDphiLSMassLow3", "Mixed event - LS mass < cut elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
a9123929 | 1655 | fOutputList->Add(fMixedDphiLSMassLow3); |
1656 | ||
b7039672 | 1657 | fMixedDphiLSMassLow4 = new TH2F("fMixedDphiLSMassLow4", "Mixed event - LS mass < cut elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
a9123929 | 1658 | fOutputList->Add(fMixedDphiLSMassLow4); |
3db57d71 | 1659 | |
b7039672 | 1660 | fDCAMetMixedPhotElecDphi = new TH2F("fDCAMetMixedPhotElecDphi", "Mixed event - Photo elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
20791315 | 1661 | fOutputList->Add(fDCAMetMixedPhotElecDphi); |
1662 | ||
b7039672 | 1663 | fDCAMetMixedPhotElecDphi1 = new TH2F("fDCAMetMixedPhotElecDphi1", "Mixed event - Photo elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
20791315 | 1664 | fOutputList->Add(fDCAMetMixedPhotElecDphi1); |
1665 | ||
b7039672 | 1666 | fDCAMetMixedPhotElecDphi2 = new TH2F("fDCAMetMixedPhotElecDphi2", "Mixed event - Photo elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
20791315 | 1667 | fOutputList->Add(fDCAMetMixedPhotElecDphi2); |
1668 | ||
b7039672 | 1669 | fDCAMetMixedPhotElecDphi3 = new TH2F("fDCAMetMixedPhotElecDphi3", "Mixed event - Photo elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
20791315 | 1670 | fOutputList->Add(fDCAMetMixedPhotElecDphi3); |
1671 | ||
b7039672 | 1672 | fDCAMetMixedPhotElecDphi4 = new TH2F("fDCAMetMixedPhotElecDphi4", "Mixed event - Photo elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
20791315 | 1673 | fOutputList->Add(fDCAMetMixedPhotElecDphi4); |
1674 | ||
b7039672 | 1675 | fDCAMetMixedDphiULSMassLow = new TH2F("fDCAMetMixedDphiULSMassLow", "Mixed event - ULS mass < cut elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
20791315 | 1676 | fOutputList->Add(fDCAMetMixedDphiULSMassLow); |
1677 | ||
b7039672 | 1678 | fDCAMetMixedDphiULSMassLow1 = new TH2F("fDCAMetMixedDphiULSMassLow1", "Mixed event - ULS mass < cut elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
20791315 | 1679 | fOutputList->Add(fDCAMetMixedDphiULSMassLow1); |
1680 | ||
b7039672 | 1681 | fDCAMetMixedDphiULSMassLow2 = new TH2F("fDCAMetMixedDphiULSMassLow2", "Mixed event - ULS mass < cut elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
20791315 | 1682 | fOutputList->Add(fDCAMetMixedDphiULSMassLow2); |
1683 | ||
b7039672 | 1684 | fDCAMetMixedDphiULSMassLow3 = new TH2F("fDCAMetMixedDphiULSMassLow3", "Mixed event - ULS mass < cut elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
20791315 | 1685 | fOutputList->Add(fDCAMetMixedDphiULSMassLow3); |
1686 | ||
b7039672 | 1687 | fDCAMetMixedDphiULSMassLow4 = new TH2F("fDCAMetMixedDphiULSMassLow4", "Mixed event - ULS mass < cut elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
20791315 | 1688 | fOutputList->Add(fDCAMetMixedDphiULSMassLow4); |
1689 | ||
b7039672 | 1690 | fDCAMetMixedDphiLSMassLow = new TH2F("fDCAMetMixedDphiLSMassLow", "Mixed event - LS mass < cut elec-had Dphi correlation",200,0,20,100,-1.57,4.71); |
20791315 | 1691 | fOutputList->Add(fDCAMetMixedDphiLSMassLow); |
1692 | ||
b7039672 | 1693 | fDCAMetMixedDphiLSMassLow1 = new TH2F("fDCAMetMixedDphiLSMassLow1", "Mixed event - LS mass < cut elec-had Dphi correlation 2<pt<4",200,0,20,100,-1.57,4.71); |
20791315 | 1694 | fOutputList->Add(fDCAMetMixedDphiLSMassLow1); |
1695 | ||
b7039672 | 1696 | fDCAMetMixedDphiLSMassLow2 = new TH2F("fDCAMetMixedDphiLSMassLow2", "Mixed event - LS mass < cut elec-had Dphi correlation 4<pt<6",200,0,20,100,-1.57,4.71); |
20791315 | 1697 | fOutputList->Add(fDCAMetMixedDphiLSMassLow2); |
1698 | ||
b7039672 | 1699 | fDCAMetMixedDphiLSMassLow3 = new TH2F("fDCAMetMixedDphiLSMassLow3", "Mixed event - LS mass < cut elec-had Dphi correlation 6<pt<8",200,0,20,100,-1.57,4.71); |
20791315 | 1700 | fOutputList->Add(fDCAMetMixedDphiLSMassLow3); |
1701 | ||
b7039672 | 1702 | fDCAMetMixedDphiLSMassLow4 = new TH2F("fDCAMetMixedDphiLSMassLow4", "Mixed event - LS mass < cut elec-had Dphi correlation 4<pt<10",200,0,20,100,-1.57,4.71); |
20791315 | 1703 | fOutputList->Add(fDCAMetMixedDphiLSMassLow4); |
1704 | ||
3c56855b | 1705 | fHadronPt = new TH1F("fHadronPt","hadron pt distribution",1000,0,100); |
1706 | fOutputList->Add(fHadronPt); | |
1707 | ||
1708 | fCentralityPass = new TH1F("fCentralityPass", "Centrality Pass", 101, -1, 100); | |
1709 | fOutputList->Add(fCentralityPass); | |
1710 | ||
1711 | fCentralityNoPass = new TH1F("fCentralityNoPass", "Centrality No Pass", 101, -1, 100); | |
1712 | fOutputList->Add(fCentralityNoPass); | |
1713 | ||
0841d73f | 1714 | fHadronDphi = new TH2F("fHadronDphi", "Hadron-had Dphi correlation",200,0,20,64,-1.57,4.71); |
1715 | fOutputList->Add(fHadronDphi); | |
1716 | ||
a9123929 | 1717 | fHadronDphi1 = new TH2F("fHadronDphi1", "Hadron-had Dphi correlation for 2<pt^{asso}<4",200,0,20,64,-1.57,4.71); |
0841d73f | 1718 | fOutputList->Add(fHadronDphi1); |
1719 | ||
a9123929 | 1720 | fHadronDphi2 = new TH2F("fHadronDphi2", "Hadron-had Dphi correlation for 4<pt^{asso}<6",200,0,20,64,-1.57,4.71); |
0841d73f | 1721 | fOutputList->Add(fHadronDphi2); |
1722 | ||
a9123929 | 1723 | fHadronDphi3 = new TH2F("fHadronDphi3", "Hadron-had Dphi correlation for 6<pt^{asso}<8",200,0,20,64,-1.57,4.71); |
0841d73f | 1724 | fOutputList->Add(fHadronDphi3); |
1725 | ||
a9123929 | 1726 | fHadronDphi4 = new TH2F("fHadronDphi4", "Hadron-had Dphi correlation for 4<pt^{asso}<10",200,0,20,64,-1.57,4.71); |
0841d73f | 1727 | fOutputList->Add(fHadronDphi4); |
1728 | ||
33369c14 | 1729 | fPiPt = new TH1F("fPiPt","Pi (-10 <TPC nsig < -3.5) pt distribution",1000,0,100); |
0841d73f | 1730 | fOutputList->Add(fPiPt); |
1731 | ||
b7039672 | 1732 | fHadronDphiNoSS = new TH2F("fHadronDphiNoSS", "Hadron-had Dphi correlation (NoSS cuts)",200,0,20,100,-1.57,4.71); |
33369c14 | 1733 | fOutputList->Add(fHadronDphiNoSS); |
1734 | ||
b7039672 | 1735 | fHadronDphiNoSS1 = new TH2F("fHadronDphiNoSS1", "Hadron-had Dphi correlation (NoSS cuts) for 2<pt^{asso}<4",200,0,20,100,-1.57,4.71); |
33369c14 | 1736 | fOutputList->Add(fHadronDphiNoSS1); |
1737 | ||
b7039672 | 1738 | fHadronDphiNoSS2 = new TH2F("fHadronDphiNoSS2", "Hadron-had Dphi correlation (NoSS cuts) for 4<pt^{asso}<6",200,0,20,100,-1.57,4.71); |
33369c14 | 1739 | fOutputList->Add(fHadronDphiNoSS2); |
1740 | ||
b7039672 | 1741 | fHadronDphiNoSS3 = new TH2F("fHadronDphiNoSS3", "Hadron-had Dphi correlation (NoSS cuts) for 6<pt^{asso}<8",200,0,20,100,-1.57,4.71); |
33369c14 | 1742 | fOutputList->Add(fHadronDphiNoSS3); |
1743 | ||
b7039672 | 1744 | fHadronDphiNoSS4 = new TH2F("fHadronDphiNoSS4", "Hadron-had Dphi correlation (NoSS cuts) for 4<pt^{asso}<10",200,0,20,100,-1.57,4.71); |
33369c14 | 1745 | fOutputList->Add(fHadronDphiNoSS4); |
1746 | ||
1747 | fPiPtNoSS = new TH1F("fPiPtNoSS","Pi (-10 <TPC nsig < -3.5) (NoSS cuts) pt distribution",1000,0,100); | |
1748 | fOutputList->Add(fPiPtNoSS); | |
1749 | ||
bf697dda | 1750 | fEovPWoSS = new TH2F("fEovPWoSS","E/p distribution without SS cuts",1000,0,50,100,0,2); |
1751 | fOutputList->Add(fEovPWoSS); | |
1752 | ||
1753 | fEovPWSS = new TH2F("fEovPWSS","E/p distribution with SS cuts",1000,0,50,100,0,2); | |
1754 | fOutputList->Add(fEovPWSS); | |
1755 | ||
1756 | fEovPHadWoSS = new TH2F("fEovPHadWoSS","E/p distribution for hadrons without SS cuts",1000,0,50,100,0,2); | |
1757 | fOutputList->Add(fEovPHadWoSS); | |
1758 | ||
1759 | fEovPHadWSS = new TH2F("fEovPHadWSS","E/p distribution for hadrons with SS cuts",1000,0,50,100,0,2); | |
1760 | fOutputList->Add(fEovPHadWSS); | |
1761 | ||
b7039672 | 1762 | fHadronDphiEta1 = new TH2F("fHadronDphiEta1", "Hadron-had Dphi correlation for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1763 | fOutputList->Add(fHadronDphiEta1); |
b7039672 | 1764 | fHadronDphiEta11 = new TH2F("fHadronDphiEta11", "Hadron-had Dphi correlation for 2<pt^{asso}<4 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1765 | fOutputList->Add(fHadronDphiEta11); |
b7039672 | 1766 | fHadronDphiEta12 = new TH2F("fHadronDphiEta12", "Hadron-had Dphi correlation for 4<pt^{asso}<6 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1767 | fOutputList->Add(fHadronDphiEta12); |
b7039672 | 1768 | fHadronDphiEta13 = new TH2F("fHadronDphiEta13", "Hadron-had Dphi correlation for 6<pt^{asso}<8 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1769 | fOutputList->Add(fHadronDphiEta13); |
b7039672 | 1770 | fHadronDphiEta14 = new TH2F("fHadronDphiEta14", "Hadron-had Dphi correlation for 8<pt^{asso}<10 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1771 | fOutputList->Add(fHadronDphiEta14); |
1772 | ||
b7039672 | 1773 | fHadronDphiNoSSEta1 = new TH2F("fHadronDphiNoSSEta1", "Hadron-had DphiNoSS correlation for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1774 | fOutputList->Add(fHadronDphiNoSSEta1); |
b7039672 | 1775 | fHadronDphiNoSSEta11 = new TH2F("fHadronDphiNoSSEta11", "Hadron-had DphiNoSS correlation for 2<pt^{asso}<4 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1776 | fOutputList->Add(fHadronDphiNoSSEta11); |
b7039672 | 1777 | fHadronDphiNoSSEta12 = new TH2F("fHadronDphiNoSSEta12", "Hadron-had DphiNoSS correlation for 4<pt^{asso}<6 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1778 | fOutputList->Add(fHadronDphiNoSSEta12); |
b7039672 | 1779 | fHadronDphiNoSSEta13 = new TH2F("fHadronDphiNoSSEta13", "Hadron-had DphiNoSS correlation for 6<pt^{asso}<8 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1780 | fOutputList->Add(fHadronDphiNoSSEta13); |
b7039672 | 1781 | fHadronDphiNoSSEta14 = new TH2F("fHadronDphiNoSSEta14", "Hadron-had DphiNoSS correlation for 4<pt^{asso}<10 for |Eta <1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1782 | fOutputList->Add(fHadronDphiNoSSEta14); |
1783 | ||
b7039672 | 1784 | fHadronDphiEta2 = new TH2F("fHadronDphiEta2", "Hadron-had Dphi correlation for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1785 | fOutputList->Add(fHadronDphiEta2); |
b7039672 | 1786 | fHadronDphiEta21 = new TH2F("fHadronDphiEta21", "Hadron-had Dphi correlation for 2<pt^{asso}<4 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1787 | fOutputList->Add(fHadronDphiEta21); |
b7039672 | 1788 | fHadronDphiEta22 = new TH2F("fHadronDphiEta22", "Hadron-had Dphi correlation for 4<pt^{asso}<6 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1789 | fOutputList->Add(fHadronDphiEta22); |
b7039672 | 1790 | fHadronDphiEta23 = new TH2F("fHadronDphiEta23", "Hadron-had Dphi correlation for 6<pt^{asso}<8 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1791 | fOutputList->Add(fHadronDphiEta23); |
b7039672 | 1792 | fHadronDphiEta24 = new TH2F("fHadronDphiEta24", "Hadron-had Dphi correlation for 4<pt^{asso}<10 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1793 | fOutputList->Add(fHadronDphiEta24); |
1794 | ||
b7039672 | 1795 | fHadronDphiNoSSEta2 = new TH2F("fHadronDphiNoSSEta2", "Hadron-had DphiNoSS correlation for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1796 | fOutputList->Add(fHadronDphiNoSSEta2); |
b7039672 | 1797 | fHadronDphiNoSSEta21 = new TH2F("fHadronDphiNoSSEta21", "Hadron-had DphiNoSS correlation for 2<pt^{asso}<4 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1798 | fOutputList->Add(fHadronDphiNoSSEta21); |
b7039672 | 1799 | fHadronDphiNoSSEta22 = new TH2F("fHadronDphiNoSSEta22", "Hadron-had DphiNoSS correlation for 4<pt^{asso}<6 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1800 | fOutputList->Add(fHadronDphiNoSSEta22); |
b7039672 | 1801 | fHadronDphiNoSSEta23 = new TH2F("fHadronDphiNoSSEta23", "Hadron-had DphiNoSS correlation for 6<pt^{asso}<8 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1802 | fOutputList->Add(fHadronDphiNoSSEta23); |
b7039672 | 1803 | fHadronDphiNoSSEta24 = new TH2F("fHadronDphiNoSSEta24", "Hadron-had DphiNoSS correlation for 4<pt^{asso}<10 for |Eta >1|",200,0,20,100,-1.57,4.71); |
0c884279 | 1804 | fOutputList->Add(fHadronDphiNoSSEta24); |
1805 | ||
1806 | ||
bf697dda | 1807 | // Int_t binsv1[6]={500,200,50,50,50,50}; //pt, TPCnsig, E/p, M20, M02, dispersion |
1808 | // Double_t xminv1[6]={0,-10,0,0,0,0}; | |
1809 | // Double_t xmaxv1[6]={25,10,2,2,2,2}; | |
1810 | // fSparseElectron = new THnSparseD ("Electron","Electron",6,binsv1,xminv1,xmaxv1); | |
1811 | // fOutputList->Add(fSparseElectron); | |
33369c14 | 1812 | |
57189f04 | 1813 | PostData(1,fOutputList); |
1814 | } | |
1815 | ||
1816 | //________________________________________________________________________ | |
1817 | void AliAnalysisTaskElecHadronCorrel::Terminate(Option_t *) | |
1818 | { | |
1819 | // Info("Terminate"); | |
520a29a8 | 1820 | AliAnalysisTaskSE::Terminate(); |
57189f04 | 1821 | } |
1822 | ||
1823 | //________________________________________________________________________ | |
1824 | Bool_t AliAnalysisTaskElecHadronCorrel::ProcessCutStep(Int_t cutStep, AliVParticle *track) | |
1825 | { | |
1826 | // Check single track cuts for a given cut step | |
1827 | const Int_t kMCOffset = AliHFEcuts::kNcutStepsMCTrack; | |
1828 | if(!fCFM->CheckParticleCuts(cutStep + kMCOffset, track)) return kFALSE; | |
1829 | return kTRUE; | |
1830 | } | |
987053ce | 1831 | |
20791315 | 1832 | //_________________________________________ |
1833 | void AliAnalysisTaskElecHadronCorrel::SelectPhotonicElectronDCAMet(Int_t itrack, AliVTrack *track, AliPIDResponse *pidResponse) | |
1834 | { | |
1835 | //Identify non-heavy flavour electrons using Invariant mass method using DCA method | |
1836 | fNonHFE = new AliSelectNonHFE(); | |
1837 | fNonHFE->SetAODanalysis(kTRUE); | |
1838 | fNonHFE->SetInvariantMassCut(fInvmassCut); | |
1839 | fNonHFE->SetAlgorithm("DCA"); //KF | |
1840 | fNonHFE->SetPIDresponse(pidResponse); | |
1841 | fNonHFE->SetTrackCuts(-3,3); | |
1842 | ||
1843 | fNonHFE->SetHistMassBack(fDCAMetInvmassLS1); | |
1844 | fNonHFE->SetHistMass(fDCAMetInvmassULS1); | |
1845 | ||
1846 | fNonHFE->FindNonHFE(itrack,track,fVevent); | |
1847 | ||
1848 | Int_t *fULSPartner = fNonHFE->GetPartnersULS(); | |
1849 | Int_t *fLSPartner = fNonHFE->GetPartnersLS(); | |
1850 | Int_t fNULS = fNonHFE->GetNULS(); | |
1851 | Int_t fNLS = fNonHFE->GetNLS(); | |
1852 | ||
1853 | if(fNonHFE->IsULS()){ | |
1854 | for(Int_t k =0; k< fNonHFE->GetNULS(); k++){ | |
1855 | fDCAMetULSElecPt->Fill(track->Pt()); | |
1856 | ||
1857 | ElectronHadCorrel(itrack,track,fDCAMetDphiULSMassLow, fDCAMetDphiULSMassLow1,fDCAMetDphiULSMassLow2,fDCAMetDphiULSMassLow3,fDCAMetDphiULSMassLow4); | |
1858 | ||
1859 | ElectronHadCorrelNoPartnerDCAMet(itrack,track, fULSPartner, fNULS,fDCAMetDphiULSMassLowNoPartner, fDCAMetDphiULSMassLowNoPartner1,fDCAMetDphiULSMassLowNoPartner2,fDCAMetDphiULSMassLowNoPartner3,fDCAMetDphiULSMassLowNoPartner4); | |
1860 | ||
1861 | MixedEvent(track,fDCAMetMixedDphiULSMassLow,fDCAMetMixedDphiULSMassLow1,fDCAMetMixedDphiULSMassLow2, fDCAMetMixedDphiULSMassLow3, fDCAMetMixedDphiULSMassLow4); | |
1862 | } | |
1863 | } | |
1864 | ||
1865 | if(fNonHFE->IsLS()){ | |
1866 | for(Int_t k=0; k < fNonHFE->GetNLS(); k++){ | |
1867 | fDCAMetLSElecPt->Fill(track->Pt()); | |
1868 | ||
1869 | ElectronHadCorrel(itrack,track,fDCAMetDphiLSMassLow,fDCAMetDphiLSMassLow1,fDCAMetDphiLSMassLow2,fDCAMetDphiLSMassLow3,fDCAMetDphiLSMassLow4); | |
1870 | ||
1871 | ElectronHadCorrelNoPartnerDCAMet(itrack,track,fLSPartner, fNLS,fDCAMetDphiLSMassLowNoPartner, fDCAMetDphiLSMassLowNoPartner1,fDCAMetDphiLSMassLowNoPartner2,fDCAMetDphiLSMassLowNoPartner3,fDCAMetDphiLSMassLowNoPartner4); | |
1872 | ||
1873 | MixedEvent(track,fDCAMetMixedDphiLSMassLow,fDCAMetMixedDphiLSMassLow1,fDCAMetMixedDphiLSMassLow2, fDCAMetMixedDphiLSMassLow3, fDCAMetMixedDphiLSMassLow4); | |
1874 | } | |
1875 | } | |
1876 | } | |
57189f04 | 1877 | //_________________________________________ |
987053ce | 1878 | void AliAnalysisTaskElecHadronCorrel::SelectPhotonicElectron(Int_t itrack, AliVTrack *track, Bool_t &fFlagPhotonicElec) |
57189f04 | 1879 | { |
1880 | //Identify non-heavy flavour electrons using Invariant mass method | |
520a29a8 | 1881 | |
57189f04 | 1882 | fTrackCuts1->SetAcceptKinkDaughters(kFALSE); |
1883 | fTrackCuts1->SetRequireTPCRefit(kTRUE); | |
28ea601a | 1884 | fTrackCuts1->SetRequireITSRefit(fAssoEleITSref); |
57189f04 | 1885 | fTrackCuts1->SetEtaRange(-0.9,0.9); |
1886 | fTrackCuts1->SetRequireSigmaToVertex(kTRUE); | |
e4b0faf2 | 1887 | fTrackCuts1->SetMaxChi2PerClusterTPC(4); |
ca49cfe0 | 1888 | fTrackCuts1->SetMinNClustersTPC(fAssoElecTPCNCls); |
987053ce | 1889 | fTrackCuts1->SetMaxDCAToVertexZ(3.2); |
1890 | fTrackCuts1->SetMaxDCAToVertexXY(2.4); | |
1891 | fTrackCuts1->SetDCAToVertex2D(kTRUE); | |
520a29a8 | 1892 | |
57189f04 | 1893 | Bool_t flagPhotonicElec = kFALSE; |
520a29a8 | 1894 | |
987053ce | 1895 | for(Int_t jTracks = 0; jTracks<fVevent->GetNumberOfTracks(); jTracks++){ |
1896 | AliVParticle* VtrackAsso = fVevent->GetTrack(jTracks); | |
1897 | if (!VtrackAsso) { | |
57189f04 | 1898 | printf("ERROR: Could not receive track %d\n", jTracks); |
1899 | continue; | |
1900 | } | |
987053ce | 1901 | |
1902 | AliVTrack *trackAsso = dynamic_cast<AliVTrack*>(VtrackAsso); | |
7e02de85 | 1903 | if(!trackAsso) continue; |
987053ce | 1904 | |
1905 | //track cuts applied | |
1906 | if(IsAODanalysis()) { | |
1907 | AliAODTrack *atrackAsso = dynamic_cast<AliAODTrack*>(VtrackAsso); | |
7e02de85 | 1908 | if(!atrackAsso) continue; |
987053ce | 1909 | if(!atrackAsso->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
ca49cfe0 | 1910 | if(atrackAsso->GetTPCNcls() < fAssoElecTPCNCls) continue; |
28ea601a | 1911 | if(!(atrackAsso->GetStatus()&AliESDtrack::kTPCrefit)) continue; |
1912 | ||
1913 | if(fAssoEleITSref){ | |
1914 | if(!(atrackAsso->GetStatus()&AliESDtrack::kITSrefit)) continue; | |
1915 | } | |
987053ce | 1916 | } |
1917 | else{ | |
1918 | AliESDtrack *etrackAsso = dynamic_cast<AliESDtrack*>(VtrackAsso); | |
7e02de85 | 1919 | if(!etrackAsso) continue; |
987053ce | 1920 | if(!fTrackCuts1->AcceptTrack(etrackAsso)) continue; |
1921 | } | |
1922 | ||
e4b0faf2 | 1923 | if(jTracks==itrack) continue; |
520a29a8 | 1924 | |
987053ce | 1925 | Double_t dEdxAsso = -999., ptAsso=-999., nsigma=-999.0; |
520a29a8 | 1926 | Double_t mass=-999., width = -999; |
57189f04 | 1927 | Bool_t fFlagLS=kFALSE, fFlagULS=kFALSE; |
520a29a8 | 1928 | |
57189f04 | 1929 | dEdxAsso = trackAsso->GetTPCsignal(); |
487ae960 | 1930 | nsigma = fPID->GetPIDResponse() ? fPID->GetPIDResponse()->NumberOfSigmasTPC(trackAsso, AliPID::kElectron) : 1000; |
57189f04 | 1931 | ptAsso = trackAsso->Pt(); |
1932 | Int_t chargeAsso = trackAsso->Charge(); | |
1933 | Int_t charge = track->Charge(); | |
520a29a8 | 1934 | |
57189f04 | 1935 | if(ptAsso <0.3) continue; |
3c56855b | 1936 | if(trackAsso->Eta()<-0.9 || trackAsso->Eta()>0.9) continue; |
487ae960 | 1937 | if(nsigma < -3 || nsigma > 3) continue; |
520a29a8 | 1938 | |
57189f04 | 1939 | Int_t fPDGe1 = 11; Int_t fPDGe2 = 11; |
1940 | if(charge>0) fPDGe1 = -11; | |
1941 | if(chargeAsso>0) fPDGe2 = -11; | |
520a29a8 | 1942 | |
57189f04 | 1943 | if(charge == chargeAsso) fFlagLS = kTRUE; |
1944 | if(charge != chargeAsso) fFlagULS = kTRUE; | |
520a29a8 | 1945 | |
20791315 | 1946 | AliKFParticle::SetField(fVevent->GetMagneticField()); |
1947 | ||
520a29a8 | 1948 | AliKFParticle ge1 = AliKFParticle(*track, fPDGe1); |
1949 | AliKFParticle ge2 = AliKFParticle(*trackAsso, fPDGe2); | |
57189f04 | 1950 | AliKFParticle recg(ge1, ge2); |
3c56855b | 1951 | |
57189f04 | 1952 | if(recg.GetNDF()<1) continue; |
1953 | Double_t chi2recg = recg.GetChi2()/recg.GetNDF(); | |
1954 | if(TMath::Sqrt(TMath::Abs(chi2recg))>3.) continue; | |
520a29a8 | 1955 | |
3c56855b | 1956 | Int_t MassCorrect; |
1957 | MassCorrect = recg.GetMass(mass,width); | |
520a29a8 | 1958 | |
487ae960 | 1959 | if(fFlagLS) { |
33369c14 | 1960 | if(track->Pt()>2)fInvmassLS1->Fill(mass); |
a9123929 | 1961 | // if(track->Pt()> 1) fInvmassLS2->Fill(mass); |
1962 | // if(track->Pt()>2) fInvmassLS3->Fill(mass); | |
1963 | // if(track->Pt()>3) fInvmassLS4->Fill(mass); | |
1964 | // if(track->Pt()>4) fInvmassLS5->Fill(mass); | |
487ae960 | 1965 | } |
1966 | if(fFlagULS) { | |
33369c14 | 1967 | if(track->Pt()>2)fInvmassULS1->Fill(mass); |
1968 | //if(track->Pt() >1) fInvmassULS2->Fill(mass); | |
1969 | //if(track->Pt() >2) fInvmassULS3->Fill(mass); | |
1970 | //if(track->Pt() >3) fInvmassULS4->Fill(mass); | |
1971 | //if(track->Pt() >4) fInvmassULS5->Fill(mass); | |
487ae960 | 1972 | } |
520a29a8 | 1973 | |
520a29a8 | 1974 | if(mass<fInvmassCut){ |
1975 | if(fFlagULS) | |
1976 | { | |
487ae960 | 1977 | ElectronHadCorrel(itrack,track,fDphiULSMassLow, fDphiULSMassLow1,fDphiULSMassLow2,fDphiULSMassLow3,fDphiULSMassLow4); |
bf697dda | 1978 | ElectronHadCorrelEtaBins(itrack,track,fDphiULSMassLowEta1, fDphiULSMassLowEta11,fDphiULSMassLowEta12,fDphiULSMassLowEta13,fDphiULSMassLowEta14,fDphiULSMassLowEta2, fDphiULSMassLowEta21,fDphiULSMassLowEta22,fDphiULSMassLowEta23,fDphiULSMassLowEta24); |
520a29a8 | 1979 | fULSElecPt->Fill(track->Pt()); |
3db57d71 | 1980 | MixedEvent(track,fMixedDphiULSMassLow,fMixedDphiULSMassLow1,fMixedDphiULSMassLow2, fMixedDphiULSMassLow3, fMixedDphiULSMassLow4); |
520a29a8 | 1981 | } |
1982 | if(fFlagLS) | |
1983 | { | |
487ae960 | 1984 | ElectronHadCorrel(itrack,track,fDphiLSMassLow,fDphiLSMassLow1,fDphiLSMassLow2,fDphiLSMassLow3,fDphiLSMassLow4); |
bf697dda | 1985 | ElectronHadCorrelEtaBins(itrack,track,fDphiLSMassLowEta1, fDphiLSMassLowEta11,fDphiLSMassLowEta12,fDphiLSMassLowEta13,fDphiLSMassLowEta14,fDphiLSMassLowEta2, fDphiLSMassLowEta21,fDphiLSMassLowEta22,fDphiLSMassLowEta23,fDphiLSMassLowEta24); |
520a29a8 | 1986 | fLSElecPt->Fill(track->Pt()); |
3db57d71 | 1987 | MixedEvent(track,fMixedDphiLSMassLow,fMixedDphiLSMassLow1,fMixedDphiLSMassLow2, fMixedDphiLSMassLow3, fMixedDphiLSMassLow4); |
520a29a8 | 1988 | } |
bf697dda | 1989 | if(fFlagLS){ |
1990 | ElectronHadCorrelNoPartner(itrack,jTracks,track,fDphiLSMassLowNoPartner, fDphiLSMassLowNoPartner1,fDphiLSMassLowNoPartner2,fDphiLSMassLowNoPartner3,fDphiLSMassLowNoPartner4); | |
1991 | ElectronHadCorrelEtaBinsNoPartner(itrack,jTracks,track,fDphiLSMassLowNoPartnerEta1, fDphiLSMassLowNoPartnerEta11,fDphiLSMassLowNoPartnerEta12,fDphiLSMassLowNoPartnerEta13,fDphiLSMassLowNoPartnerEta14,fDphiLSMassLowNoPartnerEta2, fDphiLSMassLowNoPartnerEta21,fDphiLSMassLowNoPartnerEta22,fDphiLSMassLowNoPartnerEta23,fDphiLSMassLowNoPartnerEta24); | |
1992 | } | |
1993 | if(fFlagULS) { | |
1994 | ElectronHadCorrelNoPartner(itrack,jTracks,track,fDphiULSMassLowNoPartner, fDphiULSMassLowNoPartner1,fDphiULSMassLowNoPartner2,fDphiULSMassLowNoPartner3,fDphiULSMassLowNoPartner4); | |
1995 | ElectronHadCorrelEtaBinsNoPartner(itrack,jTracks,track,fDphiULSMassLowNoPartnerEta1, fDphiULSMassLowNoPartnerEta11,fDphiULSMassLowNoPartnerEta12,fDphiULSMassLowNoPartnerEta13,fDphiULSMassLowNoPartnerEta14,fDphiULSMassLowNoPartnerEta2, fDphiULSMassLowNoPartnerEta21,fDphiULSMassLowNoPartnerEta22,fDphiULSMassLowNoPartnerEta23,fDphiULSMassLowNoPartnerEta24); | |
1996 | } | |
57189f04 | 1997 | } |
520a29a8 | 1998 | |
57189f04 | 1999 | if(mass<fInvmassCut && fFlagULS && !flagPhotonicElec){ |
2000 | flagPhotonicElec = kTRUE; | |
2001 | } | |
520a29a8 | 2002 | // } |
2003 | ||
987053ce | 2004 | } |
2005 | fFlagPhotonicElec = flagPhotonicElec; | |
57189f04 | 2006 | } |
2007 | //_________________________________________ | |
987053ce | 2008 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrel(Int_t itrack, AliVTrack *track, TH2F *DphiPt, TH2F *DphiPt1,TH2F *DphiPt2,TH2F *DphiPt3,TH2F *DphiPt4) |
57189f04 | 2009 | { |
2010 | //Construct Delta Phi between electrons and hadrons | |
520a29a8 | 2011 | |
57189f04 | 2012 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); |
2013 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2014 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2015 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2016 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
e4b0faf2 | 2017 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); |
28ea601a | 2018 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
987053ce | 2019 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2020 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2021 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
520a29a8 | 2022 | |
987053ce | 2023 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ |
2024 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2025 | if (!VtrackHad) { | |
57189f04 | 2026 | printf("ERROR: Could not receive track %d\n", ktracks); |
2027 | continue; | |
2028 | } | |
987053ce | 2029 | |
2030 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
20791315 | 2031 | if(!trackHad) continue; |
987053ce | 2032 | |
2033 | if(IsAODanalysis()) { | |
2034 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2035 | if(!atrackHad) continue; |
987053ce | 2036 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2037 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2038 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
987053ce | 2039 | } |
2040 | else{ | |
2041 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2042 | if(!etrackHad) continue; |
987053ce | 2043 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2044 | } | |
a9123929 | 2045 | |
57189f04 | 2046 | if(ktracks == itrack) continue; //do not select the same electron |
520a29a8 | 2047 | |
57189f04 | 2048 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; |
2049 | Double_t ptEle = -999; | |
2050 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2051 | Double_t pi = 3.14; | |
520a29a8 | 2052 | |
57189f04 | 2053 | dEdxHad = trackHad->GetTPCsignal(); |
2054 | ptHad = trackHad->Pt(); | |
2055 | pHad = trackHad->P(); | |
4e01b68c | 2056 | ptEle = track->Pt(); |
520a29a8 | 2057 | |
987053ce | 2058 | // if(ptHad <2) continue; |
4e01b68c | 2059 | if(ptHad > ptEle) continue; |
3c56855b | 2060 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; |
3c56855b | 2061 | |
57189f04 | 2062 | phiEle = track->Phi(); |
2063 | phiHad = trackHad->Phi(); | |
2064 | Dphi = phiEle - phiHad; | |
520a29a8 | 2065 | if (Dphi > 3*pi/2) |
57189f04 | 2066 | Dphi = Dphi - 2*pi; |
520a29a8 | 2067 | if (Dphi < -pi/2) |
57189f04 | 2068 | Dphi = Dphi + 2*pi; |
520a29a8 | 2069 | |
487ae960 | 2070 | if(ptHad>2) DphiPt->Fill(ptEle,Dphi); |
2071 | if(ptHad>2 && ptHad<4) DphiPt1->Fill(ptEle,Dphi); | |
2072 | if(ptHad>4 && ptHad<6) DphiPt2->Fill(ptEle,Dphi); | |
2073 | if(ptHad>6 && ptHad<8) DphiPt3->Fill(ptEle,Dphi); | |
a9123929 | 2074 | if(ptHad>4 && ptHad<10) DphiPt4->Fill(ptEle,Dphi); |
4e01b68c | 2075 | |
57189f04 | 2076 | } |
2077 | } | |
2078 | //_________________________________________ | |
987053ce | 2079 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrelNoPartner(Int_t itrack,Int_t jtrack, AliVTrack *track, TH2F *DphiPtNew, TH2F *DphiPtNew1,TH2F *DphiPtNew2,TH2F *DphiPtNew3,TH2F *DphiPtNew4) |
57189f04 | 2080 | { |
520a29a8 | 2081 | //Construct Delta Phi between electrons and hadrons for electrons from invariant mass calculation excluding associated track |
57189f04 | 2082 | |
520a29a8 | 2083 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); |
2084 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2085 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2086 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2087 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
e4b0faf2 | 2088 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); |
28ea601a | 2089 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
987053ce | 2090 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2091 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2092 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
520a29a8 | 2093 | |
987053ce | 2094 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ |
2095 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2096 | if (!VtrackHad) { | |
520a29a8 | 2097 | printf("ERROR: Could not receive track %d\n", ktracks); |
2098 | continue; | |
2099 | } | |
a9123929 | 2100 | |
987053ce | 2101 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); |
7e02de85 | 2102 | if(!trackHad) continue; |
520a29a8 | 2103 | |
987053ce | 2104 | if(IsAODanalysis()) { |
2105 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2106 | if(!atrackHad) continue; |
987053ce | 2107 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2108 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2109 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
987053ce | 2110 | } |
2111 | else{ | |
2112 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2113 | if(!etrackHad) continue; |
987053ce | 2114 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2115 | } | |
2116 | ||
2117 | if(ktracks == itrack || ktracks == jtrack) continue; //do not select the same electron and associated track from inv mass cal | |
520a29a8 | 2118 | |
2119 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; | |
2120 | Double_t ptEle = -999; | |
2121 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2122 | Double_t pi = 3.14; | |
2123 | ||
2124 | dEdxHad = trackHad->GetTPCsignal(); | |
2125 | ptHad = trackHad->Pt(); | |
2126 | pHad = trackHad->P(); | |
4e01b68c | 2127 | ptEle = track->Pt(); |
520a29a8 | 2128 | |
a9123929 | 2129 | // if(ptHad <2) continue; |
4e01b68c | 2130 | if(ptHad > ptEle) continue; |
3c56855b | 2131 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; |
520a29a8 | 2132 | |
2133 | phiEle = track->Phi(); | |
2134 | phiHad = trackHad->Phi(); | |
2135 | Dphi = phiEle - phiHad; | |
2136 | if (Dphi > 3*pi/2) | |
2137 | Dphi = Dphi - 2*pi; | |
2138 | if (Dphi < -pi/2) | |
2139 | Dphi = Dphi + 2*pi; | |
57189f04 | 2140 | |
487ae960 | 2141 | if(ptHad>2) DphiPtNew->Fill(ptEle,Dphi); |
2142 | if(ptHad>2 && ptHad<4) DphiPtNew1->Fill(ptEle,Dphi); | |
2143 | if(ptHad>4 && ptHad<6) DphiPtNew2->Fill(ptEle,Dphi); | |
2144 | if(ptHad>6 && ptHad<8) DphiPtNew3->Fill(ptEle,Dphi); | |
a9123929 | 2145 | if(ptHad>4 && ptHad<10) DphiPtNew4->Fill(ptEle,Dphi); |
520a29a8 | 2146 | } |
2147 | } | |
20791315 | 2148 | //_________________________________________ |
2149 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrelNoPartnerDCAMet(Int_t itrack,AliVTrack *track,Int_t * AssoList, Int_t NAsso, TH2F *DphiPtNew, TH2F *DphiPtNew1,TH2F *DphiPtNew2,TH2F *DphiPtNew3,TH2F *DphiPtNew4) | |
2150 | { | |
2151 | //Construct Delta Phi between electrons and hadrons for electrons from invariant mass calculation excluding associated track | |
2152 | ||
2153 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); | |
2154 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2155 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2156 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2157 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
2158 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); | |
2159 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); | |
2160 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); | |
2161 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2162 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
2163 | ||
2164 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ | |
2165 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2166 | if (!VtrackHad) { | |
2167 | printf("ERROR: Could not receive track %d\n", ktracks); | |
2168 | continue; | |
2169 | } | |
2170 | ||
2171 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
2172 | if(!trackHad) continue; | |
2173 | ||
2174 | if(IsAODanalysis()) { | |
2175 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
2176 | if(!atrackHad) continue; | |
2177 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; | |
2178 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
2179 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; | |
2180 | } | |
2181 | else{ | |
2182 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
2183 | if(!etrackHad) continue; | |
2184 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; | |
2185 | } | |
2186 | ||
2187 | if(ktracks == itrack ) continue; //do not select the same electron and associated track from inv mass cal | |
2188 | Bool_t IsPartner = kFALSE; | |
2189 | //Check if this track is a Non-HFE partner | |
2190 | for(Int_t i = 0; i < NAsso; i++) | |
2191 | { | |
2192 | if(AssoList[i]==ktracks) IsPartner=kTRUE; | |
2193 | } | |
2194 | if(IsPartner) continue; | |
2195 | ||
2196 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; | |
2197 | Double_t ptEle = -999; | |
2198 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2199 | Double_t pi = 3.14; | |
2200 | ||
2201 | dEdxHad = trackHad->GetTPCsignal(); | |
2202 | ptHad = trackHad->Pt(); | |
2203 | pHad = trackHad->P(); | |
2204 | ptEle = track->Pt(); | |
2205 | ||
2206 | // if(ptHad <2) continue; | |
2207 | if(ptHad > ptEle) continue; | |
2208 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; | |
2209 | ||
2210 | phiEle = track->Phi(); | |
2211 | phiHad = trackHad->Phi(); | |
2212 | Dphi = phiEle - phiHad; | |
2213 | if (Dphi > 3*pi/2) | |
2214 | Dphi = Dphi - 2*pi; | |
2215 | if (Dphi < -pi/2) | |
2216 | Dphi = Dphi + 2*pi; | |
2217 | ||
2218 | if(ptHad>2) DphiPtNew->Fill(ptEle,Dphi); | |
2219 | if(ptHad>2 && ptHad<4) DphiPtNew1->Fill(ptEle,Dphi); | |
2220 | if(ptHad>4 && ptHad<6) DphiPtNew2->Fill(ptEle,Dphi); | |
2221 | if(ptHad>6 && ptHad<8) DphiPtNew3->Fill(ptEle,Dphi); | |
2222 | if(ptHad>4 && ptHad<10) DphiPtNew4->Fill(ptEle,Dphi); | |
2223 | } | |
2224 | } | |
2225 | ||
2226 | ||
4e01b68c | 2227 | //_________________________________________ |
33369c14 | 2228 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrelEtaFarSide(Int_t itrack, AliVTrack *track, TH2F *DphiPt, TH2F *DphiPt1,TH2F *DphiPt2,TH2F *DphiPt3,TH2F *DphiPt4) |
4e01b68c | 2229 | { |
33369c14 | 2230 | //Construct Delta Phi between electrons and hadrons for 1<eta(had)<1.6 |
4e01b68c | 2231 | |
33369c14 | 2232 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); |
2233 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2234 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2235 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2236 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
2237 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); | |
28ea601a | 2238 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
33369c14 | 2239 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2240 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2241 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
4e01b68c | 2242 | |
33369c14 | 2243 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ |
2244 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2245 | if (!VtrackHad) { | |
2246 | printf("ERROR: Could not receive track %d\n", ktracks); | |
2247 | continue; | |
2248 | } | |
2249 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
20791315 | 2250 | if(!trackHad) continue; |
4e01b68c | 2251 | |
33369c14 | 2252 | if(IsAODanalysis()) { |
2253 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2254 | if(!atrackHad) continue; |
33369c14 | 2255 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2256 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2257 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
33369c14 | 2258 | } |
2259 | else{ | |
2260 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2261 | if(!etrackHad) continue; |
33369c14 | 2262 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2263 | } | |
4e01b68c | 2264 | |
33369c14 | 2265 | if(ktracks == itrack) continue; //do not select the same electron |
4e01b68c | 2266 | |
33369c14 | 2267 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; |
2268 | Double_t ptEle = -999; | |
2269 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2270 | Double_t pi = 3.14; | |
4e01b68c | 2271 | |
33369c14 | 2272 | dEdxHad = trackHad->GetTPCsignal(); |
2273 | ptHad = trackHad->Pt(); | |
2274 | pHad = trackHad->P(); | |
2275 | ptEle = track->Pt(); | |
4e01b68c | 2276 | |
33369c14 | 2277 | // if(ptHad <2) continue; |
2278 | if(ptHad > ptEle) continue; | |
2279 | if(trackHad->Eta()<1 || trackHad->Eta()>1.6) continue; | |
2280 | ||
2281 | phiEle = track->Phi(); | |
2282 | phiHad = trackHad->Phi(); | |
2283 | Dphi = phiEle - phiHad; | |
2284 | if (Dphi > 3*pi/2) | |
2285 | Dphi = Dphi - 2*pi; | |
2286 | if (Dphi < -pi/2) | |
2287 | Dphi = Dphi + 2*pi; | |
2288 | ||
2289 | if(ptHad>2) DphiPt->Fill(ptEle,Dphi); | |
2290 | if(ptHad>2 && ptHad<4) DphiPt1->Fill(ptEle,Dphi); | |
2291 | if(ptHad>4 && ptHad<6) DphiPt2->Fill(ptEle,Dphi); | |
2292 | if(ptHad>6 && ptHad<8) DphiPt3->Fill(ptEle,Dphi); | |
a9123929 | 2293 | if(ptHad>4 && ptHad<10) DphiPt4->Fill(ptEle,Dphi); |
33369c14 | 2294 | |
2295 | } | |
4e01b68c | 2296 | } |
bf697dda | 2297 | //_________________________________________ |
2298 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrelEtaBins(Int_t itrack, AliVTrack *track, TH2F *DphiPtEta1, TH2F *DphiPtEta11,TH2F *DphiPtEta12,TH2F *DphiPtEta13,TH2F *DphiPtEta14,TH2F *DphiPtEta2, TH2F *DphiPtEta21,TH2F *DphiPtEta22,TH2F *DphiPtEta23,TH2F *DphiPtEta24) | |
2299 | { | |
2300 | //Construct Delta Phi between electrons and hadrons in different Eta bins | |
2301 | ||
2302 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); | |
2303 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2304 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2305 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2306 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
2307 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); | |
28ea601a | 2308 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
bf697dda | 2309 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2310 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2311 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
2312 | ||
2313 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ | |
2314 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2315 | if (!VtrackHad) { | |
2316 | printf("ERROR: Could not receive track %d\n", ktracks); | |
2317 | continue; | |
2318 | } | |
2319 | ||
2320 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
20791315 | 2321 | if(!trackHad) continue; |
bf697dda | 2322 | |
2323 | if(IsAODanalysis()) { | |
2324 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2325 | if(!atrackHad) continue; |
bf697dda | 2326 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2327 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2328 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
bf697dda | 2329 | } |
2330 | else{ | |
2331 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2332 | if(!etrackHad) continue; |
bf697dda | 2333 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2334 | } | |
2335 | ||
2336 | if(ktracks == itrack) continue; //do not select the same electron | |
2337 | ||
2338 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; | |
2339 | Double_t ptEle = -999; | |
2340 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2341 | Double_t etaEle = -999, etaHad=-999, DEta=-999; | |
2342 | Double_t pi = 3.14; | |
2343 | ||
2344 | dEdxHad = trackHad->GetTPCsignal(); | |
2345 | ptHad = trackHad->Pt(); | |
2346 | pHad = trackHad->P(); | |
2347 | ptEle = track->Pt(); | |
2348 | ||
2349 | // if(ptHad <2) continue; | |
2350 | if(ptHad > ptEle) continue; | |
2351 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; | |
2352 | ||
2353 | phiEle = track->Phi(); | |
2354 | phiHad = trackHad->Phi(); | |
2355 | Dphi = phiEle - phiHad; | |
2356 | ||
2357 | if (Dphi > 3*pi/2) | |
2358 | Dphi = Dphi - 2*pi; | |
2359 | if (Dphi < -pi/2) | |
2360 | Dphi = Dphi + 2*pi; | |
2361 | ||
2362 | etaEle = track->Eta(); | |
2363 | etaHad = trackHad->Eta(); | |
2364 | DEta = etaEle - etaHad; | |
2365 | ||
2366 | if(TMath::Abs(DEta)<0.8){ | |
a9123929 | 2367 | if(ptHad>2) DphiPtEta1->Fill(ptEle,Dphi); |
2368 | if(ptHad>2 && ptHad<4) DphiPtEta11->Fill(ptEle,Dphi); | |
2369 | if(ptHad>4 && ptHad<6) DphiPtEta12->Fill(ptEle,Dphi); | |
2370 | if(ptHad>6 && ptHad<8) DphiPtEta13->Fill(ptEle,Dphi); | |
2371 | if(ptHad>4 && ptHad<10) DphiPtEta14->Fill(ptEle,Dphi); | |
bf697dda | 2372 | } |
2373 | ||
2374 | if(TMath::Abs(DEta)>0.8){ | |
2375 | if(ptHad>2) DphiPtEta2->Fill(ptEle,Dphi); | |
2376 | if(ptHad>2 && ptHad<4) DphiPtEta21->Fill(ptEle,Dphi); | |
2377 | if(ptHad>4 && ptHad<6) DphiPtEta22->Fill(ptEle,Dphi); | |
2378 | if(ptHad>6 && ptHad<8) DphiPtEta23->Fill(ptEle,Dphi); | |
a9123929 | 2379 | if(ptHad>4 && ptHad<10) DphiPtEta24->Fill(ptEle,Dphi); |
bf697dda | 2380 | } |
2381 | } | |
2382 | } | |
2383 | //_________________________________________ | |
2384 | void AliAnalysisTaskElecHadronCorrel::ElectronHadCorrelEtaBinsNoPartner(Int_t itrack,Int_t jtrack, AliVTrack *track, TH2F *DphiPtEta1, TH2F *DphiPtEta11,TH2F *DphiPtEta12,TH2F *DphiPtEta13,TH2F *DphiPtEta14,TH2F *DphiPtEta2, TH2F *DphiPtEta21,TH2F *DphiPtEta22,TH2F *DphiPtEta23,TH2F *DphiPtEta24) | |
2385 | { | |
2386 | //Construct Delta Phi between electrons and hadrons in different Eta bins (exclude associated track) | |
2387 | ||
2388 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); | |
2389 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2390 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2391 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2392 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
2393 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); | |
28ea601a | 2394 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
bf697dda | 2395 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2396 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2397 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
2398 | ||
2399 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ | |
2400 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2401 | if (!VtrackHad) { | |
2402 | printf("ERROR: Could not receive track %d\n", ktracks); | |
2403 | continue; | |
2404 | } | |
2405 | ||
2406 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
20791315 | 2407 | if(!trackHad) continue; |
bf697dda | 2408 | |
2409 | if(IsAODanalysis()) { | |
2410 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2411 | if(!atrackHad) continue; |
bf697dda | 2412 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2413 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2414 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
bf697dda | 2415 | } |
2416 | else{ | |
2417 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2418 | if(!etrackHad) continue; |
bf697dda | 2419 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2420 | } | |
2421 | ||
2422 | if(ktracks == itrack || ktracks == jtrack) continue; //do not select the same electron and associated track from inv mass cal | |
2423 | ||
2424 | Double_t ptHad= -999, pHad=-999., dEdxHad = -999; | |
2425 | Double_t ptEle = -999; | |
2426 | Double_t phiEle = -999, phiHad = -999, Dphi = -999; | |
2427 | Double_t etaEle = -999, etaHad=-999, DEta=-999; | |
2428 | Double_t pi = 3.14; | |
2429 | ||
2430 | dEdxHad = trackHad->GetTPCsignal(); | |
2431 | ptHad = trackHad->Pt(); | |
2432 | pHad = trackHad->P(); | |
2433 | ptEle = track->Pt(); | |
2434 | ||
2435 | // if(ptHad <2) continue; | |
2436 | if(ptHad > ptEle) continue; | |
2437 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; | |
2438 | ||
2439 | phiEle = track->Phi(); | |
2440 | phiHad = trackHad->Phi(); | |
2441 | Dphi = phiEle - phiHad; | |
2442 | ||
2443 | if (Dphi > 3*pi/2) | |
2444 | Dphi = Dphi - 2*pi; | |
2445 | if (Dphi < -pi/2) | |
2446 | Dphi = Dphi + 2*pi; | |
2447 | ||
2448 | etaEle = track->Eta(); | |
2449 | etaHad = trackHad->Eta(); | |
2450 | DEta = etaEle - etaHad; | |
2451 | ||
2452 | if(TMath::Abs(DEta)<0.8){ | |
2453 | if(ptHad>2) DphiPtEta1->Fill(ptEle,Dphi); | |
2454 | if(ptHad>2 && ptHad<4) DphiPtEta11->Fill(ptEle,Dphi); | |
2455 | if(ptHad>4 && ptHad<6) DphiPtEta12->Fill(ptEle,Dphi); | |
2456 | if(ptHad>6 && ptHad<8) DphiPtEta13->Fill(ptEle,Dphi); | |
a9123929 | 2457 | if(ptHad>4 && ptHad<10) DphiPtEta14->Fill(ptEle,Dphi); |
bf697dda | 2458 | } |
2459 | ||
2460 | if(TMath::Abs(DEta)>0.8){ | |
2461 | if(ptHad>2) DphiPtEta2->Fill(ptEle,Dphi); | |
2462 | if(ptHad>2 && ptHad<4) DphiPtEta21->Fill(ptEle,Dphi); | |
2463 | if(ptHad>4 && ptHad<6) DphiPtEta22->Fill(ptEle,Dphi); | |
2464 | if(ptHad>6 && ptHad<8) DphiPtEta23->Fill(ptEle,Dphi); | |
a9123929 | 2465 | if(ptHad>4 && ptHad<10) DphiPtEta24->Fill(ptEle,Dphi); |
bf697dda | 2466 | } |
2467 | } | |
2468 | } | |
3db57d71 | 2469 | |
33369c14 | 2470 | //_________________________________________ |
a9123929 | 2471 | void AliAnalysisTaskElecHadronCorrel::MixedEvent(AliVTrack *track, TH2F *DphiPt, TH2F *DphiPt1,TH2F *DphiPt2, TH2F *DphiPt3, TH2F *DphiPt4) |
4e01b68c | 2472 | { |
4e01b68c | 2473 | |
a9123929 | 2474 | const AliAODVertex *pVtx = fAOD->GetPrimaryVertex(); |
2475 | Double_t zVtx; | |
2476 | zVtx = pVtx->GetZ(); | |
4e01b68c | 2477 | |
7e02de85 | 2478 | AliCentrality *fCentrality1 = (AliCentrality*)fVevent->GetCentrality(); |
b8c2fb6d | 2479 | Double_t centvalue = fCentrality1->GetCentralityPercentile(fkCentralityMethod); |
3c56855b | 2480 | |
7e02de85 | 2481 | AliEventPool* pool; |
2482 | pool = fPoolMgr->GetEventPool(centvalue, zVtx); // Get the buffer associated with the current centrality and z-vtx | |
a9123929 | 2483 | if (!pool) |
7e02de85 | 2484 | { |
a9123929 | 2485 | AliFatal(Form("No pool found for centrality = %f, zVtx = %f", centvalue, zVtx)); |
7e02de85 | 2486 | return; |
2487 | } | |
a9123929 | 2488 | // pool->PrintInfo(); |
2489 | if (pool->GetCurrentNEvents() >= 5) // start mixing when 5 events are in the buffer | |
2490 | { | |
2491 | Int_t nMix = pool->GetCurrentNEvents(); | |
2492 | fNoMixedEvents->Fill(0); | |
2493 | fMixStat->Fill(pool->GetCurrentNEvents(),centvalue); | |
2494 | fMixStat1->Fill(pool->GetCurrentNEvents(),zVtx); | |
2495 | ||
2496 | // cout << "nMix = " << nMix << " tracks in pool = " << pool->NTracksInPool() << endl; | |
2497 | for (Int_t jMix=0; jMix<nMix; jMix++) // mix with each event in the buffer | |
2498 | { | |
2499 | TObjArray* bgTracks = pool->GetEvent(jMix); | |
2500 | for (Int_t i=0;i<bgTracks->GetEntriesFast(); i++) | |
2501 | { | |
2502 | AliVParticle* mixtrk = (AliVParticle*) bgTracks->At(i); | |
2503 | ||
2504 | Double_t mixtrkPhi = -999; | |
2505 | Double_t ptEle = -999; | |
2506 | Double_t phiEle = -999, Dphi = -999; | |
2507 | Double_t pi = 3.14; | |
2508 | Double_t ptmixtrk = -999; | |
2509 | ||
2510 | ptEle = track->Pt(); | |
2511 | ptmixtrk = mixtrk->Pt(); | |
2512 | if(ptmixtrk > ptEle) continue; | |
2513 | ||
2514 | mixtrkPhi = mixtrk->Phi(); | |
2515 | phiEle = track->Phi(); | |
2516 | Dphi = phiEle - mixtrkPhi; | |
2517 | ||
2518 | if (Dphi > 3*pi/2) | |
2519 | Dphi = Dphi - 2*pi; | |
2520 | if (Dphi < -pi/2) | |
2521 | Dphi = Dphi + 2*pi; | |
2522 | if(ptmixtrk>2) DphiPt->Fill(ptEle,Dphi); | |
2523 | if(ptmixtrk>2 && ptmixtrk<4) DphiPt1->Fill(ptEle,Dphi); | |
2524 | if(ptmixtrk>4 && ptmixtrk<6) DphiPt2->Fill(ptEle,Dphi); | |
2525 | if(ptmixtrk>6 && ptmixtrk<8) DphiPt3->Fill(ptEle,Dphi); | |
2526 | if(ptmixtrk>4 && ptmixtrk<10) DphiPt4->Fill(ptEle,Dphi); | |
2527 | } | |
2528 | } | |
2529 | } | |
33369c14 | 2530 | } |
3db57d71 | 2531 | |
33369c14 | 2532 | //___________________________________________ |
2533 | TObjArray* AliAnalysisTaskElecHadronCorrel::CloneAndReduceTrackList() | |
2534 | { | |
a9123929 | 2535 | // clones a track list by using AliehDPhiBasicParticle which uses much less memory (used for event mixing) |
2536 | ||
2537 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); | |
2538 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2539 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2540 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2541 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
2542 | fTrackCuts2->SetMaxChi2PerClusterTPC(3.5); | |
28ea601a | 2543 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
7e02de85 | 2544 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2545 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2546 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
a9123929 | 2547 | |
2548 | TObjArray* tracksClone = new TObjArray; | |
2549 | tracksClone->SetOwner(kTRUE); | |
2550 | ||
7e02de85 | 2551 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ |
2552 | AliVParticle* Vtrack = fVevent->GetTrack(ktracks); | |
2553 | if (!Vtrack) { | |
a9123929 | 2554 | printf("ERROR: Could not receive track %d\n", ktracks); |
2555 | continue; | |
2556 | } | |
7e02de85 | 2557 | |
2558 | AliVTrack *track = dynamic_cast<AliVTrack*>(Vtrack); | |
20791315 | 2559 | if(!track) continue; |
7e02de85 | 2560 | |
2561 | if(IsAODanalysis()) { | |
2562 | AliAODTrack *atrack = dynamic_cast<AliAODTrack*>(Vtrack); | |
2563 | if(!atrack) continue; | |
2564 | if(!atrack->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; | |
2565 | if((!(atrack->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrack->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2566 | if(atrack->GetTPCNcls() < fTPCNClsHad) continue; |
7e02de85 | 2567 | } |
2568 | else{ | |
2569 | AliESDtrack *etrack = dynamic_cast<AliESDtrack*>(Vtrack); | |
2570 | if(!etrack) continue; | |
2571 | if(!fTrackCuts2->AcceptTrack(etrack)) continue; | |
2572 | } | |
a9123929 | 2573 | |
2574 | // if(ktracks == iTrack) continue; | |
2575 | Double_t eta=-999,ptHad= -999, pHad=-999., phi=-999.0; | |
2576 | Int_t label=-9999, id=-999; | |
2577 | eta = track->Eta(); | |
2578 | ptHad = track->Pt(); | |
2579 | pHad = track->P(); | |
2580 | phi= track->Phi(); | |
2581 | label= track->GetLabel(); | |
2582 | id=track->GetID(); | |
2583 | ||
7e02de85 | 2584 | if(track->Eta()<-0.9 || track->Eta()>0.9) continue; |
a9123929 | 2585 | if(ptHad <2) continue; |
a9123929 | 2586 | |
7e02de85 | 2587 | AliVParticle* particle = (AliVParticle*) fVevent->GetTrack(ktracks); |
a9123929 | 2588 | tracksClone->Add(new AliehDPhiBasicParticle(particle->Eta(), particle->Phi(), particle->Pt(), particle->Charge())); |
2589 | ||
4e01b68c | 2590 | } |
a9123929 | 2591 | return tracksClone; |
33369c14 | 2592 | } |
3db57d71 | 2593 | |
3c56855b | 2594 | //___________________________________________ |
2595 | void AliAnalysisTaskElecHadronCorrel::HadronInfo(Int_t itrack) | |
2596 | { | |
2597 | //Hadron information | |
987053ce | 2598 | |
3c56855b | 2599 | fTrackCuts2->SetAcceptKinkDaughters(kFALSE); |
2600 | fTrackCuts2->SetRequireTPCRefit(kTRUE); | |
2601 | fTrackCuts2->SetRequireITSRefit(kTRUE); | |
2602 | fTrackCuts2->SetEtaRange(-0.9,0.9); | |
2603 | fTrackCuts2->SetRequireSigmaToVertex(kTRUE); | |
e4b0faf2 | 2604 | fTrackCuts2->SetMaxChi2PerClusterTPC(4); |
28ea601a | 2605 | fTrackCuts2->SetMinNClustersTPC(fTPCNClsHad); |
987053ce | 2606 | fTrackCuts2->SetMaxDCAToVertexZ(3.2); |
2607 | fTrackCuts2->SetMaxDCAToVertexXY(2.4); | |
2608 | fTrackCuts2->SetDCAToVertex2D(kTRUE); | |
3c56855b | 2609 | |
987053ce | 2610 | for(Int_t ktracks = 0; ktracks<fVevent->GetNumberOfTracks(); ktracks++){ |
2611 | AliVParticle* VtrackHad = fVevent->GetTrack(ktracks); | |
2612 | if (!VtrackHad) { | |
3c56855b | 2613 | printf("ERROR: Could not receive track %d\n", ktracks); |
2614 | continue; | |
2615 | } | |
987053ce | 2616 | |
2617 | AliVTrack *trackHad = dynamic_cast<AliVTrack*>(VtrackHad); | |
20791315 | 2618 | if(!trackHad) continue; |
987053ce | 2619 | |
2620 | if(IsAODanalysis()) { | |
2621 | AliAODTrack *atrackHad = dynamic_cast<AliAODTrack*>(VtrackHad); | |
7e02de85 | 2622 | if(!atrackHad) continue; |
987053ce | 2623 | if(!atrackHad->TestFilterMask(AliAODTrack::kTrkTPCOnly)) continue; |
2624 | if((!(atrackHad->GetStatus()&AliESDtrack::kITSrefit)|| (!(atrackHad->GetStatus()&AliESDtrack::kTPCrefit)))) continue; | |
28ea601a | 2625 | if(atrackHad->GetTPCNcls() < fTPCNClsHad) continue; |
987053ce | 2626 | } |
2627 | else{ | |
2628 | AliESDtrack *etrackHad = dynamic_cast<AliESDtrack*>(VtrackHad); | |
7e02de85 | 2629 | if(!etrackHad) continue; |
987053ce | 2630 | if(!fTrackCuts2->AcceptTrack(etrackHad)) continue; |
2631 | } | |
2632 | ||
3c56855b | 2633 | if(ktracks == itrack) continue; //do not select the same electron |
a9123929 | 2634 | |
3c56855b | 2635 | Double_t ptHad= -999; |
3c56855b | 2636 | ptHad = trackHad->Pt(); |
a9123929 | 2637 | |
3c56855b | 2638 | if(trackHad->Eta()<-0.9 || trackHad->Eta()>0.9) continue; |
a9123929 | 2639 | // cout << "pt had = " << ptHad <<endl; |
3c56855b | 2640 | |
3c56855b | 2641 | if(ptHad<2) continue; |
2642 | ||
2643 | fHadronPhi->Fill(trackHad->Phi()); | |
2644 | fHadronPhiPt->Fill(trackHad->Phi(),ptHad); | |
2645 | if (trackHad->Eta() >0 && trackHad->Eta() <0.9) fHadronPhiTPChalf->Fill(trackHad->Phi()); | |
2646 | ||
3c56855b | 2647 | fHadronPt->Fill(ptHad); |
2648 | } | |
2649 | } | |
2650 | //_________________________________________ | |
987053ce | 2651 | void AliAnalysisTaskElecHadronCorrel::CheckCentrality(AliVEvent* event, Bool_t ¢ralitypass) |
3c56855b | 2652 | { |
2653 | // Check if event is within the set centrality range. Falls back to V0 centrality determination if no method is set | |
2654 | if (!fkCentralityMethod) AliFatal("No centrality method set! FATAL ERROR!"); | |
2655 | fCentrality = event->GetCentrality()->GetCentralityPercentile(fkCentralityMethod); | |
a9123929 | 2656 | // cout << "Centrality evaluated-------------------------: " << fCentrality <<endl; |
3c56855b | 2657 | |
2658 | if ((fCentrality <= fCentralityMin) || (fCentrality > fCentralityMax)) | |
2659 | { | |
2660 | fCentralityNoPass->Fill(fCentrality); | |
a9123929 | 2661 | // cout << "--------------Fill no pass-------------------------"<<endl; |
3c56855b | 2662 | centralitypass = kFALSE; |
2663 | }else | |
2664 | { | |
2665 | fCentralityPass->Fill(fCentrality); | |
a9123929 | 2666 | // cout << "--------------Fill pass-------------------------"<<endl; |
3c56855b | 2667 | centralitypass = kTRUE; |
2668 | } | |
2669 | ||
2670 | } | |
2671 | //_____________________________________________________________________________ | |
2672 | void AliAnalysisTaskElecHadronCorrel::SetCentralityParameters(Double_t CentralityMin, Double_t CentralityMax, const char* CentralityMethod) | |
2673 | { | |
2674 | // Set a centrality range ]min, max] and define the method to use for centrality selection | |
2675 | fCentralityMin = CentralityMin; | |
2676 | fCentralityMax = CentralityMax; | |
2677 | fkCentralityMethod = CentralityMethod; | |
2678 | } |