]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG2/FORWARD/analysis2/OtherData.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Minor fixes
[u/mrichter/AliRoot.git] / PWG2 / FORWARD / analysis2 / OtherData.C
CommitLineData
7e4038b5 1#include <TGraphAsymmErrors.h>
2#include <TMultiGraph.h>
3#include <TStyle.h>
4#include <TMath.h>
5#include <TCanvas.h>
6#include <TLegend.h>
7
8//____________________________________________________________________
9/**
10 * @defgroup pwg2_forward_analysis_otherdata External data
11 *
12 * @ingroup pwg2_forward_analysis_scripts
13 */
14//____________________________________________________________________
15/**
16 * Values used
17 *
18 * @ingroup pwg2_forward_analysis_otherdata
19 */
20enum {
21 /** Color used for UA5 data */
22 UA5Color = kBlue+1,
23 /** Color used for CMS data */
24 CMSColor = kGreen+1,
25 /** Color used for ALICE data */
26 ALICEColor = kMagenta+1,
27 /** Marker style INEL data */
28 INELStyle = 20,
29 /** Marker style INEL>0 data */
30 INELGtStyle= 22,
31 /** Marker style NSD data */
32 NSDStyle = 21,
33 /** Colour offset for mirror data */
34 MirrorOff = 4
35};
36
37//____________________________________________________________________
38/**
39 * Set graph attributes
40 *
41 * @param g Graph
42 * @param marker Marker style
43 * @param color Marker and line color
44 * @param name Name of graph
45 * @param title Title of graph
46 *
47 * @ingroup pwg2_forward_analysis_otherdata
48 */
49void
50SetGraphAttributes(TGraph* g, Int_t marker, Int_t color,
51 const Char_t* name, const Char_t* title)
52{
53 g->SetName(name);
54 g->SetTitle(title);
55 g->SetMarkerStyle(marker);
56 g->SetMarkerColor(color);
57 g->SetLineColor(color);
58 g->SetFillColor(0);
59 g->SetFillStyle(0);
60 g->GetHistogram()->SetStats(kFALSE);
61 g->GetHistogram()->SetXTitle("#eta");
62 g->GetHistogram()->SetYTitle("#frac{1}{N} #frac{dN_{ch}}{#eta}");
63}
64
65//____________________________________________________________________
66/**
67 * Get the UA5 NSD data for pp at @f$ \sqrt{s} = 900GeV@f$
68 * p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
69 *
70 * @param mirrored Wether to produce the mirrored plot
71 *
72 * @return graph of data
73 *
74 * @ingroup pwg2_forward_analysis_otherdata
75 */
76TGraphAsymmErrors* UA5Nsd(Bool_t mirrored=false)
77{
78 //UA5 data NSD - p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
79 double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
80 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
81 4.625 };
82 double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
83 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
84 double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
85 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
86 double y[] = { 3.48, 3.38, 3.52, 3.68, 3.71, 3.86, 3.76, 3.66, 3.72,
87 3.69, 3.56, 3.41, 3.15, 3.09, 2.74, 2.73, 2.32, 1.99, 1.69 };
88 double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
89 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
90 double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
91 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
92 const int np = 20;
93 double xm[np];
94 for (int i = 0; i < np; i++) xm[i] = -x[i];
95
96 TGraphAsymmErrors* g = new TGraphAsymmErrors(19,x, y,exm,exp,eym,eyp);
97 TGraphAsymmErrors* gm = new TGraphAsymmErrors(19,xm,y,exm,exp,eym,eyp);
98 SetGraphAttributes(g, NSDStyle, UA5Color,"ua5_nsd", "UA5 NSD");
99 SetGraphAttributes(gm, NSDStyle+MirrorOff, UA5Color,"ua5_nsd_mirrored",
100 "UA5 NSD (mirrored)");
101
102 return (mirrored ? gm : g);
103}
104
105//____________________________________________________________________
106/**
107 * Get the UA5 INEL data for pp at @f$ \sqrt{s} = 900GeV@f$
108 * p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
109 *
110 * @param mirrored Wether to produce the mirrored plot
111 *
112 * @return graph of data
113 *
114 * @ingroup pwg2_forward_analysis_otherdata
115 */
116TGraphAsymmErrors* UA5Inel(Bool_t mirrored=false)
117{
118 //UA5 data INEL - p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
119 double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
120 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
121 4.625 };
122 double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
123 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
124 double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
125 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
126 double y[] = { 3.14, 3.04, 3.17, 3.33, 3.33, 3.53, 3.46, 3.41, 3.45,
127 3.39, 3.07, 3.07, 2.93, 2.93, 2.55, 2.48, 2.18, 1.91, 1.52 };
128 double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
129 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
130 double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
131 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
132 const int np = 19;
133 double xm[np];
134 for (int i = 0; i < np; i++) xm[i] = -x[i];
135 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,x, y,exm,exp,eym,eyp);
136 TGraphAsymmErrors* gm = new TGraphAsymmErrors(np,xm,y,exm,exp,eym,eyp);
137
138 SetGraphAttributes(g, INELStyle, UA5Color, "ua5_inel", "UA5 INEL");
139 SetGraphAttributes(gm, INELStyle+MirrorOff, UA5Color, "ua5_inel_mirrored",
140 "UA5 INEL (mirrored)");
141
142 return (mirrored ? gm : g);
143}
144
145//____________________________________________________________________
146/**
147 * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$ \sqrt{s}
148 * = 900GeV@f$
149 * p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
150 *
151 * @return graph of data
152 *
153 * @ingroup pwg2_forward_analysis_otherdata
154 */
155TGraphAsymmErrors* AliceCentralInel900()
156{
157 // INEL - p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
158 TGraphAsymmErrors* g = 0;
159 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3,
160 -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3 };
161 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
162 0.1, 0.1, 0.1, 0.1, 0.1 };
163 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
164 0.1, 0.1, 0.1, 0.1, 0.1 };
165 double y[] = { 3.28, 3.28, 3.22, 3.12, 3.06, 3.02, 2.98, 3.02, 3.02,
166 3.05, 3.15, 3.21, 3.26, 3.33 };
167 double eym[] = { 0.06324555320336758, 0.06324555320336758,
168 0.06324555320336758, 0.06324555320336758,
169 0.06324555320336758, 0.05385164807134505,
170 0.05385164807134505, 0.05385164807134505,
171 0.05385164807134505, 0.06324555320336758,
172 0.06324555320336758, 0.06324555320336758,
173 0.06324555320336758, 0.06324555320336758 };
174 double eyp[] = { 0.08246211251235322, 0.08246211251235322,
175 0.08246211251235322, 0.08246211251235322,
176 0.08246211251235322, 0.08246211251235322,
177 0.07280109889280519, 0.08246211251235322,
178 0.08246211251235322, 0.08246211251235322,
179 0.08246211251235322, 0.08246211251235322,
180 0.08246211251235322, 0.08246211251235322 };
181 const int np = 14;
182 for (int i = 0; i < np; i++) {
183 eym[i] += 0.02;
184 eyp[i] += 0.02;
185 }
186 g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
187 SetGraphAttributes(g, INELStyle, ALICEColor, "alice_inel",
188 "ALICE INEL (publ.)");
189
190 return g;
191}
192
193//____________________________________________________________________
194/**
195 * Get the ALICE INEL>0 data in @f$ |\eta|<1.3@f$ for pp at @f$
196 * \sqrt{s} = 900GeV@f$
197 *
198 * p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
199 *
200 * @return graph of data
201 *
202 * @ingroup pwg2_forward_analysis_otherdata
203 */
204TGraphAsymmErrors* AliceCentralInelGt900()
205{
206 // INEL>0 - p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
207 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7,
208 0.9 };
209 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
210 0.1 };
211 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
212 0.1 };
213 double y[] = { 4.0, 3.87, 3.8, 3.7, 3.67, 3.73, 3.72, 3.77, 3.92,
214 4.01 };
215 double eym[] = { 0.07615773105863909, 0.07615773105863909,
216 0.07615773105863909, 0.06324555320336758,
217 0.06324555320336758, 0.06324555320336758,
218 0.0670820393249937, 0.07615773105863909,
219 0.07615773105863909, 0.07615773105863909 };
220 double eyp[] = { 0.08544003745317531, 0.07615773105863909,
221 0.07615773105863909, 0.07280109889280519,
222 0.07280109889280519, 0.07280109889280519,
223 0.07615773105863909, 0.07615773105863909,
224 0.08544003745317531, 0.08544003745317531 };
225 const int np = 10;
226 for (int i = 0; i < np; i++) {
227 double stat = (i >= 3 && i<=5) ? 0.02 : 0.03;
228 eym[i] += stat;
229 eyp[i] += stat;
230 }
231
232 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
233 SetGraphAttributes(g, INELGtStyle, ALICEColor, "alice_inelgt900",
234 "ALICE INEL>0 (publ.)");
235 return g;
236
237}
238
239//____________________________________________________________________
240/**
241 * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
242 * \sqrt{s} = 2.36TeV@f$
243 *
244 * p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
245 *
246 * @return graph of data
247 *
248 * @ingroup pwg2_forward_analysis_otherdata
249 */
250TGraphAsymmErrors* AliceCentralInelGt2360()
251{
252 // INEL>0 - p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
253 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
254 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
255 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
256 double y[] = { 4.91, 4.76, 4.63, 4.64, 4.55, 4.55, 4.64, 4.66, 4.82, 4.88 };
257 double eym[] = { 0.08544003745317531, 0.08544003745317531,
258 0.08544003745317531, 0.08544003745317531,
259 0.08544003745317531, 0.08544003745317531,
260 0.08544003745317531, 0.08544003745317531,
261 0.08544003745317531, 0.08544003745317531 };
262 double eyp[] = { 0.11401754250991379, 0.11401754250991379,
263 0.1044030650891055, 0.1044030650891055,
264 0.1044030650891055, 0.1044030650891055,
265 0.1044030650891055, 0.1044030650891055,
266 0.11401754250991379, 0.11401754250991379 };
267 const int np = 10;
268 for (int i = 0; i < np; i++) {
269 double stat = 0.3;
270 eym[i] += stat;
271 eyp[i] += stat;
272 }
273
274 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
275 SetGraphAttributes(g, INELGtStyle, ALICEColor, "alice_inelgt2360",
276 "ALICE INEL>0 (publ.)");
277 return g;
278}
279
280//____________________________________________________________________
281/**
282 * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
283 * \sqrt{s} = 7TeV@f$
284 *
285 * p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
286 *
287 * @return graph of data
288 *
289 * @ingroup pwg2_forward_analysis_otherdata
290 */
291TGraphAsymmErrors* AliceCentralInelGt7000()
292{
293 // INEL>0 - p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
294// Plot: p7741_d6x1y1
295 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
296 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
297 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
298 double y[] = { 6.22, 6.07, 6.01, 5.84, 5.85, 5.85, 5.91, 6.01, 6.17, 6.26 };
299 double eym[] = { 0.1216552506059644, 0.1216552506059644,
300 0.1216552506059644, 0.11180339887498948,
301 0.11180339887498948, 0.11180339887498948,
302 0.11180339887498948, 0.1216552506059644,
303 0.1216552506059644, 0.1216552506059644 };
304 double eyp[] = { 0.21095023109728983, 0.21095023109728983,
305 0.2009975124224178, 0.2009975124224178,
306 0.2009975124224178, 0.2009975124224178,
307 0.2009975124224178, 0.2009975124224178,
308 0.21095023109728983, 0.21095023109728983 };
309 const int np = 10;
310 for (int i = 0; i < np; i++) {
311 double stat = 0.2;
312 eym[i] += stat;
313 eyp[i] += stat;
314 }
315
316 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
317 SetGraphAttributes(g, INELGtStyle, ALICEColor, "alice_inelgt7000",
318 "ALICE INEL>0 (publ.)");
319 return g;
320}
321
322//____________________________________________________________________
323/**
324 * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
325 * \sqrt{s} = 900GeV@f$
326 *
327 * p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
328 *
329 * @return graph of data
330 *
331 * @ingroup pwg2_forward_analysis_otherdata
332 */
333TGraphAsymmErrors* AliceCentralNsd900()
334{
335 //NSD - p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
336 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
337 0.5, 0.7, 0.9, 1.1, 1.3 };
338 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
339 0.1, 0.1, 0.1, 0.1, 0.1 };
340 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
341 0.1, 0.1, 0.1, 0.1, 0.1 };
342 double y[] = { 3.9, 3.89, 3.81, 3.7, 3.64, 3.59, 3.53, 3.58, 3.59,
343 3.61, 3.74, 3.8, 3.87, 3.95 };
344 double eym[] = { 0.13341664064126335, 0.13152946437965907,
345 0.13152946437965907, 0.1216552506059644,
346 0.1216552506059644, 0.1216552506059644,
347 0.1216552506059644, 0.1216552506059644,
348 0.1216552506059644, 0.1216552506059644,
349 0.1216552506059644, 0.13152946437965907,
350 0.13152946437965907, 0.13341664064126335 };
351 double eyp[] = { 0.13341664064126335, 0.13152946437965907,
352 0.13152946437965907, 0.1216552506059644,
353 0.1216552506059644, 0.1216552506059644,
354 0.1216552506059644, 0.1216552506059644,
355 0.1216552506059644, 0.1216552506059644,
356 0.1216552506059644, 0.13152946437965907,
357 0.13152946437965907, 0.13341664064126335 };
358 const int np = 14;
359 for (int i = 0; i < np; i++) {
360 double stat = (i == 0 || i == np-1) ? 0.03 : 0.02;
361 eym[i] += stat;
362 eyp[i] += stat;
363 }
364
365 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
366 SetGraphAttributes(g, NSDStyle, ALICEColor, "alice_nsd", "ALICE NSD (publ.)");
367
368 return g;
369}
370
371//____________________________________________________________________
372/**
373 * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$
374 * \sqrt{s} = 2.36TeV@f$
375 *
376 * p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
377 *
378 * @return graph of data
379 *
380 * @ingroup pwg2_forward_analysis_otherdata
381 */
382TGraphAsymmErrors* AliceCentralInel2360()
383{
384 // INEL - p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
385 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
386 0.5, 0.7, 0.9, 1.1, 1.3 };
387 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
388 0.1, 0.1, 0.1, 0.1, 0.1 };
389 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
390 0.1, 0.1, 0.1, 0.1, 0.1 };
391 double y[] = { 4.08, 4.01, 4.0, 3.88, 3.77, 3.8, 3.73, 3.71, 3.79,
392 3.82, 3.94, 4.02, 4.05, 4.16 };
393 double eym[] = { 0.13341664064126335, 0.12369316876852982,
394 0.12369316876852982, 0.1216552506059644,
395 0.1216552506059644, 0.1216552506059644,
396 0.1216552506059644, 0.11180339887498948,
397 0.1216552506059644, 0.1216552506059644,
398 0.12369316876852982, 0.12369316876852982,
399 0.13341664064126335, 0.13341664064126335 };
400 double eyp[] = { 0.2716615541441225, 0.2716615541441225,
401 0.2716615541441225, 0.260768096208106,
402 0.25079872407968906, 0.25079872407968906,
403 0.25079872407968906, 0.25079872407968906,
404 0.25079872407968906, 0.260768096208106,
405 0.261725046566048, 0.2716615541441225,
406 0.2716615541441225, 0.2816025568065745 };
407 const int np = 14;
408 for (int i = 0; i < np; i++) {
409 double stat = (i < 3 || i > np-1-4) ? 0.03 : 0.02;
410 eym[i] += stat;
411 eyp[i] += stat;
412 }
413
414 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
415 SetGraphAttributes(g, NSDStyle, ALICEColor, "alice_inel2360",
416 "ALICE INEL (publ.)");
417 return g;
418}
419
420//____________________________________________________________________
421/**
422 * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
423 * \sqrt{s} = 2.36TeV@f$
424 *
425 * p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
426 *
427 * @return graph of data
428 *
429 * @ingroup pwg2_forward_analysis_otherdata
430 */
431TGraphAsymmErrors* AliceCentralNsd2360()
432{
433 // NSD - p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
434 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
435 0.5, 0.7, 0.9, 1.1, 1.3 };
436 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
437 0.1, 0.1, 0.1, 0.1, 0.1 };
438 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
439 0.1, 0.1, 0.1, 0.1, 0.1 };
440 double y[] = { 4.79, 4.72, 4.7, 4.56, 4.44, 4.47, 4.39, 4.37, 4.47,
441 4.5, 4.64, 4.73, 4.76, 4.9 };
442 double eym[] = { 0.13601470508735444, 0.13341664064126335,
443 0.13341664064126335, 0.12369316876852982,
444 0.12369316876852982, 0.12369316876852982,
445 0.12369316876852982, 0.12369316876852982,
446 0.12369316876852982, 0.12369316876852982,
447 0.12369316876852982, 0.13341664064126335,
448 0.13341664064126335, 0.13341664064126335 };
449 double eyp[] = { 0.18439088914585774, 0.18248287590894657,
450 0.18248287590894657, 0.1726267650163207,
451 0.1726267650163207, 0.1726267650163207,
452 0.16278820596099708, 0.16278820596099708,
453 0.1726267650163207, 0.1726267650163207,
454 0.1726267650163207, 0.18248287590894657,
455 0.18248287590894657, 0.18248287590894657 };
456 const int np = 14;
457
458 for (int i = 0; i < np; i++) {
459 double stat = (i < 1) ? 0.03 : 0.02;
460 eym[i] += stat;
461 eyp[i] += stat;
462 }
463
464 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
465 SetGraphAttributes(g, NSDStyle, ALICEColor, "alice_nsd2360",
466 "ALICE NSD (publ.)");
467 return g;
468}
469
470
471//____________________________________________________________________
472/**
473 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
474 * \sqrt{s} = 900GeV@f$
475 *
476 * p7743_d8x1y1
477 *
478 * @return graph of data
479 *
480 * @ingroup pwg2_forward_analysis_otherdata
481 */
482TGraphAsymmErrors* CMSNsd900()
483{
484 // CMS published NSD data - p7743_d8x1y1
485 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75, 1.25, 1.75,
486 2.25 };
487 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
488 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
489 double y[] = { 3.6, 3.73, 3.62, 3.54, 3.48, 3.48, 3.54, 3.62, 3.73, 3.6 };
490 double eym[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14,0.13 };
491 double eyp[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14, 0.13 };
492 const int np = 10;
493 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
494 SetGraphAttributes(g, NSDStyle, CMSColor, "cms_nsd900", "CMS NSD");
495
496 return g;
497}
498
499
500//____________________________________________________________________
501/**
502 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
503 * \sqrt{s} = 2.36GeV@f$
504 *
505 * p7743_d8x1y2
506 *
507 * @return graph of data
508 *
509 * @ingroup pwg2_forward_analysis_otherdata
510 */
511TGraphAsymmErrors* CMSNsd2360()
512{
513 // CMS NSD 2360 - p7743_d8x1y2
514 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
515 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
516 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
517 double y[] = { 4.78, 4.81, 4.66, 4.61, 4.47, 4.47, 4.61, 4.66, 4.81, 4.78 };
518 double eym[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
519 double eyp[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
520 const int np = 10;
521 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
522 SetGraphAttributes(g, NSDStyle, CMSColor, "cms_nsd2360", "CMS NSD");
523 return g;
524}
525
526
527//____________________________________________________________________
528/**
529 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
530 * \sqrt{s} = 7TeV@f$
531 *
532 * p7838_d5x1y1
533 *
534 * @return graph of data
535 *
536 * @ingroup pwg2_forward_analysis_otherdata
537 */
538TGraphAsymmErrors* CMSNsd7000()
539{
540 // CMS NSD 7000 - Plot: p7838_d5x1y1
541 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
542 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
543 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
544 double y[] = { 6.18, 6.26, 6.14, 6.01, 5.78, 5.78, 6.01, 6.14, 6.26, 6.18 };
545 double eym[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
546 double eyp[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
547 const int np = 10;
548 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
549 SetGraphAttributes(g, NSDStyle, CMSColor, "cms_nsd7000", "CMS NSD");
550 return g;
551}
552
553//____________________________________________________________________
554/**
555 * Get a multi graph of data for a given energy and trigger type
556 *
557 * @param energy Energy in GeV (900, 2360, 7000)
558 * @param type Bit pattern of trigger type
559 * - 0x1 INEL
560 * - 0x2 INEL>0
561 * - 0x4 NSD
562 *
563 * @return A multi graph with the selected data.
564 *
565 * @ingroup pwg2_forward_analysis_otherdata
566 */
567TMultiGraph*
568GetData(Int_t energy, Int_t type)
569{
570 TMultiGraph* mp = new TMultiGraph(Form("dndeta_%dGeV_%d", energy, type),"");
571 TString tn;
572 TString en;
573 if (TMath::Abs(energy-900) < 10) {
574 en.Append(" #sqrt{s}=900GeV");
575 if (type & 0x1) {
576 tn.Append(" INEL");
577 mp->Add(UA5Inel(false));
578 mp->Add(UA5Inel(true));
579 mp->Add(AliceCentralInel900());
580 }
581 if (type & 0x4) {
582 tn.Append(" NSD");
583 mp->Add(UA5Nsd(false));
584 mp->Add(UA5Nsd(true));
585 mp->Add(AliceCentralNsd900());
586 mp->Add(CMSNsd900());
587 }
588 if (type & 0x2) {
589 tn.Append(" INEL>0");
590 mp->Add(AliceCentralInelGt900());
591 }
592 }
593 if (TMath::Abs(energy-2360) < 10) {
594 en.Append(" #sqrt{s}=2.36TeV");
595 if (type & 0x1) {
596 tn.Append(" INEL");
597 mp->Add(AliceCentralInel2360());
598 }
599 if (type & 0x4) {
600 tn.Append(" NSD");
601 mp->Add(AliceCentralNsd2360());
602 mp->Add(CMSNsd2360());
603 }
604 if (type & 0x1) {
605 tn.Append(" INEL>0");
606 mp->Add(AliceCentralInelGt2360());
607 }
608 }
609 if (TMath::Abs(energy-7000) < 10) {
610 en.Append(" #sqrt{s}=7TeV");
611 if (type & 0x1) {
612 tn.Append(" INEL");
613 }
614 if (type & 0x4) {
615 tn.Append(" NSD");
616 mp->Add(CMSNsd7000());
617 }
618 if (type & 0x1) {
619 tn.Append(" INEL>0");
620 mp->Add(AliceCentralInelGt7000());
621 }
622 }
623 mp->SetTitle(Form("1/N dN_{ch}/d#eta, pp(p#bar{p}), %s, %s",
624 en.Data(), tn.Data()));
82fc512a 625 if (!mp->GetListOfGraphs() || mp->GetListOfGraphs()->GetEntries() <= 0) {
626 delete mp;
627 mp = 0;
628 }
7e4038b5 629 return mp;
630}
631
632//____________________________________________________________________
633/**
634 * Plot external data for a given selection of energy and trigger type
635 * (see GetData)
636 *
637 * @param energy Energy in GeV
638 * @param type Trigger type bit mask
639 *
640 * @ingroup pwg2_forward_analysis_otherdata
641 */
642void
643OtherData(Int_t energy=900, Int_t type=0x1)
644{
82fc512a 645 TMultiGraph* mp = GetData(energy, type);
646 if (!mp) return;
647
7e4038b5 648 gStyle->SetTitleX(0.1);
649 gStyle->SetTitleY(1.0);
650 gStyle->SetTitleW(0.85);
651 gStyle->SetTitleH(0.05);
652 gStyle->SetTitleBorderSize(0);
653 gStyle->SetTitleTextColor(kWhite);
654 gStyle->SetTitleFillColor(kBlack);
655 gStyle->SetTitleFontSize(0.02);
656
657 gStyle->SetOptTitle(1);
658 gStyle->SetOptStat(0);
659
660 TCanvas* c = new TCanvas("c", "dN/deta", 800, 600);
661 c->SetFillColor(0);
662 c->SetBorderSize(0);
663 c->SetBorderMode(0);
664 c->SetRightMargin(0.05);
665 c->SetTopMargin(0.05);
666
82fc512a 667
7e4038b5 668 mp->SetMinimum(0);
669 mp->Draw("ap");
670 if (mp->GetXaxis())
671 mp->GetXaxis()->SetTitle("#eta");
672 if (mp->GetYaxis())
673 mp->GetYaxis()->SetTitle("#frac{1}{N} #frac{dN_{ch}}{#eta}");
674
675 TLegend* l = c->BuildLegend(0.3, 0.15, 0.7, 0.5);
676 l->SetFillColor(0);
677 l->SetBorderSize(0);
678
679 c->cd();
680}
681
682//____________________________________________________________________
683//
684// EOF
685//