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