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