]>
Commit | Line | Data |
---|---|---|
24373b38 | 1 | // S. Voloshin 03-oct-2010 |
2 | // modified from plotV2edepBevalacRhicSqrtS.C - from Art | |
3 | { | |
4 | float m_p=0.940; | |
5 | ||
6 | // gROOT->Reset(); | |
7 | ||
8 | // gROOT->SetStyle("Bold"); | |
9 | //gROOT->SetStyle("Plain"); | |
10 | // gStyle->SetTitleTextColor(kBlack); | |
11 | // gStyle->SetLabelColor(kBlack,"X"); | |
12 | // gStyle->SetLabelColor(kBlack,"Y"); | |
13 | ||
14 | int cCERES=kOrange+7; | |
15 | int cPHENIX=kYellow+3; | |
16 | int cFOPI=kGreen+2; | |
17 | int cE895=kMagenta+2; | |
18 | int cAlice=2; | |
19 | ||
20 | TCanvas *canvas = new TCanvas("v2edep","v2edep",10,10,750,600); | |
21 | canvas->cd(); | |
22 | ||
23 | hist = new TH1F("v2 vs sqrt(s)","v2 vs sqrt(s)", 1, 1., 10000.); | |
24 | hist->SetLineColor(0); | |
25 | ||
26 | TAxis *axis = hist->GetXaxis(); | |
27 | axis->SetTitle("#sqrt{s_{NN}} (GeV)"); | |
28 | //axis->CenterTitle(kTRUE); | |
29 | axis->SetTitleSize(0.05); | |
30 | axis->SetLabelSize(0.04); | |
31 | axis->SetTitleOffset(1.5); | |
32 | axis = hist->GetYaxis(); | |
33 | axis->SetTitle("v_{2}(y^{*}=0)"); | |
34 | axis->SetTitle("v_{2}"); | |
35 | axis->SetTitleOffset(1.6); | |
36 | //axis->CenterTitle(kTRUE); | |
37 | axis->SetTitleSize(0.05); | |
38 | axis->SetLabelSize(0.04); | |
39 | ||
40 | hist->SetStats(0); | |
41 | // hist->SetTitle("Elliptic Flow"); | |
42 | hist->SetTitle(""); | |
43 | hist->SetMaximum(0.08); | |
44 | hist->SetMinimum(-0.085); | |
45 | //hist->SetLabelOffset(-0.01,"X"); | |
46 | gPad->SetTopMargin(.1); | |
47 | gPad->SetBottomMargin(.2); | |
48 | gPad->SetLeftMargin(.2); | |
49 | gPad->SetRightMargin(.1); | |
50 | hist->Draw(); | |
51 | canvas->SetLogx(); | |
52 | hist->Draw(); | |
53 | ||
54 | TLine *tline=new TLine(1.,0.,10000.,0.); | |
55 | tline->SetLineWidth(1); | |
56 | tline->SetLineStyle(7); //wide dash | |
57 | tline->Draw("same"); | |
58 | ||
59 | ||
60 | ||
61 | // CERES ---------------------------------------------------------------------------- | |
62 | float xxr[3]={8.7,12.3,17.}; //sqrt(s) | |
63 | //float xxr[3]={9.57,18.7}; | |
64 | //float yyr[3]={0.028,0.031}; | |
65 | //float yyr[3]={0.028,0.031}; // taken from QM01 and INPC01 | |
66 | //float yyr[3]={0.028,0.03,0.031}; // QM01 and INPC01, 80GeV interpolated | |
67 | //float yyr[3]={0.025,0.035,0.04}; // values integrated over phase space | |
68 | float yyr[3]={0.026,0.035,0.04}; // changed to those values after request by H. Appleshaeuser | |
69 | //CERES data from yugoslavian conference | |
70 | //float xxr[3]={9.57,13.5,18.7}; | |
71 | //float yyr[3]={0.028,0.033,0.036}; | |
72 | ||
73 | float exr[3]={0.,0.,0.}; | |
74 | float exr[3]={0.,0.,0.}; | |
75 | float eyr[3]={0.,0.,0.}; | |
76 | //float eyr[3]={0.0025,0.0025,0.0025}; | |
77 | float eyr[3]={0.005,0.005,0.005}; // changed to those values after request by H. Appleshaeuser | |
78 | ||
79 | //float eyr[3]={0.004,0.005,0.004}; | |
80 | ||
81 | ceres = new TGraphAsymmErrors(3,xxr,yyr,exr,exr,eyr,eyr); | |
82 | ceres->SetTitle(""); | |
83 | ceres->SetMarkerColor(cCERES); | |
84 | ceres->SetMarkerStyle(kOpenCircle); | |
85 | ceres->SetMarkerSize(1.2); | |
86 | ceres->SetLineWidth(2); | |
87 | ceres->SetLineColor(cCERES); | |
88 | ceres->Draw("p"); | |
89 | ||
90 | // NA49 -------------------------------------------------------------------------------- | |
91 | float xxd[4]={8.7,17.,8.7,17.}; //sqrt(s) | |
92 | // float yyd[3]={0.029867,0.0324}; //mod method with Full Acceptance (40) | |
93 | //float yyd[4]={0.023676,0.0324,0.023676,0.0324}; //mod method with Cos Cut (40) | |
94 | float yyd[4]={0.023676,0.0309,0.023676,0.0309}; //mod method with Cos Cut (40) '96 and '00 data (158) | |
95 | //float yyd[4]={0.0452,0.0379,0.0452,0.0379}; //mod method with Cos Cut (40) centrality 5; | |
96 | //float yyd[4]={0.0411,0.0406,0.0411,0.0406}; //mod method with Cos Cut (40) centrality 6; | |
97 | //float yyd[4]={0.0265,0.0341,0.0265,0.0341}; //mod method with Cos Cut (40) centrality 4; | |
98 | //float yyd[4]={0.0206,0.0276,0.0206,0.0276}; //mod method with Cos Cut (40) centrality 3; | |
99 | //float yyd[4]={-0.00495,0.0189,-0.00495,0.0189}; //mod method with Cos Cut (40) centrality 2; | |
100 | //hist->SetTitle("std method: centrality 2"); | |
101 | float exd[4]={0.,0.,0.}; | |
102 | // float eyd[3]={0.001734,0.000716}; //mod method with Full Acceptance (40) | |
103 | //float eyd[4]={0.00276,0.000716,0.005,0.005}; //mod method with Cos Cut (40) | |
104 | float eyd[4]={0.00276,0.00054,0.005,0.005}; //mod method with Cos Cut (40) '96 and '00 data (158) | |
105 | data = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); | |
106 | data->SetTitle(""); | |
107 | data->SetMarkerColor(kRed); | |
108 | data->SetMarkerStyle(kFullSquare); | |
109 | data->SetMarkerSize(1.2); | |
110 | data->SetLineWidth(2); | |
111 | data->SetLineColor(kRed); | |
112 | //data->Draw("p"); | |
113 | ||
114 | // NA49 | |
115 | float xxd[4]={8.7,17.}; //sqrt(s) | |
116 | float yyd[4]={0.023676,0.0324}; //mod method with Cos Cut (40) | |
117 | float exd[4]={0.,0.,0.}; | |
118 | float eyd[4]={0.01,0.005}; //mod method with Cos Cut (40) | |
119 | dataSys = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); | |
120 | dataSys->SetTitle(""); | |
121 | //dataSys->SetMarkerColor(kRed); | |
122 | //dataSys->SetMarkerStyle(kFullSquare); | |
123 | dataSys->SetMarkerSize(1.2); | |
124 | dataSys->SetLineWidth(2); | |
125 | dataSys->SetLineColor(kRed); | |
126 | //dataSys->Draw("[]"); | |
127 | ||
128 | ||
129 | // NA49 cumulant | |
130 | float xxd[2]={7.91,15.45}; //sqrt(s) | |
131 | float yyd[2]={0.02184,0.0291746}; //order 2 | |
132 | //float yyd[2]={0.034883,0.027746}; //order 4 | |
133 | float exd[3]={0.,0.,0.}; | |
134 | //float eyd[3]={0.0105334,0.0026345}; //order 4 | |
135 | float eyd[3]={0.001173,0.000526268}; //order 2 | |
136 | cumul = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); | |
137 | cumul->SetTitle(""); | |
138 | cumul->SetMarkerColor(kRed); | |
139 | cumul->SetMarkerStyle(kFullSquare); | |
140 | cumul->SetMarkerSize(1.2); | |
141 | cumul->SetLineWidth(2); | |
142 | cumul->SetLineColor(kBlue); | |
143 | cumul->Draw("p"); | |
144 | ||
145 | // STAR ------------------------------------------------------------------------------ | |
146 | float xxd[2]={130.,200.}; //sqrt(s) | |
147 | float yyd[2]={0.0426,0.0478}; | |
148 | float exd[2]={0.}; | |
149 | float exd[2]={0.}; | |
150 | float eyd[2]={0.}; | |
151 | float eyd[2]={0.0026,0.0026}; | |
152 | star = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); | |
153 | star->SetTitle(""); | |
154 | star->SetMarkerColor(kRed); | |
155 | star->SetMarkerStyle(30); | |
156 | star->SetMarkerSize(1.6); | |
157 | star->SetLineWidth(2); | |
158 | star->SetLineColor(kRed); | |
159 | star->Draw("p"); | |
160 | ||
161 | ||
162 | // Phenix -------------------------------------------------------------------------- | |
163 | float xxd[1]={220}; //sqrt(s) | |
164 | float yyd[1]={0.054}; | |
165 | float exd[1]={0.}; | |
166 | float exd[1]={0.}; | |
167 | float eyd[1]={0.}; | |
168 | float eyd[1]={0.0041}; | |
169 | phenix = new TGraphAsymmErrors(1,xxd,yyd,exd,exd,eyd,eyd); | |
170 | phenix->SetTitle(""); | |
171 | phenix->SetMarkerColor(cPHENIX); | |
172 | phenix->SetMarkerStyle(kOpenSquare); | |
173 | phenix->SetMarkerSize(1.2); | |
174 | phenix->SetLineWidth(2); | |
175 | phenix->SetLineColor(cPHENIX); | |
176 | phenix->Draw("p"); | |
177 | ||
178 | ||
179 | // Phobos ----------------------------------------------------------------------- | |
180 | float xxd[2]={117.,180.}; //sqrt(s) | |
181 | float yyd[2]={0.048,0.051}; | |
182 | float exd[2]={0.}; | |
183 | float exd[2]={0.}; | |
184 | float eyd[2]={0.}; | |
185 | float eyd[2]={0.005,0.005}; | |
186 | phobos = new TGraphAsymmErrors(2,xxd,yyd,exd,exd,eyd,eyd); | |
187 | phobos->SetTitle(""); | |
188 | phobos->SetMarkerColor(kBlue); | |
189 | phobos->SetMarkerStyle(28); | |
190 | phobos->SetMarkerSize(1.2); | |
191 | phobos->SetLineWidth(2); | |
192 | phobos->SetLineColor(kBlue); | |
193 | phobos->Draw("p"); | |
194 | ||
195 | ||
196 | //E877 ------------------------------------------------------------------------- | |
197 | //float xxhh[4]={4.75,17.,130.,200.}; | |
198 | //float yyhh[4]={8.2,9.3,10.,9.4}; | |
199 | float xxhh[1]={4.75}; // sqrt(s) | |
200 | float yyhh[1]={0.019}; | |
201 | float exhh[1]={0.}; | |
202 | float exhh[1]={0.}; | |
203 | float eyhh[1]={0.}; | |
204 | float eyhh[1]={0.002}; | |
205 | e877 = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh,eyhh); | |
206 | e877->SetTitle(""); | |
207 | e877->SetMarkerColor(kBlue); | |
208 | e877->SetMarkerStyle(34); | |
209 | e877->SetMarkerSize(1.2); | |
210 | e877->SetLineWidth(2); | |
211 | e877->SetLineColor(kBlue); | |
212 | e877->Draw("p"); | |
213 | ||
214 | //E895 ------------------------------------------------------------------------ | |
215 | //float xxhq[4]={4.75,17.,130.,200.}; | |
216 | //float yyhq[4]={9.1.,10.,9.8,8.4}; | |
217 | float xxhq[4]={2.68,3.32,3.83,4.24}; // sqrt(s) | |
218 | float yyhq[4]={-0.05,-0.005,0.01,0.015}; | |
219 | float exhq[4]={0.,0.,0.}; | |
220 | float eyhq[4]={0.004,0.003,0.003,0.004}; | |
221 | e895 = new TGraphAsymmErrors(4,xxhq,yyhq,exhq,exhq,eyhq,eyhq); | |
222 | e895->SetTitle(""); | |
223 | e895->SetMarkerColor(cE895); | |
224 | e895->SetMarkerStyle(kFullTriangleUp); | |
225 | e895->SetMarkerSize(1.2); | |
226 | e895->SetLineWidth(2); | |
227 | e895->SetLineColor(cE895); | |
228 | e895->Draw("p"); | |
229 | //Plastic Ball ---------------------------------------------------- | |
230 | //float xxhq[4]={4.75,17.,130.,200.}; | |
231 | //float yyhq[4]={9.1.,10.,9.8,8.4}; | |
232 | float xxPB[6]={.15,.2,.4,.6,.8,1.2}; | |
233 | float yyPB[6]={-.035,-.08,-.09,-.06,-.04,-.03}; | |
234 | float exPB[6]={0.,0.,0.}; | |
235 | float eyPB[6]={0.,0.,0.}; | |
236 | plasticBall = new TGraphAsymmErrors(6,xxPB,yyPB,exPB,exPB,eyPB,eyPB); | |
237 | plasticBall->SetTitle(""); | |
238 | plasticBall->SetMarkerColor(kBlue); | |
239 | plasticBall->SetMarkerStyle(kFullTriangleDown); | |
240 | plasticBall->SetMarkerSize(1.2); | |
241 | plasticBall->SetLineWidth(1); | |
242 | plasticBall->SetLineColor(kBlue); | |
243 | //plasticBall->Draw("p"); | |
244 | ||
245 | ||
246 | //Fopi------------------------------------------------------------ | |
247 | float xxFopi[10]={0.09,0.12,0.15,0.25,0.4,0.6,0.8,1.0,1.2,1.49}; | |
248 | float yyFopi[10]={0.07456,0.02847,-0.00774,-0.05784,-0.08200,-0.07087,-0.06845,-0.06327,-0.05523,-0.04344}; | |
249 | float exFopi[10]={0.,0.,0.}; | |
250 | float eyFopi[10]={0.00746,0.00285,0.00077,0.00578,0.00656,0.00850,0.00684,0.00823,0.00828,0.00956}; | |
251 | for (int i=0;i<10;i++) xxFopi[i] = sqrt((xxFopi[i]+m_p)*m_p*2. +2.*m_p*m_p); | |
252 | ||
253 | fopi = new TGraphAsymmErrors(10,xxFopi,yyFopi,exFopi,exFopi,eyFopi,eyFopi); | |
254 | fopi->SetTitle(""); | |
255 | fopi->SetMarkerColor(cFOPI); | |
256 | fopi->SetMarkerStyle(kFullTriangleDown); | |
257 | fopi->SetMarkerSize(1.2); | |
258 | fopi->SetLineWidth(2); | |
259 | fopi->SetLineColor(cFOPI); | |
260 | fopi->Draw("p"); | |
261 | ||
262 | ||
263 | //EOS -------------------------------------------------------------------- | |
264 | //float xxhq[4]={4.75,17.,130.,200.}; | |
265 | //float yyhq[4]={9.1.,10.,9.8,8.4}; | |
266 | float xxhq[4]={2.35}; | |
267 | float yyhq[4]={-0.065}; | |
268 | float exhq[4]={0.}; | |
269 | float eyhq[4]={0.007}; | |
270 | eos = new TGraphAsymmErrors(1,xxhq,yyhq,exhq,exhq,eyhq,eyhq); | |
271 | eos->SetTitle(""); | |
272 | eos->SetMarkerColor(kRed); | |
273 | eos->SetMarkerStyle(3); | |
274 | eos->SetMarkerSize(1.2); | |
275 | eos->SetLineWidth(2); | |
276 | eos->SetLineColor(kRed); | |
277 | eos->Draw("p"); | |
278 | ||
279 | ||
280 | ||
281 | // ALICE measured v2{4} + statistical error | |
282 | float v24ALICE = 0.073; | |
283 | float v24eALICE = 0.001; | |
284 | float corrpt = 0.88; | |
285 | float xxhh[1]={2760.}; // sqrt(s) | |
286 | float yyhh[1]={v24ALICE*corrpt}; | |
287 | float exhh[1]={0.}; | |
288 | ||
289 | ||
290 | // ALICE systematic error | |
291 | // hijing 15.3 % | |
292 | // therminator 6.1% | |
293 | float v24eALICE2 = 0.004; | |
294 | float eyhh2[1]={v24eALICE2}; | |
295 | ||
296 | //draw with systematic error | |
297 | //alice = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh,eyhh); | |
298 | alice = new TGraphAsymmErrors(1,xxhh,yyhh,exhh,exhh,eyhh2,eyhh2); | |
299 | alice->SetTitle(""); | |
300 | alice->SetMarkerColor(cAlice); | |
301 | alice->SetMarkerStyle(20); | |
302 | alice->SetMarkerSize(1.2); | |
303 | alice->SetLineWidth(2); | |
304 | alice->SetLineColor(cAlice); | |
305 | alice->Draw("p"); | |
306 | ||
307 | //------------------------------------------------------------------ | |
308 | TLegend *legend = new TLegend(.68,.24,.86,.66); | |
309 | // TLegend *legend = new TLegend(.45,.23,.62,.52); | |
310 | legend->Clear(); | |
311 | legend->SetBorderSize(1); | |
312 | legend->SetFillColor(kWhite); | |
313 | legend->SetTextSize(0.028); | |
314 | //legend->AddEntry(e895,"E895","P"); | |
315 | // legend->AddEntry(e877,"E877","P"); | |
316 | // legend->AddEntry(data,"NA49 pions","P"); | |
317 | // legend->AddEntry(ceres,"CERES charged particles","P"); | |
318 | // legend->AddEntry(star,"STAR charged particles","P"); | |
319 | //legend->SetHeader("A+A:"); | |
320 | //legend->AddEntry(data,"NA49 std/mod","P"); | |
321 | legend->AddEntry(alice,"ALICE","p"); | |
322 | legend->AddEntry(star,"STAR","P"); | |
323 | legend->AddEntry(phobos,"PHOBOS","P"); | |
324 | legend->AddEntry(phenix,"PHENIX","P"); | |
325 | legend->AddEntry(cumul,"NA49","P"); | |
326 | legend->AddEntry(ceres,"CERES","P"); | |
327 | legend->AddEntry(e877,"E877","P"); | |
328 | legend->AddEntry(eos,"EOS","P"); | |
329 | legend->AddEntry(e895,"E895","P"); | |
330 | legend->AddEntry(fopi,"FOPI","P"); | |
331 | legend->SetTextFont(22); // 22 = Times New Roman (bold) | |
332 | ||
333 | legend->Draw(); | |
334 | //--------------------------- insert logo ----------- | |
335 | TImage *ps = TImage::Open("macros/ALICElogo.png"); | |
336 | float xlogo=0.75; | |
337 | float ylogo=0.76; | |
338 | float wlogo=0.08; | |
339 | TPad *aliceLogo = new TPad("aliceLogo", "aliceLogo",xlogo,ylogo,xlogo+wlogo,ylogo+wlogo); | |
340 | // aliceLogo->Draw(); | |
341 | aliceLogo->cd(); | |
342 | // ps->Draw("same"); | |
343 | //-------------------------- print ----------- | |
344 | canvas->Print("v2edep.png"); | |
345 | canvas->Print("v2edep.pdf"); | |
346 | canvas->Print("v2edep.eps"); | |
347 | ||
348 | ||
349 | } | |
350 |