]>
Commit | Line | Data |
---|---|---|
8c5ec46a | 1 | //type of analysis can be: ESD, AOD, MC, ESDMC0, ESDMC1 |
2 | //const TString type = "ESD"; | |
b25cc698 | 3 | |
4 | enum libModes {mLocal,mLocalSource}; | |
5 | //mLocal: Analyze data on your computer using aliroot | |
6 | //mLocalSource: Analyze data on your computer using root + source files | |
7 | ||
8 | //void compareFlowResults(TString type="",Int_t mode=mLocalSource) | |
9 | void compareFlowResults(TString type="ESD",Int_t mode=mLocal) | |
da24213b | 10 | { |
b25cc698 | 11 | |
12 | // load needed libraries: | |
8d5cd720 | 13 | LoadPlotLibraries(mode); |
b25cc698 | 14 | |
15 | ||
da24213b | 16 | //================================================================================== |
17 | // set here which plots will be shown by default | |
18 | //================================================================================== | |
688877b8 | 19 | Bool_t plotIntFlow = kTRUE; // integrated flow (no-name) // to be improved |
da24213b | 20 | // RP = particles used to determine the reaction plane |
688877b8 | 21 | Bool_t plotIntFlowRP = kTRUE; // integrated flow RP |
22 | Bool_t plotDiffFlowPtRP = kTRUE; // differential flow (Pt,RP) | |
23 | Bool_t plotDiffFlowEtaRP = kTRUE; // differential flow (Eta,RP) | |
24 | Bool_t plotDiffFlowPtRelativeToMCRP = kTRUE; // plot |v{MC}-v{method}/v{MC}| as a function of pt for RPs | |
da24213b | 25 | // POI = particle of interest |
688877b8 | 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) | |
da24213b | 29 | //================================================================================== |
30 | ||
b25cc698 | 31 | |
da24213b | 32 | //================================================================================== |
33 | // set here which methods will be plotted by default for differential flow (Pt,RP): | |
34 | Bool_t plotMCPtRP = kFALSE; | |
905bed86 | 35 | Bool_t plotSPPtRP = kTRUE; |
da24213b | 36 | Bool_t plotGFC2PtRP = kTRUE; |
37 | Bool_t plotGFC4PtRP = kTRUE; | |
38 | Bool_t plotGFC6PtRP = kTRUE; | |
39 | Bool_t plotGFC8PtRP = kTRUE; | |
40 | Bool_t plotQC2PtRP = kTRUE; | |
41 | Bool_t plotQC4PtRP = kTRUE; | |
688877b8 | 42 | Bool_t plotQC6PtRP = kFALSE; // not calculated yet |
43 | Bool_t plotQC8PtRP = kFALSE; // not calculated yet | |
da24213b | 44 | Bool_t plotLYZ2PtRP = kTRUE; |
45 | ||
46 | // set here which methods will be plotted by default for differential flow (Eta,RP): | |
47 | Bool_t plotMCEtaRP = kFALSE; | |
905bed86 | 48 | Bool_t plotSPEtaRP = kTRUE; |
da24213b | 49 | Bool_t plotGFC2EtaRP = kTRUE; |
50 | Bool_t plotGFC4EtaRP = kTRUE; | |
51 | Bool_t plotGFC6EtaRP = kTRUE; | |
52 | Bool_t plotGFC8EtaRP = kTRUE; | |
53 | Bool_t plotQC2EtaRP = kTRUE; | |
54 | Bool_t plotQC4EtaRP = kTRUE; | |
688877b8 | 55 | Bool_t plotQC6EtaRP = kFALSE; // not calculated yet |
56 | Bool_t plotQC8EtaRP = kFALSE; // not calculated yet | |
da24213b | 57 | Bool_t plotLYZ2EtaRP = kTRUE; |
58 | Bool_t plotLYZEPEtaRP = kTRUE; | |
59 | ||
688877b8 | 60 | // set here which methods will be plotted by default for |v{MC}-v{method}/v{MC}| as a function of pt for RPs |
61 | Bool_t plotSPRelativeToMCRP = kTRUE; | |
62 | Bool_t plotGFC2RelativeToMCRP = kTRUE; | |
63 | Bool_t plotGFC4RelativeToMCRP = kTRUE; | |
64 | Bool_t plotGFC6RelativeToMCRP = kTRUE; | |
65 | Bool_t plotGFC8RelativeToMCRP = kTRUE; | |
66 | Bool_t plotQC2RelativeToMCRP = kTRUE; | |
67 | Bool_t plotQC4RelativeToMCRP = kTRUE; | |
68 | Bool_t plotQC6RelativeToMCRP = kFALSE; // not calculated yet | |
69 | Bool_t plotQC8RelativeToMCRP = kFALSE; // not calculated yet | |
70 | Bool_t plotLYZ2RelativeToMCRP = kTRUE; | |
71 | Bool_t plotLYZEPRelativeToMCRP = kTRUE; | |
72 | ||
da24213b | 73 | // set here which methods will be plotted by default for differential flow (Pt,POI): |
74 | Bool_t plotMCPtPOI = kFALSE; | |
905bed86 | 75 | Bool_t plotSPPtPOI = kTRUE; |
da24213b | 76 | Bool_t plotGFC2PtPOI = kTRUE; |
77 | Bool_t plotGFC4PtPOI = kTRUE; | |
78 | Bool_t plotGFC6PtPOI = kTRUE; | |
79 | Bool_t plotGFC8PtPOI = kTRUE; | |
80 | Bool_t plotQC2PtPOI = kTRUE; | |
81 | Bool_t plotQC4PtPOI = kTRUE; | |
688877b8 | 82 | Bool_t plotQC6PtPOI = kFALSE; |
83 | Bool_t plotQC8PtPOI = kFALSE; | |
da24213b | 84 | Bool_t plotLYZ2PtPOI = kTRUE; |
85 | Bool_t plotLYZEPPtPOI = kTRUE; | |
86 | ||
87 | // set here which methods will be plotted by default for differential flow (Eta,POI): | |
88 | Bool_t plotMCEtaPOI = kFALSE; | |
905bed86 | 89 | Bool_t plotSPEtaPOI = kTRUE; |
da24213b | 90 | Bool_t plotGFC2EtaPOI = kTRUE; |
91 | Bool_t plotGFC4EtaPOI = kTRUE; | |
92 | Bool_t plotGFC6EtaPOI = kTRUE; | |
93 | Bool_t plotGFC8EtaPOI = kTRUE; | |
94 | Bool_t plotQC2EtaPOI = kTRUE; | |
95 | Bool_t plotQC4EtaPOI = kTRUE; | |
688877b8 | 96 | Bool_t plotQC6EtaPOI = kFALSE; // not calculated yet |
97 | Bool_t plotQC8EtaPOI = kFALSE; // not calculated yet | |
da24213b | 98 | Bool_t plotLYZ2EtaPOI = kTRUE; |
99 | Bool_t plotLYZEPEtaPOI = kTRUE; | |
100 | //================================================================================== | |
101 | ||
102 | ||
103 | //================================================================================== | |
905bed86 | 104 | // cosmetics: marker style (see TMarker) and color (see TAttFill) for each method: |
da24213b | 105 | // MC: |
905bed86 | 106 | Int_t markerStyleMC = 20; // full circle |
da24213b | 107 | Int_t markerColorMC = kRed; |
905bed86 | 108 | // SP: |
109 | Int_t markerStyleSP = 3; // star | |
110 | Int_t markerColorSP = kViolet-6; | |
da24213b | 111 | // GFC{2} |
905bed86 | 112 | Int_t markerStyleGFC2 = 21; // full square |
da24213b | 113 | Int_t markerColorGFC2 = kAzure-7; |
114 | // GFC{4} | |
905bed86 | 115 | Int_t markerStyleGFC4 = 20; // full circle |
da24213b | 116 | Int_t markerColorGFC4 = kAzure+3; |
117 | // GFC{6} | |
688877b8 | 118 | Int_t markerStyleGFC6 = 25; // open square |
da24213b | 119 | Int_t markerColorGFC6 = kAzure-7; |
120 | // GFC{8} | |
688877b8 | 121 | Int_t markerStyleGFC8 = 24; // open circle |
da24213b | 122 | Int_t markerColorGFC8 = kAzure+3; |
123 | // QC{2} | |
905bed86 | 124 | Int_t markerStyleQC2 = 21; // full square |
da24213b | 125 | Int_t markerColorQC2 = kOrange-7; |
126 | // QC{4} | |
905bed86 | 127 | Int_t markerStyleQC4 = 20; // full circle |
da24213b | 128 | Int_t markerColorQC4 = kOrange+3; |
129 | // QC{6} | |
688877b8 | 130 | Int_t markerStyleQC6 = 25; // open square |
da24213b | 131 | Int_t markerColorQC6 = kOrange-7; |
132 | // QC{8} | |
688877b8 | 133 | Int_t markerStyleQC8 = 24; // open circle |
da24213b | 134 | Int_t markerColorQC8 = kOrange+3; |
135 | // LYZ2 | |
905bed86 | 136 | Int_t markerStyleLYZ2 = 22; // full triangle |
da24213b | 137 | Int_t markerColorLYZ2 = kYellow+3; |
138 | // LYZEP | |
905bed86 | 139 | Int_t markerStyleLYZEP = 26; // open triangle |
da24213b | 140 | Int_t markerColorLYZEP = kYellow+3; |
141 | //================================================================================== | |
142 | ||
dec3fab6 | 143 | |
144 | //================================================================================== | |
145 | // set here which result goes in which bin in the plot for integrated flow (no-name) | |
146 | // MC: | |
147 | Int_t binMC = 1; | |
148 | // SP: | |
149 | Int_t binSP = 2; | |
150 | // GFC{2} | |
151 | Int_t binGFC2 = 3; | |
152 | // GFC{4} | |
153 | Int_t binGFC4 = 5; | |
154 | // GFC{6} | |
155 | Int_t binGFC6 = 7; | |
156 | // GFC{8} | |
157 | Int_t binGFC8 = 9; | |
158 | // QC{2} | |
159 | Int_t binQC2 = 4; | |
160 | // QC{4} | |
161 | Int_t binQC4 = 6; | |
162 | // QC{6} | |
163 | Int_t binQC6 = 8; | |
164 | // QC{8} | |
165 | Int_t binQC8 = 10; | |
166 | // FQD | |
167 | Int_t binFQD = 11; | |
168 | // LYZ1 | |
169 | Int_t binLYZ1 = 12; | |
170 | // LYZEP | |
171 | Int_t binLYZEP = 13; | |
172 | //================================================================================== | |
173 | ||
174 | ||
175 | //================================================================================== | |
176 | // set here which result goes in which bin in the plot for integrated flow (RP) | |
177 | // MC: | |
178 | Int_t binMCRP = 1; | |
179 | // SP: | |
180 | Int_t binSPRP = 2; | |
181 | // GFC{2} | |
182 | Int_t binGFC2RP = 3; | |
183 | // GFC{4} | |
184 | Int_t binGFC4RP = 5; | |
185 | // GFC{6} | |
186 | Int_t binGFC6RP = 7; | |
187 | // GFC{8} | |
188 | Int_t binGFC8RP = 9; | |
189 | // QC{2} | |
190 | Int_t binQC2RP = 4; | |
191 | // QC{4} | |
192 | Int_t binQC4RP = 6; | |
193 | // QC{6} | |
194 | Int_t binQC6RP = 8; | |
195 | // QC{8} | |
196 | Int_t binQC8RP = 10; | |
197 | // FQD | |
198 | Int_t binFQDRP = 11; | |
199 | // LYZ2 | |
200 | Int_t binLYZ2RP = 12; | |
201 | // LYZEP | |
202 | Int_t binLYZEPRP = 13; | |
203 | //================================================================================== | |
204 | ||
205 | ||
206 | //================================================================================== | |
207 | // set here which result goes in which bin in the plot for integrated flow (POI) | |
208 | // MC: | |
209 | Int_t binMCPOI = 1; | |
210 | // SP: | |
211 | Int_t binSPPOI = 2; | |
212 | // GFC{2} | |
213 | Int_t binGFC2POI = 3; | |
214 | // GFC{4} | |
215 | Int_t binGFC4POI = 5; | |
216 | // GFC{6} | |
217 | Int_t binGFC6POI = 7; | |
218 | // GFC{8} | |
219 | Int_t binGFC8POI = 9; | |
220 | // QC{2} | |
221 | Int_t binQC2POI = 4; | |
222 | // QC{4} | |
223 | Int_t binQC4POI = 6; | |
224 | // QC{6} | |
225 | Int_t binQC6POI = 8; | |
226 | // QC{8} | |
227 | Int_t binQC8POI = 10; | |
228 | // FQD | |
229 | Int_t binFQDPOI = 11; | |
230 | // LYZ2 | |
231 | Int_t binLYZ2POI = 12; | |
232 | // LYZEP | |
233 | Int_t binLYZEPPOI = 13; | |
234 | //================================================================================== | |
b25cc698 | 235 | |
236 | ||
8c5ec46a | 237 | //================================================================================== |
238 | // accessing output files | |
239 | //================================================================================== | |
8c5ec46a | 240 | //open the output files: |
241 | TString inputFileNameMCEP = "outputMCEPanalysis"; | |
242 | TFile* fileMCEP = NULL; | |
243 | fileMCEP = TFile::Open(((inputFileNameMCEP.Append(type)).Append(".root")).Data(), "READ"); | |
244 | ||
905bed86 | 245 | TString inputFileNameSP = "outputSPanalysis"; |
246 | TFile* fileSP = NULL; | |
247 | fileSP = TFile::Open(((inputFileNameSP.Append(type)).Append(".root")).Data(), "READ"); | |
248 | ||
8c5ec46a | 249 | TString inputFileNameLYZ1 = "outputLYZ1analysis"; |
250 | TFile* fileLYZ1 = NULL; | |
251 | fileLYZ1 = TFile::Open(((inputFileNameLYZ1.Append(type)).Append(".root")).Data(), "READ"); | |
252 | ||
253 | TString inputFileNameLYZ2 = "outputLYZ2analysis"; | |
254 | TFile* fileLYZ2 = NULL; | |
255 | fileLYZ2 = TFile::Open(((inputFileNameLYZ2.Append(type)).Append(".root")).Data(), "READ"); | |
256 | ||
257 | TString inputFileNameLYZEP = "outputLYZEPanalysis"; | |
258 | TFile* fileLYZEP = NULL; | |
259 | fileLYZEP = TFile::Open(((inputFileNameLYZEP.Append(type)).Append(".root")).Data(), "READ"); | |
260 | ||
261 | TString inputFileNameFQD = "outputFQDanalysis"; | |
262 | TFile* fileFQD = NULL; | |
263 | fileFQD = TFile::Open(((inputFileNameFQD.Append(type)).Append(".root")).Data(), "READ"); | |
264 | ||
265 | TString inputFileNameGFC = "outputGFCanalysis"; | |
266 | TFile* fileGFC = NULL; | |
267 | fileGFC = TFile::Open(((inputFileNameGFC.Append(type)).Append(".root")).Data(), "READ"); | |
268 | ||
269 | TString inputFileNameQC = "outputQCanalysis"; | |
270 | TFile* fileQC = NULL; | |
271 | fileQC = TFile::Open(((inputFileNameQC.Append(type)).Append(".root")).Data(), "READ"); | |
272 | //================================================================================== | |
1fca9c90 | 273 | |
a88e5ac2 | 274 | |
a88e5ac2 | 275 | |
9d3cbee7 | 276 | |
8c5ec46a | 277 | //================================================================================== |
278 | // cosmetics | |
279 | //================================================================================== | |
280 | //removing the title and stat. box from all histograms: | |
281 | gStyle->SetOptTitle(0); | |
282 | gStyle->SetOptStat(0); | |
283 | ||
dec3fab6 | 284 | // plot for 'no-name' integrated flow: |
285 | // choosing the style and color of mesh for MC error bands | |
8c5ec46a | 286 | Int_t meshStyle = 1001; |
dec3fab6 | 287 | Int_t meshColor = kGray; |
8c5ec46a | 288 | |
dec3fab6 | 289 | // marker style and color |
8c5ec46a | 290 | Int_t markerStyle = 21; |
dec3fab6 | 291 | Int_t markerColor = kBlack; |
292 | ||
293 | // plot for RP's integrated flow: | |
294 | // choosing the style and color of mesh for MC error bands | |
295 | Int_t meshStyleRP = 1001; | |
296 | Int_t meshColorRP = kRed-10; | |
297 | ||
298 | // marker style and color | |
299 | Int_t markerStyleRP = 21; | |
300 | Int_t markerColorRP = kRed-3; | |
301 | ||
302 | // plot for POI's integrated flow: | |
303 | // choosing the style and color of mesh for MC error bands | |
304 | Int_t meshStylePOI = 1001; | |
305 | Int_t meshColorPOI = kBlue-10; | |
306 | ||
307 | // marker style and color | |
308 | Int_t markerStylePOI = 21; | |
309 | Int_t markerColorPOI = kBlue-3; | |
310 | ||
311 | // choosing the style and color of mesh for MC error bands in the plots for diff. flow | |
312 | // plot for differential flow (Pt,RP): | |
313 | Int_t meshStyleDiffFlowPtRP = 1001; | |
314 | Int_t meshColorDiffFlowPtRP = kRed-10; | |
315 | ||
316 | // plot for differential flow (Eta,RP): | |
317 | Int_t meshStyleDiffFlowEtaRP = 1001; | |
318 | Int_t meshColorDiffFlowEtaRP = kRed-10; | |
319 | ||
320 | // plot for differential flow (Pt,POI): | |
321 | Int_t meshStyleDiffFlowPtPOI = 1001; | |
322 | Int_t meshColorDiffFlowPtPOI = kRed-10; | |
323 | ||
324 | // plot for differential flow (Eta,POI): | |
325 | Int_t meshStyleDiffFlowEtaPOI = 1001; | |
326 | Int_t meshColorDiffFlowEtaPOI = kRed-10; | |
8c5ec46a | 327 | //================================================================================== |
328 | ||
329 | ||
330 | ||
331 | ||
332 | //================================================================================== | |
dec3fab6 | 333 | // INTEGRATED FLOW (no-name, RP and POI) |
8c5ec46a | 334 | //================================================================================== |
dec3fab6 | 335 | // the number of different methods: |
905bed86 | 336 | const Int_t nMethods=13; |
8c5ec46a | 337 | |
dec3fab6 | 338 | // booking the histogram for the integrated flow results from all methods: |
8c5ec46a | 339 | TH1D* intFlowAll = new TH1D("intFlowAll","Integrated Flow",nMethods,0,nMethods); |
dec3fab6 | 340 | // intFlowAll->SetLabelSize(0.036,"X"); |
341 | // intFlowAll->SetLabelSize(0.036,"Y"); | |
8c5ec46a | 342 | intFlowAll->SetMarkerStyle(markerStyle); |
343 | intFlowAll->SetMarkerColor(markerColor); | |
dec3fab6 | 344 | (intFlowAll->GetXaxis())->SetBinLabel(binMC,"v_{2}{MC}"); |
345 | (intFlowAll->GetXaxis())->SetBinLabel(binSP,"v_{2}{SP}"); | |
346 | (intFlowAll->GetXaxis())->SetBinLabel(binGFC2,"v_{2}{2,GFC}"); | |
347 | (intFlowAll->GetXaxis())->SetBinLabel(binQC2,"v_{2}{2,QC}"); | |
348 | (intFlowAll->GetXaxis())->SetBinLabel(binGFC4,"v_{2}{4,GFC}"); | |
349 | (intFlowAll->GetXaxis())->SetBinLabel(binQC4,"v_{2}{4,QC}"); | |
350 | (intFlowAll->GetXaxis())->SetBinLabel(binGFC6,"v_{2}{6,GFC}"); | |
351 | (intFlowAll->GetXaxis())->SetBinLabel(binQC6,"v_{2}{6,QC}"); | |
352 | (intFlowAll->GetXaxis())->SetBinLabel(binGFC8,"v_{2}{8,GFC}"); | |
353 | (intFlowAll->GetXaxis())->SetBinLabel(binQC8,"v_{2}{8,QC}"); | |
354 | (intFlowAll->GetXaxis())->SetBinLabel(binFQD,"v_{2}{FQD}"); | |
355 | (intFlowAll->GetXaxis())->SetBinLabel(binLYZ1,"v_{2}{LYZ}"); | |
356 | (intFlowAll->GetXaxis())->SetBinLabel(binLYZEP,"v_{2}{LYZEP}"); | |
357 | ||
358 | // booking the graph to store flow values and errors from all methods: | |
359 | Double_t x[nMethods] = {0.}; | |
360 | for(Int_t i=0;i<nMethods;i++) | |
361 | { | |
362 | x[i]=i+0.5; | |
363 | } | |
8c5ec46a | 364 | Double_t xError[nMethods] = {0.}; |
3963f25c | 365 | Double_t flowValue[nMethods] = {0.}; |
366 | Double_t flowError[nMethods] = {0.}; | |
8c5ec46a | 367 | Double_t flowValueRP[nMethods] = {0.}; |
368 | Double_t flowErrorRP[nMethods] = {0.}; | |
369 | Double_t flowValuePOI[nMethods] = {0.}; | |
370 | Double_t flowErrorPOI[nMethods] = {0.}; | |
371 | ||
dec3fab6 | 372 | // accessing the results for integrated flow for each method: |
373 | // MCEP = Monte Carlo Event Plane | |
8c5ec46a | 374 | TList *pListMCEP = NULL; |
375 | AliFlowCommonHist *mcepCommonHist = NULL; | |
376 | AliFlowCommonHistResults *mcepCommonHistRes = NULL; | |
377 | if(fileMCEP) { | |
378 | fileMCEP->GetObject("cobjMCEP",pListMCEP); | |
379 | if(pListMCEP) { | |
dec3fab6 | 380 | mcepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListMCEP->FindObject("AliFlowCommonHistMCEP")); |
8c5ec46a | 381 | mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListMCEP->FindObject("AliFlowCommonHistResultsMCEP")); |
382 | if(mcepCommonHistRes) { | |
dec3fab6 | 383 | flowValue[binMC-1] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1); |
384 | flowError[binMC-1] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1); | |
385 | flowValueRP[binMCRP-1] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1); | |
386 | flowErrorRP[binMCRP-1] = (mcepCommonHistRes->GetHistIntFlowRP())->GetBinError(1); | |
387 | flowValuePOI[binMCPOI-1] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1); | |
388 | flowErrorPOI[binMCPOI-1] = (mcepCommonHistRes->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 389 | } |
390 | } | |
a88e5ac2 | 391 | } |
8c5ec46a | 392 | |
dec3fab6 | 393 | // SP = Scalar Product |
905bed86 | 394 | TList *pListSP = NULL; |
395 | AliFlowCommonHist *spCommonHist = NULL; | |
396 | AliFlowCommonHistResults *spCommonHistRes = NULL; | |
397 | if(fileSP) { | |
398 | fileSP->GetObject("cobjSP",pListSP); | |
399 | if(pListSP) { | |
dec3fab6 | 400 | spCommonHist = dynamic_cast<AliFlowCommonHist*> (pListSP->FindObject("AliFlowCommonHistSP")); |
905bed86 | 401 | spCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListSP->FindObject("AliFlowCommonHistResultsSP")); |
402 | if(spCommonHistRes) { | |
dec3fab6 | 403 | flowValue[binSP-1] = (spCommonHistRes->GetHistIntFlow())->GetBinContent(1); |
404 | flowError[binSP-1] = (spCommonHistRes->GetHistIntFlow())->GetBinError(1); | |
405 | flowValueRP[binSPRP-1] = (spCommonHistRes->GetHistIntFlowRP())->GetBinContent(1); | |
406 | flowErrorRP[binSPRP-1] = (spCommonHistRes->GetHistIntFlowRP())->GetBinError(1); | |
407 | flowValuePOI[binSPPOI-1] = (spCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1); | |
408 | flowErrorPOI[binSPPOI-1] = (spCommonHistRes->GetHistIntFlowPOI())->GetBinError(1); | |
905bed86 | 409 | } |
410 | } | |
411 | } | |
412 | ||
dec3fab6 | 413 | // LYZ1 = Lee-Yang Zeros (1st run) is used to get only 'no-name' integrated flow |
8c5ec46a | 414 | TList *pListLYZ1 = NULL; |
415 | AliFlowCommonHist *lyz1CommonHist = NULL; | |
416 | AliFlowCommonHistResults *lyz1CommonHistRes = NULL; | |
417 | if(fileLYZ1) { | |
418 | fileLYZ1->GetObject("cobjLYZ1",pListLYZ1); | |
419 | if(pListLYZ1) { | |
420 | lyz1CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ1->FindObject("AliFlowCommonHistLYZ1")); | |
421 | lyz1CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ1->FindObject("AliFlowCommonHistResultsLYZ1")); | |
422 | if(lyz1CommonHistRes) { | |
dec3fab6 | 423 | flowValue[binLYZ1-1] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinContent(1); |
424 | flowError[binLYZ1-1] = (lyz1CommonHistRes->GetHistIntFlow())->GetBinError(1); | |
8c5ec46a | 425 | } |
426 | } | |
a88e5ac2 | 427 | } |
8c5ec46a | 428 | |
dec3fab6 | 429 | // LYZ2 = Lee-Yang Zeros (2nd run) is used to get RP's and POI's integrated flow |
8c5ec46a | 430 | TList *pListLYZ2 = NULL; |
431 | AliFlowCommonHist *lyz2CommonHist = NULL; | |
432 | AliFlowCommonHistResults *lyz2CommonHistRes = NULL; | |
433 | if(fileLYZ2) { | |
434 | fileLYZ2->GetObject("cobjLYZ2",pListLYZ2); | |
435 | if(pListLYZ2) { | |
436 | lyz2CommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZ2->FindObject("AliFlowCommonHistLYZ2")); | |
437 | lyz2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZ2->FindObject("AliFlowCommonHistResultsLYZ2")); | |
80e93fef | 438 | if(lyz2CommonHistRes) { |
dec3fab6 | 439 | flowValueRP[binLYZ2RP-1] = (lyz2CommonHistRes->GetHistIntFlowRP())->GetBinContent(1); |
440 | flowErrorRP[binLYZ2RP-1] = (lyz2CommonHistRes->GetHistIntFlowRP())->GetBinError(1); | |
441 | flowValuePOI[binLYZ2POI-1] = (lyz2CommonHistRes->GetHistIntFlowPOI())->GetBinContent(1); | |
442 | flowErrorPOI[binLYZ2POI-1] = (lyz2CommonHistRes->GetHistIntFlowPOI())->GetBinError(1); | |
80e93fef | 443 | } |
8c5ec46a | 444 | } |
9d3cbee7 | 445 | } |
8c5ec46a | 446 | |
dec3fab6 | 447 | // LYZEP = Lee-Yang Zeros Event Plane |
8c5ec46a | 448 | TList *pListLYZEP = NULL; |
449 | AliFlowCommonHist *lyzepCommonHist = NULL; | |
450 | AliFlowCommonHistResults *lyzepCommonHistRes = NULL; | |
451 | if(fileLYZEP) { | |
452 | fileLYZEP->GetObject("cobjLYZEP",pListLYZEP); | |
453 | if(pListLYZEP) { | |
454 | lyzepCommonHist = dynamic_cast<AliFlowCommonHist*> (pListLYZEP->FindObject("AliFlowCommonHistLYZEP")); | |
455 | lyzepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListLYZEP->FindObject("AliFlowCommonHistResultsLYZEP")); | |
456 | if(lyzepCommonHistRes) { | |
dec3fab6 | 457 | flowValue[binLYZEP-1] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinContent(1); |
458 | flowError[binLYZEP-1] = (lyzepCommonHistRes->GetHistIntFlow())->GetBinError(1); | |
459 | flowValueRP[binLYZEPRP-1] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinContent(1); | |
460 | flowErrorRP[binLYZEPRP-1] = (lyzepCommonHistRes->GetHistIntFlowRP())->GetBinError(1); | |
461 | flowValuePOI[binLYZEPPOI-1] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1); | |
462 | flowErrorPOI[binLYZEPPOI-1] = (lyzepCommonHistRes->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 463 | } |
464 | } | |
9d3cbee7 | 465 | } |
9d3cbee7 | 466 | |
dec3fab6 | 467 | // FQD = Fitting q-distribution |
8c5ec46a | 468 | TList *pListFQD = NULL; |
469 | AliFlowCommonHist *fqdCommonHist = NULL; | |
470 | AliFlowCommonHistResults *fqdCommonHistRes = NULL; | |
471 | if(fileFQD) { | |
472 | fileFQD->GetObject("cobjFQD",pListFQD); | |
473 | if(pListFQD) { | |
474 | fqdCommonHist = dynamic_cast<AliFlowCommonHist*> (pListFQD->FindObject("AliFlowCommonHistFQD")); | |
475 | fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (pListFQD->FindObject("AliFlowCommonHistResultsFQD")); | |
476 | if(fqdCommonHistRes) { | |
dec3fab6 | 477 | flowValue[binFQD-1] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1); |
478 | flowError[binFQD-1] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1); | |
479 | flowValueRP[binFQDRP-1] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinContent(1); | |
480 | flowErrorRP[binFQDRP-1] = (fqdCommonHistRes->GetHistIntFlowRP())->GetBinError(1); | |
481 | flowValuePOI[binFQDPOI-1] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinContent(1); | |
482 | flowErrorPOI[binFQDPOI-1] = (fqdCommonHistRes->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 483 | } |
484 | } | |
9d3cbee7 | 485 | } |
a88e5ac2 | 486 | |
dec3fab6 | 487 | // GFC = Generating Function Cumulants |
8c5ec46a | 488 | TList *pListGFC = NULL; |
489 | AliFlowCommonHist *gfcCommonHist = NULL; | |
490 | AliFlowCommonHistResults *gfcCommonHistRes2 = NULL; | |
491 | AliFlowCommonHistResults *gfcCommonHistRes4 = NULL; | |
492 | AliFlowCommonHistResults *gfcCommonHistRes6 = NULL; | |
493 | AliFlowCommonHistResults *gfcCommonHistRes8 = NULL; | |
494 | if(fileGFC) { | |
495 | fileGFC->GetObject("cobjGFC",pListGFC); | |
496 | if(pListGFC) { | |
497 | gfcCommonHist = dynamic_cast<AliFlowCommonHist*> (pListGFC->FindObject("AliFlowCommonHistGFC")); | |
498 | gfcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC")); | |
499 | if(gfcCommonHistRes2) { | |
dec3fab6 | 500 | flowValue[binGFC2-1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinContent(1); |
501 | flowError[binGFC2-1] = (gfcCommonHistRes2->GetHistIntFlow())->GetBinError(1); | |
502 | flowValueRP[binGFC2RP-1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1); | |
503 | flowErrorRP[binGFC2RP-1] = (gfcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1); | |
504 | flowValuePOI[binGFC2POI-1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1); | |
505 | flowErrorPOI[binGFC2POI-1] = (gfcCommonHistRes2->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 506 | } |
507 | gfcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults4thOrderGFC")); | |
508 | if(gfcCommonHistRes4) { | |
dec3fab6 | 509 | flowValue[binGFC4-1] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinContent(1); |
510 | flowError[binGFC4-1] = (gfcCommonHistRes4->GetHistIntFlow())->GetBinError(1); | |
511 | flowValueRP[binGFC4RP-1] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1); | |
512 | flowErrorRP[binGFC4RP-1] = (gfcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1); | |
513 | flowValuePOI[binGFC4POI-1] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1); | |
514 | flowErrorPOI[binGFC4POI-1] = (gfcCommonHistRes4->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 515 | } |
516 | gfcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults6thOrderGFC")); | |
517 | if(gfcCommonHistRes6) { | |
dec3fab6 | 518 | flowValue[binGFC6-1] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinContent(1); |
519 | flowError[binGFC6-1] = (gfcCommonHistRes6->GetHistIntFlow())->GetBinError(1); | |
520 | flowValueRP[binGFC6RP-1] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1); | |
521 | flowErrorRP[binGFC6RP-1] = (gfcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1); | |
522 | flowValuePOI[binGFC6POI-1] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1); | |
523 | flowErrorPOI[binGFC6POI-1] = (gfcCommonHistRes6->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 524 | } |
525 | gfcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListGFC->FindObject("AliFlowCommonHistResults8thOrderGFC")); | |
526 | if(gfcCommonHistRes8) { | |
dec3fab6 | 527 | flowValue[binGFC8-1] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinContent(1); |
528 | flowError[binGFC8-1] = (gfcCommonHistRes8->GetHistIntFlow())->GetBinError(1); | |
529 | flowValueRP[binGFC8RP-1] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1); | |
530 | flowErrorRP[binGFC8RP-1] = (gfcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1); | |
531 | flowValuePOI[binGFC8POI-1] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1); | |
532 | flowErrorPOI[binGFC8POI-1] = (gfcCommonHistRes8->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 533 | } |
534 | } | |
9d3cbee7 | 535 | } |
1fca9c90 | 536 | |
8c5ec46a | 537 | //QC = Q-cumulants |
538 | TList *pListQC = NULL; | |
539 | AliFlowCommonHist *qcCommonHist2 = NULL; | |
540 | AliFlowCommonHist *qcCommonHist4 = NULL; | |
541 | AliFlowCommonHist *qcCommonHist6 = NULL; | |
542 | AliFlowCommonHist *qcCommonHist8 = NULL; | |
543 | AliFlowCommonHistResults *qcCommonHistRes2 = NULL; | |
544 | AliFlowCommonHistResults *qcCommonHistRes4 = NULL; | |
545 | AliFlowCommonHistResults *qcCommonHistRes6 = NULL; | |
546 | AliFlowCommonHistResults *qcCommonHistRes8 = NULL; | |
547 | ||
548 | if(fileQC) { | |
549 | fileQC->GetObject("cobjQC",pListQC); | |
550 | if(pListQC) { | |
551 | qcCommonHist2 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist2ndOrderQC")); | |
552 | qcCommonHistRes2 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults2ndOrderQC")); | |
553 | if(qcCommonHistRes2) { | |
dec3fab6 | 554 | flowValue[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinContent(1); |
555 | //flowError[binQC2-1] = (qcCommonHistRes2->GetHistIntFlow())->GetBinError(1); | |
556 | flowValueRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinContent(1); | |
557 | //flowErrorRP[binQC2RP-1] = (qcCommonHistRes2->GetHistIntFlowRP())->GetBinError(1); | |
558 | flowValuePOI[binQC2POI-1] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinContent(1); | |
559 | //flowErrorPOI[binQC2POI-1] = (qcCommonHistRes2->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 560 | } |
561 | qcCommonHist4 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist4thOrderQC")); | |
562 | qcCommonHistRes4 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults4thOrderQC")); | |
563 | if(qcCommonHistRes4) { | |
dec3fab6 | 564 | flowValue[binQC4-1] = (qcCommonHistRes4->GetHistIntFlow())->GetBinContent(1); |
565 | //flowError[binQC4-1] = (qcCommonHistRes4->GetHistIntFlow())->GetBinError(1); | |
566 | flowValueRP[binQC4RP-1] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinContent(1); | |
567 | //flowErrorRP[binQC4RP-1] = (qcCommonHistRes4->GetHistIntFlowRP())->GetBinError(1); | |
568 | flowValuePOI[binQC4POI-1] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinContent(1); | |
569 | //flowErrorPOI[binQC4POI-1] = (qcCommonHistRes4->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 570 | } |
571 | qcCommonHist6 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist6thOrderQC")); | |
572 | qcCommonHistRes6 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults6thOrderQC")); | |
573 | if(qcCommonHistRes6) { | |
dec3fab6 | 574 | flowValue[binQC6-1] = (qcCommonHistRes6->GetHistIntFlow())->GetBinContent(1); |
575 | //flowError[binQC6-1] = (qcCommonHistRes6->GetHistIntFlow())->GetBinError(1); | |
576 | flowValueRP[binQC6RP-1] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinContent(1); | |
577 | //flowErrorRP[binQC6RP-1] = (qcCommonHistRes6->GetHistIntFlowRP())->GetBinError(1); | |
578 | flowValuePOI[binQC6POI-1] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinContent(1); | |
579 | //flowErrorPOI[binQC6POI-1] = (qcCommonHistRes6->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 580 | } |
581 | qcCommonHist8 = dynamic_cast<AliFlowCommonHist*> (pListQC->FindObject("AliFlowCommonHist8thOrderQC")); | |
582 | qcCommonHistRes8 = dynamic_cast<AliFlowCommonHistResults*> (pListQC->FindObject("AliFlowCommonHistResults8thOrderQC")); | |
583 | if(qcCommonHistRes8) { | |
dec3fab6 | 584 | flowValue[binQC8-1] = (qcCommonHistRes8->GetHistIntFlow())->GetBinContent(1); |
585 | //flowError[binQC8-1] = (qcCommonHistRes8->GetHistIntFlow())->GetBinError(1); | |
586 | flowValueRP[binQC8RP-1] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinContent(1); | |
587 | //flowErrorRP[binQC8RP-1] = (qcCommonHistRes8->GetHistIntFlowRP())->GetBinError(1); | |
588 | flowValuePOI[binQC8POI-1] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinContent(1); | |
589 | //flowErrorPOI[binQC8POI-1] = (qcCommonHistRes8->GetHistIntFlowPOI())->GetBinError(1); | |
8c5ec46a | 590 | } |
591 | } | |
592 | } | |
593 | ||
dec3fab6 | 594 | // ranges on y-axis for 'no-name' plot: |
595 | Double_t dMax=flowValue[binMC-1]+flowError[binMC-1]; | |
596 | Double_t dMin=flowValue[binMC-1]-flowError[binMC-1]; | |
8c5ec46a | 597 | |
598 | for(Int_t i=1;i<nMethods;i++) { | |
599 | if(!(flowValue[i]==0. && flowError[i]==0.)) { | |
600 | if(dMax<flowValue[i]+flowError[i]) dMax=flowValue[i]+flowError[i]; | |
601 | if(dMin>flowValue[i]-flowError[i]) dMin=flowValue[i]-flowError[i]; | |
602 | } | |
603 | } | |
a88e5ac2 | 604 | |
dec3fab6 | 605 | // ranges on y-axis for RP plot: |
606 | Double_t dMaxRP=flowValueRP[binMCRP-1]+flowErrorRP[binMCRP-1]; | |
607 | Double_t dMinRP=flowValueRP[binMCRP-1]-flowErrorRP[binMCRP-1]; | |
8c5ec46a | 608 | |
609 | for(Int_t i=1;i<nMethods;i++) { | |
610 | if(!(flowValueRP[i]==0. && flowErrorRP[i]==0.)) { | |
611 | if(dMaxRP<flowValueRP[i]+flowErrorRP[i]) dMaxRP=flowValueRP[i]+flowErrorRP[i]; | |
612 | if(dMinRP>flowValueRP[i]-flowErrorRP[i]) dMinRP=flowValueRP[i]-flowErrorRP[i]; | |
613 | } | |
614 | } | |
1fca9c90 | 615 | |
dec3fab6 | 616 | // ranges on y-axis for POI plot: |
617 | Double_t dMaxPOI=flowValuePOI[binMCPOI-1]+flowErrorPOI[binMCPOI-1]; | |
618 | Double_t dMinPOI=flowValuePOI[binMCPOI-1]-flowErrorPOI[binMCPOI-1]; | |
a88e5ac2 | 619 | |
8c5ec46a | 620 | for(Int_t i=1;i<nMethods;i++) { |
621 | if(!(flowValuePOI[i]==0. && flowErrorPOI[i]==0.)) { | |
622 | if(dMaxPOI<flowValuePOI[i]+flowErrorPOI[i]) dMaxPOI=flowValuePOI[i]+flowErrorPOI[i]; | |
623 | if(dMinPOI>flowValuePOI[i]-flowErrorPOI[i]) dMinPOI=flowValuePOI[i]-flowErrorPOI[i]; | |
624 | } | |
625 | } | |
626 | ||
dec3fab6 | 627 | // no-name: |
8c5ec46a | 628 | TGraph* flowResults = new TGraphErrors(nMethods, x, flowValue, xError, flowError); |
8c5ec46a | 629 | flowResults->SetMarkerStyle(markerStyle); |
630 | flowResults->SetMarkerColor(markerColor); | |
1fca9c90 | 631 | |
dec3fab6 | 632 | // RP: |
3963f25c | 633 | TGraph* flowResultsRP = new TGraphErrors(nMethods, x, flowValueRP, xError, flowErrorRP); |
dec3fab6 | 634 | flowResultsRP->SetMarkerStyle(markerStyleRP); |
635 | flowResultsRP->SetMarkerColor(markerColorRP); | |
8c5ec46a | 636 | |
dec3fab6 | 637 | // POI: |
8c5ec46a | 638 | TGraph* flowResultsPOI = new TGraphErrors(nMethods, x, flowValuePOI, xError, flowErrorPOI); |
dec3fab6 | 639 | flowResultsPOI->SetMarkerStyle(markerStylePOI); |
640 | flowResultsPOI->SetMarkerColor(markerColorPOI); | |
8c5ec46a | 641 | |
642 | //----------------------------------------------------------------------------------- | |
643 | ||
644 | //---------------------------------------------------------------------------------- | |
dec3fab6 | 645 | // cosmetics: mesh for MC error bands (integrated flow) |
3963f25c | 646 | TGraph* pMesh = NULL; |
8c5ec46a | 647 | TGraph* pMeshRP = NULL; |
648 | TGraph* pMeshPOI = NULL; | |
649 | ||
650 | if(intFlowAll && mcepCommonHistRes) { | |
8c5ec46a | 651 | Int_t nPts = nMethods; |
dec3fab6 | 652 | Double_t valueMC = flowValue[binMC-1]; |
653 | Double_t errorMC = flowError[binMC-1]; | |
654 | Double_t valueMCRP = flowValueRP[binMCRP-1]; | |
655 | Double_t errorMCRP = flowErrorRP[binMCRP-1]; | |
656 | Double_t valueMCPOI = flowValuePOI[binMCPOI-1]; | |
657 | Double_t errorMCPOI = flowErrorPOI[binMCPOI-1]; | |
8c5ec46a | 658 | |
3963f25c | 659 | pMesh = new TGraph(nPts); |
8c5ec46a | 660 | pMeshRP = new TGraph(nPts); |
661 | pMeshPOI = new TGraph(nPts); | |
662 | ||
dec3fab6 | 663 | // no-name: |
8c5ec46a | 664 | pMesh->SetPoint(1,0,valueMC+errorMC); |
665 | pMesh->SetPoint(2,nPts+1,valueMC+errorMC); | |
666 | pMesh->SetPoint(3,nPts+1,valueMC-errorMC); | |
667 | pMesh->SetPoint(4,0,valueMC-errorMC); | |
3963f25c | 668 | pMesh->SetPoint(5,0,valueMC+errorMC); |
8c5ec46a | 669 | pMesh->SetFillStyle(meshStyle); |
670 | pMesh->SetFillColor(meshColor); | |
3963f25c | 671 | |
dec3fab6 | 672 | // RP: |
8c5ec46a | 673 | pMeshRP->SetPoint(1,0,valueMCRP+errorMCRP); |
674 | pMeshRP->SetPoint(2,nPts+1,valueMCRP+errorMCRP); | |
675 | pMeshRP->SetPoint(3,nPts+1,valueMCRP-errorMCRP); | |
676 | pMeshRP->SetPoint(4,0,valueMCRP-errorMCRP); | |
3963f25c | 677 | pMeshRP->SetPoint(5,0,valueMCRP+errorMCRP); |
dec3fab6 | 678 | pMeshRP->SetFillStyle(meshStyleRP); |
679 | pMeshRP->SetFillColor(meshColorRP); | |
8c5ec46a | 680 | |
dec3fab6 | 681 | // POI: |
8c5ec46a | 682 | pMeshPOI->SetPoint(1,0,valueMCPOI+errorMCPOI); |
683 | pMeshPOI->SetPoint(2,nPts+1,valueMCPOI+errorMCPOI); | |
684 | pMeshPOI->SetPoint(3,nPts+1,valueMCPOI-errorMCPOI); | |
685 | pMeshPOI->SetPoint(4,0,valueMCPOI-errorMCPOI); | |
3963f25c | 686 | pMeshPOI->SetPoint(5,0,valueMCPOI+errorMCPOI); |
dec3fab6 | 687 | pMeshPOI->SetFillStyle(meshStylePOI); |
688 | pMeshPOI->SetFillColor(meshColorPOI); | |
8c5ec46a | 689 | } |
690 | //---------------------------------------------------------------------------------- | |
1fca9c90 | 691 | |
1fca9c90 | 692 | |
8c5ec46a | 693 | //---------------------------------------------------------------------------------- |
694 | //cosmetics: text (integrated flow) | |
695 | //default text: | |
696 | TPaveText *textDefault = new TPaveText(0.05,0.77,0.95,0.90,"NDC"); | |
697 | textDefault->SetTextFont(72); | |
698 | textDefault->SetTextSize(0.08); | |
8c5ec46a | 699 | |
700 | TString *entryDefaultAvM = new TString("Average Multiplicity"); | |
701 | TString *entryDefaultAnd = new TString("and"); | |
702 | TString *entryDefaultNumOfEvts = new TString("Number of Events:"); | |
703 | ||
704 | textDefault->AddText(entryDefaultAvM->Data()); | |
705 | textDefault->AddText(entryDefaultAnd->Data()); | |
706 | textDefault->AddText(entryDefaultNumOfEvts->Data()); | |
707 | ||
dec3fab6 | 708 | // results (no-name): |
3963f25c | 709 | TPaveText *textResults = new TPaveText(0.05,0.12,0.95,0.70,"NDC"); |
710 | textResults->SetTextFont(72); | |
711 | textResults->SetTextSize(0.06); | |
8c5ec46a | 712 | |
dec3fab6 | 713 | // results (RP): |
8c5ec46a | 714 | TPaveText *textResultsRP = new TPaveText(0.05,0.12,0.95,0.70,"NDC"); |
715 | textResultsRP->SetTextFont(72); | |
716 | textResultsRP->SetTextSize(0.06); | |
717 | ||
dec3fab6 | 718 | // results (POI): |
8c5ec46a | 719 | TPaveText *textResultsPOI = new TPaveText(0.05,0.12,0.95,0.70,"NDC"); |
720 | textResultsPOI->SetTextFont(72); | |
721 | textResultsPOI->SetTextSize(0.06); | |
722 | ||
dec3fab6 | 723 | // no-name: |
8c5ec46a | 724 | TString *entryMC = new TString("MC ...... "); |
905bed86 | 725 | TString *entrySP = new TString("SP ...... "); |
8c5ec46a | 726 | TString *entryGFC = new TString("GFC ..... "); |
727 | TString *entryQC2 = new TString("QC{2} ... "); | |
728 | TString *entryQC4 = new TString("QC{4} ... "); | |
729 | TString *entryQC6 = new TString("QC{6} ... "); | |
730 | TString *entryQC8 = new TString("QC{8} ... "); | |
731 | TString *entryFQD = new TString("FQD ..... "); | |
732 | TString *entryLYZ1 = new TString("LYZ ..... "); | |
733 | TString *entryLYZEP = new TString("LYZEP ... "); | |
8c5ec46a | 734 | |
dec3fab6 | 735 | // RP: |
8c5ec46a | 736 | TString *entryMCRP = new TString("MC ...... "); |
905bed86 | 737 | TString *entrySPRP = new TString("SP ...... "); |
8c5ec46a | 738 | TString *entryGFCRP = new TString("GFC ..... "); |
739 | TString *entryQC2RP = new TString("QC{2} ... "); | |
740 | TString *entryQC4RP = new TString("QC{4} ... "); | |
741 | TString *entryQC6RP = new TString("QC{6} ... "); | |
742 | TString *entryQC8RP = new TString("QC{8} ... "); | |
743 | TString *entryFQDRP = new TString("FQD ..... "); | |
744 | TString *entryLYZ1RP = new TString("LYZ ..... "); | |
745 | TString *entryLYZEPRP = new TString("LYZEP ... "); | |
1fca9c90 | 746 | |
dec3fab6 | 747 | // POI: |
8c5ec46a | 748 | TString *entryMCPOI = new TString("MC ...... "); |
905bed86 | 749 | TString *entrySPPOI = new TString("SP ...... "); |
8c5ec46a | 750 | TString *entryGFCPOI = new TString("GFC ..... "); |
751 | TString *entryQC2POI = new TString("QC{2} ... "); | |
752 | TString *entryQC4POI = new TString("QC{4} ... "); | |
753 | TString *entryQC6POI = new TString("QC{6} ... "); | |
754 | TString *entryQC8POI = new TString("QC{8} ... "); | |
755 | TString *entryFQDPOI = new TString("FQD ..... "); | |
756 | TString *entryLYZ1POI = new TString("LYZ ..... "); | |
757 | TString *entryLYZEPPOI = new TString("LYZEP ... "); | |
758 | ||
dec3fab6 | 759 | // no-name: |
8c5ec46a | 760 | Double_t avMultMC=0.; |
761 | Long_t nEvtsMC=0; | |
905bed86 | 762 | Double_t avMultSP=0.; |
763 | Long_t nEvtsSP=0; | |
8c5ec46a | 764 | Double_t avMultGFC=0.; |
765 | Long_t nEvtsGFC=0; | |
8c5ec46a | 766 | Double_t avMultQC2=0., avMultQC4=0., avMultQC6=0., avMultQC8=0.; |
767 | Long_t nEvtsQC2=0, nEvtsQC4=0, nEvtsQC6=0, nEvtsQC8=0; | |
8c5ec46a | 768 | Double_t avMultFQD=0.; |
769 | Long_t nEvtsFQD=0; | |
8c5ec46a | 770 | Double_t avMultLYZ1=0.; |
771 | Long_t nEvtsLYZ1=0; | |
8c5ec46a | 772 | Double_t avMultLYZEP=0.; |
773 | Long_t nEvtsLYZEP=0; | |
8c5ec46a | 774 | |
dec3fab6 | 775 | // RP: |
8c5ec46a | 776 | Double_t avMultMCRP=0.; |
777 | Long_t nEvtsMCRP=0; | |
905bed86 | 778 | Double_t avMultSPRP=0.; |
779 | Long_t nEvtsSPRP=0; | |
8c5ec46a | 780 | Double_t avMultGFCRP=0.; |
781 | Long_t nEvtsGFCRP=0; | |
8c5ec46a | 782 | Double_t avMultQC2RP=0., avMultQC4RP=0., avMultQC6RP=0., avMultQC8RP=0.; |
783 | Long_t nEvtsQC2RP=0, nEvtsQC4RP=0, nEvtsQC6RP=0, nEvtsQC8RP=0; | |
8c5ec46a | 784 | Double_t avMultFQDRP=0.; |
785 | Long_t nEvtsFQDRP=0; | |
8c5ec46a | 786 | Double_t avMultLYZ1RP=0.; |
787 | Long_t nEvtsLYZ1RP=0; | |
8c5ec46a | 788 | Double_t avMultLYZEPRP=0.; |
789 | Long_t nEvtsLYZEPRP=0; | |
790 | ||
dec3fab6 | 791 | // POI: |
8c5ec46a | 792 | Double_t avMultMCPOI=0.; |
793 | Long_t nEvtsMCPOI=0; | |
905bed86 | 794 | Double_t avMultSPPOI=0.; |
795 | Long_t nEvtsSPPOI=0; | |
8c5ec46a | 796 | Double_t avMultGFCPOI=0.; |
797 | Long_t nEvtsGFCPOI=0; | |
8c5ec46a | 798 | Double_t avMultQC2POI=0., avMultQC4POI=0., avMultQC6POI=0., avMultQC8POI=0.; |
799 | Long_t nEvtsQC2POI=0, nEvtsQC4POI=0, nEvtsQC6POI=0, nEvtsQC8POI=0; | |
8c5ec46a | 800 | Double_t avMultFQDPOI=0.; |
801 | Long_t nEvtsFQDPOI=0; | |
8c5ec46a | 802 | Double_t avMultLYZ1POI=0.; |
803 | Long_t nEvtsLYZ1POI=0; | |
8c5ec46a | 804 | Double_t avMultLYZEPPOI=0.; |
805 | Long_t nEvtsLYZEPPOI=0; | |
806 | ||
dec3fab6 | 807 | // MC: |
8c5ec46a | 808 | if(mcepCommonHist) { |
6a1a854d | 809 | avMultMC = (mcepCommonHist->GetHistMultRP())->GetMean(); |
810 | nEvtsMC = (mcepCommonHist->GetHistMultRP())->GetEntries(); | |
811 | avMultMCRP = (mcepCommonHist->GetHistMultRP())->GetMean(); | |
812 | nEvtsMCRP = (mcepCommonHist->GetHistMultRP())->GetEntries(); | |
813 | avMultMCPOI = (mcepCommonHist->GetHistMultPOI())->GetMean(); | |
814 | nEvtsMCPOI = (mcepCommonHist->GetHistMultPOI())->GetEntries(); | |
8c5ec46a | 815 | } |
3963f25c | 816 | |
8c5ec46a | 817 | if(entryMC) { |
818 | entryMC->Append("M = "); | |
819 | (*entryMC)+=(Long_t)avMultMC; | |
820 | entryMC->Append(", N = "); | |
821 | (*entryMC)+=(Long_t)nEvtsMC; | |
822 | } | |
3963f25c | 823 | |
8c5ec46a | 824 | if(entryMCRP) { |
825 | entryMCRP->Append("M = "); | |
826 | (*entryMCRP)+=(Long_t)avMultMCRP; | |
827 | entryMCRP->Append(", N = "); | |
828 | (*entryMCRP)+=(Long_t)nEvtsMCRP; | |
829 | } | |
3963f25c | 830 | |
831 | if(entryMCPOI) { | |
8c5ec46a | 832 | entryMCPOI->Append("M = "); |
833 | (*entryMCPOI)+=(Long_t)avMultMCPOI; | |
834 | entryMCPOI->Append(", N = "); | |
835 | (*entryMCPOI)+=(Long_t)nEvtsMCPOI; | |
3963f25c | 836 | } |
1fca9c90 | 837 | |
dec3fab6 | 838 | // SP: |
905bed86 | 839 | if(spCommonHist) { |
6a1a854d | 840 | avMultSP = (spCommonHist->GetHistMultRP())->GetMean(); |
841 | nEvtsSP = (spCommonHist->GetHistMultRP())->GetEntries(); | |
842 | avMultSPRP = (spCommonHist->GetHistMultRP())->GetMean(); | |
843 | nEvtsSPRP = (spCommonHist->GetHistMultRP())->GetEntries(); | |
844 | avMultSPPOI = (spCommonHist->GetHistMultPOI())->GetMean(); | |
845 | nEvtsSPPOI = (spCommonHist->GetHistMultPOI())->GetEntries(); | |
905bed86 | 846 | } |
847 | ||
848 | if(entrySP) { | |
849 | entrySP->Append("M = "); | |
850 | (*entrySP)+=(Long_t)avMultSP; | |
851 | entrySP->Append(", N = "); | |
852 | (*entrySP)+=(Long_t)nEvtsSP; | |
853 | } | |
854 | ||
855 | if(entrySPRP) { | |
856 | entrySPRP->Append("M = "); | |
857 | (*entrySPRP)+=(Long_t)avMultSPRP; | |
858 | entrySPRP->Append(", N = "); | |
859 | (*entrySPRP)+=(Long_t)nEvtsSPRP; | |
860 | } | |
861 | ||
862 | if(entrySPPOI) { | |
863 | entrySPPOI->Append("M = "); | |
864 | (*entrySPPOI)+=(Long_t)avMultSPPOI; | |
865 | entrySPPOI->Append(", N = "); | |
866 | (*entrySPPOI)+=(Long_t)nEvtsSPPOI; | |
867 | } | |
868 | ||
dec3fab6 | 869 | // GFC: |
8c5ec46a | 870 | if(gfcCommonHist) { |
6a1a854d | 871 | avMultGFC = (gfcCommonHist->GetHistMultRP())->GetMean(); |
872 | nEvtsGFC = (gfcCommonHist->GetHistMultRP())->GetEntries(); | |
873 | avMultGFCRP = (gfcCommonHist->GetHistMultRP())->GetMean(); | |
874 | nEvtsGFCRP = (gfcCommonHist->GetHistMultRP())->GetEntries(); | |
875 | avMultGFCPOI = (gfcCommonHist->GetHistMultPOI())->GetMean(); | |
876 | nEvtsGFCPOI = (gfcCommonHist->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 877 | } |
3963f25c | 878 | |
8c5ec46a | 879 | if(entryGFC) { |
880 | entryGFC->Append("M = "); | |
881 | (*entryGFC)+=(Long_t)avMultGFC; | |
882 | entryGFC->Append(", N = "); | |
883 | (*entryGFC)+=(Long_t)nEvtsGFC; | |
a88e5ac2 | 884 | } |
3963f25c | 885 | |
8c5ec46a | 886 | if(entryGFCRP) { |
887 | entryGFCRP->Append("M = "); | |
888 | (*entryGFCRP)+=(Long_t)avMultGFCRP; | |
889 | entryGFCRP->Append(", N = "); | |
890 | (*entryGFCRP)+=(Long_t)nEvtsGFCRP; | |
1fca9c90 | 891 | } |
8c5ec46a | 892 | if(entryGFCPOI) { |
893 | entryGFCPOI->Append("M = "); | |
894 | (*entryGFCPOI)+=(Long_t)avMultGFCPOI; | |
895 | entryGFCPOI->Append(", N = "); | |
896 | (*entryGFCPOI)+=(Long_t)nEvtsGFCPOI; | |
1fca9c90 | 897 | } |
a88e5ac2 | 898 | |
dec3fab6 | 899 | // QC: |
8c5ec46a | 900 | if(qcCommonHist2) { |
6a1a854d | 901 | avMultQC2 = (qcCommonHist2->GetHistMultRP())->GetMean(); |
902 | nEvtsQC2 = (qcCommonHist2->GetHistMultRP())->GetEntries(); | |
903 | avMultQC2RP = (qcCommonHist2->GetHistMultRP())->GetMean(); | |
904 | nEvtsQC2RP = (qcCommonHist2->GetHistMultRP())->GetEntries(); | |
905 | avMultQC2POI = (qcCommonHist2->GetHistMultPOI())->GetMean(); | |
906 | nEvtsQC2POI = (qcCommonHist2->GetHistMultPOI())->GetEntries(); | |
1fca9c90 | 907 | } |
3963f25c | 908 | |
a88e5ac2 | 909 | if(entryQC2) |
910 | { | |
70288839 | 911 | entryQC2->Append("M = "); |
912 | (*entryQC2)+=(Long_t)avMultQC2; | |
913 | entryQC2->Append(", N = "); | |
914 | (*entryQC2)+=(Long_t)nEvtsQC2; | |
70288839 | 915 | } |
3963f25c | 916 | |
1fca9c90 | 917 | if(entryQC2RP) |
918 | { | |
919 | entryQC2RP->Append("M = "); | |
920 | (*entryQC2RP)+=(Long_t)avMultQC2RP; | |
921 | entryQC2RP->Append(", N = "); | |
922 | (*entryQC2RP)+=(Long_t)nEvtsQC2RP; | |
923 | } | |
3963f25c | 924 | |
1fca9c90 | 925 | if(entryQC2POI) |
926 | { | |
927 | entryQC2POI->Append("M = "); | |
928 | (*entryQC2POI)+=(Long_t)avMultQC2POI; | |
929 | entryQC2POI->Append(", N = "); | |
930 | (*entryQC2POI)+=(Long_t)nEvtsQC2POI; | |
931 | } | |
932 | ||
70288839 | 933 | if(qcCommonHist4) |
934 | { | |
6a1a854d | 935 | avMultQC4 = (qcCommonHist4->GetHistMultRP())->GetMean(); |
936 | nEvtsQC4 = (qcCommonHist4->GetHistMultRP())->GetEntries(); | |
937 | avMultQC4RP = (qcCommonHist4->GetHistMultRP())->GetMean(); | |
938 | nEvtsQC4RP = (qcCommonHist4->GetHistMultRP())->GetEntries(); | |
939 | avMultQC4POI = (qcCommonHist4->GetHistMultPOI())->GetMean(); | |
940 | nEvtsQC4POI = (qcCommonHist4->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 941 | } |
3963f25c | 942 | |
a88e5ac2 | 943 | if(entryQC4) |
944 | { | |
70288839 | 945 | entryQC4->Append("M = "); |
946 | (*entryQC4)+=(Long_t)avMultQC4; | |
947 | entryQC4->Append(", N = "); | |
948 | (*entryQC4)+=(Long_t)nEvtsQC4; | |
70288839 | 949 | } |
3963f25c | 950 | |
1fca9c90 | 951 | if(entryQC4RP) |
952 | { | |
953 | entryQC4RP->Append("M = "); | |
954 | (*entryQC4RP)+=(Long_t)avMultQC4RP; | |
955 | entryQC4RP->Append(", N = "); | |
956 | (*entryQC4RP)+=(Long_t)nEvtsQC4RP; | |
957 | } | |
3963f25c | 958 | |
1fca9c90 | 959 | if(entryQC4POI) |
960 | { | |
961 | entryQC4POI->Append("M = "); | |
962 | (*entryQC4POI)+=(Long_t)avMultQC4POI; | |
963 | entryQC4POI->Append(", N = "); | |
964 | (*entryQC4POI)+=(Long_t)nEvtsQC4POI; | |
965 | } | |
966 | ||
a88e5ac2 | 967 | if(qcCommonHist6) |
968 | { | |
6a1a854d | 969 | avMultQC6 = (qcCommonHist6->GetHistMultRP())->GetMean(); |
970 | nEvtsQC6 = (qcCommonHist6->GetHistMultRP())->GetEntries(); | |
971 | avMultQC6RP = (qcCommonHist6->GetHistMultRP())->GetMean(); | |
972 | nEvtsQC6RP = (qcCommonHist6->GetHistMultRP())->GetEntries(); | |
973 | avMultQC6POI = (qcCommonHist6->GetHistMultPOI())->GetMean(); | |
974 | nEvtsQC6POI = (qcCommonHist6->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 975 | } |
3963f25c | 976 | |
a88e5ac2 | 977 | if(entryQC6) |
978 | { | |
979 | entryQC6->Append("M = "); | |
980 | (*entryQC6)+=(Long_t)avMultQC6; | |
981 | entryQC6->Append(", N = "); | |
982 | (*entryQC6)+=(Long_t)nEvtsQC6; | |
983 | } | |
3963f25c | 984 | |
1fca9c90 | 985 | if(entryQC6RP) |
986 | { | |
987 | entryQC6RP->Append("M = "); | |
988 | (*entryQC6RP)+=(Long_t)avMultQC6RP; | |
989 | entryQC6RP->Append(", N = "); | |
990 | (*entryQC6RP)+=(Long_t)nEvtsQC6RP; | |
991 | } | |
3963f25c | 992 | |
1fca9c90 | 993 | if(entryQC6POI) |
994 | { | |
995 | entryQC6POI->Append("M = "); | |
996 | (*entryQC6POI)+=(Long_t)avMultQC6POI; | |
997 | entryQC6POI->Append(", N = "); | |
998 | (*entryQC6POI)+=(Long_t)nEvtsQC6POI; | |
999 | } | |
1000 | ||
a88e5ac2 | 1001 | if(qcCommonHist8) |
1002 | { | |
6a1a854d | 1003 | avMultQC8 = (qcCommonHist8->GetHistMultRP())->GetMean(); |
1004 | nEvtsQC8 = (qcCommonHist8->GetHistMultRP())->GetEntries(); | |
1005 | avMultQC8RP = (qcCommonHist8->GetHistMultRP())->GetMean(); | |
1006 | nEvtsQC8RP = (qcCommonHist8->GetHistMultRP())->GetEntries(); | |
1007 | avMultQC8POI = (qcCommonHist8->GetHistMultPOI())->GetMean(); | |
1008 | nEvtsQC8POI = (qcCommonHist8->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 1009 | } |
3963f25c | 1010 | |
a88e5ac2 | 1011 | if(entryQC8) |
1012 | { | |
1013 | entryQC8->Append("M = "); | |
1014 | (*entryQC8)+=(Long_t)avMultQC8; | |
1015 | entryQC8->Append(", N = "); | |
1016 | (*entryQC8)+=(Long_t)nEvtsQC8; | |
1017 | } | |
3963f25c | 1018 | |
1fca9c90 | 1019 | if(entryQC8RP) |
1020 | { | |
1021 | entryQC8RP->Append("M = "); | |
1022 | (*entryQC8RP)+=(Long_t)avMultQC8RP; | |
1023 | entryQC8RP->Append(", N = "); | |
1024 | (*entryQC8RP)+=(Long_t)nEvtsQC8RP; | |
1025 | } | |
3963f25c | 1026 | |
1fca9c90 | 1027 | if(entryQC8POI) |
1028 | { | |
1029 | entryQC8POI->Append("M = "); | |
1030 | (*entryQC8POI)+=(Long_t)avMultQC8POI; | |
1031 | entryQC8POI->Append(", N = "); | |
1032 | (*entryQC8POI)+=(Long_t)nEvtsQC8POI; | |
1033 | } | |
1034 | ||
dec3fab6 | 1035 | // FQD: |
a88e5ac2 | 1036 | if(fqdCommonHist) |
1037 | { | |
6a1a854d | 1038 | avMultFQD = (fqdCommonHist->GetHistMultRP())->GetMean(); |
1039 | nEvtsFQD = (fqdCommonHist->GetHistMultRP())->GetEntries(); | |
1040 | avMultFQDRP = (fqdCommonHist->GetHistMultRP())->GetMean(); | |
1041 | nEvtsFQDRP = (fqdCommonHist->GetHistMultRP())->GetEntries(); | |
1042 | avMultFQDPOI = (fqdCommonHist->GetHistMultPOI())->GetMean(); | |
1043 | nEvtsFQDPOI = (fqdCommonHist->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 1044 | } |
3963f25c | 1045 | |
a88e5ac2 | 1046 | if(entryFQD) |
1047 | { | |
1048 | entryFQD->Append("M = "); | |
1049 | (*entryFQD)+=(Long_t)avMultFQD; | |
1050 | entryFQD->Append(", N = "); | |
1051 | (*entryFQD)+=(Long_t)nEvtsFQD; | |
1052 | } | |
3963f25c | 1053 | |
1fca9c90 | 1054 | if(entryFQDRP) |
1055 | { | |
1056 | entryFQDRP->Append("M = "); | |
1057 | (*entryFQDRP)+=(Long_t)avMultFQDRP; | |
1058 | entryFQDRP->Append(", N = "); | |
1059 | (*entryFQDRP)+=(Long_t)nEvtsFQDRP; | |
1060 | } | |
3963f25c | 1061 | |
1fca9c90 | 1062 | if(entryFQDPOI) |
1063 | { | |
1064 | entryFQDPOI->Append("M = "); | |
1065 | (*entryFQDPOI)+=(Long_t)avMultFQDPOI; | |
1066 | entryFQDPOI->Append(", N = "); | |
1067 | (*entryFQDPOI)+=(Long_t)nEvtsFQDPOI; | |
1068 | } | |
1069 | ||
dec3fab6 | 1070 | // LYZ1: |
a88e5ac2 | 1071 | if(lyz1CommonHist) |
1072 | { | |
6a1a854d | 1073 | avMultLYZ1 = (lyz1CommonHist->GetHistMultRP())->GetMean(); |
1074 | nEvtsLYZ1 = (lyz1CommonHist->GetHistMultRP())->GetEntries(); | |
1075 | avMultLYZ1RP = (lyz1CommonHist->GetHistMultRP())->GetMean(); | |
1076 | nEvtsLYZ1RP = (lyz1CommonHist->GetHistMultRP())->GetEntries(); | |
1077 | avMultLYZ1POI = (lyz1CommonHist->GetHistMultPOI())->GetMean(); | |
1078 | nEvtsLYZ1POI = (lyz1CommonHist->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 1079 | } |
3963f25c | 1080 | |
a88e5ac2 | 1081 | if(entryLYZ1) |
1082 | { | |
1083 | entryLYZ1->Append("M = "); | |
1084 | (*entryLYZ1)+=(Long_t)avMultLYZ1; | |
1085 | entryLYZ1->Append(", N = "); | |
1086 | (*entryLYZ1)+=(Long_t)nEvtsLYZ1; | |
1087 | } | |
3963f25c | 1088 | |
1fca9c90 | 1089 | if(entryLYZ1RP) |
1090 | { | |
1091 | entryLYZ1RP->Append("M = "); | |
1092 | (*entryLYZ1RP)+=(Long_t)avMultLYZ1RP; | |
1093 | entryLYZ1RP->Append(", N = "); | |
1094 | (*entryLYZ1RP)+=(Long_t)nEvtsLYZ1RP; | |
1095 | } | |
3963f25c | 1096 | |
1fca9c90 | 1097 | if(entryLYZ1POI) |
1098 | { | |
1099 | entryLYZ1POI->Append("M = "); | |
1100 | (*entryLYZ1POI)+=(Long_t)avMultLYZ1POI; | |
1101 | entryLYZ1POI->Append(", N = "); | |
1102 | (*entryLYZ1POI)+=(Long_t)nEvtsLYZ1POI; | |
1103 | } | |
70288839 | 1104 | |
dec3fab6 | 1105 | // LYZEP: |
a88e5ac2 | 1106 | if(lyzepCommonHist) |
70288839 | 1107 | { |
6a1a854d | 1108 | avMultLYZEP = (lyzepCommonHist->GetHistMultRP())->GetMean(); |
1109 | nEvtsLYZEP = (lyzepCommonHist->GetHistMultRP())->GetEntries(); | |
1110 | avMultLYZEPRP = (lyzepCommonHist->GetHistMultRP())->GetMean(); | |
1111 | nEvtsLYZEPRP = (lyzepCommonHist->GetHistMultRP())->GetEntries(); | |
1112 | avMultLYZEPPOI = (lyzepCommonHist->GetHistMultPOI())->GetMean(); | |
1113 | nEvtsLYZEPPOI = (lyzepCommonHist->GetHistMultPOI())->GetEntries(); | |
70288839 | 1114 | } |
3963f25c | 1115 | |
a88e5ac2 | 1116 | if(entryLYZEP) |
70288839 | 1117 | { |
a88e5ac2 | 1118 | entryLYZEP->Append("M = "); |
1119 | (*entryLYZEP)+=(Long_t)avMultLYZEP; | |
1120 | entryLYZEP->Append(", N = "); | |
1121 | (*entryLYZEP)+=(Long_t)nEvtsLYZEP; | |
70288839 | 1122 | } |
3963f25c | 1123 | |
1fca9c90 | 1124 | if(entryLYZEPRP) |
1125 | { | |
1126 | entryLYZEPRP->Append("M = "); | |
1127 | (*entryLYZEPRP)+=(Long_t)avMultLYZEPRP; | |
1128 | entryLYZEPRP->Append(", N = "); | |
1129 | (*entryLYZEPRP)+=(Long_t)nEvtsLYZEPRP; | |
1130 | } | |
3963f25c | 1131 | |
1fca9c90 | 1132 | if(entryLYZEPPOI) |
1133 | { | |
1134 | entryLYZEPPOI->Append("M = "); | |
1135 | (*entryLYZEPPOI)+=(Long_t)avMultLYZEPPOI; | |
1136 | entryLYZEPPOI->Append(", N = "); | |
1137 | (*entryLYZEPPOI)+=(Long_t)nEvtsLYZEPPOI; | |
1138 | } | |
1139 | ||
dec3fab6 | 1140 | // no-name: |
a88e5ac2 | 1141 | if(textResults) |
70288839 | 1142 | { |
a88e5ac2 | 1143 | textResults->AddText(entryMC->Data()); |
905bed86 | 1144 | textResults->AddText(entrySP->Data()); |
a88e5ac2 | 1145 | textResults->AddText(entryGFC->Data()); |
1146 | textResults->AddText(entryQC2->Data()); | |
1147 | textResults->AddText(entryQC4->Data()); | |
1148 | textResults->AddText(entryQC6->Data()); | |
1149 | textResults->AddText(entryQC8->Data()); | |
1150 | textResults->AddText(entryFQD->Data()); | |
1151 | textResults->AddText(entryLYZ1->Data()); | |
1152 | textResults->AddText(entryLYZEP->Data()); | |
70288839 | 1153 | } |
3963f25c | 1154 | |
dec3fab6 | 1155 | // RP: |
1fca9c90 | 1156 | if(textResultsRP) |
1157 | { | |
1158 | textResultsRP->AddText(entryMCRP->Data()); | |
905bed86 | 1159 | textResultsRP->AddText(entrySPRP->Data()); |
1fca9c90 | 1160 | textResultsRP->AddText(entryGFCRP->Data()); |
1161 | textResultsRP->AddText(entryQC2RP->Data()); | |
1162 | textResultsRP->AddText(entryQC4RP->Data()); | |
1163 | textResultsRP->AddText(entryQC6RP->Data()); | |
1164 | textResultsRP->AddText(entryQC8RP->Data()); | |
1165 | textResultsRP->AddText(entryFQDRP->Data()); | |
1166 | textResultsRP->AddText(entryLYZ1RP->Data()); | |
1167 | textResultsRP->AddText(entryLYZEPRP->Data()); | |
1168 | } | |
dec3fab6 | 1169 | |
1170 | // POI: | |
1fca9c90 | 1171 | if(textResultsPOI) |
1172 | { | |
1173 | textResultsPOI->AddText(entryMCPOI->Data()); | |
905bed86 | 1174 | textResultsPOI->AddText(entrySPPOI->Data()); |
1fca9c90 | 1175 | textResultsPOI->AddText(entryGFCPOI->Data()); |
1176 | textResultsPOI->AddText(entryQC2POI->Data()); | |
1177 | textResultsPOI->AddText(entryQC4POI->Data()); | |
1178 | textResultsPOI->AddText(entryQC6POI->Data()); | |
1179 | textResultsPOI->AddText(entryQC8POI->Data()); | |
1180 | textResultsPOI->AddText(entryFQDPOI->Data()); | |
1181 | textResultsPOI->AddText(entryLYZ1POI->Data()); | |
1182 | textResultsPOI->AddText(entryLYZEPPOI->Data()); | |
1183 | } | |
a88e5ac2 | 1184 | //---------------------------------------------------------------------------------- |
70288839 | 1185 | |
3963f25c | 1186 | |
1187 | ||
1188 | ||
70288839 | 1189 | //---------------------------------------------------------------------------------- |
da24213b | 1190 | // final drawing for integrated flow (no-name): |
1191 | if(plotIntFlow) | |
1192 | { | |
1193 | TCanvas* intFlowAllCanvas = new TCanvas("Integrated Flow","Integrated Flow",1000,600); | |
6f62c2a7 | 1194 | |
da24213b | 1195 | intFlowAllCanvas->Divide(2,1); |
9d3cbee7 | 1196 | |
da24213b | 1197 | // 1st pad is for plot: |
1198 | (intFlowAllCanvas->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
1199 | ||
1200 | if(intFlowAll) | |
a88e5ac2 | 1201 | { |
da24213b | 1202 | if(dMin>0. && dMax>0.) |
1203 | { | |
1204 | (intFlowAll->GetYaxis())->SetRangeUser(0.9744*dMin,1.0144*dMax); | |
1205 | } else if(dMin<0. && dMax>0.) | |
1206 | { | |
1207 | if(!(-1.*dMin<4.*dMax)) | |
1208 | { | |
1209 | (intFlowAll->GetYaxis())->SetRangeUser(1.0266*dMin,1.0144*dMax); | |
1210 | } else {(intFlowAll->GetYaxis())->SetRangeUser(1.1266*dMin,1.0144*dMax);} | |
1211 | } else if(dMin<0. && dMax<0.) | |
1212 | { | |
1213 | (intFlowAll->GetYaxis())->SetRangeUser(1.0266*dMin,0.9866*dMax); | |
1214 | } | |
1215 | intFlowAll->Draw("E1"); | |
1216 | } | |
e83922f2 | 1217 | |
da24213b | 1218 | if(pMesh) pMesh->Draw("LFSAME"); |
1219 | ||
1220 | if(flowResults) flowResults->Draw("PSAME"); | |
a88e5ac2 | 1221 | |
da24213b | 1222 | // 2nd pad is for legend: |
1223 | (intFlowAllCanvas->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
a88e5ac2 | 1224 | |
dec3fab6 | 1225 | if(textDefault) textDefault->Draw(); |
1226 | ||
1227 | if(textResults) textResults->Draw(); | |
1228 | ||
da24213b | 1229 | }// end of if(plotIntFlow) |
a88e5ac2 | 1230 | //---------------------------------------------------------------------------------- |
1fca9c90 | 1231 | |
e83922f2 | 1232 | |
1233 | ||
e83922f2 | 1234 | |
1fca9c90 | 1235 | //---------------------------------------------------------------------------------- |
1236 | //final drawing for integrated flow of RP (i.e. of particles used to determine the reaction plane): | |
e83922f2 | 1237 | if(plotIntFlowRP) |
1238 | { | |
1239 | TCanvas* intFlowAllCanvasRP = new TCanvas("Integrated Flow RP","Integrated Flow RP",1000,600); | |
1fca9c90 | 1240 | |
e83922f2 | 1241 | intFlowAllCanvasRP->Divide(2,1); |
1fca9c90 | 1242 | |
e83922f2 | 1243 | //1st pad is for plot: |
1244 | (intFlowAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
1fca9c90 | 1245 | |
dec3fab6 | 1246 | TH1D *intFlowAllRP = new TH1D(*intFlowAll); |
1247 | intFlowAllRP->SetMarkerStyle(markerStyleRP); | |
1248 | intFlowAllRP->SetMarkerColor(markerColorRP); | |
1249 | (intFlowAllRP->GetXaxis())->SetBinLabel(binMCRP,"v_{2}{MC}"); | |
1250 | (intFlowAllRP->GetXaxis())->SetBinLabel(binSPRP,"v_{2}{SP}"); | |
1251 | (intFlowAllRP->GetXaxis())->SetBinLabel(binGFC2RP,"v_{2}{2,GFC}"); | |
1252 | (intFlowAllRP->GetXaxis())->SetBinLabel(binQC2RP,"v_{2}{2,QC}"); | |
1253 | (intFlowAllRP->GetXaxis())->SetBinLabel(binGFC4RP,"v_{2}{4,GFC}"); | |
1254 | (intFlowAllRP->GetXaxis())->SetBinLabel(binQC4RP,"v_{2}{4,QC}"); | |
1255 | (intFlowAllRP->GetXaxis())->SetBinLabel(binGFC6RP,"v_{2}{6,GFC}"); | |
1256 | (intFlowAllRP->GetXaxis())->SetBinLabel(binQC6RP,"v_{2}{6,QC}"); | |
1257 | (intFlowAllRP->GetXaxis())->SetBinLabel(binGFC8RP,"v_{2}{8,GFC}"); | |
1258 | (intFlowAllRP->GetXaxis())->SetBinLabel(binQC8RP,"v_{2}{8,QC}"); | |
1259 | (intFlowAllRP->GetXaxis())->SetBinLabel(binFQDRP,"v_{2}{FQD}"); | |
1260 | (intFlowAllRP->GetXaxis())->SetBinLabel(binLYZ2RP,"v_{2}{LYZ}"); | |
1261 | (intFlowAllRP->GetXaxis())->SetBinLabel(binLYZEPRP,"v_{2}{LYZEP}"); | |
694e0378 | 1262 | |
e83922f2 | 1263 | if(intFlowAllRP) |
1264 | { | |
1265 | if(dMinRP>0. && dMaxRP>0.) | |
1266 | { | |
1267 | (intFlowAllRP->GetYaxis())->SetRangeUser(0.9744*dMinRP,1.0144*dMaxRP); | |
1268 | } else if(dMinRP<0. && dMaxRP>0.) | |
1269 | { | |
1270 | if(!(-1.*dMinRP<4.*dMaxRP)) | |
1271 | { | |
1272 | (intFlowAllRP->GetYaxis())->SetRangeUser(1.0266*dMinRP,1.0144*dMaxRP); | |
1273 | } else {(intFlowAllRP->GetYaxis())->SetRangeUser(1.1266*dMinRP,1.0144*dMaxRP);} | |
1274 | } else if(dMinRP<0. && dMaxRP<0.) | |
1275 | { | |
1276 | (intFlowAllRP->GetYaxis())->SetRangeUser(1.0266*dMinRP,0.9866*dMaxRP); | |
1277 | } | |
1278 | intFlowAllRP->Draw("E1"); | |
1279 | } | |
694e0378 | 1280 | |
e83922f2 | 1281 | if(pMeshRP) pMeshRP->Draw("LFSAME"); |
1282 | ||
1283 | if(flowResultsRP) flowResultsRP->Draw("PSAME"); | |
1fca9c90 | 1284 | |
e83922f2 | 1285 | //2nd pad is for legend: |
1286 | (intFlowAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
1287 | ||
1288 | if(textDefault) textDefault->Draw(); | |
1fca9c90 | 1289 | |
e83922f2 | 1290 | if(textResultsRP) textResultsRP->Draw(); |
1291 | }//end of if(plotIntFlowRP} | |
1fca9c90 | 1292 | //---------------------------------------------------------------------------------- |
1293 | ||
3963f25c | 1294 | |
1295 | ||
1296 | ||
1fca9c90 | 1297 | //---------------------------------------------------------------------------------- |
1298 | //final drawing for integrated flow of POI (i.e. of particles of interest): | |
e83922f2 | 1299 | if(plotIntFlowPOI) |
1300 | { | |
1301 | TCanvas* intFlowAllCanvasPOI = new TCanvas("Integrated Flow POI","Integrated Flow POI",1000,600); | |
dec3fab6 | 1302 | |
e83922f2 | 1303 | intFlowAllCanvasPOI->Divide(2,1); |
1fca9c90 | 1304 | |
e83922f2 | 1305 | //1st pad is for plot: |
1306 | (intFlowAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
1307 | ||
694e0378 | 1308 | TH1D *intFlowAllPOI = new TH1D(*intFlowAll); |
dec3fab6 | 1309 | intFlowAllPOI->SetMarkerStyle(markerStylePOI); |
1310 | intFlowAllPOI->SetMarkerColor(markerColorPOI); | |
1311 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binMCPOI,"v_{2}{MC}"); | |
1312 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binSPPOI,"v_{2}{SP}"); | |
1313 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binGFC2POI,"v_{2}{2,GFC}"); | |
1314 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binQC2POI,"v_{2}{2,QC}"); | |
1315 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binGFC4POI,"v_{2}{4,GFC}"); | |
1316 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binQC4POI,"v_{2}{4,QC}"); | |
1317 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binGFC6POI,"v_{2}{6,GFC}"); | |
1318 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binQC6POI,"v_{2}{6,QC}"); | |
1319 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binGFC8POI,"v_{2}{8,GFC}"); | |
1320 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binQC8POI,"v_{2}{8,QC}"); | |
1321 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binFQDPOI,"v_{2}{FQD}"); | |
1322 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binLYZ2POI,"v_{2}{LYZ}"); | |
1323 | (intFlowAllPOI->GetXaxis())->SetBinLabel(binLYZEPPOI,"v_{2}{LYZEP}"); | |
694e0378 | 1324 | |
e83922f2 | 1325 | if(intFlowAllPOI) |
1fca9c90 | 1326 | { |
e83922f2 | 1327 | if(dMinPOI>0. && dMaxPOI>0.) |
1328 | { | |
1329 | (intFlowAllPOI->GetYaxis())->SetRangeUser(0.9744*dMinPOI,1.0144*dMaxPOI); | |
1330 | } else if(dMinPOI<0. && dMaxPOI>0.) | |
1331 | { | |
1332 | if(!(-1.*dMinPOI<4.*dMaxPOI)) | |
1333 | { | |
1334 | (intFlowAllPOI->GetYaxis())->SetRangeUser(1.0266*dMinPOI,1.0144*dMaxPOI); | |
1335 | } else {(intFlowAllPOI->GetYaxis())->SetRangeUser(1.1266*dMinPOI,1.0144*dMaxPOI);} | |
1336 | } else if(dMinPOI<0. && dMaxPOI<0.) | |
1337 | { | |
1338 | (intFlowAllPOI->GetYaxis())->SetRangeUser(1.0266*dMinPOI,0.9866*dMaxPOI); | |
1339 | } | |
1340 | intFlowAllPOI->Draw("E1"); | |
1341 | } | |
1fca9c90 | 1342 | |
e83922f2 | 1343 | if(pMeshPOI) pMeshPOI->Draw("LFSAME"); |
1344 | ||
1345 | if(flowResultsPOI) flowResultsPOI->Draw("PSAME"); | |
1fca9c90 | 1346 | |
e83922f2 | 1347 | //2nd pad is for legend: |
1348 | (intFlowAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
1349 | ||
1350 | if(textDefault) textDefault->Draw(); | |
1fca9c90 | 1351 | |
e83922f2 | 1352 | if(textResultsPOI) textResultsPOI->Draw(); |
1353 | }// end of if(plotIntFlowPOI) | |
1fca9c90 | 1354 | //---------------------------------------------------------------------------------- |
1355 | ||
9d3cbee7 | 1356 | //================================================================================== |
e83922f2 | 1357 | |
1358 | ||
1359 | ||
9d3cbee7 | 1360 | |
1361 | //================================================================================== | |
1362 | // DIFFERENTIAL FLOW | |
1363 | //================================================================================== | |
70288839 | 1364 | Int_t iNbinsPt = AliFlowCommonConstants::GetNbinsPt(); |
1365 | Double_t dPtMin = AliFlowCommonConstants::GetPtMin(); | |
1366 | Double_t dPtMax = AliFlowCommonConstants::GetPtMax(); | |
1367 | ||
1fca9c90 | 1368 | Int_t iNbinsEta = AliFlowCommonConstants::GetNbinsEta(); |
1369 | Double_t dEtaMin = AliFlowCommonConstants::GetEtaMin(); | |
1370 | Double_t dEtaMax = AliFlowCommonConstants::GetEtaMax(); | |
1371 | ||
a88e5ac2 | 1372 | //---------------------------------------------------------------------------------- |
1fca9c90 | 1373 | //cosmetics: the style histogram for differential flow (pt): |
1374 | TH1D *styleHistPt = new TH1D("styleHistPt","styleHistPt",iNbinsPt,dPtMin,dPtMax); | |
1375 | styleHistPt->SetTitle("Differential Flow"); | |
1376 | styleHistPt->SetXTitle("p_{t} [GeV]"); | |
1377 | styleHistPt->SetYTitle("v_{n}"); | |
1378 | ||
1379 | //cosmetics: the style histogram for differential flow (eta): | |
1380 | TH1D *styleHistEta = new TH1D("styleHistEta","styleHistEta",iNbinsEta,dEtaMin,dEtaMax); | |
1381 | styleHistEta->SetTitle("Differential Flow"); | |
1382 | styleHistEta->SetXTitle("#eta"); | |
1383 | styleHistEta->SetYTitle("v_{n}"); | |
a88e5ac2 | 1384 | //---------------------------------------------------------------------------------- |
9d3cbee7 | 1385 | |
1fca9c90 | 1386 | |
1fca9c90 | 1387 | |
1fca9c90 | 1388 | //---------------------------------------------------------------------------------- |
1389 | //RP: | |
1390 | //cosmetics: Monte Carlo error bands for differential flow (Pt) | |
1391 | TGraph* pMeshDiffFlowPtRP = NULL; | |
1392 | if(mcepCommonHistRes) | |
1393 | { | |
7d02d77f | 1394 | Int_t nBinsDiffFlowPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetNbinsX(); |
1fca9c90 | 1395 | Double_t binWidthPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinWidth(1);//assuming that all bins have the same width |
7d02d77f | 1396 | |
1397 | //counting the non-empty bins: | |
1398 | Int_t nNonEmptyBinsDiffFlowPtRP=0; | |
1399 | for(Int_t i=1;i<nBinsDiffFlowPtRP+1;i++) | |
1400 | { | |
1401 | if(!(mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinError(i)==0.0)) | |
1402 | { | |
1403 | nNonEmptyBinsDiffFlowPtRP++; | |
1404 | } | |
1405 | } | |
1fca9c90 | 1406 | |
7d02d77f | 1407 | pMeshDiffFlowPtRP = new TGraph(2*nNonEmptyBinsDiffFlowPtRP+1); |
1fca9c90 | 1408 | |
7d02d77f | 1409 | Double_t valueMCPtRP=0.,errorMCPtRP=0.; |
1410 | Int_t countDiffFlowPtRP=1; | |
1411 | Double_t xFirstDiffFlowPtRP=0.,yUpFirstDiffFlowPtRP=0.;//needed to close up the mesh | |
1412 | for(Int_t i=1;i<nBinsDiffFlowPtRP+1;i++) | |
1fca9c90 | 1413 | { |
7d02d77f | 1414 | //setting up the upper limit of the mesh: |
1fca9c90 | 1415 | valueMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(i); |
7d02d77f | 1416 | errorMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinError(i); |
1417 | if(!(errorMCPtRP==0.0)) | |
1418 | { | |
1419 | pMeshDiffFlowPtRP->SetPoint(countDiffFlowPtRP++,(i-0.5)*binWidthPtRP+dPtMin,valueMCPtRP+errorMCPtRP); | |
1420 | if(xFirstDiffFlowPtRP==0.) | |
1421 | { | |
1422 | xFirstDiffFlowPtRP=(i-0.5)*binWidthPtRP+dPtMin; | |
1423 | yUpFirstDiffFlowPtRP=valueMCPtRP+errorMCPtRP; | |
1424 | } | |
1425 | } | |
1426 | } | |
1427 | for(Int_t i=nBinsDiffFlowPtRP+1;i<2*nBinsDiffFlowPtRP+1;i++) | |
1428 | { | |
1429 | //setting up the lower limit of the mesh: | |
1430 | valueMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(2*nBinsDiffFlowPtRP+1-i); | |
1431 | errorMCPtRP = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinError(2*nBinsDiffFlowPtRP+1-i); | |
1432 | if(!(errorMCPtRP==0.0)) | |
1433 | { | |
1434 | pMeshDiffFlowPtRP->SetPoint(countDiffFlowPtRP++,(2*nBinsDiffFlowPtRP-i+0.5)*binWidthPtRP+dPtMin,valueMCPtRP-errorMCPtRP); | |
1435 | } | |
1436 | } | |
1437 | //closing the mesh area: | |
1438 | pMeshDiffFlowPtRP->SetPoint(2*nNonEmptyBinsDiffFlowPtRP+1,xFirstDiffFlowPtRP,yUpFirstDiffFlowPtRP); | |
1439 | ||
1440 | //setting the mesh style and color: | |
dec3fab6 | 1441 | pMeshDiffFlowPtRP->SetFillStyle(meshStyleDiffFlowPtRP); |
1442 | pMeshDiffFlowPtRP->SetFillColor(meshColorDiffFlowPtRP); | |
1fca9c90 | 1443 | } |
7d02d77f | 1444 | |
1fca9c90 | 1445 | //cosmetics: Monte Carlo error bands for differential flow (Eta) |
1446 | TGraph* pMeshDiffFlowEtaRP = NULL; | |
1447 | if(mcepCommonHistRes) | |
1448 | { | |
7d02d77f | 1449 | Int_t nBinsDiffFlowEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetNbinsX(); |
1fca9c90 | 1450 | Double_t binWidthEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinWidth(1);//assuming that all bins have the same width |
e83922f2 | 1451 | |
7d02d77f | 1452 | //counting the non-empty bins: |
1453 | Int_t nNonEmptyBinsDiffFlowEtaRP=0; | |
1454 | for(Int_t i=1;i<nBinsDiffFlowEtaRP+1;i++) | |
1fca9c90 | 1455 | { |
7d02d77f | 1456 | if(!(mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinError(i)==0.0)) |
1457 | { | |
1458 | nNonEmptyBinsDiffFlowEtaRP++; | |
1459 | } | |
1fca9c90 | 1460 | } |
7d02d77f | 1461 | |
1462 | pMeshDiffFlowEtaRP = new TGraph(2*nNonEmptyBinsDiffFlowEtaRP+1); | |
1463 | ||
1464 | Double_t valueMCEtaRP=0.,errorMCEtaRP=0.; | |
1465 | Int_t countDiffFlowEtaRP=1; | |
1466 | Double_t xFirstDiffFlowEtaRP=0.,yUpFirstDiffFlowEtaRP=0.;//needed to close up the mesh | |
1467 | for(Int_t i=1;i<nBinsDiffFlowEtaRP+1;i++) | |
1fca9c90 | 1468 | { |
7d02d77f | 1469 | //setting up the upper limit of the mesh: |
1470 | valueMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinContent(i); | |
1471 | errorMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinError(i); | |
1472 | if(!(errorMCEtaRP==0.0)) | |
1473 | { | |
1474 | pMeshDiffFlowEtaRP->SetPoint(countDiffFlowEtaRP++,(i-0.5)*binWidthEtaRP+dEtaMin,valueMCEtaRP+errorMCEtaRP); | |
1475 | if(xFirstDiffFlowEtaRP==0.) | |
1476 | { | |
1477 | xFirstDiffFlowEtaRP=(i-0.5)*binWidthEtaRP+dEtaMin; | |
1478 | yUpFirstDiffFlowEtaRP=valueMCEtaRP+errorMCEtaRP; | |
1479 | } | |
1480 | } | |
1481 | } | |
1482 | for(Int_t i=nBinsDiffFlowEtaRP+1;i<2*nBinsDiffFlowEtaRP+1;i++) | |
1483 | { | |
1484 | //setting up the lower limit of the mesh: | |
1485 | valueMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinContent(2*nBinsDiffFlowEtaRP+1-i); | |
1486 | errorMCEtaRP = (mcepCommonHistRes->GetHistDiffFlowEtaRP())->GetBinError(2*nBinsDiffFlowEtaRP+1-i); | |
1487 | if(!(errorMCEtaRP==0.0)) | |
1488 | { | |
1489 | pMeshDiffFlowEtaRP->SetPoint(countDiffFlowEtaRP++,(2*nBinsDiffFlowEtaRP-i+0.5)*binWidthEtaRP+dEtaMin,valueMCEtaRP-errorMCEtaRP); | |
1490 | } | |
1491 | } | |
1492 | //closing the mesh area: | |
1493 | pMeshDiffFlowEtaRP->SetPoint(2*nNonEmptyBinsDiffFlowEtaRP+1,xFirstDiffFlowEtaRP,yUpFirstDiffFlowEtaRP); | |
1494 | ||
1495 | //setting the mesh style and color: | |
dec3fab6 | 1496 | pMeshDiffFlowEtaRP->SetFillStyle(meshStyleDiffFlowEtaRP); |
1497 | pMeshDiffFlowEtaRP->SetFillColor(meshColorDiffFlowEtaRP); | |
1fca9c90 | 1498 | } |
1499 | //---------------------------------------------------------------------------------- | |
1500 | ||
1501 | ||
7d02d77f | 1502 | |
1503 | ||
1fca9c90 | 1504 | //---------------------------------------------------------------------------------- |
1505 | //POI: | |
1506 | //cosmetics: Monte Carlo error bands for differential flow (Pt) | |
1507 | TGraph* pMeshDiffFlowPtPOI = NULL; | |
1508 | if(mcepCommonHistRes) | |
1509 | { | |
7d02d77f | 1510 | Int_t nBinsDiffFlowPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetNbinsX(); |
1fca9c90 | 1511 | Double_t binWidthPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinWidth(1);//assuming that all bins have the same width |
7d02d77f | 1512 | |
1513 | //counting the non-empty bins: | |
1514 | Int_t nNonEmptyBinsDiffFlowPtPOI=0; | |
1515 | for(Int_t i=1;i<nBinsDiffFlowPtPOI+1;i++) | |
1516 | { | |
1517 | if(!(mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinError(i)==0.0)) | |
1518 | { | |
1519 | nNonEmptyBinsDiffFlowPtPOI++; | |
1520 | } | |
1521 | } | |
1fca9c90 | 1522 | |
7d02d77f | 1523 | pMeshDiffFlowPtPOI = new TGraph(2*nNonEmptyBinsDiffFlowPtPOI+1); |
1fca9c90 | 1524 | |
7d02d77f | 1525 | Double_t valueMCPtPOI=0.,errorMCPtPOI=0.; |
1526 | Int_t countDiffFlowPtPOI=1; | |
1527 | Double_t xFirstDiffFlowPtPOI=0.,yUpFirstDiffFlowPtPOI=0.;//needed to close up the mesh | |
1528 | for(Int_t i=1;i<nBinsDiffFlowPtPOI+1;i++) | |
1fca9c90 | 1529 | { |
7d02d77f | 1530 | //setting up the upper limit of the mesh: |
1fca9c90 | 1531 | valueMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinContent(i); |
7d02d77f | 1532 | errorMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinError(i); |
1533 | if(!(errorMCPtPOI==0.0)) | |
1534 | { | |
1535 | pMeshDiffFlowPtPOI->SetPoint(countDiffFlowPtPOI++,(i-0.5)*binWidthPtPOI+dPtMin,valueMCPtPOI+errorMCPtPOI); | |
1536 | if(xFirstDiffFlowPtPOI==0.) | |
1537 | { | |
1538 | xFirstDiffFlowPtPOI=(i-0.5)*binWidthPtPOI+dPtMin; | |
1539 | yUpFirstDiffFlowPtPOI=valueMCPtPOI+errorMCPtPOI; | |
1540 | } | |
1541 | } | |
1542 | } | |
1543 | for(Int_t i=nBinsDiffFlowPtPOI+1;i<2*nBinsDiffFlowPtPOI+1;i++) | |
1544 | { | |
1545 | //setting up the lower limit of the mesh: | |
1546 | valueMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinContent(2*nBinsDiffFlowPtPOI+1-i); | |
1547 | errorMCPtPOI = (mcepCommonHistRes->GetHistDiffFlowPtPOI())->GetBinError(2*nBinsDiffFlowPtPOI+1-i); | |
1548 | if(!(errorMCPtPOI==0.0)) | |
1549 | { | |
1550 | pMeshDiffFlowPtPOI->SetPoint(countDiffFlowPtPOI++,(2*nBinsDiffFlowPtPOI-i+0.5)*binWidthPtPOI+dPtMin,valueMCPtPOI-errorMCPtPOI); | |
1551 | } | |
1552 | } | |
1553 | //closing the mesh area: | |
1554 | pMeshDiffFlowPtPOI->SetPoint(2*nNonEmptyBinsDiffFlowPtPOI+1,xFirstDiffFlowPtPOI,yUpFirstDiffFlowPtPOI); | |
1555 | ||
1556 | //setting the mesh style and color: | |
dec3fab6 | 1557 | pMeshDiffFlowPtPOI->SetFillStyle(meshStyleDiffFlowPtPOI); |
1558 | pMeshDiffFlowPtPOI->SetFillColor(meshColorDiffFlowPtPOI); | |
1fca9c90 | 1559 | } |
7d02d77f | 1560 | |
1fca9c90 | 1561 | //cosmetics: Monte Carlo error bands for differential flow (Eta) |
1562 | TGraph* pMeshDiffFlowEtaPOI = NULL; | |
1563 | if(mcepCommonHistRes) | |
1564 | { | |
7d02d77f | 1565 | Int_t nBinsDiffFlowEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetNbinsX(); |
1fca9c90 | 1566 | Double_t binWidthEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinWidth(1);//assuming that all bins have the same width |
7d02d77f | 1567 | |
1568 | //counting the non-empty bins: | |
1569 | Int_t nNonEmptyBinsDiffFlowEtaPOI=0; | |
1570 | for(Int_t i=1;i<nBinsDiffFlowEtaPOI+1;i++) | |
1571 | { | |
1572 | if(!(mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinError(i)==0.0)) | |
1573 | { | |
1574 | nNonEmptyBinsDiffFlowEtaPOI++; | |
1575 | } | |
1576 | } | |
1fca9c90 | 1577 | |
7d02d77f | 1578 | pMeshDiffFlowEtaPOI = new TGraph(2*nNonEmptyBinsDiffFlowEtaPOI+1); |
a88e5ac2 | 1579 | |
7d02d77f | 1580 | Double_t valueMCEtaPOI=0.,errorMCEtaPOI=0.; |
1581 | Int_t countDiffFlowEtaPOI=1; | |
1582 | Double_t xFirstDiffFlowEtaPOI=0.,yUpFirstDiffFlowEtaPOI=0.;//needed to close up the mesh | |
1583 | for(Int_t i=1;i<nBinsDiffFlowEtaPOI+1;i++) | |
1fca9c90 | 1584 | { |
7d02d77f | 1585 | //setting up the upper limit of the mesh: |
1fca9c90 | 1586 | valueMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinContent(i); |
7d02d77f | 1587 | errorMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinError(i); |
1588 | if(!(errorMCEtaPOI==0.0)) | |
1589 | { | |
1590 | pMeshDiffFlowEtaPOI->SetPoint(countDiffFlowEtaPOI++,(i-0.5)*binWidthEtaPOI+dEtaMin,valueMCEtaPOI+errorMCEtaPOI); | |
1591 | if(xFirstDiffFlowEtaPOI==0.) | |
1592 | { | |
1593 | xFirstDiffFlowEtaPOI=(i-0.5)*binWidthEtaPOI+dEtaMin; | |
1594 | yUpFirstDiffFlowEtaPOI=valueMCEtaPOI+errorMCEtaPOI; | |
1595 | } | |
1596 | } | |
1597 | } | |
1598 | for(Int_t i=nBinsDiffFlowEtaPOI+1;i<2*nBinsDiffFlowEtaPOI+1;i++) | |
1599 | { | |
1600 | //setting up the lower limit of the mesh: | |
1601 | valueMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinContent(2*nBinsDiffFlowEtaPOI+1-i); | |
1602 | errorMCEtaPOI = (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->GetBinError(2*nBinsDiffFlowEtaPOI+1-i); | |
1603 | if(!(errorMCEtaPOI==0.0)) | |
1604 | { | |
1605 | pMeshDiffFlowEtaPOI->SetPoint(countDiffFlowEtaPOI++,(2*nBinsDiffFlowEtaPOI-i+0.5)*binWidthEtaPOI+dEtaMin,valueMCEtaPOI-errorMCEtaPOI); | |
1606 | } | |
1607 | } | |
1608 | //closing the mesh area: | |
1609 | pMeshDiffFlowEtaPOI->SetPoint(2*nNonEmptyBinsDiffFlowEtaPOI+1,xFirstDiffFlowEtaPOI,yUpFirstDiffFlowEtaPOI); | |
1610 | ||
1611 | //setting the mesh style and color: | |
dec3fab6 | 1612 | pMeshDiffFlowEtaPOI->SetFillStyle(meshStyleDiffFlowEtaPOI); |
1613 | pMeshDiffFlowEtaPOI->SetFillColor(meshColorDiffFlowEtaPOI); | |
7d02d77f | 1614 | } |
1fca9c90 | 1615 | //---------------------------------------------------------------------------------- |
1616 | ||
70288839 | 1617 | //MCEP = Monte Carlo Event Plane |
1fca9c90 | 1618 | Double_t avMultDiffFlowMCRP=0.; |
1619 | Double_t nEvtsDiffFlowMCRP=0; | |
1620 | Double_t avMultDiffFlowMCPOI=0.; | |
1621 | Double_t nEvtsDiffFlowMCPOI=0; | |
a88e5ac2 | 1622 | if(fileMCEP) |
9d3cbee7 | 1623 | { |
70288839 | 1624 | if(mcepCommonHistRes) |
1625 | { | |
7d02d77f | 1626 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorMC); |
1627 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleMC); | |
1628 | (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorMC); | |
1629 | (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleMC); | |
1630 | (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorMC); | |
1631 | (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleMC); | |
1632 | (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorMC); | |
1633 | (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleMC); | |
a88e5ac2 | 1634 | } |
1635 | if(mcepCommonHist) | |
1636 | { | |
6a1a854d | 1637 | avMultDiffFlowMCRP = (mcepCommonHist->GetHistMultRP())->GetMean(); |
1638 | nEvtsDiffFlowMCRP = (mcepCommonHist->GetHistMultRP())->GetEntries(); | |
1639 | avMultDiffFlowMCPOI = (mcepCommonHist->GetHistMultPOI())->GetMean(); | |
1640 | nEvtsDiffFlowMCPOI = (mcepCommonHist->GetHistMultPOI())->GetEntries(); | |
70288839 | 1641 | } |
1642 | } | |
905bed86 | 1643 | |
1644 | //SP = Scalar Product | |
1645 | Double_t avMultDiffFlowSPRP=0.; | |
1646 | Double_t nEvtsDiffFlowSPRP=0; | |
1647 | Double_t avMultDiffFlowSPPOI=0.; | |
1648 | Double_t nEvtsDiffFlowSPPOI=0; | |
1649 | if(fileSP) | |
1650 | { | |
1651 | if(spCommonHistRes) | |
1652 | { | |
1653 | (spCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorSP); | |
1654 | (spCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleSP); | |
1655 | (spCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorSP); | |
1656 | (spCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleSP); | |
1657 | (spCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorSP); | |
1658 | (spCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleSP); | |
1659 | (spCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorSP); | |
1660 | (spCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleSP); | |
1661 | } | |
1662 | if(spCommonHist) | |
1663 | { | |
6a1a854d | 1664 | avMultDiffFlowSPRP = (spCommonHist->GetHistMultRP())->GetMean(); |
dec3fab6 | 1665 | nEvtsDiffFlowSPRP = (spCommonHist->GetHistMultRP())->GetEntries(); |
6a1a854d | 1666 | avMultDiffFlowSPPOI = (spCommonHist->GetHistMultPOI())->GetMean(); |
dec3fab6 | 1667 | nEvtsDiffFlowSPPOI = (spCommonHist->GetHistMultPOI())->GetEntries(); |
905bed86 | 1668 | } |
1669 | } | |
a88e5ac2 | 1670 | |
9d3cbee7 | 1671 | //GFC = Generating Function Cumulants |
1fca9c90 | 1672 | Double_t avMultDiffFlowGFC=0.;//to be removed |
1673 | Double_t nEvtsDiffFlowGFC=0.;//to be removed | |
1674 | Double_t avMultDiffFlowGFCRP=0.; | |
1675 | Double_t nEvtsDiffFlowGFCRP=0.; | |
1676 | Double_t avMultDiffFlowGFCPOI=0.; | |
1677 | Double_t nEvtsDiffFlowGFCPOI=0.; | |
a88e5ac2 | 1678 | if(fileGFC) |
9d3cbee7 | 1679 | { |
a88e5ac2 | 1680 | if(gfcCommonHistRes2) |
9d3cbee7 | 1681 | { |
7d02d77f | 1682 | (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC2); |
1683 | (gfcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC2); | |
1684 | (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC2); | |
1685 | (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleGFC2); | |
1686 | (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorGFC2); | |
1687 | (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleGFC2); | |
1688 | (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorGFC2); | |
1689 | (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleGFC2); | |
9d3cbee7 | 1690 | } |
a88e5ac2 | 1691 | if(gfcCommonHistRes4) |
9d3cbee7 | 1692 | { |
7d02d77f | 1693 | (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC4); |
1694 | (gfcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC4); | |
1695 | (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC4); | |
1696 | (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleGFC4); | |
1697 | (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorGFC4); | |
1698 | (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleGFC4); | |
1699 | (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorGFC4); | |
1700 | (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleGFC4); | |
1fca9c90 | 1701 | } |
1702 | if(gfcCommonHistRes6) | |
1703 | { | |
7d02d77f | 1704 | (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC6); |
1705 | (gfcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC6); | |
1706 | (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC6); | |
1707 | (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleGFC6); | |
1708 | (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorGFC6); | |
1709 | (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleGFC6); | |
1710 | (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorGFC6); | |
1711 | (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleGFC6); | |
1fca9c90 | 1712 | } |
1713 | if(gfcCommonHistRes8) | |
1714 | { | |
7d02d77f | 1715 | (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorGFC8); |
1716 | (gfcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleGFC8); | |
1717 | (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorGFC8); | |
1718 | (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleGFC8); | |
1719 | (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorGFC8); | |
1720 | (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleGFC8); | |
1721 | (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorGFC8); | |
1722 | (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleGFC8); | |
9d3cbee7 | 1723 | } |
a88e5ac2 | 1724 | if(gfcCommonHist) |
1725 | { | |
6a1a854d | 1726 | avMultDiffFlowGFCRP = (gfcCommonHist->GetHistMultRP())->GetMean(); |
1727 | nEvtsDiffFlowGFCRP = (gfcCommonHist->GetHistMultRP())->GetEntries(); | |
1728 | avMultDiffFlowGFCPOI = (gfcCommonHist->GetHistMultPOI())->GetMean(); | |
1729 | nEvtsDiffFlowGFCPOI = (gfcCommonHist->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 1730 | } |
9d3cbee7 | 1731 | } |
a88e5ac2 | 1732 | |
9d3cbee7 | 1733 | //QC = Q-cumulants |
1fca9c90 | 1734 | Double_t avMultDiffFlowQC2RP=0.; |
1735 | Double_t nEvtsDiffFlowQC2RP=0.; | |
1736 | Double_t avMultDiffFlowQC2POI=0.; | |
1737 | Double_t nEvtsDiffFlowQC2POI=0.; | |
1738 | Double_t avMultDiffFlowQC4RP=0.; | |
1739 | Double_t nEvtsDiffFlowQC4RP=0.; | |
1740 | Double_t avMultDiffFlowQC4POI=0.; | |
1741 | Double_t nEvtsDiffFlowQC4POI=0.; | |
e83922f2 | 1742 | Double_t avMultDiffFlowQC6RP=0.; |
1743 | Double_t nEvtsDiffFlowQC6RP=0.; | |
1744 | Double_t avMultDiffFlowQC6POI=0.; | |
1745 | Double_t nEvtsDiffFlowQC6POI=0.; | |
1746 | Double_t avMultDiffFlowQC8RP=0.; | |
1747 | Double_t nEvtsDiffFlowQC8RP=0.; | |
1748 | Double_t avMultDiffFlowQC8POI=0.; | |
1749 | Double_t nEvtsDiffFlowQC8POI=0.; | |
1fca9c90 | 1750 | |
a88e5ac2 | 1751 | if(fileQC) |
9d3cbee7 | 1752 | { |
a88e5ac2 | 1753 | //QC{2} |
70288839 | 1754 | if(qcCommonHistRes2) |
9d3cbee7 | 1755 | { |
7d02d77f | 1756 | (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC2); |
1757 | (qcCommonHistRes2->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC2); | |
1758 | (qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC2); | |
1759 | (qcCommonHistRes2->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleQC2); | |
1760 | (qcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorQC2); | |
1761 | (qcCommonHistRes2->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleQC2); | |
1762 | (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorQC2); | |
1763 | (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleQC2); | |
9d3cbee7 | 1764 | } |
a88e5ac2 | 1765 | if(qcCommonHist2) |
1766 | { | |
6a1a854d | 1767 | avMultDiffFlowQC2RP = (qcCommonHist2->GetHistMultRP())->GetMean(); |
1768 | nEvtsDiffFlowQC2RP = (qcCommonHist2->GetHistMultRP())->GetEntries(); | |
1769 | avMultDiffFlowQC2POI = (qcCommonHist2->GetHistMultPOI())->GetMean(); | |
1770 | nEvtsDiffFlowQC2POI = (qcCommonHist2->GetHistMultPOI())->GetEntries(); | |
a88e5ac2 | 1771 | } |
1772 | //QC{4} | |
70288839 | 1773 | if(qcCommonHistRes4) |
9d3cbee7 | 1774 | { |
7d02d77f | 1775 | (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC4); |
1776 | (qcCommonHistRes4->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC4); | |
1777 | (qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC4); | |
1778 | (qcCommonHistRes4->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleQC4); | |
1779 | (qcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorQC4); | |
1780 | (qcCommonHistRes4->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleQC4); | |
1781 | (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorQC4); | |
1782 | (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleQC4); | |
a88e5ac2 | 1783 | } |
1784 | if(qcCommonHist4) | |
1785 | { | |
6a1a854d | 1786 | avMultDiffFlowQC4RP = (qcCommonHist4->GetHistMultRP())->GetMean(); |
1787 | nEvtsDiffFlowQC4RP = (qcCommonHist4->GetHistMultRP())->GetEntries(); | |
1788 | avMultDiffFlowQC4POI = (qcCommonHist4->GetHistMultPOI())->GetMean(); | |
1789 | nEvtsDiffFlowQC4POI = (qcCommonHist4->GetHistMultPOI())->GetEntries(); | |
9d3cbee7 | 1790 | } |
e83922f2 | 1791 | //QC{6} |
1792 | if(qcCommonHistRes6) | |
1793 | { | |
7d02d77f | 1794 | (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC6); |
1795 | (qcCommonHistRes6->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC6); | |
1796 | (qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC6); | |
1797 | (qcCommonHistRes6->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleQC6); | |
1798 | (qcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorQC6); | |
1799 | (qcCommonHistRes6->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleQC6); | |
1800 | (qcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorQC6); | |
1801 | (qcCommonHistRes6->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleQC6); | |
e83922f2 | 1802 | } |
1803 | if(qcCommonHist6) | |
1804 | { | |
6a1a854d | 1805 | avMultDiffFlowQC6RP = (qcCommonHist6->GetHistMultRP())->GetMean(); |
1806 | nEvtsDiffFlowQC6RP = (qcCommonHist6->GetHistMultRP())->GetEntries(); | |
1807 | avMultDiffFlowQC6POI = (qcCommonHist6->GetHistMultPOI())->GetMean(); | |
1808 | nEvtsDiffFlowQC6POI = (qcCommonHist6->GetHistMultPOI())->GetEntries(); | |
e83922f2 | 1809 | } |
1810 | //QC{8} | |
1811 | if(qcCommonHistRes8) | |
1812 | { | |
7d02d77f | 1813 | (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorQC8); |
1814 | (qcCommonHistRes8->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleQC8); | |
1815 | (qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorQC8); | |
1816 | (qcCommonHistRes8->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleQC8); | |
1817 | (qcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorQC8); | |
1818 | (qcCommonHistRes8->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleQC8); | |
1819 | (qcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorQC8); | |
1820 | (qcCommonHistRes8->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleQC8); | |
e83922f2 | 1821 | } |
1822 | if(qcCommonHist8) | |
1823 | { | |
6a1a854d | 1824 | avMultDiffFlowQC8RP = (qcCommonHist8->GetHistMultRP())->GetMean(); |
1825 | nEvtsDiffFlowQC8RP = (qcCommonHist8->GetHistMultRP())->GetEntries(); | |
1826 | avMultDiffFlowQC8POI = (qcCommonHist8->GetHistMultPOI())->GetMean(); | |
1827 | nEvtsDiffFlowQC8POI = (qcCommonHist8->GetHistMultPOI())->GetEntries(); | |
e83922f2 | 1828 | } |
9d3cbee7 | 1829 | } |
70288839 | 1830 | |
a88e5ac2 | 1831 | //LYZ2 = Lee-Yang Zeros (2nd run) |
1fca9c90 | 1832 | Double_t avMultDiffFlowLYZ2RP=0.; |
1833 | Double_t nEvtsDiffFlowLYZ2RP=0; | |
1834 | Double_t avMultDiffFlowLYZ2POI=0.; | |
1835 | Double_t nEvtsDiffFlowLYZ2POI=0; | |
a88e5ac2 | 1836 | if(fileLYZ2) |
9d3cbee7 | 1837 | { |
a88e5ac2 | 1838 | if(lyz2CommonHistRes) |
9d3cbee7 | 1839 | { |
7d02d77f | 1840 | (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZ2); |
1841 | (lyz2CommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZ2); | |
1842 | (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZ2); | |
1843 | (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleLYZ2); | |
1844 | (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorLYZ2); | |
1845 | (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleLYZ2); | |
1846 | (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorLYZ2); | |
1847 | (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleLYZ2); | |
a88e5ac2 | 1848 | } |
1849 | if(lyz2CommonHist) | |
9d3cbee7 | 1850 | { |
6a1a854d | 1851 | avMultDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultRP())->GetMean(); |
1852 | nEvtsDiffFlowLYZ2RP = (lyz2CommonHist->GetHistMultRP())->GetEntries(); | |
1853 | avMultDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultPOI())->GetMean(); | |
1854 | nEvtsDiffFlowLYZ2POI = (lyz2CommonHist->GetHistMultPOI())->GetEntries(); | |
9d3cbee7 | 1855 | } |
1856 | } | |
70288839 | 1857 | |
e83922f2 | 1858 | //LYZEP = Lee-Yang Zeros Event Plane |
e83922f2 | 1859 | Double_t avMultDiffFlowLYZEPRP=0.; |
1860 | Double_t nEvtsDiffFlowLYZEPRP=0; | |
1861 | Double_t avMultDiffFlowLYZEPPOI=0.; | |
1862 | Double_t nEvtsDiffFlowLYZEPPOI=0; | |
1863 | if(fileLYZEP) | |
1864 | { | |
1865 | if(lyzepCommonHistRes) | |
1866 | { | |
7d02d77f | 1867 | (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerColor(markerColorLYZEP); |
1868 | (lyzepCommonHistRes->GetHistDiffFlowPtRP())->SetMarkerStyle(markerStyleLYZEP); | |
1869 | (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerColor(markerColorLYZEP); | |
1870 | (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->SetMarkerStyle(markerStyleLYZEP); | |
1871 | (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerColor(markerColorLYZEP); | |
1872 | (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->SetMarkerStyle(markerStyleLYZEP); | |
1873 | (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerColor(markerColorLYZEP); | |
1874 | (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->SetMarkerStyle(markerStyleLYZEP); | |
e83922f2 | 1875 | } |
1876 | if(lyzepCommonHist) | |
1877 | { | |
6a1a854d | 1878 | avMultDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultRP())->GetMean(); |
1879 | nEvtsDiffFlowLYZEPRP = (lyzepCommonHist->GetHistMultRP())->GetEntries(); | |
1880 | avMultDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultPOI())->GetMean(); | |
1881 | nEvtsDiffFlowLYZEPPOI = (lyzepCommonHist->GetHistMultPOI())->GetEntries(); | |
e83922f2 | 1882 | } |
1883 | } | |
1884 | ||
1fca9c90 | 1885 | |
1fca9c90 | 1886 | //---------------------------------------------------------------------------------- |
1887 | //final drawing for differential flow (Pt, RP): | |
e83922f2 | 1888 | Bool_t plotLYZEPPtRP = kTRUE; |
1889 | if(plotDiffFlowPtRP) | |
1890 | { | |
1891 | TCanvas* diffFlowPtAllCanvasRP = new TCanvas("Differential Flow (Pt) of RP","Differential Flow (Pt) of RP ",1000,600); | |
1892 | ||
1893 | diffFlowPtAllCanvasRP->Divide(2,1); | |
1fca9c90 | 1894 | |
e83922f2 | 1895 | //1st pad is for plot: |
1896 | (diffFlowPtAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
6f62c2a7 | 1897 | |
e83922f2 | 1898 | if(styleHistPt) |
1899 | { | |
1900 | (styleHistPt->GetYaxis())->SetRangeUser(-0.3,1.0); | |
1901 | styleHistPt->Draw(); | |
1902 | } | |
1903 | if(pMeshDiffFlowPtRP) | |
1904 | { | |
1905 | pMeshDiffFlowPtRP->Draw("LFSAME"); | |
1906 | } | |
1fca9c90 | 1907 | |
e83922f2 | 1908 | //MC |
1909 | if(plotMCPtRP && mcepCommonHistRes) | |
1910 | { | |
1911 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1912 | } | |
905bed86 | 1913 | //SP |
1914 | if(plotSPPtRP && spCommonHistRes) | |
1915 | { | |
1916 | (spCommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1917 | } | |
e83922f2 | 1918 | //GFC |
7d02d77f | 1919 | if(plotGFC2PtRP && gfcCommonHistRes2)Pt |
e83922f2 | 1920 | { |
1921 | (gfcCommonHistRes2->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1922 | } | |
1923 | if(plotGFC4PtRP && gfcCommonHistRes4) | |
1924 | { | |
1925 | (gfcCommonHistRes4->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1926 | } | |
1927 | if(plotGFC6PtRP && gfcCommonHistRes6) | |
1928 | { | |
1929 | (gfcCommonHistRes6->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1930 | } | |
1931 | if(plotGFC8PtRP && gfcCommonHistRes8) | |
1932 | { | |
1933 | (gfcCommonHistRes8->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1934 | } | |
1935 | //QC | |
1936 | if(plotQC2PtRP && qcCommonHistRes2) | |
1937 | { | |
1938 | (qcCommonHistRes2->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1939 | } | |
1940 | if(plotQC4PtRP && qcCommonHistRes4) | |
1941 | { | |
1942 | (qcCommonHistRes4->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1943 | } | |
1944 | if(plotQC6PtRP && qcCommonHistRes6) | |
1945 | { | |
7d02d77f | 1946 | //(qcCommonHistRes6->GetHistDiffFlowPtRP())->Draw("E1PSAME"); |
e83922f2 | 1947 | } |
1948 | if(plotQC8PtRP && qcCommonHistRes8) | |
1949 | { | |
7d02d77f | 1950 | //(qcCommonHistRes8->GetHistDiffFlowPtRP())->Draw("E1PSAME"); |
e83922f2 | 1951 | } |
1952 | //LYZ2 | |
1953 | if(plotLYZ2PtRP && lyz2CommonHistRes) | |
1954 | { | |
1955 | (lyz2CommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1956 | } | |
1957 | //LYZEP | |
1958 | if(plotLYZEPPtRP && lyzepCommonHistRes) | |
1959 | { | |
1960 | (lyzepCommonHistRes->GetHistDiffFlowPtRP())->Draw("E1PSAME"); | |
1961 | } | |
1fca9c90 | 1962 | |
e83922f2 | 1963 | //2nd pad is for legend: |
dec3fab6 | 1964 | (diffFlowPtAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0); |
1965 | ||
1966 | TLegend* legendDiffFlowPtRP = new TLegend(0.02,0.12,0.97,0.70); | |
e83922f2 | 1967 | legendDiffFlowPtRP->SetTextFont(72); |
1968 | legendDiffFlowPtRP->SetTextSize(0.06); | |
1969 | ||
7d02d77f | 1970 | //legend's entries:Pt |
e83922f2 | 1971 | TString *entryDiffMCPtRP = new TString("MC ....... "); |
905bed86 | 1972 | TString *entryDiffSPPtRP = new TString("SP ....... "); |
e83922f2 | 1973 | TString *entryDiffGFC2PtRP = new TString("GFC{2} ... "); |
1974 | TString *entryDiffGFC4PtRP = new TString("GFC{4} ... "); | |
1975 | TString *entryDiffGFC6PtRP = new TString("GFC{6} ... "); | |
1976 | TString *entryDiffGFC8PtRP = new TString("GFC{8} ... "); | |
1977 | TString *entryDiffQC2PtRP = new TString("QC{2} .... "); | |
1978 | TString *entryDiffQC4PtRP = new TString("QC{4} .... "); | |
1979 | TString *entryDiffQC6PtRP = new TString("QC{6} .... "); | |
1980 | TString *entryDiffQC8PtRP = new TString("QC{8} .... "); | |
1981 | TString *entryDiffLYZ2PtRP = new TString("LYZ ...... "); | |
1982 | TString *entryDiffLYZEPPtRP = new TString("LYZEP ... "); | |
1983 | ||
1984 | //MC | |
1985 | if(mcepCommonHistRes) | |
1986 | { | |
dec3fab6 | 1987 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillStyle(meshStyleDiffFlowPtRP); |
1988 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillColor(meshColorDiffFlowPtRP); | |
e83922f2 | 1989 | entryDiffMCPtRP->Append("M = "); |
1990 | (*entryDiffMCPtRP)+=(Long_t)avMultDiffFlowMCRP; | |
1991 | entryDiffMCPtRP->Append(", N = "); | |
1992 | (*entryDiffMCPtRP)+=(Long_t)nEvtsDiffFlowMCRP; | |
1993 | legendDiffFlowPtRP->AddEntry(mcepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffMCPtRP->Data(),"f"); | |
1994 | } | |
905bed86 | 1995 | |
1996 | //SP | |
1997 | if(spCommonHistRes) | |
1998 | { | |
905bed86 | 1999 | entryDiffSPPtRP->Append("M = "); |
2000 | (*entryDiffSPPtRP)+=(Long_t)avMultDiffFlowSPRP; | |
2001 | entryDiffSPPtRP->Append(", N = "); | |
2002 | (*entryDiffSPPtRP)+=(Long_t)nEvtsDiffFlowSPRP; | |
2003 | legendDiffFlowPtRP->AddEntry(spCommonHistRes->GetHistDiffFlowPtRP(),entryDiffSPPtRP->Data(),"p"); | |
2004 | } | |
1fca9c90 | 2005 | |
e83922f2 | 2006 | //GFC |
2007 | if(plotGFC2PtRP && gfcCommonHistRes2) | |
2008 | { | |
2009 | entryDiffGFC2PtRP->Append("M = "); | |
2010 | (*entryDiffGFC2PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2011 | entryDiffGFC2PtRP->Append(", N = "); | |
2012 | (*entryDiffGFC2PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2013 | legendDiffFlowPtRP->AddEntry(gfcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffGFC2PtRP->Data(),"p"); | |
2014 | } | |
2015 | if(plotGFC4PtRP && gfcCommonHistRes4) | |
2016 | { | |
2017 | entryDiffGFC4PtRP->Append("M = "); | |
2018 | (*entryDiffGFC4PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2019 | entryDiffGFC4PtRP->Append(", N = "); | |
2020 | (*entryDiffGFC4PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2021 | legendDiffFlowPtRP->AddEntry(gfcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffGFC4PtRP->Data(),"p"); | |
2022 | } | |
2023 | if(plotGFC6PtRP && gfcCommonHistRes6) | |
2024 | { | |
2025 | entryDiffGFC6PtRP->Append("M = "); | |
2026 | (*entryDiffGFC6PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2027 | entryDiffGFC6PtRP->Append(", N = "); | |
2028 | (*entryDiffGFC6PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2029 | legendDiffFlowPtRP->AddEntry(gfcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffGFC6PtRP->Data(),"p"); | |
2030 | } | |
2031 | if(plotGFC8PtRP && gfcCommonHistRes8) | |
2032 | { | |
2033 | entryDiffGFC8PtRP->Append("M = "); | |
2034 | (*entryDiffGFC8PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2035 | entryDiffGFC8PtRP->Append(", N = "); | |
2036 | (*entryDiffGFC8PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2037 | legendDiffFlowPtRP->AddEntry(gfcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffGFC8PtRP->Data(),"p"); | |
2038 | } | |
2039 | ||
2040 | //QC | |
2041 | if(plotQC2PtRP && qcCommonHistRes2) | |
2042 | { | |
2043 | entryDiffQC2PtRP->Append("M = "); | |
2044 | (*entryDiffQC2PtRP)+=(Long_t)avMultDiffFlowQC2RP; | |
2045 | entryDiffQC2PtRP->Append(", N = "); | |
2046 | (*entryDiffQC2PtRP)+=(Long_t)nEvtsDiffFlowQC2RP; | |
2047 | legendDiffFlowPtRP->AddEntry(qcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffQC2PtRP->Data(),"p"); | |
2048 | } | |
2049 | if(plotQC4PtRP && qcCommonHistRes4) | |
2050 | { | |
2051 | entryDiffQC4PtRP->Append("M = "); | |
2052 | (*entryDiffQC4PtRP)+=(Long_t)avMultDiffFlowQC4RP; | |
2053 | entryDiffQC4PtRP->Append(", N = "); | |
2054 | (*entryDiffQC4PtRP)+=(Long_t)nEvtsDiffFlowQC4RP; | |
2055 | legendDiffFlowPtRP->AddEntry(qcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffQC4PtRP->Data(),"p"); | |
2056 | } | |
2057 | if(plotQC6PtRP && qcCommonHistRes6) | |
2058 | { | |
2059 | entryDiffQC6PtRP->Append("M = "); | |
2060 | (*entryDiffQC6PtRP)+=(Long_t)avMultDiffFlowQC6RP; | |
2061 | entryDiffQC6PtRP->Append(", N = "); | |
2062 | (*entryDiffQC6PtRP)+=(Long_t)nEvtsDiffFlowQC6RP; | |
2063 | legendDiffFlowPtRP->AddEntry(qcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffQC6PtRP->Data(),"p"); | |
2064 | } | |
2065 | if(plotQC8PtRP && qcCommonHistRes8) | |
2066 | { | |
2067 | entryDiffQC8PtRP->Append("M = "); | |
2068 | (*entryDiffQC8PtRP)+=(Long_t)avMultDiffFlowQC8RP; | |
2069 | entryDiffQC8PtRP->Append(", N = "); | |
2070 | (*entryDiffQC8PtRP)+=(Long_t)nEvtsDiffFlowQC8RP; | |
2071 | legendDiffFlowPtRP->AddEntry(qcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffQC8PtRP->Data(),"p"); | |
2072 | } | |
2073 | ||
2074 | //LYZ2 | |
2075 | if(plotLYZ2PtRP && lyz2CommonHistRes) | |
2076 | { | |
2077 | entryDiffLYZ2PtRP->Append("M = "); | |
2078 | (*entryDiffLYZ2PtRP)+=(Long_t)avMultDiffFlowLYZ2RP; | |
2079 | entryDiffLYZ2PtRP->Append(", N = "); | |
2080 | (*entryDiffLYZ2PtRP)+=(Long_t)nEvtsDiffFlowLYZ2RP; | |
2081 | legendDiffFlowPtRP->AddEntry(lyz2CommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZ2PtRP->Data(),"p"); | |
2082 | } | |
2083 | ||
2084 | //LYZEP | |
2085 | if(plotLYZEPPtRP && lyzepCommonHistRes) | |
2086 | { | |
2087 | entryDiffLYZEPPtRP->Append("M = "); | |
2088 | (*entryDiffLYZEPPtRP)+=(Long_t)avMultDiffFlowLYZEPRP; | |
2089 | entryDiffLYZEPPtRP->Append(", N = "); | |
2090 | (*entryDiffLYZEPPtRP)+=(Long_t)nEvtsDiffFlowLYZEPRP; | |
2091 | legendDiffFlowPtRP->AddEntry(lyzepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZEPPtRP->Data(),"p"); | |
2092 | } | |
1fca9c90 | 2093 | |
dec3fab6 | 2094 | //drawing finally the legend in the 2nd pad: |
2095 | if(textDefault) textDefault->Draw(); | |
2096 | ||
e83922f2 | 2097 | if(legendDiffFlowPtRP) |
2098 | { | |
2099 | legendDiffFlowPtRP->SetMargin(0.15); | |
2100 | legendDiffFlowPtRP->Draw(); | |
2101 | } | |
2102 | }// end of if(plotDiffFlowPtRP) | |
1fca9c90 | 2103 | //---------------------------------------------------------------------------------- |
2104 | ||
2105 | ||
2106 | //---------------------------------------------------------------------------------- | |
2107 | //final drawing for differential flow (Eta, RP): | |
e83922f2 | 2108 | if(plotDiffFlowEtaRP) |
2109 | { | |
2110 | TCanvas* diffFlowEtaAllCanvasRP = new TCanvas("Differential Flow (Eta) of RP","Differential Flow (Eta) of RP ",1000,600); | |
2111 | ||
2112 | diffFlowEtaAllCanvasRP->Divide(2,1); | |
2113 | ||
2114 | //1st pad is for plot: | |
2115 | (diffFlowEtaAllCanvasRP->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
2116 | ||
2117 | if(styleHistEta) | |
2118 | { | |
2119 | (styleHistEta->GetYaxis())->SetRangeUser(-0.3,1.0); | |
2120 | styleHistEta->Draw(); | |
2121 | } | |
2122 | if(pMeshDiffFlowEtaRP) | |
2123 | { | |
2124 | pMeshDiffFlowEtaRP->Draw("LFSAME"); | |
2125 | } | |
1fca9c90 | 2126 | |
e83922f2 | 2127 | //MC |
2128 | if(plotMCEtaRP && mcepCommonHistRes) | |
2129 | { | |
2130 | (mcepCommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2131 | } | |
905bed86 | 2132 | //SP |
2133 | if(plotSPEtaRP && spCommonHistRes) | |
2134 | { | |
2135 | (spCommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2136 | } | |
e83922f2 | 2137 | //GFC |
2138 | if(plotGFC2EtaRP && gfcCommonHistRes2) | |
2139 | { | |
2140 | (gfcCommonHistRes2->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2141 | } | |
2142 | if(plotGFC4EtaRP && gfcCommonHistRes4) | |
2143 | { | |
2144 | (gfcCommonHistRes4->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2145 | } | |
2146 | if(plotGFC6EtaRP && gfcCommonHistRes6) | |
2147 | { | |
2148 | (gfcCommonHistRes6->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2149 | } | |
2150 | if(plotGFC8EtaRP && gfcCommonHistRes8) | |
2151 | { | |
2152 | (gfcCommonHistRes8->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2153 | } | |
2154 | //QC | |
2155 | if(plotQC2EtaRP && qcCommonHistRes2) | |
2156 | { | |
2157 | (qcCommonHistRes2->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2158 | } | |
2159 | if(plotQC4EtaRP && qcCommonHistRes4) | |
2160 | { | |
2161 | (qcCommonHistRes4->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2162 | } | |
2163 | if(plotQC6EtaRP && qcCommonHistRes6) | |
2164 | { | |
7d02d77f | 2165 | //(qcCommonHistRes6->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); |
e83922f2 | 2166 | } |
2167 | if(plotQC8EtaRP && qcCommonHistRes8) | |
2168 | { | |
7d02d77f | 2169 | //(qcCommonHistRes8->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); |
e83922f2 | 2170 | } |
2171 | //LYZ2 | |
2172 | if(plotLYZ2EtaRP && lyz2CommonHistRes) | |
2173 | { | |
2174 | (lyz2CommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2175 | } | |
2176 | //LYZEP | |
2177 | if(plotLYZEPEtaRP && lyzepCommonHistRes) | |
2178 | { | |
2179 | (lyzepCommonHistRes->GetHistDiffFlowEtaRP())->Draw("E1PSAME"); | |
2180 | } | |
1fca9c90 | 2181 | |
e83922f2 | 2182 | //2nd pad is for legend: |
2183 | (diffFlowEtaAllCanvasRP->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
2184 | ||
dec3fab6 | 2185 | TLegend* legendDiffFlowEtaRP = new TLegend(0.02,0.12,0.97,0.70); |
e83922f2 | 2186 | legendDiffFlowEtaRP->SetTextFont(72); |
2187 | legendDiffFlowEtaRP->SetTextSize(0.06); | |
2188 | ||
2189 | //legend's entries: | |
2190 | TString *entryDiffMCEtaRP = new TString("MC ....... "); | |
905bed86 | 2191 | TString *entryDiffSPEtaRP = new TString("SP ....... "); |
e83922f2 | 2192 | TString *entryDiffGFC2EtaRP = new TString("GFC{2} ... "); |
2193 | TString *entryDiffGFC4EtaRP = new TString("GFC{4} ... "); | |
2194 | TString *entryDiffGFC6EtaRP = new TString("GFC{6} ... "); | |
2195 | TString *entryDiffGFC8EtaRP = new TString("GFC{8} ... "); | |
2196 | TString *entryDiffQC2EtaRP = new TString("QC{2} .... "); | |
2197 | TString *entryDiffQC4EtaRP = new TString("QC{4} .... "); | |
2198 | TString *entryDiffQC6EtaRP = new TString("QC{6} .... "); | |
2199 | TString *entryDiffQC8EtaRP = new TString("QC{8} .... "); | |
2200 | TString *entryDiffLYZ2EtaRP = new TString("LYZ ...... "); | |
2201 | TString *entryDiffLYZEPEtaRP = new TString("LYZEP ... "); | |
2202 | ||
2203 | //MC | |
2204 | if(mcepCommonHistRes) | |
2205 | { | |
dec3fab6 | 2206 | (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetFillStyle(meshStyleDiffFlowEtaRP); |
2207 | (mcepCommonHistRes->GetHistDiffFlowEtaRP())->SetFillColor(meshColorDiffFlowEtaRP); | |
e83922f2 | 2208 | entryDiffMCEtaRP->Append("M = "); |
2209 | (*entryDiffMCEtaRP)+=(Long_t)avMultDiffFlowMCRP; | |
2210 | entryDiffMCEtaRP->Append(", N = "); | |
2211 | (*entryDiffMCEtaRP)+=(Long_t)nEvtsDiffFlowMCRP; | |
2212 | legendDiffFlowEtaRP->AddEntry(mcepCommonHistRes->GetHistDiffFlowEtaRP(),entryDiffMCEtaRP->Data(),"f"); | |
2213 | } | |
905bed86 | 2214 | |
2215 | //SP | |
2216 | if(spCommonHistRes) | |
2217 | { | |
905bed86 | 2218 | entryDiffSPEtaRP->Append("M = "); |
2219 | (*entryDiffSPEtaRP)+=(Long_t)avMultDiffFlowSPRP; | |
2220 | entryDiffSPEtaRP->Append(", N = "); | |
2221 | (*entryDiffSPEtaRP)+=(Long_t)nEvtsDiffFlowSPRP; | |
2222 | legendDiffFlowEtaRP->AddEntry(spCommonHistRes->GetHistDiffFlowEtaRP(),entryDiffSPEtaRP->Data(),"p"); | |
2223 | } | |
1fca9c90 | 2224 | |
e83922f2 | 2225 | //GFC |
2226 | if(plotGFC2EtaRP && gfcCommonHistRes2) | |
2227 | { | |
2228 | entryDiffGFC2EtaRP->Append("M = "); | |
2229 | (*entryDiffGFC2EtaRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2230 | entryDiffGFC2EtaRP->Append(", N = "); | |
2231 | (*entryDiffGFC2EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2232 | legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes2->GetHistDiffFlowEtaRP(),entryDiffGFC2EtaRP->Data(),"p"); | |
2233 | } | |
2234 | if(plotGFC4EtaRP && gfcCommonHistRes4) | |
2235 | { | |
2236 | entryDiffGFC4EtaRP->Append("M = "); | |
2237 | (*entryDiffGFC4EtaRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2238 | entryDiffGFC4EtaRP->Append(", N = "); | |
2239 | (*entryDiffGFC4EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2240 | legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes4->GetHistDiffFlowEtaRP(),entryDiffGFC4EtaRP->Data(),"p"); | |
2241 | } | |
2242 | if(plotGFC6EtaRP && gfcCommonHistRes6) | |
2243 | { | |
2244 | entryDiffGFC6EtaRP->Append("M = "); | |
2245 | (*entryDiffGFC6EtaRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2246 | entryDiffGFC6EtaRP->Append(", N = "); | |
2247 | (*entryDiffGFC6EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2248 | legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes6->GetHistDiffFlowEtaRP(),entryDiffGFC6EtaRP->Data(),"p"); | |
2249 | } | |
2250 | if(plotGFC8EtaRP && gfcCommonHistRes8) | |
2251 | { | |
2252 | entryDiffGFC8EtaRP->Append("M = "); | |
2253 | (*entryDiffGFC8EtaRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2254 | entryDiffGFC8EtaRP->Append(", N = "); | |
2255 | (*entryDiffGFC8EtaRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2256 | legendDiffFlowEtaRP->AddEntry(gfcCommonHistRes8->GetHistDiffFlowEtaRP(),entryDiffGFC8EtaRP->Data(),"p"); | |
2257 | } | |
2258 | ||
2259 | //QC | |
2260 | if(plotQC2EtaRP && qcCommonHistRes2) | |
2261 | { | |
2262 | entryDiffQC2EtaRP->Append("M = "); | |
2263 | (*entryDiffQC2EtaRP)+=(Long_t)avMultDiffFlowQC2RP; | |
2264 | entryDiffQC2EtaRP->Append(", N = "); | |
2265 | (*entryDiffQC2EtaRP)+=(Long_t)nEvtsDiffFlowQC2RP; | |
2266 | legendDiffFlowEtaRP->AddEntry(qcCommonHistRes2->GetHistDiffFlowEtaRP(),entryDiffQC2EtaRP->Data(),"p"); | |
2267 | } | |
2268 | if(plotQC4EtaRP && qcCommonHistRes4) | |
2269 | { | |
2270 | entryDiffQC4EtaRP->Append("M = "); | |
2271 | (*entryDiffQC4EtaRP)+=(Long_t)avMultDiffFlowQC4RP; | |
2272 | entryDiffQC4EtaRP->Append(", N = "); | |
2273 | (*entryDiffQC4EtaRP)+=(Long_t)nEvtsDiffFlowQC4RP; | |
2274 | legendDiffFlowEtaRP->AddEntry(qcCommonHistRes4->GetHistDiffFlowEtaRP(),entryDiffQC4EtaRP->Data(),"p"); | |
2275 | } | |
2276 | if(plotQC6EtaRP && qcCommonHistRes6) | |
2277 | { | |
2278 | entryDiffQC6EtaRP->Append("M = "); | |
2279 | (*entryDiffQC6EtaRP)+=(Long_t)avMultDiffFlowQC6RP; | |
2280 | entryDiffQC6EtaRP->Append(", N = "); | |
2281 | (*entryDiffQC6EtaRP)+=(Long_t)nEvtsDiffFlowQC6RP; | |
2282 | legendDiffFlowEtaRP->AddEntry(qcCommonHistRes6->GetHistDiffFlowEtaRP(),entryDiffQC6EtaRP->Data(),"p"); | |
2283 | } | |
2284 | if(plotQC8EtaRP && qcCommonHistRes8) | |
2285 | { | |
2286 | entryDiffQC8EtaRP->Append("M = "); | |
2287 | (*entryDiffQC8EtaRP)+=(Long_t)avMultDiffFlowQC8RP; | |
2288 | entryDiffQC8EtaRP->Append(", N = "); | |
2289 | (*entryDiffQC8EtaRP)+=(Long_t)nEvtsDiffFlowQC8RP; | |
2290 | legendDiffFlowEtaRP->AddEntry(qcCommonHistRes8->GetHistDiffFlowEtaRP(),entryDiffQC8EtaRP->Data(),"p"); | |
2291 | } | |
1fca9c90 | 2292 | |
e83922f2 | 2293 | //LYZ2 |
2294 | if(plotLYZ2EtaRP && lyz2CommonHistRes) | |
2295 | { | |
2296 | entryDiffLYZ2EtaRP->Append("M = "); | |
2297 | (*entryDiffLYZ2EtaRP)+=(Long_t)avMultDiffFlowLYZ2RP; | |
2298 | entryDiffLYZ2EtaRP->Append(", N = "); | |
2299 | (*entryDiffLYZ2EtaRP)+=(Long_t)nEvtsDiffFlowLYZ2RP; | |
2300 | legendDiffFlowEtaRP->AddEntry(lyz2CommonHistRes->GetHistDiffFlowEtaRP(),entryDiffLYZ2EtaRP->Data(),"p"); | |
2301 | } | |
2302 | ||
2303 | //LYZEP | |
2304 | if(plotLYZEPEtaRP && lyzepCommonHistRes) | |
2305 | { | |
2306 | entryDiffLYZEPEtaRP->Append("M = "); | |
2307 | (*entryDiffLYZEPEtaRP)+=(Long_t)avMultDiffFlowLYZEPRP; | |
2308 | entryDiffLYZEPEtaRP->Append(", N = "); | |
2309 | (*entryDiffLYZEPEtaRP)+=(Long_t)nEvtsDiffFlowLYZEPRP; | |
2310 | legendDiffFlowEtaRP->AddEntry(lyzepCommonHistRes->GetHistDiffFlowEtaRP(),entryDiffLYZEPEtaRP->Data(),"p"); | |
2311 | } | |
1fca9c90 | 2312 | |
dec3fab6 | 2313 | //drawing finally the legend in the 2nd pad: |
2314 | if(textDefault) textDefault->Draw(); | |
2315 | ||
e83922f2 | 2316 | if(legendDiffFlowEtaRP) |
2317 | { | |
2318 | legendDiffFlowEtaRP->SetMargin(0.15); | |
2319 | legendDiffFlowEtaRP->Draw(); | |
2320 | } | |
2321 | }// end of if(plotDiffFlowEtaRP) | |
1fca9c90 | 2322 | //---------------------------------------------------------------------------------- |
688877b8 | 2323 | |
2324 | ||
2325 | ||
2326 | ||
2327 | ||
2328 | ||
2329 | ||
2330 | ||
2331 | ||
2332 | ||
2333 | ||
2334 | ||
2335 | ||
2336 | ||
2337 | ||
2338 | ||
2339 | ||
2340 | ||
2341 | ||
2342 | ||
2343 | ||
2344 | ||
2345 | ||
2346 | ||
2347 | ||
2348 | ||
2349 | ||
2350 | ||
2351 | ||
2352 | ||
2353 | ||
2354 | ||
2355 | ||
2356 | ||
2357 | //---------------------------------------------------------------------------------- | |
2358 | // final drawing for plot |(v{method}-v{MC})/v{MC}| as a function of pt for RPs | |
2359 | if(plotDiffFlowPtRelativeToMCRP) | |
2360 | { | |
2361 | TCanvas* diffFlowPtRelativeToMCRP = new TCanvas("Differential Flow (Pt) of RPs relative to MC","Differential Flow (Pt) of RPs relative to MC",1000,600); | |
2362 | ||
2363 | diffFlowPtRelativeToMCRP->Divide(2,1); | |
2364 | ||
2365 | //1st pad is for plot: | |
2366 | (diffFlowPtRelativeToMCRP->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
2367 | ||
2368 | if(styleHistPt) | |
2369 | { | |
2370 | (styleHistPt->GetYaxis())->SetRangeUser(-0.3,1.0); | |
2371 | styleHistPt->Draw(); | |
2372 | } | |
2373 | ||
2374 | TH1D *spDiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2375 | TH1D *gfc2DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2376 | TH1D *gfc4DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2377 | TH1D *gfc6DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2378 | TH1D *gfc8DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2379 | TH1D *qc2DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2380 | TH1D *qc4DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2381 | TH1D *qc6DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2382 | TH1D *qc8DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2383 | TH1D *lyz2DiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2384 | TH1D *lyzepDiffFlowPtRelativeToMCRP = new TH1D("","",iNbinsPt,dPtMin,dPtMax); | |
2385 | ||
2386 | if(mcepCommonHistRes->GetHistDiffFlowPtRP()) | |
2387 | { | |
2388 | for(Int_t p=1;p<iNbinsPt+1;p++) | |
2389 | { | |
2390 | Double_t dvnMC = (mcepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2391 | ||
2392 | // SP: | |
2393 | Double_t dvnSP = 0.; | |
2394 | if(spCommonHistRes && spCommonHistRes->GetHistDiffFlowPtRP()) | |
2395 | { | |
2396 | dvnSP = (spCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2397 | } | |
2398 | if(dvnMC!=0. && dvnSP!=0.) | |
2399 | { | |
2400 | spDiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnSP-dvnMC)/dvnMC)); | |
2401 | } | |
2402 | ||
2403 | // GFC{2}: | |
2404 | Double_t dvnGFC2 = 0.; | |
2405 | if(gfcCommonHistRes2 && gfcCommonHistRes2->GetHistDiffFlowPtRP()) | |
2406 | { | |
2407 | dvnGFC2 = (gfcCommonHistRes2->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2408 | } | |
2409 | if(dvnMC!=0. && dvnGFC2!=0.) | |
2410 | { | |
2411 | gfc2DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnGFC2-dvnMC)/dvnMC)); | |
2412 | } | |
2413 | ||
2414 | // GFC{4}: | |
2415 | Double_t dvnGFC4 = 0.; | |
2416 | if(gfcCommonHistRes4 && gfcCommonHistRes4->GetHistDiffFlowPtRP()) | |
2417 | { | |
2418 | dvnGFC4 = (gfcCommonHistRes4->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2419 | } | |
2420 | if(dvnMC!=0. && dvnGFC4!=0.) | |
2421 | { | |
2422 | gfc4DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnGFC4-dvnMC)/dvnMC)); | |
2423 | } | |
2424 | ||
2425 | // GFC{6}: | |
2426 | Double_t dvnGFC6 = 0.; | |
2427 | if(gfcCommonHistRes6 && gfcCommonHistRes6->GetHistDiffFlowPtRP()) | |
2428 | { | |
2429 | dvnGFC6 = (gfcCommonHistRes6->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2430 | } | |
2431 | if(dvnMC!=0. && dvnGFC6!=0.) | |
2432 | { | |
2433 | gfc6DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnGFC6-dvnMC)/dvnMC)); | |
2434 | } | |
2435 | ||
2436 | // GFC{8}: | |
2437 | Double_t dvnGFC8 = 0.; | |
2438 | if(gfcCommonHistRes8 && gfcCommonHistRes8->GetHistDiffFlowPtRP()) | |
2439 | { | |
2440 | dvnGFC8 = (gfcCommonHistRes8->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2441 | } | |
2442 | if(dvnMC!=0. && dvnGFC8!=0.) | |
2443 | { | |
2444 | gfc8DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnGFC8-dvnMC)/dvnMC)); | |
2445 | } | |
2446 | ||
2447 | // QC{2}: | |
2448 | Double_t dvnQC2 = 0.; | |
2449 | if(qcCommonHistRes2 && qcCommonHistRes2->GetHistDiffFlowPtRP()) | |
2450 | { | |
2451 | dvnQC2 = (qcCommonHistRes2->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2452 | } | |
2453 | if(dvnMC!=0. && dvnQC2!=0.) | |
2454 | { | |
2455 | qc2DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnQC2-dvnMC)/dvnMC)); | |
2456 | } | |
2457 | ||
2458 | // QC{4}: | |
2459 | Double_t dvnQC4 = 0.; | |
2460 | if(qcCommonHistRes4 && qcCommonHistRes4->GetHistDiffFlowPtRP()) | |
2461 | { | |
2462 | dvnQC4 = (qcCommonHistRes4->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2463 | } | |
2464 | if(dvnMC!=0. && dvnQC4!=0.) | |
2465 | { | |
2466 | qc4DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnQC4-dvnMC)/dvnMC)); | |
2467 | } | |
2468 | ||
2469 | // QC{6}: | |
2470 | Double_t dvnQC6 = 0.; | |
2471 | if(qcCommonHistRes6 && qcCommonHistRes6->GetHistDiffFlowPtRP()) | |
2472 | { | |
2473 | dvnQC6 = (qcCommonHistRes6->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2474 | } | |
2475 | if(dvnMC!=0. && dvnQC6!=0.) | |
2476 | { | |
2477 | qc6DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnQC6-dvnMC)/dvnMC)); | |
2478 | } | |
2479 | ||
2480 | // QC{8}: | |
2481 | Double_t dvnQC8 = 0.; | |
2482 | if(qcCommonHistRes8 && qcCommonHistRes8->GetHistDiffFlowPtRP()) | |
2483 | { | |
2484 | dvnQC8 = (qcCommonHistRes8->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2485 | } | |
2486 | if(dvnMC!=0. && dvnQC8!=0.) | |
2487 | { | |
2488 | qc8DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnQC8-dvnMC)/dvnMC)); | |
2489 | } | |
2490 | ||
2491 | // LYZ2: | |
2492 | Double_t dvnLYZ2 = 0.; | |
2493 | if(lyz2CommonHistRes && lyz2CommonHistRes->GetHistDiffFlowPtRP()) | |
2494 | { | |
2495 | dvnLYZ2 = (lyz2CommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2496 | } | |
2497 | if(dvnMC!=0. && dvnLYZ2!=0.) | |
2498 | { | |
2499 | lyz2DiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnLYZ2-dvnMC)/dvnMC)); | |
2500 | } | |
2501 | ||
2502 | // LYZEP: | |
2503 | Double_t dvnLYZEP = 0.; | |
2504 | if(lyzepCommonHistRes && lyzepCommonHistRes->GetHistDiffFlowPtRP()) | |
2505 | { | |
2506 | dvnLYZEP = (lyzepCommonHistRes->GetHistDiffFlowPtRP())->GetBinContent(p); | |
2507 | } | |
2508 | if(dvnMC!=0. && dvnLYZEP!=0.) | |
2509 | { | |
2510 | lyzepDiffFlowPtRelativeToMCRP->SetBinContent(p,TMath::Abs((dvnLYZEP-dvnMC)/dvnMC)); | |
2511 | } | |
2512 | ||
2513 | } // end of for(Int_t p=1;p<iNbinsPt+1;p++) | |
2514 | } // end of if(mcepCommonHistRes->GetHistDiffFlowPtRP()) | |
2515 | ||
2516 | ||
2517 | // final drawings: | |
2518 | spDiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorSP); | |
2519 | spDiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleSP); | |
2520 | if(plotSPRelativeToMCRP && spCommonHistRes) spDiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2521 | ||
2522 | gfc2DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorGFC2); | |
2523 | gfc2DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleGFC2); | |
2524 | if(plotGFC2RelativeToMCRP && gfcCommonHistRes2) gfc2DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2525 | ||
2526 | gfc4DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorGFC4); | |
2527 | gfc4DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleGFC4); | |
2528 | if(plotGFC4RelativeToMCRP && gfcCommonHistRes4) gfc4DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2529 | ||
2530 | gfc6DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorGFC6); | |
2531 | gfc6DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleGFC6); | |
2532 | if(plotGFC6RelativeToMCRP && gfcCommonHistRes6) gfc6DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2533 | ||
2534 | gfc8DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorGFC8); | |
2535 | gfc8DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleGFC8); | |
2536 | if(plotGFC8RelativeToMCRP && gfcCommonHistRes8) gfc8DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2537 | ||
2538 | qc2DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorQC2); | |
2539 | qc2DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleQC2); | |
2540 | if(plotQC2RelativeToMCRP && qcCommonHistRes2) qc2DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2541 | ||
2542 | qc4DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorQC4); | |
2543 | qc4DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleQC4); | |
2544 | if(plotQC4RelativeToMCRP && qcCommonHistRes4) qc4DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2545 | ||
2546 | qc6DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorQC6); | |
2547 | qc6DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleQC6); | |
2548 | if(plotQC6RelativeToMCRP && qcCommonHistRes6) qc6DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2549 | ||
2550 | qc8DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorQC8); | |
2551 | qc8DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleQC8); | |
2552 | if(plotQC8RelativeToMCRP && qcCommonHistRes8) qc8DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2553 | ||
2554 | lyz2DiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorLYZ2); | |
2555 | lyz2DiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleLYZ2); | |
2556 | if(plotLYZ2RelativeToMCRP && lyz2CommonHistRes) lyz2DiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2557 | ||
2558 | lyzepDiffFlowPtRelativeToMCRP->SetMarkerColor(markerColorLYZEP); | |
2559 | lyzepDiffFlowPtRelativeToMCRP->SetMarkerStyle(markerStyleLYZEP); | |
2560 | if(plotLYZEPRelativeToMCRP && lyzepCommonHistRes) lyzepDiffFlowPtRelativeToMCRP->Draw("PSAME"); | |
2561 | ||
2562 | //2nd pad is for legend: | |
2563 | (diffFlowPtRelativeToMCRP->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
2564 | ||
2565 | TLegend* legendDiffFlowPtRP = new TLegend(0.02,0.12,0.97,0.70); | |
2566 | legendDiffFlowPtRP->SetTextFont(72); | |
2567 | legendDiffFlowPtRP->SetTextSize(0.06); | |
2568 | ||
2569 | //legend's entries:Pt | |
2570 | TString *entryDiffMCPtRP = new TString("MC ....... "); | |
2571 | TString *entryDiffSPPtRP = new TString("SP ....... "); | |
2572 | TString *entryDiffGFC2PtRP = new TString("GFC{2} ... "); | |
2573 | TString *entryDiffGFC4PtRP = new TString("GFC{4} ... "); | |
2574 | TString *entryDiffGFC6PtRP = new TString("GFC{6} ... "); | |
2575 | TString *entryDiffGFC8PtRP = new TString("GFC{8} ... "); | |
2576 | TString *entryDiffQC2PtRP = new TString("QC{2} .... "); | |
2577 | TString *entryDiffQC4PtRP = new TString("QC{4} .... "); | |
2578 | TString *entryDiffQC6PtRP = new TString("QC{6} .... "); | |
2579 | TString *entryDiffQC8PtRP = new TString("QC{8} .... "); | |
2580 | TString *entryDiffLYZ2PtRP = new TString("LYZ ...... "); | |
2581 | TString *entryDiffLYZEPPtRP = new TString("LYZEP ... "); | |
2582 | ||
2583 | //MC | |
2584 | if(mcepCommonHistRes) | |
2585 | { | |
2586 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillStyle(meshStyleDiffFlowPtRP); | |
2587 | (mcepCommonHistRes->GetHistDiffFlowPtRP())->SetFillColor(meshColorDiffFlowPtRP); | |
2588 | entryDiffMCPtRP->Append("M = "); | |
2589 | (*entryDiffMCPtRP)+=(Long_t)avMultDiffFlowMCRP; | |
2590 | entryDiffMCPtRP->Append(", N = "); | |
2591 | (*entryDiffMCPtRP)+=(Long_t)nEvtsDiffFlowMCRP; | |
2592 | //legendDiffFlowPtRP->AddEntry(mcepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffMCPtRP->Data(),"f"); | |
2593 | } | |
2594 | ||
2595 | //SP | |
2596 | if(spCommonHistRes) | |
2597 | { | |
2598 | entryDiffSPPtRP->Append("M = "); | |
2599 | (*entryDiffSPPtRP)+=(Long_t)avMultDiffFlowSPRP; | |
2600 | entryDiffSPPtRP->Append(", N = "); | |
2601 | (*entryDiffSPPtRP)+=(Long_t)nEvtsDiffFlowSPRP; | |
2602 | if(plotSPRelativeToMCRP) legendDiffFlowPtRP->AddEntry(spCommonHistRes->GetHistDiffFlowPtRP(),entryDiffSPPtRP->Data(),"p"); | |
2603 | } | |
2604 | ||
2605 | //GFC | |
2606 | if(plotGFC2PtRP && gfcCommonHistRes2) | |
2607 | { | |
2608 | entryDiffGFC2PtRP->Append("M = "); | |
2609 | (*entryDiffGFC2PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2610 | entryDiffGFC2PtRP->Append(", N = "); | |
2611 | (*entryDiffGFC2PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2612 | if(plotGFC2RelativeToMCRP) legendDiffFlowPtRP->AddEntry(gfcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffGFC2PtRP->Data(),"p"); | |
2613 | } | |
2614 | if(plotGFC4PtRP && gfcCommonHistRes4) | |
2615 | { | |
2616 | entryDiffGFC4PtRP->Append("M = "); | |
2617 | (*entryDiffGFC4PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2618 | entryDiffGFC4PtRP->Append(", N = "); | |
2619 | (*entryDiffGFC4PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2620 | if(plotGFC4RelativeToMCRP) legendDiffFlowPtRP->AddEntry(gfcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffGFC4PtRP->Data(),"p"); | |
2621 | } | |
2622 | if(plotGFC6PtRP && gfcCommonHistRes6) | |
2623 | { | |
2624 | entryDiffGFC6PtRP->Append("M = "); | |
2625 | (*entryDiffGFC6PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2626 | entryDiffGFC6PtRP->Append(", N = "); | |
2627 | (*entryDiffGFC6PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2628 | if(plotGFC6RelativeToMCRP) legendDiffFlowPtRP->AddEntry(gfcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffGFC6PtRP->Data(),"p"); | |
2629 | } | |
2630 | if(plotGFC8PtRP && gfcCommonHistRes8) | |
2631 | { | |
2632 | entryDiffGFC8PtRP->Append("M = "); | |
2633 | (*entryDiffGFC8PtRP)+=(Long_t)avMultDiffFlowGFCRP; | |
2634 | entryDiffGFC8PtRP->Append(", N = "); | |
2635 | (*entryDiffGFC8PtRP)+=(Long_t)nEvtsDiffFlowGFCRP; | |
2636 | if(plotGFC8RelativeToMCRP) legendDiffFlowPtRP->AddEntry(gfcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffGFC8PtRP->Data(),"p"); | |
2637 | } | |
2638 | ||
2639 | //QC | |
2640 | if(plotQC2PtRP && qcCommonHistRes2) | |
2641 | { | |
2642 | entryDiffQC2PtRP->Append("M = "); | |
2643 | (*entryDiffQC2PtRP)+=(Long_t)avMultDiffFlowQC2RP; | |
2644 | entryDiffQC2PtRP->Append(", N = "); | |
2645 | (*entryDiffQC2PtRP)+=(Long_t)nEvtsDiffFlowQC2RP; | |
2646 | if(plotQC2RelativeToMCRP) legendDiffFlowPtRP->AddEntry(qcCommonHistRes2->GetHistDiffFlowPtRP(),entryDiffQC2PtRP->Data(),"p"); | |
2647 | } | |
2648 | if(plotQC4PtRP && qcCommonHistRes4) | |
2649 | { | |
2650 | entryDiffQC4PtRP->Append("M = "); | |
2651 | (*entryDiffQC4PtRP)+=(Long_t)avMultDiffFlowQC4RP; | |
2652 | entryDiffQC4PtRP->Append(", N = "); | |
2653 | (*entryDiffQC4PtRP)+=(Long_t)nEvtsDiffFlowQC4RP; | |
2654 | if(plotQC4RelativeToMCRP) legendDiffFlowPtRP->AddEntry(qcCommonHistRes4->GetHistDiffFlowPtRP(),entryDiffQC4PtRP->Data(),"p"); | |
2655 | } | |
2656 | if(plotQC6PtRP && qcCommonHistRes6) | |
2657 | { | |
2658 | entryDiffQC6PtRP->Append("M = "); | |
2659 | (*entryDiffQC6PtRP)+=(Long_t)avMultDiffFlowQC6RP; | |
2660 | entryDiffQC6PtRP->Append(", N = "); | |
2661 | (*entryDiffQC6PtRP)+=(Long_t)nEvtsDiffFlowQC6RP; | |
2662 | if(plotQC6RelativeToMCRP) legendDiffFlowPtRP->AddEntry(qcCommonHistRes6->GetHistDiffFlowPtRP(),entryDiffQC6PtRP->Data(),"p"); | |
2663 | } | |
2664 | if(plotQC8PtRP && qcCommonHistRes8) | |
2665 | { | |
2666 | entryDiffQC8PtRP->Append("M = "); | |
2667 | (*entryDiffQC8PtRP)+=(Long_t)avMultDiffFlowQC8RP; | |
2668 | entryDiffQC8PtRP->Append(", N = "); | |
2669 | (*entryDiffQC8PtRP)+=(Long_t)nEvtsDiffFlowQC8RP; | |
2670 | if(plotQC8RelativeToMCRP) legendDiffFlowPtRP->AddEntry(qcCommonHistRes8->GetHistDiffFlowPtRP(),entryDiffQC8PtRP->Data(),"p"); | |
2671 | } | |
2672 | ||
2673 | //LYZ2 | |
2674 | if(plotLYZ2PtRP && lyz2CommonHistRes) | |
2675 | { | |
2676 | entryDiffLYZ2PtRP->Append("M = "); | |
2677 | (*entryDiffLYZ2PtRP)+=(Long_t)avMultDiffFlowLYZ2RP; | |
2678 | entryDiffLYZ2PtRP->Append(", N = "); | |
2679 | (*entryDiffLYZ2PtRP)+=(Long_t)nEvtsDiffFlowLYZ2RP; | |
2680 | if(plotLYZ2RelativeToMCRP) legendDiffFlowPtRP->AddEntry(lyz2CommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZ2PtRP->Data(),"p"); | |
2681 | } | |
2682 | ||
2683 | //LYZEP | |
2684 | if(plotLYZEPPtRP && lyzepCommonHistRes) | |
2685 | { | |
2686 | entryDiffLYZEPPtRP->Append("M = "); | |
2687 | (*entryDiffLYZEPPtRP)+=(Long_t)avMultDiffFlowLYZEPRP; | |
2688 | entryDiffLYZEPPtRP->Append(", N = "); | |
2689 | (*entryDiffLYZEPPtRP)+=(Long_t)nEvtsDiffFlowLYZEPRP; | |
2690 | if(plotLYZEPRelativeToMCRP) legendDiffFlowPtRP->AddEntry(lyzepCommonHistRes->GetHistDiffFlowPtRP(),entryDiffLYZEPPtRP->Data(),"p"); | |
2691 | } | |
2692 | ||
2693 | //drawing finally the legend in the 2nd pad: | |
2694 | if(textDefault) textDefault->Draw(); | |
2695 | ||
2696 | if(legendDiffFlowPtRP) | |
2697 | { | |
2698 | legendDiffFlowPtRP->SetMargin(0.15); | |
2699 | legendDiffFlowPtRP->Draw(); | |
2700 | } | |
2701 | } | |
2702 | //---------------------------------------------------------------------------------- | |
1fca9c90 | 2703 | |
2704 | //---------------------------------------------------------------------------------- | |
2705 | //final drawing for differential flow (Pt, POI): | |
e83922f2 | 2706 | if(plotDiffFlowPtPOI) |
2707 | { | |
2708 | TCanvas* diffFlowPtAllCanvasPOI = new TCanvas("Differential Flow (Pt) of POI","Differential Flow (Pt) of POI ",1000,600); | |
2709 | ||
2710 | diffFlowPtAllCanvasPOI->Divide(2,1); | |
2711 | ||
2712 | //1st pad is for plot: | |
2713 | (diffFlowPtAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
2714 | ||
2715 | if(styleHistPt) | |
2716 | { | |
2717 | (styleHistPt->GetYaxis())->SetRangeUser(-0.3,1.0); | |
2718 | styleHistPt->Draw(); | |
2719 | } | |
2720 | if(pMeshDiffFlowPtPOI) | |
2721 | { | |
2722 | pMeshDiffFlowPtPOI->Draw("LFSAME"); | |
2723 | } | |
1fca9c90 | 2724 | |
e83922f2 | 2725 | //MC |
2726 | if(plotMCPtPOI && mcepCommonHistRes) | |
905bed86 | 2727 | { |
e83922f2 | 2728 | (mcepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); |
2729 | } | |
905bed86 | 2730 | //SP |
2731 | if(plotSPPtPOI && spCommonHistRes) | |
2732 | { | |
2733 | (spCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2734 | } | |
e83922f2 | 2735 | //GFC |
2736 | if(plotGFC2PtPOI && gfcCommonHistRes2) | |
2737 | { | |
2738 | (gfcCommonHistRes2->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2739 | } | |
2740 | if(plotGFC4PtPOI && gfcCommonHistRes4) | |
2741 | { | |
2742 | (gfcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2743 | } | |
2744 | if(plotGFC6PtPOI && gfcCommonHistRes6) | |
2745 | { | |
2746 | (gfcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2747 | } | |
2748 | if(plotGFC8PtPOI && gfcCommonHistRes8) | |
2749 | { | |
2750 | (gfcCommonHistRes8->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2751 | } | |
2752 | //QC | |
2753 | if(plotQC2PtPOI && qcCommonHistRes2) | |
2754 | { | |
2755 | (qcCommonHistRes2->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2756 | } | |
2757 | if(plotQC4PtPOI && qcCommonHistRes4) | |
2758 | { | |
2759 | (qcCommonHistRes4->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2760 | } | |
2761 | if(plotQC6PtPOI && qcCommonHistRes6) | |
3963f25c | 2762 | { kTRUE; |
7d02d77f | 2763 | //(qcCommonHistRes6->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); |
e83922f2 | 2764 | } |
2765 | if(plotQC8PtPOI && qcCommonHistRes8) | |
2766 | { | |
7d02d77f | 2767 | //(qcCommonHistRes8->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); |
e83922f2 | 2768 | } |
2769 | //LYZ2 | |
2770 | if(plotLYZ2PtPOI && lyz2CommonHistRes) | |
2771 | { | |
2772 | (lyz2CommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2773 | } | |
2774 | //LYZEP | |
2775 | if(plotLYZEPPtPOI && lyzepCommonHistRes) | |
2776 | { | |
2777 | (lyzepCommonHistRes->GetHistDiffFlowPtPOI())->Draw("E1PSAME"); | |
2778 | } | |
1fca9c90 | 2779 | |
e83922f2 | 2780 | //2nd pad is for legend: |
2781 | (diffFlowPtAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
2782 | ||
dec3fab6 | 2783 | TLegend* legendDiffFlowPtPOI = new TLegend(0.02,0.12,0.97,0.70); |
e83922f2 | 2784 | legendDiffFlowPtPOI->SetTextFont(72); |
2785 | legendDiffFlowPtPOI->SetTextSize(0.06); | |
2786 | ||
2787 | //legend's entries: | |
2788 | TString *entryDiffMCPtPOI = new TString("MC ....... "); | |
905bed86 | 2789 | TString *entryDiffSPPtPOI = new TString("SP ....... "); |
e83922f2 | 2790 | TString *entryDiffGFC2PtPOI = new TString("GFC{2} ... "); |
2791 | TString *entryDiffGFC4PtPOI = new TString("GFC{4} ... "); | |
2792 | TString *entryDiffGFC6PtPOI = new TString("GFC{6} ... "); | |
2793 | TString *entryDiffGFC8PtPOI = new TString("GFC{8} ... "); | |
2794 | TString *entryDiffQC2PtPOI = new TString("QC{2} .... "); | |
2795 | TString *entryDiffQC4PtPOI = new TString("QC{4} .... "); | |
2796 | TString *entryDiffQC6PtPOI = new TString("QC{6} .... "); | |
2797 | TString *entryDiffQC8PtPOI = new TString("QC{8} .... "); | |
2798 | TString *entryDiffLYZ2PtPOI = new TString("LYZ ...... "); | |
2799 | TString *entryDiffLYZEPPtPOI = new TString("LYZEP ... "); | |
2800 | ||
2801 | //MC | |
2802 | if(mcepCommonHistRes) | |
2803 | { | |
dec3fab6 | 2804 | (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetFillStyle(meshStyleDiffFlowPtPOI); |
2805 | (mcepCommonHistRes->GetHistDiffFlowPtPOI())->SetFillColor(meshColorDiffFlowPtPOI); | |
e83922f2 | 2806 | entryDiffMCPtPOI->Append("M = "); |
2807 | (*entryDiffMCPtPOI)+=(Long_t)avMultDiffFlowMCPOI; | |
2808 | entryDiffMCPtPOI->Append(", N = "); | |
2809 | (*entryDiffMCPtPOI)+=(Long_t)nEvtsDiffFlowMCPOI; | |
2810 | legendDiffFlowPtPOI->AddEntry(mcepCommonHistRes->GetHistDiffFlowPtPOI(),entryDiffMCPtPOI->Data(),"f"); | |
2811 | } | |
905bed86 | 2812 | |
2813 | //SP | |
2814 | if(spCommonHistRes) | |
2815 | { | |
905bed86 | 2816 | entryDiffSPPtPOI->Append("M = "); |
2817 | (*entryDiffSPPtPOI)+=(Long_t)avMultDiffFlowSPPOI; | |
2818 | entryDiffSPPtPOI->Append(", N = "); | |
2819 | (*entryDiffSPPtPOI)+=(Long_t)nEvtsDiffFlowSPPOI; | |
2820 | legendDiffFlowPtPOI->AddEntry(spCommonHistRes->GetHistDiffFlowPtPOI(),entryDiffSPPtPOI->Data(),"p"); | |
2821 | } | |
1fca9c90 | 2822 | |
e83922f2 | 2823 | //GFC |
2824 | if(plotGFC2PtPOI && gfcCommonHistRes2) | |
2825 | { | |
2826 | entryDiffGFC2PtPOI->Append("M = "); | |
2827 | (*entryDiffGFC2PtPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
2828 | entryDiffGFC2PtPOI->Append(", N = "); | |
2829 | (*entryDiffGFC2PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
2830 | legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes2->GetHistDiffFlowPtPOI(),entryDiffGFC2PtPOI->Data(),"p"); | |
2831 | } | |
2832 | if(plotGFC4PtPOI && gfcCommonHistRes4) | |
2833 | { | |
2834 | entryDiffGFC4PtPOI->Append("M = "); | |
2835 | (*entryDiffGFC4PtPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
2836 | entryDiffGFC4PtPOI->Append(", N = "); | |
2837 | (*entryDiffGFC4PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
2838 | legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes4->GetHistDiffFlowPtPOI(),entryDiffGFC4PtPOI->Data(),"p"); | |
2839 | } | |
2840 | if(plotGFC6PtPOI && gfcCommonHistRes6) | |
2841 | { | |
2842 | entryDiffGFC6PtPOI->Append("M = "); | |
2843 | (*entryDiffGFC6PtPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
2844 | entryDiffGFC6PtPOI->Append(", N = "); | |
2845 | (*entryDiffGFC6PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
2846 | legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes6->GetHistDiffFlowPtPOI(),entryDiffGFC6PtPOI->Data(),"p"); | |
2847 | } | |
2848 | if(plotGFC8PtPOI && gfcCommonHistRes8) | |
2849 | { | |
2850 | entryDiffGFC8PtPOI->Append("M = "); | |
2851 | (*entryDiffGFC8PtPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
2852 | entryDiffGFC8PtPOI->Append(", N = "); | |
2853 | (*entryDiffGFC8PtPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
2854 | legendDiffFlowPtPOI->AddEntry(gfcCommonHistRes8->GetHistDiffFlowPtPOI(),entryDiffGFC8PtPOI->Data(),"p"); | |
2855 | } | |
2856 | ||
2857 | //QC | |
2858 | if(plotQC2PtPOI && qcCommonHistRes2) | |
2859 | { | |
2860 | entryDiffQC2PtPOI->Append("M = "); | |
2861 | (*entryDiffQC2PtPOI)+=(Long_t)avMultDiffFlowQC2POI; | |
2862 | entryDiffQC2PtPOI->Append(", N = "); | |
2863 | (*entryDiffQC2PtPOI)+=(Long_t)nEvtsDiffFlowQC2POI; | |
2864 | legendDiffFlowPtPOI->AddEntry(qcCommonHistRes2->GetHistDiffFlowPtPOI(),entryDiffQC2PtPOI->Data(),"p"); | |
2865 | } | |
2866 | if(plotQC4PtPOI && qcCommonHistRes4) | |
2867 | { | |
2868 | entryDiffQC4PtPOI->Append("M = "); | |
2869 | (*entryDiffQC4PtPOI)+=(Long_t)avMultDiffFlowQC4POI; | |
2870 | entryDiffQC4PtPOI->Append(", N = "); | |
2871 | (*entryDiffQC4PtPOI)+=(Long_t)nEvtsDiffFlowQC4POI; | |
2872 | legendDiffFlowPtPOI->AddEntry(qcCommonHistRes4->GetHistDiffFlowPtPOI(),entryDiffQC4PtPOI->Data(),"p"); | |
2873 | } | |
2874 | if(plotQC6PtPOI && qcCommonHistRes6) | |
2875 | { | |
2876 | entryDiffQC6PtPOI->Append("M = "); | |
2877 | (*entryDiffQC6PtPOI)+=(Long_t)avMultDiffFlowQC6POI; | |
2878 | entryDiffQC6PtPOI->Append(", N = "); | |
2879 | (*entryDiffQC6PtPOI)+=(Long_t)nEvtsDiffFlowQC6POI; | |
2880 | legendDiffFlowPtPOI->AddEntry(qcCommonHistRes6->GetHistDiffFlowPtPOI(),entryDiffQC6PtPOI->Data(),"p"); | |
2881 | } | |
2882 | if(plotQC8PtPOI && qcCommonHistRes8) | |
2883 | { | |
2884 | entryDiffQC8PtPOI->Append("M = "); | |
2885 | (*entryDiffQC8PtPOI)+=(Long_t)avMultDiffFlowQC8POI; | |
2886 | entryDiffQC8PtPOI->Append(", N = "); | |
2887 | (*entryDiffQC8PtPOI)+=(Long_t)nEvtsDiffFlowQC8POI; | |
2888 | legendDiffFlowPtPOI->AddEntry(qcCommonHistRes8->GetHistDiffFlowPtPOI(),entryDiffQC8PtPOI->Data(),"p"); | |
2889 | } | |
1fca9c90 | 2890 | |
e83922f2 | 2891 | //LYZ2 |
2892 | if(plotLYZ2PtPOI && lyz2CommonHistRes) | |
2893 | { | |
2894 | entryDiffLYZ2PtPOI->Append("M = "); | |
2895 | (*entryDiffLYZ2PtPOI)+=(Long_t)avMultDiffFlowLYZ2POI; | |
2896 | entryDiffLYZ2PtPOI->Append(", N = "); | |
2897 | (*entryDiffLYZ2PtPOI)+=(Long_t)nEvtsDiffFlowLYZ2POI; | |
2898 | legendDiffFlowPtPOI->AddEntry(lyz2CommonHistRes->GetHistDiffFlowPtPOI(),entryDiffLYZ2PtPOI->Data(),"p"); | |
2899 | } | |
2900 | ||
2901 | //LYZEP | |
2902 | if(plotLYZEPPtPOI && lyzepCommonHistRes) | |
2903 | { | |
2904 | entryDiffLYZEPPtPOI->Append("M = "); | |
2905 | (*entryDiffLYZEPPtPOI)+=(Long_t)avMultDiffFlowLYZEPPOI; | |
2906 | entryDiffLYZEPPtPOI->Append(", N = "); | |
2907 | (*entryDiffLYZEPPtPOI)+=(Long_t)nEvtsDiffFlowLYZEPPOI; | |
2908 | legendDiffFlowPtPOI->AddEntry(lyzepCommonHistRes->GetHistDiffFlowPtPOI(),entryDiffLYZEPPtPOI->Data(),"p"); | |
2909 | } | |
1fca9c90 | 2910 | |
dec3fab6 | 2911 | //drawing finally the legend in the 2nd pad: |
2912 | if(textDefault) textDefault->Draw(); | |
2913 | ||
e83922f2 | 2914 | if(legendDiffFlowPtPOI) |
2915 | { | |
2916 | legendDiffFlowPtPOI->SetMargin(0.15); | |
2917 | legendDiffFlowPtPOI->Draw(); | |
2918 | } | |
2919 | }//end of if(plotDiffFlowPtPOI) | |
1fca9c90 | 2920 | //---------------------------------------------------------------------------------- |
688877b8 | 2921 | |
1fca9c90 | 2922 | |
2923 | //---------------------------------------------------------------------------------- | |
2924 | //final drawing for differential flow (Eta, POI): | |
e83922f2 | 2925 | if(plotDiffFlowEtaPOI) |
2926 | { | |
2927 | TCanvas* diffFlowEtaAllCanvasPOI = new TCanvas("Differential Flow (Eta) of POI","Differential Flow (Eta) of POI ",1000,600); | |
2928 | ||
2929 | diffFlowEtaAllCanvasPOI->Divide(2,1); | |
2930 | ||
2931 | //1st pad is for plot: | |
2932 | (diffFlowEtaAllCanvasPOI->cd(1))->SetPad(0.0,0.0,0.75,1.0); | |
1fca9c90 | 2933 | |
e83922f2 | 2934 | if(styleHistEta) |
2935 | { | |
2936 | (styleHistEta->GetYaxis())->SetRangeUser(-0.3,1.0); | |
2937 | styleHistEta->Draw(); | |
2938 | } | |
2939 | if(pMeshDiffFlowEtaPOI) | |
2940 | { | |
2941 | pMeshDiffFlowEtaPOI->Draw("LFSAME"); | |
2942 | } | |
1fca9c90 | 2943 | |
e83922f2 | 2944 | //MC |
2945 | if(plotMCEtaPOI && mcepCommonHistRes) | |
2946 | { | |
2947 | (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2948 | } | |
905bed86 | 2949 | //SP |
2950 | if(plotSPEtaPOI && spCommonHistRes) | |
2951 | { | |
2952 | (spCommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2953 | } | |
e83922f2 | 2954 | //GFC |
2955 | if(plotGFC2EtaPOI && gfcCommonHistRes2) | |
2956 | { | |
2957 | (gfcCommonHistRes2->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2958 | } | |
2959 | if(plotGFC4EtaPOI && gfcCommonHistRes4) | |
2960 | { | |
2961 | (gfcCommonHistRes4->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2962 | } | |
2963 | if(plotGFC6EtaPOI && gfcCommonHistRes6) | |
2964 | { | |
2965 | (gfcCommonHistRes6->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2966 | } | |
2967 | if(plotGFC8EtaPOI && gfcCommonHistRes8) | |
2968 | { | |
2969 | (gfcCommonHistRes8->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2970 | } | |
2971 | //QC | |
2972 | if(plotQC2EtaPOI && qcCommonHistRes2) | |
2973 | { | |
2974 | (qcCommonHistRes2->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2975 | } | |
2976 | if(plotQC4EtaPOI && qcCommonHistRes4) | |
2977 | { | |
2978 | (qcCommonHistRes4->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2979 | } | |
2980 | if(plotQC6EtaPOI && qcCommonHistRes6) | |
2981 | { | |
7d02d77f | 2982 | //(qcCommonHistRes6->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); |
e83922f2 | 2983 | } |
2984 | if(plotQC8EtaPOI && qcCommonHistRes8) | |
2985 | { | |
7d02d77f | 2986 | //(qcCommonHistRes8->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); |
e83922f2 | 2987 | } |
2988 | //LYZ2 | |
2989 | if(plotLYZ2EtaPOI && lyz2CommonHistRes) | |
2990 | { | |
2991 | (lyz2CommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2992 | } | |
2993 | //LYZEP | |
2994 | if(plotLYZEPEtaPOI && lyzepCommonHistRes) | |
2995 | { | |
2996 | (lyzepCommonHistRes->GetHistDiffFlowEtaPOI())->Draw("E1PSAME"); | |
2997 | } | |
1fca9c90 | 2998 | |
e83922f2 | 2999 | //2nd pad is for legend: |
3000 | (diffFlowEtaAllCanvasPOI->cd(2))->SetPad(0.75,0.0,1.0,1.0); | |
3001 | ||
dec3fab6 | 3002 | TLegend* legendDiffFlowEtaPOI = new TLegend(0.02,0.12,0.97,0.70); |
e83922f2 | 3003 | legendDiffFlowEtaPOI->SetTextFont(72); |
3004 | legendDiffFlowEtaPOI->SetTextSize(0.06); | |
3005 | ||
3006 | //legend's entries: | |
3007 | TString *entryDiffMCEtaPOI = new TString("MC ....... "); | |
905bed86 | 3008 | TString *entryDiffSPEtaPOI = new TString("SP ....... "); |
e83922f2 | 3009 | TString *entryDiffGFC2EtaPOI = new TString("GFC{2} ... "); |
3010 | TString *entryDiffGFC4EtaPOI = new TString("GFC{4} ... "); | |
3011 | TString *entryDiffGFC6EtaPOI = new TString("GFC{6} ... "); | |
3012 | TString *entryDiffGFC8EtaPOI = new TString("GFC{8} ... "); | |
3013 | TString *entryDiffQC2EtaPOI = new TString("QC{2} .... "); | |
3014 | TString *entryDiffQC4EtaPOI = new TString("QC{4} .... "); | |
3015 | TString *entryDiffQC6EtaPOI = new TString("QC{6} .... "); | |
3016 | TString *entryDiffQC8EtaPOI = new TString("QC{8} .... "); | |
3017 | TString *entryDiffLYZ2EtaPOI = new TString("LYZ ...... "); | |
3018 | TString *entryDiffLYZEPEtaPOI = new TString("LYZEP ... "); | |
3019 | ||
3020 | //MC | |
3021 | if(mcepCommonHistRes) | |
3022 | { | |
dec3fab6 | 3023 | (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetFillStyle(meshStyleDiffFlowEtaPOI); |
3024 | (mcepCommonHistRes->GetHistDiffFlowEtaPOI())->SetFillColor(meshColorDiffFlowEtaPOI); | |
e83922f2 | 3025 | entryDiffMCEtaPOI->Append("M = "); |
3026 | (*entryDiffMCEtaPOI)+=(Long_t)avMultDiffFlowMCPOI; | |
3027 | entryDiffMCEtaPOI->Append(", N = "); | |
3028 | (*entryDiffMCEtaPOI)+=(Long_t)nEvtsDiffFlowMCPOI; | |
3029 | legendDiffFlowEtaPOI->AddEntry(mcepCommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffMCEtaPOI->Data(),"f"); | |
3030 | } | |
905bed86 | 3031 | |
3032 | //SP | |
3033 | if(spCommonHistRes) | |
3034 | { | |
905bed86 | 3035 | entryDiffSPEtaPOI->Append("M = "); |
3036 | (*entryDiffSPEtaPOI)+=(Long_t)avMultDiffFlowSPPOI; | |
3037 | entryDiffSPEtaPOI->Append(", N = "); | |
3038 | (*entryDiffSPEtaPOI)+=(Long_t)nEvtsDiffFlowSPPOI; | |
3039 | legendDiffFlowEtaPOI->AddEntry(spCommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffSPEtaPOI->Data(),"p"); | |
3040 | } | |
1fca9c90 | 3041 | |
e83922f2 | 3042 | //GFC |
3043 | if(plotGFC2EtaPOI && gfcCommonHistRes2) | |
3044 | { | |
3045 | entryDiffGFC2EtaPOI->Append("M = "); | |
3046 | (*entryDiffGFC2EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
3047 | entryDiffGFC2EtaPOI->Append(", N = "); | |
3048 | (*entryDiffGFC2EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
3049 | legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes2->GetHistDiffFlowEtaPOI(),entryDiffGFC2EtaPOI->Data(),"p"); | |
3050 | } | |
3051 | if(plotGFC4EtaPOI && gfcCommonHistRes4) | |
3052 | { | |
3053 | entryDiffGFC4EtaPOI->Append("M = "); | |
3054 | (*entryDiffGFC4EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
3055 | entryDiffGFC4EtaPOI->Append(", N = "); | |
3056 | (*entryDiffGFC4EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
3057 | legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes4->GetHistDiffFlowEtaPOI(),entryDiffGFC4EtaPOI->Data(),"p"); | |
3058 | } | |
3059 | if(plotGFC6EtaPOI && gfcCommonHistRes6) | |
3060 | { | |
3061 | entryDiffGFC6EtaPOI->Append("M = "); | |
3062 | (*entryDiffGFC6EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
3063 | entryDiffGFC6EtaPOI->Append(", N = "); | |
3064 | (*entryDiffGFC6EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
3065 | legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes6->GetHistDiffFlowEtaPOI(),entryDiffGFC6EtaPOI->Data(),"p"); | |
3066 | } | |
3067 | if(plotGFC8EtaPOI && gfcCommonHistRes8) | |
3068 | { | |
3069 | entryDiffGFC8EtaPOI->Append("M = "); | |
3070 | (*entryDiffGFC8EtaPOI)+=(Long_t)avMultDiffFlowGFCPOI; | |
3071 | entryDiffGFC8EtaPOI->Append(", N = "); | |
3072 | (*entryDiffGFC8EtaPOI)+=(Long_t)nEvtsDiffFlowGFCPOI; | |
3073 | legendDiffFlowEtaPOI->AddEntry(gfcCommonHistRes8->GetHistDiffFlowEtaPOI(),entryDiffGFC8EtaPOI->Data(),"p"); | |
3074 | } | |
3075 | ||
3076 | //QC | |
3077 | if(plotQC2EtaPOI && qcCommonHistRes2) | |
3078 | { | |
3079 | entryDiffQC2EtaPOI->Append("M = "); | |
3080 | (*entryDiffQC2EtaPOI)+=(Long_t)avMultDiffFlowQC2POI; | |
3081 | entryDiffQC2EtaPOI->Append(", N = "); | |
3082 | (*entryDiffQC2EtaPOI)+=(Long_t)nEvtsDiffFlowQC2POI; | |
3083 | legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes2->GetHistDiffFlowEtaPOI(),entryDiffQC2EtaPOI->Data(),"p"); | |
3084 | } | |
3085 | if(plotQC4EtaPOI && qcCommonHistRes4) | |
3086 | { | |
3087 | entryDiffQC4EtaPOI->Append("M = "); | |
3088 | (*entryDiffQC4EtaPOI)+=(Long_t)avMultDiffFlowQC4POI; | |
3089 | entryDiffQC4EtaPOI->Append(", N = "); | |
3090 | (*entryDiffQC4EtaPOI)+=(Long_t)nEvtsDiffFlowQC4POI; | |
3091 | legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes4->GetHistDiffFlowEtaPOI(),entryDiffQC4EtaPOI->Data(),"p"); | |
3092 | } | |
3093 | if(plotQC6EtaPOI && qcCommonHistRes6) | |
3094 | { | |
3095 | entryDiffQC6EtaPOI->Append("M = "); | |
3096 | (*entryDiffQC6EtaPOI)+=(Long_t)avMultDiffFlowQC6POI; | |
3097 | entryDiffQC6EtaPOI->Append(", N = "); | |
3098 | (*entryDiffQC6EtaPOI)+=(Long_t)nEvtsDiffFlowQC6POI; | |
3099 | legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes6->GetHistDiffFlowEtaPOI(),entryDiffQC6EtaPOI->Data(),"p"); | |
3100 | } | |
3101 | if(plotQC8EtaPOI && qcCommonHistRes8) | |
3102 | { | |
3103 | entryDiffQC8EtaPOI->Append("M = "); | |
3104 | (*entryDiffQC8EtaPOI)+=(Long_t)avMultDiffFlowQC8POI; | |
3105 | entryDiffQC8EtaPOI->Append(", N = "); | |
3106 | (*entryDiffQC8EtaPOI)+=(Long_t)nEvtsDiffFlowQC8POI; | |
3107 | legendDiffFlowEtaPOI->AddEntry(qcCommonHistRes8->GetHistDiffFlowEtaPOI(),entryDiffQC8EtaPOI->Data(),"p"); | |
3108 | } | |
1fca9c90 | 3109 | |
e83922f2 | 3110 | //LYZ2 |
3111 | if(plotLYZ2EtaPOI && lyz2CommonHistRes) | |
3112 | { | |
3113 | entryDiffLYZ2EtaPOI->Append("M = "); | |
3114 | (*entryDiffLYZ2EtaPOI)+=(Long_t)avMultDiffFlowLYZ2POI; | |
3115 | entryDiffLYZ2EtaPOI->Append(", N = "); | |
3116 | (*entryDiffLYZ2EtaPOI)+=(Long_t)nEvtsDiffFlowLYZ2POI; | |
3117 | legendDiffFlowEtaPOI->AddEntry(lyz2CommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffLYZ2EtaPOI->Data(),"p"); | |
3118 | } | |
3119 | ||
3120 | //LYZEP | |
3121 | if(plotLYZEPEtaPOI && lyzepCommonHistRes) | |
3122 | { | |
3123 | entryDiffLYZEPEtaPOI->Append("M = "); | |
3124 | (*entryDiffLYZEPEtaPOI)+=(Long_t)avMultDiffFlowLYZEPPOI; | |
3125 | entryDiffLYZEPEtaPOI->Append(", N = "); | |
3126 | (*entryDiffLYZEPEtaPOI)+=(Long_t)nEvtsDiffFlowLYZEPPOI; | |
3127 | legendDiffFlowEtaPOI->AddEntry(lyzepCommonHistRes->GetHistDiffFlowEtaPOI(),entryDiffLYZEPEtaPOI->Data(),"p"); | |
3128 | } | |
1fca9c90 | 3129 | |
dec3fab6 | 3130 | //drawing finally the legend in the 2nd pad: |
3131 | if(textDefault) textDefault->Draw(); | |
3132 | ||
e83922f2 | 3133 | if(legendDiffFlowEtaPOI) |
3134 | { | |
3135 | legendDiffFlowEtaPOI->SetMargin(0.15); | |
3136 | legendDiffFlowEtaPOI->Draw(); | |
3137 | } | |
3138 | }//end of if(plotDiffFlowEtaPOI) | |
1fca9c90 | 3139 | //---------------------------------------------------------------------------------- |
3140 | ||
3141 | ||
a88e5ac2 | 3142 | //===================================================================================== |
3143 | ||
70288839 | 3144 | } |
b25cc698 | 3145 | |
8d5cd720 | 3146 | void LoadPlotLibraries(const libModes mode) { |
b25cc698 | 3147 | |
3148 | //-------------------------------------- | |
3149 | // Load the needed libraries most of them already loaded by aliroot | |
3150 | //-------------------------------------- | |
3151 | gSystem->Load("libTree.so"); | |
3152 | gSystem->Load("libGeom.so"); | |
3153 | gSystem->Load("libVMC.so"); | |
3154 | gSystem->Load("libXMLIO.so"); | |
3155 | gSystem->Load("libPhysics.so"); | |
3156 | ||
3157 | //---------------------------------------------------------- | |
3158 | // >>>>>>>>>>> Local mode <<<<<<<<<<<<<< | |
3159 | //---------------------------------------------------------- | |
3160 | if (mode==mLocal) { | |
3161 | //-------------------------------------------------------- | |
3162 | // If you want to use already compiled libraries | |
3163 | // in the aliroot distribution | |
3164 | //-------------------------------------------------------- | |
3165 | ||
3166 | //================================================================================== | |
3167 | //load needed libraries: | |
3168 | gSystem->AddIncludePath("-I$ROOTSYS/include"); | |
3169 | gSystem->Load("libTree.so"); | |
3170 | ||
3171 | // for AliRoot | |
3172 | gSystem->AddIncludePath("-I$ALICE_ROOT/include"); | |
3173 | gSystem->Load("libANALYSIS.so"); | |
3174 | gSystem->Load("libPWG2flowCommon.so"); | |
3175 | cerr<<"libPWG2flowCommon.so loaded ..."<<endl; | |
3176 | ||
3177 | } | |
3178 | ||
3179 | else if (mode==mLocalSource) { | |
3180 | ||
3181 | // In root inline compile | |
d90df6c3 | 3182 | |
3183 | // Constants | |
3184 | gROOT->LoadMacro("AliFlowCommon/AliFlowCommonConstants.cxx+"); | |
3185 | gROOT->LoadMacro("AliFlowCommon/AliFlowLYZConstants.cxx+"); | |
3186 | gROOT->LoadMacro("AliFlowCommon/AliFlowCumuConstants.cxx+"); | |
3187 | ||
3188 | // Flow event | |
3189 | gROOT->LoadMacro("AliFlowCommon/AliFlowVector.cxx+"); | |
3190 | gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimple.cxx+"); | |
3191 | gROOT->LoadMacro("AliFlowCommon/AliFlowEventSimple.cxx+"); | |
3192 | ||
3193 | // Cuts | |
3194 | gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimpleCuts.cxx+"); | |
b25cc698 | 3195 | |
3196 | // Output histosgrams | |
3197 | gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHist.cxx+"); | |
3198 | gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHistResults.cxx+"); | |
3199 | gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist1.cxx+"); | |
3200 | gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist2.cxx+"); | |
3201 | ||
3202 | cout << "finished loading macros!" << endl; | |
3203 | ||
3204 | } | |
3205 | ||
3206 | } | |
3207 | ||
3208 |