]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG2/FLOW/macros/CompareFlowResults.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Additional protection (macosx with gfortran)
[u/mrichter/AliRoot.git] / PWG2 / FLOW / macros / CompareFlowResults.C
CommitLineData
8c5ec46a 1//type of analysis can be: ESD, AOD, MC, ESDMC0, ESDMC1
2//const TString type = "ESD";
694e0378 3void CompareFlowResults(TString type="ESD")
9d3cbee7 4{
8c5ec46a 5 gSystem->AddIncludePath("-I$ALICE_ROOT/include");
6 gSystem->AddIncludePath("-I$ROOTSYS/include");
7
8 //load needed libraries:
9 gSystem->Load("libTree.so");
10 gSystem->Load("libANALYSIS.so");
11 gSystem->Load("libPWG2flow.so");
12 cerr<<"libPWG2flow.so loaded ..."<<endl;
13 cout<<endl;
14
15
16
17
18 //==================================================================================
19 // set here which plots will be shown by default
20 //==================================================================================
21 //RP = particles used to determine the reaction plane
694e0378 22 Bool_t plotIntFlowRP = kFALSE; //integrated flow RP
8c5ec46a 23 Bool_t plotDiffFlowPtRP = kFALSE; //differential flow (Pt,RP)
24 Bool_t plotDiffFlowEtaRP = kFALSE; //differential flow (Eta,RP)
25 //POI = particle of interest
694e0378 26 Bool_t plotIntFlowPOI = kTRUE; //integrated flow POI
27 Bool_t plotDiffFlowPtPOI = kTRUE; //differential flow (Pt,POI)
28 Bool_t plotDiffFlowEtaPOI = kTRUE; //differential flow (Eta,POI)
8c5ec46a 29 //==================================================================================
30
31
32
33
34 //==================================================================================
35 // accessing output files
36 //==================================================================================
37
38 //open the output files:
39 TString inputFileNameMCEP = "outputMCEPanalysis";
40 TFile* fileMCEP = NULL;
41 fileMCEP = TFile::Open(((inputFileNameMCEP.Append(type)).Append(".root")).Data(), "READ");
42
43 TString inputFileNameLYZ1 = "outputLYZ1analysis";
44 TFile* fileLYZ1 = NULL;
45 fileLYZ1 = TFile::Open(((inputFileNameLYZ1.Append(type)).Append(".root")).Data(), "READ");
46
47 TString inputFileNameLYZ2 = "outputLYZ2analysis";
48 TFile* fileLYZ2 = NULL;
49 fileLYZ2 = TFile::Open(((inputFileNameLYZ2.Append(type)).Append(".root")).Data(), "READ");
50
51 TString inputFileNameLYZEP = "outputLYZEPanalysis";
52 TFile* fileLYZEP = NULL;
53 fileLYZEP = TFile::Open(((inputFileNameLYZEP.Append(type)).Append(".root")).Data(), "READ");
54
55 TString inputFileNameFQD = "outputFQDanalysis";
56 TFile* fileFQD = NULL;
57 fileFQD = TFile::Open(((inputFileNameFQD.Append(type)).Append(".root")).Data(), "READ");
58
59 TString inputFileNameGFC = "outputGFCanalysis";
60 TFile* fileGFC = NULL;
61 fileGFC = TFile::Open(((inputFileNameGFC.Append(type)).Append(".root")).Data(), "READ");
62
63 TString inputFileNameQC = "outputQCanalysis";
64 TFile* fileQC = NULL;
65 fileQC = TFile::Open(((inputFileNameQC.Append(type)).Append(".root")).Data(), "READ");
66 //==================================================================================
1fca9c90 67
a88e5ac2 68
a88e5ac2 69
9d3cbee7 70
8c5ec46a 71 //==================================================================================
72 // cosmetics
73 //==================================================================================
74 //removing the title and stat. box from all histograms:
75 gStyle->SetOptTitle(0);
76 gStyle->SetOptStat(0);
77
78 //choosing the style and color of mesh for MC error bands:
79 Int_t meshStyle = 1001;
80 Int_t meshColor = kRed-10;
81
82 //marker style and color (int. flow)
83 Int_t markerStyle = 21;
84 Int_t markerColor = kRed-3;
85 //==================================================================================
86
87
88
89
90 //==================================================================================
91 // INTEGRATED FLOW
92 //==================================================================================
93 //the number of different methods:
94 const Int_t nMethods=12;
95
96 //booking the histogram for the integrated flow results from all methods:
97 TH1D* intFlowAll = new TH1D("intFlowAll","Integrated Flow",nMethods,0,nMethods);
98 //intFlowAll->SetLabelSize(0.036,"X");
99 //intFlowAll->SetLabelSize(0.036,"Y");
100 intFlowAll->SetMarkerStyle(markerStyle);
101 intFlowAll->SetMarkerColor(markerColor);
102 (intFlowAll->GetXaxis())->SetBinLabel(1,"v_{2}{MC}");
103 (intFlowAll->GetXaxis())->SetBinLabel(2,"v_{2}{2,GFC}");
104 (intFlowAll->GetXaxis())->SetBinLabel(3,"v_{2}{2,QC}");
105 (intFlowAll->GetXaxis())->SetBinLabel(4,"v_{2}{4,GFC}");
106 (intFlowAll->GetXaxis())->SetBinLabel(5,"v_{2}{4,QC}");
107 (intFlowAll->GetXaxis())->SetBinLabel(6,"v_{2}{6,GFC}");
108 (intFlowAll->GetXaxis())->SetBinLabel(7,"v_{2}{6,QC}");
109 (intFlowAll->GetXaxis())->SetBinLabel(8,"v_{2}{8,GFC}");
110 (intFlowAll->GetXaxis())->SetBinLabel(9,"v_{2}{8,QC}");
111 (intFlowAll->GetXaxis())->SetBinLabel(10,"v_{2}{FQD}");
112 (intFlowAll->GetXaxis())->SetBinLabel(11,"v_{2}{LYZ}");
113 (intFlowAll->GetXaxis())->SetBinLabel(12,"v_{2}{LYZEP}");
114
115 //booking the graph to store flow values and errors from all methods:
116 Double_t x[nMethods] = {0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, 9.5, 10.5, 11.5};
117 Double_t xError[nMethods] = {0.};
118 Double_t flowValue[nMethods] = {0.};//to be removed
119 Double_t flowError[nMethods] = {0.};//to be removed
120 Double_t flowValueRP[nMethods] = {0.};
121 Double_t flowErrorRP[nMethods] = {0.};
122 Double_t flowValuePOI[nMethods] = {0.};
123 Double_t flowErrorPOI[nMethods] = {0.};
124
125 //accessing the results for each method:
126 //MCEP = Monte Carlo Event Plane
127 TList *pListMCEP = NULL;
128 AliFlowCommonHist *mcepCommonHist = NULL;
129 AliFlowCommonHistResults *mcepCommonHistRes = NULL;
130 if(fileMCEP) {
131 fileMCEP->GetObject("cobjMCEP",pListMCEP);
132 if(pListMCEP) {
133 mcepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListMCEP->FindObject("AliFlowCommonHistMCEP"));
134 mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListMCEP->FindObject("AliFlowCommonHistResultsMCEP"));
135 if(mcepCommonHistRes) {
136 flowValue[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
137 flowError[0] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
138 flowValueRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
139 flowErrorRP[0] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
140 flowValuePOI[0] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
141 flowErrorPOI[0] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinError(1);
142 }
143 }
a88e5ac2 144 }
8c5ec46a 145
146 //LYZ1 = Lee-Yang Zeros (1st run)
147 TList *pListLYZ1 = NULL;
148 AliFlowCommonHist *lyz1CommonHist = NULL;
149 AliFlowCommonHistResults *lyz1CommonHistRes = NULL;
150 if(fileLYZ1) {
151 fileLYZ1->GetObject("cobjLYZ1",pListLYZ1);
152 if(pListLYZ1) {
153 lyz1CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ1->FindObject("AliFlowCommonHistLYZ1"));
154 lyz1CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ1->FindObject("AliFlowCommonHistResultsLYZ1"));
155 if(lyz1CommonHistRes) {
156 flowValue[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
157 flowError[10] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
158 flowValueRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
159 flowErrorRP[10] = (lyz1CommonHistRes->GetHistIntFlowRP())->GetBinError(1);
160 flowValuePOI[10] = (lyz1CommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
161 flowErrorPOI[10] = (lyz1CommonHistRes->GetHistIntFlowPOI())->GetBinError(1);
162 }
163 }
a88e5ac2 164 }
8c5ec46a 165
166 //LYZ2 = Lee-Yang Zeros (2nd run) (needed only for differential flow)
167 TList *pListLYZ2 = NULL;
168 AliFlowCommonHist *lyz2CommonHist = NULL;
169 AliFlowCommonHistResults *lyz2CommonHistRes = NULL;
170 if(fileLYZ2) {
171 fileLYZ2->GetObject("cobjLYZ2",pListLYZ2);
172 if(pListLYZ2) {
173 lyz2CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ2->FindObject("AliFlowCommonHistLYZ2"));
174 lyz2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ2->FindObject("AliFlowCommonHistResultsLYZ2"));
175 }
9d3cbee7 176 }
8c5ec46a 177
178 //LYZEP = Lee-Yang Zeros Event Plane
179 TList *pListLYZEP = NULL;
180 AliFlowCommonHist *lyzepCommonHist = NULL;
181 AliFlowCommonHistResults *lyzepCommonHistRes = NULL;
182 if(fileLYZEP) {
183 fileLYZEP->GetObject("cobjLYZEP",pListLYZEP);
184 if(pListLYZEP) {
185 lyzepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZEP->FindObject("AliFlowCommonHistLYZEP"));
186 lyzepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZEP->FindObject("AliFlowCommonHistResultsLYZEP"));
187 if(lyzepCommonHistRes) {
188 flowValue[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
189 //flowError[11] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
190 flowValueRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
191 //flowErrorRP[11] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
192 flowValuePOI[11] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
193 //flowErrorPOI[11] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinError(1);
194 }
195 }
9d3cbee7 196 }
9d3cbee7 197
8c5ec46a 198 //FQD = Fitting q-distribution
199 TList *pListFQD = NULL;
200 AliFlowCommonHist *fqdCommonHist = NULL;
201 AliFlowCommonHistResults *fqdCommonHistRes = NULL;
202 if(fileFQD) {
203 fileFQD->GetObject("cobjFQD",pListFQD);
204 if(pListFQD) {
205 fqdCommonHist = dynamic_cast<AliFlowCommonHist*> (pListFQD->FindObject("AliFlowCommonHistFQD"));
206 fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListFQD->FindObject("AliFlowCommonHistResultsFQD"));
207 if(fqdCommonHistRes) {
208 flowValue[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);//to be removed
209 flowError[9] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);//to be removed
210 flowValueRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinContent(1);
211 flowErrorRP[9] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinError(1);
212 flowValuePOI[9] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1);
213 flowErrorPOI[9] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinError(1);
214 }
215 }
9d3cbee7 216 }
a88e5ac2 217
8c5ec46a 218 //GFC = Generating Function Cumulants
219 TList *pListGFC = NULL;
220 AliFlowCommonHist *gfcCommonHist = NULL;
221 AliFlowCommonHistResults *gfcCommonHistRes2 = NULL;
222 AliFlowCommonHistResults *gfcCommonHistRes4 = NULL;
223 AliFlowCommonHistResults *gfcCommonHistRes6 = NULL;
224 AliFlowCommonHistResults *gfcCommonHistRes8 = NULL;
225 if(fileGFC) {
226 fileGFC->GetObject("cobjGFC",pListGFC);
227 if(pListGFC) {
228 gfcCommonHist = dynamic_cast<AliFlowCommonHist*> (pListGFC->FindObject("AliFlowCommonHistGFC"));
229 gfcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC"));
230 if(gfcCommonHistRes2) {
231 flowValue[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
232 flowError[1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
233 flowValueRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
234 flowErrorRP[1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
235 flowValuePOI[1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
236 flowErrorPOI[1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinError(1);
237 }
238 gfcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults4thOrderGFC"));
239 if(gfcCommonHistRes4) {
240 flowValue[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
241 flowError[3] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
242 flowValueRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
243 flowErrorRP[3] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
244 flowValuePOI[3] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
245 flowErrorPOI[3] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinError(1);
246 }
247 gfcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults6thOrderGFC"));
248 if(gfcCommonHistRes6) {
249 flowValue[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
250 flowError[5] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
251 flowValueRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
252 flowErrorRP[5] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
253 flowValuePOI[5] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
254 flowErrorPOI[5] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinError(1);
255 }
256 gfcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults8thOrderGFC"));
257 if(gfcCommonHistRes8) {
258 flowValue[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
259 flowError[7] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
260 flowValueRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
261 flowErrorRP[7] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
262 flowValuePOI[7] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
263 flowErrorPOI[7] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinError(1);
264 }
265 }
9d3cbee7 266 }
1fca9c90 267
8c5ec46a 268 //QC = Q-cumulants
269 TList *pListQC = NULL;
270 AliFlowCommonHist *qcCommonHist2 = NULL;
271 AliFlowCommonHist *qcCommonHist4 = NULL;
272 AliFlowCommonHist *qcCommonHist6 = NULL;
273 AliFlowCommonHist *qcCommonHist8 = NULL;
274 AliFlowCommonHistResults *qcCommonHistRes2 = NULL;
275 AliFlowCommonHistResults *qcCommonHistRes4 = NULL;
276 AliFlowCommonHistResults *qcCommonHistRes6 = NULL;
277 AliFlowCommonHistResults *qcCommonHistRes8 = NULL;
278
279 if(fileQC) {
280 fileQC->GetObject("cobjQC",pListQC);
281 if(pListQC) {
282 qcCommonHist2 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist2ndOrderQC"));
283 qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
284 if(qcCommonHistRes2) {
285 flowValue[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1);//to be removed
286 //flowError[2] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1);//to be removed
287 flowValueRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1);
288 //flowErrorRP[2] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1);
289 flowValuePOI[2] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1);
290 //flowErrorPOI[2] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinError(1);
291 }
292 qcCommonHist4 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist4thOrderQC"));
293 qcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults4thOrderQC"));
294 if(qcCommonHistRes4) {
295 flowValue[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1);//to be removed
296 //flowError[4] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1);//to be removed
297 flowValueRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1);
298 //flowErrorRP[4] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1);
299 flowValuePOI[4] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1);
300 //flowErrorPOI[4] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinError(1);
301 }
302 qcCommonHist6 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist6thOrderQC"));
303 qcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults6thOrderQC"));
304 if(qcCommonHistRes6) {
305 flowValue[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1);//to be removed
306 //flowError[6] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1);//to be removed
307 flowValueRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1);
308 //flowErrorRP[6] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1);
309 flowValuePOI[6] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1);
310 //flowErrorPOI[6] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinError(1);
311 }
312 qcCommonHist8 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist8thOrderQC"));
313 qcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults8thOrderQC"));
314 if(qcCommonHistRes8) {
315 flowValue[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1);//to be removed
316 //flowError[8] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1);//to be removed
317 flowValueRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1);
318 //flowErrorRP[8] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1);
319 flowValuePOI[8] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1);
320 //flowErrorPOI[8] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinError(1);
321 }
322 }
323 }
324
325 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
326 // !!!! to be removed !!!!
327 Double_t dMax=flowValue[0]+flowError[0];
328 Double_t dMin=flowValue[0]-flowError[0];
329
330 for(Int_t i=1;i<nMethods;i++) {
331 if(!(flowValue[i]==0. && flowError[i]==0.)) {
332 if(dMax<flowValue[i]+flowError[i]) dMax=flowValue[i]+flowError[i];
333 if(dMin>flowValue[i]-flowError[i]) dMin=flowValue[i]-flowError[i];
334 }
335 }
336 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
a88e5ac2 337
1fca9c90 338 //RP:
8c5ec46a 339 Double_t dMaxRP=flowValueRP[0]+flowErrorRP[0];
340 Double_t dMinRP=flowValueRP[0]-flowErrorRP[0];
341
342 for(Int_t i=1;i<nMethods;i++) {
343 if(!(flowValueRP[i]==0. && flowErrorRP[i]==0.)) {
344 if(dMaxRP<flowValueRP[i]+flowErrorRP[i]) dMaxRP=flowValueRP[i]+flowErrorRP[i];
345 if(dMinRP>flowValueRP[i]-flowErrorRP[i]) dMinRP=flowValueRP[i]-flowErrorRP[i];
346 }
347 }
1fca9c90 348
349 //POI:
8c5ec46a 350 Double_t dMaxPOI=flowValuePOI[0]+flowErrorPOI[0];
351 Double_t dMinPOI=flowValuePOI[0]-flowErrorPOI[0];
a88e5ac2 352
8c5ec46a 353 for(Int_t i=1;i<nMethods;i++) {
354 if(!(flowValuePOI[i]==0. && flowErrorPOI[i]==0.)) {
355 if(dMaxPOI<flowValuePOI[i]+flowErrorPOI[i]) dMaxPOI=flowValuePOI[i]+flowErrorPOI[i];
356 if(dMinPOI>flowValuePOI[i]-flowErrorPOI[i]) dMinPOI=flowValuePOI[i]-flowErrorPOI[i];
357 }
358 }
359
360 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
361 // !!!! to be removed !!!!
362 TGraph* flowResults = new TGraphErrors(nMethods, x, flowValue, xError, flowError);
363
364 flowResults->SetMarkerStyle(markerStyle);
365 flowResults->SetMarkerColor(markerColor);
366 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
1fca9c90 367
8c5ec46a 368 //RP:
369 TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP);
370
371 flowResultsRP->SetMarkerStyle(markerStyle);
372 flowResultsRP->SetMarkerColor(markerColor);
373
374 //POI:
375 TGraph* flowResultsPOI = new TGraphErrors(nMethods, x, flowValuePOI, xError, flowErrorPOI);
376
377 flowResultsPOI->SetMarkerStyle(markerStyle);
378 flowResultsPOI->SetMarkerColor(markerColor);
379
380 //-----------------------------------------------------------------------------------
381
382 //----------------------------------------------------------------------------------
383 //cosmetics: mesh for MC error bands (integrated flow)
384 TGraph* pMesh = NULL;//to be removed
385 TGraph* pMeshRP = NULL;
386 TGraph* pMeshPOI = NULL;
387
388 if(intFlowAll && mcepCommonHistRes) {
389 //Double_t valueMC = intFlowAll->GetBinContent(1);
390 //Double_t errorMC = intFlowAll->GetBinError(1);
391 //Int_t nPts = intFlowAll->GetNbinsX();
392
393 Int_t nPts = nMethods;
394 Double_t valueMC = flowValue[0];//to be removed
395 Double_t errorMC = flowError[0];//to be removed
396 Double_t valueMCRP = flowValueRP[0];
397 Double_t errorMCRP = flowErrorRP[0];
398 Double_t valueMCPOI = flowValuePOI[0];
399 Double_t errorMCPOI = flowErrorPOI[0];
400
401 pMesh = new TGraph(nPts);//to be removed
402 pMeshRP = new TGraph(nPts);
403 pMeshPOI = new TGraph(nPts);
404
405 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
406 // !!!! to be removed !!!!
407 pMesh->SetPoint(1,0,valueMC+errorMC);
408 pMesh->SetPoint(2,nPts+1,valueMC+errorMC);
409 pMesh->SetPoint(3,nPts+1,valueMC-errorMC);
410 pMesh->SetPoint(4,0,valueMC-errorMC);
411 pMesh->SetPoint(5,0,valueMC+errorMC);
412
413 pMesh->SetFillStyle(meshStyle);
414 pMesh->SetFillColor(meshColor);
415 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
416
417 //RP:
418 pMeshRP->SetPoint(1,0,valueMCRP+errorMCRP);
419 pMeshRP->SetPoint(2,nPts+1,valueMCRP+errorMCRP);
420 pMeshRP->SetPoint(3,nPts+1,valueMCRP-errorMCRP);
421 pMeshRP->SetPoint(4,0,valueMCRP-errorMCRP);
422 pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP);
423
424 pMeshRP->SetFillStyle(meshStyle);
425 pMeshRP->SetFillColor(meshColor);
426
427 //POI:
428 pMeshPOI->SetPoint(1,0,valueMCPOI+errorMCPOI);
429 pMeshPOI->SetPoint(2,nPts+1,valueMCPOI+errorMCPOI);
430 pMeshPOI->SetPoint(3,nPts+1,valueMCPOI-errorMCPOI);
431 pMeshPOI->SetPoint(4,0,valueMCPOI-errorMCPOI);
432 pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI);
433
434 pMeshPOI->SetFillStyle(meshStyle);
435 pMeshPOI->SetFillColor(meshColor);
436 }
437 //----------------------------------------------------------------------------------
1fca9c90 438
1fca9c90 439
8c5ec46a 440 //----------------------------------------------------------------------------------
441 //cosmetics: text (integrated flow)
442 //default text:
443 TPaveText *textDefault = new TPaveText(0.05,0.77,0.95,0.90,"NDC");
444 textDefault->SetTextFont(72);
445 textDefault->SetTextSize(0.08);
446 //textDefault->SetLineColor(kFALSE);
447 //textDefault->SetShadowColor(kFALSE);
448
449 TString *entryDefaultAvM = new TString("Average Multiplicity");
450 TString *entryDefaultAnd = new TString("and");
451 TString *entryDefaultNumOfEvts = new TString("Number of Events:");
452
453 textDefault->AddText(entryDefaultAvM->Data());
454 textDefault->AddText(entryDefaultAnd->Data());
455 textDefault->AddText(entryDefaultNumOfEvts->Data());
456
457 //results:
458 TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC");//to be removed
459 textResults->SetTextFont(72);//to be removed
460 textResults->SetTextSize(0.06);//to be removed
461 //textResults->SetLineColor(kFALSE);
462 //textResults->SetShadowColor(kFALSE);
463
464 //results (RP):
465 TPaveText *textResultsRP = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
466 textResultsRP->SetTextFont(72);
467 textResultsRP->SetTextSize(0.06);
468
469 //results (POI):
470 TPaveText *textResultsPOI = new TPaveText(0.05,0.12,0.95,0.70,"NDC");
471 textResultsPOI->SetTextFont(72);
472 textResultsPOI->SetTextSize(0.06);
473
474 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
475 // !!!! to be removed !!!!!!
476 TString *entryMC = new TString("MC ...... ");
477 TString *entryGFC = new TString("GFC ..... ");
478 TString *entryQC2 = new TString("QC{2} ... ");
479 TString *entryQC4 = new TString("QC{4} ... ");
480 TString *entryQC6 = new TString("QC{6} ... ");
481 TString *entryQC8 = new TString("QC{8} ... ");
482 TString *entryFQD = new TString("FQD ..... ");
483 TString *entryLYZ1 = new TString("LYZ ..... ");
484 TString *entryLYZEP = new TString("LYZEP ... ");
485 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
486
487 //RP:
488 TString *entryMCRP = new TString("MC ...... ");
489 TString *entryGFCRP = new TString("GFC ..... ");
490 TString *entryQC2RP = new TString("QC{2} ... ");
491 TString *entryQC4RP = new TString("QC{4} ... ");
492 TString *entryQC6RP = new TString("QC{6} ... ");
493 TString *entryQC8RP = new TString("QC{8} ... ");
494 TString *entryFQDRP = new TString("FQD ..... ");
495 TString *entryLYZ1RP = new TString("LYZ ..... ");
496 TString *entryLYZEPRP = new TString("LYZEP ... ");
1fca9c90 497
8c5ec46a 498 //POI:
499 TString *entryMCPOI = new TString("MC ...... ");
500 TString *entryGFCPOI = new TString("GFC ..... ");
501 TString *entryQC2POI = new TString("QC{2} ... ");
502 TString *entryQC4POI = new TString("QC{4} ... ");
503 TString *entryQC6POI = new TString("QC{6} ... ");
504 TString *entryQC8POI = new TString("QC{8} ... ");
505 TString *entryFQDPOI = new TString("FQD ..... ");
506 TString *entryLYZ1POI = new TString("LYZ ..... ");
507 TString *entryLYZEPPOI = new TString("LYZEP ... ");
508
509 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
510 // !!!! to be removed !!!!
511 Double_t avMultMC=0.;
512 Long_t nEvtsMC=0;
513
514 Double_t avMultGFC=0.;
515 Long_t nEvtsGFC=0;
516
517 Double_t avMultQC2=0., avMultQC4=0., avMultQC6=0., avMultQC8=0.;
518 Long_t nEvtsQC2=0, nEvtsQC4=0, nEvtsQC6=0, nEvtsQC8=0;
519
520 Double_t avMultFQD=0.;
521 Long_t nEvtsFQD=0;
522
523 Double_t avMultLYZ1=0.;
524 Long_t nEvtsLYZ1=0;
525
526 Double_t avMultLYZEP=0.;
527 Long_t nEvtsLYZEP=0;
528 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
529
530 //RP:
531 Double_t avMultMCRP=0.;
532 Long_t nEvtsMCRP=0;
533
534 Double_t avMultGFCRP=0.;
535 Long_t nEvtsGFCRP=0;
536
537 Double_t avMultQC2RP=0., avMultQC4RP=0., avMultQC6RP=0., avMultQC8RP=0.;
538 Long_t nEvtsQC2RP=0, nEvtsQC4RP=0, nEvtsQC6RP=0, nEvtsQC8RP=0;
539
540 Double_t avMultFQDRP=0.;
541 Long_t nEvtsFQDRP=0;
542
543 Double_t avMultLYZ1RP=0.;
544 Long_t nEvtsLYZ1RP=0;
545
546 Double_t avMultLYZEPRP=0.;
547 Long_t nEvtsLYZEPRP=0;
548
549 //POI:
550 Double_t avMultMCPOI=0.;
551 Long_t nEvtsMCPOI=0;
552
553 Double_t avMultGFCPOI=0.;
554 Long_t nEvtsGFCPOI=0;
555
556 Double_t avMultQC2POI=0., avMultQC4POI=0., avMultQC6POI=0., avMultQC8POI=0.;
557 Long_t nEvtsQC2POI=0, nEvtsQC4POI=0, nEvtsQC6POI=0, nEvtsQC8POI=0;
558
559 Double_t avMultFQDPOI=0.;
560 Long_t nEvtsFQDPOI=0;
561
562 Double_t avMultLYZ1POI=0.;
563 Long_t nEvtsLYZ1POI=0;
564
565 Double_t avMultLYZEPPOI=0.;
566 Long_t nEvtsLYZEPPOI=0;
567
568 //MC:
569 if(mcepCommonHist) {
570 avMultMC = (mcepCommonHist->GetHistMultInt())->GetMean();//to be removed
571 nEvtsMC = (mcepCommonHist->GetHistMultInt())->GetEntries();//to be removed
572 avMultMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
573 nEvtsMCRP = (mcepCommonHist->GetHistMultInt())->GetEntries();
574 avMultMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
575 nEvtsMCPOI = (mcepCommonHist->GetHistMultDiff())->GetEntries();
576 }
577 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
578 // !!!! to be removed !!!!
579 if(entryMC) {
580 entryMC->Append("M = ");
581 (*entryMC)+=(Long_t)avMultMC;
582 entryMC->Append(", N = ");
583 (*entryMC)+=(Long_t)nEvtsMC;
584 }
1fca9c90 585 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
8c5ec46a 586 if(entryMCRP) {
587 entryMCRP->Append("M = ");
588 (*entryMCRP)+=(Long_t)avMultMCRP;
589 entryMCRP->Append(", N = ");
590 (*entryMCRP)+=(Long_t)nEvtsMCRP;
591 }
592 if(entryMCPOI) {
593 entryMCPOI->Append("M = ");
594 (*entryMCPOI)+=(Long_t)avMultMCPOI;
595 entryMCPOI->Append(", N = ");
596 (*entryMCPOI)+=(Long_t)nEvtsMCPOI;
1fca9c90 597 }
598
a88e5ac2 599 //GFC:
8c5ec46a 600 if(gfcCommonHist) {
601 avMultGFC = (gfcCommonHist->GetHistMultInt())->GetMean();//to be removed
602 nEvtsGFC = (gfcCommonHist->GetHistMultInt())->GetEntries();//to be removed
603 avMultGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();
604 nEvtsGFCRP = (gfcCommonHist->GetHistMultInt())->GetEntries();
605 avMultGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
606 nEvtsGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetEntries();
a88e5ac2 607 }
1fca9c90 608 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
609 // !!!! to be removed !!!!
8c5ec46a 610 if(entryGFC) {
611 entryGFC->Append("M = ");
612 (*entryGFC)+=(Long_t)avMultGFC;
613 entryGFC->Append(", N = ");
614 (*entryGFC)+=(Long_t)nEvtsGFC;
a88e5ac2 615 }
1fca9c90 616 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
8c5ec46a 617 if(entryGFCRP) {
618 entryGFCRP->Append("M = ");
619 (*entryGFCRP)+=(Long_t)avMultGFCRP;
620 entryGFCRP->Append(", N = ");
621 (*entryGFCRP)+=(Long_t)nEvtsGFCRP;
1fca9c90 622 }
8c5ec46a 623 if(entryGFCPOI) {
624 entryGFCPOI->Append("M = ");
625 (*entryGFCPOI)+=(Long_t)avMultGFCPOI;
626 entryGFCPOI->Append(", N = ");
627 (*entryGFCPOI)+=(Long_t)nEvtsGFCPOI;
1fca9c90 628 }
a88e5ac2 629
630 //QC:
8c5ec46a 631 if(qcCommonHist2) {
632 avMultQC2 = (qcCommonHist2->GetHistMultInt())->GetMean();//to be removed
633 nEvtsQC2 = (qcCommonHist2->GetHistMultInt())->GetEntries();//to be removed
634 avMultQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
635 nEvtsQC2RP = (qcCommonHist2->GetHistMultInt())->GetEntries();
636 avMultQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
637 nEvtsQC2POI = (qcCommonHist2->GetHistMultDiff())->GetEntries();
1fca9c90 638 }
639 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
640 // !!!! to be removed !!!!
a88e5ac2 641 if(entryQC2)
642 {
70288839 643 entryQC2->Append("M = ");
644 (*entryQC2)+=(Long_t)avMultQC2;
645 entryQC2->Append(", N = ");
646 (*entryQC2)+=(Long_t)nEvtsQC2;
70288839 647 }
1fca9c90 648 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
649 if(entryQC2RP)
650 {
651 entryQC2RP->Append("M = ");
652 (*entryQC2RP)+=(Long_t)avMultQC2RP;
653 entryQC2RP->Append(", N = ");
654 (*entryQC2RP)+=(Long_t)nEvtsQC2RP;
655 }
656 if(entryQC2POI)
657 {
658 entryQC2POI->Append("M = ");
659 (*entryQC2POI)+=(Long_t)avMultQC2POI;
660 entryQC2POI->Append(", N = ");
661 (*entryQC2POI)+=(Long_t)nEvtsQC2POI;
662 }
663
70288839 664 if(qcCommonHist4)
665 {
1fca9c90 666 avMultQC4 = (qcCommonHist4->GetHistMultInt())->GetMean();//to be removed
667 nEvtsQC4 = (qcCommonHist4->GetHistMultInt())->GetEntries();//to be removed
668 avMultQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
669 nEvtsQC4RP = (qcCommonHist4->GetHistMultInt())->GetEntries();
670 avMultQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
671 nEvtsQC4POI = (qcCommonHist4->GetHistMultDiff())->GetEntries();
a88e5ac2 672 }
1fca9c90 673 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
674 // !!!! to be removed !!!!
a88e5ac2 675 if(entryQC4)
676 {
70288839 677 entryQC4->Append("M = ");
678 (*entryQC4)+=(Long_t)avMultQC4;
679 entryQC4->Append(", N = ");
680 (*entryQC4)+=(Long_t)nEvtsQC4;
70288839 681 }
1fca9c90 682 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
683 if(entryQC4RP)
684 {
685 entryQC4RP->Append("M = ");
686 (*entryQC4RP)+=(Long_t)avMultQC4RP;
687 entryQC4RP->Append(", N = ");
688 (*entryQC4RP)+=(Long_t)nEvtsQC4RP;
689 }
690 if(entryQC4POI)
691 {
692 entryQC4POI->Append("M = ");
693 (*entryQC4POI)+=(Long_t)avMultQC4POI;
694 entryQC4POI->Append(", N = ");
695 (*entryQC4POI)+=(Long_t)nEvtsQC4POI;
696 }
697
a88e5ac2 698 if(qcCommonHist6)
699 {
1fca9c90 700 avMultQC6 = (qcCommonHist6->GetHistMultInt())->GetMean();//to be removed
701 nEvtsQC6 = (qcCommonHist6->GetHistMultInt())->GetEntries();//to be removed
702 avMultQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
703 nEvtsQC6RP = (qcCommonHist6->GetHistMultInt())->GetEntries();
704 avMultQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
705 nEvtsQC6POI = (qcCommonHist6->GetHistMultDiff())->GetEntries();
a88e5ac2 706 }
1fca9c90 707 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
708 // !!!! to be removed !!!!
a88e5ac2 709 if(entryQC6)
710 {
711 entryQC6->Append("M = ");
712 (*entryQC6)+=(Long_t)avMultQC6;
713 entryQC6->Append(", N = ");
714 (*entryQC6)+=(Long_t)nEvtsQC6;
715 }
1fca9c90 716 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
717 if(entryQC6RP)
718 {
719 entryQC6RP->Append("M = ");
720 (*entryQC6RP)+=(Long_t)avMultQC6RP;
721 entryQC6RP->Append(", N = ");
722 (*entryQC6RP)+=(Long_t)nEvtsQC6RP;
723 }
724 if(entryQC6POI)
725 {
726 entryQC6POI->Append("M = ");
727 (*entryQC6POI)+=(Long_t)avMultQC6POI;
728 entryQC6POI->Append(", N = ");
729 (*entryQC6POI)+=(Long_t)nEvtsQC6POI;
730 }
731
a88e5ac2 732 if(qcCommonHist8)
733 {
1fca9c90 734 avMultQC8 = (qcCommonHist8->GetHistMultInt())->GetMean();//to be removed
735 nEvtsQC8 = (qcCommonHist8->GetHistMultInt())->GetEntries();//to be removed
736 avMultQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
737 nEvtsQC8RP = (qcCommonHist8->GetHistMultInt())->GetEntries();
738 avMultQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
739 nEvtsQC8POI = (qcCommonHist8->GetHistMultDiff())->GetEntries();
a88e5ac2 740 }
1fca9c90 741 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
742 // !!!! to be removed !!!!
a88e5ac2 743 if(entryQC8)
744 {
745 entryQC8->Append("M = ");
746 (*entryQC8)+=(Long_t)avMultQC8;
747 entryQC8->Append(", N = ");
748 (*entryQC8)+=(Long_t)nEvtsQC8;
749 }
1fca9c90 750 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
751 if(entryQC8RP)
752 {
753 entryQC8RP->Append("M = ");
754 (*entryQC8RP)+=(Long_t)avMultQC8RP;
755 entryQC8RP->Append(", N = ");
756 (*entryQC8RP)+=(Long_t)nEvtsQC8RP;
757 }
758 if(entryQC8POI)
759 {
760 entryQC8POI->Append("M = ");
761 (*entryQC8POI)+=(Long_t)avMultQC8POI;
762 entryQC8POI->Append(", N = ");
763 (*entryQC8POI)+=(Long_t)nEvtsQC8POI;
764 }
765
a88e5ac2 766 //FQD:
767 if(fqdCommonHist)
768 {
1fca9c90 769 avMultFQD = (fqdCommonHist->GetHistMultInt())->GetMean();//to be removed
770 nEvtsFQD = (fqdCommonHist->GetHistMultInt())->GetEntries();//to be removed
771 avMultFQDRP = (fqdCommonHist->GetHistMultInt())->GetMean();
772 nEvtsFQDRP = (fqdCommonHist->GetHistMultInt())->GetEntries();
773 avMultFQDPOI = (fqdCommonHist->GetHistMultDiff())->GetMean();
774 nEvtsFQDPOI = (fqdCommonHist->GetHistMultDiff())->GetEntries();
a88e5ac2 775 }
1fca9c90 776 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
777 // !!!! to be removed !!!!
a88e5ac2 778 if(entryFQD)
779 {
780 entryFQD->Append("M = ");
781 (*entryFQD)+=(Long_t)avMultFQD;
782 entryFQD->Append(", N = ");
783 (*entryFQD)+=(Long_t)nEvtsFQD;
784 }
1fca9c90 785 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
786 if(entryFQDRP)
787 {
788 entryFQDRP->Append("M = ");
789 (*entryFQDRP)+=(Long_t)avMultFQDRP;
790 entryFQDRP->Append(", N = ");
791 (*entryFQDRP)+=(Long_t)nEvtsFQDRP;
792 }
793 if(entryFQDPOI)
794 {
795 entryFQDPOI->Append("M = ");
796 (*entryFQDPOI)+=(Long_t)avMultFQDPOI;
797 entryFQDPOI->Append(", N = ");
798 (*entryFQDPOI)+=(Long_t)nEvtsFQDPOI;
799 }
800
a88e5ac2 801 //LYZ1:
802 if(lyz1CommonHist)
803 {
1fca9c90 804 avMultLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetMean();//to be removed
805 nEvtsLYZ1 = (lyz1CommonHist->GetHistMultInt())->GetEntries();//to be removed
806 avMultLYZ1RP = (lyz1CommonHist->GetHistMultInt())->GetMean();
807 nEvtsLYZ1RP = (lyz1CommonHist->GetHistMultInt())->GetEntries();
808 avMultLYZ1POI = (lyz1CommonHist->GetHistMultDiff())->GetMean();
809 nEvtsLYZ1POI = (lyz1CommonHist->GetHistMultDiff())->GetEntries();
a88e5ac2 810 }
1fca9c90 811 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
812 // !!!! to be removed !!!!
a88e5ac2 813 if(entryLYZ1)
814 {
815 entryLYZ1->Append("M = ");
816 (*entryLYZ1)+=(Long_t)avMultLYZ1;
817 entryLYZ1->Append(", N = ");
818 (*entryLYZ1)+=(Long_t)nEvtsLYZ1;
819 }
1fca9c90 820 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
821 if(entryLYZ1RP)
822 {
823 entryLYZ1RP->Append("M = ");
824 (*entryLYZ1RP)+=(Long_t)avMultLYZ1RP;
825 entryLYZ1RP->Append(", N = ");
826 (*entryLYZ1RP)+=(Long_t)nEvtsLYZ1RP;
827 }
828 if(entryLYZ1POI)
829 {
830 entryLYZ1POI->Append("M = ");
831 (*entryLYZ1POI)+=(Long_t)avMultLYZ1POI;
832 entryLYZ1POI->Append(", N = ");
833 (*entryLYZ1POI)+=(Long_t)nEvtsLYZ1POI;
834 }
70288839 835
a88e5ac2 836 //LYZEP:
837 if(lyzepCommonHist)
70288839 838 {
e83922f2 839 avMultLYZEP = (lyzepCommonHist->GetHistMultInt())->GetMean();//to be removed
1fca9c90 840 nEvtsLYZEP = (lyzepCommonHist->GetHistMultInt())->GetEntries();//to be removed
e83922f2 841 avMultLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
1fca9c90 842 nEvtsLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetEntries();
e83922f2 843 avMultLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
1fca9c90 844 nEvtsLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetEntries();
70288839 845 }
1fca9c90 846 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
847 // !!!! to be removed !!!!
a88e5ac2 848 if(entryLYZEP)
70288839 849 {
a88e5ac2 850 entryLYZEP->Append("M = ");
851 (*entryLYZEP)+=(Long_t)avMultLYZEP;
852 entryLYZEP->Append(", N = ");
853 (*entryLYZEP)+=(Long_t)nEvtsLYZEP;
70288839 854 }
1fca9c90 855 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
856 if(entryLYZEPRP)
857 {
858 entryLYZEPRP->Append("M = ");
859 (*entryLYZEPRP)+=(Long_t)avMultLYZEPRP;
860 entryLYZEPRP->Append(", N = ");
861 (*entryLYZEPRP)+=(Long_t)nEvtsLYZEPRP;
862 }
863 if(entryLYZEPPOI)
864 {
865 entryLYZEPPOI->Append("M = ");
866 (*entryLYZEPPOI)+=(Long_t)avMultLYZEPPOI;
867 entryLYZEPPOI->Append(", N = ");
868 (*entryLYZEPPOI)+=(Long_t)nEvtsLYZEPPOI;
869 }
870
871 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
872 // !!!! to be removed !!!!
a88e5ac2 873 if(textResults)
70288839 874 {
a88e5ac2 875 textResults->AddText(entryMC->Data());
876 textResults->AddText(entryGFC->Data());
877 textResults->AddText(entryQC2->Data());
878 textResults->AddText(entryQC4->Data());
879 textResults->AddText(entryQC6->Data());
880 textResults->AddText(entryQC8->Data());
881 textResults->AddText(entryFQD->Data());
882 textResults->AddText(entryLYZ1->Data());
883 textResults->AddText(entryLYZEP->Data());
70288839 884 }
1fca9c90 885 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxx
886 //RP:
887 if(textResultsRP)
888 {
889 textResultsRP->AddText(entryMCRP->Data());
890 textResultsRP->AddText(entryGFCRP->Data());
891 textResultsRP->AddText(entryQC2RP->Data());
892 textResultsRP->AddText(entryQC4RP->Data());
893 textResultsRP->AddText(entryQC6RP->Data());
894 textResultsRP->AddText(entryQC8RP->Data());
895 textResultsRP->AddText(entryFQDRP->Data());
896 textResultsRP->AddText(entryLYZ1RP->Data());
897 textResultsRP->AddText(entryLYZEPRP->Data());
898 }
899 //POI:
900 if(textResultsPOI)
901 {
902 textResultsPOI->AddText(entryMCPOI->Data());
903 textResultsPOI->AddText(entryGFCPOI->Data());
904 textResultsPOI->AddText(entryQC2POI->Data());
905 textResultsPOI->AddText(entryQC4POI->Data());
906 textResultsPOI->AddText(entryQC6POI->Data());
907 textResultsPOI->AddText(entryQC8POI->Data());
908 textResultsPOI->AddText(entryFQDPOI->Data());
909 textResultsPOI->AddText(entryLYZ1POI->Data());
910 textResultsPOI->AddText(entryLYZEPPOI->Data());
911 }
a88e5ac2 912 //----------------------------------------------------------------------------------
70288839 913
1fca9c90 914 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
915 // !!!! to be removed !!!!
70288839 916 //----------------------------------------------------------------------------------
a88e5ac2 917 //final drawing for integrated flow:
918 TCanvas* intFlowAllCanvas = new TCanvas("Integrated Flow","Integrated Flow",1000,600);
6f62c2a7 919
a88e5ac2 920 intFlowAllCanvas->Divide(2,1);
6f62c2a7 921
a88e5ac2 922 //1st pad is for plot:
923 (intFlowAllCanvas->cd(1))->SetPad(0.0,0.0,0.75,1.0);
9d3cbee7 924
a88e5ac2 925 if(intFlowAll)
926 {
927 if(dMin>0. && dMax>0.)
928 {
929 (intFlowAll->GetYaxis())->SetRangeUser(0.9744*dMin,1.0144*dMax);
930 } else if(dMin<0. && dMax>0.)
931 {
932 if(!(-1.*dMin<4.*dMax))
933 {
934 (intFlowAll->GetYaxis())->SetRangeUser(1.0266*dMin,1.0144*dMax);
935 } else {(intFlowAll->GetYaxis())->SetRangeUser(1.1266*dMin,1.0144*dMax);}
936 } else if(dMin<0. && dMax<0.)
937 {
938 (intFlowAll->GetYaxis())->SetRangeUser(1.0266*dMin,0.9866*dMax);
939 }
940 intFlowAll->Draw("E1");
e83922f2 941 }
942
a88e5ac2 943 if(pMesh) pMesh->Draw("LFSAME");
944
945 if(flowResults) flowResults->Draw("PSAME");
946
947 //2nd pad is for legend:
948 (intFlowAllCanvas->cd(2))->SetPad(0.75,0.0,1.0,1.0);
949
950 if(textDefault)
951 {
952 textDefault->Draw();
953 textResults->Draw();
954 }
e83922f2 955
a88e5ac2 956 //----------------------------------------------------------------------------------
1fca9c90 957 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
958
e83922f2 959
960
e83922f2 961
1fca9c90 962 //----------------------------------------------------------------------------------
963 //final drawing for integrated flow of RP (i.e. of particles used to determine the reaction plane):
e83922f2 964 if(plotIntFlowRP)
965 {
966 TCanvas* intFlowAllCanvasRP = new TCanvas("Integrated Flow RP","Integrated Flow RP",1000,600);
1fca9c90 967
e83922f2 968 intFlowAllCanvasRP->Divide(2,1);
1fca9c90 969
e83922f2 970 //1st pad is for plot:
971 (intFlowAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0);
1fca9c90 972
694e0378 973 TH1D *intFlowAllRP = new TH1D(*intFlowAll);
974
e83922f2 975 if(intFlowAllRP)
976 {
977 if(dMinRP>0. && dMaxRP>0.)
978 {
979 (intFlowAllRP->GetYaxis())->SetRangeUser(0.9744*dMinRP,1.0144*dMaxRP);
980 } else if(dMinRP<0. && dMaxRP>0.)
981 {
982 if(!(-1.*dMinRP<4.*dMaxRP))
983 {
984 (intFlowAllRP->GetYaxis())->SetRangeUser(1.0266*dMinRP,1.0144*dMaxRP);
985 } else {(intFlowAllRP->GetYaxis())->SetRangeUser(1.1266*dMinRP,1.0144*dMaxRP);}
986 } else if(dMinRP<0. && dMaxRP<0.)
987 {
988 (intFlowAllRP->GetYaxis())->SetRangeUser(1.0266*dMinRP,0.9866*dMaxRP);
989 }
990 intFlowAllRP->Draw("E1");
991 }
694e0378 992
e83922f2 993 if(pMeshRP) pMeshRP->Draw("LFSAME");
994
995 if(flowResultsRP) flowResultsRP->Draw("PSAME");
1fca9c90 996
e83922f2 997 //2nd pad is for legend:
998 (intFlowAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0);
999
1000 if(textDefault) textDefault->Draw();
1fca9c90 1001
e83922f2 1002 if(textResultsRP) textResultsRP->Draw();
1003 }//end of if(plotIntFlowRP}
1fca9c90 1004 //----------------------------------------------------------------------------------
1005
1006 //----------------------------------------------------------------------------------
1007 //final drawing for integrated flow of POI (i.e. of particles of interest):
e83922f2 1008 if(plotIntFlowPOI)
1009 {
1010 TCanvas* intFlowAllCanvasPOI = new TCanvas("Integrated Flow POI","Integrated Flow POI",1000,600);
1fca9c90 1011
e83922f2 1012 intFlowAllCanvasPOI->Divide(2,1);
1fca9c90 1013
e83922f2 1014 //1st pad is for plot:
1015 (intFlowAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0);
1016
694e0378 1017 TH1D *intFlowAllPOI = new TH1D(*intFlowAll);
1018
e83922f2 1019 if(intFlowAllPOI)
1fca9c90 1020 {
e83922f2 1021 if(dMinPOI>0. && dMaxPOI>0.)
1022 {
1023 (intFlowAllPOI->GetYaxis())->SetRangeUser(0.9744*dMinPOI,1.0144*dMaxPOI);
1024 } else if(dMinPOI<0. && dMaxPOI>0.)
1025 {
1026 if(!(-1.*dMinPOI<4.*dMaxPOI))
1027 {
1028 (intFlowAllPOI->GetYaxis())->SetRangeUser(1.0266*dMinPOI,1.0144*dMaxPOI);
1029 } else {(intFlowAllPOI->GetYaxis())->SetRangeUser(1.1266*dMinPOI,1.0144*dMaxPOI);}
1030 } else if(dMinPOI<0. && dMaxPOI<0.)
1031 {
1032 (intFlowAllPOI->GetYaxis())->SetRangeUser(1.0266*dMinPOI,0.9866*dMaxPOI);
1033 }
1034 intFlowAllPOI->Draw("E1");
1035 }
1fca9c90 1036
e83922f2 1037 if(pMeshPOI) pMeshPOI->Draw("LFSAME");
1038
1039 if(flowResultsPOI) flowResultsPOI->Draw("PSAME");
1fca9c90 1040
e83922f2 1041 //2nd pad is for legend:
1042 (intFlowAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0);
1043
1044 if(textDefault) textDefault->Draw();
1fca9c90 1045
e83922f2 1046 if(textResultsPOI) textResultsPOI->Draw();
1047 }// end of if(plotIntFlowPOI)
1fca9c90 1048 //----------------------------------------------------------------------------------
1049
9d3cbee7 1050 //==================================================================================
e83922f2 1051
1052
1053
9d3cbee7 1054
1055 //==================================================================================
1056 // DIFFERENTIAL FLOW
1057 //==================================================================================
70288839 1058 Int_t iNbinsPt = AliFlowCommonConstants::GetNbinsPt();
1059 Double_t dPtMin = AliFlowCommonConstants::GetPtMin();
1060 Double_t dPtMax = AliFlowCommonConstants::GetPtMax();
1061
1fca9c90 1062 Int_t iNbinsEta = AliFlowCommonConstants::GetNbinsEta();
1063 Double_t dEtaMin = AliFlowCommonConstants::GetEtaMin();
1064 Double_t dEtaMax = AliFlowCommonConstants::GetEtaMax();
1065
a88e5ac2 1066 //----------------------------------------------------------------------------------
1fca9c90 1067 //cosmetics: the style histogram for differential flow (pt):
1068 TH1D *styleHistPt = new TH1D("styleHistPt","styleHistPt",iNbinsPt,dPtMin,dPtMax);
1069 styleHistPt->SetTitle("Differential Flow");
1070 styleHistPt->SetXTitle("p_{t} [GeV]");
1071 styleHistPt->SetYTitle("v_{n}");
1072
1073 //cosmetics: the style histogram for differential flow (eta):
1074 TH1D *styleHistEta = new TH1D("styleHistEta","styleHistEta",iNbinsEta,dEtaMin,dEtaMax);
1075 styleHistEta->SetTitle("Differential Flow");
1076 styleHistEta->SetXTitle("#eta");
1077 styleHistEta->SetYTitle("v_{n}");
a88e5ac2 1078 //----------------------------------------------------------------------------------
9d3cbee7 1079
1fca9c90 1080
1081
1082 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
1083 // !!!! to be removed !!!!
a88e5ac2 1084 //----------------------------------------------------------------------------------
1085 //cosmetics: Monte Carlo error bands for differential flow
1086 TGraph* pMeshDiffFlow = NULL;
1087 if(mcepCommonHistRes)
1088 {
1089 Int_t nPtsDF = (mcepCommonHistRes->GetHistDiffFlow())->GetNbinsX();
1090 Double_t binWidth = (mcepCommonHistRes->GetHistDiffFlow())->GetBinWidth(1);//assuming that all bins have the same width
1091
1fca9c90 1092 pMeshDiffFlow = new TGraph(2*nPtsDF+1);
a88e5ac2 1093
1094 Double_t valueMC=0., errorMC=0.;
1095 for(Int_t i=1;i<nPtsDF+1;i++)
1096 {
1097 valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(i);
1098 errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(i);
1099 pMeshDiffFlow->SetPoint(i,(i-0.5)*binWidth,valueMC+errorMC);
1100 }
1101 for(Int_t i=nPtsDF+1;i<2*nPtsDF+1;i++)
1102 {
1103 valueMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinContent(2*nPtsDF+1-i);
1104 errorMC = (mcepCommonHistRes->GetHistDiffFlow())->GetBinError(2*nPtsDF+1-i);
1105 pMeshDiffFlow->SetPoint(i,(2*nPtsDF-i+0.5)*binWidth,valueMC-errorMC);
1106 }
1107 pMeshDiffFlow->SetPoint(2*nPtsDF+1,0.5*binWidth,valueMC+errorMC);
1108 pMeshDiffFlow->SetFillStyle(meshStyle);
1109 pMeshDiffFlow->SetFillColor(meshColor);
1110 }
1111 //----------------------------------------------------------------------------------
1fca9c90 1112 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
1113
1fca9c90 1114 //----------------------------------------------------------------------------------
1115 //RP:
1116 //cosmetics: Monte Carlo error bands for differential flow (Pt)
1117 TGraph* pMeshDiffFlowPtRP = NULL;
1118 if(mcepCommonHistRes)
1119 {
1120 Int_t nPtsDFPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetNbinsX();
1121 Double_t binWidthPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinWidth(1);//assuming that all bins have the same width
1122
1123 pMeshDiffFlowPtRP = new TGraph(2*nPtsDFPtRP+1);
1124
1125 Double_t valueMCPtRP=0., errorMCPtRP=0.;
1126 for(Int_t i=1;i<nPtsDFPtRP+1;i++)
1127 {
1128 valueMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(i);
1129 errorMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinError(i);
1130 pMeshDiffFlowPtRP->SetPoint(i,(i-0.5)*binWidthPtRP,valueMCPtRP+errorMCPtRP);
1131 }
1132 for(Int_t i=nPtsDFPtRP+1;i<2*nPtsDFPtRP+1;i++)
1133 {
1134 valueMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(2*nPtsDFPtRP+1-i);
1135 errorMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinError(2*nPtsDFPtRP+1-i);
1136 pMeshDiffFlowPtRP->SetPoint(i,(2*nPtsDFPtRP-i+0.5)*binWidthPtRP,valueMCPtRP-errorMCPtRP);
1137 }
1138 pMeshDiffFlowPtRP->SetPoint(2*nPtsDFPtRP+1,0.5*binWidthPtRP,valueMCPtRP+errorMCPtRP);
1139 pMeshDiffFlowPtRP->SetFillStyle(meshStyle);
1140 pMeshDiffFlowPtRP->SetFillColor(meshColor);
1141 }
1142
1143 //cosmetics: Monte Carlo error bands for differential flow (Eta)
1144 TGraph* pMeshDiffFlowEtaRP = NULL;
1145 if(mcepCommonHistRes)
1146 {
1147 Int_t nPtsDFEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetNbinsX();
1148 Double_t binWidthEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinWidth(1);//assuming that all bins have the same width
1149
1150 pMeshDiffFlowEtaRP = new TGraph(2*nPtsDFEtaRP+1);
e83922f2 1151
1fca9c90 1152 Double_t valueMCEtaRP=0., errorMCEtaRP=0.;
1153 for(Int_t i=1;i<nPtsDFEtaRP+1;i++)
1154 {
1155 valueMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinContent(i);
e83922f2 1156 errorMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinError(i);
1157 pMeshDiffFlowEtaRP->SetPoint(i,(i-0.5)*binWidthEtaRP+dEtaMin,valueMCEtaRP+errorMCEtaRP);
1fca9c90 1158 }
1159 for(Int_t i=nPtsDFEtaRP+1;i<2*nPtsDFEtaRP+1;i++)
1160 {
1161 valueMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinContent(2*nPtsDFEtaRP+1-i);
1162 errorMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinError(2*nPtsDFEtaRP+1-i);
e83922f2 1163 pMeshDiffFlowEtaRP->SetPoint(i,(2*nPtsDFEtaRP-i+0.5)*binWidthEtaRP+dEtaMin,valueMCEtaRP-errorMCEtaRP);
1fca9c90 1164 }
e83922f2 1165 pMeshDiffFlowEtaRP->SetPoint(2*nPtsDFEtaRP+1,0.5*binWidthEtaRP+dEtaMin,valueMCEtaRP+errorMCEtaRP);
1fca9c90 1166 pMeshDiffFlowEtaRP->SetFillStyle(meshStyle);
1167 pMeshDiffFlowEtaRP->SetFillColor(meshColor);
1168 }
1169 //----------------------------------------------------------------------------------
1170
1171
1172 //----------------------------------------------------------------------------------
1173 //POI:
1174 //cosmetics: Monte Carlo error bands for differential flow (Pt)
1175 TGraph* pMeshDiffFlowPtPOI = NULL;
1176 if(mcepCommonHistRes)
1177 {
1178 Int_t nPtsDFPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetNbinsX();
1179 Double_t binWidthPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinWidth(1);//assuming that all bins have the same width
1180
1181 pMeshDiffFlowPtPOI = new TGraph(2*nPtsDFPtPOI+1);
1182
1183 Double_t valueMCPtPOI=0., errorMCPtPOI=0.;
1184 for(Int_t i=1;i<nPtsDFPtPOI+1;i++)
1185 {
1186 valueMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinContent(i);
1187 errorMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinError(i);
1188 pMeshDiffFlowPtPOI->SetPoint(i,(i-0.5)*binWidthPtPOI,valueMCPtPOI+errorMCPtPOI);
1189 }
1190 for(Int_t i=nPtsDFPtPOI+1;i<2*nPtsDFPtPOI+1;i++)
1191 {
1192 valueMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinContent(2*nPtsDFPtPOI+1-i);
1193 errorMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinError(2*nPtsDFPtPOI+1-i);
1194 pMeshDiffFlowPtPOI->SetPoint(i,(2*nPtsDFPtPOI-i+0.5)*binWidthPtPOI,valueMCPtPOI-errorMCPtPOI);
1195 }
1196 pMeshDiffFlowPtPOI->SetPoint(2*nPtsDFPtPOI+1,0.5*binWidthPtPOI,valueMCPtPOI+errorMCPtPOI);
1197 pMeshDiffFlowPtPOI->SetFillStyle(meshStyle);
1198 pMeshDiffFlowPtPOI->SetFillColor(meshColor);
1199 }
1200
1201 //cosmetics: Monte Carlo error bands for differential flow (Eta)
1202 TGraph* pMeshDiffFlowEtaPOI = NULL;
1203 if(mcepCommonHistRes)
1204 {
1205 Int_t nPtsDFEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetNbinsX();
1206 Double_t binWidthEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinWidth(1);//assuming that all bins have the same width
1207
1208 pMeshDiffFlowEtaPOI = new TGraph(2*nPtsDFEtaPOI+1);
a88e5ac2 1209
1fca9c90 1210 Double_t valueMCEtaPOI=0., errorMCEtaPOI=0.;
1211 for(Int_t i=1;i<nPtsDFEtaPOI+1;i++)
1212 {
1213 valueMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinContent(i);
1214 errorMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinError(i);
e83922f2 1215 pMeshDiffFlowEtaPOI->SetPoint(i,(i-0.5)*binWidthEtaPOI+dEtaMin,valueMCEtaPOI+errorMCEtaPOI);
1fca9c90 1216 }
1217 for(Int_t i=nPtsDFEtaPOI+1;i<2*nPtsDFEtaPOI+1;i++)
1218 {
1219 valueMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinContent(2*nPtsDFEtaPOI+1-i);
1220 errorMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinError(2*nPtsDFEtaPOI+1-i);
e83922f2 1221 pMeshDiffFlowEtaPOI->SetPoint(i,(2*nPtsDFEtaPOI-i+0.5)*binWidthEtaPOI+dEtaMin,valueMCEtaPOI-errorMCEtaPOI);
1fca9c90 1222 }
e83922f2 1223 pMeshDiffFlowEtaPOI->SetPoint(2*nPtsDFEtaPOI+1,0.5*binWidthEtaPOI+dEtaMin,valueMCEtaPOI+errorMCEtaPOI);
1fca9c90 1224 pMeshDiffFlowEtaPOI->SetFillStyle(meshStyle);
1225 pMeshDiffFlowEtaPOI->SetFillColor(meshColor);
1226 }
1227 //----------------------------------------------------------------------------------
1228
70288839 1229 //MCEP = Monte Carlo Event Plane
1fca9c90 1230 Double_t avMultDiffFlowMC=0.;//to be removed
1231 Double_t nEvtsDiffFlowMC=0;//to be removed
1232 Double_t avMultDiffFlowMCRP=0.;
1233 Double_t nEvtsDiffFlowMCRP=0;
1234 Double_t avMultDiffFlowMCPOI=0.;
1235 Double_t nEvtsDiffFlowMCPOI=0;
a88e5ac2 1236 if(fileMCEP)
9d3cbee7 1237 {
70288839 1238 if(mcepCommonHistRes)
1239 {
1fca9c90 1240 (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(2);//to be removed
1241 (mcepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(20);//to be removed
1242 (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(2);
1243 (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(20);
1244 (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(2);
1245 (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(20);
1246 (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(2);
1247 (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(20);
1248 (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(2);
1249 (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(20);
a88e5ac2 1250 }
1251 if(mcepCommonHist)
1252 {
1fca9c90 1253 avMultDiffFlowMC = (mcepCommonHist->GetHistMultDiff())->GetMean();//to be removed
1254 nEvtsDiffFlowMC = (mcepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
1255 avMultDiffFlowMCRP = (mcepCommonHist->GetHistMultInt())->GetMean();
1256 nEvtsDiffFlowMCRP = (mcepCommonHist->GetHistMultInt())->GetEntries();
1257 avMultDiffFlowMCPOI = (mcepCommonHist->GetHistMultDiff())->GetMean();
1258 nEvtsDiffFlowMCPOI = (mcepCommonHist->GetHistMultDiff())->GetEntries();
70288839 1259 }
1260 }
a88e5ac2 1261
9d3cbee7 1262 //GFC = Generating Function Cumulants
1fca9c90 1263 Double_t avMultDiffFlowGFC=0.;//to be removed
1264 Double_t nEvtsDiffFlowGFC=0.;//to be removed
1265 Double_t avMultDiffFlowGFCRP=0.;
1266 Double_t nEvtsDiffFlowGFCRP=0.;
1267 Double_t avMultDiffFlowGFCPOI=0.;
1268 Double_t nEvtsDiffFlowGFCPOI=0.;
a88e5ac2 1269 if(fileGFC)
9d3cbee7 1270 {
a88e5ac2 1271 if(gfcCommonHistRes2)
9d3cbee7 1272 {
e83922f2 1273 (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(kViolet+7);//to be removed
1fca9c90 1274 (gfcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(20);//to be removed
e83922f2 1275 (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(kViolet+7);
1fca9c90 1276 (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(20);
e83922f2 1277 (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(kViolet+7);
1fca9c90 1278 (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerStyle(20);
e83922f2 1279 (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerColor(kViolet+7);
1fca9c90 1280 (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerStyle(20);
e83922f2 1281 (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerColor(kViolet+7);
1fca9c90 1282 (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerStyle(20);
9d3cbee7 1283 }
a88e5ac2 1284 if(gfcCommonHistRes4)
9d3cbee7 1285 {
e83922f2 1286 (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(kViolet-3);//to be removed
1fca9c90 1287 (gfcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(21);//to be removed
e83922f2 1288 (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(kViolet-3);
1fca9c90 1289 (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(21);
e83922f2 1290 (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(kViolet-3);
1fca9c90 1291 (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerStyle(21);
e83922f2 1292 (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerColor(kViolet-3);
1fca9c90 1293 (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerStyle(21);
e83922f2 1294 (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerColor(kViolet-3);
1fca9c90 1295 (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerStyle(21);
1296 }
1297 if(gfcCommonHistRes6)
1298 {
e83922f2 1299 (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(kViolet+7);//to be removed
1fca9c90 1300 (gfcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(24);//to be removed
e83922f2 1301 (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(kViolet+7);
1fca9c90 1302 (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(24);
e83922f2 1303 (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(kViolet+7);
1fca9c90 1304 (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerStyle(24);
e83922f2 1305 (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerColor(kViolet+7);
1fca9c90 1306 (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerStyle(24);
e83922f2 1307 (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerColor(kViolet+7);
1fca9c90 1308 (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerStyle(24);
1309 }
1310 if(gfcCommonHistRes8)
1311 {
e83922f2 1312 (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(kViolet-3);//to be removed
1fca9c90 1313 (gfcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(25);//to be removed
e83922f2 1314 (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(kViolet-3);
1fca9c90 1315 (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(25);
e83922f2 1316 (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(kViolet-3);
1fca9c90 1317 (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerStyle(25);
e83922f2 1318 (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerColor(kViolet-3);
1fca9c90 1319 (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerStyle(25);
e83922f2 1320 (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerColor(kViolet-3);
1fca9c90 1321 (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerStyle(25);
9d3cbee7 1322 }
a88e5ac2 1323 if(gfcCommonHist)
1324 {
1fca9c90 1325 avMultDiffFlowGFC = (gfcCommonHist->GetHistMultDiff())->GetMean();//to be removed
1326 nEvtsDiffFlowGFC = (gfcCommonHist->GetHistMultDiff())->GetEntries();//to be removed
1327 avMultDiffFlowGFCRP = (gfcCommonHist->GetHistMultInt())->GetMean();
1328 nEvtsDiffFlowGFCRP = (gfcCommonHist->GetHistMultInt())->GetEntries();
1329 avMultDiffFlowGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetMean();
1330 nEvtsDiffFlowGFCPOI = (gfcCommonHist->GetHistMultDiff())->GetEntries();
a88e5ac2 1331 }
9d3cbee7 1332 }
a88e5ac2 1333
9d3cbee7 1334 //QC = Q-cumulants
1fca9c90 1335 Double_t avMultDiffFlowQC2=0., avMultDiffFlowQC4=0.;//to be removed
1336 Double_t nEvtsDiffFlowQC2=0., nEvtsDiffFlowQC4=0.;//to be removed
e83922f2 1337 Double_t avMultDiffFlowQC6=0., avMultDiffFlowQC8=0.;//to be removed
1338 Double_t nEvtsDiffFlowQC6=0., nEvtsDiffFlowQC8=0.;//to be removed
1fca9c90 1339 Double_t avMultDiffFlowQC2RP=0.;
1340 Double_t nEvtsDiffFlowQC2RP=0.;
1341 Double_t avMultDiffFlowQC2POI=0.;
1342 Double_t nEvtsDiffFlowQC2POI=0.;
1343 Double_t avMultDiffFlowQC4RP=0.;
1344 Double_t nEvtsDiffFlowQC4RP=0.;
1345 Double_t avMultDiffFlowQC4POI=0.;
1346 Double_t nEvtsDiffFlowQC4POI=0.;
e83922f2 1347 Double_t avMultDiffFlowQC6RP=0.;
1348 Double_t nEvtsDiffFlowQC6RP=0.;
1349 Double_t avMultDiffFlowQC6POI=0.;
1350 Double_t nEvtsDiffFlowQC6POI=0.;
1351 Double_t avMultDiffFlowQC8RP=0.;
1352 Double_t nEvtsDiffFlowQC8RP=0.;
1353 Double_t avMultDiffFlowQC8POI=0.;
1354 Double_t nEvtsDiffFlowQC8POI=0.;
1fca9c90 1355
a88e5ac2 1356 if(fileQC)
9d3cbee7 1357 {
a88e5ac2 1358 //QC{2}
70288839 1359 if(qcCommonHistRes2)
9d3cbee7 1360 {
1fca9c90 1361 (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerColor(kOrange+3);//to be removed
1362 (qcCommonHistRes2->GetHistDiffFlow())->SetMarkerStyle(20);//to be removed
1363 (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(kOrange+3);
1364 (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(20);
1365 (qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(kOrange+3);
1366 (qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerStyle(20);
1367 (qcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerColor(kOrange+3);
1368 (qcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerStyle(20);
1369 (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerColor(kOrange+3);
1370 (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerStyle(20);
9d3cbee7 1371 }
a88e5ac2 1372 if(qcCommonHist2)
1373 {
1fca9c90 1374 avMultDiffFlowQC2 = (qcCommonHist2->GetHistMultDiff())->GetMean();//to be removed
1375 nEvtsDiffFlowQC2 = (qcCommonHist2->GetHistMultDiff())->GetEntries();//to be removed
1376 avMultDiffFlowQC2RP = (qcCommonHist2->GetHistMultInt())->GetMean();
1377 nEvtsDiffFlowQC2RP = (qcCommonHist2->GetHistMultInt())->GetEntries();
1378 avMultDiffFlowQC2POI = (qcCommonHist2->GetHistMultDiff())->GetMean();
1379 nEvtsDiffFlowQC2POI = (qcCommonHist2->GetHistMultDiff())->GetEntries();
a88e5ac2 1380 }
1381 //QC{4}
70288839 1382 if(qcCommonHistRes4)
9d3cbee7 1383 {
1fca9c90 1384 (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerColor(kOrange-6);//to be removed
1385 (qcCommonHistRes4->GetHistDiffFlow())->SetMarkerStyle(21);//to be removed
1386 (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(kOrange-6);
1387 (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(21);
1388 (qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(kOrange-6);
1389 (qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerStyle(21);
1390 (qcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerColor(kOrange-6);
1391 (qcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerStyle(21);
1392 (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerColor(kOrange-6);
1393 (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerStyle(21);
a88e5ac2 1394 }
1395 if(qcCommonHist4)
1396 {
1fca9c90 1397 avMultDiffFlowQC4 = (qcCommonHist4->GetHistMultDiff())->GetMean();//to be removed
1398 nEvtsDiffFlowQC4 = (qcCommonHist4->GetHistMultDiff())->GetEntries();//to be removed
1399 avMultDiffFlowQC4RP = (qcCommonHist4->GetHistMultInt())->GetMean();
1400 nEvtsDiffFlowQC4RP = (qcCommonHist4->GetHistMultInt())->GetEntries();
1401 avMultDiffFlowQC4POI = (qcCommonHist4->GetHistMultDiff())->GetMean();
1402 nEvtsDiffFlowQC4POI = (qcCommonHist4->GetHistMultDiff())->GetEntries();
9d3cbee7 1403 }
e83922f2 1404 //QC{6}
1405 if(qcCommonHistRes6)
1406 {
1407 (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerColor(kOrange+3);//to be removed
1408 (qcCommonHistRes6->GetHistDiffFlow())->SetMarkerStyle(24);//to be removed
1409 (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(kOrange+3);
1410 (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(24);
1411 (qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(kOrange+3);
1412 (qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerStyle(24);
1413 (qcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerColor(kOrange+3);
1414 (qcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerStyle(24);
1415 (qcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerColor(kOrange+3);
1416 (qcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerStyle(24);
1417 }
1418 if(qcCommonHist6)
1419 {
1420 avMultDiffFlowQC6 = (qcCommonHist6->GetHistMultDiff())->GetMean();//to be removed
1421 nEvtsDiffFlowQC6 = (qcCommonHist6->GetHistMultDiff())->GetEntries();//to be removed
1422 avMultDiffFlowQC6RP = (qcCommonHist6->GetHistMultInt())->GetMean();
1423 nEvtsDiffFlowQC6RP = (qcCommonHist6->GetHistMultInt())->GetEntries();
1424 avMultDiffFlowQC6POI = (qcCommonHist6->GetHistMultDiff())->GetMean();
1425 nEvtsDiffFlowQC6POI = (qcCommonHist6->GetHistMultDiff())->GetEntries();
1426 }
1427 //QC{8}
1428 if(qcCommonHistRes8)
1429 {
1430 (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerColor(kOrange-6);//to be removed
1431 (qcCommonHistRes8->GetHistDiffFlow())->SetMarkerStyle(25);//to be removed
1432 (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(kOrange-6);
1433 (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(25);
1434 (qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(kOrange-6);
1435 (qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerStyle(25);
1436 (qcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerColor(kOrange-6);
1437 (qcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerStyle(25);
1438 (qcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerColor(kOrange-6);
1439 (qcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerStyle(25);
1440 }
1441 if(qcCommonHist8)
1442 {
1443 avMultDiffFlowQC8 = (qcCommonHist8->GetHistMultDiff())->GetMean();//to be removed
1444 nEvtsDiffFlowQC8 = (qcCommonHist8->GetHistMultDiff())->GetEntries();//to be removed
1445 avMultDiffFlowQC8RP = (qcCommonHist8->GetHistMultInt())->GetMean();
1446 nEvtsDiffFlowQC8RP = (qcCommonHist8->GetHistMultInt())->GetEntries();
1447 avMultDiffFlowQC8POI = (qcCommonHist8->GetHistMultDiff())->GetMean();
1448 nEvtsDiffFlowQC8POI = (qcCommonHist8->GetHistMultDiff())->GetEntries();
1449 }
9d3cbee7 1450 }
70288839 1451
a88e5ac2 1452 //LYZ2 = Lee-Yang Zeros (2nd run)
1fca9c90 1453 Double_t avMultDiffFlowLYZ2=0.;//to be removed
1454 Double_t nEvtsDiffFlowLYZ2=0;//to be removed
1455 Double_t avMultDiffFlowLYZ2RP=0.;
1456 Double_t nEvtsDiffFlowLYZ2RP=0;
1457 Double_t avMultDiffFlowLYZ2POI=0.;
1458 Double_t nEvtsDiffFlowLYZ2POI=0;
a88e5ac2 1459 if(fileLYZ2)
9d3cbee7 1460 {
a88e5ac2 1461 if(lyz2CommonHistRes)
9d3cbee7 1462 {
1fca9c90 1463 (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerColor(kGreen+3);//to be removed
1464 (lyz2CommonHistRes->GetHistDiffFlow())->SetMarkerStyle(22);//to be removed
1465 (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(kGreen+3);
1466 (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(22);
1467 (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(kGreen+3);
1468 (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(22);
1469 (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(kGreen+3);
1470 (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(22);
1471 (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(kGreen+3);
1472 (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(22);
a88e5ac2 1473 }
1474 if(lyz2CommonHist)
9d3cbee7 1475 {
1fca9c90 1476 avMultDiffFlowLYZ2 = (lyz2CommonHist->GetHistMultDiff())->GetMean();//to be removed
1477 nEvtsDiffFlowLYZ2 = (lyz2CommonHist->GetHistMultDiff())->GetEntries();//to be removed
1478 avMultDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultInt())->GetMean();
1479 nEvtsDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultInt())->GetEntries();
1480 avMultDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultDiff())->GetMean();
1481 nEvtsDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultDiff())->GetEntries();
9d3cbee7 1482 }
1483 }
70288839 1484
e83922f2 1485 //LYZEP = Lee-Yang Zeros Event Plane
1486 Double_t avMultDiffFlowLYZEP=0.;//to be removed
1487 Double_t nEvtsDiffFlowLYZEP=0;//to be removed
1488 Double_t avMultDiffFlowLYZEPRP=0.;
1489 Double_t nEvtsDiffFlowLYZEPRP=0;
1490 Double_t avMultDiffFlowLYZEPPOI=0.;
1491 Double_t nEvtsDiffFlowLYZEPPOI=0;
1492 if(fileLYZEP)
1493 {
1494 if(lyzepCommonHistRes)
1495 {
1496 (lyzepCommonHistRes->GetHistDiffFlow())->Scale(0.01);//to be improved
1497 (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerColor(kGreen+3);//to be removed
1498 (lyzepCommonHistRes->GetHistDiffFlow())->SetMarkerStyle(26);//to be removed
1499 (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(kGreen+3);
1500 (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(26);
1501 (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(kGreen+3);
1502 (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(26);
1503 (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(kGreen+3);
1504 (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(26);
1505 (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(kGreen+3);
1506 (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(26);
1507 }
1508 if(lyzepCommonHist)
1509 {
1510 avMultDiffFlowLYZEP = (lyzepCommonHist->GetHistMultDiff())->GetMean();//to be removed
1511 nEvtsDiffFlowLYZEP = (lyzepCommonHist->GetHistMultDiff())->GetEntries();//to be removed
1512 avMultDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetMean();
1513 nEvtsDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultInt())->GetEntries();
1514 avMultDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetMean();
1515 nEvtsDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultDiff())->GetEntries();
1516 }
1517 }
1518
1fca9c90 1519 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
1520 // !!!! to be removed !!!!
a88e5ac2 1521 //----------------------------------------------------------------------------------
1522 //final drawing for differential flow:
1523 TCanvas* diffFlowAllCanvas = new TCanvas("Differential Flow","Differential Flow",1000,600);
70288839 1524
a88e5ac2 1525 diffFlowAllCanvas->Divide(2,1);
70288839 1526
a88e5ac2 1527 //1st pad is for plot:
1528 (diffFlowAllCanvas->cd(1))->SetPad(0.0,0.0,0.75,1.0);
1fca9c90 1529
1530 if(styleHistPt)
70288839 1531 {
1fca9c90 1532 styleHistPt->Draw();
70288839 1533 }
1fca9c90 1534
a88e5ac2 1535 if(pMeshDiffFlow)
70288839 1536 {
a88e5ac2 1537 pMeshDiffFlow->Draw("LFSAME");
70288839 1538 }
a88e5ac2 1539 //MC
1540 if(mcepCommonHistRes)
1541 {
1fca9c90 1542 (mcepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
a88e5ac2 1543 }
1544 //GFC
1545 if(gfcCommonHistRes2)
1546 {
1547 (gfcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME");
1548 }
1549 if(gfcCommonHistRes4)
1550 {
1551 (gfcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME");
1552 }
1fca9c90 1553 if(gfcCommonHistRes6)
1554 {
1555 (gfcCommonHistRes6->GetHistDiffFlow())->Draw("E1PSAME");
1556 }
1557 if(gfcCommonHistRes8)
1558 {
1559 (gfcCommonHistRes8->GetHistDiffFlow())->Draw("E1PSAME");
1560 }
a88e5ac2 1561 //QC
1562 if(qcCommonHistRes2)
1563 {
1564 (qcCommonHistRes2->GetHistDiffFlow())->Draw("E1PSAME");
1565 }
1566 if(qcCommonHistRes4)
1567 {
1568 (qcCommonHistRes4->GetHistDiffFlow())->Draw("E1PSAME");
1569 }
1570 //LYZ2
1571 if(lyz2CommonHistRes)
1572 {
1573 (lyz2CommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
70288839 1574 }
e83922f2 1575 //LYZEP
1576 if(lyzepCommonHistRes)
1577 {
1578 (lyzepCommonHistRes->GetHistDiffFlow())->Draw("E1PSAME");
1579 }
70288839 1580
a88e5ac2 1581 //2nd pad is for legend:
1582 (diffFlowAllCanvas->cd(2))->SetPad(0.75,0.0,1.0,1.0);
1583
1584 TLegend* legendDiffFlow = new TLegend(0.02,0.25,0.97,0.75);
1585 legendDiffFlow->SetTextFont(72);
1586 legendDiffFlow->SetTextSize(0.06);
1587
1588 //legend's entries:
1589 TString *entryDiffMC = new TString("MC ....... ");
1590 TString *entryDiffGFC2 = new TString("GFC{2} ... ");
1591 TString *entryDiffGFC4 = new TString("GFC{4} ... ");
1fca9c90 1592 TString *entryDiffGFC6 = new TString("GFC{6} ... ");
1593 TString *entryDiffGFC8 = new TString("GFC{8} ... ");
a88e5ac2 1594 TString *entryDiffQC2 = new TString("QC{2} .... ");
1595 TString *entryDiffQC4 = new TString("QC{4} .... ");
1596 TString *entryDiffLYZ2 = new TString("LYZ ...... ");
e83922f2 1597 TString *entryDiffLYZEP = new TString("LYZEP ... ");
a88e5ac2 1598
1599 //MC
70288839 1600 if(mcepCommonHistRes)
6f62c2a7 1601 {
a88e5ac2 1602 (mcepCommonHistRes->GetHistDiffFlow())->SetFillStyle(meshStyle);
1603 (mcepCommonHistRes->GetHistDiffFlow())->SetFillColor(meshColor);
1604 entryDiffMC->Append("M = ");
1605 (*entryDiffMC)+=(Long_t)avMultDiffFlowMC;
1606 entryDiffMC->Append(", N = ");
1607 (*entryDiffMC)+=(Long_t)nEvtsDiffFlowMC;
1608 legendDiffFlow->AddEntry(mcepCommonHistRes->GetHistDiffFlow(),entryDiffMC->Data(),"f");
1609 }
6f62c2a7 1610
a88e5ac2 1611 //GFC
1612 if(gfcCommonHistRes2)
1613 {
1614 entryDiffGFC2->Append("M = ");
1615 (*entryDiffGFC2)+=(Long_t)avMultDiffFlowGFC;
1616 entryDiffGFC2->Append(", N = ");
1617 (*entryDiffGFC2)+=(Long_t)nEvtsDiffFlowGFC;
1618 legendDiffFlow->AddEntry(gfcCommonHistRes2->GetHistDiffFlow(),entryDiffGFC2->Data(),"p");
1619 }
1620 if(gfcCommonHistRes4)
1621 {
1622 entryDiffGFC4->Append("M = ");
1623 (*entryDiffGFC4)+=(Long_t)avMultDiffFlowGFC;
1624 entryDiffGFC4->Append(", N = ");
1625 (*entryDiffGFC4)+=(Long_t)nEvtsDiffFlowGFC;
1626 legendDiffFlow->AddEntry(gfcCommonHistRes4->GetHistDiffFlow(),entryDiffGFC4->Data(),"p");
1627 }
1fca9c90 1628 if(gfcCommonHistRes6)
1629 {
1630 entryDiffGFC6->Append("M = ");
1631 (*entryDiffGFC6)+=(Long_t)avMultDiffFlowGFC;
1632 entryDiffGFC6->Append(", N = ");
1633 (*entryDiffGFC6)+=(Long_t)nEvtsDiffFlowGFC;
1634 legendDiffFlow->AddEntry(gfcCommonHistRes6->GetHistDiffFlow(),entryDiffGFC6->Data(),"p");
1635 }
1636 if(gfcCommonHistRes8)
1637 {
1638 entryDiffGFC8->Append("M = ");
1639 (*entryDiffGFC8)+=(Long_t)avMultDiffFlowGFC;
1640 entryDiffGFC8->Append(", N = ");
1641 (*entryDiffGFC8)+=(Long_t)nEvtsDiffFlowGFC;
1642 legendDiffFlow->AddEntry(gfcCommonHistRes8->GetHistDiffFlow(),entryDiffGFC8->Data(),"p");
1643 }
a88e5ac2 1644 //QC
1645 if(qcCommonHistRes2)
1646 {
1647 entryDiffQC2->Append("M = ");
1648 (*entryDiffQC2)+=(Long_t)avMultDiffFlowQC2;
1649 entryDiffQC2->Append(", N = ");
1650 (*entryDiffQC2)+=(Long_t)nEvtsDiffFlowQC2;
1651 legendDiffFlow->AddEntry(qcCommonHistRes2->GetHistDiffFlow(),entryDiffQC2->Data(),"p");
1652 }
1653 if(qcCommonHistRes4)
1654 {
1655 entryDiffQC4->Append("M = ");
1656 (*entryDiffQC4)+=(Long_t)avMultDiffFlowQC4;
1657 entryDiffQC4->Append(", N = ");
1658 (*entryDiffQC4)+=(Long_t)nEvtsDiffFlowQC4;
1659 legendDiffFlow->AddEntry(qcCommonHistRes4->GetHistDiffFlow(),entryDiffQC4->Data(),"p");
1660 }
6f62c2a7 1661
a88e5ac2 1662 //LYZ
1663 if(lyz2CommonHistRes)
1664 {
1665 entryDiffLYZ2->Append("M = ");
1666 (*entryDiffLYZ2)+=(Long_t)avMultDiffFlowLYZ2;
1667 entryDiffLYZ2->Append(", N = ");
1668 (*entryDiffLYZ2)+=(Long_t)nEvtsDiffFlowLYZ2;
1669 legendDiffFlow->AddEntry(lyz2CommonHistRes->GetHistDiffFlow(),entryDiffLYZ2->Data(),"p");
1670 }
e83922f2 1671
1672 //LYZEP
1673 if(lyzepCommonHistRes)
1674 {
1675 entryDiffLYZEP->Append("M = ");
1676 (*entryDiffLYZEP)+=(Long_t)avMultDiffFlowLYZEP;
1677 entryDiffLYZEP->Append(", N = ");
1678 (*entryDiffLYZEP)+=(Long_t)nEvtsDiffFlowLYZEP;
1679 legendDiffFlow->AddEntry(lyzepCommonHistRes->GetHistDiffFlow(),entryDiffLYZEP->Data(),"p");
1680 }
1681
a88e5ac2 1682
1683 //drawing finally the legend in the 2nd pad:
1684 if(legendDiffFlow)
1685 {
1fca9c90 1686 legendDiffFlow->SetMargin(0.15);
a88e5ac2 1687 legendDiffFlow->Draw();
1688 }
1fca9c90 1689 //----------------------------------------------------------------------------------
1690 //xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
1691
1692
1693 //----------------------------------------------------------------------------------
1694 //final drawing for differential flow (Pt, RP):
e83922f2 1695 //set here the results of which methods will be plotted by default:
1696 Bool_t plotMCPtRP = kFALSE;
1697 Bool_t plotGFC2PtRP = kTRUE;
1698 Bool_t plotGFC4PtRP = kTRUE;
1699 Bool_t plotGFC6PtRP = kTRUE;
1700 Bool_t plotGFC8PtRP = kTRUE;
1701 Bool_t plotQC2PtRP = kTRUE;
1702 Bool_t plotQC4PtRP = kTRUE;
1703 Bool_t plotQC6PtRP = kTRUE;
1704 Bool_t plotQC8PtRP = kTRUE;
1705 Bool_t plotLYZ2PtRP = kTRUE;
1706 Bool_t plotLYZEPPtRP = kTRUE;
1707 if(plotDiffFlowPtRP)
1708 {
1709 TCanvas* diffFlowPtAllCanvasRP = new TCanvas("Differential Flow (Pt) of RP","Differential Flow (Pt) of RP ",1000,600);
1710
1711 diffFlowPtAllCanvasRP->Divide(2,1);
1fca9c90 1712
e83922f2 1713 //1st pad is for plot:
1714 (diffFlowPtAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0);
6f62c2a7 1715
e83922f2 1716 if(styleHistPt)
1717 {
1718 (styleHistPt->GetYaxis())->SetRangeUser(-0.3,1.0);
1719 styleHistPt->Draw();
1720 }
1721 if(pMeshDiffFlowPtRP)
1722 {
1723 pMeshDiffFlowPtRP->Draw("LFSAME");
1724 }
1fca9c90 1725
e83922f2 1726 //MC
1727 if(plotMCPtRP && mcepCommonHistRes)
1728 {
1729 (mcepCommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1730 }
1731 //GFC
1732 if(plotGFC2PtRP && gfcCommonHistRes2)
1733 {
1734 (gfcCommonHistRes2->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1735 }
1736 if(plotGFC4PtRP && gfcCommonHistRes4)
1737 {
1738 (gfcCommonHistRes4->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1739 }
1740 if(plotGFC6PtRP && gfcCommonHistRes6)
1741 {
1742 (gfcCommonHistRes6->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1743 }
1744 if(plotGFC8PtRP && gfcCommonHistRes8)
1745 {
1746 (gfcCommonHistRes8->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1747 }
1748 //QC
1749 if(plotQC2PtRP && qcCommonHistRes2)
1750 {
1751 (qcCommonHistRes2->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1752 }
1753 if(plotQC4PtRP && qcCommonHistRes4)
1754 {
1755 (qcCommonHistRes4->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1756 }
1757 if(plotQC6PtRP && qcCommonHistRes6)
1758 {
1759 (qcCommonHistRes6->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1760 }
1761 if(plotQC8PtRP && qcCommonHistRes8)
1762 {
1763 (qcCommonHistRes8->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1764 }
1765 //LYZ2
1766 if(plotLYZ2PtRP && lyz2CommonHistRes)
1767 {
1768 (lyz2CommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1769 }
1770 //LYZEP
1771 if(plotLYZEPPtRP && lyzepCommonHistRes)
1772 {
1773 (lyzepCommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME");
1774 }
1fca9c90 1775
e83922f2 1776 //2nd pad is for legend:
1777 (diffFlowPtAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0);
1778
1779 TLegend* legendDiffFlowPtRP = new TLegend(0.02,0.25,0.97,0.75);
1780 legendDiffFlowPtRP->SetTextFont(72);
1781 legendDiffFlowPtRP->SetTextSize(0.06);
1782
1783 //legend's entries:
1784 TString *entryDiffMCPtRP = new TString("MC ....... ");
1785 TString *entryDiffGFC2PtRP = new TString("GFC{2} ... ");
1786 TString *entryDiffGFC4PtRP = new TString("GFC{4} ... ");
1787 TString *entryDiffGFC6PtRP = new TString("GFC{6} ... ");
1788 TString *entryDiffGFC8PtRP = new TString("GFC{8} ... ");
1789 TString *entryDiffQC2PtRP = new TString("QC{2} .... ");
1790 TString *entryDiffQC4PtRP = new TString("QC{4} .... ");
1791 TString *entryDiffQC6PtRP = new TString("QC{6} .... ");
1792 TString *entryDiffQC8PtRP = new TString("QC{8} .... ");
1793 TString *entryDiffLYZ2PtRP = new TString("LYZ ...... ");
1794 TString *entryDiffLYZEPPtRP = new TString("LYZEP ... ");
1795
1796 //MC
1797 if(mcepCommonHistRes)
1798 {
1799 (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillStyle(meshStyle);
1800 (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillColor(meshColor);
1801 entryDiffMCPtRP->Append("M = ");
1802 (*entryDiffMCPtRP)+=(Long_t)avMultDiffFlowMCRP;
1803 entryDiffMCPtRP->Append(", N = ");
1804 (*entryDiffMCPtRP)+=(Long_t)nEvtsDiffFlowMCRP;
1805 legendDiffFlowPtRP->AddEntry(mcepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffMCPtRP->Data(),"f");
1806 }
1fca9c90 1807
e83922f2 1808 //GFC
1809 if(plotGFC2PtRP && gfcCommonHistRes2)
1810 {
1811 entryDiffGFC2PtRP->Append("M = ");
1812 (*entryDiffGFC2PtRP)+=(Long_t)avMultDiffFlowGFCRP;
1813 entryDiffGFC2PtRP->Append(", N = ");
1814 (*entryDiffGFC2PtRP)+=(Long_t)nEvtsDiffFlowGFCRP;
1815 legendDiffFlowPtRP->AddEntry(gfcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffGFC2PtRP->Data(),"p");
1816 }
1817 if(plotGFC4PtRP && gfcCommonHistRes4)
1818 {
1819 entryDiffGFC4PtRP->Append("M = ");
1820 (*entryDiffGFC4PtRP)+=(Long_t)avMultDiffFlowGFCRP;
1821 entryDiffGFC4PtRP->Append(", N = ");
1822 (*entryDiffGFC4PtRP)+=(Long_t)nEvtsDiffFlowGFCRP;
1823 legendDiffFlowPtRP->AddEntry(gfcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffGFC4PtRP->Data(),"p");
1824 }
1825 if(plotGFC6PtRP && gfcCommonHistRes6)
1826 {
1827 entryDiffGFC6PtRP->Append("M = ");
1828 (*entryDiffGFC6PtRP)+=(Long_t)avMultDiffFlowGFCRP;
1829 entryDiffGFC6PtRP->Append(", N = ");
1830 (*entryDiffGFC6PtRP)+=(Long_t)nEvtsDiffFlowGFCRP;
1831 legendDiffFlowPtRP->AddEntry(gfcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffGFC6PtRP->Data(),"p");
1832 }
1833 if(plotGFC8PtRP && gfcCommonHistRes8)
1834 {
1835 entryDiffGFC8PtRP->Append("M = ");
1836 (*entryDiffGFC8PtRP)+=(Long_t)avMultDiffFlowGFCRP;
1837 entryDiffGFC8PtRP->Append(", N = ");
1838 (*entryDiffGFC8PtRP)+=(Long_t)nEvtsDiffFlowGFCRP;
1839 legendDiffFlowPtRP->AddEntry(gfcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffGFC8PtRP->Data(),"p");
1840 }
1841
1842 //QC
1843 if(plotQC2PtRP && qcCommonHistRes2)
1844 {
1845 entryDiffQC2PtRP->Append("M = ");
1846 (*entryDiffQC2PtRP)+=(Long_t)avMultDiffFlowQC2RP;
1847 entryDiffQC2PtRP->Append(", N = ");
1848 (*entryDiffQC2PtRP)+=(Long_t)nEvtsDiffFlowQC2RP;
1849 legendDiffFlowPtRP->AddEntry(qcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffQC2PtRP->Data(),"p");
1850 }
1851 if(plotQC4PtRP && qcCommonHistRes4)
1852 {
1853 entryDiffQC4PtRP->Append("M = ");
1854 (*entryDiffQC4PtRP)+=(Long_t)avMultDiffFlowQC4RP;
1855 entryDiffQC4PtRP->Append(", N = ");
1856 (*entryDiffQC4PtRP)+=(Long_t)nEvtsDiffFlowQC4RP;
1857 legendDiffFlowPtRP->AddEntry(qcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffQC4PtRP->Data(),"p");
1858 }
1859 if(plotQC6PtRP && qcCommonHistRes6)
1860 {
1861 entryDiffQC6PtRP->Append("M = ");
1862 (*entryDiffQC6PtRP)+=(Long_t)avMultDiffFlowQC6RP;
1863 entryDiffQC6PtRP->Append(", N = ");
1864 (*entryDiffQC6PtRP)+=(Long_t)nEvtsDiffFlowQC6RP;
1865 legendDiffFlowPtRP->AddEntry(qcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffQC6PtRP->Data(),"p");
1866 }
1867 if(plotQC8PtRP && qcCommonHistRes8)
1868 {
1869 entryDiffQC8PtRP->Append("M = ");
1870 (*entryDiffQC8PtRP)+=(Long_t)avMultDiffFlowQC8RP;
1871 entryDiffQC8PtRP->Append(", N = ");
1872 (*entryDiffQC8PtRP)+=(Long_t)nEvtsDiffFlowQC8RP;
1873 legendDiffFlowPtRP->AddEntry(qcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffQC8PtRP->Data(),"p");
1874 }
1875
1876 //LYZ2
1877 if(plotLYZ2PtRP && lyz2CommonHistRes)
1878 {
1879 entryDiffLYZ2PtRP->Append("M = ");
1880 (*entryDiffLYZ2PtRP)+=(Long_t)avMultDiffFlowLYZ2RP;
1881 entryDiffLYZ2PtRP->Append(", N = ");
1882 (*entryDiffLYZ2PtRP)+=(Long_t)nEvtsDiffFlowLYZ2RP;
1883 legendDiffFlowPtRP->AddEntry(lyz2CommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZ2PtRP->Data(),"p");
1884 }
1885
1886 //LYZEP
1887 if(plotLYZEPPtRP && lyzepCommonHistRes)
1888 {
1889 entryDiffLYZEPPtRP->Append("M = ");
1890 (*entryDiffLYZEPPtRP)+=(Long_t)avMultDiffFlowLYZEPRP;
1891 entryDiffLYZEPPtRP->Append(", N = ");
1892 (*entryDiffLYZEPPtRP)+=(Long_t)nEvtsDiffFlowLYZEPRP;
1893 legendDiffFlowPtRP->AddEntry(lyzepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZEPPtRP->Data(),"p");
1894 }
1fca9c90 1895
e83922f2 1896 //drawing finally the legend in the 2nd pad:
1897 if(legendDiffFlowPtRP)
1898 {
1899 legendDiffFlowPtRP->SetMargin(0.15);
1900 legendDiffFlowPtRP->Draw();
1901 }
1902 }// end of if(plotDiffFlowPtRP)
1fca9c90 1903 //----------------------------------------------------------------------------------
1904
1905
1906 //----------------------------------------------------------------------------------
1907 //final drawing for differential flow (Eta, RP):
e83922f2 1908 //set here the results of which methods will be plotted by default:
1909 Bool_t plotMCEtaRP = kFALSE;
1910 Bool_t plotGFC2EtaRP = kTRUE;
1911 Bool_t plotGFC4EtaRP = kTRUE;
1912 Bool_t plotGFC6EtaRP = kTRUE;
1913 Bool_t plotGFC8EtaRP = kTRUE;
1914 Bool_t plotQC2EtaRP = kTRUE;
1915 Bool_t plotQC4EtaRP = kTRUE;
1916 Bool_t plotQC6EtaRP = kTRUE;
1917 Bool_t plotQC8EtaRP = kTRUE;
1918 Bool_t plotLYZ2EtaRP = kTRUE;
1919 Bool_t plotLYZEPEtaRP = kTRUE;
1920 if(plotDiffFlowEtaRP)
1921 {
1922 TCanvas* diffFlowEtaAllCanvasRP = new TCanvas("Differential Flow (Eta) of RP","Differential Flow (Eta) of RP ",1000,600);
1923
1924 diffFlowEtaAllCanvasRP->Divide(2,1);
1925
1926 //1st pad is for plot:
1927 (diffFlowEtaAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0);
1928
1929 if(styleHistEta)
1930 {
1931 (styleHistEta->GetYaxis())->SetRangeUser(-0.3,1.0);
1932 styleHistEta->Draw();
1933 }
1934 if(pMeshDiffFlowEtaRP)
1935 {
1936 pMeshDiffFlowEtaRP->Draw("LFSAME");
1937 }
1fca9c90 1938
e83922f2 1939 //MC
1940 if(plotMCEtaRP && mcepCommonHistRes)
1941 {
1942 (mcepCommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1943 }
1944 //GFC
1945 if(plotGFC2EtaRP && gfcCommonHistRes2)
1946 {
1947 (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1948 }
1949 if(plotGFC4EtaRP && gfcCommonHistRes4)
1950 {
1951 (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1952 }
1953 if(plotGFC6EtaRP && gfcCommonHistRes6)
1954 {
1955 (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1956 }
1957 if(plotGFC8EtaRP && gfcCommonHistRes8)
1958 {
1959 (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1960 }
1961 //QC
1962 if(plotQC2EtaRP && qcCommonHistRes2)
1963 {
1964 (qcCommonHistRes2->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1965 }
1966 if(plotQC4EtaRP && qcCommonHistRes4)
1967 {
1968 (qcCommonHistRes4->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1969 }
1970 if(plotQC6EtaRP && qcCommonHistRes6)
1971 {
1972 (qcCommonHistRes6->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1973 }
1974 if(plotQC8EtaRP && qcCommonHistRes8)
1975 {
1976 (qcCommonHistRes8->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1977 }
1978 //LYZ2
1979 if(plotLYZ2EtaRP && lyz2CommonHistRes)
1980 {
1981 (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1982 }
1983 //LYZEP
1984 if(plotLYZEPEtaRP && lyzepCommonHistRes)
1985 {
1986 (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME");
1987 }
1fca9c90 1988
e83922f2 1989 //2nd pad is for legend:
1990 (diffFlowEtaAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0);
1991
1992 TLegend* legendDiffFlowEtaRP = new TLegend(0.02,0.25,0.97,0.75);
1993 legendDiffFlowEtaRP->SetTextFont(72);
1994 legendDiffFlowEtaRP->SetTextSize(0.06);
1995
1996 //legend's entries:
1997 TString *entryDiffMCEtaRP = new TString("MC ....... ");
1998 TString *entryDiffGFC2EtaRP = new TString("GFC{2} ... ");
1999 TString *entryDiffGFC4EtaRP = new TString("GFC{4} ... ");
2000 TString *entryDiffGFC6EtaRP = new TString("GFC{6} ... ");
2001 TString *entryDiffGFC8EtaRP = new TString("GFC{8} ... ");
2002 TString *entryDiffQC2EtaRP = new TString("QC{2} .... ");
2003 TString *entryDiffQC4EtaRP = new TString("QC{4} .... ");
2004 TString *entryDiffQC6EtaRP = new TString("QC{6} .... ");
2005 TString *entryDiffQC8EtaRP = new TString("QC{8} .... ");
2006 TString *entryDiffLYZ2EtaRP = new TString("LYZ ...... ");
2007 TString *entryDiffLYZEPEtaRP = new TString("LYZEP ... ");
2008
2009 //MC
2010 if(mcepCommonHistRes)
2011 {
2012 (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetFillStyle(meshStyle);
2013 (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetFillColor(meshColor);
2014 entryDiffMCEtaRP->Append("M = ");
2015 (*entryDiffMCEtaRP)+=(Long_t)avMultDiffFlowMCRP;
2016 entryDiffMCEtaRP->Append(", N = ");
2017 (*entryDiffMCEtaRP)+=(Long_t)nEvtsDiffFlowMCRP;
2018 legendDiffFlowEtaRP->AddEntry(mcepCommonHistRes->GetHistDiffFlowEtaRP(),entryDiffMCEtaRP->Data(),"f");
2019 }
1fca9c90 2020
e83922f2 2021 //GFC
2022 if(plotGFC2EtaRP && gfcCommonHistRes2)
2023 {
2024 entryDiffGFC2EtaRP->Append("M = ");
2025 (*entryDiffGFC2EtaRP)+=(Long_t)avMultDiffFlowGFCRP;
2026 entryDiffGFC2EtaRP->Append(", N = ");
2027 (*entryDiffGFC2EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP;
2028 legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes2->GetHistDiffFlowEtaRP(),entryDiffGFC2EtaRP->Data(),"p");
2029 }
2030 if(plotGFC4EtaRP && gfcCommonHistRes4)
2031 {
2032 entryDiffGFC4EtaRP->Append("M = ");
2033 (*entryDiffGFC4EtaRP)+=(Long_t)avMultDiffFlowGFCRP;
2034 entryDiffGFC4EtaRP->Append(", N = ");
2035 (*entryDiffGFC4EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP;
2036 legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes4->GetHistDiffFlowEtaRP(),entryDiffGFC4EtaRP->Data(),"p");
2037 }
2038 if(plotGFC6EtaRP && gfcCommonHistRes6)
2039 {
2040 entryDiffGFC6EtaRP->Append("M = ");
2041 (*entryDiffGFC6EtaRP)+=(Long_t)avMultDiffFlowGFCRP;
2042 entryDiffGFC6EtaRP->Append(", N = ");
2043 (*entryDiffGFC6EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP;
2044 legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes6->GetHistDiffFlowEtaRP(),entryDiffGFC6EtaRP->Data(),"p");
2045 }
2046 if(plotGFC8EtaRP && gfcCommonHistRes8)
2047 {
2048 entryDiffGFC8EtaRP->Append("M = ");
2049 (*entryDiffGFC8EtaRP)+=(Long_t)avMultDiffFlowGFCRP;
2050 entryDiffGFC8EtaRP->Append(", N = ");
2051 (*entryDiffGFC8EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP;
2052 legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes8->GetHistDiffFlowEtaRP(),entryDiffGFC8EtaRP->Data(),"p");
2053 }
2054
2055 //QC
2056 if(plotQC2EtaRP && qcCommonHistRes2)
2057 {
2058 entryDiffQC2EtaRP->Append("M = ");
2059 (*entryDiffQC2EtaRP)+=(Long_t)avMultDiffFlowQC2RP;
2060 entryDiffQC2EtaRP->Append(", N = ");
2061 (*entryDiffQC2EtaRP)+=(Long_t)nEvtsDiffFlowQC2RP;
2062 legendDiffFlowEtaRP->AddEntry(qcCommonHistRes2->GetHistDiffFlowEtaRP(),entryDiffQC2EtaRP->Data(),"p");
2063 }
2064 if(plotQC4EtaRP && qcCommonHistRes4)
2065 {
2066 entryDiffQC4EtaRP->Append("M = ");
2067 (*entryDiffQC4EtaRP)+=(Long_t)avMultDiffFlowQC4RP;
2068 entryDiffQC4EtaRP->Append(", N = ");
2069 (*entryDiffQC4EtaRP)+=(Long_t)nEvtsDiffFlowQC4RP;
2070 legendDiffFlowEtaRP->AddEntry(qcCommonHistRes4->GetHistDiffFlowEtaRP(),entryDiffQC4EtaRP->Data(),"p");
2071 }
2072 if(plotQC6EtaRP && qcCommonHistRes6)
2073 {
2074 entryDiffQC6EtaRP->Append("M = ");
2075 (*entryDiffQC6EtaRP)+=(Long_t)avMultDiffFlowQC6RP;
2076 entryDiffQC6EtaRP->Append(", N = ");
2077 (*entryDiffQC6EtaRP)+=(Long_t)nEvtsDiffFlowQC6RP;
2078 legendDiffFlowEtaRP->AddEntry(qcCommonHistRes6->GetHistDiffFlowEtaRP(),entryDiffQC6EtaRP->Data(),"p");
2079 }
2080 if(plotQC8EtaRP && qcCommonHistRes8)
2081 {
2082 entryDiffQC8EtaRP->Append("M = ");
2083 (*entryDiffQC8EtaRP)+=(Long_t)avMultDiffFlowQC8RP;
2084 entryDiffQC8EtaRP->Append(", N = ");
2085 (*entryDiffQC8EtaRP)+=(Long_t)nEvtsDiffFlowQC8RP;
2086 legendDiffFlowEtaRP->AddEntry(qcCommonHistRes8->GetHistDiffFlowEtaRP(),entryDiffQC8EtaRP->Data(),"p");
2087 }
1fca9c90 2088
e83922f2 2089 //LYZ2
2090 if(plotLYZ2EtaRP && lyz2CommonHistRes)
2091 {
2092 entryDiffLYZ2EtaRP->Append("M = ");
2093 (*entryDiffLYZ2EtaRP)+=(Long_t)avMultDiffFlowLYZ2RP;
2094 entryDiffLYZ2EtaRP->Append(", N = ");
2095 (*entryDiffLYZ2EtaRP)+=(Long_t)nEvtsDiffFlowLYZ2RP;
2096 legendDiffFlowEtaRP->AddEntry(lyz2CommonHistRes->GetHistDiffFlowEtaRP(),entryDiffLYZ2EtaRP->Data(),"p");
2097 }
2098
2099 //LYZEP
2100 if(plotLYZEPEtaRP && lyzepCommonHistRes)
2101 {
2102 entryDiffLYZEPEtaRP->Append("M = ");
2103 (*entryDiffLYZEPEtaRP)+=(Long_t)avMultDiffFlowLYZEPRP;
2104 entryDiffLYZEPEtaRP->Append(", N = ");
2105 (*entryDiffLYZEPEtaRP)+=(Long_t)nEvtsDiffFlowLYZEPRP;
2106 legendDiffFlowEtaRP->AddEntry(lyzepCommonHistRes->GetHistDiffFlowEtaRP(),entryDiffLYZEPEtaRP->Data(),"p");
2107 }
1fca9c90 2108
e83922f2 2109 //drawing finally the legend in the 2nd pad:
2110 if(legendDiffFlowEtaRP)
2111 {
2112 legendDiffFlowEtaRP->SetMargin(0.15);
2113 legendDiffFlowEtaRP->Draw();
2114 }
2115 }// end of if(plotDiffFlowEtaRP)
1fca9c90 2116 //----------------------------------------------------------------------------------
2117
2118 //----------------------------------------------------------------------------------
2119 //final drawing for differential flow (Pt, POI):
e83922f2 2120 //set here the results of which methods will be plotted by default:
2121 Bool_t plotMCPtPOI = kFALSE;
2122 Bool_t plotGFC2PtPOI = kTRUE;
2123 Bool_t plotGFC4PtPOI = kTRUE;
2124 Bool_t plotGFC6PtPOI = kTRUE;
2125 Bool_t plotGFC8PtPOI = kTRUE;
2126 Bool_t plotQC2PtPOI = kTRUE;
2127 Bool_t plotQC4PtPOI = kTRUE;
2128 Bool_t plotQC6PtPOI = kTRUE;
2129 Bool_t plotQC8PtPOI = kTRUE;
2130 Bool_t plotLYZ2PtPOI = kTRUE;
2131 Bool_t plotLYZEPPtPOI = kTRUE;
2132 if(plotDiffFlowPtPOI)
2133 {
2134 TCanvas* diffFlowPtAllCanvasPOI = new TCanvas("Differential Flow (Pt) of POI","Differential Flow (Pt) of POI ",1000,600);
2135
2136 diffFlowPtAllCanvasPOI->Divide(2,1);
2137
2138 //1st pad is for plot:
2139 (diffFlowPtAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0);
2140
2141 if(styleHistPt)
2142 {
2143 (styleHistPt->GetYaxis())->SetRangeUser(-0.3,1.0);
2144 styleHistPt->Draw();
2145 }
2146 if(pMeshDiffFlowPtPOI)
2147 {
2148 pMeshDiffFlowPtPOI->Draw("LFSAME");
2149 }
1fca9c90 2150
e83922f2 2151 //MC
2152 if(plotMCPtPOI && mcepCommonHistRes)
2153 {
2154 (mcepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2155 }
2156 //GFC
2157 if(plotGFC2PtPOI && gfcCommonHistRes2)
2158 {
2159 (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2160 }
2161 if(plotGFC4PtPOI && gfcCommonHistRes4)
2162 {
2163 (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2164 }
2165 if(plotGFC6PtPOI && gfcCommonHistRes6)
2166 {
2167 (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2168 }
2169 if(plotGFC8PtPOI && gfcCommonHistRes8)
2170 {
2171 (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2172 }
2173 //QC
2174 if(plotQC2PtPOI && qcCommonHistRes2)
2175 {
2176 (qcCommonHistRes2->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2177 }
2178 if(plotQC4PtPOI && qcCommonHistRes4)
2179 {
2180 (qcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2181 }
2182 if(plotQC6PtPOI && qcCommonHistRes6)
2183 {
2184 (qcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2185 }
2186 if(plotQC8PtPOI && qcCommonHistRes8)
2187 {
2188 (qcCommonHistRes8->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2189 }
2190 //LYZ2
2191 if(plotLYZ2PtPOI && lyz2CommonHistRes)
2192 {
2193 (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2194 }
2195 //LYZEP
2196 if(plotLYZEPPtPOI && lyzepCommonHistRes)
2197 {
2198 (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME");
2199 }
1fca9c90 2200
e83922f2 2201 //2nd pad is for legend:
2202 (diffFlowPtAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0);
2203
2204 TLegend* legendDiffFlowPtPOI = new TLegend(0.02,0.25,0.97,0.75);
2205 legendDiffFlowPtPOI->SetTextFont(72);
2206 legendDiffFlowPtPOI->SetTextSize(0.06);
2207
2208 //legend's entries:
2209 TString *entryDiffMCPtPOI = new TString("MC ....... ");
2210 TString *entryDiffGFC2PtPOI = new TString("GFC{2} ... ");
2211 TString *entryDiffGFC4PtPOI = new TString("GFC{4} ... ");
2212 TString *entryDiffGFC6PtPOI = new TString("GFC{6} ... ");
2213 TString *entryDiffGFC8PtPOI = new TString("GFC{8} ... ");
2214 TString *entryDiffQC2PtPOI = new TString("QC{2} .... ");
2215 TString *entryDiffQC4PtPOI = new TString("QC{4} .... ");
2216 TString *entryDiffQC6PtPOI = new TString("QC{6} .... ");
2217 TString *entryDiffQC8PtPOI = new TString("QC{8} .... ");
2218 TString *entryDiffLYZ2PtPOI = new TString("LYZ ...... ");
2219 TString *entryDiffLYZEPPtPOI = new TString("LYZEP ... ");
2220
2221 //MC
2222 if(mcepCommonHistRes)
2223 {
2224 (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetFillStyle(meshStyle);
2225 (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetFillColor(meshColor);
2226 entryDiffMCPtPOI->Append("M = ");
2227 (*entryDiffMCPtPOI)+=(Long_t)avMultDiffFlowMCPOI;
2228 entryDiffMCPtPOI->Append(", N = ");
2229 (*entryDiffMCPtPOI)+=(Long_t)nEvtsDiffFlowMCPOI;
2230 legendDiffFlowPtPOI->AddEntry(mcepCommonHistRes->GetHistDiffFlowPtPOI(),entryDiffMCPtPOI->Data(),"f");
2231 }
1fca9c90 2232
e83922f2 2233 //GFC
2234 if(plotGFC2PtPOI && gfcCommonHistRes2)
2235 {
2236 entryDiffGFC2PtPOI->Append("M = ");
2237 (*entryDiffGFC2PtPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2238 entryDiffGFC2PtPOI->Append(", N = ");
2239 (*entryDiffGFC2PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2240 legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes2->GetHistDiffFlowPtPOI(),entryDiffGFC2PtPOI->Data(),"p");
2241 }
2242 if(plotGFC4PtPOI && gfcCommonHistRes4)
2243 {
2244 entryDiffGFC4PtPOI->Append("M = ");
2245 (*entryDiffGFC4PtPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2246 entryDiffGFC4PtPOI->Append(", N = ");
2247 (*entryDiffGFC4PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2248 legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes4->GetHistDiffFlowPtPOI(),entryDiffGFC4PtPOI->Data(),"p");
2249 }
2250 if(plotGFC6PtPOI && gfcCommonHistRes6)
2251 {
2252 entryDiffGFC6PtPOI->Append("M = ");
2253 (*entryDiffGFC6PtPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2254 entryDiffGFC6PtPOI->Append(", N = ");
2255 (*entryDiffGFC6PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2256 legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes6->GetHistDiffFlowPtPOI(),entryDiffGFC6PtPOI->Data(),"p");
2257 }
2258 if(plotGFC8PtPOI && gfcCommonHistRes8)
2259 {
2260 entryDiffGFC8PtPOI->Append("M = ");
2261 (*entryDiffGFC8PtPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2262 entryDiffGFC8PtPOI->Append(", N = ");
2263 (*entryDiffGFC8PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2264 legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes8->GetHistDiffFlowPtPOI(),entryDiffGFC8PtPOI->Data(),"p");
2265 }
2266
2267 //QC
2268 if(plotQC2PtPOI && qcCommonHistRes2)
2269 {
2270 entryDiffQC2PtPOI->Append("M = ");
2271 (*entryDiffQC2PtPOI)+=(Long_t)avMultDiffFlowQC2POI;
2272 entryDiffQC2PtPOI->Append(", N = ");
2273 (*entryDiffQC2PtPOI)+=(Long_t)nEvtsDiffFlowQC2POI;
2274 legendDiffFlowPtPOI->AddEntry(qcCommonHistRes2->GetHistDiffFlowPtPOI(),entryDiffQC2PtPOI->Data(),"p");
2275 }
2276 if(plotQC4PtPOI && qcCommonHistRes4)
2277 {
2278 entryDiffQC4PtPOI->Append("M = ");
2279 (*entryDiffQC4PtPOI)+=(Long_t)avMultDiffFlowQC4POI;
2280 entryDiffQC4PtPOI->Append(", N = ");
2281 (*entryDiffQC4PtPOI)+=(Long_t)nEvtsDiffFlowQC4POI;
2282 legendDiffFlowPtPOI->AddEntry(qcCommonHistRes4->GetHistDiffFlowPtPOI(),entryDiffQC4PtPOI->Data(),"p");
2283 }
2284 if(plotQC6PtPOI && qcCommonHistRes6)
2285 {
2286 entryDiffQC6PtPOI->Append("M = ");
2287 (*entryDiffQC6PtPOI)+=(Long_t)avMultDiffFlowQC6POI;
2288 entryDiffQC6PtPOI->Append(", N = ");
2289 (*entryDiffQC6PtPOI)+=(Long_t)nEvtsDiffFlowQC6POI;
2290 legendDiffFlowPtPOI->AddEntry(qcCommonHistRes6->GetHistDiffFlowPtPOI(),entryDiffQC6PtPOI->Data(),"p");
2291 }
2292 if(plotQC8PtPOI && qcCommonHistRes8)
2293 {
2294 entryDiffQC8PtPOI->Append("M = ");
2295 (*entryDiffQC8PtPOI)+=(Long_t)avMultDiffFlowQC8POI;
2296 entryDiffQC8PtPOI->Append(", N = ");
2297 (*entryDiffQC8PtPOI)+=(Long_t)nEvtsDiffFlowQC8POI;
2298 legendDiffFlowPtPOI->AddEntry(qcCommonHistRes8->GetHistDiffFlowPtPOI(),entryDiffQC8PtPOI->Data(),"p");
2299 }
1fca9c90 2300
e83922f2 2301 //LYZ2
2302 if(plotLYZ2PtPOI && lyz2CommonHistRes)
2303 {
2304 entryDiffLYZ2PtPOI->Append("M = ");
2305 (*entryDiffLYZ2PtPOI)+=(Long_t)avMultDiffFlowLYZ2POI;
2306 entryDiffLYZ2PtPOI->Append(", N = ");
2307 (*entryDiffLYZ2PtPOI)+=(Long_t)nEvtsDiffFlowLYZ2POI;
2308 legendDiffFlowPtPOI->AddEntry(lyz2CommonHistRes->GetHistDiffFlowPtPOI(),entryDiffLYZ2PtPOI->Data(),"p");
2309 }
2310
2311 //LYZEP
2312 if(plotLYZEPPtPOI && lyzepCommonHistRes)
2313 {
2314 entryDiffLYZEPPtPOI->Append("M = ");
2315 (*entryDiffLYZEPPtPOI)+=(Long_t)avMultDiffFlowLYZEPPOI;
2316 entryDiffLYZEPPtPOI->Append(", N = ");
2317 (*entryDiffLYZEPPtPOI)+=(Long_t)nEvtsDiffFlowLYZEPPOI;
2318 legendDiffFlowPtPOI->AddEntry(lyzepCommonHistRes->GetHistDiffFlowPtPOI(),entryDiffLYZEPPtPOI->Data(),"p");
2319 }
1fca9c90 2320
e83922f2 2321 //drawing finally the legend in the 2nd pad:
2322 if(legendDiffFlowPtPOI)
2323 {
2324 legendDiffFlowPtPOI->SetMargin(0.15);
2325 legendDiffFlowPtPOI->Draw();
2326 }
2327 }//end of if(plotDiffFlowPtPOI)
1fca9c90 2328 //----------------------------------------------------------------------------------
2329
2330
2331 //----------------------------------------------------------------------------------
2332 //final drawing for differential flow (Eta, POI):
e83922f2 2333 //set here the results of which methods will be plotted by default:
2334 Bool_t plotMCEtaPOI = kFALSE;
2335 Bool_t plotGFC2EtaPOI = kTRUE;
2336 Bool_t plotGFC4EtaPOI = kTRUE;
2337 Bool_t plotGFC6EtaPOI = kTRUE;
2338 Bool_t plotGFC8EtaPOI = kTRUE;
2339 Bool_t plotQC2EtaPOI = kTRUE;
2340 Bool_t plotQC4EtaPOI = kTRUE;
2341 Bool_t plotQC6EtaPOI = kTRUE;
2342 Bool_t plotQC8EtaPOI = kTRUE;
2343 Bool_t plotLYZ2EtaPOI = kTRUE;
2344 Bool_t plotLYZEPEtaPOI = kTRUE;
2345 if(plotDiffFlowEtaPOI)
2346 {
2347 TCanvas* diffFlowEtaAllCanvasPOI = new TCanvas("Differential Flow (Eta) of POI","Differential Flow (Eta) of POI ",1000,600);
2348
2349 diffFlowEtaAllCanvasPOI->Divide(2,1);
2350
2351 //1st pad is for plot:
2352 (diffFlowEtaAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0);
1fca9c90 2353
e83922f2 2354 if(styleHistEta)
2355 {
2356 (styleHistEta->GetYaxis())->SetRangeUser(-0.3,1.0);
2357 styleHistEta->Draw();
2358 }
2359 if(pMeshDiffFlowEtaPOI)
2360 {
2361 pMeshDiffFlowEtaPOI->Draw("LFSAME");
2362 }
1fca9c90 2363
e83922f2 2364 //MC
2365 if(plotMCEtaPOI && mcepCommonHistRes)
2366 {
2367 (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2368 }
2369 //GFC
2370 if(plotGFC2EtaPOI && gfcCommonHistRes2)
2371 {
2372 (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2373 }
2374 if(plotGFC4EtaPOI && gfcCommonHistRes4)
2375 {
2376 (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2377 }
2378 if(plotGFC6EtaPOI && gfcCommonHistRes6)
2379 {
2380 (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2381 }
2382 if(plotGFC8EtaPOI && gfcCommonHistRes8)
2383 {
2384 (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2385 }
2386 //QC
2387 if(plotQC2EtaPOI && qcCommonHistRes2)
2388 {
2389 (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2390 }
2391 if(plotQC4EtaPOI && qcCommonHistRes4)
2392 {
2393 (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2394 }
2395 if(plotQC6EtaPOI && qcCommonHistRes6)
2396 {
2397 (qcCommonHistRes6->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2398 }
2399 if(plotQC8EtaPOI && qcCommonHistRes8)
2400 {
2401 (qcCommonHistRes8->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2402 }
2403 //LYZ2
2404 if(plotLYZ2EtaPOI && lyz2CommonHistRes)
2405 {
2406 (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2407 }
2408 //LYZEP
2409 if(plotLYZEPEtaPOI && lyzepCommonHistRes)
2410 {
2411 (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME");
2412 }
1fca9c90 2413
e83922f2 2414 //2nd pad is for legend:
2415 (diffFlowEtaAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0);
2416
2417 TLegend* legendDiffFlowEtaPOI = new TLegend(0.02,0.25,0.97,0.75);
2418 legendDiffFlowEtaPOI->SetTextFont(72);
2419 legendDiffFlowEtaPOI->SetTextSize(0.06);
2420
2421 //legend's entries:
2422 TString *entryDiffMCEtaPOI = new TString("MC ....... ");
2423 TString *entryDiffGFC2EtaPOI = new TString("GFC{2} ... ");
2424 TString *entryDiffGFC4EtaPOI = new TString("GFC{4} ... ");
2425 TString *entryDiffGFC6EtaPOI = new TString("GFC{6} ... ");
2426 TString *entryDiffGFC8EtaPOI = new TString("GFC{8} ... ");
2427 TString *entryDiffQC2EtaPOI = new TString("QC{2} .... ");
2428 TString *entryDiffQC4EtaPOI = new TString("QC{4} .... ");
2429 TString *entryDiffQC6EtaPOI = new TString("QC{6} .... ");
2430 TString *entryDiffQC8EtaPOI = new TString("QC{8} .... ");
2431 TString *entryDiffLYZ2EtaPOI = new TString("LYZ ...... ");
2432 TString *entryDiffLYZEPEtaPOI = new TString("LYZEP ... ");
2433
2434 //MC
2435 if(mcepCommonHistRes)
2436 {
2437 (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetFillStyle(meshStyle);
2438 (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetFillColor(meshColor);
2439 entryDiffMCEtaPOI->Append("M = ");
2440 (*entryDiffMCEtaPOI)+=(Long_t)avMultDiffFlowMCPOI;
2441 entryDiffMCEtaPOI->Append(", N = ");
2442 (*entryDiffMCEtaPOI)+=(Long_t)nEvtsDiffFlowMCPOI;
2443 legendDiffFlowEtaPOI->AddEntry(mcepCommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffMCEtaPOI->Data(),"f");
2444 }
1fca9c90 2445
e83922f2 2446 //GFC
2447 if(plotGFC2EtaPOI && gfcCommonHistRes2)
2448 {
2449 entryDiffGFC2EtaPOI->Append("M = ");
2450 (*entryDiffGFC2EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2451 entryDiffGFC2EtaPOI->Append(", N = ");
2452 (*entryDiffGFC2EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2453 legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes2->GetHistDiffFlowEtaPOI(),entryDiffGFC2EtaPOI->Data(),"p");
2454 }
2455 if(plotGFC4EtaPOI && gfcCommonHistRes4)
2456 {
2457 entryDiffGFC4EtaPOI->Append("M = ");
2458 (*entryDiffGFC4EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2459 entryDiffGFC4EtaPOI->Append(", N = ");
2460 (*entryDiffGFC4EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2461 legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes4->GetHistDiffFlowEtaPOI(),entryDiffGFC4EtaPOI->Data(),"p");
2462 }
2463 if(plotGFC6EtaPOI && gfcCommonHistRes6)
2464 {
2465 entryDiffGFC6EtaPOI->Append("M = ");
2466 (*entryDiffGFC6EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2467 entryDiffGFC6EtaPOI->Append(", N = ");
2468 (*entryDiffGFC6EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2469 legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes6->GetHistDiffFlowEtaPOI(),entryDiffGFC6EtaPOI->Data(),"p");
2470 }
2471 if(plotGFC8EtaPOI && gfcCommonHistRes8)
2472 {
2473 entryDiffGFC8EtaPOI->Append("M = ");
2474 (*entryDiffGFC8EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI;
2475 entryDiffGFC8EtaPOI->Append(", N = ");
2476 (*entryDiffGFC8EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI;
2477 legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes8->GetHistDiffFlowEtaPOI(),entryDiffGFC8EtaPOI->Data(),"p");
2478 }
2479
2480 //QC
2481 if(plotQC2EtaPOI && qcCommonHistRes2)
2482 {
2483 entryDiffQC2EtaPOI->Append("M = ");
2484 (*entryDiffQC2EtaPOI)+=(Long_t)avMultDiffFlowQC2POI;
2485 entryDiffQC2EtaPOI->Append(", N = ");
2486 (*entryDiffQC2EtaPOI)+=(Long_t)nEvtsDiffFlowQC2POI;
2487 legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes2->GetHistDiffFlowEtaPOI(),entryDiffQC2EtaPOI->Data(),"p");
2488 }
2489 if(plotQC4EtaPOI && qcCommonHistRes4)
2490 {
2491 entryDiffQC4EtaPOI->Append("M = ");
2492 (*entryDiffQC4EtaPOI)+=(Long_t)avMultDiffFlowQC4POI;
2493 entryDiffQC4EtaPOI->Append(", N = ");
2494 (*entryDiffQC4EtaPOI)+=(Long_t)nEvtsDiffFlowQC4POI;
2495 legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes4->GetHistDiffFlowEtaPOI(),entryDiffQC4EtaPOI->Data(),"p");
2496 }
2497 if(plotQC6EtaPOI && qcCommonHistRes6)
2498 {
2499 entryDiffQC6EtaPOI->Append("M = ");
2500 (*entryDiffQC6EtaPOI)+=(Long_t)avMultDiffFlowQC6POI;
2501 entryDiffQC6EtaPOI->Append(", N = ");
2502 (*entryDiffQC6EtaPOI)+=(Long_t)nEvtsDiffFlowQC6POI;
2503 legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes6->GetHistDiffFlowEtaPOI(),entryDiffQC6EtaPOI->Data(),"p");
2504 }
2505 if(plotQC8EtaPOI && qcCommonHistRes8)
2506 {
2507 entryDiffQC8EtaPOI->Append("M = ");
2508 (*entryDiffQC8EtaPOI)+=(Long_t)avMultDiffFlowQC8POI;
2509 entryDiffQC8EtaPOI->Append(", N = ");
2510 (*entryDiffQC8EtaPOI)+=(Long_t)nEvtsDiffFlowQC8POI;
2511 legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes8->GetHistDiffFlowEtaPOI(),entryDiffQC8EtaPOI->Data(),"p");
2512 }
1fca9c90 2513
e83922f2 2514 //LYZ2
2515 if(plotLYZ2EtaPOI && lyz2CommonHistRes)
2516 {
2517 entryDiffLYZ2EtaPOI->Append("M = ");
2518 (*entryDiffLYZ2EtaPOI)+=(Long_t)avMultDiffFlowLYZ2POI;
2519 entryDiffLYZ2EtaPOI->Append(", N = ");
2520 (*entryDiffLYZ2EtaPOI)+=(Long_t)nEvtsDiffFlowLYZ2POI;
2521 legendDiffFlowEtaPOI->AddEntry(lyz2CommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffLYZ2EtaPOI->Data(),"p");
2522 }
2523
2524 //LYZEP
2525 if(plotLYZEPEtaPOI && lyzepCommonHistRes)
2526 {
2527 entryDiffLYZEPEtaPOI->Append("M = ");
2528 (*entryDiffLYZEPEtaPOI)+=(Long_t)avMultDiffFlowLYZEPPOI;
2529 entryDiffLYZEPEtaPOI->Append(", N = ");
2530 (*entryDiffLYZEPEtaPOI)+=(Long_t)nEvtsDiffFlowLYZEPPOI;
2531 legendDiffFlowEtaPOI->AddEntry(lyzepCommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffLYZEPEtaPOI->Data(),"p");
2532 }
1fca9c90 2533
e83922f2 2534 //drawing finally the legend in the 2nd pad:
2535 if(legendDiffFlowEtaPOI)
2536 {
2537 legendDiffFlowEtaPOI->SetMargin(0.15);
2538 legendDiffFlowEtaPOI->Draw();
2539 }
2540 }//end of if(plotDiffFlowEtaPOI)
1fca9c90 2541 //----------------------------------------------------------------------------------
2542
2543
a88e5ac2 2544 //=====================================================================================
2545
70288839 2546}