]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGCF/FEMTOSCOPY/macros/Plot_plots.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
(For the train) small change in the macro with parameters.
[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Plot_plots.C
CommitLineData
9b479142 1#include <math.h>
2#include <time.h>
3#include <stdio.h>
4#include <stdlib.h>
5#include <Riostream.h>
6
7#include "TVector2.h"
8#include "TFile.h"
9#include "TString.h"
10#include "TF1.h"
11#include "TH1.h"
12#include "TH2.h"
13#include "TH3.h"
14#include "TProfile.h"
15#include "TProfile2D.h"
16#include "TMath.h"
17#include "TText.h"
18#include "TRandom3.h"
19#include "TArray.h"
20#include "TLegend.h"
21#include "TStyle.h"
22#include "TMinuit.h"
23
24#define BohrR 1963.6885
25#define FmToGeV 0.19733 // conversion to fm
26#define PI 3.1415926
27#define masspiC 0.1395702 // pi+ mass (GeV/c^2)
28
29using namespace std;
30
9b34d0fd 31bool TherminatorC2=kFALSE;
32const int BOI_1=0;// centrality bin (0-9)
9b479142 33const int BOI_2=9;// centrality of second bin for C2 fit parameter plot only
9b479142 34const int ChProdBOI=0;// 0=SameCharge, 1=MixedCharge
35const int KT3Bin=0;// Kt3 bin. 0=low Kt3 bin. 1=high Kt3 bin
36//
37const int CoulChoice=0;// 0 for GRS (default), 1 for Omega0
38
39double MIN[10]={0.97, 0.97, 0.97, 0.97, 0.97, 0.97, 0.97, 0.97, 0.97, 0.97};// C2 y-axis min
3e748846 40double MAX[10]={1.39, 1.38, 1.38, 1.38, 1.38, 1.38, 1.38, 1.38, 1.38, 1.53};// C2 y-axis max
9b479142 41//
9b34d0fd 42const int KTBINS = 6;
43int KTINDEX;
44bool ChargeConstraint=kFALSE;
45bool LinkRadii=kFALSE;
46//
47int TextFont=42;// 63, or 42
3e748846 48float SizeLabel=0.1;// 20(63 font), 0.08(42 font)
49float SizeLegend=0.1;// .08
50float SizeTitle=0.12;//
9b34d0fd 51float SizeSpecif=0.075;//
52float SF1=2/3.*0.95;
53float SF2=1/2.*0.95;
54
3e748846 55double RightMargin=0.004;// 0.002
9b479142 56//
57double Chi2_C2global;
58double NFitPoints_C2global;
59TH1D *C2_ss[KTBINS][10];
60TH1D *C2_os[KTBINS][10];
61TH1D *C2_ss_MomSys[KTBINS][10];
62TH1D *C2_ss_KSys[KTBINS][10];
63TH1D *C2_os_MomSys[KTBINS][10];
64TH1D *C2_os_KSys[KTBINS][10];
9b34d0fd 65TH1D *hfitC2ss_noG[KTBINS][10];
66TH1D *hfitC2os_noG[KTBINS][10];
9b479142 67TF1 *fitC2ss_noG[KTBINS][10];
68TF1 *fitC2os_noG[KTBINS][10];
69TF1 *fitC2ss_yesG[KTBINS][10];
70TF1 *fitC2os_yesG[KTBINS][10];
71TF1 *fitC2ss_noGEWfromTherm[KTBINS][10];
72TF1 *fitC2os_noGEWfromTherm[KTBINS][10];
73TF1 *fitC2ss_yesGEWfromTherm[KTBINS][10];
74TF1 *fitC2os_yesGEWfromTherm[KTBINS][10];
b4f433d0 75TF1 *fitC2ss_noGEW[KTBINS][10];
76TF1 *fitC2os_noGEW[KTBINS][10];
9b34d0fd 77TH1D *hfitC2ss_noGEW[KTBINS][10];
78TH1D *hfitC2os_noGEW[KTBINS][10];
9b479142 79TH1D *K2_ss[KTBINS][10];
80TH1D *K2_os[KTBINS][10];
81TH1D *C2Therm_ss[KTBINS][10];
82TH1D *C2Therm_os[KTBINS][10];
83TF1 *fitC2ss_Therm[KTBINS][10];
84TF1 *fitC2os_Therm[KTBINS][10];
85TF1 *fitC2ss_ThermGaus[KTBINS][10];
86TF1 *fitC2os_ThermGaus[KTBINS][10];
87
88
89void DrawALICELogo(Bool_t, Float_t, Float_t, Float_t, Float_t);
90
91void Plot_plots(){
92
93 gStyle->SetOptStat(0);
94 gStyle->SetOptDate(0);
95
96
97 TFile *files_2_noG[KTBINS][10];
98 TFile *files_2_yesG[KTBINS][10];
99 TFile *files_2_noGEWfromTherm[KTBINS][10];
100 TFile *files_2_yesGEWfromTherm[KTBINS][10];
b4f433d0 101 TFile *files_2_noGEW[KTBINS][10];
9b479142 102 TFile *files_3[2][2][2][2][10];// SC/MC, +/-, GRS/Omega0, Therm/Gauss, MBINS
b4f433d0 103 //
9b479142 104 double intercept1[10]={0};
105 double intercept1_e[10]={0};
106 double intercept2[10]={0};
107 double intercept2_e[10]={0};
108 double ChiNDF1[10]={0};
109 double ChiNDF2[10]={0};
110 double cent_mean[10]={0};
111 double cent_mean_e[10]={0};
112
113
114 TH1D *Chi2NDF_yesG[KTBINS][10];
115 TH1D *Chi2NDF_yesGEWfromTherm[KTBINS][10];
116 TH1D *Chi2NDF_noG[KTBINS][10];
117 TH1D *Chi2NDF_noGEWfromTherm[KTBINS][10];
118 TH1D *Cterm1;
119 TH1D *C3[2][2][2][10];// SC/MC, GRS/Omega0, Therm/Gauss, cb
120 TH1D *c3[2][2][2][10];// SC/MC, GRS/Omega0, Therm/Gauss, cb
121 TH1D *r3_Q3[2][2][10];// GRS/Omega0, Therm/Gauss, cb
122 TH1D *C3EW[2][2][10];// GRS/Omega0, Therm/Gauss, cb
123 TH1D *K3[2][2][2][10];// SC/MC, GRS/Omega0, Therm/Gauss, cb
124 //
125 TH1D *ParHisto_ch[4][10];
126 TH1D *ParHisto_coh[4][10];
127 TH1D *ParHisto_chEWfromTherm[4][10];
128 TH1D *ParHisto_cohEWfromTherm[4][10];
b4f433d0 129 TH1D *ParHisto_chEW[4][10];
130
9b479142 131 for(int ii=0; ii<10; ii++){
132 for(int par=1; par<=4; par++){
133 TString *name_ch = new TString("ParHisto_ch");
134 *name_ch += ii; *name_ch += par;
135 TString *name_coh = new TString("ParHisto_coh");
136 *name_coh += ii; *name_coh += par;
b4f433d0 137 TString *name_EWfromThermch = new TString("ParHisto_EWfromThermch");
138 *name_EWfromThermch += ii; *name_EWfromThermch += par;
9b479142 139 TString *name_EWch = new TString("ParHisto_EWch");
140 *name_EWch += ii; *name_EWch += par;
141 TString *name_EWcoh = new TString("ParHisto_EWcoh");
142 *name_EWcoh += ii; *name_EWcoh += par;
143
144 int MStyle_ch, MStyle_coh;
145 if(ii==BOI_1) {MStyle_ch = 20; MStyle_coh = 24;}
146 else {MStyle_ch = 22; MStyle_coh = 26;}
147
148 ParHisto_ch[par-1][ii] = new TH1D(name_ch->Data(),"",10,0,1);
149 ParHisto_ch[par-1][ii]->SetMarkerStyle(MStyle_ch);
150 ParHisto_ch[par-1][ii]->SetMarkerColor(1);
151 ParHisto_ch[par-1][ii]->SetLineColor(1);
152 ParHisto_ch[par-1][ii]->SetMarkerSize(1.5);
153 ParHisto_coh[par-1][ii] = new TH1D(name_coh->Data(),"",10,0,1);
154 ParHisto_coh[par-1][ii]->SetMarkerStyle(MStyle_coh);
155 ParHisto_coh[par-1][ii]->SetMarkerColor(4);
156 ParHisto_coh[par-1][ii]->SetLineColor(4);
157 ParHisto_coh[par-1][ii]->SetMarkerSize(1.5);
b4f433d0 158 ParHisto_chEWfromTherm[par-1][ii] = new TH1D(name_EWfromThermch->Data(),"",10,0,1);
9b479142 159 ParHisto_chEWfromTherm[par-1][ii]->SetMarkerStyle(MStyle_ch);
160 ParHisto_chEWfromTherm[par-1][ii]->SetMarkerColor(1);
161 ParHisto_chEWfromTherm[par-1][ii]->SetLineColor(1);
162 ParHisto_chEWfromTherm[par-1][ii]->SetMarkerSize(1.5);
b4f433d0 163 ParHisto_chEW[par-1][ii] = new TH1D(name_EWch->Data(),"",10,0,1);
164 ParHisto_chEW[par-1][ii]->SetMarkerStyle(MStyle_coh);
165 ParHisto_chEW[par-1][ii]->SetMarkerColor(4);
166 ParHisto_chEW[par-1][ii]->SetLineColor(4);
167 ParHisto_chEW[par-1][ii]->SetMarkerSize(1.5);
9b479142 168 ParHisto_cohEWfromTherm[par-1][ii] = new TH1D(name_EWcoh->Data(),"",10,0,1);
169 ParHisto_cohEWfromTherm[par-1][ii]->SetMarkerStyle(MStyle_coh);
170 ParHisto_cohEWfromTherm[par-1][ii]->SetMarkerColor(4);
171 ParHisto_cohEWfromTherm[par-1][ii]->SetLineColor(4);
172 ParHisto_cohEWfromTherm[par-1][ii]->SetMarkerSize(1.5);
173 }
174 }
b4f433d0 175 //////////////////////////////
176
177 // Start File access
9b479142 178 for(int cb=0; cb<10; cb++){
179 for(int ChComb=0; ChComb<2; ChComb++) {// SC or MC
b4f433d0 180 for(int ch=0; ch<2; ch++) {// - or +
181 for(int KT3=0; KT3<2; KT3++) {// Kt3 bin
182 TString *name3 = new TString("OutFiles/OutFile");
183 if(ChComb==0) name3->Append("SC");
184 else name3->Append("MC");
185 if(ch==0) name3->Append("Neg");
186 else name3->Append("Pos");
187 TString *name3_Omega0=new TString(name3->Data());
188 name3->Append("NoGEWGRS");
189 name3_Omega0->Append("NoGEWOmega0");
190 name3->Append("Kt3_"); name3_Omega0->Append("Kt3_");
191 *name3 += KT3+1; *name3_Omega0 += KT3+1;
192
193 name3->Append("_Kt10_M");
194 name3_Omega0->Append("_Kt10_M");
195 if(cb<10) {*name3 += cb; *name3_Omega0 += cb;}
196 else {*name3 += 0; *name3_Omega0 += 0;}
197 name3->Append(".root");
198 name3_Omega0->Append(".root");
199 files_3[ChComb][ch][0][KT3][cb] = new TFile(name3->Data(),"READ");
200 files_3[ChComb][ch][1][KT3][cb] = new TFile(name3_Omega0->Data(),"READ");
201 ///////////////////////////////
202 for(int Coul=0; Coul<2; Coul++){
9b479142 203 if(ch==0) {
204 C3[ChComb][Coul][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("C3");
b4f433d0 205 C3[ChComb][Coul][KT3][cb]->SetDirectory(0);
9b479142 206 c3[ChComb][Coul][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("c3");
b4f433d0 207 c3[ChComb][Coul][KT3][cb]->SetDirectory(0);
9b479142 208 if(ChComb==0) C3EW[Coul][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("C3_EWexpectation");
209 if(ChComb==0) r3_Q3[Coul][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("r3_Q3");
b4f433d0 210 C3EW[Coul][KT3][cb]->SetDirectory(0);
211 r3_Q3[Coul][KT3][cb]->SetDirectory(0);
212
9b479142 213 if(Coul==1){
214 K3[ChComb][0][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("Coul_GRiverside");
b4f433d0 215 K3[ChComb][0][KT3][cb]->SetDirectory(0);
9b479142 216 K3[ChComb][1][KT3][cb]=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("Coul_Omega0");
b4f433d0 217 K3[ChComb][1][KT3][cb]->SetDirectory(0);
9b479142 218 }
219 }else{
220 TH1D *tempC3=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("C3");
221 TH1D *tempc3=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("c3");
222 TH1D *tempC3EW=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("C3_EWexpectation");
223 TH1D *tempr3_Q3;
224 if(ChComb==0) tempr3_Q3=(TH1D*)files_3[ChComb][ch][Coul][KT3][cb]->Get("r3_Q3");
225 for(int bin=1; bin<=tempC3->GetNbinsX(); bin++){
226 double valueC3 = (C3[ChComb][Coul][KT3][cb]->GetBinContent(bin) + tempC3->GetBinContent(bin))/2.;
227 double valueC3_e = sqrt(pow(C3[ChComb][Coul][KT3][cb]->GetBinError(bin),2)+pow(tempC3->GetBinError(bin),2))/sqrt(2.);
228 double valuec3 = (c3[ChComb][Coul][KT3][cb]->GetBinContent(bin) + tempc3->GetBinContent(bin))/2.;
229 double valuec3_e = sqrt(pow(c3[ChComb][Coul][KT3][cb]->GetBinError(bin),2)+pow(tempc3->GetBinError(bin),2))/sqrt(2.);
230 C3[ChComb][Coul][KT3][cb]->SetBinContent(bin, valueC3);
231 C3[ChComb][Coul][KT3][cb]->SetBinError(bin, valueC3_e);
232 c3[ChComb][Coul][KT3][cb]->SetBinContent(bin, valuec3);
233 c3[ChComb][Coul][KT3][cb]->SetBinError(bin, valuec3_e);
234 if(Coul==0) {c3[ChComb][Coul][KT3][cb]->SetMarkerStyle(20); c3[ChComb][Coul][KT3][cb]->SetMarkerColor(4); c3[ChComb][Coul][KT3][cb]->SetLineColor(4);}
235 else {c3[ChComb][Coul][KT3][cb]->SetMarkerStyle(22); c3[ChComb][Coul][KT3][cb]->SetMarkerColor(2); c3[ChComb][Coul][KT3][cb]->SetLineColor(2);}
3e748846 236 c3[ChComb][Coul][KT3][cb]->GetYaxis()->SetTitleOffset(1.2);
9b479142 237 if(ChComb==1){
238 c3[ChComb][Coul][KT3][cb]->SetMinimum(0.96);// 0.99
239 c3[ChComb][Coul][KT3][cb]->SetMaximum(1.02);// 1.02
240 c3[ChComb][Coul][KT3][cb]->GetXaxis()->SetRangeUser(0,0.13);
3e748846 241 c3[ChComb][Coul][KT3][cb]->GetYaxis()->SetTitleOffset(1.2);
9b479142 242 c3[ChComb][Coul][KT3][cb]->GetYaxis()->SetTitleSize(.05);
243 }
244
245 if(ChComb==0) {// only same-charge for r3
246 double valuer3_Q3 = (r3_Q3[Coul][KT3][cb]->GetBinContent(bin) + tempr3_Q3->GetBinContent(bin))/2.;
247 double valuer3_Q3_e = sqrt(pow(r3_Q3[Coul][KT3][cb]->GetBinError(bin),2) + pow(tempr3_Q3->GetBinError(bin),2))/sqrt(2.);
248 double valueC3EW = (C3EW[Coul][KT3][cb]->GetBinContent(bin) + tempC3EW->GetBinContent(bin))/2.;
249 double valueC3EW_e = sqrt(pow(C3EW[Coul][KT3][cb]->GetBinError(bin),2)+pow(tempC3EW->GetBinError(bin),2))/sqrt(2.);
250 C3EW[Coul][KT3][cb]->SetBinContent(bin, valueC3EW);
251 C3EW[Coul][KT3][cb]->SetBinError(bin, valueC3EW_e);
252 r3_Q3[Coul][KT3][cb]->SetBinContent(bin, valuer3_Q3);
253 r3_Q3[Coul][KT3][cb]->SetBinError(bin, valuer3_Q3_e);
254 r3_Q3[Coul][KT3][cb]->GetXaxis()->SetRangeUser(0,0.14);
255 if(cb<2) r3_Q3[Coul][KT3][cb]->SetMinimum(-3);
256 else r3_Q3[Coul][KT3][cb]->SetMinimum(-0.5);
257 r3_Q3[Coul][KT3][cb]->SetMaximum(2.6);
258 if(Coul==0) {r3_Q3[Coul][KT3][cb]->SetMarkerStyle(20); r3_Q3[Coul][KT3][cb]->SetMarkerColor(4); r3_Q3[Coul][KT3][cb]->SetLineColor(4);}
259 else {r3_Q3[Coul][KT3][cb]->SetMarkerStyle(22); r3_Q3[Coul][KT3][cb]->SetMarkerColor(2); r3_Q3[Coul][KT3][cb]->SetLineColor(2);}
260 }
261 }
9b479142 262 }
b4f433d0 263 files_3[ChComb][ch][Coul][KT3][cb]->Close();
264 }// Coul
265 }
266 }
267 }
268
269
9b479142 270 for(int kt=0; kt<KTBINS; kt++){
271
b4f433d0 272 TString *name = new TString("OutFiles/OutFileSCPosNoG");
273 TString *nameEWfromTherm = new TString(name->Data()); nameEWfromTherm->Append("EWfromThermGRS");
274 TString *nameEW = new TString(name->Data()); nameEW->Append("EWGRS");
275 name->Append("GRS");
276 name->Append("Kt3_1"); nameEWfromTherm->Append("Kt3_1"); nameEW->Append("Kt3_1");
277 name->Append("_Kt"); nameEWfromTherm->Append("_Kt"); nameEW->Append("_Kt");
278 *name += kt+1; *nameEWfromTherm += kt+1; *nameEW += kt+1;
279 name->Append("_M"); nameEWfromTherm->Append("_M"); nameEW->Append("_M");
280 if(cb<10) {*name += cb; *nameEWfromTherm += cb; *nameEW += cb;}
281 else {*name += 0; *nameEWfromTherm += 0; *nameEW += 0;}
282 name->Append(".root"); nameEWfromTherm->Append(".root"); nameEW->Append(".root");
283 files_2_noG[kt][cb] = new TFile(name->Data(),"READ");
284 files_2_noGEWfromTherm[kt][cb] = new TFile(nameEWfromTherm->Data(),"READ");
285 files_2_noGEW[kt][cb] = new TFile(nameEW->Data(),"READ");
286 //
287
288 //
289 TString *nameYesG = new TString("OutFiles/OutFileSCPosYesG");
290 TString *nameYesGEWfromTherm = new TString(nameYesG->Data()); nameYesGEWfromTherm->Append("EWfromThermGRS");
291 nameYesG->Append("GRS");
292 nameYesG->Append("Kt3_1"); nameYesGEWfromTherm->Append("Kt3_1");
293 nameYesG->Append("_Kt"); nameYesGEWfromTherm->Append("_Kt");
294 *nameYesG += kt+1; *nameYesGEWfromTherm += kt+1;
295 nameYesG->Append("_M"); nameYesGEWfromTherm->Append("_M");
296 if(cb<10) {*nameYesG += cb; *nameYesGEWfromTherm += cb;}
297 else {*nameYesG += 0; *nameYesGEWfromTherm += 0;}
298 nameYesG->Append(".root"); nameYesGEWfromTherm->Append(".root");
299 files_2_yesG[kt][cb] = new TFile(nameYesG->Data(),"READ");
300 files_2_yesGEWfromTherm[kt][cb] = new TFile(nameYesGEWfromTherm->Data(),"READ");
301 //////////////////////////////////////////
302
9b479142 303 TString *Centname = new TString("Centrality ");
304 *Centname += (cb)*5;
305 Centname->Append("-");
306 *Centname += (cb+1)*5;
307 Centname->Append("%");
308
309
b4f433d0 310 C2_ss[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_ss"); C2_ss[kt][cb]->SetDirectory(0);
311 C2_os[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_os"); C2_os[kt][cb]->SetDirectory(0);
312 C2_ss_MomSys[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_ss_Momsys"); C2_ss_MomSys[kt][cb]->SetDirectory(0);
313 C2_os_MomSys[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_os_Momsys"); C2_os_MomSys[kt][cb]->SetDirectory(0);
314 C2_ss_KSys[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_ss_Ksys"); C2_ss_KSys[kt][cb]->SetDirectory(0);
315 C2_os_KSys[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2_os_Ksys"); C2_os_KSys[kt][cb]->SetDirectory(0);
9b34d0fd 316 fitC2ss_noG[kt][cb] = (TF1*)files_2_noG[kt][cb]->Get("fitC2ss");
317 fitC2os_noG[kt][cb] = (TF1*)files_2_noG[kt][cb]->Get("fitC2os");
318 hfitC2ss_noG[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("fitC2ss_h"); hfitC2ss_noG[kt][cb]->SetDirectory(0);
319 hfitC2os_noG[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("fitC2os_h"); hfitC2os_noG[kt][cb]->SetDirectory(0);
320 fitC2ss_yesG[kt][cb] = (TF1*)files_2_yesG[kt][cb]->Get("fitC2ss");
321 fitC2os_yesG[kt][cb] = (TF1*)files_2_yesG[kt][cb]->Get("fitC2os");
322 fitC2ss_noGEW[kt][cb] = (TF1*)files_2_noGEW[kt][cb]->Get("fitC2ss");
323 fitC2os_noGEW[kt][cb] = (TF1*)files_2_noGEW[kt][cb]->Get("fitC2os");
324 hfitC2ss_noGEW[kt][cb] = (TH1D*)files_2_noGEW[kt][cb]->Get("fitC2ss_h"); hfitC2ss_noGEW[kt][cb]->SetDirectory(0);
325 hfitC2os_noGEW[kt][cb] = (TH1D*)files_2_noGEW[kt][cb]->Get("fitC2os_h"); hfitC2os_noGEW[kt][cb]->SetDirectory(0);
326 fitC2ss_noGEWfromTherm[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2ss");
327 fitC2os_noGEWfromTherm[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2os");
328 fitC2ss_yesGEWfromTherm[kt][cb] = (TF1*)files_2_yesGEWfromTherm[kt][cb]->Get("fitC2ss");
329 fitC2os_yesGEWfromTherm[kt][cb] = (TF1*)files_2_yesGEWfromTherm[kt][cb]->Get("fitC2os");
330 fitC2ss_Therm[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2ssTherm");
331 fitC2os_Therm[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2osTherm");
332 fitC2ss_ThermGaus[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2ssThermGaus");
333 fitC2os_ThermGaus[kt][cb] = (TF1*)files_2_noGEWfromTherm[kt][cb]->Get("fitC2osThermGaus");
b4f433d0 334 K2_ss[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("K2_ss"); K2_ss[kt][cb]->SetDirectory(0);
335 K2_os[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("K2_os"); K2_os[kt][cb]->SetDirectory(0);
336 C2Therm_ss[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2Therm_ss"); C2Therm_ss[kt][cb]->SetDirectory(0);
337 C2Therm_os[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("C2Therm_os"); C2Therm_os[kt][cb]->SetDirectory(0);
338 Chi2NDF_noG[kt][cb] = (TH1D*)files_2_noG[kt][cb]->Get("ChiSquaredNDF"); Chi2NDF_noG[kt][cb]->SetDirectory(0);
339 Chi2NDF_noGEWfromTherm[kt][cb] = (TH1D*)files_2_noGEWfromTherm[kt][cb]->Get("ChiSquaredNDF"); Chi2NDF_noGEWfromTherm[kt][cb]->SetDirectory(0);
340 Chi2NDF_yesG[kt][cb] = (TH1D*)files_2_yesG[kt][cb]->Get("ChiSquaredNDF"); Chi2NDF_yesG[kt][cb]->SetDirectory(0);
341 Chi2NDF_yesGEWfromTherm[kt][cb] = (TH1D*)files_2_yesGEWfromTherm[kt][cb]->Get("ChiSquaredNDF"); Chi2NDF_yesGEWfromTherm[kt][cb]->SetDirectory(0);
9b34d0fd 342
9b479142 343 //
344 for(int par=1; par<=4; par++){
345 //if(par!=4 && par!=2)
346 ParHisto_ch[par-1][cb]->SetBinContent(kt+3, fitC2ss_noG[kt][cb]->GetParameter(par));
347 ParHisto_ch[par-1][cb]->SetBinError(kt+3, fitC2ss_noG[kt][cb]->GetParError(par));
348 ParHisto_coh[par-1][cb]->SetBinContent(kt+3, fitC2ss_yesG[kt][cb]->GetParameter(par));
349 ParHisto_coh[par-1][cb]->SetBinError(kt+3, fitC2ss_yesG[kt][cb]->GetParError(par));
350 //if(par!=4 && par!=2)
351 ParHisto_chEWfromTherm[par-1][cb]->SetBinContent(kt+3, fitC2ss_noGEWfromTherm[kt][cb]->GetParameter(par));
352 ParHisto_chEWfromTherm[par-1][cb]->SetBinError(kt+3, fitC2ss_noGEWfromTherm[kt][cb]->GetParError(par));
353 ParHisto_cohEWfromTherm[par-1][cb]->SetBinContent(kt+3, fitC2ss_yesGEWfromTherm[kt][cb]->GetParameter(par));
354 ParHisto_cohEWfromTherm[par-1][cb]->SetBinError(kt+3, fitC2ss_yesGEWfromTherm[kt][cb]->GetParError(par));
b4f433d0 355 ParHisto_chEW[par-1][cb]->SetBinContent(kt+3, fitC2ss_noGEW[kt][cb]->GetParameter(par));
356 ParHisto_chEW[par-1][cb]->SetBinError(kt+3, fitC2ss_noGEW[kt][cb]->GetParError(par));
9b479142 357 }
b4f433d0 358 //////////////////////////////////////
359 files_2_noG[kt][cb]->Close();
360 files_2_noGEWfromTherm[kt][cb]->Close();
361 files_2_noGEW[kt][cb]->Close();
362 files_2_yesG[kt][cb]->Close();
363 files_2_yesGEWfromTherm[kt][cb]->Close();
364 }// kt loop
365 }
9b479142 366
367
b4f433d0 368 TF1 *Unity = new TF1("Unity","1",0,100);
369 Unity->SetLineStyle(2);
370 Unity->SetLineColor(1);
371
372
373
9b479142 374
375 // merge C3 histogram centralities
376 TH1D *C3merged[2][2][2][6];// SC/MC, GRS/Omega0, Therm/Gauss, cb_merged
377 TH1D *c3merged[2][2][2][6];// SC/MC, GRS/Omega0, Therm/Gauss, cb_merged
378 TH1D *r3merged[2][2][6];// GRS/Omega0, Therm/Gauss, cb_merged
379 for(int Coul=0; Coul<2; Coul++){// GRS or Omega0
380 for(int ChComb=0; ChComb<2; ChComb++) {
381 for(int KT3=0; KT3<2; KT3++) {// Therminator or Gaussian FSI source
382 for(int cb=0; cb<6; cb++){
383 C3merged[ChComb][Coul][KT3][cb]=(TH1D*)C3[ChComb][Coul][KT3][cb]->Clone();
384 c3merged[ChComb][Coul][KT3][cb]=(TH1D*)c3[ChComb][Coul][KT3][cb]->Clone();
385 if(ChComb==0) r3merged[Coul][KT3][cb]=(TH1D*)r3_Q3[Coul][KT3][cb]->Clone();
386 }
387 }
388 }
389 }
390
391 for(int MergedBin=0; MergedBin<4; MergedBin++){
392 for(int Coul=0; Coul<2; Coul++){// GRS or Omega0
393 for(int ChComb=0; ChComb<2; ChComb++) {
394 for(int KT3=0; KT3<2; KT3++) {// Therminator or Gaussian FSI source
395 for(int bin=1; bin<=C3[0][Coul][KT3][0]->GetNbinsX(); bin++){
396 int startbin = MergedBin*2 + 2;
397
398 double valueC3 = C3[ChComb][Coul][KT3][startbin]->GetBinContent(bin) + C3[ChComb][Coul][KT3][startbin+1]->GetBinContent(bin);
399 valueC3 /= 2.;
400 double valueC3_e = pow(C3[ChComb][Coul][KT3][startbin]->GetBinError(bin),2) + pow(C3[ChComb][Coul][KT3][startbin+1]->GetBinError(bin),2);
401 valueC3_e = sqrt(valueC3_e);
402 valueC3_e /= 2.;
403 double valuec3 = c3[ChComb][Coul][KT3][startbin]->GetBinContent(bin) + c3[ChComb][Coul][KT3][startbin+1]->GetBinContent(bin);
404 valuec3 /= 2.;
405 double valuec3_e = pow(c3[ChComb][Coul][KT3][startbin]->GetBinError(bin),2) + pow(c3[ChComb][Coul][KT3][startbin+1]->GetBinError(bin),2);
406 valuec3_e = sqrt(valuec3_e);
407 valuec3_e /= 2.;
408 //
409 C3merged[ChComb][Coul][KT3][MergedBin+2]->SetBinContent(bin, valueC3);
410 C3merged[ChComb][Coul][KT3][MergedBin+2]->SetBinError(bin, valueC3_e);
411 c3merged[ChComb][Coul][KT3][MergedBin+2]->SetBinContent(bin, valuec3);
412 c3merged[ChComb][Coul][KT3][MergedBin+2]->SetBinError(bin, valuec3_e);
413 if(ChComb==0){
414 double valueC3EW = C3EW[Coul][KT3][startbin]->GetBinContent(bin) + C3EW[Coul][KT3][startbin+1]->GetBinContent(bin);
415 valueC3EW /= 2.;
416 double valueC3EW_e = pow(C3EW[Coul][KT3][startbin]->GetBinError(bin),2) + pow(C3EW[Coul][KT3][startbin+1]->GetBinError(bin),2);
417 valueC3EW_e = sqrt(valueC3EW_e);
418 valueC3EW_e /= 2.;
419 double valuer3_Q3 = r3_Q3[Coul][KT3][startbin]->GetBinContent(bin) + r3_Q3[Coul][KT3][startbin+1]->GetBinContent(bin);
420 valuer3_Q3 /= 2.;
421 double valuer3_Q3_e = pow(r3_Q3[Coul][KT3][startbin]->GetBinError(bin),2) + pow(r3_Q3[Coul][KT3][startbin+1]->GetBinError(bin),2);
422 valuer3_Q3_e = sqrt(valuer3_Q3_e);
423 valuer3_Q3_e /= 2.;
424 r3merged[Coul][KT3][MergedBin+2]->SetBinContent(bin, valuer3_Q3);
425 r3merged[Coul][KT3][MergedBin+2]->SetBinError(bin, valuer3_Q3_e);
426 C3EW[Coul][KT3][MergedBin+2]->SetBinContent(bin, valueC3EW);
427 C3EW[Coul][KT3][MergedBin+2]->SetBinError(bin, valueC3EW_e);
428 }
429 //
430 }
431 }// KT3
432 }
433 }
434 }// Merged Bin
435
436
437
438 //////////////////////////////////////////////////////////////////////////
439 // Draw K3 factors
440 TCanvas *canK3 = new TCanvas("canK3", "canK3",10,0,600,900);// 11,53,700,500
441 canK3->SetHighLightColor(2);
442 //canK3->Range(-0.7838115,-1.033258,0.7838115,1.033258);
443 gStyle->SetOptFit(0111);
3e748846 444 canK3->SetFillColor(0);//10
9b479142 445 canK3->SetBorderMode(0);
446 canK3->SetBorderSize(2);
447 canK3->SetFrameFillColor(0);
448 canK3->SetFrameBorderMode(0);
449 canK3->SetFrameBorderMode(0);
450
451 TPad *padK3 = new TPad("padK3","padK3",0.0,0.0,1.,1.);
452 gPad->SetGridx(0);
453 gPad->SetGridy(0);
454 gPad->SetTickx();
455 gPad->SetTicky();
456 padK3->SetTopMargin(0.02);//0.05
457 padK3->SetRightMargin(0.01);//1e-2
9b34d0fd 458 padK3->SetBottomMargin(0.07);//0.12
9b479142 459 padK3->Divide(1,2,0,0);
460 padK3->Draw();
461 padK3->cd(1);
3e748846 462 gPad->SetLeftMargin(0.14);
463 gPad->SetRightMargin(0.03);
9b34d0fd 464 //double Dim1=gPad->GetAbsHNDC();
465 //cout<<gPad->GetAbsHNDC()<<endl;
466 //cout<<gPad->GetAspectRatio()<<endl;
9b479142 467 TLegend *legendK3 = new TLegend(.45,.05,.95,.3,NULL,"brNDC");
9b34d0fd 468 legendK3->SetBorderSize(0);
469 legendK3->SetTextSize(SizeLabel*SF1);// *2/3.*0.96
9b479142 470 legendK3->SetFillColor(0);
9b479142 471
472 K3[0][1][KT3Bin][BOI_1]->SetMinimum(0.48);
473 K3[0][1][KT3Bin][BOI_1]->SetMaximum(1.33);
474
475 K3[0][1][KT3Bin][BOI_1]->SetMarkerStyle(20);
476 K3[0][1][KT3Bin][BOI_1]->SetMarkerColor(2);
477 K3[0][0][KT3Bin][BOI_1]->SetMarkerStyle(24);
478 K3[0][0][KT3Bin][BOI_1]->SetMarkerColor(2);
479 K3[1][1][KT3Bin][BOI_1]->SetMarkerStyle(22);
480 K3[1][1][KT3Bin][BOI_1]->SetMarkerColor(4);
481 K3[1][0][KT3Bin][BOI_1]->SetMarkerStyle(26);
482 K3[1][0][KT3Bin][BOI_1]->SetMarkerColor(4);
483 //
9b34d0fd 484 K3[0][1][KT3Bin][BOI_1]->GetXaxis()->SetLabelFont(TextFont); K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetLabelFont(TextFont);
485 K3[0][1][KT3Bin][BOI_1]->GetXaxis()->SetLabelSize(SizeLabel*SF1); K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetLabelSize(SizeLabel*SF1);
486 K3[0][1][KT3Bin][BOI_1]->GetXaxis()->SetNdivisions(808); K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetNdivisions(909);
487 TGaxis::SetMaxDigits(3);
488 K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetTitleFont(TextFont);
3e748846 489 K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetTitleSize(SizeTitle*SF1);
490 K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetTitle("#font[12]{K_{3}}");
491 K3[0][1][KT3Bin][BOI_1]->GetYaxis()->SetTitleOffset(0.8);
492
9b479142 493 //
494 K3[0][1][KT3Bin][BOI_1]->Draw();
495 K3[0][0][KT3Bin][BOI_1]->Draw("same");
496 K3[1][1][KT3Bin][BOI_1]->Draw("same");
497 K3[1][0][KT3Bin][BOI_1]->Draw("same");
498 Unity->Draw("same");
499 legendK3->AddEntry(K3[0][1][KT3Bin][BOI_1],"same-charge, #Omega_{0}","p");
500 legendK3->AddEntry(K3[0][0][KT3Bin][BOI_1],"same-charge, GRS","p");
501 legendK3->AddEntry(K3[1][1][KT3Bin][BOI_1],"mixed-charge, #Omega_{0}","p");
502 legendK3->AddEntry(K3[1][0][KT3Bin][BOI_1],"mixed-charge, GRS","p");
503 legendK3->Draw("same");
3e748846 504 //TLatex *K3Label = new TLatex(0.02,1.25,"K_{3}");// -0.011,0.92 (ss), 1.26 (os)
505 //K3Label->SetTextFont(TextFont);
506 //K3Label->SetTextSize(SizeTitle*SF1);
507 //K3Label->SetTextAngle(90);
9b34d0fd 508 //K3Label->Draw();
9b479142 509 //
510 padK3->cd(2);
3e748846 511 double SF_correction=0.97;
512 gPad->SetLeftMargin(0.14);
513 gPad->SetRightMargin(0.03);
514 gPad->SetBottomMargin(0.16);
9b34d0fd 515 //double Dim2=gPad->GetAbsHNDC();
516 TLegend *legendK3comp = new TLegend(.45,.8,.95,.95,NULL,"brNDC");
517 legendK3comp->SetBorderSize(0);
3e748846 518 legendK3comp->SetTextSize(SizeLabel*SF1*SF_correction);// small .03; large .06
9b34d0fd 519 legendK3comp->SetFillColor(0);
520
9b479142 521 TH1D *K3sc_compOmega0 = (TH1D*)K3[0][1][KT3Bin][BOI_1]->Clone();
522 TH1D *K3sc_compGRS = (TH1D*)K3[0][0][KT3Bin][BOI_1]->Clone();
523 TH1D *K3mc_compOmega0 = (TH1D*)K3[1][1][KT3Bin][BOI_1]->Clone();
524 TH1D *K3mc_compGRS = (TH1D*)K3[1][0][KT3Bin][BOI_1]->Clone();
525 for(int ii=0; ii<K3sc_compOmega0->GetNbinsX(); ii++){
526 K3sc_compOmega0->SetBinContent(ii+1, K3sc_compOmega0->GetBinContent(ii+1)-1.0);
527 K3sc_compGRS->SetBinContent(ii+1, K3sc_compGRS->GetBinContent(ii+1)-1.0);
528 K3mc_compOmega0->SetBinContent(ii+1, K3mc_compOmega0->GetBinContent(ii+1)-1.0);
529 K3mc_compGRS->SetBinContent(ii+1, K3mc_compGRS->GetBinContent(ii+1)-1.0);
530 }
531 K3sc_compOmega0->Add(K3sc_compGRS,-1);
532 K3sc_compOmega0->Divide(K3sc_compGRS);
533 K3mc_compOmega0->Add(K3mc_compGRS,-1);
534 K3mc_compOmega0->Divide(K3mc_compGRS);
535 K3sc_compOmega0->SetMarkerStyle(20);
536 K3sc_compOmega0->SetMarkerColor(2);
537 K3sc_compOmega0->SetLineColor(2);
538 K3mc_compOmega0->SetMarkerStyle(24);
539 K3mc_compOmega0->SetMarkerColor(4);
540 K3mc_compOmega0->SetLineColor(4);
541 K3sc_compOmega0->SetMinimum(-0.12);// -0.021
542 K3sc_compOmega0->SetMaximum(0.12);// 0.021
543 K3sc_compOmega0->SetBinContent(1,-100); K3mc_compOmega0->SetBinContent(1,-100);
9b34d0fd 544 K3sc_compOmega0->GetYaxis()->SetNdivisions(404);
3e748846 545 K3sc_compOmega0->GetYaxis()->SetTitleFont(TextFont); K3sc_compOmega0->GetYaxis()->SetTitleSize(SizeTitle*SF1*SF_correction);
9b34d0fd 546 //K3sc_compOmega0->GetXaxis()->SetTitleFont(TextFont); K3sc_compOmega0->GetXaxis()->SetTitleSize(SizeTitle*2/3.*.96);
547 //K3sc_compOmega0->GetXaxis()->SetTitle("Q_{3} (GeV/c)");
3e748846 548 K3sc_compOmega0->GetYaxis()->SetTitle("#Delta#font[12]{K_{3} / (K_{3}}(#Omega_{0})-1)");
549 K3sc_compOmega0->GetXaxis()->SetTitleOffset(1.2); K3sc_compOmega0->GetYaxis()->SetTitleOffset(0.8);
9b479142 550 K3sc_compOmega0->Draw();
551 K3mc_compOmega0->Draw("same");
552 TF1 *Zero=new TF1("Zero","0",0,1);
553 Zero->SetLineStyle(2); Zero->SetLineColor(1);
554 Zero->Draw("same");
555 legendK3comp->AddEntry(K3sc_compOmega0,"same-charge","p");
556 legendK3comp->AddEntry(K3mc_compOmega0,"mixed-charge","p");
557 legendK3comp->Draw("same");
3e748846 558 //TLatex *RatioLabel = new TLatex(-.011,0.02,"#Delta K_{3}/(K_{3}(#Omega_{0})-1)");// -0.011,0.04
559 //RatioLabel->SetTextFont(TextFont);
560 //RatioLabel->SetTextSize(SizeTitle*SF1);
561 //RatioLabel->SetTextAngle(90);
9b34d0fd 562 //RatioLabel->Draw();
563
3e748846 564 TLatex *Q3Label = new TLatex(.071,-0.155,"#font[12]{Q}_{3} (GeV/#font[12]{c})");//.065,-0.147
9b34d0fd 565 Q3Label->SetTextFont(TextFont);
566 Q3Label->SetTextSize(SizeTitle*SF1);// 0.08
9b479142 567 Q3Label->Draw();
568
569
570
571
572
573 //////////////////////////////////////////
574 // print global fit values
575
576 // kappa3 and kappa4
b4f433d0 577 cout<<"ALICE EW"<<endl;
578 cout<<"kappa3"<<endl;
579 for(int cb=0; cb<10; cb++){
580 if(cb!=0 && cb!=9) continue;
581 TString *Cname = new TString("$");
582 *Cname += cb*5; Cname->Append("-");
583 *Cname += (cb+1)*5; Cname->Append("\\%$");
584 cout<<Cname->Data()<<" & ";
585 cout.precision(2);
586 for(int kt=0; kt<KTBINS; kt++) {
587 cout<<fitC2ss_noGEW[kt][cb]->GetParameter(5)<<" ";
588 if((kt+1)!=KTBINS) cout<<"& ";
589 else cout<<" \\\\\ \\hline";
590 }
591 cout<<endl;
592 }
593 //
594 cout<<"kappa4"<<endl;
595 for(int cb=0; cb<10; cb++){
596 if(cb!=0 && cb!=9) continue;
597 TString *Cname = new TString("$");
598 *Cname += cb*5; Cname->Append("-");
599 *Cname += (cb+1)*5; Cname->Append("\\%$");
600 cout<<Cname->Data()<<" & ";
601 cout.precision(2);
602 for(int kt=0; kt<KTBINS; kt++) {
603 cout<<fitC2ss_noGEW[kt][cb]->GetParameter(6)<<" ";
604 if((kt+1)!=KTBINS) cout<<"& ";
605 else cout<<" \\\\\ \\hline";
606 }
607 cout<<endl;
608 }
609 // Therminator
610 cout<<"Therminator EW"<<endl;
9b479142 611 cout<<"kappa3"<<endl;
612 for(int cb=0; cb<10; cb++){
613 if(cb!=0 && cb!=9) continue;
614 TString *Cname = new TString("$");
615 *Cname += cb*5; Cname->Append("-");
616 *Cname += (cb+1)*5; Cname->Append("\\%$");
617 cout<<Cname->Data()<<" & ";
618 cout.precision(2);
619 for(int kt=0; kt<KTBINS; kt++) {
620 cout<<fitC2ss_noGEWfromTherm[kt][cb]->GetParameter(5)<<" ";
621 if((kt+1)!=KTBINS) cout<<"& ";
622 else cout<<" \\\\\ \\hline";
623 }
624 cout<<endl;
625 }
626 //
627 cout<<"kappa4"<<endl;
628 for(int cb=0; cb<10; cb++){
629 if(cb!=0 && cb!=9) continue;
630 TString *Cname = new TString("$");
631 *Cname += cb*5; Cname->Append("-");
632 *Cname += (cb+1)*5; Cname->Append("\\%$");
633 cout<<Cname->Data()<<" & ";
634 cout.precision(2);
635 for(int kt=0; kt<KTBINS; kt++) {
636 cout<<fitC2ss_noGEWfromTherm[kt][cb]->GetParameter(6)<<" ";
637 if((kt+1)!=KTBINS) cout<<"& ";
638 else cout<<" \\\\\ \\hline";
639 }
640 cout<<endl;
641 }
642
b4f433d0 643
644
645
9b479142 646 // Chi2/NDF
b4f433d0 647 /*cout<<"Chi2/NDF"<<endl;
9b479142 648 cout<<"fits with Coherence"<<endl;
649 for(int cb=0; cb<10; cb++){
650 if(cb!=0 && cb!=9) continue;
651 TString *Cname = new TString("$");
652 *Cname += cb*5; Cname->Append("-");
653 *Cname += (cb+1)*5; Cname->Append("\\%$");
654 cout<<Cname->Data()<<" & ";
655 cout.precision(2);
656 for(int kt=0; kt<KTBINS; kt++) {
657 //cout<<fitC2ss_yesGEWfromTherm[kt][cb]->GetParameter(4)<<" ";
658 //cout<<fitC2ss_yesG[kt][cb]->GetParameter(4)<<" ";
659 //cout<<fitC2ss_noGEWfromTherm[kt][cb]->GetParameter(3)<<" ";
660 cout<<Chi2NDF_yesG[kt][cb]->GetBinContent(1)<<" ";// Chi^2/NDF
661 if((kt+1)!=KTBINS) cout<<"& ";
662 else cout<<" \\\\\ \\hline";
663 }
664 cout<<endl;
665 }
666 //
667 cout<<"fits with Therminator EW"<<endl;
668 for(int cb=0; cb<10; cb++){
669 if(cb!=0 && cb!=9) continue;
670 TString *Cname = new TString("$");
671 *Cname += cb*5; Cname->Append("-");
672 *Cname += (cb+1)*5; Cname->Append("\\%$");
673 cout<<Cname->Data()<<" & ";
674 cout.precision(2);
675 for(int kt=0; kt<KTBINS; kt++) {
676 cout<<Chi2NDF_noGEWfromTherm[kt][cb]->GetBinContent(1)<<" ";// Chi^2/NDF
677 if((kt+1)!=KTBINS) cout<<"& ";
678 else cout<<" \\\\\ \\hline";
679 }
680 cout<<endl;
681 }
b4f433d0 682 */
9b479142 683
684
9b34d0fd 685
9b479142 686 /////////////////////////////////////////////
687 // Draw C2 in 6 kt bins
688 TCanvas *canC2 = new TCanvas("canC2", "canC2",10,0,600,900);// 11,53,700,500
689 canC2->SetHighLightColor(2);
690 //canC2->Range(-0.7838115,-1.033258,0.7838115,1.033258);
691 gStyle->SetOptFit(0111);
3e748846 692 canC2->SetFillColor(0);//10
9b479142 693 canC2->SetBorderMode(0);
694 canC2->SetBorderSize(2);
695 canC2->SetFrameFillColor(0);
696 canC2->SetFrameBorderMode(0);
697 canC2->SetFrameBorderMode(0);
698
699
3e748846 700 TPad *pad1 = new TPad("pad1","pad1",0.06,0.06,1.,1.);//0.05,0.05,1.,1.
9b479142 701 gPad->SetGridx(0);
702 gPad->SetGridy(0);
703 gPad->SetTickx();
704 gPad->SetTicky();
705 pad1->SetTopMargin(0.02);//0.05
706 pad1->SetLeftMargin(0.12);//
707 pad1->SetRightMargin(0.01);//3e-2
3e748846 708 pad1->SetBottomMargin(0.08);//0.12
9b479142 709 pad1->Divide(2,3,0,0);
710 pad1->Draw();
9b34d0fd 711
9b479142 712
3e748846 713 TLegend *legend1 = new TLegend(.2,.75,.99,.95,NULL,"brNDC");//.45 or .4 for x1
9b34d0fd 714 legend1->SetBorderSize(0);
715 //legend1->SetTextSize(.08);// small .03; large .036
716 legend1->SetTextFont(TextFont);
717 legend1->SetTextSize(SizeLegend);
9b479142 718 legend1->SetFillColor(0);
719 //
3e748846 720 TLegend *legend2 = new TLegend(.35,.74, .99,.98,NULL,"brNDC");//.45 or .4 for x1
9b34d0fd 721 legend2->SetBorderSize(0);
722 //legend2->SetTextSize(.08);// small .03; large .036
723 legend2->SetTextFont(TextFont);
724 legend2->SetTextSize(SizeLegend);
9b479142 725 legend2->SetFillColor(0);
726
727
9b34d0fd 728 TGaxis::SetMaxDigits(4);
9b479142 729 /////////////////////
730 for(int kt=0; kt<KTBINS; kt++) {
731 pad1->cd(kt+1);
732 gPad->SetRightMargin(RightMargin);
3e748846 733 if(kt%2==0) gPad->SetRightMargin(0.001);
734 if(kt==4) SF_correction=0.92;
735 else if(kt==5) SF_correction=1.015;
736 else SF_correction=1.0;
9b34d0fd 737 //if(kt==0) Dim1=gPad->GetAbsWNDC()*gPad->GetAbsHNDC();
738 //Dim2=gPad->GetAbsWNDC()*gPad->GetAbsHNDC();
9b479142 739 C2_ss[kt][BOI_1]->SetMinimum(MIN[BOI_1]);
740 C2_ss[kt][BOI_1]->SetMaximum(MAX[BOI_1]);
9b34d0fd 741 C2_ss[kt][BOI_1]->GetXaxis()->SetRangeUser(0,0.061);
9b479142 742 C2_ss[kt][BOI_1]->GetXaxis()->SetTitleOffset(0.85);
9b34d0fd 743 C2_ss[kt][BOI_1]->GetYaxis()->SetTitleOffset(0.85);// 0.85
9b479142 744 C2_ss[kt][BOI_1]->GetXaxis()->SetTitle(""); C2_ss[kt][BOI_1]->GetYaxis()->SetTitle("");
9b34d0fd 745 C2_ss[kt][BOI_1]->GetXaxis()->SetTitleFont(TextFont); C2_ss[kt][BOI_1]->GetYaxis()->SetTitleFont(TextFont);
746 C2_ss[kt][BOI_1]->GetXaxis()->SetLabelFont(TextFont); C2_ss[kt][BOI_1]->GetYaxis()->SetLabelFont(TextFont);
3e748846 747 C2_ss[kt][BOI_1]->GetXaxis()->SetLabelSize(SizeLabel*SF_correction); C2_ss[kt][BOI_1]->GetYaxis()->SetLabelSize(SizeLabel*SF_correction);
748 if(kt%2==0) C2_ss[kt][BOI_1]->GetYaxis()->SetLabelSize(SizeLabel*SF_correction);
9b479142 749 else {C2_ss[kt][BOI_1]->GetYaxis()->SetLabelSize(.0);}
750 if(kt<KTBINS-2) C2_ss[kt][BOI_1]->GetXaxis()->SetLabelSize(.0);
3e748846 751 C2_ss[kt][BOI_1]->GetXaxis()->SetNdivisions(603); C2_ss[kt][BOI_1]->GetYaxis()->SetNdivisions(606);
9b34d0fd 752
9b479142 753 //
754 C2Therm_ss[kt][BOI_1]->SetMinimum(MIN[BOI_1]);
755 C2Therm_ss[kt][BOI_1]->SetMaximum(MAX[BOI_1]);
756 C2Therm_ss[kt][BOI_1]->GetXaxis()->SetRangeUser(0,0.061);
757 C2Therm_ss[kt][BOI_1]->GetXaxis()->SetTitleOffset(0.85);
758 C2Therm_ss[kt][BOI_1]->GetYaxis()->SetTitleOffset(0.85);
759 C2Therm_ss[kt][BOI_1]->GetXaxis()->SetTitle(""); C2Therm_ss[kt][BOI_1]->GetYaxis()->SetTitle("");
9b34d0fd 760 C2Therm_ss[kt][BOI_1]->GetXaxis()->SetLabelFont(TextFont); C2Therm_ss[kt][BOI_1]->GetYaxis()->SetLabelFont(TextFont);
3e748846 761 if(kt%2==0) C2Therm_ss[kt][BOI_1]->GetYaxis()->SetLabelSize(SizeLabel*SF_correction);
9b479142 762 else {C2Therm_ss[kt][BOI_1]->GetYaxis()->SetLabelSize(.0);}
763 if(kt<KTBINS-2) C2Therm_ss[kt][BOI_1]->GetXaxis()->SetLabelSize(.0);
3e748846 764 else {C2Therm_ss[kt][BOI_1]->GetXaxis()->SetLabelSize(SizeLabel*SF_correction);}
9b34d0fd 765 C2Therm_ss[kt][BOI_1]->GetXaxis()->SetNdivisions(606); C2Therm_ss[kt][BOI_1]->GetYaxis()->SetNdivisions(606);
9b479142 766
767 if(BOI_1==9 && kt==5) {// Point with large statistical error
768 C2_ss[kt][BOI_1]->SetBinContent(3,-100);
769 C2_ss_MomSys[kt][BOI_1]->SetBinContent(3,-100);
770 }
771
772 if(TherminatorC2==kFALSE){
773 C2_ss[kt][BOI_1]->DrawCopy();
9b34d0fd 774 C2_ss_MomSys[kt][BOI_1]->SetFillColor(kRed-9); C2_os_MomSys[kt][BOI_1]->SetFillColor(kBlue-9);
9b479142 775 C2_ss_MomSys[kt][BOI_1]->DrawCopy("E2 same");
776 C2_os_MomSys[kt][BOI_1]->DrawCopy("E2 same");
777 C2_ss[kt][BOI_1]->DrawCopy("same");
778 C2_os[kt][BOI_1]->DrawCopy("same");
9b34d0fd 779 hfitC2ss_noG[kt][BOI_1]->SetLineStyle(3); fitC2ss_noG[kt][BOI_1]->SetLineStyle(3);
780 hfitC2os_noG[kt][BOI_1]->SetLineStyle(3); fitC2os_noG[kt][BOI_1]->SetLineStyle(3);
781 hfitC2ss_noG[kt][BOI_1]->DrawCopy("C same");
782 hfitC2os_noG[kt][BOI_1]->DrawCopy("C same");
783 hfitC2ss_noGEW[kt][BOI_1]->DrawCopy("C same");
784 hfitC2os_noGEW[kt][BOI_1]->DrawCopy("C same");
b4f433d0 785 //fitC2ss_yesG[kt][BOI_1]->DrawCopy("same");
786 //fitC2os_yesG[kt][BOI_1]->DrawCopy("same");
787
9b479142 788 }else{// Therminator
789 C2Therm_ss[kt][BOI_1]->DrawCopy();
790 C2Therm_os[kt][BOI_1]->DrawCopy("same");
791 fitC2ss_Therm[kt][BOI_1]->DrawCopy("same");
792 fitC2os_Therm[kt][BOI_1]->DrawCopy("same");
793 fitC2ss_ThermGaus[kt][BOI_1]->SetLineStyle(3);
794 fitC2os_ThermGaus[kt][BOI_1]->SetLineStyle(3);
795 fitC2ss_ThermGaus[kt][BOI_1]->DrawCopy("same");
796 fitC2os_ThermGaus[kt][BOI_1]->DrawCopy("same");
797 }
798 //
799 if(kt==0){
800 TLatex *Specif;
3e748846 801 if(!TherminatorC2) Specif = new TLatex(0.007,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.91,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");// ALICE specifications
9b34d0fd 802 else Specif = new TLatex(0.0015,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.91,"Therminator");// Therminator
803 Specif->SetTextFont(TextFont);
804 Specif->SetTextSize(SizeSpecif*.9);
9b479142 805 TString *CentName;
806 if(!TherminatorC2) {CentName = new TString(); *CentName += BOI_1*5; CentName->Append("-"); *CentName += (BOI_1+1)*5; CentName->Append("%");}
807 else {
808 if(BOI_1==0) CentName = new TString("0-5%");
809 if(BOI_1==9) CentName = new TString("40-50%");
810 }
9b34d0fd 811 TLatex *Specif2 = new TLatex(0.027,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.8,CentName->Data());
812 Specif2->SetTextFont(TextFont);
3e748846 813 Specif2->SetTextSize(SizeSpecif);
9b479142 814 Specif->Draw("same");
815 Specif2->Draw("same");
816 }
817 if(kt==1){
818 legend1->AddEntry(C2_ss[kt][BOI_1],"same-charge","p");
819 legend1->AddEntry(C2_os[kt][BOI_1],"mixed-charge","p");
820 legend1->Draw("same");
821 }
822 if(kt==2){
823 fitC2ss_noG[kt][BOI_1]->SetLineColor(1);
824 fitC2ss_yesG[kt][BOI_1]->SetLineColor(1);
b4f433d0 825 fitC2ss_noGEW[kt][BOI_1]->SetLineColor(1);
9b479142 826 fitC2ss_Therm[kt][BOI_1]->SetLineColor(1);
827 fitC2ss_ThermGaus[kt][BOI_1]->SetLineColor(1);
828 if(!TherminatorC2) {
b4f433d0 829 legend2->AddEntry(fitC2ss_noG[kt][BOI_1],"Gauss","l");
830 legend2->AddEntry(fitC2ss_noGEW[kt][BOI_1],"Edgeworth","l");
831 //legend2->AddEntry(fitC2ss_yesG[kt][BOI_1],"Gauss, G#neq0","l");
9b479142 832 }else {
b4f433d0 833 legend2->AddEntry(fitC2ss_ThermGaus[kt][BOI_1],"Gauss","l");// Therminator
834 legend2->AddEntry(fitC2ss_Therm[kt][BOI_1],"Edgeworth","l");// Therminator
9b479142 835 }
836 legend2->Draw("same");
837 }
838 TLatex *KtLabel;
3e748846 839 if(kt==0) KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.2<#font[12]{k}_{T}<0.3 GeV/#font[12]{c}");// 0.003,.988
840 else if(kt==1) KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.3<#font[12]{k}_{T}<0.4 GeV/#font[12]{c}");
841 else if(kt==2) KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.4<#font[12]{k}_{T}<0.5 GeV/#font[12]{c}");
842 else if(kt==3) KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.5<#font[12]{k}_{T}<0.6 GeV/#font[12]{c}");
843 else if(kt==4) KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.6<#font[12]{k}_{T}<0.7 GeV/#font[12]{c}");
844 else KtLabel = new TLatex(0.028,MIN[BOI_1]+(MAX[BOI_1]-MIN[BOI_1])*0.68,"0.7<#font[12]{k}_{T}<0.8 GeV/#font[12]{c}");
9b34d0fd 845 KtLabel->SetTextFont(TextFont);// 23
3e748846 846 KtLabel->SetTextSize(SizeSpecif*SF_correction);// 24
9b479142 847 KtLabel->Draw("same");
848 }
849
850 canC2->cd();
9b34d0fd 851
9b479142 852 TPad *pad1_2 = new TPad("pad1_2","pad1_2",0.0,0.0,1.,1.);
853 pad1_2->SetFillStyle(0);
854 pad1_2->Draw();
855 pad1_2->cd();
3e748846 856 TLatex *C2Label = new TLatex(.04,.94,"#font[12]{C}_{2}");// 0.04,.91
9b34d0fd 857 C2Label->SetTextFont(TextFont);
858 C2Label->SetTextSize(SizeTitle*SF2);
9b479142 859 C2Label->SetTextAngle(90);
860 C2Label->Draw();
3e748846 861 TLatex *Q2Label = new TLatex(.75,.015,"#font[12]{q} (GeV/#font[12]{c})");// .7,.02
9b34d0fd 862 Q2Label->SetTextFont(TextFont);
863 Q2Label->SetTextSize(SizeTitle*SF2);
9b479142 864 Q2Label->Draw();
3e748846 865 TBox *CoverUp = new TBox(0.554,0.05,0.58,.083);
866 CoverUp->SetFillColor(0);
9b479142 867 CoverUp->Draw();
b4f433d0 868
9b479142 869
870
871
872 ////////////////////////////////////////////////////
873 ////////////////////////////////////////////////////
874 // C2 fit parameters
b4f433d0 875 TCanvas *can2 = new TCanvas("can2", "can2",680,0,600,900);// 680,0,600,900 with G
9b479142 876 can2->SetHighLightColor(2);
9b479142 877 gStyle->SetOptFit(0111);
3e748846 878 can2->SetFillColor(0);//10
9b479142 879 can2->SetBorderMode(0);
880 can2->SetBorderSize(2);
881 can2->SetFrameFillColor(0);
882 can2->SetFrameBorderMode(0);
883 can2->SetFrameBorderMode(0);
884 can2->cd();
885
886 TPad *pad2 = new TPad("pad2","pad2",0.,0.,1.,1.);
887 gPad->SetGridx(0);
888 gPad->SetGridy(0);
889 gPad->SetTickx();
890 gPad->SetTicky();
891 pad2->SetTopMargin(0.02);
892 pad2->SetLeftMargin(0.0);
893 pad2->SetRightMargin(0.0);
9b34d0fd 894 pad2->SetBottomMargin(0.07);
b4f433d0 895 pad2->Divide(1,2,0,0);// 1,3,0,0 with G
9b479142 896 pad2->Draw();
897 pad2->cd();
3e748846 898 TLegend *legend3 = new TLegend(.44,.66, .96,.98,NULL,"brNDC");// .56,.6, .96,.98
9b34d0fd 899 legend3->SetBorderSize(0);
900 legend3->SetTextFont(TextFont);
901 legend3->SetTextSize(SizeLegend*SF1);
9b479142 902 legend3->SetFillColor(0);
903
3e748846 904 double LowerLimits[4]={0.46, 0.32, 4.8, 0.08};// 0.46, 0.32, 4.8, 0.08
905 double UpperLimits[4]={0.87, 0.63, 13.2, 1.28};// 0.84, 0.63, 13.2, 1.28
9b479142 906 for(int par=1; par<=3; par++){
b4f433d0 907 //pad2->cd(par);
908 if(par!=1 && par!=3) continue;
909 if(par==1) pad2->cd(1);
910 else pad2->cd(2);
3e748846 911 gPad->SetRightMargin(0.03); gPad->SetLeftMargin(0.14);
912 if(par==3) gPad->SetBottomMargin(0.16);// 0.22 with G
b4f433d0 913 // settings applied to ParHisto_coh[par-1][BOI_1] to include G.
914 ParHisto_ch[par-1][BOI_1]->SetMinimum(LowerLimits[par-1]);
915 ParHisto_ch[par-1][BOI_1]->SetMaximum(UpperLimits[par-1]);
9b34d0fd 916 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetTitleFont(TextFont); ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitleFont(TextFont);
917 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetLabelFont(TextFont); ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetLabelFont(TextFont);
9b479142 918 if(par==3){
9b34d0fd 919 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetLabelSize(SizeLabel*SF1);
920 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetTitleSize(SizeTitle*SF1);
921 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetLabelSize(SizeLabel*SF1);
9b479142 922 }else{
b4f433d0 923 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetLabelSize(0);
924 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetTitleSize(0);
925 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetLabelSize(0);
9b479142 926 }
3e748846 927 if(par==1) SF_correction=1.0;
928 else SF_correction=0.97;
929 ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitleSize(SizeTitle*SF1*SF_correction);
930 ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetLabelSize(SizeLabel*SF1*SF_correction);
931 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetTitle("#font[12]{k}_{T} (GeV/#font[12]{c})");
b4f433d0 932 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetRangeUser(0.2,0.79);
933 if(par==1) ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitle("#lambda");
3e748846 934 if(par==2) ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitle("#font[12]{G}");
935 if(par==3) ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitle("#font[12]{R_{ch}} (fm)");
936 if(par==4) ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitle("#font[12]{R_{coh}} (fm)");
937 ParHisto_ch[par-1][BOI_1]->GetXaxis()->SetTitleOffset(0.9);// 2.0
938 ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetTitleOffset(0.8);// 1.4
939 ParHisto_ch[par-1][BOI_1]->GetYaxis()->SetNdivisions(505);
9b479142 940
941 // ParN=2 (lambda), 3(G), 4(Rch), 5(Rcoh)
b4f433d0 942 ParHisto_ch[par-1][BOI_1]->Draw();
b4f433d0 943
944
9b34d0fd 945 TH1D *ParHisto1_Sys=(TH1D*)ParHisto_ch[par-1][BOI_1]->Clone();
946 TH1D *ParHisto2_Sys=(TH1D*)ParHisto_ch[par-1][BOI_2]->Clone();
947 TH1D *ParHisto3_Sys=(TH1D*)ParHisto_chEW[par-1][BOI_1]->Clone();
948 TH1D *ParHisto4_Sys=(TH1D*)ParHisto_chEW[par-1][BOI_2]->Clone();
949 TF1 *RSys1=new TF1("RSys1","0.3",0,1);// M0, R Gauss Fit range syst in fm.
950 TF1 *RSys2=new TF1("RSys2","0.2",0,1);// M9, R Gauss Fit range syst in fm.
951 TF1 *RSys3=new TF1("RSys3","0.2",0,1);// M0, R EW Fit range syst in fm.
952 TF1 *RSys4=new TF1("RSys4","0.2",0,1);// M9, R EW Fit range syst in fm.
953 TF1 *LamSys1=new TF1("LamSys1","0.015",0,1);// M0, Lam Gauss Fit range syst in fm.
954 TF1 *LamSys2=new TF1("LamSys2","0.015",0,1);// M9, Lam Gauss Fit range syst in fm.
955 TF1 *LamSys3=new TF1("LamSys3","0.01",0,1);// M0, Lam EW Fit range syst in fm.
956 TF1 *LamSys4=new TF1("LamSys4","0.01",0,1);// M9, Lam EW Fit range syst in fm.
957
958 for(int bin=1; bin<=ParHisto1_Sys->GetNbinsX(); bin++){
959 double ktvalue=ParHisto1_Sys->GetXaxis()->GetBinCenter(bin);
960 if(par==1){// lambda Sys
961 ParHisto1_Sys->SetBinError(bin, LamSys1->Eval(ktvalue));
962 ParHisto2_Sys->SetBinError(bin, LamSys2->Eval(ktvalue));
963 ParHisto3_Sys->SetBinError(bin, LamSys3->Eval(ktvalue));
964 ParHisto4_Sys->SetBinError(bin, LamSys4->Eval(ktvalue));
965 }
966 if(par==3){// R sys
967 ParHisto1_Sys->SetBinError(bin, RSys1->Eval(ktvalue));
968 ParHisto2_Sys->SetBinError(bin, RSys2->Eval(ktvalue));
969 ParHisto3_Sys->SetBinError(bin, RSys3->Eval(ktvalue));
970 ParHisto4_Sys->SetBinError(bin, RSys4->Eval(ktvalue));
971 }
b4f433d0 972 }
9b34d0fd 973 ParHisto1_Sys->SetMarkerSize(0); ParHisto1_Sys->SetFillStyle(1000); ParHisto1_Sys->SetFillColor(kGray);
974 ParHisto2_Sys->SetMarkerSize(0); ParHisto2_Sys->SetFillStyle(1000); ParHisto2_Sys->SetFillColor(kGray);
975 ParHisto3_Sys->SetMarkerSize(0); ParHisto3_Sys->SetFillStyle(1000); ParHisto3_Sys->SetFillColor(kBlue-10);
976 ParHisto4_Sys->SetMarkerSize(0); ParHisto4_Sys->SetFillStyle(1000); ParHisto4_Sys->SetFillColor(kBlue-10);
977
978 ParHisto1_Sys->Draw("E2 same"); ParHisto2_Sys->Draw("E2 same"); ParHisto3_Sys->Draw("E2 same"); ParHisto4_Sys->Draw("E2 same");
979 ParHisto_ch[par-1][BOI_1]->Draw("same");
980 ParHisto_ch[par-1][BOI_2]->Draw("same");
981 ParHisto_chEW[par-1][BOI_1]->Draw("same");
982 ParHisto_chEW[par-1][BOI_2]->Draw("same");
983
9b479142 984 if(par==1){
b4f433d0 985 //legend3->AddEntry(ParHisto_chEWfromTherm[par-1][BOI_1],"E_{w}, 0-5%","p"); legend3->AddEntry(ParHisto_chEWfromTherm[par-1][BOI_2],"E_{w}, 45-50%","p");
986 //legend3->AddEntry(ParHisto_coh[par-1][BOI_1],"G#neq0, 0-5%","p"); legend3->AddEntry(ParHisto_coh[par-1][BOI_2],"G#neq0, 45-50%","p");
987 legend3->AddEntry(ParHisto_ch[par-1][BOI_1],"Gauss, 0-5%","p"); legend3->AddEntry(ParHisto_ch[par-1][BOI_2],"Gauss, 45-50%","p");
3e748846 988 legend3->AddEntry(ParHisto_chEW[par-1][BOI_1],"Edgeworth, 0-5%","p"); legend3->AddEntry(ParHisto_chEW[par-1][BOI_2],"Edgeworth, 45-50%","p");
9b479142 989 }
990
b4f433d0 991 if(par==3) legend3->Draw("same");// par==2 with G
9b479142 992 if(par==1){
3e748846 993 TLatex *Specif_2 = new TLatex(0.35,LowerLimits[par-1]+(UpperLimits[par-1]-LowerLimits[par-1])*.94,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");
9b34d0fd 994 Specif_2->SetTextFont(TextFont);
995 Specif_2->SetTextSize(SizeSpecif*SF1);
996 Specif_2->Draw();
9b479142 997 }
998 }
999
1000 ////////////////////////////////////////////////////
1001 ////////////////////////////////////////////////////
1002 // C2(+-) Therminator comparisons
3e748846 1003
9b34d0fd 1004 TCanvas *can3 = new TCanvas("can3", "can3",1800,0,600,600);// 11,53,700,500
9b479142 1005 can3->SetHighLightColor(2);
1006 gStyle->SetOptFit(0111);
3e748846 1007 can3->SetFillColor(0);//10
9b479142 1008 can3->SetBorderMode(0);
1009 can3->SetBorderSize(2);
1010 can3->SetFrameFillColor(0);
1011 can3->SetFrameBorderMode(0);
1012 can3->SetFrameBorderMode(0);
1013 can3->cd();
1014 TPad *pad3 = new TPad("pad3","pad3",0.0,0.0,1.,1.);
1015 gPad->SetGridx(0);
1016 gPad->SetGridy(0);
1017 gPad->SetTickx();
1018 gPad->SetTicky();
1019 pad3->SetTopMargin(0.02);//0.05
1020 pad3->SetRightMargin(0.01);//3e-2
9b34d0fd 1021 pad3->SetBottomMargin(0.07);//0.12
1022 //pad3->Divide(1,2,0,0);
9b479142 1023 pad3->Draw();
1024 pad3->cd();
3e748846 1025 TLegend *legend4 = new TLegend(.24,.78, .85,.97,NULL,"brNDC");//.45 or .4 for x1
9b34d0fd 1026 legend4->SetBorderSize(0);
9b479142 1027 legend4->SetFillColor(0);
9b34d0fd 1028
3e748846 1029 gPad->SetRightMargin(0.01); gPad->SetLeftMargin(0.16); gPad->SetBottomMargin(0.11); gPad->SetTopMargin(0.01);
9b34d0fd 1030
9b479142 1031 TH1D *C2_os_ktcompare[2][KTBINS];
1032 TH1D *C2Therm_os_ktcompare[2][KTBINS];
9b34d0fd 1033 float MIN_ktcomp=0.985, MAX_ktcomp=1.045;// MIN_ktcomp=0.92, MAX_ktcomp=2.04 (cumulant comparison)
9b479142 1034 for(int type=0; type<2; type++){
9b34d0fd 1035 //pad3->cd(type+1);
1036 //gPad->SetRightMargin(RightMargin); gPad->SetLeftMargin(0.12);
1037 //gPad->SetBottomMargin(0.145);
9b479142 1038 for(int kt=0; kt<KTBINS; kt++){
1039 C2_os_ktcompare[type][kt] = (TH1D*)C2_os[kt][BOI_1]->Clone();
1040 TH1D *FSIbase;
1041 if(type==0) {
1042 FSIbase = (TH1D*)C2Therm_os[0][BOI_1]->Clone();// diluted
1043 C2Therm_os_ktcompare[type][kt] = (TH1D*)C2Therm_os[kt][BOI_1]->Clone();// diluted
1044 }else{
1045 FSIbase = (TH1D*)K2_os[0][BOI_1]->Clone();// undiluted
1046 C2Therm_os_ktcompare[type][kt] = (TH1D*)K2_os[kt][BOI_1]->Clone();// undiluted
1047 }
1048 //
9b34d0fd 1049
9b479142 1050 for(int qbin=2; qbin<=20; qbin++){
1051 if(C2_os[0][BOI_1]->GetBinContent(qbin)<=0) continue;
1052 if(C2Therm_os[0][BOI_1]->GetBinContent(qbin)<=0) continue;
9b34d0fd 1053
1054 double value = C2_os_ktcompare[type][kt]->GetBinContent(qbin)/C2_os[0][BOI_1]->GetBinContent(qbin);
1055 double value_e = pow(C2_os_ktcompare[type][kt]->GetBinError(qbin)/C2_os[0][BOI_1]->GetBinContent(qbin),2);
1056 value_e += pow(C2_os[0][BOI_1]->GetBinError(qbin) * C2_os_ktcompare[type][kt]->GetBinContent(qbin)/pow(C2_os[0][BOI_1]->GetBinContent(qbin),2),2);
1057 value_e = sqrt(value_e);
1058 double valueC2Therm = C2Therm_os_ktcompare[type][kt]->GetBinContent(qbin)/FSIbase->GetBinContent(qbin);
1059 double valueC2Therm_e = pow(C2Therm_os_ktcompare[type][kt]->GetBinError(qbin)/FSIbase->GetBinContent(qbin),2);
1060 valueC2Therm_e += pow(FSIbase->GetBinError(qbin) * C2Therm_os_ktcompare[type][kt]->GetBinContent(qbin)/pow(FSIbase->GetBinContent(qbin),2),2);
1061 valueC2Therm_e = sqrt(valueC2Therm_e);
1062 // cumulant comparison
3e748846 1063 //double value = (C2_os_ktcompare[type][kt]->GetBinContent(qbin)-1.)/(C2_os[0][BOI_1]->GetBinContent(qbin)-1.);
1064 //double value_e = pow(C2_os_ktcompare[type][kt]->GetBinError(qbin)/(C2_os[0][BOI_1]->GetBinContent(qbin)-1.),2);
1065 //value_e += pow(C2_os[0][BOI_1]->GetBinError(qbin) * (C2_os_ktcompare[type][kt]->GetBinContent(qbin)-1.)/pow((C2_os[0][BOI_1]->GetBinContent(qbin)-1.),2),2);
1066 //value_e = sqrt(value_e);
1067 //double valueC2Therm = (C2Therm_os_ktcompare[type][kt]->GetBinContent(qbin)-1.)/(FSIbase->GetBinContent(qbin)-1.);
1068 //double valueC2Therm_e = pow(C2Therm_os_ktcompare[type][kt]->GetBinError(qbin)/(FSIbase->GetBinContent(qbin)-1.),2);
1069 //valueC2Therm_e += pow(FSIbase->GetBinError(qbin) * (C2Therm_os_ktcompare[type][kt]->GetBinContent(qbin)-1.)/pow((FSIbase->GetBinContent(qbin)-1.),2),2);
1070 //valueC2Therm_e = sqrt(valueC2Therm_e);
1071
9b479142 1072
1073 //
1074 C2_os_ktcompare[type][kt]->SetBinContent(qbin, value);
1075 C2_os_ktcompare[type][kt]->SetBinError(qbin, value_e);
1076 C2Therm_os_ktcompare[type][kt]->SetBinContent(qbin, valueC2Therm);
b4f433d0 1077 if(type==0) C2Therm_os_ktcompare[type][kt]->SetBinError(qbin, valueC2Therm_e);
1078 else C2Therm_os_ktcompare[type][kt]->SetBinError(qbin, C2Therm_os_ktcompare[0][kt]->GetBinError(qbin));
9b479142 1079 if(kt==0) {C2_os_ktcompare[type][kt]->SetBinError(qbin, 0);}
1080 }
1081 C2_os_ktcompare[type][kt]->SetBinContent(1,-10); C2Therm_os_ktcompare[type][kt]->SetBinContent(1,-10);
9b34d0fd 1082 C2_os_ktcompare[type][kt]->SetMinimum(MIN_ktcomp);
1083 C2_os_ktcompare[type][kt]->SetMaximum(MAX_ktcomp);
9b479142 1084 C2_os_ktcompare[type][kt]->GetXaxis()->SetRangeUser(0,0.09);
1085 int MarkerStyle=20;
1086 if(kt==1) MarkerStyle=24;
1087 if(kt==5) MarkerStyle=21;
1088 C2_os_ktcompare[type][kt]->SetMarkerStyle(MarkerStyle);
b4f433d0 1089 C2Therm_os_ktcompare[type][kt]->SetMarkerStyle(24);// MarkerStyle
1090 C2Therm_os_ktcompare[type][kt]->SetMarkerColor(1); C2Therm_os_ktcompare[type][kt]->SetLineColor(1);
1091 C2Therm_os_ktcompare[type][kt]->SetMarkerSize(C2_os_ktcompare[type][kt]->GetMarkerSize());
9b479142 1092 C2_os_ktcompare[type][kt]->GetXaxis()->SetTitleOffset(20); C2_os_ktcompare[type][kt]->GetYaxis()->SetTitleOffset(20);
1093
9b34d0fd 1094 C2_os_ktcompare[type][kt]->GetXaxis()->SetNdivisions(606);
1095
b4f433d0 1096 //if(kt!=5 && kt!=1 ) continue;// old way
1097 if(kt!=5) continue;
1098 if(kt==5) {
9b34d0fd 1099 C2_os_ktcompare[type][kt]->GetXaxis()->SetTitleFont(TextFont); C2_os_ktcompare[type][kt]->GetYaxis()->SetTitleFont(TextFont);
1100 C2_os_ktcompare[type][kt]->GetXaxis()->SetTitleSize(SizeTitle*SF2);
1101 C2_os_ktcompare[type][kt]->GetYaxis()->SetTitleSize(SizeTitle*SF2);
3e748846 1102 C2_os_ktcompare[type][kt]->GetXaxis()->SetTitle("#font[12]{q} (GeV/#font[12]{c})");
1103 C2_os_ktcompare[type][kt]->GetXaxis()->SetTitleOffset(0.8);
9b34d0fd 1104 C2_os_ktcompare[type][kt]->GetXaxis()->SetLabelFont(TextFont); C2_os_ktcompare[type][kt]->GetYaxis()->SetLabelFont(TextFont);
1105 C2_os_ktcompare[type][kt]->GetXaxis()->SetLabelSize(SizeLabel*SF2);
1106 C2_os_ktcompare[type][kt]->GetYaxis()->SetLabelSize(SizeLabel*SF2);
1107 if(type==1) C2Therm_os_ktcompare[type][kt]->SetMarkerStyle(25);
1108
1109 if(type==0) C2_os_ktcompare[type][kt]->DrawCopy();
b4f433d0 1110 C2Therm_os_ktcompare[type][kt]->DrawCopy("same");
1111 }
1112 //else {
1113 //C2_os_ktcompare[type][kt]->DrawCopy("same");
1114 //}
3e748846 1115 TString *ktname = new TString("#font[12]{k}_{T,");
b4f433d0 1116 *ktname += kt+1;
3e748846 1117 ktname->Append("}/#font[12]{k}_{T,1}");
b4f433d0 1118 TString *nameReal=new TString(ktname->Data());
1119 nameReal->Append(", ALICE");
1120 TString *nameTherm=new TString(ktname->Data());
1121 nameTherm->Append(", Therminator");
1122 if(type==0) {
9b34d0fd 1123 nameTherm->Append(" (diluted)");
b4f433d0 1124 legend4->AddEntry(C2_os_ktcompare[type][kt],nameReal->Data(),"p");
1125 legend4->AddEntry(C2Therm_os_ktcompare[type][kt],nameTherm->Data(),"p");
9b34d0fd 1126 }else {
1127 TString *nameTherm2=new TString(nameTherm->Data());
3e748846 1128 //nameTherm2->Append(", r*<80 fm");
9b34d0fd 1129 nameTherm2->Append(" (undiluted)");
1130 legend4->AddEntry(C2Therm_os_ktcompare[type][kt],nameTherm2->Data(),"p");
9b479142 1131 }
9b34d0fd 1132
9b479142 1133 }
1134
9b34d0fd 1135 if(type==1) {
3e748846 1136 TLatex *Specif_3 = new TLatex(0.4,0.15,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");
9b34d0fd 1137 Specif_3->SetNDC();
1138 Specif_3->SetTextFont(TextFont);
1139 Specif_3->SetTextSize(SizeSpecif*SF2);
3e748846 1140 TLatex *Specif2_3 = new TLatex(0.55,0.22,"0-5%");
9b34d0fd 1141 Specif2_3->SetNDC();
1142 Specif2_3->SetTextFont(TextFont);
1143 Specif2_3->SetTextSize(SizeLegend*SF2);
1144 Specif_3->Draw();
1145 Specif2_3->Draw();
1146 legend4->SetTextFont(TextFont);
1147 legend4->SetTextSize(SizeLegend*SF2);
1148 legend4->Draw("same");
9b479142 1149 }
9b34d0fd 1150
9b479142 1151 TF1 *Unity = new TF1("Unity","1",0,100);
1152 Unity->SetLineStyle(2);
1153 Unity->Draw("same");
1154
1155 }
1156 can3->cd();
1157 TPad *pad3_2 = new TPad("pad3_2","pad3_2",0.0,0.0,1.,1.);
1158 pad3_2->SetFillStyle(0);
1159 pad3_2->Draw();
1160 pad3_2->cd();
3e748846 1161 TLatex *DthLabel = new TLatex(.044,.58,"#font[12]{C}_{2}^{+-}(#font[12]{k}_{T,6}) / #font[12]{C}_{2}^{+-}(#font[12]{k}_{T,1})");// .04,.96, "D^{+-}_{Th}"
9b34d0fd 1162 DthLabel->SetTextFont(TextFont);
1163 DthLabel->SetTextSize(SizeTitle*SF2);
9b479142 1164 DthLabel->SetTextAngle(90);
1165 DthLabel->Draw();
3e748846 1166
9b479142 1167
1168 ////////////////////////////////////////////////////
1169 ////////////////////////////////////////////////////
1170 // Three-pion correlations
1171
1172
1173 TCanvas *can4 = new TCanvas("can4", "can4",10,300,600,900);// 11,53,700,500
1174 can4->SetHighLightColor(2);
1175 gStyle->SetOptFit(0111);
3e748846 1176 can4->SetFillColor(0);//10
9b479142 1177 can4->SetBorderMode(0);
1178 can4->SetBorderSize(2);
1179 can4->SetFrameFillColor(0);
1180 can4->SetFrameBorderMode(0);
1181 can4->SetFrameBorderMode(0);
1182
1183
3e748846 1184 TPad *pad4 = new TPad("pad4","pad4",0.06,0.06,1.,1.);
9b479142 1185 gPad->SetGridx(0);
1186 gPad->SetGridy(0);
1187 gPad->SetTickx();
1188 gPad->SetTicky();
1189 pad4->SetTopMargin(0.02);//0.05
1190 pad4->SetLeftMargin(0.12);
1191 pad4->SetRightMargin(0.01);//3e-2
3e748846 1192 pad4->SetBottomMargin(0.08);//0.12
9b479142 1193 pad4->Divide(2,3,0,0);
1194 pad4->Draw();
1195
1196
b4f433d0 1197 TLegend *legend5 = new TLegend(.55,.6,.99,.99,NULL,"brNDC");//.55,.6,.99,.95
9b34d0fd 1198 legend5->SetBorderSize(0);
1199 legend5->SetTextFont(TextFont);
1200 legend5->SetTextSize(SizeLegend);
9b479142 1201 legend5->SetFillColor(0);
1202
1203
1204 TF1 *BaseLine = new TF1("BaseLine","1.0",0,1);
1205 BaseLine->SetLineStyle(3);
1206 BaseLine->SetLineColor(1);
1207
1208
1209 TF1 *GaussFit[6];
1210 for(int ii=0; ii<6; ii++) {
1211 TString *nameFit=new TString("GaussFit_");
1212 *nameFit += ii;
b4f433d0 1213 GaussFit[ii] = new TF1(nameFit->Data(),"[0]+[1]*exp(-pow([2]*x,2))",0,0.2);
1214 GaussFit[ii]->SetParameter(0, 1.);
1215 GaussFit[ii]->SetParameter(1, 2.);
1216 GaussFit[ii]->SetParameter(2, 10./FmToGeV);
9b479142 1217 }
1218
b4f433d0 1219 TF1 *C3_SysFit=new TF1("C3_SysFit","sqrt(pow(pol3,2)+pow(0.003,2))",0,0.15);// MRC + PID
9b479142 1220 C3_SysFit->FixParameter(0, -0.01354);
1221 C3_SysFit->FixParameter(1, 0.2948);
1222 C3_SysFit->FixParameter(2, -2.237);
1223 C3_SysFit->FixParameter(3, 5.697);
1224 TH1D *C3Sys[6];// cb
1225 TH1D *c3Sys[6];// cb
b4f433d0 1226 TF1 *c3_Coul_SysFit=new TF1("c3_Coul_SysFit","pol1",0,0.15);
9b479142 1227 TH1D *c3Sys_Coul[6];
1228 float C3MAX=2.3;
1229 if(ChProdBOI==1) C3MAX=2.1;
1230
1231 for(int cb=0; cb<6; cb++) {
1232 pad4->cd(cb+1);
1233 gPad->SetRightMargin(RightMargin);
3e748846 1234 if(cb%2==0) gPad->SetRightMargin(0.001);
1235 if(cb==4) SF_correction=0.92;
1236 else if(cb==5) SF_correction=1.015;
1237 else SF_correction=1.0;
9b34d0fd 1238 //if(cb==0) Dim1=gPad->GetAbsWNDC()*gPad->GetAbsHNDC()*gPad->PixeltoX(10)*gPad->PixeltoY(10);
1239 //Dim2=gPad->GetAbsWNDC()*gPad->GetAbsHNDC()*gPad->PixeltoX(10)*gPad->PixeltoY(10);
9b479142 1240
1241 C3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerSize(1.5); C3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerSize(1.5); c3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerSize(1.5); c3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerSize(1.5);
1242 C3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerStyle(20); C3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerStyle(21);
1243 c3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerStyle(24); c3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerStyle(25);
1244 C3EW[CoulChoice][KT3Bin][cb]->SetMarkerStyle(22);
1245 C3EW[CoulChoice][KT3Bin][cb]->SetMarkerColor(1);
1246 C3EW[CoulChoice][KT3Bin][cb]->SetLineColor(1);
1247 C3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerColor(2); C3merged[0][CoulChoice][KT3Bin][cb]->SetLineColor(2); c3merged[0][CoulChoice][KT3Bin][cb]->SetMarkerColor(2); c3merged[0][CoulChoice][KT3Bin][cb]->SetLineColor(2);
1248 C3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerColor(4); C3merged[1][CoulChoice][KT3Bin][cb]->SetLineColor(4); c3merged[1][CoulChoice][KT3Bin][cb]->SetMarkerColor(4); c3merged[1][CoulChoice][KT3Bin][cb]->SetLineColor(4);
1249 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetTitle(""); C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetYaxis()->SetTitle("");
9b34d0fd 1250 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelFont(TextFont);
1251 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelFont(TextFont);
3e748846 1252 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelSize(SizeLabel*SF_correction);
1253 if(cb%2==0) C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelSize(SizeLabel*SF_correction);
9b479142 1254 else {C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelSize(.0);}
1255 if(cb<4) C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelSize(.0);
9b34d0fd 1256 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetRangeUser(0,0.11);
9b479142 1257 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->SetMinimum(0.9);//0.9
1258 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->SetMaximum(C3MAX);
1259 if(ChProdBOI==0) c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->Fit(GaussFit[cb],"IMEQN","",0.02,0.1);
1260 //
1261 if(cb>2 && ChProdBOI==0) {
1262 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->SetBinContent(2,-100);
1263 c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->SetBinContent(2,-100);
1264 }
9b34d0fd 1265 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetXaxis()->SetNdivisions(404);
3e748846 1266 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetYaxis()->SetNdivisions(505);
9b34d0fd 1267 //TGaxis::SetMaxDigits(3);
9b479142 1268 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->DrawCopy();
1269 c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->DrawCopy("same");
1270
1271 C3Sys[cb] = (TH1D*)C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->Clone();
1272 c3Sys[cb] = (TH1D*)c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->Clone();
1273 c3Sys_Coul[cb] = (TH1D*)c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->Clone();// Coulomb choice
1274 c3Sys_Coul[cb]->Add(c3merged[ChProdBOI][1-int(CoulChoice)][KT3Bin][cb],-1);// Alternate Coulomb
1275 c3Sys_Coul[cb]->Fit(c3_Coul_SysFit,"MENQ","",0,0.1);
b4f433d0 1276 for(int ii=1; ii<=11; ii++){
9b479142 1277 double Q3=(ii-0.5)*0.01;
b4f433d0 1278 C3Sys[cb]->SetBinError(ii, C3_SysFit->Eval(Q3) * C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetBinContent(ii));// MRC + PID
9b479142 1279 //
b4f433d0 1280 double err = pow(0.2 * C3_SysFit->Eval(Q3) * c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetBinContent(ii),2);// MRC + PID (1/5 that of C3)
9b479142 1281 err += pow(c3_Coul_SysFit->Eval(Q3),2);// Coulomb
1282 c3Sys[cb]->SetBinError(ii, sqrt(err));
1283 }
1284 C3Sys[cb]->SetMarkerSize(0); c3Sys[cb]->SetMarkerSize(0);
1285 C3Sys[cb]->SetFillStyle(1000); c3Sys[cb]->SetFillStyle(1000);
9b34d0fd 1286 if(ChProdBOI==0) { C3Sys[cb]->SetFillColor(kRed-9); c3Sys[cb]->SetFillColor(kRed-9);}
1287 else {C3Sys[cb]->SetFillColor(kBlue-9); c3Sys[cb]->SetFillColor(kBlue-9);}
9b479142 1288 C3Sys[cb]->Draw("E2 same");
1289 c3Sys[cb]->Draw("E2 same");
1290 C3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->DrawCopy("same");// redraw
1291 c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->DrawCopy("same");// redraw
1292 BaseLine->Draw("same");
1293 //ResidueFit->Draw("same");
1294
1295
1296 if(cb==0){
1297 if(ChProdBOI==0){
3e748846 1298 legend5->AddEntry(C3merged[ChProdBOI][CoulChoice][KT3Bin][cb],"#font[12]{C}_{3}^{#pm#pm#pm}","p");
1299 legend5->AddEntry(c3merged[ChProdBOI][CoulChoice][KT3Bin][cb],"#font[12]{#bf{c}}_{3}^{#pm#pm#pm}","p");
9b479142 1300 }else{
3e748846 1301 legend5->AddEntry(C3merged[ChProdBOI][CoulChoice][KT3Bin][cb],"#font[12]{C}_{3}^{#pm#pm#mp}","p");
1302 legend5->AddEntry(c3merged[ChProdBOI][CoulChoice][KT3Bin][cb],"#font[12]{#bf{c}}_{3}^{#pm#pm#mp}","p");
9b479142 1303 }
1304 legend5->Draw("same");
1305 }
3e748846 1306 TLatex *Specif_3 = new TLatex(0.007,0.9+(C3MAX-0.9)*0.91,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");// ALICE specifications
9b34d0fd 1307 Specif_3->SetTextFont(TextFont);
3e748846 1308 Specif_3->SetTextSize(SizeSpecif*SF_correction);
9b34d0fd 1309 if(cb==3) Specif_3->Draw("same");
9b479142 1310
1311 TString *KTString = new TString("");
3e748846 1312 if(KT3Bin==0) KTString->Append("0.16<#font[12]{K}_{T,3}<0.3 GeV/#font[12]{c}");
1313 else KTString->Append("0.3<#font[12]{K}_{T,3}<1.0 GeV/#font[12]{c}");
9b34d0fd 1314 TLatex *Kt3Name_1 = new TLatex(0.016,0.9+(C3MAX-0.9)*0.91,KTString->Data());
1315 Kt3Name_1->SetTextFont(TextFont);
3e748846 1316 Kt3Name_1->SetTextSize(SizeLabel*SF_correction);
9b34d0fd 1317 if(cb==1) Kt3Name_1->Draw("same");
1318
1319 TLatex *CentLabel_1;
1320 if(cb==0) CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"0-5%");// ALICE specifications
1321 else if(cb==1) CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"5-10%");
1322 else if(cb==2) CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"10-20%");
1323 else if(cb==3) CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"20-30%");
1324 else if(cb==4) CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"30-40%");
1325 else CentLabel_1 = new TLatex(0.07,C3MAX-0.75,"40-50%");
9b479142 1326
9b34d0fd 1327 CentLabel_1->SetTextFont(TextFont);
3e748846 1328 CentLabel_1->SetTextSize(SizeLabel*SF_correction);//.1
9b34d0fd 1329 CentLabel_1->Draw("same");
9b479142 1330
1331 if(ChProdBOI==1){
1332 TF1 *ResidueFit=new TF1("ResidueFit","pol0",0,0.1);
1333 c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->Fit(ResidueFit,"NME","",0,0.1);
1334 double Chi2=0;
1335 for(int q3bin=2; q3bin<=10; q3bin++){
1336 Chi2 += pow(c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetBinContent(q3bin)-1.0,2)/pow(c3merged[ChProdBOI][CoulChoice][KT3Bin][cb]->GetBinError(q3bin),2);
1337 }
1338 cout<<"Chi2/NDF MC from baseline: "<<Chi2/9.<<endl;
1339 }
1340
1341 }
1342
1343 can4->cd();
1344 TPad *pad4_2 = new TPad("pad4_2","pad4_2",0.0,0.0,1.,1.);
1345 pad4_2->SetFillStyle(0);
1346 pad4_2->Draw();
1347 pad4_2->cd();
3e748846 1348 TLatex *C3Label = new TLatex(.04,.86,"#font[12]{C}_{3} or #font[12]{#bf{c}}_{3}");// 0.05,0.9
9b34d0fd 1349 C3Label->SetTextFont(TextFont);
1350 C3Label->SetTextSize(SizeTitle*SF2);
9b479142 1351 C3Label->SetTextAngle(90);
1352 C3Label->Draw();
3e748846 1353 TLatex *Q3Label = new TLatex(.72,.02,"#font[12]{Q}_{3} (GeV/#font[12]{c})");// 0.05,0.9
9b34d0fd 1354 Q3Label->SetTextFont(TextFont);
1355 Q3Label->SetTextSize(SizeTitle*SF2);
9b479142 1356 Q3Label->Draw();
1357 CoverUp->Draw();
9b479142 1358
1359 if(ChProdBOI==1) return;
1360
1361 /////////////////////////////////////////////////////////////////////////////
1362 /////////////////////////////////////////////////////////////////////////////
1363 // r3
1364
1365 TF1 *ChaoticLimit = new TF1("ChaoticLimit","2.0",0,1);
1366 ChaoticLimit->SetLineStyle(3);
1367
1368 TF1 *QuarticFit[4][6];// Stat/Det/Coul1/Coul2, cb
1369 TF1 *QuadraticFit[4][6];// Stat/Det/Coul1/Coul2, cb
1370 TF1 *QuadQuarticFit[4][6];// Stat/Det/Coul1/Coul2, cb
1371 for(int ii=0; ii<6; ii++){
1372 for(int jj=0; jj<4; jj++){
1373 TString *nameFit=new TString("Quartic_");
1374 *nameFit += ii;
1375 nameFit->Append("_E");
1376 *nameFit += jj;
1377 QuarticFit[jj][ii] = new TF1(nameFit->Data(),"[0]*(1-[1]*pow(x,4))",0,0.2);
1378 QuarticFit[jj][ii]->SetParameter(0,2.0);
1379 QuarticFit[jj][ii]->SetParameter(1,0);
1380 QuarticFit[jj][ii]->SetLineColor(4);
1381 TString *nameFit2=new TString("Quadratic_");
1382 *nameFit2 += ii;
1383 nameFit2->Append("_E");
1384 *nameFit2 += jj;
1385 QuadraticFit[jj][ii] = new TF1(nameFit2->Data(),"[0]*(1-[1]*pow(x,2))",0,0.2);
1386 QuadraticFit[jj][ii]->SetParameter(0,2.0);
1387 QuadraticFit[jj][ii]->SetParameter(1,0);
1388 QuadraticFit[jj][ii]->SetLineColor(1);
1389 QuadraticFit[jj][ii]->SetLineStyle(4);
1390 TString *nameFit3=new TString("QuadQuartic_");
1391 *nameFit3 += ii;
1392 nameFit3->Append("_E");
1393 *nameFit3 += jj;
1394 QuadQuarticFit[jj][ii] = new TF1(nameFit3->Data(),"[0]*(1-[1]*pow(x,2)-[2]*pow(x,4))",0,0.2);
1395 QuadQuarticFit[jj][ii]->SetParameter(0,2.0);
1396 QuadQuarticFit[jj][ii]->SetParameter(1,50);
1397 QuadQuarticFit[jj][ii]->SetParameter(2,5000);
1398 QuadQuarticFit[jj][ii]->SetLineColor(6);
1399 QuadQuarticFit[jj][ii]->SetLineStyle(3);
1400
1401 }
1402 }
1403
1404 TCanvas *can5 = new TCanvas("can5", "can5",680,300,600,900);// 11,53,700,500
1405 can5->SetHighLightColor(2);
1406 gStyle->SetOptFit(0111);
3e748846 1407 can5->SetFillColor(0);//10
9b479142 1408 can5->SetBorderMode(0);
1409 can5->SetBorderSize(2);
1410 can5->SetFrameFillColor(0);
1411 can5->SetFrameBorderMode(0);
1412 can5->SetFrameBorderMode(0);
1413
1414
3e748846 1415 TPad *pad5 = new TPad("pad5","pad5",0.06,0.06,1.,1.);
9b479142 1416 gPad->SetGridx(0);
1417 gPad->SetGridy(0);
1418 gPad->SetTickx();
1419 gPad->SetTicky();
1420 pad5->SetTopMargin(0.02);//0.05
1421 pad5->SetLeftMargin(0.12);//.14 for wide title, .10 for narrow title, 0.08 for DeltaQ
1422 pad5->SetRightMargin(0.01);//3e-2
3e748846 1423 pad5->SetBottomMargin(0.08);//0.12
9b479142 1424 pad5->Divide(2,3,0,0);
1425 pad5->Draw();
1426
3e748846 1427 TLegend *legend6 = new TLegend(.18,.72, .53,.98,NULL,"brNDC");// .25,.16, .6,.46
9b34d0fd 1428 legend6->SetBorderSize(0);
1429 legend6->SetTextFont(TextFont);
1430 legend6->SetTextSize(SizeLegend);
9b479142 1431 legend6->SetFillColor(0);
1432
1433
1434 double LambdaSysPar0[6]={0.0143, 0.014, 0.01294, 0.00204, 0.001227, 0.00999};
1435 double LambdaSysPar1[6]={-0.857, -0.741, -0.6655, -0.3516, -0.2327, -0.3325};
3e748846 1436 double r3MIN = 1.45;// 0.3
1437 double r3MAX = 2.45;// 2.68
1438 r3merged[CoulChoice][1][5]->SetBinContent(3,100); r3merged[CoulChoice][1][5]->SetBinError(3,1000);// remove high stat error point
1439
9b479142 1440 for(int cb=0; cb<6; cb++) {
1441 pad5->cd(cb+1);
1442 gPad->SetRightMargin(RightMargin);
3e748846 1443 if(cb%2==0) gPad->SetRightMargin(0.001);
1444 if(cb==4) SF_correction=0.92;
1445 else if(cb==5) SF_correction=1.015;
1446 else SF_correction=1.0;
9b479142 1447 r3merged[CoulChoice][KT3Bin][cb]->SetMarkerSize(1.5);
1448 r3merged[CoulChoice][KT3Bin][cb]->SetMarkerStyle(20);
1449 r3merged[CoulChoice][KT3Bin][cb]->SetMarkerColor(2); r3merged[CoulChoice][KT3Bin][cb]->SetLineColor(2);
1450 r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetTitle("");
1451 r3merged[CoulChoice][KT3Bin][cb]->GetYaxis()->SetTitle("");
9b34d0fd 1452 r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelFont(TextFont);
1453 r3merged[CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelFont(TextFont);
3e748846 1454 r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelSize(SizeLabel*SF_correction);
1455 if(cb%2==0) r3merged[CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelSize(SizeLabel*SF_correction);
9b479142 1456 else {r3merged[CoulChoice][KT3Bin][cb]->GetYaxis()->SetLabelSize(.0);}
1457 if(cb<4) r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetLabelSize(.0);
3e748846 1458 r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetRangeUser(0,0.071);// was 0,0.11
b4f433d0 1459 r3merged[CoulChoice][KT3Bin][cb]->SetMinimum(r3MIN);
1460 r3merged[CoulChoice][KT3Bin][cb]->SetMaximum(r3MAX);
9b479142 1461 if(cb>2 || KT3Bin==1) r3merged[CoulChoice][KT3Bin][cb]->SetBinContent(2,-100);// 10-20 MeV bin is insignificant
3e748846 1462 r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->SetNdivisions(502);// was 404
1463 r3merged[CoulChoice][KT3Bin][cb]->GetYaxis()->SetNdivisions(505);// was absent
9b479142 1464 r3merged[CoulChoice][KT3Bin][cb]->DrawCopy();
9b34d0fd 1465 //ChaoticLimit->Draw("same");
9b479142 1466
1467 ///////////////
1468 // Systematics
b4f433d0 1469 TF1 *ResidueFit=new TF1("ResidueFit","[0]+[1]*exp(-[2]*x)",0,0.2);
9b479142 1470 ResidueFit->SetParameter(0,1);
1471 ResidueFit->SetParameter(1,0.001);
1472 ResidueFit->SetParameter(2,1);
1473 ResidueFit->SetParLimits(0,.9,1.1); ResidueFit->SetParLimits(1,-.2,.2); ResidueFit->SetParLimits(2,-40,200);
b4f433d0 1474 c3merged[1][CoulChoice][KT3Bin][cb]->Fit(ResidueFit,"IMNQ","",0.01,0.14);// Remove "N" to see residue fits
1475 //cout<<setprecision(5)<<ResidueFit->Eval(0.08)<<endl;
9b479142 1476 TH1D *r3MethSys = (TH1D*)r3merged[CoulChoice][KT3Bin][cb]->Clone();
1477 r3MethSys->SetBinContent(1,-10);
1478 r3MethSys->SetMarkerSize(0);
1479 r3MethSys->SetFillStyle(1000);
1480 r3MethSys->SetFillColor(kRed-10);
1481 TH1D *r3DetSys = (TH1D*)r3merged[CoulChoice][KT3Bin][cb]->Clone();
1482 r3DetSys->SetBinContent(1,-10);
1483 r3DetSys->SetBinContent(1,-10);
1484 r3DetSys->SetMarkerSize(0);
1485 r3DetSys->SetFillStyle(1000);
1486 r3DetSys->SetFillColor(kGray+2);
1487 TH1D *r3MethDetSys = (TH1D*)r3merged[CoulChoice][KT3Bin][cb]->Clone();// for fitting systematics
1488 TH1D *r3Coul1 = (TH1D*)r3merged[CoulChoice][KT3Bin][cb]->Clone();
1489 r3Coul1->SetMarkerSize(0);
1490 TH1D *r3Coul2 = (TH1D*)r3merged[1-int(CoulChoice)][KT3Bin][cb]->Clone();
1491
b4f433d0 1492 TF1 *r3LambdaSysFit = new TF1("r3LambdaSysFit","pol1",0,0.15);
9b479142 1493 r3LambdaSysFit->FixParameter(0,LambdaSysPar0[cb]);
1494 r3LambdaSysFit->FixParameter(1,LambdaSysPar1[cb]);
1495 const int Nb = 20;
1496 double x[Nb]={0};
1497 double y[Nb]={0};
1498 double exL[Nb]={0};
1499 double exH[Nb]={0};
1500 double eyL[Nb]={0};
1501 double eyH[Nb]={0};
1502 for(int ii=0; ii<Nb; ii++){
1503 x[ii] = r3merged[CoulChoice][KT3Bin][cb]->GetXaxis()->GetBinCenter(ii+1);
1504 y[ii] = r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(ii+1);
1505 double ey = r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(ii+1)-r3merged[1-int(CoulChoice)][KT3Bin][cb]->GetBinContent(ii+1);
1506 if(ey>0) eyL[ii] = fabs(ey);
1507 else eyH[ii] = fabs(ey);
1508 }
1509
1510 TGraphAsymmErrors *r3CoulSys = new TGraphAsymmErrors(Nb,x,y,exL,exH,eyL,eyH);
1511 r3CoulSys->SetLineStyle(2);
1512 for(int bin=1; bin<=11; bin++){
1513 double Q3 = (bin-0.5)*0.01;
b4f433d0 1514 double SysMeth=0;
1515 if(C3merged[1][CoulChoice][KT3Bin][cb]->GetBinContent(bin)>0) {
3e748846 1516 //SysMeth = pow(2 * (fabs(ResidueFit->Eval(Q3)-1.0) + 0.002)/(c3merged[0][CoulChoice][KT3Bin][cb]->GetBinContent(bin)-1.0),2);// mixed-charge baseline (new). Includes 0.002 as allowed variation through lambda powers
1517 SysMeth = pow((fabs(ResidueFit->Eval(Q3)-1.0)+0.002)/(c3merged[0][CoulChoice][KT3Bin][cb]->GetBinContent(bin)-1.0) * r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(bin),2);// factional error determined from MC residue fit divided by SC cumulant. Error = Fraction * r3
1518 }
9b479142 1519 SysMeth += pow(r3LambdaSysFit->Eval(Q3,2)*r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(bin),2);// Lambda 0.7 vs 0.6
1520 double SysDet = pow(0.01*r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(bin),2);// MRC
1521 SysDet += pow(0.01*r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(bin),2);// PID
1522
1523 r3MethSys->SetBinError(bin, sqrt(SysMeth));
1524 r3DetSys->SetBinError(bin, sqrt(SysDet));
9b34d0fd 1525 r3MethDetSys->SetBinError(bin, sqrt(pow(r3merged[CoulChoice][KT3Bin][cb]->GetBinError(bin),2) + SysMeth + SysDet));// was SysDet + SysDet by mistake
9b479142 1526 r3Coul1->SetBinError(bin, sqrt(pow(r3merged[CoulChoice][KT3Bin][cb]->GetBinError(bin),2) + SysMeth + SysDet));
1527 r3Coul2->SetBinError(bin, sqrt(pow(r3merged[CoulChoice][KT3Bin][cb]->GetBinError(bin),2) + SysMeth + SysDet));
1528 }
1529
1530 r3MethSys->Draw("E2 same");
1531 r3DetSys->Draw("E2 same");
1532 r3CoulSys->Draw("P");
1533 //
3e748846 1534 float MinQ3Fit=0.01;//0.01
1535 float MaxQ3Fit=0.08;//0.1
9b479142 1536 if(r3merged[CoulChoice][KT3Bin][cb]->GetBinContent(2)<-10) MinQ3Fit=0.02;
1537
9b34d0fd 1538 r3merged[CoulChoice][KT3Bin][cb]->Fit(QuarticFit[0][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Stat only
1539 r3MethDetSys->Fit(QuarticFit[1][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Meth + Det Sys
1540 r3Coul1->Fit(QuarticFit[2][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// GRS
1541 r3Coul2->Fit(QuarticFit[3][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Omega0
9b479142 1542 //
9b34d0fd 1543 r3merged[CoulChoice][KT3Bin][cb]->Fit(QuadraticFit[0][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Stat only
1544 r3MethDetSys->Fit(QuadraticFit[1][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Meth + Det Sys
1545 r3Coul1->Fit(QuadraticFit[2][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// GRS
1546 r3Coul2->Fit(QuadraticFit[3][cb],"IMENQ","",MinQ3Fit,MaxQ3Fit);// Omega0
9b479142 1547 //
3e748846 1548 //if(cb==1){
1549 //for(int i=0; i<10; i++){
1550 // cout<<r3MethDetSys->GetBinError(i+1)<<endl;
1551 //}
1552 //}
9b479142 1553 //
1554 QuarticFit[2][cb]->Draw("same");
1555 QuadraticFit[2][cb]->Draw("same");
1556 r3merged[CoulChoice][KT3Bin][cb]->Draw("same");
1557 //QuadQuarticFit[1][cb]->Draw("same");
1558 //cout<<"Quartic Chi^2/NDF = "<<QuarticFit[cb]->GetChisquare()/QuarticFit[cb]->GetNDF()<<endl;
1559 //cout<<"Quadratic Chi^2/NDF = "<<QuadraticFit[cb]->GetChisquare()/QuadraticFit[cb]->GetNDF()<<endl;
1560
3e748846 1561 TLatex *Specif_4 = new TLatex(0.007,r3MIN+(r3MAX-r3MIN)*0.91,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");// 0.005,2.1
9b34d0fd 1562 Specif_4->SetTextFont(TextFont);
3e748846 1563 Specif_4->SetTextSize(SizeSpecif*SF_correction);
9b34d0fd 1564 if(cb==3) Specif_4->Draw("same");
b4f433d0 1565
9b479142 1566
1567 TString *KTString = new TString("");
3e748846 1568 if(KT3Bin==0) KTString->Append("0.16<#font[12]{K}_{T,3}<0.3 GeV/#font[12]{c}");
1569 else KTString->Append("0.3<#font[12]{K}_{T,3}<1.0 GeV/#font[12]{c}");
1570 TLatex *Kt3Name_2 = new TLatex(0.002,r3MIN+(r3MAX-r3MIN)*0.91,KTString->Data());//
9b34d0fd 1571 Kt3Name_2->SetTextFont(TextFont);
3e748846 1572 Kt3Name_2->SetTextSize(SizeLabel*SF_correction);
9b34d0fd 1573 if(cb==1) Kt3Name_2->Draw("same");
3e748846 1574
9b34d0fd 1575 TLatex *CentLabel_2;
1576 //TString *CentLabel2;
3e748846 1577 if(cb==0) {CentLabel_2 = new TLatex(0.03,1.52,"0-5%"); }//CentLabel2 = new TString("0-5\\%");}// was 0.04,0.4
1578 else if(cb==1) {CentLabel_2 = new TLatex(0.03,1.52,"5-10%");}// CentLabel2 = new TString("5-10\\%");}
1579 else if(cb==2) {CentLabel_2 = new TLatex(0.03,1.52,"10-20%");}// CentLabel2 = new TString("10-20\\%");}
1580 else if(cb==3) {CentLabel_2 = new TLatex(0.03,1.52,"20-30%");}// CentLabel2 = new TString("20-30\\%");}
1581 else if(cb==4) {CentLabel_2 = new TLatex(0.03,1.52,"30-40%");}// CentLabel2 = new TString("30-40\\%");}
1582 else {CentLabel_2 = new TLatex(0.03,1.52,"40-50%");}// CentLabel2 = new TString("40-50\\%");}
9b479142 1583
9b34d0fd 1584 CentLabel_2->SetTextFont(TextFont);
3e748846 1585 CentLabel_2->SetTextSize(SizeLabel*SF_correction);
9b34d0fd 1586 CentLabel_2->Draw("same");
9b479142 1587
1588 if(cb==0){
1589 legend6->AddEntry(QuarticFit[1][cb],"Quartic","l");
1590 legend6->AddEntry(QuadraticFit[1][cb],"Quadratic","l");
1591 legend6->Draw("same");
1592 }
1593
9b34d0fd 1594 ChaoticLimit->Draw("same");
1595 }// cb bin
9b479142 1596 can5->cd();
1597 TPad *pad5_2 = new TPad("pad5_2","pad5_2",0.0,0.0,1.,1.);
1598 pad5_2->SetFillStyle(0);
1599 pad5_2->Draw();
1600 pad5_2->cd();
3e748846 1601 TLatex *r3Label = new TLatex(.035,.96,"#font[12]{r}_{3}");// 0.035,0.88
9b34d0fd 1602 r3Label->SetTextFont(TextFont);
3e748846 1603 r3Label->SetTextSize(SizeTitle*SF2);
9b479142 1604 r3Label->SetTextAngle(90);
1605 r3Label->Draw();
9b479142 1606 CoverUp->Draw();
b4f433d0 1607 Q3Label->Draw();
1608
1609
9b479142 1610 ////////////////////////////
1611 // Print Quartic r3 fit values
1612 double I_quartic_Avg=0, I_quadratic_Avg=0;
1613 double I_quartic_AvgErr1=0, I_quadratic_AvgErr1=0;
1614 double I_quartic_AvgErr2=0, I_quadratic_AvgErr2=0;
1615 double a_quartic_Avg=0, a_quadratic_Avg=0;
1616 double a_quartic_AvgErr1=0, a_quadratic_AvgErr1=0;
1617 double a_quartic_AvgErr2=0, a_quadratic_AvgErr2=0;
3e748846 1618
9b479142 1619 cout<<"Quartic"<<endl;
1620 for(int cb=0; cb<6; cb++) {
1621 TString *CentLabel2;
1622 if(cb==0) {CentLabel2 = new TString("0-5\\%");}
1623 else if(cb==1) {CentLabel2 = new TString("5-10\\%");}
1624 else if(cb==2) {CentLabel2 = new TString("10-20\\%");}
1625 else if(cb==3) {CentLabel2 = new TString("20-30\\%");}
1626 else if(cb==4) {CentLabel2 = new TString("30-40\\%");}
1627 else {CentLabel2 = new TString("40-50\\%");}
1628
1629 double SF=1000.;
3e748846 1630 double I_Sys1=fabs(QuarticFit[1][cb]->GetParameter(0)-QuarticFit[0][cb]->GetParameter(0));// Meth+Det Sys
1631 double I_Sys2=fabs(QuarticFit[2][cb]->GetParameter(0)-QuarticFit[3][cb]->GetParameter(0));// Coulomb Sys
9b479142 1632 double I_SysTotal = sqrt(pow(I_Sys1,2)+pow(I_Sys2,2));
1633 double a_Sys1=fabs(QuarticFit[1][cb]->GetParameter(1)-QuarticFit[0][cb]->GetParameter(1));
1634 double a_Sys2=fabs(QuarticFit[2][cb]->GetParameter(1)-QuarticFit[3][cb]->GetParameter(1));
1635 a_Sys1 /= SF; a_Sys2 /= SF;
1636 double a_SysTotal = sqrt(pow(a_Sys1,2)+pow(a_Sys2,2));
1637 I_quartic_Avg += QuarticFit[1][cb]->GetParameter(0);// [1] has both Stat+Sys errors
1638 I_quartic_AvgErr1 += pow(QuarticFit[0][cb]->GetParError(0),2);// [0] has only Stat errors
3e748846 1639 I_quartic_AvgErr2 += I_SysTotal;
9b479142 1640 a_quartic_Avg += QuarticFit[1][cb]->GetParameter(1)/SF;
1641 a_quartic_AvgErr1 += pow(QuarticFit[0][cb]->GetParError(1)/SF,2);
9b34d0fd 1642 a_quartic_AvgErr2 += sqrt(pow(a_Sys1,2) + pow(a_Sys2,2));
9b479142 1643
1644 cout.precision(3);
b4f433d0 1645 cout<<"$"<<setprecision(2)<<fixed<<CentLabel2->Data()<<"$ & $"<<QuarticFit[1][cb]->GetParameter(0)<<" \\pm "<<QuarticFit[0][cb]->GetParError(0)<<" \\pm "<<I_SysTotal<<"$ & $"<<setprecision(1)<<fixed<<QuarticFit[1][cb]->GetParameter(1)/SF<<" \\pm "<<QuarticFit[0][cb]->GetParError(1)/SF<<" \\pm "<<a_SysTotal<<"$ \\\\ \\hline"<<endl;
9b479142 1646 }
9b34d0fd 1647 I_quartic_Avg /= 6.; I_quartic_AvgErr1 = sqrt(I_quartic_AvgErr1)/6.; I_quartic_AvgErr2 = I_quartic_AvgErr2/6.;
1648 a_quartic_Avg /= 6.; a_quartic_AvgErr1 = sqrt(a_quartic_AvgErr1)/6.; a_quartic_AvgErr2 = a_quartic_AvgErr2/6.;
9b479142 1649 cout<<"$"<<setprecision(2)<<fixed<<"0-50\\%"<<"$ & $"<<I_quartic_Avg<<" \\pm "<<I_quartic_AvgErr1<<" \\pm "<<I_quartic_AvgErr2<<"$ & $"<<setprecision(1)<<fixed<<a_quartic_Avg<<" \\pm "<<a_quartic_AvgErr1<<" \\pm "<<a_quartic_AvgErr2<<"$ \\\\ \\hline"<<endl;
1650
3e748846 1651 /*TCanvas *can5_5 = new TCanvas("can5_5", "can5_5",680,300,600,600);// 11,53,700,500
1652 can5_5->SetHighLightColor(2);
1653 gStyle->SetOptFit(0111);
1654 can5_5->SetFillColor(0);//10
1655 can5_5->SetBorderMode(0);
1656 can5_5->SetBorderSize(2);
1657 can5_5->SetFrameFillColor(0);
1658 can5_5->SetFrameBorderMode(0);
1659 can5_5->SetFrameBorderMode(0);
1660 can5_5->Draw();
1661 double Iold_K1[6]={1.84,1.84,1.82,1.84,1.79,1.73};
1662 double Iold_K2[6]={1.94,1.92,2.02,1.98,2.01,2.03};
1663 double aold_K1[6]={13.8,10.6,6.9,6.7,3.7,4.9};
1664 double aold_K2[6]={8.7,4.6,4.4,1.8,-0.2,-1.8};
1665
1666 TH1D *ratioMuonCorr=new TH1D("ratioMuonCorr","",6,0.5,6.5);
1667 ratioMuonCorr->GetXaxis()->SetBinLabel(1,"0-5%");
1668 ratioMuonCorr->GetXaxis()->SetBinLabel(2,"5-10%");
1669 ratioMuonCorr->GetXaxis()->SetBinLabel(3,"10-20%");
1670 ratioMuonCorr->GetXaxis()->SetBinLabel(4,"20-30%");
1671 ratioMuonCorr->GetXaxis()->SetBinLabel(5,"30-40%");
1672 ratioMuonCorr->GetXaxis()->SetBinLabel(6,"40-50%");
1673 ratioMuonCorr->GetYaxis()->SetTitle("New/Old"); ratioMuonCorr->GetYaxis()->SetTitleOffset(1.8);
1674 ratioMuonCorr->SetTitle("Quartic coefficient");
1675 for(int i=0; i<6; i++){
1676 ratioMuonCorr->SetBinContent(i+1, QuarticFit[1][i]->GetParameter(1)/SF/aold_K1[i]);
1677 ratioMuonCorr->SetBinError(i+1, 1000./QuarticFit[1][i]->GetParameter(1));// 0.015/1.8 for I_K1, 0.03/1.8 for I_K2,
1678 }
1679 ratioMuonCorr->Draw();
1680 */
1681
9b479142 1682 ////////////////////////////
1683 cout<<"Quadratic"<<endl;
1684 for(int cb=0; cb<6; cb++) {
1685 TString *CentLabel2;
1686 if(cb==0) {CentLabel2 = new TString("0-5\\%");}
1687 else if(cb==1) {CentLabel2 = new TString("5-10\\%");}
1688 else if(cb==2) {CentLabel2 = new TString("10-20\\%");}
1689 else if(cb==3) {CentLabel2 = new TString("20-30\\%");}
1690 else if(cb==4) {CentLabel2 = new TString("30-40\\%");}
1691 else {CentLabel2 = new TString("40-50\\%");}
1692 double SF=10.;
3e748846 1693 cout.precision(4);
1694 double I_Sys1 = fabs(QuadraticFit[1][cb]->GetParameter(0)-QuadraticFit[0][cb]->GetParameter(0));// Meth+Det Sys
1695 double I_Sys2 = fabs(QuadraticFit[2][cb]->GetParameter(0)-QuadraticFit[3][cb]->GetParameter(0));// Coulomb Sys
9b479142 1696 double I_SysTotal = sqrt(pow(I_Sys1,2)+pow(I_Sys2,2));
1697 double a_Sys1=fabs(QuadraticFit[1][cb]->GetParameter(1)-QuadraticFit[0][cb]->GetParameter(1));
1698 double a_Sys2 = fabs(QuadraticFit[2][cb]->GetParameter(1)-QuadraticFit[3][cb]->GetParameter(1));
1699 a_Sys1 /= SF; a_Sys2 /= SF;
1700 double a_SysTotal = sqrt(pow(a_Sys1,2)+pow(a_Sys2,2));
1701 I_quadratic_Avg += QuadraticFit[1][cb]->GetParameter(0);
1702 I_quadratic_AvgErr1 += pow(QuadraticFit[0][cb]->GetParError(0),2);
3e748846 1703 I_quadratic_AvgErr2 += I_SysTotal;
9b479142 1704 a_quadratic_Avg += QuadraticFit[1][cb]->GetParameter(1)/SF;
1705 a_quadratic_AvgErr1 += pow(QuadraticFit[0][cb]->GetParError(1)/SF,2);
9b34d0fd 1706 a_quadratic_AvgErr2 += sqrt(pow(a_Sys1,2) + pow(a_Sys2,2));
9b479142 1707 cout.precision(3);
1708 cout<<"$"<<setprecision(2)<<fixed<<CentLabel2->Data()<<"$ & $"<<QuadraticFit[1][cb]->GetParameter(0)<<" \\pm "<<QuadraticFit[0][cb]->GetParError(0)<<" \\pm "<<I_SysTotal<<"$ & $"<<setprecision(1)<<fixed<<QuadraticFit[1][cb]->GetParameter(1)/SF<<" \\pm "<<QuadraticFit[0][cb]->GetParError(1)/SF<<" \\pm "<<a_SysTotal<<"$ \\\\ \\hline"<<endl;
1709 }
9b34d0fd 1710 I_quadratic_Avg /= 6.; I_quadratic_AvgErr1 = sqrt(I_quadratic_AvgErr1)/6.; I_quadratic_AvgErr2 = I_quadratic_AvgErr2/6.;
1711 a_quadratic_Avg /= 6.; a_quadratic_AvgErr1 = sqrt(a_quadratic_AvgErr1)/6.; a_quadratic_AvgErr2 = a_quadratic_AvgErr2/6.;
9b479142 1712 cout<<"$"<<setprecision(2)<<fixed<<"0-50\\%"<<"$ & $"<<I_quadratic_Avg<<" \\pm "<<I_quadratic_AvgErr1<<" \\pm "<<I_quadratic_AvgErr2<<"$ & $"<<setprecision(1)<<fixed<<a_quadratic_Avg<<" \\pm "<<a_quadratic_AvgErr1<<" \\pm "<<a_quadratic_AvgErr2<<"$ \\\\ \\hline"<<endl;
3e748846 1713
9b479142 1714
1715
1716
3e748846 1717
9b479142 1718 ////////////////////////////////////////////////////////////////////////////////
3e748846 1719 // Coulomb Variation
9b479142 1720
1721 int CBOI_SS=1;// 1=5-10%
1722 TCanvas *can6 = new TCanvas("can6", "can6",1800,300,600,900);// 680,300,600,600
1723 can6->SetHighLightColor(2);
1724 gStyle->SetOptFit(0111);
3e748846 1725 can6->SetFillColor(0);//10
9b479142 1726 can6->SetBorderMode(0);
1727 can6->SetBorderSize(2);
1728 can6->SetFrameFillColor(0);
1729 can6->SetFrameBorderMode(0);
1730 can6->SetFrameBorderMode(0);
1731
1732 TPad *pad6 = new TPad("pad6","pad6",0.0,0.0,1.,1.);
1733 gPad->SetGridx(0);
1734 gPad->SetGridy(0);
1735 gPad->SetTickx();
1736 gPad->SetTicky();
9b34d0fd 1737 pad6->SetBottomMargin(0.07);//0.12
9b479142 1738 pad6->SetTopMargin(0.02);//0.05
1739 pad6->SetRightMargin(0.01);//1e-2
1740 pad6->Divide(1,2,0,0);
1741 pad6->Draw();
1742 pad6->cd(1);
3e748846 1743 TLegend *legend7 = new TLegend(.8,.05, .98,.2,NULL,"brNDC");//c3
9b34d0fd 1744 legend7->SetBorderSize(0);
1745 legend7->SetTextFont(TextFont);
1746 legend7->SetTextSize(SizeLegend*SF1);
9b479142 1747 legend7->SetFillColor(0);
1748
1749 //
1750 r3merged[1][KT3Bin][CBOI_SS]->SetMarkerStyle(24); r3merged[1][KT3Bin][CBOI_SS]->SetMarkerColor(2); r3merged[1][KT3Bin][CBOI_SS]->SetLineColor(2);// Omega0, Kt3_1
1751 r3merged[0][KT3Bin][CBOI_SS]->SetMinimum(1.1); r3merged[0][KT3Bin][CBOI_SS]->SetMaximum(2.35);
1752 c3merged[1][1][KT3Bin][CBOI_SS]->SetMarkerStyle(21); c3merged[1][1][KT3Bin][CBOI_SS]->SetMarkerColor(4); c3merged[1][1][KT3Bin][CBOI_SS]->SetLineColor(4);// Omega0, Kt3_1
1753 c3merged[1][0][KT3Bin][CBOI_SS]->SetMinimum(0.995); c3merged[1][0][KT3Bin][CBOI_SS]->SetMaximum(1.065);
1754 //
3e748846 1755 gPad->SetRightMargin(RightMargin); gPad->SetLeftMargin(0.13);
9b34d0fd 1756 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelFont(TextFont); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelFont(TextFont);
1757 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelSize(0); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelSize(SizeLabel*SF1);
9b479142 1758 //
3e748846 1759 r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleOffset(0.8);
9b34d0fd 1760 r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleFont(TextFont);
1761 r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleSize(SizeTitle*SF1);
3e748846 1762 r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitle("#font[12]{r}_{3}");
1763 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetRangeUser(0,0.08);
9b479142 1764 //
3e748846 1765 r3merged[0][KT3Bin][CBOI_SS]->SetMinimum(0.9); r3merged[0][KT3Bin][CBOI_SS]->SetMaximum(2.4);
9b34d0fd 1766 r3merged[0][KT3Bin][CBOI_SS]->SetMarkerSize(1.5);
1767 r3merged[1][KT3Bin][CBOI_SS]->SetMarkerSize(1.7);
9b479142 1768 r3merged[0][KT3Bin][CBOI_SS]->DrawCopy();
1769 r3merged[1][KT3Bin][CBOI_SS]->DrawCopy("same");
1770 ChaoticLimit->Draw("same");
1771 legend7->AddEntry(r3merged[0][KT3Bin][CBOI_SS],"GRS","p");
1772 legend7->AddEntry(r3merged[1][KT3Bin][CBOI_SS],"#Omega_{0}","p");
1773 legend7->Draw("same");
3e748846 1774 TLatex *Specif_5 = new TLatex(0.36,0.94,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");
9b34d0fd 1775 Specif_5->SetNDC();
1776 Specif_5->SetTextFont(TextFont);
1777 Specif_5->SetTextSize(SizeSpecif*SF1);
3e748846 1778 TLatex *Specif2_5 = new TLatex(0.5,0.1,"5-10%");
9b34d0fd 1779 Specif2_5->SetNDC();
1780 Specif2_5->SetTextFont(TextFont);
3e748846 1781 Specif2_5->SetTextSize(SizeLabel*SF1);
9b34d0fd 1782 Specif_5->Draw();
1783 Specif2_5->Draw();
9b479142 1784 TString *KTString = new TString("");
3e748846 1785 if(KT3Bin==0) KTString->Append("0.16<#font[12]{K}_{T,3}<0.3 GeV/#font[12]{c}");
1786 else KTString->Append("0.3<#font[12]{K}_{T,3}<1.0 GeV/#font[12]{c}");
1787 TLatex *Kt3Name_3 = new TLatex(0.36,0.22,KTString->Data());
9b34d0fd 1788 Kt3Name_3->SetNDC();
1789 Kt3Name_3->SetTextFont(TextFont);
3e748846 1790 Kt3Name_3->SetTextSize(SizeLabel*SF1);
9b34d0fd 1791 Kt3Name_3->Draw("same");
9b479142 1792
1793 //
1794 //
1795 pad6->cd(2);
3e748846 1796 SF_correction=0.97;
1797 gPad->SetRightMargin(RightMargin); gPad->SetLeftMargin(0.13); gPad->SetBottomMargin(0.16);
9b34d0fd 1798 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelFont(TextFont); c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelFont(TextFont);
3e748846 1799 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelSize(SizeLabel*SF1*SF_correction); c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelSize(SizeLabel*SF1*SF_correction);
9b34d0fd 1800 c3merged[1][0][KT3Bin][CBOI_SS]->SetMarkerStyle(21);
1801 c3merged[1][1][KT3Bin][CBOI_SS]->SetMarkerStyle(25);
1802 c3merged[1][0][KT3Bin][CBOI_SS]->SetMarkerSize(1.5);
1803 c3merged[1][1][KT3Bin][CBOI_SS]->SetMarkerSize(1.7);
9b479142 1804 //
3e748846 1805 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleOffset(0.9); c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleOffset(0.8);
9b34d0fd 1806 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleFont(TextFont); c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleFont(TextFont);
1807 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleSize(SizeTitle*SF1); c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleSize(SizeTitle*SF1);
3e748846 1808 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/#font[12]{c})");
1809 c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitle("#font[12]{#bf{c}}_{3}^{#pm#pm#mp}");
1810 c3merged[1][0][KT3Bin][CBOI_SS]->GetYaxis()->SetNdivisions(504);
1811 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetNdivisions(503);
9b479142 1812 //
3e748846 1813 c3merged[1][0][KT3Bin][CBOI_SS]->GetXaxis()->SetRangeUser(0,0.08);
1814 c3merged[1][0][KT3Bin][CBOI_SS]->SetMaximum(1.08);
1815 TLegend *legend8 = new TLegend(.8,.8, .98,.95,NULL,"brNDC");//c3
9b34d0fd 1816 legend8->SetBorderSize(0);
1817 legend8->SetTextFont(TextFont);
1818 legend8->SetTextSize(SizeLegend*SF1);
9b479142 1819 legend8->SetFillColor(0);
1820 c3merged[1][0][KT3Bin][CBOI_SS]->DrawCopy();
1821 c3merged[1][1][KT3Bin][CBOI_SS]->DrawCopy("same");
1822 legend8->AddEntry(c3merged[1][0][KT3Bin][CBOI_SS],"GRS","p");
1823 legend8->AddEntry(c3merged[1][1][KT3Bin][CBOI_SS],"#Omega_{0}","p");
1824 legend8->Draw("same");
1825 BaseLine->Draw("same");
1826
1827
9b34d0fd 1828 //can6->SaveAs("SaveFigs/Fig_r3c3_CoulVar.eps");
9b479142 1829 ///////////////////////////////////////////////////
3e748846 1830 // Lambda variation
9b479142 1831
1832 TCanvas *can7 = new TCanvas("can7", "can7",1800,300,600,600);// 680,300,600,600
1833 can7->SetHighLightColor(2);
1834 gStyle->SetOptFit(0111);
3e748846 1835 can7->SetFillColor(0);//10
9b479142 1836 can7->SetBorderMode(0);
1837 can7->SetBorderSize(2);
1838 can7->SetFrameFillColor(0);
1839 can7->SetFrameBorderMode(0);
1840 can7->SetFrameBorderMode(0);
1841
1842 TPad *pad7 = new TPad("pad7","pad7",0.0,0.0,1.,1.);
1843 gPad->SetGridx(0);
1844 gPad->SetGridy(0);
9b34d0fd 1845 //gPad->SetTickx();
1846 //gPad->SetTicky();
1847 pad7->SetBottomMargin(0.07);//0.12
9b479142 1848 pad7->SetTopMargin(0.02);//0.05
1849 pad7->SetRightMargin(0.01);//1e-2
1850 pad7->Draw();
1851 TLegend *legend9 = new TLegend(.75,.7, .95,.85,NULL,"brNDC");//c3
9b34d0fd 1852 legend9->SetBorderSize(0);
1853 legend9->SetTextFont(TextFont);
1854 legend9->SetTextSize(SizeLegend*SF2);
9b479142 1855 legend9->SetFillColor(0);
1856
3e748846 1857 gPad->SetRightMargin(0.01); gPad->SetLeftMargin(0.13); gPad->SetBottomMargin(0.13); gPad->SetTopMargin(0.01);
1858 r3merged[0][KT3Bin][CBOI_SS]->SetMinimum(r3MIN); r3merged[0][KT3Bin][CBOI_SS]->SetMaximum(r3MAX);
9b34d0fd 1859 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelFont(TextFont); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelFont(TextFont);
1860 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetLabelSize(SizeLabel*SF2); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetLabelSize(SizeLabel*SF2);
9b479142 1861 //
9b34d0fd 1862 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleFont(TextFont); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleFont(TextFont);
1863 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleSize(SizeTitle*SF2); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleSize(SizeTitle*SF2);
3e748846 1864 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitle("#font[12]{Q}_{3} (GeV/#font[12]{c})"); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitle("#font[12]{r}_{3}");
1865 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetTitleOffset(0.9); r3merged[0][KT3Bin][CBOI_SS]->GetYaxis()->SetTitleOffset(1.05);
9b479142 1866 TH1D *r3LambdaVaried = (TH1D*)r3merged[0][KT3Bin][CBOI_SS]->Clone();
1867 r3LambdaVaried->SetMarkerStyle(25);
1868 TF1 *r3LambdaSysFit = new TF1("r3LambdaSysFit","pol1",0,0.1);
1869 r3LambdaSysFit->FixParameter(0,LambdaSysPar0[CBOI_SS]);
1870 r3LambdaSysFit->FixParameter(1,LambdaSysPar1[CBOI_SS]);
1871 for(int bin=1; bin<=11; bin++){
1872 double Q3 = (bin-0.5)*0.01;
1873 r3LambdaVaried->SetBinContent(bin, (1-r3LambdaSysFit->Eval(Q3)) * r3LambdaVaried->GetBinContent(bin));
1874 }
3e748846 1875 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetNdivisions(503);
1876 r3merged[0][KT3Bin][CBOI_SS]->GetXaxis()->SetRangeUser(0,0.08);
9b34d0fd 1877
9b479142 1878 r3merged[0][KT3Bin][CBOI_SS]->DrawCopy();
1879 r3LambdaVaried->DrawCopy("same");
1880 ChaoticLimit->Draw("same");
1881 legend9->AddEntry(r3merged[0][KT3Bin][CBOI_SS],"#lambda=0.7","p");
1882 legend9->AddEntry(r3LambdaVaried,"#lambda=0.6","p");
1883 legend9->Draw("same");
9b34d0fd 1884
3e748846 1885 TLatex *Specif_6 = new TLatex(0.36,0.94,"ALICE Pb-Pb #\sqrt{#font[12]{s}_{NN}}=2.76 TeV");
9b34d0fd 1886 Specif_6->SetNDC();
1887 Specif_6->SetTextFont(TextFont);
1888 Specif_6->SetTextSize(SizeSpecif*SF2);
3e748846 1889 TLatex *Specif2_6 = new TLatex(0.5,0.2,"5-10%");
9b34d0fd 1890 Specif2_6->SetNDC();
1891 Specif2_6->SetTextFont(TextFont);
1892 Specif2_6->SetTextSize(SizeLegend*SF2);
1893 Specif_6->Draw();
1894 Specif2_6->Draw();
1895 TString *KTString = new TString("");
3e748846 1896 if(KT3Bin==0) KTString->Append("0.16<#font[12]{K}_{T,3}<0.3 GeV/#font[12]{c}");
1897 else KTString->Append("0.3<#font[12]{K}_{T,3}<1.0 GeV/#font[12]{c}");
1898 TLatex *Kt3Name_4 = new TLatex(0.36,0.28,KTString->Data());
9b34d0fd 1899 Kt3Name_4->SetNDC();
1900 Kt3Name_4->SetTextFont(TextFont);
3e748846 1901 Kt3Name_4->SetTextSize(SizeLabel*SF2);
9b34d0fd 1902 Kt3Name_4->Draw("same");
1903
1904 //can7->SaveAs("SaveFigs/Fig_r3_LamVar.eps");
9b479142 1905
1906 //TString *OutName=new TString("r3_M");
1907 //TString *OutName=new TString("c3mc_M");
1908 //TString *OutName=new TString("c3mc_GaussK3_M");
1909 //*OutName += BOI_1;
1910 //OutName->Append(".png");
1911 //can6->SaveAs(OutName->Data());
3e748846 1912
9b479142 1913
1914
9b479142 1915
1916 /*
1917 // digitized from Wilson's Thesis
1918 double Wilson_Q3[18]={0};
1919 double Wilson_Q3_e[18]={0};
1920 double Wilson_C3[18]={2.02, 1.345, 1.51, 1.514, 1.51, 1.453, 1.394, 1.329, 1.277, 1.228, 1.188, 1.16, 1.132, 1.116, 1.092, 1.076, 1.064, 1.052};
1921 double Wilson_C3_eh[18]={2.97, 1.445, 1.54, 1.53, 1.518, 0,0,0,0,0,0,0,0,0,0,0,0,0};
1922 double Wilson_C3_e[18]={0};
1923 for(int i=0; i<18; i++){
1924 Wilson_Q3[i] = (i+0.5)*0.005 + .010;
1925 Wilson_C3_e[i] = Wilson_C3_eh[i]-Wilson_C3[i];
1926 if(Wilson_C3_e[i] > 1) Wilson_C3_e[i]=0;
1927 }
1928 TGraphErrors *gr_Wilson = new TGraphErrors(18, Wilson_Q3, Wilson_C3, Wilson_Q3_e, Wilson_C3_e);
1929 gr_Wilson->SetMarkerStyle(20);
1930 gr_Wilson->SetMarkerColor(2);
1931 gr_Wilson->SetLineColor(2);
1932 //gr_Wilson->Draw("P");
1933 //legend1->AddEntry(gr_Wilson,"STAR 0-12%: 130 GeV","p");
1934
1935 // digitized from Phys. Rev. Lett. 87 (2001) 82301
1936 double STAR_130GeV_C2[20]={1.21, 1.211, 1.182, 1.135, 1.095, 1.064, 1.0389, 1.023, 1.014, 1.0047, 1.0016, 1.001, 1,1,1,1,1,1,1,1};
1937 double STAR_130GeV_C2_e[20]={0};
1938 double STAR_130GeV_C2_qinv[20]={0};
1939 for(int i=0; i<20; i++) STAR_130GeV_C2_qinv[i] = (i+0.5)*0.005;
1940 TGraphErrors *gr_STAR_C2_130 = new TGraphErrors(20,STAR_130GeV_C2_qinv, STAR_130GeV_C2, 0,0);
1941 gr_STAR_C2_130->SetMarkerStyle(20);
1942 gr_STAR_C2_130->SetMarkerColor(2);
1943 //gr_STAR_C2_130->Draw("P");
1944 //legend1->AddEntry(gr_STAR_C2_130,"STAR 0-12%: 130 GeV","p");
1945
1946 TF1 *line=new TF1("line","3",0,100);
1947 line->SetLineStyle(1);
1948 //line->Draw("same");
1949 //legend1->AddEntry(line,"Chaotic limit #font[12]{c}_{3}","l");
1950 TF1 *line2=new TF1("line2","6",0,100);
1951 line2->SetLineStyle(1);
1952 line2->SetLineColor(4);
1953 //line2->Draw("same");
1954 //legend1->AddEntry(line2,"Chaotic limit C_{3}","l");
1955 TF1 *line3=new TF1("line3","1",0,100);
1956 line3->SetLineStyle(1);
1957 //line3->Draw("same");
1958 //legend1->AddEntry(line3,"Chaotic limit","l");
1959 */
1960 /*
1961 //TLatex *tex = new TLatex(0.001,0.09,"q_{ij} > 0.01 GeV/c");//1.54,2300
1962 //TLatex *tex = new TLatex(0.05,4,"#\lambda=0.4#pm0.04, R=7#pm0.5 fm");// for C3
1963 //TLatex *tex = new TLatex(0.05,1.17,"#\lambda=0.4#pm0.04, R=7#pm0.5 fm");// for C2(+-) 0-5%
1964 //TLatex *tex = new TLatex(0.05,1.17,"#\lambda=0.4#pm0.04, R=4#pm0.5 fm");// for C2(+-) 45-50%
1965 TLatex *tex = new TLatex(0.0005,-.9,"#\lambda=0.4#pm0.04, R=7#pm0.5 fm");// for r3
1966
1967 tex->SetTextSize(.04);
1968 //tex->Draw("same");
1969
1970 TLatex *tex2 = new TLatex(0.002,5.75,"Chaotic Limit C_{3}");
1971 tex2->SetTextSize(.04);
1972 tex2->SetTextColor(4);
1973 //tex2->Draw("same");
1974
1975 TLatex *tex3 = new TLatex(0.002,2.75,"Chaotic Limit #font[12]{c}_{3}");
1976 tex3->SetTextSize(.04);
1977 tex3->SetTextColor(2);
1978 //tex3->Draw("same");
1979
1980 TLatex *tex4 = new TLatex(0.0005,-0.5,"Spherical Gaussian Coulomb WF");
1981 tex4->SetTextSize(.05);
1982 tex4->SetTextColor();
1983 //tex4->Draw("same");
1984
3e748846 1985 TLatex *tex5 = new TLatex(0.0005,.5,"Pb-Pb #sqrt{#font[12]{s}_{NN}}=2.76 TeV: 0-10% Centrality");
9b479142 1986 tex5->SetTextSize(.04);
1987 tex5->SetTextColor();
1988 //tex5->Draw("same");
1989 */
1990 //
1991 /*
1992 TF1 *C2noConstraint=new TF1("C2noConstraint","2-pow(x,2)",0,1);
1993 C2noConstraint->SetLineColor(2);
1994 TF1 *C2Constraint=new TF1("C2Constraint","2-4/5.*pow(x,2)",0,1);
1995 C2Constraint->SetLineColor(4);
1996 C2Constraint->GetXaxis()->SetTitle("G");
1997 C2Constraint->GetYaxis()->SetTitle("C_{2}^{++} Intercept");
1998 C2Constraint->GetXaxis()->SetTitleSize(0.06);
1999 C2Constraint->GetYaxis()->SetTitleSize(0.06);
2000 C2Constraint->GetXaxis()->SetTitleOffset(.8);
2001 C2Constraint->GetYaxis()->SetTitleOffset(.7);
2002 C2Constraint->SetTitle("");
2003 TF1 *C3noConstraint=new TF1("C3noConstraint","6-9*pow(x,2)+4*pow(x,3)",0,1);
2004 C3noConstraint->SetLineColor(2);
2005 //C3noConstraint->SetLineStyle(2);
2006 TF1 *C3Constraint=new TF1("C3Constraint","6-36/5.*pow(x,2)+96/35.*pow(x,3)",0,1);
2007 C3Constraint->SetLineColor(4);
2008 //C3Constraint->SetLineStyle(2);
2009 C3Constraint->GetXaxis()->SetTitle("G");
2010 C3Constraint->GetYaxis()->SetTitle("C_{3}^{+++} Intercept");
2011 C3Constraint->GetXaxis()->SetTitleSize(0.06);
2012 C3Constraint->GetYaxis()->SetTitleSize(0.06);
2013 C3Constraint->GetXaxis()->SetTitleOffset(.8);
2014 C3Constraint->GetYaxis()->SetTitleOffset(.7);
2015 C3Constraint->SetTitle("");
2016 C3Constraint->Draw();
2017 legend1->AddEntry(C3Constraint,"C_{3}^{+++} Charge-Constrained","l");
2018 C3noConstraint->Draw("same");
2019 legend1->AddEntry(C3noConstraint,"C_{3}^{+++} Charge-Unconstrained","l");
2020 //C2Constraint->Draw();
2021 //legend1->AddEntry(C2Constraint,"C_{2}^{++} Charge-Constrained","l");
2022 //C2noConstraint->Draw("same");
2023 //legend1->AddEntry(C2noConstraint,"C_{2}^{++} Charge-Unconstrained","l");
2024 legend1->SetTextSize(.03);
2025 */
9b34d0fd 2026
2027 /*TF1 *r3noConstraint = new TF1("r3noConstraint","(2-6*pow(x,2)+4*pow(x,3))/pow(1-pow(x,2),1.5)",0,1);
9b479142 2028 r3noConstraint->SetLineColor(2);
2029 r3noConstraint->GetXaxis()->SetTitle("G");
2030 r3noConstraint->GetYaxis()->SetTitle("r_{3} Intercept");
2031 r3noConstraint->GetXaxis()->SetTitleSize(0.06);
2032 r3noConstraint->GetYaxis()->SetTitleSize(0.06);
2033 r3noConstraint->GetXaxis()->SetTitleOffset(.8);
2034 r3noConstraint->GetYaxis()->SetTitleOffset(.7);
2035 r3noConstraint->SetTitle("");
9b34d0fd 2036 TF1 *r3Constraint = new TF1("r3Constraint","(2-24/5.*pow(x,2)+96/35.*pow(x,3))/pow(1-4/5.*pow(x,2),1.5)",0,1);
9b479142 2037 r3Constraint->SetLineColor(4);
2038 r3Constraint->GetXaxis()->SetTitle("G");
2039 r3Constraint->GetYaxis()->SetTitle("r_{3} Intercept");
2040 r3Constraint->GetXaxis()->SetTitleSize(0.06);
2041 r3Constraint->GetYaxis()->SetTitleSize(0.06);
2042 r3Constraint->GetXaxis()->SetTitleOffset(.8);
2043 r3Constraint->GetYaxis()->SetTitleOffset(.7);
2044 r3Constraint->SetTitle("");
9b34d0fd 2045 r3Constraint->Draw();
9b479142 2046 //legend1->AddEntry(r3Constraint,"Charge-Constrained","l");
9b34d0fd 2047 r3noConstraint->Draw("same");
9b479142 2048 //legend1->AddEntry(r3noConstraint,"Charge-Unconstrained","l");
9b34d0fd 2049 //cout<<setprecision(5)<<r3noConstraint->Eval(0.3)<<" "<<r3Constraint->Eval(0.3)<<endl;
9b479142 2050 */
2051 //legend1->Draw("same");
2052
2053 //DrawALICELogo(kTRUE, 0.82,0.24,1.12,0.44);// C3(+++)
2054 //DrawALICELogo(kTRUE, 0.32,0.47,0.65,0.7);// C2(+-)
2055 //DrawALICELogo(kTRUE, 0.3,0.18,0.63,0.36);// r3
2056 //can1->SaveAs("test.eps");
2057
2058}
2059//____________________________________________________________________________________________________
2060void DrawALICELogo(Bool_t prel, Float_t x1, Float_t y1, Float_t x2, Float_t y2)
2061{
2062 // revision on July 24th or 25th 2012
2063 // correct for aspect ratio of figure plus aspect ratio of pad (coordinates are NDC!)
2064 x2 = x1 + (y2 - y1) * (466. / 523) * gPad->GetWh() * gPad->GetHNDC() / (gPad->GetWNDC() * gPad->GetWw());
2065 // Printf("%f %f %f %f", x1, x2, y1, y2);
2066
2067 TPad *myPadLogo = new TPad("myPadLogo", "Pad for ALICE Logo", x1, y1, x2, y2);
2068 myPadLogo->SetLeftMargin(0);
2069 myPadLogo->SetTopMargin(0);
2070 myPadLogo->SetRightMargin(0);
2071 myPadLogo->SetBottomMargin(0);
2072 myPadLogo->Draw();
2073 myPadLogo->cd();
2074 TASImage *myAliceLogo = new TASImage((prel) ? "~/Pictures/alice_logo_preliminary.eps" : "~/Pictures/alice_logo_performance.eps");
2075 myAliceLogo->Draw();
2076}
2077