23 #include "TObjArray.h" 24 #include "TMatrixDSym.h" 28 #include "Riostream.h" 31 #include "AliAnalysisManager.h" 32 #include "AliInputEventHandler.h" 35 #include "AliAODZDC.h" 36 #include "AliAODVZERO.h" 37 #include "AliAODVertex.h" 40 #include "AliESDEvent.h" 41 #include "AliAODHeader.h" 42 #include "AliAODEvent.h" 43 #include "AliAODTrack.h" 46 #include "AliPIDResponse.h" 47 #include "AliPIDCombined.h" 50 #include "AliVEventHandler.h" 51 #include "AliVEvent.h" 52 #include "AliVTrack.h" 53 #include "AliVParticle.h" 54 #include "AliCentrality.h" 57 #include "AliMultSelection.h" 58 #include "AliAnalysisUtils.h" 59 #include "AliPhysicsSelection.h" 61 #include "AliAnalysisTaskSE.h" 86 fHistTaskConfigParameters(NULL),
87 fHistPileUpCount(NULL),
88 fHistMultSelPUCount(NULL),
89 fHistEtaPtBefore(NULL),
90 fHistEtaPtAfter(NULL),
91 fHistTPCvsGlobalMultBefore(NULL),
92 fHistTPCvsGlobalMultAfter(NULL),
93 fHistTPCdEdxvsPBefore(NULL),
94 fHistTPCdEdxvsPAfter(NULL),
95 fHistTOFBetavsPBefore(NULL),
96 fHistTOFBetavsPAfter(NULL),
97 fHistTOFMassvsPtBefore(NULL),
98 fHistTOFMatchCount(NULL),
99 fHistTPCVsESDTrkBefore(NULL),
100 fHistTPCVsESDTrkAfter(NULL),
101 fHistTPConlyVsCL1Before(NULL),
102 fHistTPConlyVsV0MBefore(NULL),
103 fHistTPConlyVsCL1After(NULL),
104 fHistTPConlyVsV0MAfter(NULL),
105 fHistGlobalVsV0MBefore(NULL),
106 fHistGlobalVsV0MAfter(NULL),
107 fHistRawVsCorrMultFB(NULL),
108 hCentvsTPCmultCuts(NULL),
109 fHV0AEventPlaneVsCent(NULL),
110 fHV0CEventPlaneVsCent(NULL),
111 fHTPCAEventPlaneVsCent(NULL),
112 fHTPCCEventPlaneVsCent(NULL),
113 fHTPCEventPlaneVsCent(NULL),
114 fV0MultChVsRun(NULL),
115 fCentDistBefore(NULL),
116 fCentDistAfter(NULL),
118 fHAvgerageQnV0A(NULL),
119 fHAvgerageQnV0C(NULL),
120 fHCentWeightForRun(NULL),
121 fQAEtaPhiAfterNUA(NULL),
122 fQAEtaPhiAfterNUAPion(NULL),
123 fQAEtaPhiAfterNUAKaon(NULL),
124 fQAEtaPhiAfterNUAProton(NULL),
125 fV0AQ2xVsCentRun(NULL),
126 fV0AQ2yVsCentRun(NULL),
127 fV0CQ2xVsCentRun(NULL),
128 fV0CQ2yVsCentRun(NULL),
129 fV0AQ3xVsCentRun(NULL),
130 fV0AQ3yVsCentRun(NULL),
131 fV0CQ3xVsCentRun(NULL),
132 fV0CQ3yVsCentRun(NULL),
133 fTPCAQ2xVsCentRun(NULL),
134 fTPCAQ2yVsCentRun(NULL),
135 fTPCCQ2xVsCentRun(NULL),
136 fTPCCQ2yVsCentRun(NULL),
137 fTPCAQ3xVsCentRun(NULL),
138 fTPCAQ3yVsCentRun(NULL),
139 fTPCCQ3xVsCentRun(NULL),
140 fTPCCQ3yVsCentRun(NULL),
141 fTPCAQ4xVsCentRun(NULL),
142 fTPCAQ4yVsCentRun(NULL),
143 fTPCCQ4xVsCentRun(NULL),
144 fTPCCQ4yVsCentRun(NULL),
145 fTPCFQ2xVsCentRun(NULL),
146 fTPCFQ2yVsCentRun(NULL),
163 fCentralityPercentMin(0),
164 fCentralityPercentMax(90),
165 fPileUpSlopeParm(3.43),
166 fPileUpConstParm(43),
169 bApplyMCcorr(kFALSE),
170 bV0MGainCorr(kFALSE),
171 bSkipPileUpCut(kFALSE),
172 bFillNUAHistPID(kFALSE),
173 bUseKinkTracks(kFALSE),
174 sPathOfMCFile("/alien"),
176 sCentrEstimator("V0M"),
177 fHistEventCount(NULL)
179 for(
int i=0;i<3;i++){
180 fHistPtwithTPCNsigma[i]=NULL;
181 fHistPtwithTOFmasscut[i]=NULL;
182 fHistPtwithTOFSignal[i]=NULL;
183 fHistTOFnSigmavsPtAfter[i]=NULL;
184 fHistTPCnSigmavsPtAfter[i]=NULL;
185 fHistTPCTOFnSigmavsPtAfter[i]=NULL;
186 fHistTPCdEdxvsPtPIDAfter[i]=NULL;
188 for(
int i=0;i<5;i++){
189 fHCorrectNUApos[i] = NULL;
190 fHCorrectNUAneg[i] = NULL;
192 for(
int i=0;i<5;i++){
193 fHCorrectNUAposPion[i] = NULL;
194 fHCorrectNUAnegPion[i] = NULL;
195 fHCorrectNUAposKaon[i] = NULL;
196 fHCorrectNUAnegKaon[i] = NULL;
197 fHCorrectNUAposProton[i] = NULL;
198 fHCorrectNUAnegProton[i] = NULL;
201 for(
int i=0;i<2;i++){
202 for(
int j=0;j<4;j++){
203 fHist_Corr3p_EP_Norm_PN[i][j] = NULL;
204 fHist_Corr3p_EP_Norm_PP[i][j] = NULL;
205 fHist_Corr3p_EP_Norm_NN[i][j] = NULL;
207 for(
int j=0;j<4;j++) {
208 fHist_Reso2n_EP_Norm_Det[i][j] = NULL;
211 for(
int j=0;j<4;j++){
212 fHist_Corr3p_Pion_EP_Norm_PN[i][j] = NULL;
213 fHist_Corr3p_Pion_EP_Norm_PP[i][j] = NULL;
214 fHist_Corr3p_Pion_EP_Norm_NN[i][j] = NULL;
215 fHist_Corr3p_Kaon_EP_Norm_PN[i][j] = NULL;
216 fHist_Corr3p_Kaon_EP_Norm_PP[i][j] = NULL;
217 fHist_Corr3p_Kaon_EP_Norm_NN[i][j] = NULL;
218 fHist_Corr3p_Proton_EP_Norm_PN[i][j] = NULL;
219 fHist_Corr3p_Proton_EP_Norm_PP[i][j] = NULL;
220 fHist_Corr3p_Proton_EP_Norm_NN[i][j] = NULL;
224 for(
int i=0;i<2;i++){
225 for(
int j=0;j<4;j++){
226 fHist_Corr3p_EP_Refm_PN[i][j] = NULL;
227 fHist_Corr3p_EP_Refm_PP[i][j] = NULL;
228 fHist_Corr3p_EP_Refm_NN[i][j] = NULL;
230 for(
int j=0;j<4;j++) {
231 fHist_Reso2n_EP_Refm_Det[i][j] = NULL;
234 for(
int j=0;j<4;j++){
235 fHist_Corr3p_Pion_EP_Refm_PN[i][j] = NULL;
236 fHist_Corr3p_Pion_EP_Refm_PP[i][j] = NULL;
237 fHist_Corr3p_Pion_EP_Refm_NN[i][j] = NULL;
238 fHist_Corr3p_Kaon_EP_Refm_PN[i][j] = NULL;
239 fHist_Corr3p_Kaon_EP_Refm_PP[i][j] = NULL;
240 fHist_Corr3p_Kaon_EP_Refm_NN[i][j] = NULL;
241 fHist_Corr3p_Proton_EP_Refm_PN[i][j] = NULL;
242 fHist_Corr3p_Proton_EP_Refm_PP[i][j] = NULL;
243 fHist_Corr3p_Proton_EP_Refm_NN[i][j] = NULL;
261 for(
int i=0;i<2;i++){
262 for(
int j=0;j<6;j++){
263 fHist_Corr3p_pTSum_EP_V0A_PN[i][j] = NULL;
264 fHist_Corr3p_pTSum_EP_V0A_PP[i][j] = NULL;
265 fHist_Corr3p_pTSum_EP_V0A_NN[i][j] = NULL;
266 fHist_Corr3p_pTSum_EP_V0C_PN[i][j] = NULL;
267 fHist_Corr3p_pTSum_EP_V0C_PP[i][j] = NULL;
268 fHist_Corr3p_pTSum_EP_V0C_NN[i][j] = NULL;
270 fHist_Corr3p_pTDiff_EP_V0A_PN[i][j] = NULL;
271 fHist_Corr3p_pTDiff_EP_V0A_PP[i][j] = NULL;
272 fHist_Corr3p_pTDiff_EP_V0A_NN[i][j] = NULL;
273 fHist_Corr3p_pTDiff_EP_V0C_PN[i][j] = NULL;
274 fHist_Corr3p_pTDiff_EP_V0C_PP[i][j] = NULL;
275 fHist_Corr3p_pTDiff_EP_V0C_NN[i][j] = NULL;
277 fHist_Corr3p_EtaDiff_EP_V0A_PN[i][j] = NULL;
278 fHist_Corr3p_EtaDiff_EP_V0A_PP[i][j] = NULL;
279 fHist_Corr3p_EtaDiff_EP_V0A_NN[i][j] = NULL;
280 fHist_Corr3p_EtaDiff_EP_V0C_PN[i][j] = NULL;
281 fHist_Corr3p_EtaDiff_EP_V0C_PP[i][j] = NULL;
282 fHist_Corr3p_EtaDiff_EP_V0C_NN[i][j] = NULL;
286 for(
int i=0;i<2;i++){
287 for(
int j=0;j<6;j++){
289 fHist_Corr3p_Pion_pTSum_EP_V0A_PN[i][j] = NULL;
290 fHist_Corr3p_Pion_pTSum_EP_V0A_PP[i][j] = NULL;
291 fHist_Corr3p_Pion_pTSum_EP_V0A_NN[i][j] = NULL;
292 fHist_Corr3p_Pion_pTSum_EP_V0C_PN[i][j] = NULL;
293 fHist_Corr3p_Pion_pTSum_EP_V0C_PP[i][j] = NULL;
294 fHist_Corr3p_Pion_pTSum_EP_V0C_NN[i][j] = NULL;
296 fHist_Corr3p_Pion_pTDiff_EP_V0A_PN[i][j] = NULL;
297 fHist_Corr3p_Pion_pTDiff_EP_V0A_PP[i][j] = NULL;
298 fHist_Corr3p_Pion_pTDiff_EP_V0A_NN[i][j] = NULL;
299 fHist_Corr3p_Pion_pTDiff_EP_V0C_PN[i][j] = NULL;
300 fHist_Corr3p_Pion_pTDiff_EP_V0C_PP[i][j] = NULL;
301 fHist_Corr3p_Pion_pTDiff_EP_V0C_NN[i][j] = NULL;
303 fHist_Corr3p_Pion_EtaDiff_EP_V0A_PN[i][j] = NULL;
304 fHist_Corr3p_Pion_EtaDiff_EP_V0A_PP[i][j] = NULL;
305 fHist_Corr3p_Pion_EtaDiff_EP_V0A_NN[i][j] = NULL;
306 fHist_Corr3p_Pion_EtaDiff_EP_V0C_PN[i][j] = NULL;
307 fHist_Corr3p_Pion_EtaDiff_EP_V0C_PP[i][j] = NULL;
308 fHist_Corr3p_Pion_EtaDiff_EP_V0C_NN[i][j] = NULL;
310 fHist_Corr3p_Kaon_pTSum_EP_V0A_PN[i][j] = NULL;
311 fHist_Corr3p_Kaon_pTSum_EP_V0A_PP[i][j] = NULL;
312 fHist_Corr3p_Kaon_pTSum_EP_V0A_NN[i][j] = NULL;
313 fHist_Corr3p_Kaon_pTSum_EP_V0C_PN[i][j] = NULL;
314 fHist_Corr3p_Kaon_pTSum_EP_V0C_PP[i][j] = NULL;
315 fHist_Corr3p_Kaon_pTSum_EP_V0C_NN[i][j] = NULL;
317 fHist_Corr3p_Kaon_pTDiff_EP_V0A_PN[i][j] = NULL;
318 fHist_Corr3p_Kaon_pTDiff_EP_V0A_PP[i][j] = NULL;
319 fHist_Corr3p_Kaon_pTDiff_EP_V0A_NN[i][j] = NULL;
320 fHist_Corr3p_Kaon_pTDiff_EP_V0C_PN[i][j] = NULL;
321 fHist_Corr3p_Kaon_pTDiff_EP_V0C_PP[i][j] = NULL;
322 fHist_Corr3p_Kaon_pTDiff_EP_V0C_NN[i][j] = NULL;
324 fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PN[i][j] = NULL;
325 fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PP[i][j] = NULL;
326 fHist_Corr3p_Kaon_EtaDiff_EP_V0A_NN[i][j] = NULL;
327 fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PN[i][j] = NULL;
328 fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PP[i][j] = NULL;
329 fHist_Corr3p_Kaon_EtaDiff_EP_V0C_NN[i][j] = NULL;
331 fHist_Corr3p_Proton_pTSum_EP_V0A_PN[i][j] = NULL;
332 fHist_Corr3p_Proton_pTSum_EP_V0A_PP[i][j] = NULL;
333 fHist_Corr3p_Proton_pTSum_EP_V0A_NN[i][j] = NULL;
334 fHist_Corr3p_Proton_pTSum_EP_V0C_PN[i][j] = NULL;
335 fHist_Corr3p_Proton_pTSum_EP_V0C_PP[i][j] = NULL;
336 fHist_Corr3p_Proton_pTSum_EP_V0C_NN[i][j] = NULL;
338 fHist_Corr3p_Proton_pTDiff_EP_V0A_PN[i][j] = NULL;
339 fHist_Corr3p_Proton_pTDiff_EP_V0A_PP[i][j] = NULL;
340 fHist_Corr3p_Proton_pTDiff_EP_V0A_NN[i][j] = NULL;
341 fHist_Corr3p_Proton_pTDiff_EP_V0C_PN[i][j] = NULL;
342 fHist_Corr3p_Proton_pTDiff_EP_V0C_PP[i][j] = NULL;
343 fHist_Corr3p_Proton_pTDiff_EP_V0C_NN[i][j] = NULL;
345 fHist_Corr3p_Proton_EtaDiff_EP_V0A_PN[i][j] = NULL;
346 fHist_Corr3p_Proton_EtaDiff_EP_V0A_PP[i][j] = NULL;
347 fHist_Corr3p_Proton_EtaDiff_EP_V0A_NN[i][j] = NULL;
348 fHist_Corr3p_Proton_EtaDiff_EP_V0C_PN[i][j] = NULL;
349 fHist_Corr3p_Proton_EtaDiff_EP_V0C_PP[i][j] = NULL;
350 fHist_Corr3p_Proton_EtaDiff_EP_V0C_NN[i][j] = NULL;
356 for(
int i=0;i<4;i++){
357 for(
int j=0;j<5;j++){
358 fHist3DEtaPhiVz_Pos_Run[i][j]=NULL;
359 fHist3DEtaPhiVz_Neg_Run[i][j]=NULL;
362 for(
int i=0;i<10;i++){
363 fFB_Efficiency_Cent[i] = NULL;
364 fFB_Efficiency_Pion_Cent[i] = NULL;
365 fFB_Efficiency_Kaon_Cent[i] = NULL;
367 fFB_Efficiency_Proton_Pos_Cent[i] = NULL;
368 fFB_Efficiency_Proton_Neg_Cent[i] = NULL;
372 DefineInput(0,TChain::Class());
373 DefineOutput(1,TList::Class());
382 fMultSelection(NULL),
389 fHistTaskConfigParameters(NULL),
390 fHistPileUpCount(NULL),
391 fHistMultSelPUCount(NULL),
392 fHistEtaPtBefore(NULL),
393 fHistEtaPtAfter(NULL),
394 fHistTPCvsGlobalMultBefore(NULL),
395 fHistTPCvsGlobalMultAfter(NULL),
396 fHistTPCdEdxvsPBefore(NULL),
397 fHistTPCdEdxvsPAfter(NULL),
398 fHistTOFBetavsPBefore(NULL),
399 fHistTOFBetavsPAfter(NULL),
400 fHistTOFMassvsPtBefore(NULL),
401 fHistTOFMatchCount(NULL),
402 fHistTPCVsESDTrkBefore(NULL),
403 fHistTPCVsESDTrkAfter(NULL),
404 fHistTPConlyVsCL1Before(NULL),
405 fHistTPConlyVsV0MBefore(NULL),
406 fHistTPConlyVsCL1After(NULL),
407 fHistTPConlyVsV0MAfter(NULL),
408 fHistGlobalVsV0MBefore(NULL),
409 fHistGlobalVsV0MAfter(NULL),
410 fHistRawVsCorrMultFB(NULL),
411 hCentvsTPCmultCuts(NULL),
412 fHV0AEventPlaneVsCent(NULL),
413 fHV0CEventPlaneVsCent(NULL),
414 fHTPCAEventPlaneVsCent(NULL),
415 fHTPCCEventPlaneVsCent(NULL),
416 fHTPCEventPlaneVsCent(NULL),
417 fV0MultChVsRun(NULL),
418 fCentDistBefore(NULL),
419 fCentDistAfter(NULL),
421 fHAvgerageQnV0A(NULL),
422 fHAvgerageQnV0C(NULL),
423 fHCentWeightForRun(NULL),
424 fQAEtaPhiAfterNUA(NULL),
425 fQAEtaPhiAfterNUAPion(NULL),
426 fQAEtaPhiAfterNUAKaon(NULL),
427 fQAEtaPhiAfterNUAProton(NULL),
428 fV0AQ2xVsCentRun(NULL),
429 fV0AQ2yVsCentRun(NULL),
430 fV0CQ2xVsCentRun(NULL),
431 fV0CQ2yVsCentRun(NULL),
432 fV0AQ3xVsCentRun(NULL),
433 fV0AQ3yVsCentRun(NULL),
434 fV0CQ3xVsCentRun(NULL),
435 fV0CQ3yVsCentRun(NULL),
436 fTPCAQ2xVsCentRun(NULL),
437 fTPCAQ2yVsCentRun(NULL),
438 fTPCCQ2xVsCentRun(NULL),
439 fTPCCQ2yVsCentRun(NULL),
440 fTPCAQ3xVsCentRun(NULL),
441 fTPCAQ3yVsCentRun(NULL),
442 fTPCCQ3xVsCentRun(NULL),
443 fTPCCQ3yVsCentRun(NULL),
444 fTPCAQ4xVsCentRun(NULL),
445 fTPCAQ4yVsCentRun(NULL),
446 fTPCCQ4xVsCentRun(NULL),
447 fTPCCQ4yVsCentRun(NULL),
448 fTPCFQ2xVsCentRun(NULL),
449 fTPCFQ2yVsCentRun(NULL),
466 fCentralityPercentMin(0),
467 fCentralityPercentMax(90),
468 fPileUpSlopeParm(3.43),
469 fPileUpConstParm(43),
472 bApplyMCcorr(kFALSE),
473 bV0MGainCorr(kFALSE),
474 bSkipPileUpCut(kFALSE),
475 bFillNUAHistPID(kFALSE),
476 bUseKinkTracks(kFALSE),
477 sPathOfMCFile(
"/alien"),
479 sCentrEstimator(
"V0M"),
480 fHistEventCount(NULL)
482 for(
int i=0;i<3;i++){
491 for(
int i=0;i<5;i++){
495 for(
int i=0;i<5;i++){
504 for(
int i=0;i<2;i++){
505 for(
int j=0;j<4;j++){
510 for(
int j=0;j<4;j++) {
514 for(
int j=0;j<4;j++){
527 for(
int i=0;i<2;i++){
528 for(
int j=0;j<4;j++){
533 for(
int j=0;j<4;j++) {
537 for(
int j=0;j<4;j++){
570 for(
int i=0;i<2;i++){
571 for(
int j=0;j<6;j++){
595 for(
int i=0;i<2;i++){
596 for(
int j=0;j<6;j++){
665 for(
int i=0;i<4;i++){
666 for(
int j=0;j<5;j++){
671 for(
int i=0;i<10;i++){
734 AliInputEventHandler *inputHandler=
dynamic_cast<AliInputEventHandler*
>(mgr->GetInputEventHandler());
735 if (!inputHandler) { printf(
"\n ...Input handler missing!!!...\n");
return; }
752 TGrid::Connect(
"alien://");
776 Int_t gMaxTPCFB1mult = 0;
777 Int_t gMaxGlobalmult = 0;
778 Int_t gMaxTPCcorrmult = 0;
779 Int_t gMaxESDtracks = 0;
780 Int_t nBinRefMult = 200;
781 Int_t nRefMultMax = 4000;
786 gMaxGlobalmult = 200;
787 gMaxTPCFB1mult = 200;
788 gMaxTPCcorrmult = 500;
789 gMaxESDtracks = 1000;
795 gMaxGlobalmult = 400;
796 gMaxTPCFB1mult = 400;
797 gMaxTPCcorrmult = 500;
798 gMaxESDtracks = 2000;
804 gMaxGlobalmult = 4000;
805 gMaxTPCFB1mult = 4000;
806 gMaxTPCcorrmult = 5000;
807 gMaxESDtracks = 20000;
816 fHistTPCVsESDTrkBefore =
new TH2F(
"fHistTPCVsESDTrkBefore",
"Before; TPC1; ESD trk",500,0,gMaxTPCcorrmult,200,0,gMaxESDtracks);
818 fHistTPCVsESDTrkAfter =
new TH2F(
"fHistTPCVsESDTrkAfter",
" After; TPC1; ESD trk",500,0,gMaxTPCcorrmult,200,0,gMaxESDtracks);
841 fHistRawVsCorrMultFB =
new TH2F(
"fHistRawVsCorrMultFB",Form(
"FB%d;Mult_{raw};Mult_{corr}",
fFilterBit),gMaxTPCFB1mult,0,gMaxTPCFB1mult,gMaxTPCcorrmult,0,gMaxTPCcorrmult);
892 fCentDistBefore =
new TH1F(
"fCentDistBefore",
"no Cut; Cent (%); Events ",100,0,100);
895 fCentDistAfter =
new TH1F(
"fCentDistAfter",
"with Cut; Cent (%); Events ",100,0,100);
910 Double_t centRange[11] = {0,5,10,20,30,40,50,60,70,80,90};
914 for(
int i=0;i<2;i++){
916 for(
int j=0;j<4;j++){
918 fHist_Corr3p_EP_Norm_PN[i][j] =
new TProfile(Form(
"fHist_Corr3p_EP_Norm_PosNeg_Mag%d_Det%d",i,j+1),Form(
"US, #Psi_{2} %d",j),10,centRange,
"");
921 fHist_Corr3p_EP_Norm_PP[i][j] =
new TProfile(Form(
"fHist_Corr3p_EP_Norm_PosPos_Mag%d_Det%d",i,j+1),Form(
"P-P, #Psi_{2} %d",j),10,centRange,
"");
924 fHist_Corr3p_EP_Norm_NN[i][j] =
new TProfile(Form(
"fHist_Corr3p_EP_Norm_NegNeg_Mag%d_Det%d",i,j+1),Form(
"N-N, #Psi_{2}, %d",j),10,centRange,
"");
929 for(
int j=0;j<4;j++){
931 fHist_Reso2n_EP_Norm_Det[i][j] =
new TProfile(Form(
"fHist_Reso2n_EP_Norm_Mag%d_DetComb%d",i,j+1),
"Event plane Resolution",10,centRange,
"");
936 for(
int j=0;j<4;j++){
938 fHist_Corr3p_Pion_EP_Norm_PN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Pion_EP_Norm_PosNeg_Mag%d_Det%d",i,j+1),Form(
"US, #Psi_{2} %d",j),10,centRange,
"");
941 fHist_Corr3p_Pion_EP_Norm_PP[i][j] =
new TProfile(Form(
"fHist_Corr3p_Pion_EP_Norm_PosPos_Mag%d_Det%d",i,j+1),Form(
"P-P, #Psi_{2} %d",j),10,centRange,
"");
944 fHist_Corr3p_Pion_EP_Norm_NN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Pion_EP_Norm_NegNeg_Mag%d_Det%d",i,j+1),Form(
"N-N, #Psi_{2}, %d",j),10,centRange,
"");
948 fHist_Corr3p_Kaon_EP_Norm_PN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Kaon_EP_Norm_PosNeg_Mag%d_Det%d",i,j+1),Form(
"US, #Psi_{2} %d",j),10,centRange,
"");
951 fHist_Corr3p_Kaon_EP_Norm_PP[i][j] =
new TProfile(Form(
"fHist_Corr3p_Kaon_EP_Norm_PosPos_Mag%d_Det%d",i,j+1),Form(
"P-P, #Psi_{2} %d",j),10,centRange,
"");
954 fHist_Corr3p_Kaon_EP_Norm_NN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Kaon_EP_Norm_NegNeg_Mag%d_Det%d",i,j+1),Form(
"N-N, #Psi_{2}, %d",j),10,centRange,
"");
958 fHist_Corr3p_Proton_EP_Norm_PN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Proton_EP_Norm_PosNeg_Mag%d_Det%d",i,j+1),Form(
"US, #Psi_{2} %d",j),10,centRange,
"");
961 fHist_Corr3p_Proton_EP_Norm_PP[i][j] =
new TProfile(Form(
"fHist_Corr3p_Proton_EP_Norm_PosPos_Mag%d_Det%d",i,j+1),Form(
"P-P, #Psi_{2} %d",j),10,centRange,
"");
964 fHist_Corr3p_Proton_EP_Norm_NN[i][j] =
new TProfile(Form(
"fHist_Corr3p_Proton_EP_Norm_NegNeg_Mag%d_Det%d",i,j+1),Form(
"N-N, #Psi_{2}, %d",j),10,centRange,
"");
1048 Double_t pTRange[21] = {0.0,0.2,0.4,0.6,0.8,1.0,1.2,1.4,1.6,1.8,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0};
1052 for(
Int_t i=0;i<2;i++){
1053 for(
Int_t j=0;j<6;j++){
1054 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0A_PN_Mag%d_Cent%d",i,j);
1055 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1060 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0A_PP_Mag%d_Cent%d",i,j);
1061 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1066 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0A_NN_Mag%d_Cent%d",i,j);
1067 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1072 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0C_PN_Mag%d_Cent%d",i,j);
1073 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1078 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0C_PP_Mag%d_Cent%d",i,j);
1079 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1084 sprintf(name,
"fHist_Corr3p_pTSum_EP_V0C_NN_Mag%d_Cent%d",i,j);
1085 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1092 for(
Int_t i=0;i<2;i++){
1093 for(
Int_t j=0;j<6;j++){
1094 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1095 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1100 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1101 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1106 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1107 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1112 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1113 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1118 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1119 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1124 sprintf(name,
"fHist_Corr3p_pTDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1125 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1135 for(
Int_t i=0;i<2;i++){
1136 for(
Int_t j=0;j<6;j++){
1137 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1138 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1143 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1144 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1149 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1150 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1155 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1156 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1161 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1162 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1167 sprintf(name,
"fHist_Corr3p_EtaDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1168 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1176 for(
Int_t i=0;i<2;i++){
1177 for(
Int_t j=0;j<6;j++){
1178 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0A_PN_Mag%d_Cent%d",i,j);
1179 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1184 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0A_PP_Mag%d_Cent%d",i,j);
1185 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1190 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0A_NN_Mag%d_Cent%d",i,j);
1191 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1196 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0C_PN_Mag%d_Cent%d",i,j);
1197 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1202 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0C_PP_Mag%d_Cent%d",i,j);
1203 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1208 sprintf(name,
"fHist_Corr3p_Pion_pTSum_EP_V0C_NN_Mag%d_Cent%d",i,j);
1209 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1216 for(
Int_t i=0;i<2;i++){
1217 for(
Int_t j=0;j<6;j++){
1218 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1219 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1224 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1225 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1230 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1231 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1236 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1237 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1242 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1243 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1248 sprintf(name,
"fHist_Corr3p_Pion_pTDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1249 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1258 for(
Int_t i=0;i<2;i++){
1259 for(
Int_t j=0;j<6;j++){
1260 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1261 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1266 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1267 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1272 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1273 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1278 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1279 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1284 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1285 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1290 sprintf(name,
"fHist_Corr3p_Pion_EtaDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1291 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1299 for(
Int_t i=0;i<2;i++){
1300 for(
Int_t j=0;j<6;j++){
1301 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0A_PN_Mag%d_Cent%d",i,j);
1302 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1307 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0A_PP_Mag%d_Cent%d",i,j);
1308 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1313 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0A_NN_Mag%d_Cent%d",i,j);
1314 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1319 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0C_PN_Mag%d_Cent%d",i,j);
1320 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1325 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0C_PP_Mag%d_Cent%d",i,j);
1326 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1331 sprintf(name,
"fHist_Corr3p_Kaon_pTSum_EP_V0C_NN_Mag%d_Cent%d",i,j);
1332 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1339 for(
Int_t i=0;i<2;i++){
1340 for(
Int_t j=0;j<6;j++){
1341 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1342 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1347 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1348 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1353 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1354 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1359 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1360 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1365 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1366 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1371 sprintf(name,
"fHist_Corr3p_Kaon_pTDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1372 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1381 for(
Int_t i=0;i<2;i++){
1382 for(
Int_t j=0;j<6;j++){
1383 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1384 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1389 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1390 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1395 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1396 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1401 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1402 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1407 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1408 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1413 sprintf(name,
"fHist_Corr3p_Kaon_EtaDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1414 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1422 for(
Int_t i=0;i<2;i++){
1423 for(
Int_t j=0;j<6;j++){
1424 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0A_PN_Mag%d_Cent%d",i,j);
1425 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1430 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0A_PP_Mag%d_Cent%d",i,j);
1431 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1436 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0A_NN_Mag%d_Cent%d",i,j);
1437 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1442 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0C_PN_Mag%d_Cent%d",i,j);
1443 sprintf(title,
"PN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1448 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0C_PP_Mag%d_Cent%d",i,j);
1449 sprintf(title,
"PP 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1454 sprintf(name,
"fHist_Corr3p_Proton_pTSum_EP_V0C_NN_Mag%d_Cent%d",i,j);
1455 sprintf(title,
"NN 3p vs (pT1+pT2)/2, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1462 for(
Int_t i=0;i<2;i++){
1463 for(
Int_t j=0;j<6;j++){
1464 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1465 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1470 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1471 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1476 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1477 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1482 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1483 sprintf(title,
"PN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1488 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1489 sprintf(title,
"PP 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1494 sprintf(name,
"fHist_Corr3p_Proton_pTDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1495 sprintf(title,
"NN 3p vs |pT1-pT2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1504 for(
Int_t i=0;i<2;i++){
1505 for(
Int_t j=0;j<6;j++){
1506 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0A_PN_Mag%d_Cent%d",i,j);
1507 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1512 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0A_PP_Mag%d_Cent%d",i,j);
1513 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1518 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0A_NN_Mag%d_Cent%d",i,j);
1519 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1524 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0C_PN_Mag%d_Cent%d",i,j);
1525 sprintf(title,
"PN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1530 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0C_PP_Mag%d_Cent%d",i,j);
1531 sprintf(title,
"PP 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1536 sprintf(name,
"fHist_Corr3p_Proton_EtaDiff_EP_V0C_NN_Mag%d_Cent%d",i,j);
1537 sprintf(title,
"NN 3p vs |Eta1-Eta2|, Cent %2.0f-%2.0f",centRange[i],centRange[i+1]);
1557 fListNUACalib->SetName(
"fListNUACalib");
1558 fListNUACalib->SetOwner(kTRUE);
1582 fQAEtaPhiAfterNUA =
new TH2F(
"fQAEtaPhiAfterNUA",
"eta vs phi with NUA corr",50,0,6.283185,16,-0.8,0.8);
1611 fV0AQ2xVsCentRun =
new TProfile(
"fV0ACos2nVsCentRun",Form(
"<Cos2> vs cent (%d)",isCorr),90,0,90,
"");
1613 fV0AQ2yVsCentRun =
new TProfile(
"fV0ASin2nVsCentRun",Form(
"<Sin2> vs cent (%d)",isCorr),90,0,90,
"");
1615 fV0CQ2xVsCentRun =
new TProfile(
"fV0CCos2nVsCentRun",Form(
"<Cos2> vs cent (%d)",isCorr),90,0,90,
"");
1617 fV0CQ2yVsCentRun =
new TProfile(
"fV0CSin2nVsCentRun",Form(
"<Sin2> vs cent (%d)",isCorr),90,0,90,
"");
1620 fV0AQ3xVsCentRun =
new TProfile(
"fV0ACos3nVsCentRun",Form(
"<Cos3> vs cent (%d)",isCorr),90,0,90,
"");
1622 fV0AQ3yVsCentRun =
new TProfile(
"fV0ASin3nVsCentRun",Form(
"<Sin3> vs cent (%d)",isCorr),90,0,90,
"");
1624 fV0CQ3xVsCentRun =
new TProfile(
"fV0CCos3nVsCentRun",Form(
"<Cos3> vs cent (%d)",isCorr),90,0,90,
"");
1626 fV0CQ3yVsCentRun =
new TProfile(
"fV0CSin3nVsCentRun",Form(
"<Sin3> vs cent (%d)",isCorr),90,0,90,
"");
1631 cout<<
"\n =========> NUA file found for NUA correction <========== \n";
1635 fTPCAQ2xVsCentRun =
new TProfile(
"fTPCACos2nVsCentRun",Form(
"<Cos2> vs cent (%d)",isCorr),90,0,90,
"");
1637 fTPCAQ2yVsCentRun =
new TProfile(
"fTPCASin2nVsCentRun",Form(
"<Sin2> vs cent (%d)",isCorr),90,0,90,
"");
1639 fTPCCQ2xVsCentRun =
new TProfile(
"fTPCCCos2nVsCentRun",Form(
"<Cos2> vs cent (%d)",isCorr),90,0,90,
"");
1641 fTPCCQ2yVsCentRun =
new TProfile(
"fTPCCSin2nVsCentRun",Form(
"<Sin2> vs cent (%d)",isCorr),90,0,90,
"");
1644 fTPCAQ3xVsCentRun =
new TProfile(
"fTPCACos3nVsCentRun",Form(
"<Cos3> vs cent (%d)",isCorr),90,0,90,
"");
1646 fTPCAQ3yVsCentRun =
new TProfile(
"fTPCASin3nVsCentRun",Form(
"<Sin3> vs cent (%d)",isCorr),90,0,90,
"");
1648 fTPCCQ3xVsCentRun =
new TProfile(
"fTPCCCos3nVsCentRun",Form(
"<Cos3> vs cent (%d)",isCorr),90,0,90,
"");
1650 fTPCCQ3yVsCentRun =
new TProfile(
"fTPCCSin3nVsCentRun",Form(
"<Sin3> vs cent (%d)",isCorr),90,0,90,
"");
1653 fTPCAQ4xVsCentRun =
new TProfile(
"fTPCACos4nVsCentRun",Form(
"<Cos4> vs cent (%d)",isCorr),90,0,90,
"");
1655 fTPCAQ4yVsCentRun =
new TProfile(
"fTPCASin4nVsCentRun",Form(
"<Sin4> vs cent (%d)",isCorr),90,0,90,
"");
1657 fTPCCQ4xVsCentRun =
new TProfile(
"fTPCCCos4nVsCentRun",Form(
"<Cos4> vs cent (%d)",isCorr),90,0,90,
"");
1659 fTPCCQ4yVsCentRun =
new TProfile(
"fTPCCSin4nVsCentRun",Form(
"<Sin4> vs cent (%d)",isCorr),90,0,90,
"");
1663 fTPCFQ2xVsCentRun =
new TProfile(
"fTPCFQ2xVsCentRun",Form(
"<Cos2> vs cent (%d)",isCorr),90,0,90,
"");
1665 fTPCFQ2yVsCentRun =
new TProfile(
"fTPCFQ2yVsCentRun",Form(
"<Sin2> vs cent (%d)",isCorr),90,0,90,
"");
1674 Int_t gCentForNUA[6] = {0,5,10,20,40,90};
1678 for(
int i=0;i<4;i++){
1679 for(
int j=0;j<5;j++){
1680 sprintf(name,
"fHistEtaPhiVz_%d_Pos_Cent%d_Run%d",i,j,1);
1681 sprintf(title,
"eta,phi,Vz %dPos, Cent%d-%d, FB %d",i,gCentForNUA[j],gCentForNUA[j+1],
fFilterBit);
1682 fHist3DEtaPhiVz_Pos_Run[i][j] =
new TH3F(name,title,10,-10,10,50,0,6.283185,16,-0.8,0.8);
1685 sprintf(name,
"fHistEtaPhiVz_%d_Neg_Cent%d_Run%d",i,j,1);
1686 sprintf(title,
"eta,phi,Vz %dNeg, Cent%d-%d, FB %d",i,gCentForNUA[j],gCentForNUA[j+1],
fFilterBit);
1687 fHist3DEtaPhiVz_Neg_Run[i][j] =
new TH3F(name,title,10,-10,10,50,0,6.283185,16,-0.8,0.8);
1734 if(!(
fESD ||
fAOD)){ printf(
"ERROR: fESD & fAOD not available\n");
return; }
1736 fVevent =
dynamic_cast<AliVEvent*
>(InputEvent());
1738 if (!
fVevent) { printf(
"ERROR: fVevent not available\n");
return; }
1751 AliInputEventHandler *inputHandler=
dynamic_cast<AliInputEventHandler*
>(mgr->GetInputEventHandler());
1752 if(inputHandler)
fPIDResponse=inputHandler->GetPIDResponse();
1754 printf(
"\n\n...... PIDResponse object not found..... \n\n");
1764 const AliVVertex *pVtx =
fVevent->GetPrimaryVertex();
1767 pVtxZ = pVtx->GetZ();
1771 if(pVtxZ<fMinVzCut || pVtxZ>
fMaxVzCut )
return;
1786 AliCentrality* Alicentr = ((AliVAODHeader*)
fAOD->GetHeader())->GetCentralityP();
1788 fMultSelection = (AliMultSelection*) InputEvent()->FindListObject(
"MultSelection");
1789 if(!
fMultSelection) { printf(
"\n\n **WARNING** \n::UserExec() AliMultSelection object not found.\n\n"); exit(1); }
1795 centrV0M = Alicentr->GetCentralityPercentile(
"V0M");
1796 centrCL1 = Alicentr->GetCentralityPercentile(
"CL1");
1799 centrality = centrV0M;
1802 centrality = centrCL1;
1805 centrality = Alicentr->GetCentralityPercentile(
"V0C");
1808 centrality = Alicentr->GetCentralityPercentile(
"V0A");
1811 centrality = Alicentr->GetCentralityPercentile(
"TRK");
1821 centrality = centrV0M;
1824 centrality = centrCL1;
1875 Int_t ntracks=
fAOD->GetNumberOfTracks();
1876 if(ntracks<2)
return;
1886 Int_t cent10bin = -1;
1890 if(centrality<5.0) {
1893 else if(centrality>=5.0 && centrality<10){
1896 else if(centrality>=10.0) {
1897 cent10bin = abs(centrality/10.0)+1;
1905 if(centrality<5.0) {
1908 else if(centrality>=5.0 && centrality<10){
1911 else if(centrality>=10.0 && centrality<20) {
1914 else if(centrality>=20 && centrality<40){
1917 else if(centrality>=40){
1929 const Int_t QAindex = (fMagField > 0) ? 1 : 0;
1996 Int_t isPileup =
fAOD->IsPileupFromSPD(3);
2006 else if(((AliAODHeader*)
fAOD->GetHeader())->GetRefMultiplicityComb08() < 0) {
2010 else if(
fAOD->IsIncompleteDAQ()) {
2014 else if(fabs(centrV0M-centrCL1)> 5.0) {
2021 const AliAODVertex* vtTrc =
fAOD->GetPrimaryVertex();
2022 const AliAODVertex* vtSPD =
fAOD->GetPrimaryVertexSPD();
2024 if(vtTrc->GetNContributors() < 2 || vtSPD->GetNContributors()<1) {
2029 double covTrc[6], covSPD[6];
2030 vtTrc->GetCovarianceMatrix(covTrc);
2031 vtSPD->GetCovarianceMatrix(covSPD);
2033 double dz = vtTrc->GetZ() - vtSPD->GetZ();
2035 double errTot = TMath::Sqrt(covTrc[5]+covSPD[5]);
2036 double errTrc = TMath::Sqrt(covTrc[5]);
2037 double nsigTot = dz/errTot;
2038 double nsigTrc = dz/errTrc;
2040 if(TMath::Abs(dz)>0.2 || TMath::Abs(nsigTot)>10 || TMath::Abs(nsigTrc)>20) {
2047 Float_t nSigPionTPC[40000] = {0.,};
2048 Float_t nSigKaonTPC[40000] = {0.,};
2049 Float_t nSigProtonTPC[40000] = {0.,};
2050 Float_t nSigPionTOF[40000] = {0.,};
2051 Float_t nSigKaonTOF[40000] = {0.,};
2052 Float_t nSigProtonTOF[40000] = {0.,};
2054 if(ntracks > 40000)
return;
2084 Int_t multEtaNeg, multEtaPos, multEtaFull;
2085 Double_t SumWEtaNeg, SumWEtaPos, SumWEtaFull;
2106 for(
Int_t iTrack = 0; iTrack < ntracks; iTrack++) {
2108 AliAODTrack* AODtrack =
dynamic_cast<AliAODTrack*
>(
fVevent->GetTrack(iTrack));
2109 if(!AODtrack)
continue;
2123 nSigPionTPC[iTrack] =
fPIDResponse->NumberOfSigmasTPC(AODtrack, AliPID::kPion);
2124 nSigKaonTPC[iTrack] =
fPIDResponse->NumberOfSigmasTPC(AODtrack, AliPID::kKaon);
2127 nSigPionTOF[iTrack] =
fPIDResponse->NumberOfSigmasTOF(AODtrack, AliPID::kPion);
2128 nSigKaonTOF[iTrack] =
fPIDResponse->NumberOfSigmasTOF(AODtrack, AliPID::kKaon);
2134 ptTrk = AODtrack->Pt();
2135 etaTrk = AODtrack->Eta();
2136 Chi2Trk = AODtrack->Chi2perNDF();
2137 nClustTPC = AODtrack->GetTPCNcls();
2147 phiTrk = AODtrack->Phi();
2148 dEdx = AODtrack->GetDetPid()->GetTPCsignal();
2149 ChTrk = AODtrack->Charge();
2162 RefMultCorrFB += ptWgtMC;
2184 sumTPCQn2x[0] += WgtNUA*TMath::Cos(
gPsiN*phiTrk);
2185 sumTPCQn2y[0] += WgtNUA*TMath::Sin(
gPsiN*phiTrk);
2186 sumTPCQn3x[0] += WgtNUA*TMath::Cos(3*phiTrk);
2187 sumTPCQn3y[0] += WgtNUA*TMath::Sin(3*phiTrk);
2188 sumTPCQn4x[0] += WgtNUA*TMath::Cos(4*phiTrk);
2189 sumTPCQn4y[0] += WgtNUA*TMath::Sin(4*phiTrk);
2191 SumWEtaNeg += WgtNUA;
2193 else if(etaTrk > 0.05){
2194 sumTPCQn2x[1] += WgtNUA*TMath::Cos(
gPsiN*phiTrk);
2195 sumTPCQn2y[1] += WgtNUA*TMath::Sin(
gPsiN*phiTrk);
2196 sumTPCQn3x[1] += WgtNUA*TMath::Cos(3*phiTrk);
2197 sumTPCQn3y[1] += WgtNUA*TMath::Sin(3*phiTrk);
2198 sumTPCQn4x[1] += WgtNUA*TMath::Cos(4*phiTrk);
2199 sumTPCQn4y[1] += WgtNUA*TMath::Sin(4*phiTrk);
2201 SumWEtaPos += WgtNUA;
2203 sumTPCQn2x[3] += WgtNUA*TMath::Cos(
gPsiN*phiTrk);
2204 sumTPCQn2y[3] += WgtNUA*TMath::Sin(
gPsiN*phiTrk);
2205 sumTPCQn3x[3] += WgtNUA*TMath::Cos(3*phiTrk);
2206 sumTPCQn3y[3] += WgtNUA*TMath::Sin(3*phiTrk);
2207 sumTPCQn4x[3] += WgtNUA*TMath::Cos(4*phiTrk);
2208 sumTPCQn4y[3] += WgtNUA*TMath::Sin(4*phiTrk);
2210 SumWEtaFull += WgtNUA;
2215 if(AODtrack->TestFilterBit(128)) multTPCAll++;
2217 if(!(AODtrack->TestFilterBit(1)))
continue;
2219 dEdx = AODtrack->GetDetPid()->GetTPCsignal();
2221 if((ptTrk < 0.2) || (ptTrk > 5.0) || (TMath::Abs(etaTrk) > 0.8) || (nClustTPC < 70) || (dEdx < 10.0) || (Chi2Trk < 0.1))
continue;
2223 if(AODtrack->GetDetPid() && Chi2Trk > 0.2) multTPC++;
2225 if(!AODtrack->TestFilterBit(16) || AODtrack->Chi2perNDF() < 0.1)
continue;
2227 Double_t bCov[3] = {-99., -99., -99.};
2228 AliAODTrack copy(*AODtrack);
2229 if(copy.PropagateToDCA(
fVevent->GetPrimaryVertex(),
fVevent->GetMagneticField(), 100., b, bCov) && TMath::Abs(b[0]) < 0.3 && TMath::Abs(b[1]) < 0.3) multGlobal++;
2238 Int_t multEsd = ((AliAODHeader*)
fAOD->GetHeader())->GetNumberOfESDTracks();
2247 else if(bIsPileup==kFALSE) {
2293 Bool_t bIsOutLier=kFALSE;
2295 if(multTPC < (-20.0+1.15*multGlobal) || multTPC > (200.+1.45*multGlobal)) { bIsOutLier = kTRUE;}
2367 if(multEtaNeg<2 || multEtaPos<2)
return;
2393 Int_t icentV0Qn = centrCL1;
2396 if(icentV0Qn>90)
return;
2406 const AliAODVZERO *fAODV0 =
fAOD->GetVZEROData();
2425 for(
int iV0 = 0; iV0 < 64; iV0++) {
2427 fMultv0 = fAODV0->GetMultiplicity(iV0);
2430 fMultv0 = fMultv0 *
fHCorrectV0M->GetBinContent(iV0+1);
2437 phiV0 = TMath::PiOver4()*(0.5 + iV0 % 8);
2440 Qxcn2 += TMath::Cos(2*phiV0) * fMultv0;
2441 Qycn2 += TMath::Sin(2*phiV0) * fMultv0;
2442 Qxcn3 += TMath::Cos(3*phiV0) * fMultv0;
2443 Qycn3 += TMath::Sin(3*phiV0) * fMultv0;
2447 Qxan2 += TMath::Cos(2*phiV0) * fMultv0;
2448 Qyan2 += TMath::Sin(2*phiV0) * fMultv0;
2449 Qxan3 += TMath::Cos(3*phiV0) * fMultv0;
2450 Qyan3 += TMath::Sin(3*phiV0) * fMultv0;
2457 QxanCor = Qxan3/sumMa;
2458 QyanCor = Qyan3/sumMa;
2459 QxcnCor = Qxcn3/sumMc;
2460 QycnCor = Qycn3/sumMc;
2476 QxanCor = Qxan2/sumMa;
2477 QyanCor = Qyan2/sumMa;
2478 QxcnCor = Qxcn2/sumMc;
2479 QycnCor = Qycn2/sumMc;
2509 PsiNV0C = 1.0/
gPsiN*( TMath::ATan2(QycnCor,QxcnCor) + TMath::Pi() );
2510 PsiNV0A = 1.0/
gPsiN*( TMath::ATan2(QyanCor,QxanCor) + TMath::Pi() );
2513 PsiNV0C = 1.0/
gPsiN*TMath::ATan2(QycnCor,QxcnCor) ;
2514 if(PsiNV0C<0.) PsiNV0C += 2*TMath::Pi()/
gPsiN;
2516 PsiNV0A = 1.0/
gPsiN*TMath::ATan2(QyanCor,QxanCor) ;
2517 if(PsiNV0A<0.) PsiNV0A += 2*TMath::Pi()/
gPsiN;
2562 PDGmassProton *= PDGmassProton;
2563 PDGmassPion *= PDGmassPion;
2564 PDGmassKaon *= PDGmassKaon;
2566 Double_t dEdx1, dEdx2, dPhi1, dPhi2, dPt1, dPt2;
2567 Double_t dEta1, dEta2, ptw1, ptw2, deltaPhi;
2571 Double_t nSigTOFproton,nSigTPCproton;
2572 Double_t nSigTOFpion2, nSigTPCpion2;
2573 Double_t nSigTOFkaon2, nSigTPCkaon2;
2574 Double_t nSigTOFproton2,nSigTPCproton2;
2583 Int_t ptBin,gCharge1,gCharge2;
2595 Int_t skipPairHBT = 0;
2601 Double_t Chi2Trk1,ptwPion1,ptwKaon1,ptwProton1;
2602 Double_t Chi2Trk2,ptwPion2,ptwKaon2,ptwProton2;
2603 Double_t wNUAPion1,wNUAKaon1,wNUAProton1;
2604 Double_t wNUAPion2,wNUAKaon2,wNUAProton2;
2613 Bool_t isProton1 = kFALSE;
2616 Bool_t isProton2 = kFALSE;
2618 Int_t multPOI1st = 0;
2619 Int_t multPOI2nd = 0;
2702 for(
Int_t itrack = 0; itrack < ntracks; itrack++) {
2704 AliAODTrack *track=
dynamic_cast<AliAODTrack*
>(
fVevent->GetTrack(itrack));
2705 if(!track)
continue;
2707 if(!track->TestFilterBit(
fFilterBit))
continue;
2713 dPhi1 = track->Phi();
2714 dEta1 = track->Eta();
2715 gCharge1 = track->Charge();
2716 dEdx1 = track->GetDetPid()->GetTPCsignal();
2717 Chi2Trk1 = track->Chi2perNDF();
2779 nSigTPCpion = nSigPionTPC[itrack];
2780 nSigTPCkaon = nSigKaonTPC[itrack];
2781 nSigTPCproton = nSigProtonTPC[itrack];
2782 nSigTOFpion = nSigPionTOF[itrack];
2783 nSigTOFkaon = nSigKaonTOF[itrack];
2784 nSigTOFproton = nSigProtonTOF[itrack];
2823 if(dPt1<=0.6 && TMath::Abs(nSigTPCpion)<=3.0){
2826 else if(dPt1>0.6 && dPt1<=2.0 && TMath::Abs(nSigTPCpion)<=3.0 && TMath::Abs(nSigTOFpion)<=2.0 ){
2830 if(dPt1<=0.45 && TMath::Abs(nSigTPCkaon)<=3.0){
2833 else if(dPt1>0.45 && dPt1<=2.0 && TMath::Abs(nSigTPCkaon)<=3.0 && TMath::Abs(nSigTOFkaon)<=2.0){
2837 if(dPt1<=0.8 && TMath::Abs(nSigTPCproton)<=3.0){
2839 if(dPt1<0.4) isProton1 = kFALSE;
2841 else if(dPt1>0.8 && dPt1<=2.0 && TMath::Abs(nSigTPCproton)<=3.0 && TMath::Abs(nSigTOFproton)<=3.0){
2962 if(w1NUA==0) w1NUA = 1.0;
2967 sumQxTPCneg = sumTPCQn2x[0];
2968 sumQyTPCneg = sumTPCQn2y[0];
2969 sumQxTPCpos = sumTPCQn2x[1];
2970 sumQyTPCpos = sumTPCQn2y[1];
2971 SumWgtNeg = SumWEtaNeg;
2972 SumWgtPos = SumWEtaPos;
2977 sumQxTPCneg -= w1NUA*TMath::Cos(
gPsiN*dPhi1);
2978 sumQyTPCneg -= w1NUA*TMath::Sin(
gPsiN*dPhi1);
2981 else if(dEta1 > 0.05){
2982 sumQxTPCpos -= w1NUA*TMath::Cos(
gPsiN*dPhi1);
2983 sumQyTPCpos -= w1NUA*TMath::Sin(
gPsiN*dPhi1);
3002 for(
Int_t jtrack = 0; jtrack < ntracks; jtrack++) {
3004 if(jtrack==itrack)
continue;
3008 AliAODTrack *track2=
dynamic_cast<AliAODTrack*
>(
fVevent->GetTrack(jtrack));
3009 if(!track2)
continue;
3011 if(!(track2->TestFilterBit(
fFilterBit)))
continue;
3013 dPt2 = track2->Pt();
3014 dPhi2 = track2->Phi();
3015 dEta2 = track2->Eta();
3016 gCharge2= track2->Charge();
3017 dEdx2 = track2->GetDetPid()->GetTPCsignal();
3018 Chi2Trk2 = track2->Chi2perNDF();
3049 nSigTPCpion2 = nSigPionTPC[jtrack];
3050 nSigTPCkaon2 = nSigKaonTPC[jtrack];
3051 nSigTPCproton2 = nSigProtonTPC[jtrack];
3052 nSigTOFpion2 = nSigPionTOF[jtrack];
3053 nSigTOFkaon2 = nSigKaonTOF[jtrack];
3054 nSigTOFproton2 = nSigProtonTOF[jtrack];
3067 if(dPt2<=0.6 && TMath::Abs(nSigTPCpion2)<=3.0){
3070 else if(dPt2>0.6 && dPt2<=2.0 && TMath::Abs(nSigTPCpion2)<=3.0 && TMath::Abs(nSigTOFpion2)<=2.0 ){
3074 if(dPt2<=0.45 && TMath::Abs(nSigTPCkaon2)<=3.0){
3077 else if(dPt2>0.45 && dPt2<=2.0 && TMath::Abs(nSigTPCkaon2)<=3.0 && TMath::Abs(nSigTOFkaon2)<=2.0){
3081 if(dPt2<=0.8 && TMath::Abs(nSigTPCproton2)<=3.0){
3083 if(dPt2<0.4) isProton2 = kFALSE;
3085 else if(dPt2>0.8 && dPt2<=2.0 && TMath::Abs(nSigTPCproton2)<=3.0 && TMath::Abs(nSigTOFproton2)<=3.0){
3156 WgtEPPion = ptwPion1*ptwPion2*wNUAPion1*wNUAPion2;
3182 WgtEPKaon = ptwKaon1*ptwKaon2*wNUAKaon1*wNUAKaon2;
3208 WgtEPProton = ptwProton1*ptwProton2*wNUAProton1*wNUAProton2;
3242 if(w2NUA==0) w2NUA = 1.0;
3256 sumQxTPCneg2 = sumQxTPCneg;
3257 sumQyTPCneg2 = sumQyTPCneg;
3258 sumQxTPCpos2 = sumQxTPCpos;
3259 sumQyTPCpos2 = sumQyTPCpos;
3260 SumWgtNeg2 = SumWgtNeg;
3261 SumWgtPos2 = SumWgtPos;
3267 sumQxTPCneg2 -= w2NUA*TMath::Cos(
gPsiN*dPhi2);
3268 sumQyTPCneg2 -= w2NUA*TMath::Sin(
gPsiN*dPhi2);
3269 SumWgtNeg2 -= w2NUA;
3271 else if(dEta2 > 0.05){
3272 sumQxTPCpos2 -= w2NUA*TMath::Cos(
gPsiN*dPhi2);
3273 sumQyTPCpos2 -= w2NUA*TMath::Sin(
gPsiN*dPhi2);
3274 SumWgtPos2 -= w2NUA;
3277 if(SumWgtNeg2>0 && SumWgtPos2>0){
3278 sumQyTPCneg2 = sumQyTPCneg2/SumWgtNeg2;
3279 sumQxTPCneg2 = sumQxTPCneg2/SumWgtNeg2;
3281 sumQyTPCpos2 = sumQyTPCpos2/SumWgtPos2;
3282 sumQxTPCpos2 = sumQxTPCpos2/SumWgtPos2;
3287 PsiNTPCC = (1.0/
gPsiN)*( TMath::ATan2(sumQyTPCneg2,sumQxTPCneg2) ) ;
3288 if(PsiNTPCC<0.) PsiNTPCC += 2*TMath::Pi()/
gPsiN;
3290 PsiNTPCA = (1.0/
gPsiN)*( TMath::ATan2(sumQyTPCpos2,sumQxTPCpos2) ) ;
3291 if(PsiNTPCA<0.) PsiNTPCA += 2*TMath::Pi()/
gPsiN;
3302 WgtEP = ptw1*ptw2*w1NUA*w2NUA;
3306 deltaPhi = dPhi1 - dPhi2;
3310 if(gCharge1!=gCharge2) {
3311 fHist_Corr3p_EP_Norm_PN[QAindex][0]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0A), WgtEP*fWgtCent);
3312 fHist_Corr3p_EP_Norm_PN[QAindex][1]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0C), WgtEP*fWgtCent);
3313 fHist_Corr3p_EP_Norm_PN[QAindex][2]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCA),WgtEP*fWgtCent);
3314 fHist_Corr3p_EP_Norm_PN[QAindex][3]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCC),WgtEP*fWgtCent);
3332 if(isPion1 && isPion2){
3355 if(isKaon1 && isKaon2){
3378 if(isProton1 && isProton2){
3410 else if(gCharge1>0 && gCharge2>0 && skipPairHBT==0) {
3411 fHist_Corr3p_EP_Norm_PP[QAindex][0]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0A), WgtEP*fWgtCent);
3412 fHist_Corr3p_EP_Norm_PP[QAindex][1]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0C), WgtEP*fWgtCent);
3413 fHist_Corr3p_EP_Norm_PP[QAindex][2]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCA),WgtEP*fWgtCent);
3414 fHist_Corr3p_EP_Norm_PP[QAindex][3]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCC),WgtEP*fWgtCent);
3433 if(isPion1 && isPion2){
3456 if(isKaon1 && isKaon2){
3479 if(isProton1 && isProton2){
3511 else if(gCharge1<0 && gCharge2<0 && skipPairHBT==0){
3512 fHist_Corr3p_EP_Norm_NN[QAindex][0]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0A), WgtEP*fWgtCent);
3513 fHist_Corr3p_EP_Norm_NN[QAindex][1]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNV0C), WgtEP*fWgtCent);
3514 fHist_Corr3p_EP_Norm_NN[QAindex][2]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCA),WgtEP*fWgtCent);
3515 fHist_Corr3p_EP_Norm_NN[QAindex][3]->Fill(EvtCent, TMath::Cos(n*dPhi1 + m*dPhi2 - p*PsiNTPCC),WgtEP*fWgtCent);
3533 if(isPion1 && isPion2){
3556 if(isKaon1 && isKaon2){
3579 if(isProton1 && isProton2){
3617 if(TMath::Abs(pVtxZ) < 8.0 && gCharge1 > 0){
3681 else if(gCharge1<0){
3691 else if(gCharge1<0){
3701 else if(gCharge1<0){
3710 else if(gCharge1<0){
3854 sumTPCQn2x[0] = sumTPCQn2x[0]/SumWEtaNeg;
3855 sumTPCQn2y[0] = sumTPCQn2y[0]/SumWEtaNeg;
3856 sumTPCQn3x[0] = sumTPCQn3x[0]/SumWEtaNeg;
3857 sumTPCQn3y[0] = sumTPCQn3y[0]/SumWEtaNeg;
3858 sumTPCQn4x[0] = sumTPCQn4x[0]/SumWEtaNeg;
3859 sumTPCQn4y[0] = sumTPCQn4y[0]/SumWEtaNeg;
3861 sumTPCQn2x[1] = sumTPCQn2x[1]/SumWEtaPos;
3862 sumTPCQn2y[1] = sumTPCQn2y[1]/SumWEtaPos;
3863 sumTPCQn3x[1] = sumTPCQn3x[1]/SumWEtaPos;
3864 sumTPCQn3y[1] = sumTPCQn3y[1]/SumWEtaPos;
3865 sumTPCQn4x[1] = sumTPCQn4x[1]/SumWEtaPos;
3866 sumTPCQn4y[1] = sumTPCQn4y[1]/SumWEtaPos;
3868 sumTPCQn2x[3] = sumTPCQn2x[3]/SumWEtaFull;
3869 sumTPCQn2y[3] = sumTPCQn2y[3]/SumWEtaFull;
3870 sumTPCQn3x[3] = sumTPCQn3x[3]/SumWEtaFull;
3871 sumTPCQn3y[3] = sumTPCQn3y[3]/SumWEtaFull;
3872 sumTPCQn4x[3] = sumTPCQn4x[3]/SumWEtaFull;
3873 sumTPCQn4y[3] = sumTPCQn4y[3]/SumWEtaFull;
3878 PsiNTPCC = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[0],sumTPCQn2x[0]) + TMath::Pi() ) ;
3879 PsiNTPCA = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[1],sumTPCQn2x[1]) + TMath::Pi() ) ;
3880 PsiNTPCF = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[3],sumTPCQn2x[3]) + TMath::Pi() ) ;
3883 PsiNTPCC = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[0],sumTPCQn2x[0]) ) ;
3884 if(PsiNTPCC<0.) PsiNTPCC += 2*TMath::Pi()/
gPsiN;
3885 PsiNTPCA = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[1],sumTPCQn2x[1]) ) ;
3886 if(PsiNTPCA<0.) PsiNTPCA += 2*TMath::Pi()/
gPsiN;
3887 PsiNTPCF = (1.0/
gPsiN)*( TMath::ATan2(sumTPCQn2y[3],sumTPCQn2x[3]) ) ;
3888 if(PsiNTPCF<0.) PsiNTPCF += 2*TMath::Pi()/
gPsiN;
4011 AliDebug(2,
"\n\n ::GetWDist => One of vertices is not valid\n\n");
4014 static TMatrixDSym vVb(3);
4016 double dx = v0->GetX()-v1->GetX();
4017 double dy = v0->GetY()-v1->GetY();
4018 double dz = v0->GetZ()-v1->GetZ();
4019 double cov0[6],cov1[6];
4020 v0->GetCovarianceMatrix(cov0);
4021 v1->GetCovarianceMatrix(cov1);
4022 vVb(0,0) = cov0[0]+cov1[0];
4023 vVb(1,1) = cov0[2]+cov1[2];
4024 vVb(2,2) = cov0[5]+cov1[5];
4025 vVb(1,0) = vVb(0,1) = cov0[1]+cov1[1];
4026 vVb(0,2) = vVb(1,2) = vVb(2,0) = vVb(2,1) = 0.;
4028 if (!vVb.IsValid()) {
4029 AliDebug(2,
"Singular Matrix\n");
4032 dist = vVb(0,0)*dx*dx + vVb(1,1)*dy*dy + vVb(2,2)*dz*dz
4033 + 2*vVb(0,1)*dx*dy + 2*vVb(0,2)*dx*dz + 2*vVb(1,2)*dy*dz;
4034 return dist>0 ? TMath::Sqrt(dist) : -1;
4040 const int kMinPlpContrib = 5;
4041 const double kMaxPlpChi2 = 5.0;
4042 const double kMinWDist = 15;
4044 const AliVVertex* vtPrm = 0;
4045 const AliVVertex* vtPlp = 0;
4049 if(!(nPlp=faod->GetNumberOfPileupVerticesTracks()))
4052 vtPrm = faod->GetPrimaryVertex();
4053 if(vtPrm == faod->GetPrimaryVertexSPD())
4058 for(
int ipl=0;ipl<nPlp;ipl++) {
4059 vtPlp = (
const AliVVertex*)faod->GetPileupVertexTracks(ipl);
4060 if (vtPlp->GetNContributors() < kMinPlpContrib)
continue;
4061 if (vtPlp->GetChi2perNDF() > kMaxPlpChi2)
continue;
4066 double wDst =
GetWDist(vtPrm,vtPlp);
4067 if (wDst<kMinWDist)
continue;
4086 cout<<
"\n =========> Info: Using MC efficiency correction Map <=========== "<<endl;
4087 for(
int i=0;i<10;i++) {
4091 for(
int i=0;i<10;i++) {
4114 std::cout<<
"\n\n !!!!**** Warning: FB Efficiency File/List not found Use Weight = 1.0 *****\n\n"<<std::endl;
4129 else if(fCent>=5.0 && fCent<10){
4132 else if(fCent>=10.0) {
4133 cIndex = abs(fCent/10.0)+1;
4140 Float_t fMeanTPC[100] = {2902.95,2758.33,2642.78,2536.67,2435.37,2340.06,2248.44,2163.71,2080.49,2001.54,1925.86,1852.64,1781.97,1715.56,1650.53,1587.23,1527.51,1468.19,1412.73,1357.86,1305.35,1254.33,1205.57,1157.28,1111.53,1066.42,1023.15,981.594,940.795,901.766,863.651,826.183,790.53,756.358,722.654,690.513,659.443,628.807,599.748,571.664,544.446,518.042,492.369,468.072,444.694,422.487,400.104,379.129,359.147,339.62,320.817,302.788,285.791,269.015,253.688,238.671,224.039,209.932,196.915,184.647,172.76,161.381,150.395,140.288,131.033,121.58,113.112,104.938,97.3078,90.2178,83.5974,77.2645,70.7126,65.4424,60.1404,55.5644,50.8314,46.3761,43.024,38.625,35.3435,32.2304,29.4192,26.821,24.3303,21.9332,19.4215,16.7163,14.9414,13.1092,0.};
4142 Float_t fSigmaTPC[100] = {122.209,107.901,103.452,100.498,97.7403,94.7845,93.2543,90.0548,88.1106,85.7382,84.0812,82.2978,80.3817,78.6002,77.3448,75.5086,73.6842,71.9733,70.3447,69.1999,67.878,66.3511,65.0406,63.4866,62.4409,60.7899,59.1328,58.426,56.8618,55.8871,54.1031,53.4959,52.0482,51.0441,49.6218,48.7646,47.5166,46.5247,45.0727,44.4311,43.4531,42.0404,41.0238,40.1384,39.2588,38.2461,36.5951,36.0552,35.3727,33.7883,32.7167,32.4486,31.3709,30.3444,29.505,28.5139,27.4471,26.5359,25.9506,25.127,24.3797,23.2985,22.279,21.4698,20.781,20.8193,19.9509,18.8036,17.9145,16.961,16.7375,15.852,14.9324,14.7663,13.5969,13.4533,12.3067,12.7835,11.7283,10.6758,10.6676,10.6492,9.04614,8.89065,8.66093,8.50997,7.98812,6.91087,7.12045,7.29593,0.};
4144 for(
int i=0;i<90;i++) {
4198 fHistPileUpCount =
new TH1F(
"fHistPileUpCount",
"fHistPileUpCount", 15, 0., 15.);
4207 fHistPileUpCount->GetXaxis()->SetBinLabel(8,Form(
"multESDTPCDif>%d",puConst));
4208 fHistPileUpCount->GetXaxis()->SetBinLabel(9,Form(
"multGlobTPCDif>%d",puConst));
4213 fHistMultSelPUCount =
new TH1F(
"fHistMultSelPileUpCount",
"no of PU Event from MultSelTask", 10, 0., 10);
4264 cout<<
"\n =========== Info:: Using Centrality wgts for run = "<<run<<
"============"<<endl;
4267 cout<<
"\n =========== Info:: No Centrality wgt. correction.!! for run = "<<run<<
"============"<<endl;
4271 cout<<
"\n =========== Info:: Setting up V0 gain correction for run = "<<run<<
"============"<<endl;
4275 cout<<
"\n =========== Info:: No V0 gain correction..!!! for run = "<<run<<
"============"<<endl;
4278 cout<<
"\n =========== Info:: Setting up V0 <Qn> recentering for run = "<<run<<
"============"<<endl;
4282 cout<<
"\n =========== Info:: No V0 <Qn> recentering..!!! for run = "<<run<<
"============"<<endl;
4286 cout<<
"\n ======== !!!! Error:: List For V0-Gain and Qn not found for run "<<run<<
"============"<<endl;
4301 for(
int i=0;i<5;i++){
4306 cout<<
"\n=========== Info:: Setting up NUA corrections for run "<<run<<
"============"<<endl;
4311 printf(
"\n ******** Warning: No NUA Correction for Charge Particle in run %d, Use Wgt = 1.0 ********* \n",run);
4325 for(
int i=0;i<5;i++){
4339 cout<<
"\n=========== Info:: Setting up --> PID NUA corrections for run = "<<run<<
"============"<<endl;
4342 cout<<
"\n=========== WARNING :: PID NUA corrections NOT found for run = "<<run<<
"============"<<endl;
4346 printf(
"\n ******** Error/Warning: No NUA Correction found for PID for run %d, Use PID Wgt = 1.0 ********* \n",run);
TProfile * fHist_Corr3p_pTSum_EP_V0C_PP[2][6]
TProfile * fTPCFQ2xVsCentRun
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_EP_Refm_NN[2][4]
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0C_PP[2][6]
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_Proton_EP_Refm_PN[2][4]
Double_t fCentralityPercentMin
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0A_PP[2][6]
TProfile * fTPCCQ2yVsCentRun
TProfile * fHist_Corr3p_Kaon_EP_Refm_NN[2][4]
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0C_PP[2][6]
Int_t GetCentralityScaled0to10(Double_t fCent)
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PN[2][6]
void GetV0MCorrectionHist(Int_t run)
TProfile * fV0AQ2xVsCentRun
TH2D * fHAvgerageQnV0A
To read Gain Correction file.
TProfile * fHist_Corr3p_pTDiff_EP_V0C_PP[2][6]
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_Kaon_EP_Norm_PP[2][4]
AliAnalysisUtils * fAnalysisUtil
Double_t fPileUpConstParm
TProfile * fHist_Corr3p_pTSum_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_pTSum_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_Proton_EP_Norm_NN[2][4]
TH2F * fHistTPConlyVsCL1Before
Double_t fPileUpSlopeParm
TProfile * fV0AQ2yVsCentRun
TH3F * fHist3DEtaPhiVz_Pos_Run[4][5]
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0A_PP[2][6]
AliMultSelection * fMultSelection
PID response Handler.
TH2F * fHistRawVsCorrMultFB
TProfile * fV0CQ3yVsCentRun
TProfile * fHist_Corr3p_Kaon_EP_Norm_PN[2][4]
TProfile * fTPCCQ3xVsCentRun
TProfile * fHist_Corr3p_EtaDiff_EP_V0A_PN[2][6]
TH1F * fHistPtwithTPCNsigma[3]
last in the list
TProfile * fHist_Reso2n_EP_Norm_Det[2][4]
TProfile * fV0CQ2yVsCentRun
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0C_PP[2][6]
TH1D * fFB_Efficiency_Proton_Neg_Cent[10]
TH1D * fFB_Efficiency_Pion_Cent[10]
for charge
TH2F * fQAEtaPhiAfterNUAPion
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0A_PP[2][6]
TH2F * fHistTPCVsESDTrkBefore
TProfile * fV0CQ3xVsCentRun
virtual void UserCreateOutputObjects()
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0C_NN[2][6]
TH2F * fHTPCAEventPlaneVsCent
TH2F * fHistEtaPtAfter
Eta-Pt acceptance.
TProfile * fTPCCQ4yVsCentRun
TProfile * fHist_Corr3p_EP_Norm_NN[2][4]
TProfile * fHist_Corr3p_pTDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_EtaDiff_EP_V0C_PP[2][6]
TProfile * fTPCCQ2xVsCentRun
TH3D * fHCorrectNUAnegProton[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0C_PP[2][6]
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0A_PN[2][6]
TH3D * fHCorrectNUAposPion[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0C_PP[2][6]
TH2F * fHistGlobalVsV0MBefore
TH1F * fCentDistAfter
without PileUp cut
TProfile * fHist_Corr3p_Pion_EP_Norm_PP[2][4]
void GetNUACorrectionHist(Int_t run)
TH3D * fHCorrectNUAnegPion[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_EtaDiff_EP_V0A_NN[2][6]
AliPIDResponse * fPIDResponse
aod
TProfile * fHist_Corr3p_pTSum_EP_V0C_PN[2][6]
TH1D * fHCentWeightForRun
V0C Average <Qn>, n=2,3.
TProfile * fHist_Corr3p_EtaDiff_EP_V0A_PP[2][6]
TH2F * fHV0CEventPlaneVsCent
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0C_PP[2][6]
TH2F * fHistTPConlyVsCL1After
TProfile * fTPCAQ3yVsCentRun
TProfile * fTPCCQ4xVsCentRun
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0C_NN[2][6]
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0A_PP[2][6]
TProfile * fTPCAQ3xVsCentRun
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0C_NN[2][6]
TH3F * fHist3DEtaPhiVz_Neg_Run[4][5]
4 particle 5 centrality bin
TProfile * fHist_Corr3p_Proton_EP_Norm_PP[2][4]
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0A_PP[2][6]
TProfile * fTPCAQ2xVsCentRun
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0A_NN[2][6]
virtual ~AliAnalysisTaskCMEV0PID()
TProfile * fHist_Corr3p_Kaon_EP_Refm_PN[2][4]
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PP[2][6]
TH2F * fHistTPCVsESDTrkAfter
TH2F * fHistTPCnSigmavsPtAfter[3]
Double_t fCentralityPercentMax
TH2F * fHistTOFMatchCount
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0C_PP[2][6]
TH1D * fHCorrectV0M
with PileUp cut
TProfile * fTPCAQ4yVsCentRun
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0A_NN[2][6]
TH2F * fHistTPCdEdxvsPtPIDAfter[3]
TProfile * fV0CQ2xVsCentRun
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0C_NN[2][6]
TH2F * fHistTPCvsGlobalMultBefore
Eta-Pt acceptance.
TH3D * fHCorrectNUAposKaon[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_Pion_EP_Norm_PN[2][4]
TProfile * fTPCAQ4xVsCentRun
TH2F * fQAEtaPhiAfterNUAProton
TProfile * fHist_Corr3p_Kaon_EP_Norm_NN[2][4]
TProfile * fHist_Corr3p_Pion_EP_Norm_NN[2][4]
TProfile * fHist_Corr3p_Pion_EP_Refm_PN[2][4]
TH2F * fHTPCCEventPlaneVsCent
void SetupMCcorrectionMap()
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0A_PN[2][6]
TH2F * hCentvsTPCmultCuts
TProfile * fHist_Corr3p_pTDiff_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_Proton_EP_Norm_PN[2][4]
TProfile * fTPCAQ2yVsCentRun
TH2F * fHistTOFnSigmavsPtAfter[3]
TH1F * fCentDistBefore
To fill VOM multiplicity.
TList * fListHist
Event selection.
virtual void UserExec(Option_t *)
TH1F * fHistPileUpCount
Task input parameters FB / cut values etc.
Bool_t PileUpMultiVertex(const AliAODEvent *faod)
TH2F * fHTPCEventPlaneVsCent
TH2F * fHistGlobalVsV0MAfter
TH3D * fHCorrectNUAnegKaon[5]
5 centrality bin, read NUA from file
TH3D * fHCorrectNUAneg[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_Kaon_EP_Refm_PP[2][4]
double GetWDist(const AliVVertex *v0, const AliVVertex *v1)
TH3D * fHCorrectNUApos[5]
TH1D * fFB_Efficiency_Cent[10]
4 particle 5 centrality bin
TH3F * fHistTPCTOFnSigmavsPtAfter[3]
TProfile * fHist_Corr3p_EP_Norm_PN[2][4]
5 centrality bin, read NUA from file
TProfile * fV0AQ3xVsCentRun
TH1D * fFB_Efficiency_Kaon_Cent[10]
AliAnalysisTaskCMEV0PID()
TProfile * fHist_Corr3p_EtaDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_pTSum_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0C_NN[2][6]
void SetupEventAndTaskConfigInfo()
TProfile * fHist_Corr3p_pTDiff_EP_V0C_NN[2][6]
TH2F * fHistTPCvsGlobalMultAfter
TH1F * fHistPtwithTOFmasscut[3]
void SetUpCentralityOutlierCut()
TH2F * fQAEtaPhiAfterNUAKaon
TProfile * fTPCCQ3yVsCentRun
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Pion_EP_Refm_PP[2][4]
TProfile * fHist_Corr3p_Pion_pTSum_EP_V0A_NN[2][6]
TH2F * fQAEtaPhiAfterNUA
Event weights for non-flat centrality.
TProfile * fHist_Corr3p_Kaon_pTDiff_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0A_PN[2][6]
TH3D * fHCorrectNUAposProton[5]
5 centrality bin, read NUA from file
TProfile * fHist_Corr3p_EP_Refm_PN[2][4]
TH1F * fHistMultSelPUCount
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Pion_pTDiff_EP_V0C_PN[2][6]
TH1D * fFB_Efficiency_Proton_Pos_Cent[10]
TProfile * fHist_Corr3p_pTDiff_EP_V0A_PN[2][6]
TProfile * fHist_Corr3p_pTDiff_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Pion_EP_Refm_NN[2][4]
TH2F * fV0MultChVsRun
Full Event plane.
TProfile * fHist_Corr3p_Proton_pTSum_EP_V0C_PP[2][6]
TProfile * fHist_Corr3p_Kaon_EtaDiff_EP_V0C_PN[2][6]
TH1F * fHistPtwithTOFSignal[3]
TProfile * fV0AQ3yVsCentRun
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Proton_EP_Refm_PP[2][4]
TH2F * fHistTPConlyVsV0MAfter
TProfile * fHist_Corr3p_EP_Refm_PP[2][4]
TProfile * fTPCFQ2yVsCentRun
TProfile * fHist_Reso2n_EP_Refm_Det[2][4]
TProfile * fHist_Corr3p_Kaon_pTSum_EP_V0A_NN[2][6]
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Pion_EtaDiff_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_pTSum_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_Proton_EP_Refm_NN[2][4]
TProfile * fHist_Corr3p_Proton_EtaDiff_EP_V0C_PN[2][6]
TProfile * fHist_Corr3p_Proton_pTDiff_EP_V0A_PP[2][6]
TProfile * fHist_Corr3p_EtaDiff_EP_V0C_NN[2][6]
TH2F * fHV0AEventPlaneVsCent
TH1F * fHistTaskConfigParameters
TH2D * fHAvgerageQnV0C
V0A Average <Qn>, n=2,3.
TProfile * fHist_Corr3p_EP_Norm_PP[2][4]
TH2F * fHistTPConlyVsV0MBefore