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