]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGCF/FLOW/papers/arXiv_1207.0900/drawPaperFigure1.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGCF / FLOW / papers / arXiv_1207.0900 / drawPaperFigure1.C
CommitLineData
24373b38 1//=======================================================================//
2//Macro to draw the main results of the MH analysis:
3//i) drawCentralityDependence3P:: Draws the <cos(psi1 + psi2 - 2psi3)
4// vs centrality percentile
5//ii) drawNpartDependence3P:: Draws the <cos(psi1 + psi2 - 2psi3)
6// vs Nparticipants
7//iii) drawNpartDependenceScaled3P:: Draws the Npart*<cos(psi1 + psi2 - 2psi3)
8// vs Nparticipants
9//iv) drawCentralityDependenceRP:: Draws the <cos(psi1 + psi2 - 2PsiRP)
10// vs centrality percentile
11//v) drawNpartDependenceRP:: Draws the <cos(psi1 + psi2 - 2PsiRP)
12// vs Nparticipants
13//vi) drawNpartDependenceScaledRP:: Draws the Npart*<cos(psi1 + psi2 - 2PsiRP)
14// vs Nparticipants
15//=======================================================================//
16
17// #include "SetFlowStyle.C"
18
19Bool_t gPreliminary = kFALSE;
20
21float myMarkerSize = 2.0;
22static int myDarkRed = TColor::GetColor(128,0,0);
23static int myLightRed = TColor::GetColor(128,0,0);
24static int myBlue = 9;
25static int myGreen = kGreen+3;
26
27//+++++++++++++++++++++GLOBAL VARIABLES+++++++++++++++++++++//
28const Int_t nCentralityBins = 9;
29TString strCentralityBins[nCentralityBins] = {"0-5","5-10","10-20",
30 "20-30","30-40","40-50",
31 "50-60","60-70","70-80"};
32//Double_t gCentralityPercentile[nCentralityBins] = {75.,65.,55.,45.,35.,25.,15.,7.5,2.5};
33//Double_t gCentralityPercentile[nCentralityBins] = {2.5,7.5,15.,25.,35.,45.,55.,65.,75.};
34//Double_t gCentralityPercentileError[nCentralityBins] = {0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5};
35Double_t gCentralityPercentileSystematicError[nCentralityBins] = {1.5,1.5,3.0,3.0,3.0,3.0,3.0,3.0,3.0};
36
37Double_t gCentralityPercentileError[nCentralityBins] = {0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5};
38Double_t gCentralityPercentileVZERO[nCentralityBins] = {3.,8.,16.,26.,36.,46.,56.,66.,76.};
39Double_t gCentralityPercentileZDC[nCentralityBins] = {4.25,9.25,18.,28.,38.,48.,58.,68.,78.};
40Double_t gCentralityPercentile[nCentralityBins] = {0.5,5.5,12.,22.,32.,42.,52.,62.,74.};
41Double_t gCentralityPercentileTPC[nCentralityBins] = {1.75,6.75,14.,24.,34.,44.,54.,64.,74.};
42
43//================================ALICE================================//
44//================================Panos-VZERO================================//
45Double_t g3pCorrelatorPanosCorrectedPlusMinus[nCentralityBins] = {1.09401e-07,-1.09401e-06,-1.41212e-05,-2.49265e-05,-6.4138e-05,-2.57942e-05,0.000107605,0.000314927,0.000989089};
46Double_t g3pCorrelatorPanosCorrectedPlusMinusError[nCentralityBins] = {4.05503e-06,6.16342e-06,9.18652e-06,1.42343e-05,2.32727e-05,4.07319e-05,8.00524e-05,0.00012171819,0.0002171819};
47Double_t g3pCorrelatorPanosCorrectedSameCharge[nCentralityBins] = {-1.09401e-07,-9.82253e-06,-4.64251e-05,-9.80566e-05,-0.000147225,-0.000248039,-0.000366027,-0.000370332,-0.000393409};
48Double_t g3pCorrelatorPanosCorrectedSameChargeError[nCentralityBins] = {4.10413e-06,6.23425e-06,9.33561e-06,1.4481e-05,2.3662e-05,4.13592e-05,8.14661e-05,0.000176888,0.000376888};
49//================================Panos-VZERO================================//
50
51//================================Ilya-ZDC================================//
52//Double_t g3pCorrelatorIlyaCorrectedPlusMinus[nCentralityBins] = {4.297763e-05,-2.408162e-05,-1.637881e-05,-4.876254e-06,-2.316564e-05,-7.809004e-05,6.207905e-06,-0.0004647401,-0.0001632693};
53Double_t g3pCorrelatorIlyaCorrectedPlusMinus[nCentralityBins] = {4.297763e+15,-2.408162e-05,-1.637881e-05,-4.876254e-06,-2.316564e-05,-7.809004e-05,6.207905e-06,-0.0004647401,-0.0001632693};
54Double_t g3pCorrelatorIlyaCorrectedPlusMinusError[nCentralityBins] = {3.5165e-05,2.527027e-05,1.044865e-05,1.15894e-05,1.778335e-05,3.373117e-05,7.60854e-05,0.000225265,0.0008951721};
55//Double_t g3pCorrelatorIlyaCorrectedSameCharge[nCentralityBins] = {-5.446542e-05,-3.563536e-05,-5.561554e-05,-9.352835e-05,-0.0001713429,-0.0003199126,-0.0003791019,-0.0006610145,-0.00121334};
56Double_t g3pCorrelatorIlyaCorrectedSameCharge[nCentralityBins] = {-5.446542e+15,-3.563536e-05,-5.561554e-05,-9.352835e-05,-0.0001713429,-0.0003199126,-0.0003791019,-0.0006610145,-0.00121334};
57Double_t g3pCorrelatorIlyaCorrectedSameChargeError[nCentralityBins] = {3.556981e-05,2.425196e-05,1.01993e-05,1.134225e-05,1.746077e-05,3.325377e-05,7.538573e-05,0.0002243383,0.0008956909};
58//================================Ilya-ZDC================================//
59
60//==============================Alexandru-TPC==============================//
61Double_t g3pCorrelatorAlexandruCorrectedPlusMinus[nCentralityBins] = {-6.41377e-06,1.87768e-06,-5.07504e-06,-1.13659e-05,-1.03051e-05,-3.0122e-05,-9.96488e-06,0.000246629,0.00108092};
62Double_t g3pCorrelatorAlexandruCorrectedPlusMinusError[nCentralityBins] = {3.47328e-06,3.67446e-06,3.31881e-06,4.77789e-06,7.23694e-06,1.17143e-05,2.16551e-05,4.5914e-05,0.000123041};
63Double_t g3pCorrelatorAlexandruCorrectedSameCharge[nCentralityBins] = {-1.03186e-05,-2.34825e-05,-5.40217e-05,-9.6981e-05,-0.000159935,-0.000268521,-0.000371931,-0.000445415,-0.000214926};
64Double_t g3pCorrelatorAlexandruCorrectedSameChargeError[nCentralityBins] = {3.47395e-06,3.67532e-06,3.31981e-06,4.77997e-06,7.24162e-06,1.17259e-05,2.16903e-05,4.60457e-05,0.000123736};
65//==============================Alexandru-TPC==============================//
66
67
68//================================ALICE================================//
69//<cos*cos>
70Double_t gCosCosALICEDataSameCharge[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
71Double_t gCosCosALICEDataSameChargeError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
72Double_t gCosCosALICEDataSameChargeSystematicError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
73Double_t gCosCosALICEDataPlusMinus[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
74Double_t gCosCosALICEDataPlusMinusError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
75Double_t gCosCosALICEDataPlusMinusSystematicError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
76
77//<sin*sin>
78Double_t gSinSinALICEDataSameCharge[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
79Double_t gSinSinALICEDataSameChargeError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
80Double_t gSinSinALICEDataSameChargeSystematicError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
81Double_t gSinSinALICEDataPlusMinus[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
82Double_t gSinSinALICEDataPlusMinusError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
83Double_t gSinSinALICEDataPlusMinusSystematicError[nCentralityBins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.};
84//================================ALICE================================//
85
86//================================ALICE================================//
87Double_t g2pCorrelatorALICEDataPlusMinus[nCentralityBins] = {0.000456523, 0.000560778, 0.000731902, 0.00102103, 0.00143198, 0.00212054, 0.00330624, 0.00544651, 0.00919686};
88Double_t g2pCorrelatorALICEDataPlusMinusError[nCentralityBins] = {2.30045e-06,
89 2.81061e-06,
90 2.64844e-06,
91 3.9621e-06,
92 6.07271e-06,
93 9.91051e-06,
94 1.75048e-05,
95 3.37665e-05,
96 7.36336e-05};
97Double_t g2pCorrelatorALICEDataPlusMinusSystematicError[nCentralityBins] = {1.35E-07,
98 2.16E-07,
99 1.55E-07,
100 3.98E-07,
101 1.54E-06,
102 1.49E-06,
103 2.90E-06,
104 2.60E-06,
105 2.27E-05};
106
107Double_t g2pCorrelatorALICEDataSameCharge[nCentralityBins] = {0.000139404,
108 0.000177643,
109 0.000241073,
110 0.000356669,
111 0.000501588,
112 0.000756712,
113 0.00116784,
114 0.00209485,
115 0.00392973};
116Double_t g2pCorrelatorALICEDataSameChargeError[nCentralityBins] = {1.63072e-06,
117 1.99034e-06,
118 1.87724e-06,
119 2.80842e-06,
120 4.29531e-06,
121 7.01318e-06,
122 1.24019e-05,
123 2.39888e-05,
124 5.25678e-05};
125Double_t g2pCorrelatorALICEDataSameChargeSystematicError[nCentralityBins] = {1.38541E-07,
126 2.82213E-07,
127 6.48121E-07,
128 1.35494E-06,
129 2.78063E-06,
130 6.51713E-06,
131 7.92961E-06,
132 5.92337E-06,
133 6.47308E-05};
134//================================ALICE================================//
135Double_t g3pCorrelatorALICEDataSameCharge[nCentralityBins] = {-1.60885e-05,-2.33387e-05,-5.49677e-05,-0.000105864,-0.000166243,-0.000276424,-0.00037915,-0.000425016,-0.000293409};
136Double_t g3pCorrelatorALICEDataSameChargeError[nCentralityBins] = {3.22359e-06,
1373.45502e-06,3.26195e-06,4.85736e-06,7.60362e-06,1.31918e-05,2.61394e-05,6.2836e-05,0.000196171};
138Double_t g3pCorrelatorALICEDataSameChargeSystematicError[nCentralityBins] = {3.58357E-05,1.66674E-05,1.27416E-05,1.98483E-05,3.97548E-05,8.97482E-05,0.00010582,
1390.000233554,0.00031035126};
140
141Double_t g3pCorrelatorALICEDataPlusMinus[nCentralityBins] = {-3.25362e-07,
142 3.86547e-06,
143 -3.41652e-06,
144 -1.40672e-05,
145 -2.64535e-05,
146 -3.86616e-05,
147 2.17488e-05,
148 0.000248817,
149 0.00105866};
150Double_t g3pCorrelatorALICEDataPlusMinusError[nCentralityBins] = {4.40278e-06,
151 4.87877e-06,
152 4.57899e-06,
153 6.81769e-06,
154 1.06697e-05,
155 1.83207e-05,
156 3.62493e-05,
157 8.68075e-05,
158 0.000279468};
159Double_t g3pCorrelatorALICEDataPlusMinusSystematicError[nCentralityBins] = {3.55452E-05, 1.8482E-05,
160 6.31586E-06,
161 1.20367E-05,
162 3.66273E-05,
163 4.12716E-05,
164 8.08092E-05,
165 0.000158143,
166 0.00031038405};
167
168//================================STAR================================//
169Double_t g2pCorrelatorSTARPlusMinus[nCentralityBins] = {0.000345933,
170 0.000426757,
171 0.000546398,
172 0.000743832,
173 0.00102597,
174 0.00140117,
175 0.00195394,
176 0.0027375,
177 0.00363606};
178Double_t g2pCorrelatorSTARPlusMinusError[nCentralityBins] = {2.06688e-06,
179 2.5451e-06,
180 2.43789e-06,
181 3.24192e-06,
182 4.82798e-06,
183 7.43894e-06,
184 1.22702e-05,
185 2.25468e-05,
186 4.74887e-05};
187
188Double_t g2pCorrelatorSTARPlusPlus[nCentralityBins] = {-0.000162705,
189 -0.000184491,
190 -0.000201566,
191 -0.00022444,
192 -0.000262823,
193 -0.000328381,
194 -0.000364104,
195 -0.000444511,
196 -0.000558773};
197Double_t g2pCorrelatorSTARPlusPlusError[nCentralityBins] = {2.85059e-06,
198 3.50367e-06,
199 3.35781e-06,
200 4.46332e-06,
201 6.67312e-06,
202 1.03454e-05,
203 1.72965e-05,
204 3.21068e-05,
205 7.06495e-05};
206//================================STAR================================//
207
208//================================HIJING================================//
209Double_t g2pCorrelatorHIJINGPlusMinus[nCentralityBins] = {-10.,0.00041372,0.000590745,0.000879697,0.00133003,0.00214314,0.00357934,0.00626945,0.0111598};
210Double_t g2pCorrelatorHIJINGPlusMinusError[nCentralityBins] = {0.,1.0236e-06,8.68066e-07,1.33362e-06,2.56541e-06,3.58639e-06,6.34634e-06,1.22583e-05,2.49233e-05};
211
212Double_t g2pCorrelatorHIJINGSameCharge[nCentralityBins] = {-10.,0.00038675,0.000515641,0.000779283,0.001169,0.00190631,0.00321236,0.00564697,0.0104565};
213Double_t g2pCorrelatorHIJINGSameChargeError[nCentralityBins] = {0.,1.73833e-06,1.87754e-06,2.89076e-06,5.52952e-06,7.75963e-06,1.37514e-05,2.6596e-05,5.41739e-05};
214//================================HIJING================================//
215//+++++++++++++++++++++END OF VARIABLES+++++++++++++++++++++//
216
217//_____________________________________________________//
218void drawPaperFigure1() {
219 //Draws the <cos*cos> and <sin*sin> vs centrality percentile
220 //gROOT->LoadMacro("SetFlowStyle.C");
221 //SetFlowStyle();
222 gROOT->LoadMacro("SetPlotStyle.C");
223 SetPlotStyle();
224
225 TGaxis::SetMaxDigits(5);
226
227 Double_t arrOpp[4] = {g3pCorrelatorALICEDataPlusMinus[0],g3pCorrelatorPanosCorrectedPlusMinus[0],g3pCorrelatorIlyaCorrectedPlusMinus[0],g3pCorrelatorAlexandruCorrectedPlusMinus[0]};
228 Double_t arrSame[4] = {g3pCorrelatorALICEDataSameCharge[0],g3pCorrelatorPanosCorrectedSameCharge[0],g3pCorrelatorIlyaCorrectedSameCharge[0],g3pCorrelatorAlexandruCorrectedSameCharge[0]};
229 Printf("(Same charge) Systematic: %3.10lf",TMath::Abs(arrSame[TMath::LocMax(4,arrSame)]-arrSame[TMath::LocMin(3,arrSame)])/2.);
230 Printf("(Opp charge) Systematic: %3.10lf",TMath::Abs(arrOpp[TMath::LocMax(4,arrOpp)]-arrOpp[TMath::LocMin(3,arrOpp)])/2.);
231
232 //Calculate the coscos and sinsin terms and their errors
233 for(Int_t iBin = 0; iBin < nCentralityBins; iBin++) {
234 //coscos (same charge)
235 gCosCosALICEDataSameCharge[iBin] = 0.5*(g2pCorrelatorALICEDataSameCharge[iBin] + g3pCorrelatorALICEDataSameCharge[iBin]);
236 gCosCosALICEDataSameChargeError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataSameChargeError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataSameChargeError[iBin],2));
237 gCosCosALICEDataSameChargeSystematicError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataSameChargeSystematicError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataSameChargeSystematicError[iBin],2));
238
239 //coscos (opposite charge)
240 gCosCosALICEDataPlusMinus[iBin] = 0.5*(g2pCorrelatorALICEDataPlusMinus[iBin] + g3pCorrelatorALICEDataPlusMinus[iBin]);
241 gCosCosALICEDataPlusMinusError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataPlusMinusError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataPlusMinusError[iBin],2));
242 gCosCosALICEDataPlusMinusSystematicError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataPlusMinusSystematicError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataPlusMinusSystematicError[iBin],2));
243
244 //sinsin (same charge)
245 gSinSinALICEDataSameCharge[iBin] = 0.5*(g2pCorrelatorALICEDataSameCharge[iBin] - g3pCorrelatorALICEDataSameCharge[iBin]);
246 gSinSinALICEDataSameChargeError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataSameChargeError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataSameChargeError[iBin],2));
247 gSinSinALICEDataSameChargeSystematicError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataSameChargeSystematicError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataSameChargeSystematicError[iBin],2));
248
249 //sinsin (opposite charge)
250 gSinSinALICEDataPlusMinus[iBin] = 0.5*(g2pCorrelatorALICEDataPlusMinus[iBin] - g3pCorrelatorALICEDataPlusMinus[iBin]);
251 gSinSinALICEDataPlusMinusError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataPlusMinusError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataPlusMinusError[iBin],2));
252 gSinSinALICEDataPlusMinusSystematicError[iBin] = 0.5*TMath::Sqrt(TMath::Power(g2pCorrelatorALICEDataPlusMinusSystematicError[iBin],2) + TMath::Power(g3pCorrelatorALICEDataPlusMinusSystematicError[iBin],2));
253
254 }
255
256 //================================================//
257 //(+-)
258 TGraphErrors *grALICEDataPlusMinus = new TGraphErrors(nCentralityBins,
259 gCentralityPercentile,
260 g3pCorrelatorALICEDataPlusMinus,
261 gCentralityPercentileError,
262 g3pCorrelatorALICEDataPlusMinusError);
263 myTGraphSetUp(grALICEDataPlusMinus,24,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myDarkRed);
264
265 TGraphErrors *grAlexandruCorrectedPlusMinus = new TGraphErrors(nCentralityBins,
266 gCentralityPercentileTPC,
267 g3pCorrelatorAlexandruCorrectedPlusMinus,
268 gCentralityPercentileError,
269 g3pCorrelatorAlexandruCorrectedPlusMinusError);
270 myTGraphSetUp(grAlexandruCorrectedPlusMinus,25,myGreen,myMarkerSize,1,myGreen,2,1001,myGreen);
271
272 TGraphErrors *grPanosCorrectedPlusMinus = new TGraphErrors(nCentralityBins,
273 gCentralityPercentileVZERO,
274 g3pCorrelatorPanosCorrectedPlusMinus,
275 gCentralityPercentileError,
276 g3pCorrelatorPanosCorrectedPlusMinusError);
277 myTGraphSetUp(grPanosCorrectedPlusMinus,26,myBlue,myMarkerSize,1,myBlue,2,1001,myBlue);
278
279 TGraphErrors *grIlyaCorrectedPlusMinus = new TGraphErrors(nCentralityBins,
280 gCentralityPercentileZDC,
281 g3pCorrelatorIlyaCorrectedPlusMinus,
282 gCentralityPercentileError,
283 g3pCorrelatorIlyaCorrectedPlusMinusError);
284 myTGraphSetUp(grIlyaCorrectedPlusMinus,30,1,myMarkerSize,1,1,2,1001,1);
285
286 //================================================//
287 //(--)&(++)
288 TGraphErrors *grALICEDataSameCharge = new TGraphErrors(nCentralityBins,
289 gCentralityPercentile,
290 g3pCorrelatorALICEDataSameCharge,
291 gCentralityPercentileError,
292 g3pCorrelatorALICEDataSameChargeError);
293 myTGraphSetUp(grALICEDataSameCharge,20,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myDarkRed);
294
295 TGraphErrors *grAlexandruCorrectedSameCharge = new TGraphErrors(nCentralityBins,
296 gCentralityPercentileTPC,
297 g3pCorrelatorAlexandruCorrectedSameCharge,
298 gCentralityPercentileError,
299 g3pCorrelatorAlexandruCorrectedSameChargeError);
300 myTGraphSetUp(grAlexandruCorrectedSameCharge,21,myGreen,myMarkerSize,1,myGreen,2,1001,myGreen);
301
302 TGraphErrors *grPanosCorrectedSameCharge = new TGraphErrors(nCentralityBins,
303 gCentralityPercentileVZERO,
304 g3pCorrelatorPanosCorrectedSameCharge,
305 gCentralityPercentileError,
306 g3pCorrelatorPanosCorrectedSameChargeError);
307 myTGraphSetUp(grPanosCorrectedSameCharge,22,myBlue,myMarkerSize,1,myBlue,2,1001,myBlue);
308
309 TGraphErrors *grIlyaCorrectedSameCharge = new TGraphErrors(nCentralityBins,
310 gCentralityPercentileZDC,
311 g3pCorrelatorIlyaCorrectedSameCharge,
312 gCentralityPercentileError,
313 g3pCorrelatorIlyaCorrectedSameChargeError);
314 myTGraphSetUp(grIlyaCorrectedSameCharge,29,1,myMarkerSize,1,1,2,1001,1);
315
316 //================================================//
317 //STAR
318 TGraphErrors *gr2pSTARDataPlusMinus = new TGraphErrors(nCentralityBins,
319 gCentralityPercentileVZERO,
320 g2pCorrelatorSTARPlusMinus,
321 gCentralityPercentileError,
322 g2pCorrelatorSTARPlusMinusError);
323 myTGraphSetUp(gr2pSTARDataPlusMinus,30,myGreen,myMarkerSize+0.5,1,myGreen,2,1001,myBlue);
324
325 TGraphErrors *gr2pSTARDataPlusPlus = new TGraphErrors(nCentralityBins,
326 gCentralityPercentileVZERO,
327 g2pCorrelatorSTARPlusPlus,
328 gCentralityPercentileError,
329 g2pCorrelatorSTARPlusPlusError);
330 myTGraphSetUp(gr2pSTARDataPlusPlus,29,myGreen,myMarkerSize+0.5,1,myDarkRed,2,1001,myGreen);
331 //================================================//
332
333 //================================================//
334 //(+-)
335 TGraphErrors *gr2pALICEDataPlusMinus = new TGraphErrors(nCentralityBins-1,
336 gCentralityPercentileVZERO,
337 g2pCorrelatorALICEDataPlusMinus,
338 gCentralityPercentileError,
339 g2pCorrelatorALICEDataPlusMinusError);
340 myTGraphSetUp(gr2pALICEDataPlusMinus,24,myDarkRed,myMarkerSize,1,myDarkRed,5,1001,myBlue);
341
342 TGraphErrors *gr2pALICEDataPlusMinusSystematic = new TGraphErrors(nCentralityBins-1,
343 gCentralityPercentileVZERO,
344 g2pCorrelatorALICEDataPlusMinus,
345 gCentralityPercentileSystematicError,
346 g2pCorrelatorALICEDataPlusMinusSystematicError);
347 myTGraphSetUp(gr2pALICEDataPlusMinusSystematic,24,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myBlue);
348
349 //(++)&(--)
350 TGraphErrors *gr2pALICEDataSameCharge = new TGraphErrors(nCentralityBins-1,
351 gCentralityPercentileVZERO,
352 g2pCorrelatorALICEDataSameCharge,
353 gCentralityPercentileError,
354 g2pCorrelatorALICEDataSameChargeError);
355 myTGraphSetUp(gr2pALICEDataSameCharge,20,myDarkRed,myMarkerSize,1,myDarkRed,5,1001,myDarkRed);
356
357 TGraphErrors *gr2pALICEDataSameChargeSystematic = new TGraphErrors(nCentralityBins-1,
358 gCentralityPercentileVZERO,
359 g2pCorrelatorALICEDataSameCharge,
360 gCentralityPercentileSystematicError,
361 g2pCorrelatorALICEDataSameChargeSystematicError);
362 myTGraphSetUp(gr2pALICEDataSameChargeSystematic,20,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myDarkRed);
363 //================================================//
364
365 //================================================//
366 //(+-)
367 TGraphErrors *grCosCosALICEDataPlusMinus = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gCosCosALICEDataPlusMinus,gCentralityPercentileError,gCosCosALICEDataPlusMinusError);
368 myTGraphSetUp(grCosCosALICEDataPlusMinus,24,myDarkRed,myMarkerSize,1,myDarkRed,5,1001,myBlue);
369
370 TGraphErrors *grCosCosALICEDataPlusMinusSystematic = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gCosCosALICEDataPlusMinus,gCentralityPercentileError,gCosCosALICEDataPlusMinusSystematicError);
371 myTGraphSetUp(grCosCosALICEDataPlusMinusSystematic,24,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myBlue);
372
373 TGraphErrors *grSinSinALICEDataPlusMinus = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gSinSinALICEDataPlusMinus,gCentralityPercentileError,gSinSinALICEDataPlusMinusError);
374 myTGraphSetUp(grSinSinALICEDataPlusMinus,25,myBlue,myMarkerSize,1,myBlue,5,1001,myBlue);
375
376 TGraphErrors *grSinSinALICEDataPlusMinusSystematic = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gSinSinALICEDataPlusMinus,gCentralityPercentileError,gSinSinALICEDataPlusMinusSystematicError);
377 myTGraphSetUp(grSinSinALICEDataPlusMinusSystematic,25,myBlue,myMarkerSize,1,myBlue,2,1001,myBlue);
378
379
380 //================================================//
381 //(++)&(--)
382 TGraphErrors *grCosCosALICEDataSameCharge = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gCosCosALICEDataSameCharge,gCentralityPercentileError,gCosCosALICEDataSameChargeError);
383 myTGraphSetUp(grCosCosALICEDataSameCharge,20,myDarkRed,myMarkerSize,1,myDarkRed,5,1001,myBlue);
384
385 TGraphErrors *grCosCosALICEDataSameChargeSystematic = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gCosCosALICEDataSameCharge,gCentralityPercentileError,gCosCosALICEDataSameChargeSystematicError);
386 myTGraphSetUp(grCosCosALICEDataSameChargeSystematic,20,myDarkRed,myMarkerSize,1,myDarkRed,2,1001,myBlue);
387
388 TGraphErrors *grSinSinALICEDataSameCharge = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gSinSinALICEDataSameCharge,gCentralityPercentileError,gSinSinALICEDataSameChargeError);
389 myTGraphSetUp(grSinSinALICEDataSameCharge,21,myBlue,myMarkerSize,1,myBlue,5,1001,myBlue);
390
391 TGraphErrors *grSinSinALICEDataSameChargeSystematic = new TGraphErrors(nCentralityBins,gCentralityPercentileVZERO,gSinSinALICEDataSameCharge,gCentralityPercentileError,gSinSinALICEDataSameChargeSystematicError);
392 myTGraphSetUp(grSinSinALICEDataSameChargeSystematic,21,myBlue,myMarkerSize,1,myBlue,2,1001,myBlue);
393
394 //================================================//
395 //HIJING (+-)
396 TGraphErrors *gr2pHIJINGPlusMinus = new TGraphErrors(nCentralityBins-1,
397 gCentralityPercentileVZERO,
398 g2pCorrelatorHIJINGPlusMinus,
399 gCentralityPercentileError,
400 g2pCorrelatorHIJINGPlusMinusError);
401 myTGraphSetUp(gr2pHIJINGPlusMinus,26,myBlue,myMarkerSize,1,myBlue,5,1001,myBlue);
402
403 //(++)&(--)
404 TGraphErrors *gr2pHIJINGSameCharge = new TGraphErrors(nCentralityBins-1,
405 gCentralityPercentileVZERO,
406 g2pCorrelatorHIJINGSameCharge,
407 gCentralityPercentileError,
408 g2pCorrelatorHIJINGSameChargeError);
409 myTGraphSetUp(gr2pHIJINGSameCharge,22,myBlue,myMarkerSize,1,myBlue,5,1001,myBlue);
410 //================================================//
411
412 //_____________________________________________________//
413 //Draw the results
414
415 //====================================//
416 //<cos(psi1+psi2-2phi3)> vs centrality
417 //TH2F *gEmpty1 = new TH2F("gEmpty1",";centrality, %;#LT cos(#Delta #phi_{#alpha}) cos(#Delta #phi_{#beta}) #GT", nCentralityBins,0,80,1000,-1.5e-03,0.01);
418 TH2F *gEmpty1 = new TH2F("gEmpty1",";centrality, %;", nCentralityBins,0,72,1000,-1.5e-03,0.015);
419 gEmpty1->SetStats(kFALSE);
420 gEmpty1->GetYaxis()->SetTitleSize(0.07);
421 gEmpty1->GetYaxis()->SetTitleOffset(0.95);
422 gEmpty1->GetYaxis()->SetNdivisions(10);
423 gEmpty1->GetXaxis()->SetNdivisions(10);
424
425 TF1 *f1 = new TF1("f1","0",0,1000);
426 f1->SetLineColor(1); f1->SetLineStyle(1); f1->SetLineWidth(1);
427
428 //================================================//
429 TCanvas *c1 = new TCanvas("c1","Centrality dependence",0,0,700,1200);
430 c1->SetFillColor(10); c1->SetHighLightColor(10);
431 c1->Divide(1,3,0.99,0.0,10);
432 //TPad *myPad1 = new TPad("myPad1","myPad1",0,0.667,1,1);
433 //c1->cd();
434 //myPadSetUp(myPad1,0.13,0.07,0.04,0.0,10);
435 //myPad1->Draw();
436 c1->cd(1)->SetLeftMargin(0.19);
437 c1->cd(1)->SetTopMargin(0.083);
438 c1->cd(1)->SetRightMargin(0.01);
439 //myPad1->cd();
440 gEmpty1->GetYaxis()->SetLabelSize(0.075);
441 gEmpty1->GetXaxis()->SetLabelSize(0.075);
442 gEmpty1->GetYaxis()->SetTitleSize(0.095);
443 gEmpty1->GetYaxis()->SetTitleOffset(1.01);
444 gEmpty1->GetXaxis()->SetTitleSize(0.075);
445 gEmpty1->GetYaxis()->SetNdivisions(4);
446 gEmpty1->GetXaxis()->SetNdivisions(0);
447 gEmpty1->GetYaxis()->SetRangeUser(-9.e-04,5.e-04);
448 gEmpty1->GetYaxis()->SetTitle("#LT cos(#phi_{#alpha} + #phi_{#beta} - 2#Psi_{RP}) #GT");
449 //gEmpty1->GetYaxis()->SetNdivisions(3);
450 gEmpty1->DrawCopy();
451 f1->Draw("same");
452 grAlexandruCorrectedPlusMinus->Draw("P,Z");
453 grAlexandruCorrectedSameCharge->Draw("P,Z");
454 grIlyaCorrectedPlusMinus->Draw("P,Z");
455 grIlyaCorrectedSameCharge->Draw("eZ,P");
456 grPanosCorrectedPlusMinus->Draw("eZ,P");
457 grPanosCorrectedSameCharge->Draw("eZ,P");
458 grALICEDataPlusMinus->Draw("eZ,P");
459 grALICEDataSameCharge->Draw("eZ,P");
460
461 TLegend *legend1 = new TLegend(0.23,0.03,0.85,0.35,"","brNDC");
462 myLegendSetUp(legend1,0.06);
463 legend1->SetNColumns(2);
464 legend1->AddEntry(grALICEDataSameCharge," ","lp");
465 legend1->AddEntry(grALICEDataPlusMinus," TPC (cumulants) ","lp");
466 legend1->AddEntry(grAlexandruCorrectedSameCharge," ","lp");
467 legend1->AddEntry(grAlexandruCorrectedPlusMinus," TPC","lp");
468 legend1->AddEntry(grPanosCorrectedSameCharge," ","lp");
469 legend1->AddEntry(grPanosCorrectedPlusMinus," VZERO","lp");
470 legend1->AddEntry(grIlyaCorrectedSameCharge," ","lp");
471 legend1->AddEntry(grIlyaCorrectedPlusMinus," ZDC","lp");
472 legend1->Draw();
473 TLatex *myText1 = new TLatex();
474 myText1->SetNDC();
475 myText1->SetTextSize(0.07);
476 myText1->SetTextColor(1);
477 myText1->DrawLatex(0.23,0.37,"same");
478 myText1->SetTextColor(1);
479 myText1->DrawLatex(0.35,0.37,"opp.");
480 myText1->SetTextColor(1);
481 myText1->DrawLatex(0.37,0.85,"ALICE Pb-Pb @ #sqrt{s_{NN}} = 2.76 TeV");
482 myText1->DrawLatex(0.285,0.85,"(a)");
483
484 //================================================//
485 //TPad *myPad2 = new TPad("myPad2","myPad2",0,0.333,1,0.667);
486 //c1->cd();
487 //myPadSetUp(myPad2,0.13,0.00,0.04,0.00,10);
488 //myPad2->Draw();
489 //c1->cd();
490 //myPad2->cd();
491 c1->cd(2)->SetLeftMargin(0.19); c1->cd(2)->SetRightMargin(0.01);
492 gEmpty1->GetYaxis()->SetLabelSize(0.075);
493 gEmpty1->GetXaxis()->SetLabelSize(0.075);
494 gEmpty1->GetYaxis()->SetTitleSize(0.095);
495 gEmpty1->GetXaxis()->SetTitleSize(0.075);
496 gEmpty1->GetXaxis()->SetNdivisions(0);
497 gEmpty1->GetYaxis()->SetNdivisions(5);
498 gEmpty1->GetYaxis()->SetTitle("#LT cos(#phi_{#alpha}-#phi_{#beta}) #GT");
499 gEmpty1->GetYaxis()->CenterTitle();
500 gEmpty1->GetYaxis()->SetRangeUser(-1.5e-03,0.007);
501 gEmpty1->DrawCopy();
502 f1->Draw("same");
503 gr2pSTARDataPlusMinus->Draw("P,eZ");
504 gr2pSTARDataPlusPlus->Draw("P,eZ");
505 gr2pALICEDataPlusMinusSystematic->Draw("E3");
506 gr2pALICEDataPlusMinusSystematic->Draw("L,X0,same");
507 gr2pALICEDataPlusMinus->Draw("P,eZ");
508 gr2pALICEDataSameChargeSystematic->Draw("E3");
509 gr2pALICEDataSameChargeSystematic->Draw("L,X0,same");
510 gr2pALICEDataSameCharge->Draw("P,eZ");
511
512 //gr2pHIJINGPlusMinus->Draw("P,eZ");
513 //gr2pHIJINGSameCharge->Draw("P,eZ");
514
515 TLatex *myText2 = new TLatex();
516 myText2->SetNDC();
517 myText2->SetTextSize(0.07);
518 myText2->SetTextColor(1);
519 myText2->DrawLatex(0.22,0.86,"same");
520 myText2->SetTextColor(1);
521 myText2->DrawLatex(0.325,0.86,"opp.");
522 myText2->DrawLatex(0.285,0.93,"(b)");
523
524 TLegend *legend2 = new TLegend(0.22,0.63,0.89,0.83,"","brNDC");
525 myLegendSetUp(legend2,0.06);
526 legend2->SetNColumns(2);
527 legend2->AddEntry(gr2pALICEDataSameCharge," ","PL");
528 legend2->AddEntry(gr2pALICEDataPlusMinus," ALICE Pb-Pb @ #sqrt{s_{NN}} = 2.76 TeV","PL");
529 legend2->AddEntry(gr2pSTARDataPlusPlus," ","P");
530 legend2->AddEntry(gr2pSTARDataPlusMinus," STAR Au-Au @ #sqrt{s_{NN}} = 0.2 TeV","P");
531 legend2->Draw();
532
533 //================================================//
534 //cos*cos & sin*sin
535 //TPad *myPad3 = new TPad("myPad3","myPad3",0,0,1,0.333);
536 //c1->cd();
537 //myPadSetUp(myPad3,0.13,0.00,0.04,0.15,10);
538 //myPad3->Draw();
539 //c1->cd();
540 //myPad3->cd();
541 //gEmpty1->GetYaxis()->SetTitle("#LT cos(#phi_{#alpha}-#phi_{#beta}) #GT #pm #LT cos(#phi_{#alpha} + #phi_{#beta} - 2#Psi_{RP}) #GT");
542 c1->cd(3)->SetLeftMargin(0.19); c1->cd(3)->SetRightMargin(0.01);
543 gEmpty1->GetYaxis()->SetLabelSize(0.065);
544 gEmpty1->GetXaxis()->SetLabelSize(0.065);
545 gEmpty1->GetYaxis()->SetTitleSize(0.065);
546 gEmpty1->GetXaxis()->SetTitleSize(0.075);
547 gEmpty1->GetXaxis()->SetNdivisions(10);
548 gEmpty1->GetYaxis()->SetTitle("");
549 gEmpty1->GetYaxis()->SetRangeUser(0.0,0.0035);
550 gEmpty1->GetYaxis()->SetNdivisions(4);
551 gEmpty1->DrawCopy();
552 f1->Draw("same");
553 gStyle->SetErrorX(0);
554
555 grCosCosALICEDataSameCharge->Draw("P,eZ");
556 grCosCosALICEDataSameChargeSystematic->Draw("Z");
557 grCosCosALICEDataPlusMinus->Draw("P,eZ");
558 grCosCosALICEDataPlusMinusSystematic->Draw("Z");
559 grSinSinALICEDataSameCharge->Draw("P,eZ");
560 grSinSinALICEDataSameChargeSystematic->Draw("Z");
561 grSinSinALICEDataPlusMinus->Draw("P,eZ");
562 grSinSinALICEDataPlusMinusSystematic->Draw("Z");
563
564 TLatex *myText3 = new TLatex();
565 myText3->SetNDC();
566 myText3->SetTextSize(0.055);
567 myText3->SetTextColor(1);
568 myText3->DrawLatex(0.22,0.86,"same");
569 myText3->SetTextColor(1);
570 myText3->DrawLatex(0.325,0.86,"opp.");
571 myText3->DrawLatex(0.285,0.93,"(c)");
572 myText3->DrawLatex(0.37,0.93,"ALICE Pb-Pb @ #sqrt{s_{NN}} = 2.76 TeV");
573
574 TLegend *legend3 = new TLegend(0.22,0.63,0.89,0.83,"","brNDC");
575 myLegendSetUp(legend3,0.06);
576 legend3->SetNColumns(2);
577 legend3->AddEntry(grCosCosALICEDataSameCharge," ","P");
578 legend3->AddEntry(grCosCosALICEDataPlusMinus," #LT cos(#Delta #phi_{#alpha}) cos(#Delta #phi_{#beta}) #GT","P");
579 legend3->AddEntry(grSinSinALICEDataSameCharge," ","P");
580 legend3->AddEntry(grSinSinALICEDataPlusMinus," #LT sin(#Delta #phi_{#alpha}) sin(#Delta #phi_{#beta}) #GT","P");
581 legend3->Draw();
582
583 if(gPreliminary) {
584 TLatex *alice = new TLatex(0.23,0.38,"Preliminary");
585 alice->SetNDC();
586 alice->SetTextColor(myDarkRed);
587 alice->SetTextSize(0.035);
588 alice->SetLineWidth(2);
589 alice->Draw();
590
591 TPad *myPadLogo = new TPad("myPadLogo",
592 "Pad for ALICE Logo",0.26,0.41,0.36,0.51);
593 //myPadLogo->SetFillColor(2); // color to first figure out where is the pad then comment !
594 myPadSetUp(myPadLogo,0,0,0,0);
595 myPadLogo->Draw();
596 myPadLogo->cd();
597 TASImage *myAliceLogo = new TASImage("alice_logo_transparent.png");
598 myAliceLogo->Draw();
599 }
600
601 c1->SaveAs("figure1.eps");
602 c1->SaveAs("figure1.pdf");
603 c1->SaveAs("figure1.png");
604}
605
606void SetPadSetUp(TPad *currentPad,
607 float currentLeft=0.11, float currentTop=0.04,
608 float currentRight=0.04, float currentBottom=0.15){
609 currentPad->SetLeftMargin(currentLeft);
610 currentPad->SetTopMargin(currentTop);
611 currentPad->SetRightMargin(currentRight);
612 currentPad->SetBottomMargin(currentBottom);
613 currentPad->SetFillColor(10);
614
615 return;
616}
617
618void myTGraphSetUp(TGraphErrors *currentGraph=0,
619 int myMarkerStyle=8,
620 int myMarkerColor=1,
621 float myMarkerSize=1,
622 int myLineStyle=1,
623 int myLineColor=1,
624 float myLineWidth=1,
625 int myFillStyle =1001,
626 int myFillColor =1) {
627 currentGraph->SetMarkerStyle(myMarkerStyle);
628 currentGraph->SetMarkerColor(myMarkerColor);
629 currentGraph->SetMarkerSize(myMarkerSize);
630 currentGraph->SetLineColor(myLineColor);
631 currentGraph->SetLineStyle(myLineStyle);
632 currentGraph->SetLineWidth(myLineWidth);
633 currentGraph->SetFillStyle(myFillStyle);
634 currentGraph->SetFillColor(myFillColor);
635// currentGraph->Set();
636}
637
638void myLegendSetUp(TLegend *currentLegend=0,float currentTextSize=0.07){
639 currentLegend->SetTextFont(42);
640 currentLegend->SetBorderSize(0);
641 currentLegend->SetFillStyle(0);
642 currentLegend->SetFillColor(0);
643 currentLegend->SetMargin(0.25);
644 currentLegend->SetTextSize(currentTextSize);
645 currentLegend->SetEntrySeparation(0.5);
646 return;
647}