Fix some documentation issues
[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / scripts / OtherPNchData.C
CommitLineData
281a2bf8 1#include <TGraphAsymmErrors.h>
2
671df6c9 3/**
4 * Class that holds information on other @f$ P(N_{ch})@f$ results
5 *
6 * @ingroup pwglf_forward_multdist
7 */
281a2bf8 8struct OtherPNch
9{
10 /* ===================================================================
11 *
12 * CMS data
13 */
14 /**
c8b1a7db 15 * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le0.5@f$
281a2bf8 16 * http://hepdata.cedar.ac.uk/view/p8068/d2
17 *
18 * @return p8068_d2x1y1
19 */
20 TGraphAsymmErrors* CMSsqrts900eta05()
21 {
6b079439 22 //printf("CMS sqrt(s)=900GeV |eta|<05");
281a2bf8 23 double x[] = { 0.0, 1.0,
24 2.0, 3.0,
25 4.0, 5.0,
26 6.0, 7.0,
27 8.0, 9.0,
28 10.0, 11.0,
29 12.0, 13.0,
30 14.0, 15.0,
31 16.0, 17.0,
32 18.0, 19.0,
33 20.0, 22.0,
34 26.0 };
35 double xem[] = { 0.5, 0.5,
36 0.5, 0.5,
37 0.5, 0.5,
38 0.5, 0.5,
39 0.5, 0.5,
40 0.5, 0.5,
41 0.5, 0.5,
42 0.5, 0.5,
43 0.5, 0.5,
44 0.5, 0.5,
45 0.5, 1.5,
46 2.5 };
47 double xep[] = { 0.5, 0.5,
48 0.5, 0.5,
49 0.5, 0.5,
50 0.5, 0.5,
51 0.5, 0.5,
52 0.5, 0.5,
53 0.5, 0.5,
54 0.5, 0.5,
55 0.5, 0.5,
56 0.5, 0.5,
57 0.5, 1.5,
58 2.5 };
59 double y[] = { 0.193, 0.1504,
60 0.13507, 0.11647,
61 0.09566, 0.07567,
62 0.05817, 0.04427,
63 0.0337, 0.02572,
64 0.01948, 0.01457,
65 0.01069, 0.00769,
66 0.0055, 0.00396,
67 0.00289, 0.002112,
68 0.001544, 0.001118,
69 7.52E-4, 3.35E-4,
70 9.8E-5 };
71 double yem[] = { 0.02189155088156159, 0.012180722474467595,
72 0.007909551188278638, 0.006072832946821442,
73 0.004838605170914444, 0.003828276374558138,
74 0.003016836091006603, 0.002365797962633327,
75 0.0018847015678881366, 0.0015610893632332517,
76 0.001308166656049603, 0.0011140017953306899,
77 9.219544457292888E-4, 7.203471385380802E-4,
78 5.508175741568165E-4, 4.3011626335213136E-4,
79 3.3837848631377264E-4, 2.746961958236772E-4,
80 2.2377220560203628E-4, 1.9070920271449933E-4,
81 1.6011246047700347E-4, 8.134494452638099E-5,
82 2.5059928172283337E-5 };
83 double yep[] = { 0.031663227883461285, 0.015861904047118684,
84 0.010948538715280683, 0.008490583018850945,
85 0.007046453008429135, 0.005675746646917919,
86 0.004395201929377079, 0.003365070578754627,
87 0.0026049952015310893, 0.0020535335400231475,
88 0.0016522106403240478, 0.0013497407158413798,
89 0.001118033988749895, 9.060905032059435E-4,
90 7.061161377563892E-4, 5.457105459856901E-4,
91 4.148493702538308E-4, 3.257759966602819E-4,
92 2.6689698387205504E-4, 2.292705825002414E-4,
93 1.904205871222962E-4, 9.646242791885347E-5,
94 3.417601498127012E-5 };
95 int np = 23;
96 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
97 x, y,
98 xem,
99 xep,
100 yem,
101 yep);
102 g->SetName("/HepData/8068/d2x1y1");
103 g->SetTitle("/HepData/8068/d2x1y1");
104
105 return g;
106 }
107 /**
c8b1a7db 108 * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.0@f$
281a2bf8 109 * http://hepdata.cedar.ac.uk/view/p8068/d3
110 *
111 * @return p8068_d3x1y1
112 */
113 TGraphAsymmErrors* CMSsqrts900eta10()
114 {
6b079439 115 //printf("CMS sqrt(s)=900GeV |eta|<10");
281a2bf8 116
117 double x[] = { 0.0, 1.0,
118 2.0, 3.0,
119 4.0, 5.0,
120 6.0, 7.0,
121 8.0, 9.0,
122 10.0, 11.0,
123 12.0, 13.0,
124 14.0, 15.0,
125 16.0, 17.0,
126 18.0, 19.0,
127 20.0, 21.0,
128 22.0, 23.0,
129 24.0, 25.0,
130 26.0, 27.0,
131 28.0, 29.0,
132 30.0, 31.0,
133 32.0, 33.0,
134 34.0, 36.0,
135 40.0, 45.0,
136 50.0, 55.0 };
137 double xem[] = { 0.5, 0.5,
138 0.5, 0.5,
139 0.5, 0.5,
140 0.5, 0.5,
141 0.5, 0.5,
142 0.5, 0.5,
143 0.5, 0.5,
144 0.5, 0.5,
145 0.5, 0.5,
146 0.5, 0.5,
147 0.5, 0.5,
148 0.5, 0.5,
149 0.5, 0.5,
150 0.5, 0.5,
151 0.5, 0.5,
152 0.5, 0.5,
153 0.5, 0.5,
154 0.5, 1.5,
155 2.5, 2.5,
156 2.5, 2.5 };
157 double xep[] = { 0.5, 0.5,
158 0.5, 0.5,
159 0.5, 0.5,
160 0.5, 0.5,
161 0.5, 0.5,
162 0.5, 0.5,
163 0.5, 0.5,
164 0.5, 0.5,
165 0.5, 0.5,
166 0.5, 0.5,
167 0.5, 0.5,
168 0.5, 0.5,
169 0.5, 0.5,
170 0.5, 0.5,
171 0.5, 0.5,
172 0.5, 0.5,
173 0.5, 0.5,
174 0.5, 1.5,
175 2.5, 2.5,
176 2.5, 2.5 };
177 double y[] = { 0.1076, 0.0655,
178 0.07574, 0.08168,
179 0.0813, 0.07622,
180 0.06843, 0.06034,
181 0.05259, 0.0458,
182 0.03987, 0.03472,
183 0.03035, 0.02661,
184 0.02328, 0.02018,
185 0.01728, 0.0147,
186 0.0125, 0.01067,
187 0.00916, 0.00786,
188 0.00671, 0.0057,
189 0.00478, 0.00394,
190 0.00319, 0.00251,
191 0.00196, 0.00156,
192 0.001301, 0.001121,
193 9.84E-4, 8.25E-4,
194 6.39E-4, 3.78E-4,
195 1.96E-4, 6.0E-5,
196 1.25E-5, 1.6E-6 };
197 double yem[] = { 0.015802847844613326, 0.009464142856064674,
198 0.008214195030555337, 0.00676426640516176,
199 0.0058916636020736966, 0.005200240379059414,
200 0.004441733895676326, 0.003671511950137164,
201 0.0029149442533262966, 0.0024150569351466646,
202 0.0020554318281081475, 0.001767059704707229,
203 0.0015500322577288513, 0.0013800362314084368,
204 0.0012684636376341263, 0.0011763077828527702,
205 0.001084158659975559, 9.82344135219425E-4,
206 8.805679985100526E-4, 7.692203845452875E-4,
207 6.705221845696084E-4, 5.913543776789007E-4,
208 5.375872022286245E-4, 5.154609587543949E-4,
209 5.124451190127582E-4, 4.904079934095691E-4,
210 4.4922154890432404E-4, 3.7E-4,
211 3.0083217912982643E-4, 2.4166091947189146E-4,
212 2.1108529081866412E-4, 1.8252944967867517E-4,
213 1.662077013859466E-4, 1.4663219291819924E-4,
214 1.2229881438509532E-4, 7.826237921249263E-5,
215 4.74236228055175E-5, 2.1095023109728988E-5,
216 6.161980201201558E-6, 1.3601470508735443E-6 };
217 double yep[] = { 0.02725949375905576, 0.011154371340420759,
218 0.009740395269186974, 0.008248078564126314,
219 0.007285581651453781, 0.006454339625399333,
220 0.005595158621522718, 0.004804039966528172,
221 0.00405504623894722, 0.0034544174617437305,
222 0.0029644561052577585, 0.002545584412271571,
223 0.0021980445855350615, 0.0019087430418995638,
224 0.0016690416411821486, 0.001468911161370898,
225 0.0013179529581893276, 0.0011866338946785568,
226 0.001065129100156408, 9.339164844888434E-4,
227 8.052328855678958E-4, 6.868041933477111E-4,
228 5.94810894318522E-4, 5.440588203494177E-4,
229 5.220153254455275E-4, 5.192301994298867E-4,
230 4.972926703662542E-4, 4.368065933568311E-4,
231 3.5735136770411276E-4, 2.8792360097775935E-4,
232 2.4200206610688264E-4, 2.0310588371585893E-4,
233 1.670748335327616E-4, 1.364587849865299E-4,
234 1.1970797801316335E-4, 9.470480452437458E-5,
235 5.758472019555187E-5, 2.195449840010015E-5,
236 7.481310045707236E-6, 1.7029386365926402E-6 };
237 int np = 40;
238 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
239 x,
240 y,
241 xem,
242 xep,
243 yem,
244 yep);
245 g->SetName("/HepData/8068/d3x1y1");
246 g->SetTitle("/HepData/8068/d3x1y1");
247 return g;
248 }
249 /**
c8b1a7db 250 * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.5@f$
281a2bf8 251 * http://hepdata.cedar.ac.uk/view/p8068/d4
252 *
253 * @return p8068_d4x1y1
254 */
255 TGraphAsymmErrors* CMSsqrts900eta15()
256 {
6b079439 257 //printf("CMS sqrt(s)=900GeV |eta|<15");
258 double x[] = { 0.0, 1.0, // 2
259 2.0, 3.0, // 4
260 4.0, 5.0, // 6
261 6.0, 7.0, // 8
262 8.0, 9.0, // 10
263 10.0, 11.0, // 12
264 12.0, 13.0, // 14
265 14.0, 15.0, // 16
266 16.0, 17.0, // 18
267 18.0, 19.0, // 20
268 20.0, 21.0, // 22
269 22.0, 23.0, // 24
270 24.0, 25.0, // 26
271 26.0, 27.0, // 28
272 28.0, 29.0, // 30
273 30.0, 31.0, // 32
274 32.0, 33.0, // 34
275 34.0, 35.0, // 36
276 36.0, 37.0, // 38
277 38.0, 39.0, // 40
278 40.0, 41.0, // 42
279 42.0, 43.0, // 44
280 44.0, 45.0, // 46
281 46.5, 49.0, // 48
282 53.0, 58.0, // 50
283 63.0, 68.0 }; // 52
284 double xem[] = { 0.5, 0.5, // 2
285 0.5, 0.5, // 4
286 0.5, 0.5, // 6
287 0.5, 0.5, // 8
288 0.5, 0.5, // 10
289 0.5, 0.5, // 12
290 0.5, 0.5, // 14
291 0.5, 0.5, // 16
292 0.5, 0.5, // 18
293 0.5, 0.5, // 20
294 0.5, 0.5, // 22
295 0.5, 0.5, // 24
296 0.5, 0.5, // 26
297 0.5, 0.5, // 28
298 0.5, 0.5, // 30
299 0.5, 0.5, // 32
300 0.5, 0.5, // 34
301 0.5, 0.5, // 36
302 0.5, 0.5, // 38
303 0.5, 0.5, // 40
304 0.5, 0.5, // 42
305 0.5, 0.5, // 44
306 0.5, 0.5, // 46
307 1.0, 1.5, // 48
308 2.5, 2.5, // 50
309 2.5, 2.5 }; // 52
281a2bf8 310 double xep[] = { 0.5, 0.5,
311 0.5, 0.5,
312 0.5, 0.5,
313 0.5, 0.5,
314 0.5, 0.5,
315 0.5, 0.5,
316 0.5, 0.5,
317 0.5, 0.5,
318 0.5, 0.5,
319 0.5, 0.5,
320 0.5, 0.5,
321 0.5, 0.5,
322 0.5, 0.5,
323 0.5, 0.5,
324 0.5, 0.5,
325 0.5, 0.5,
326 0.5, 0.5,
327 0.5, 0.5,
328 0.5, 0.5,
329 0.5, 0.5,
330 0.5, 0.5,
331 0.5, 0.5,
332 0.5, 0.5,
333 1.0, 1.5,
334 2.5, 2.5,
335 2.5, 2.5 };
336 double y[] = { 0.0753, 0.03557,
337 0.04347, 0.05068,
338 0.05566, 0.05786,
339 0.05835, 0.05688,
340 0.05348, 0.0489,
341 0.04422, 0.04002,
342 0.03647, 0.03337,
343 0.03055, 0.02797,
344 0.02556, 0.02329,
345 0.02122, 0.01936,
346 0.01765, 0.01604,
347 0.01451, 0.01304,
348 0.01168, 0.01045,
349 0.00937, 0.0084,
350 0.00751, 0.00667,
351 0.00592, 0.00526,
352 0.00466, 0.00413,
353 0.00362, 0.00316,
354 0.00272, 0.00233,
355 0.00203, 0.00178,
356 0.00156, 0.00136,
357 0.0012, 0.001053,
358 9.13E-4, 7.84E-4,
359 5.85E-4, 3.73E-4,
360 2.13E-4, 9.6E-5,
361 3.6E-5, 7.4E-6 };
362 double yem[] = { 0.014524806367039803, 0.0061507154055443014,
363 0.006900384047283165, 0.0068944397886992964,
364 0.00601252858621063, 0.004686747272896204,
365 0.0042880881520789655, 0.0039679339712248235,
366 0.0035288099977187778, 0.0031197756329582422,
367 0.0027127108213003464, 0.0023229722340140013,
368 0.002033937068839643, 0.0017866449003649271,
369 0.001591257364476281, 0.0014504137340772805,
370 0.0013412307780542466, 0.0012490796611905903,
371 0.0011594826432508594, 0.0010893117092916976,
372 0.0010384603988597735, 9.972462083156796E-4,
373 9.588013350011566E-4, 8.984430978086481E-4,
374 8.160882305241266E-4, 7.178439941937245E-4,
375 6.390618123468182E-4, 5.88727441181401E-4,
376 5.571355310873647E-4, 5.255473337388365E-4,
377 4.846648326421054E-4, 4.6615448083226656E-4,
378 4.4384682042344296E-4, 4.3081318457076036E-4,
379 4.0853396431630995E-4, 3.992492955535426E-4,
380 3.676955262170047E-4, 3.2695565448543625E-4,
381 2.9068883707497264E-4, 2.505992817228334E-4,
382 2.1954498400100151E-4, 1.8601075237738274E-4,
383 1.8601075237738274E-4, 1.7042300314218148E-4,
384 1.7998055450520203E-4, 1.6058953888718904E-4,
385 1.4985659811966905E-4, 8.561541917201597E-5,
386 6.324555320336759E-5, 4.657252408878007E-5,
387 2.4166091947189146E-5, 5.590169943749474E-6 };
388 double yep[] = { 0.03405304685340212, 0.00724989655098609,
389 0.007817985674072318, 0.007751728839426725,
390 0.006919566460407762, 0.00571192611997039,
391 0.005232991496266739, 0.004813003220443552,
392 0.00435332057170156, 0.0038838769290491172,
393 0.003445692963686695, 0.0030549959083442323,
394 0.0027252339349127445, 0.0024367396249907374,
395 0.002189794510907359, 0.001969568480657629,
396 0.0017906702655709678, 0.0016395731151735808,
397 0.0015002999700059986, 0.0013710215169719256,
398 0.0012712198865656563, 0.0011812705024675763,
399 0.0011035397591387453, 0.0010139033484509261,
400 9.121403400793104E-4, 8.132035415564789E-4,
401 7.242237223399962E-4, 6.545991139621257E-4,
402 6.135144660071188E-4, 5.724508712544685E-4,
403 5.692099788303082E-4, 5.126402247190518E-4,
404 4.531004303683677E-4, 4.4011362169330773E-4,
405 4.1785164831552356E-4, 4.0853396431630995E-4,
406 3.8626415831655934E-4, 3.54682957019364E-4,
407 3.1780497164141406E-4, 2.7730849247724094E-4,
408 2.459674775249769E-4, 2.3706539182259396E-4,
409 2.1954498400100151E-4, 2.0068881383873889E-4,
410 1.5864425612041553E-4, 1.6058953888718904E-4,
411 1.4338758663147937E-4, 1.0245974819410792E-4,
412 6.609841147864296E-5, 4.657252408878007E-5,
413 2.5079872407968906E-5, 5.941380311005179E-6 };
414 int np = 52;
415 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
416 x,
417 y,
418 xem,
419 xep,
420 yem,
421 yep);
422 g->SetName("/HepData/8068/d4x1y1");
423 g->SetTitle("/HepData/8068/d4x1y1");
424 return g;
425 }
281a2bf8 426 /**
c8b1a7db 427 * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le2.0@f$
281a2bf8 428 * http://hepdata.cedar.ac.uk/view/p8068/d6
429 *
430 * @return p8068_d6x1y1
431 */
432 TGraphAsymmErrors* CMSsqrts900eta20()
433 {
6b079439 434 //printf("CMS sqrt(s)=900GeV |eta|<20");
281a2bf8 435 double x[] = { 0.0, 1.0,
436 2.0, 3.0,
437 4.0, 5.0,
438 6.0, 7.0,
439 8.0, 9.0,
440 10.0, 11.0,
441 12.0, 13.0,
442 14.0, 15.0,
443 16.0, 17.0,
444 18.0, 19.0,
445 20.0, 21.0,
446 22.0, 23.0,
447 24.0, 25.0,
448 26.0, 27.0,
449 28.0, 29.0,
450 30.0, 31.0,
451 32.0, 33.0,
452 34.0, 35.0,
453 36.0, 37.0,
454 38.0, 39.0,
455 40.0, 41.0,
456 42.0, 43.0,
457 44.0, 45.0,
458 46.0, 47.0,
459 48.0, 49.0,
460 50.0, 51.0,
461 52.0, 53.0,
462 54.0, 55.0,
463 56.0, 57.0,
464 58.0, 59.0,
465 60.0, 61.0,
466 62.5, 64.5,
467 67.5, 74.5,
468 84.5, 94.5 };
469 double xem[] = { 0.5, 0.5,
470 0.5, 0.5,
471 0.5, 0.5,
472 0.5, 0.5,
473 0.5, 0.5,
474 0.5, 0.5,
475 0.5, 0.5,
476 0.5, 0.5,
477 0.5, 0.5,
478 0.5, 0.5,
479 0.5, 0.5,
480 0.5, 0.5,
481 0.5, 0.5,
482 0.5, 0.5,
483 0.5, 0.5,
484 0.5, 0.5,
485 0.5, 0.5,
486 0.5, 0.5,
487 0.5, 0.5,
488 0.5, 0.5,
489 0.5, 0.5,
490 0.5, 0.5,
491 0.5, 0.5,
492 0.5, 0.5,
493 0.5, 0.5,
494 0.5, 0.5,
495 0.5, 0.5,
496 0.5, 0.5,
497 0.5, 0.5,
498 0.5, 0.5,
499 0.5, 0.5,
500 1.0, 1.0,
501 2.0, 5.0,
502 5.0, 5.0 };
503 double xep[] = { 0.5, 0.5,
504 0.5, 0.5,
505 0.5, 0.5,
506 0.5, 0.5,
507 0.5, 0.5,
508 0.5, 0.5,
509 0.5, 0.5,
510 0.5, 0.5,
511 0.5, 0.5,
512 0.5, 0.5,
513 0.5, 0.5,
514 0.5, 0.5,
515 0.5, 0.5,
516 0.5, 0.5,
517 0.5, 0.5,
518 0.5, 0.5,
519 0.5, 0.5,
520 0.5, 0.5,
521 0.5, 0.5,
522 0.5, 0.5,
523 0.5, 0.5,
524 0.5, 0.5,
525 0.5, 0.5,
526 0.5, 0.5,
527 0.5, 0.5,
528 0.5, 0.5,
529 0.5, 0.5,
530 0.5, 0.5,
531 0.5, 0.5,
532 0.5, 0.5,
533 0.5, 0.5,
534 1.0, 1.0,
535 2.0, 5.0,
536 5.0, 5.0 };
537 double y[] = { 0.0494, 0.01789,
538 0.01895, 0.02174,
539 0.0255, 0.02961,
540 0.0335, 0.03646,
541 0.03826, 0.03906,
542 0.0392, 0.03868,
543 0.03753, 0.03596,
544 0.0341, 0.03206,
545 0.03, 0.02802,
546 0.02616, 0.02448,
547 0.023, 0.0217,
548 0.02058, 0.01959,
549 0.01866, 0.01776,
550 0.01685, 0.01595,
551 0.01504, 0.01415,
552 0.01326, 0.01241,
553 0.01161, 0.01087,
554 0.01015, 0.00948,
555 0.00887, 0.00829,
556 0.00777, 0.00734,
557 0.00681, 0.00624,
558 0.00578, 0.00535,
559 0.00494, 0.00457,
560 0.00422, 0.00388,
561 0.00356, 0.00326,
562 0.00298, 0.00274,
563 0.00253, 0.00234,
564 0.00212, 0.00189,
565 0.00167, 0.0015,
566 0.00136, 0.00123,
567 0.00112, 0.001041,
568 9.17E-4, 7.64E-4,
569 4.92E-4, 2.45E-4,
570 9.9E-5, 1.64E-5 };
571 double yem[] = { 0.031045933711196384, 0.0025678006153126453,
572 0.003808017857100988, 0.004928792955683978,
573 0.005475372133471842, 0.005444713031923721,
574 0.004928792955683978, 0.0041622710147226115,
575 0.003405994715204356, 0.0030985964564621835,
576 0.0029811742652854096, 0.002811547616527239,
577 0.0026404734423962684, 0.002492308167141455,
578 0.0023525518060183073, 0.0021654098919142305,
579 0.001967053634245899, 0.0017886866690396059,
580 0.0016107451691686056, 0.0014552319402761885,
581 0.001326989073052224, 0.001225275479229059,
582 0.001164817582284883, 0.0011168258592994702,
583 0.0010756393447619887, 0.001027861858422619,
584 9.7082439194738E-4, 9.202716990106781E-4,
585 8.825531145489205E-4, 8.509406559801923E-4,
586 8.228000972289684E-4, 8.099382692526635E-4,
587 7.912016177940992E-4, 7.725283166331187E-4,
588 7.316419889536139E-4, 6.815423684555495E-4,
589 6.407807737440318E-4, 5.950630218724736E-4,
590 5.636488268416782E-4, 5.636488268416782E-4,
591 5.322593352868506E-4, 4.919349550499538E-4,
592 4.7853944456021594E-4, 4.7423622805517506E-4,
593 4.651881339845203E-4, 4.428317965096906E-4,
594 4.248529157249601E-4, 3.9824615503479755E-4,
595 3.671511950137164E-4, 3.3615472627943223E-4,
596 3.138470965295043E-4, 3.138470965295043E-4,
597 3.176476034853718E-4, 3.3541019662496844E-4,
598 3.3105890714493693E-4, 2.996664812754339E-4,
599 2.641968962724581E-4, 2.3853720883753126E-4,
600 2.1633307652783938E-4, 1.9416487838947602E-4,
601 1.6401219466856724E-4, 1.5181897114655994E-4,
602 1.5966527487215246E-4, 1.3961733416735901E-4,
603 9.740636529508736E-5, 6.800735254367721E-5,
604 3.733630940518893E-5, 9.437160589923222E-6 };
605 double yep[] = { 0.027627884464793896, 0.0038148263394288343,
606 0.004395600072800072, 0.005347317084295637,
607 0.005874325493194942, 0.0058935897380119695,
608 0.005446999173857106, 0.004729587719875803,
609 0.004031935019317648, 0.0037237615390892046,
610 0.0035759474269066098, 0.00339607125955861,
611 0.0032050585018061684, 0.003026549190084311,
612 0.0028468403537957655, 0.00266865134478073,
613 0.002469412885687608, 0.0022703523955544874,
614 0.0020813697413001857, 0.0019043371550227129,
615 0.00175524927004685, 0.001643471934655411,
616 0.0015533190271158081, 0.001475127113166862,
617 0.0014046351839534703, 0.0013267252918370102,
618 0.0012490796611905903, 0.0011691449867317568,
619 0.0011208925015361642, 0.0011472575996697516,
620 9.932774033471214E-4, 9.52102935611481E-4,
621 9.139474820797965E-4, 8.75956619930462E-4,
622 8.34865258589672E-4, 7.749193506423748E-4,
623 7.244998274671983E-4, 6.685057965343307E-4,
624 6.092618484691126E-4, 5.727128425310541E-4,
625 5.412947441089743E-4, 5.280151512977634E-4,
626 5.147815070493501E-4, 5.015974481593781E-4,
627 4.924428900898052E-4, 4.792702786528704E-4,
628 4.7010637094172637E-4, 4.428317965096906E-4,
629 4.204759208325728E-4, 3.8910152916687437E-4,
630 3.6674241641784496E-4, 3.488552708502481E-4,
631 3.2649655434629013E-4, 3.0886890422961E-4,
632 3.04138126514911E-4, 2.996664812754339E-4,
633 2.9068883707497264E-4, 2.641968962724581E-4,
634 2.418677324489565E-4, 2.1095023109728988E-4,
635 1.972308292331602E-4, 1.739482681718907E-4,
636 1.4475496537252186E-4, 1.2083045973594572E-4,
637 1.097679370308106E-4, 6.992138442565336E-5,
638 3.54682957019364E-5, 9.078546139112804E-6 };
639 int np = 68;
640 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
641 x,
642 y,
643 xem,
644 xep,
645 yem,
646 yep);
647 g->SetName("/HepData/8068/d6x1y1");
648 g->SetTitle("/HepData/8068/d6x1y1");
649 return g;
650 }
6b079439 651 /**
c8b1a7db 652 * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le2.4@f$
6b079439 653 * http://hepdata.cedar.ac.uk/view/p8068/d5
654 *
655 * @return p8068_d5x1y1
656 */
657 TGraphAsymmErrors* CMSsqrts900eta24()
658 {
659 //printf("CMS sqrt(s)=900GeV |eta|<24");
660 double x[] = { 0.0, 1.0,
661 2.0, 3.0,
662 4.0, 5.0,
663 6.0, 7.0,
664 8.0, 9.0,
665 10.0, 11.0,
666 12.0, 13.0,
667 14.0, 15.0,
668 16.0, 17.0,
669 18.0, 19.0,
670 20.0, 21.0,
671 22.0, 23.0,
672 24.0, 25.0,
673 26.0, 27.0,
674 28.0, 29.0,
675 30.0, 31.0,
676 32.0, 33.0,
677 34.0, 35.0,
678 36.0, 37.0,
679 38.0, 39.0,
680 40.0, 41.0,
681 42.0, 43.0,
682 44.0, 45.0,
683 46.0, 47.0,
684 48.0, 49.0,
685 50.0, 51.0,
686 52.0, 53.0,
687 54.0, 55.0,
688 56.0, 57.5,
689 60.0, 66.5,
690 76.5, 86.5 };
691 double xem[] = { 0.5, 0.5,
692 0.5, 0.5,
693 0.5, 0.5,
694 0.5, 0.5,
695 0.5, 0.5,
696 0.5, 0.5,
697 0.5, 0.5,
698 0.5, 0.5,
699 0.5, 0.5,
700 0.5, 0.5,
701 0.5, 0.5,
702 0.5, 0.5,
703 0.5, 0.5,
704 0.5, 0.5,
705 0.5, 0.5,
706 0.5, 0.5,
707 0.5, 0.5,
708 0.5, 0.5,
709 0.5, 0.5,
710 0.5, 0.5,
711 0.5, 0.5,
712 0.5, 0.5,
713 0.5, 0.5,
714 0.5, 0.5,
715 0.5, 0.5,
716 0.5, 0.5,
717 0.5, 0.5,
718 0.5, 0.5,
719 0.5, 1.0,
720 1.5, 5.0,
721 5.0, 5.0 };
722 double xep[] = { 0.5, 0.5,
723 0.5, 0.5,
724 0.5, 0.5,
725 0.5, 0.5,
726 0.5, 0.5,
727 0.5, 0.5,
728 0.5, 0.5,
729 0.5, 0.5,
730 0.5, 0.5,
731 0.5, 0.5,
732 0.5, 0.5,
733 0.5, 0.5,
734 0.5, 0.5,
735 0.5, 0.5,
736 0.5, 0.5,
737 0.5, 0.5,
738 0.5, 0.5,
739 0.5, 0.5,
740 0.5, 0.5,
741 0.5, 0.5,
742 0.5, 0.5,
743 0.5, 0.5,
744 0.5, 0.5,
745 0.5, 0.5,
746 0.5, 0.5,
747 0.5, 0.5,
748 0.5, 0.5,
749 0.5, 0.5,
750 0.5, 1.0,
751 1.5, 5.0,
752 5.0, 5.0 };
753 double y[] = { 0.0561, 0.02321,
754 0.02645, 0.031,
755 0.0362, 0.04081,
756 0.04402, 0.04548,
757 0.04592, 0.04515,
758 0.04342, 0.04101,
759 0.03836, 0.03573,
760 0.03316, 0.03079,
761 0.02868, 0.02679,
762 0.02506, 0.02343,
763 0.02187, 0.02039,
764 0.01901, 0.0178,
765 0.01668, 0.01556,
766 0.01445, 0.01334,
767 0.01229, 0.01132,
768 0.01045, 0.00967,
769 0.00896, 0.00829,
770 0.00763, 0.00699,
771 0.00637, 0.00578,
772 0.00527, 0.00482,
773 0.00441, 0.00412,
774 0.00363, 0.00328,
775 0.00292, 0.00261,
776 0.00233, 0.00209,
777 0.00185, 0.00164,
778 0.00146, 0.00131,
779 0.00117, 0.00106,
780 9.7E-4, 8.9E-4,
781 8.2E-4, 6.63E-4,
782 4.53E-4, 2.12E-4,
783 7.7E-5, 1.12E-5 };
784 double yem[] = { 0.01767766952966369, 0.003660737630587584,
785 0.005079094407470687, 0.006069810540700591,
786 0.006147235150862541, 0.005458580401532985,
787 0.004473712105176193, 0.0035970126494078384,
788 0.003456848275524976, 0.003228761372415125,
789 0.0029797483115189447, 0.002739361239413305,
790 0.0025401181074902798, 0.002313114783144148,
791 0.0020573040611440983, 0.001820109886792553,
792 0.0016324827717314506, 0.001484318025222358,
793 0.0013727709204379296, 0.001287866452703851,
794 0.001217538500417954, 0.0011279184367674819,
795 0.0010412012293500234, 9.741663102366044E-4,
796 9.360555539069249E-4, 9.139474820797965E-4,
797 9.013878188659973E-4, 8.697700845625814E-4,
798 8.190848551890091E-4, 7.433034373659253E-4,
799 6.741661516273269E-4, 6.239390995922599E-4,
800 5.961543424315552E-4, 5.8309518948453E-4,
801 5.738466694161429E-4, 5.515432893255072E-4,
802 5.292447448959697E-4, 4.976946855251721E-4,
803 4.609772228646444E-4, 4.295346318982906E-4,
804 4.0249223594996216E-4, 4.386342439892262E-4,
805 3.712142238654117E-4, 3.6674241641784496E-4,
806 3.6249137920783716E-4, 3.5341194094144583E-4,
807 3.3105890714493693E-4, 3.130495168499705E-4,
808 2.8178005607210744E-4, 2.5553864678361276E-4,
809 2.3323807579381201E-4, 2.1633307652783938E-4,
810 1.8601075237738274E-4, 1.6401219466856724E-4,
811 1.562049935181331E-4, 1.562049935181331E-4,
812 1.5480633061990717E-4, 1.395886814895821E-4,
813 1.0756393447619885E-4, 6.676076692189808E-5,
814 2.9966648127543395E-5, 6.551335741663679E-6 };
815 double yep[] = { 0.09323352401362935, 0.004710169848317574,
816 0.005736898116578331, 0.006618164398078972,
817 0.006715772479767314, 0.006096630216767292,
818 0.00519042387479096, 0.004362247127341595,
819 0.004172217635742412, 0.003913693907295255,
820 0.003644283194264683, 0.0033638668225719043,
821 0.0031145144083789367, 0.0028765604460883488,
822 0.002629087294100369, 0.0023903974564912843,
823 0.002181604913819182, 0.0020026232796010335,
824 0.0018514858897652987, 0.0017280335644888382,
825 0.00160822883943797, 0.0014787156589419076,
826 0.0013612494260788505, 0.0012539936203984452,
827 0.0011763077828527702, 0.0011157060544785082,
828 0.001084158659975559, 0.001033247308247159,
829 9.630160954002794E-4, 9.244457799135652E-4,
830 8.35224520712844E-4, 6.989277502002622E-4,
831 6.519202405202648E-4, 6.203224967708329E-4,
832 6.10982814815605E-4, 5.88727441181401E-4,
833 5.664803615307418E-4, 5.348831648126533E-4,
834 4.976946855251721E-4, 4.477722635447622E-4,
835 4.204759208325728E-4, 3.9357337308308855E-4,
836 3.981205847478877E-4, 3.8483762809787716E-4,
837 3.807886552931954E-4, 3.716180835212409E-4,
838 3.584689665786984E-4, 3.4928498393145964E-4,
839 3.1780497164141406E-4, 2.9068883707497264E-4,
840 2.505992817228334E-4, 2.2472205054244234E-4,
841 2.0248456731316588E-4, 1.8027756377319944E-4,
842 1.6401219466856724E-4, 1.6401219466856724E-4,
843 1.5640012787718557E-4, 1.4223923509355636E-4,
844 1.0288342918079667E-4, 6.86804193347711E-5,
845 3.361547262794322E-5, 7.071067811865475E-6 };
846 int np = 62;
847 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
848 x,
849 y,
850 xem,
851 xep,
852 yem,
853 yep);
854 g->SetName("/HepData/8068/d5x1y1");
855 g->SetTitle("/HepData/8068/d5x1y1");
856 return g;
857 }
281a2bf8 858 /**
859 * Get CMS data
860 *
861 * @param eta Limits on eta
862 * @param sNN Center of mass energy
863 *
864 * @return Graph or null
865 */
866 TGraphAsymmErrors* GetCMS(Double_t eta, UShort_t sNN)
867 {
868 Double_t aEta = TMath::Abs(eta);
869 TGraphAsymmErrors* g = 0;
6b079439 870 //printf("CMS %dGeV |eta|<%f ->", sNN, aEta);
281a2bf8 871 switch (sNN) {
872 case 900:
873 if (aEta <= 0.51) g = CMSsqrts900eta05();
874 else if (aEta <= 1.01) g = CMSsqrts900eta10();
875 else if (aEta <= 1.51) g = CMSsqrts900eta15();
876 else if (aEta <= 2.01) g = CMSsqrts900eta20();
877 else if (aEta <= 2.41) g = CMSsqrts900eta24();
878 break;
879 default:
880 return g;
881 }
882 return g;
883 }
884
885 /* ===================================================================
886 *
887 * ALICE data
888 */
889 /**
c8b1a7db 890 * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le0.5@f$
281a2bf8 891 * http://hepdata.cedar.ac.uk/view/p7742/d8
892 *
893 * @return p7742_d8x1y1
894 */
895 TGraphAsymmErrors* ALICEsqrts900eta05()
896 {
6b079439 897 //printf("ALICE sqrt(s)=900GeV |eta|<05");
281a2bf8 898 double x[] = { 0.0, 1.0,
899 2.0, 3.0,
900 4.0, 5.0,
901 6.0, 7.0,
902 8.0, 9.0,
903 10.0, 11.0,
904 12.0, 13.0,
905 14.0, 15.0,
906 16.0, 17.0,
907 18.0, 19.0,
908 20.0, 21.5,
909 23.5, 25.5 };
910 double xem[] = { 0.5, 0.5,
911 0.5, 0.5,
912 0.5, 0.5,
913 0.5, 0.5,
914 0.5, 0.5,
915 0.5, 0.5,
916 0.5, 0.5,
917 0.5, 0.5,
918 0.5, 0.5,
919 0.5, 0.5,
920 0.5, 1.0,
921 1.0, 1.0 };
922 double xep[] = { 0.5, 0.5,
923 0.5, 0.5,
924 0.5, 0.5,
925 0.5, 0.5,
926 0.5, 0.5,
927 0.5, 0.5,
928 0.5, 0.5,
929 0.5, 0.5,
930 0.5, 0.5,
931 0.5, 0.5,
932 0.5, 1.0,
933 1.0, 1.0 };
934 double y[] = { 0.179767, 0.155985,
935 0.140806, 0.120933,
936 0.097273, 0.075258,
937 0.057395, 0.043794,
938 0.033212, 0.024939,
939 0.018642, 0.013927,
940 0.010337, 0.007693,
941 0.005665, 0.004134,
942 0.002994, 0.002144,
943 0.001527, 0.001093,
944 8.03E-4, 4.82E-4,
945 2.56E-4, 1.02E-4 };
946 double yem[] = { 0.028627623478032542, 0.006826681770816624,
947 0.004736912179891031, 0.005286866936097409,
948 0.004605827830043151, 0.003753897441326814,
949 0.003172958871463669, 0.0023132784095305087,
950 0.0017387768689512751, 0.0012885437516824954,
951 9.891637882575362E-4, 7.64476945368531E-4,
952 5.769662035162892E-4, 4.500277769204919E-4,
953 3.7107950630558943E-4, 3.1458703088334716E-4,
954 2.7004073766748605E-4, 2.210701246211256E-4,
955 1.7819090885900997E-4, 1.3576450198781714E-4,
956 1.2447088012864694E-4, 9.762171889492625E-5,
957 6.095900261651268E-5, 3.4828149534536E-5 };
958 double yep[] = { 0.028627623478032542, 0.006826681770816624,
959 0.004736912179891031, 0.005286866936097409,
960 0.004605827830043151, 0.003753897441326814,
961 0.003172958871463669, 0.0023132784095305087,
962 0.0017387768689512751, 0.0012885437516824954,
963 9.891637882575362E-4, 7.64476945368531E-4,
964 5.769662035162892E-4, 4.500277769204919E-4,
965 3.7107950630558943E-4, 3.1458703088334716E-4,
966 2.7004073766748605E-4, 2.210701246211256E-4,
967 1.7819090885900997E-4, 1.3576450198781714E-4,
968 1.2447088012864694E-4, 9.762171889492625E-5,
969 6.095900261651268E-5, 3.4828149534536E-5 };
970 int np = 24;
971 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
972 x,
973 y,
974 xem,
975 xep,
976 yem,
977 yep);
978 g->SetName("/HepData/7742/d8x1y1");
979 g->SetTitle("/HepData/7742/d8x1y1");
980 return g;
981 }
982 /**
c8b1a7db 983 * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.0@f$
281a2bf8 984 * http://hepdata.cedar.ac.uk/view/p7742/d9
985 *
986 * @return p7742_d9x1y1
987 */
988 TGraphAsymmErrors* ALICEsqrts900eta10()
989 {
6b079439 990 //printf("ALICE sqrt(s)=900GeV |eta|<10");
281a2bf8 991 double x[] = { 0.0, 1.0,
992 2.0, 3.0,
993 4.0, 5.0,
994 6.0, 7.0,
995 8.0, 9.0,
996 10.0, 11.0,
997 12.0, 13.0,
998 14.0, 15.0,
999 16.0, 17.0,
1000 18.0, 19.0,
1001 20.0, 21.0,
1002 22.0, 23.0,
1003 24.0, 25.0,
1004 26.0, 27.0,
1005 28.0, 29.0,
1006 30.0, 31.0,
1007 32.0, 33.0,
1008 34.0, 35.0,
1009 36.0, 37.0,
1010 38.0, 39.0,
1011 40.0, 41.5 };
1012 double xem[] = { 0.5, 0.5,
1013 0.5, 0.5,
1014 0.5, 0.5,
1015 0.5, 0.5,
1016 0.5, 0.5,
1017 0.5, 0.5,
1018 0.5, 0.5,
1019 0.5, 0.5,
1020 0.5, 0.5,
1021 0.5, 0.5,
1022 0.5, 0.5,
1023 0.5, 0.5,
1024 0.5, 0.5,
1025 0.5, 0.5,
1026 0.5, 0.5,
1027 0.5, 0.5,
1028 0.5, 0.5,
1029 0.5, 0.5,
1030 0.5, 0.5,
1031 0.5, 0.5,
1032 0.5, 1.0 };
1033 double xep[] = { 0.5, 0.5,
1034 0.5, 0.5,
1035 0.5, 0.5,
1036 0.5, 0.5,
1037 0.5, 0.5,
1038 0.5, 0.5,
1039 0.5, 0.5,
1040 0.5, 0.5,
1041 0.5, 0.5,
1042 0.5, 0.5,
1043 0.5, 0.5,
1044 0.5, 0.5,
1045 0.5, 0.5,
1046 0.5, 0.5,
1047 0.5, 0.5,
1048 0.5, 0.5,
1049 0.5, 0.5,
1050 0.5, 0.5,
1051 0.5, 0.5,
1052 0.5, 0.5,
1053 0.5, 1.0 };
1054 double y[] = { 0.092396, 0.076087,
1055 0.073233, 0.078506,
1056 0.080885, 0.077084,
1057 0.070058, 0.062258,
1058 0.054437, 0.047195,
1059 0.040978, 0.035352,
1060 0.03048, 0.026336,
1061 0.022788, 0.019725,
1062 0.017071, 0.014674,
1063 0.012509, 0.010593,
1064 0.008935, 0.007527,
1065 0.006361, 0.005416,
1066 0.004665, 0.004041,
1067 0.003519, 0.00306,
1068 0.002637, 0.002234,
1069 0.001856, 0.001512,
1070 0.001213, 9.62E-4,
1071 7.61E-4, 6.02E-4,
1072 4.84E-4, 3.93E-4,
1073 3.27E-4, 2.76E-4,
1074 2.37E-4, 1.7E-4 };
1075 double yem[] = { 0.021888152137629163, 0.0066076761421849355,
1076 0.004451948000594796, 0.0034909631908686747,
1077 0.002952765483407038, 0.0028889466938661224,
1078 0.0029957812002881653, 0.0032637882897026274,
1079 0.002922475833946279, 0.0026227742945209753,
1080 0.0023487871338203465, 0.002091141554271255,
1081 0.0017347841364273539, 0.0014762689456870655,
1082 0.0012594490859101847, 0.0010506721658062519,
1083 9.019456746389995E-4, 7.836938177630343E-4,
1084 6.432239112470866E-4, 5.679656679765071E-4,
1085 4.967293025381128E-4, 4.475187146924696E-4,
1086 4.040519768544636E-4, 3.3646693745448454E-4,
1087 2.9206163733020465E-4, 2.820017730440715E-4,
1088 2.6982401672201087E-4, 2.5470178640912594E-4,
1089 2.3519353732617736E-4, 2.1914607000811127E-4,
1090 1.9727392123643713E-4, 1.7402586014727813E-4,
1091 1.5060212481900778E-4, 1.3009611831257687E-4,
1092 1.1269871339105873E-4, 9.693812459502196E-5,
1093 8.48528137423857E-5, 7.427651041883969E-5,
1094 6.726812023536856E-5, 6.030754513325841E-5,
1095 5.544366510251645E-5, 4.601086828130936E-5 };
1096 double yep[] = { 0.021888152137629163, 0.0066076761421849355,
1097 0.004451948000594796, 0.0034909631908686747,
1098 0.002952765483407038, 0.0028889466938661224,
1099 0.0029957812002881653, 0.0032637882897026274,
1100 0.002922475833946279, 0.0026227742945209753,
1101 0.0023487871338203465, 0.002091141554271255,
1102 0.0017347841364273539, 0.0014762689456870655,
1103 0.0012594490859101847, 0.0010506721658062519,
1104 9.019456746389995E-4, 7.836938177630343E-4,
1105 6.432239112470866E-4, 5.679656679765071E-4,
1106 4.967293025381128E-4, 4.475187146924696E-4,
1107 4.040519768544636E-4, 3.3646693745448454E-4,
1108 2.9206163733020465E-4, 2.820017730440715E-4,
1109 2.6982401672201087E-4, 2.5470178640912594E-4,
1110 2.3519353732617736E-4, 2.1914607000811127E-4,
1111 1.9727392123643713E-4, 1.7402586014727813E-4,
1112 1.5060212481900778E-4, 1.3009611831257687E-4,
1113 1.1269871339105873E-4, 9.693812459502196E-5,
1114 8.48528137423857E-5, 7.427651041883969E-5,
1115 6.726812023536856E-5, 6.030754513325841E-5,
1116 5.544366510251645E-5, 4.601086828130936E-5 };
1117 int np = 42;
1118 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
1119 x,
1120 y,
1121 xem,
1122 xep,
1123 yem,
1124 yep);
1125 g->SetName("/HepData/7742/d9x1y1");
1126 g->SetTitle("/HepData/7742/d9x1y1");
1127 return g;
1128 }
1129
1130 /**
c8b1a7db 1131 * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.3@f$
281a2bf8 1132 * http://hepdata.cedar.ac.uk/view/p7742/d10
1133 *
1134 * @return p7742_d10x1y1
1135 */
1136 TGraphAsymmErrors* ALICEsqrts900eta13()
1137 {
6b079439 1138 //printf("ALICE sqrt(s)=900GeV |eta|<13");
281a2bf8 1139 double x[] = { 0.0, 1.0,
1140 2.0, 3.0,
1141 4.0, 5.0,
1142 6.0, 7.0,
1143 8.0, 9.0,
1144 10.0, 11.0,
1145 12.0, 13.0,
1146 14.0, 15.0,
1147 16.0, 17.0,
1148 18.0, 19.0,
1149 20.0, 21.0,
1150 22.0, 23.0,
1151 24.0, 25.0,
1152 26.0, 27.0,
1153 28.0, 29.0,
1154 30.0, 31.0,
1155 32.0, 33.0,
1156 34.0, 35.0,
1157 36.0, 37.0,
1158 38.0, 39.0,
1159 40.0, 41.5,
1160 43.5, 45.5,
1161 47.5, 49.5,
1162 51.5, 53.5 };
1163 double xem[] = { 0.5, 0.5,
1164 0.5, 0.5,
1165 0.5, 0.5,
1166 0.5, 0.5,
1167 0.5, 0.5,
1168 0.5, 0.5,
1169 0.5, 0.5,
1170 0.5, 0.5,
1171 0.5, 0.5,
1172 0.5, 0.5,
1173 0.5, 0.5,
1174 0.5, 0.5,
1175 0.5, 0.5,
1176 0.5, 0.5,
1177 0.5, 0.5,
1178 0.5, 0.5,
1179 0.5, 0.5,
1180 0.5, 0.5,
1181 0.5, 0.5,
1182 0.5, 0.5,
1183 0.5, 1.0,
1184 1.0, 1.0,
1185 1.0, 1.0,
1186 1.0, 1.0 };
1187 double xep[] = { 0.5, 0.5,
1188 0.5, 0.5,
1189 0.5, 0.5,
1190 0.5, 0.5,
1191 0.5, 0.5,
1192 0.5, 0.5,
1193 0.5, 0.5,
1194 0.5, 0.5,
1195 0.5, 0.5,
1196 0.5, 0.5,
1197 0.5, 0.5,
1198 0.5, 0.5,
1199 0.5, 0.5,
1200 0.5, 0.5,
1201 0.5, 0.5,
1202 0.5, 0.5,
1203 0.5, 0.5,
1204 0.5, 0.5,
1205 0.5, 0.5,
1206 0.5, 0.5,
1207 0.5, 1.0,
1208 1.0, 1.0,
1209 1.0, 1.0,
1210 1.0, 1.0 };
1211 double y[] = { 0.068695, 0.056382,
1212 0.049828, 0.056116,
1213 0.063345, 0.066198,
1214 0.065598, 0.061117,
1215 0.055042, 0.04927,
1216 0.044132, 0.039615,
1217 0.035906, 0.03262,
1218 0.029489, 0.026657,
1219 0.023839, 0.021302,
1220 0.018884, 0.016732,
1221 0.014707, 0.012944,
1222 0.01137, 0.009982,
1223 0.008753, 0.007661,
1224 0.006711, 0.005891,
1225 0.005179, 0.004572,
1226 0.004064, 0.003636,
1227 0.003278, 0.002968,
1228 0.002677, 0.002392,
1229 0.002111, 0.001829,
1230 0.001561, 0.001319,
1231 0.001115, 7.93E-4,
1232 5.24E-4, 3.51E-4,
1233 2.2E-4, 1.56E-4,
1234 1.27E-4, 8.5E-5 };
1235 double yem[] = { 0.021322725646595934, 0.007630831737104416,
1236 0.004038016592338372, 0.0036354708085748677,
1237 0.0035003511252444373, 0.0030936373413831173,
1238 0.0026374775828431228, 0.002920906195001818,
1239 0.00276560481631053, 0.002864668392676542,
1240 0.0026843833556330957, 0.0022856071403458645,
1241 0.002325611317481922, 0.002006865217198205,
1242 0.0017203360136903488, 0.0014989629748596194,
1243 0.0012394615766533467, 0.0011637718848640398,
1244 9.151568171630477E-4, 9.078199160626517E-4,
1245 7.98723982361867E-4, 6.967302203866285E-4,
1246 6.308502199413107E-4, 5.846990678973245E-4,
1247 5.4516236113657E-4, 5.460009157501477E-4,
1248 4.930608481719068E-4, 4.3796118549478783E-4,
1249 3.4394330928221294E-4, 3.2323366161339073E-4,
1250 3.0869078379504626E-4, 2.9443165590676553E-4,
1251 2.756320010448714E-4, 2.630304164920856E-4,
1252 2.5628889948649746E-4, 2.4142286552851616E-4,
1253 2.245551157288562E-4, 2.0678491240900531E-4,
1254 1.8828170383762728E-4, 1.7255723688098393E-4,
1255 1.5911316727411343E-4, 1.3023056476879765E-4,
1256 1.0761505470890214E-4, 8.62670273047588E-5,
1257 6.129437168288782E-5, 6.037383539249432E-5,
1258 3.535533905932738E-5, 3.8078865529319545E-5 };
1259 double yep[] = { 0.021322725646595934, 0.007630831737104416,
1260 0.004038016592338372, 0.0036354708085748677,
1261 0.0035003511252444373, 0.0030936373413831173,
1262 0.0026374775828431228, 0.002920906195001818,
1263 0.00276560481631053, 0.002864668392676542,
1264 0.0026843833556330957, 0.0022856071403458645,
1265 0.002325611317481922, 0.002006865217198205,
1266 0.0017203360136903488, 0.0014989629748596194,
1267 0.0012394615766533467, 0.0011637718848640398,
1268 9.151568171630477E-4, 9.078199160626517E-4,
1269 7.98723982361867E-4, 6.967302203866285E-4,
1270 6.308502199413107E-4, 5.846990678973245E-4,
1271 5.4516236113657E-4, 5.460009157501477E-4,
1272 4.930608481719068E-4, 4.3796118549478783E-4,
1273 3.4394330928221294E-4, 3.2323366161339073E-4,
1274 3.0869078379504626E-4, 2.9443165590676553E-4,
1275 2.756320010448714E-4, 2.630304164920856E-4,
1276 2.5628889948649746E-4, 2.4142286552851616E-4,
1277 2.245551157288562E-4, 2.0678491240900531E-4,
1278 1.8828170383762728E-4, 1.7255723688098393E-4,
1279 1.5911316727411343E-4, 1.3023056476879765E-4,
1280 1.0761505470890214E-4, 8.62670273047588E-5,
1281 6.129437168288782E-5, 6.037383539249432E-5,
1282 3.535533905932738E-5, 3.8078865529319545E-5 };
1283 int np = 48;
1284 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,
1285 x,
1286 y,
1287 xem,
1288 xep,
1289 yem,
1290 yep);
1291 g->SetName("/HepData/7742/d10x1y1");
1292 g->SetTitle("/HepData/7742/d10x1y1");
1293 return g;
1294 }
1295 /**
1296 * Get ALICE data
1297 *
1298 * @param eta Eta limit
1299 * @param sNN Center of mass energy
1300 *
1301 * @return
1302 */
1303 TGraphAsymmErrors* GetALICE(Double_t eta, UShort_t sNN)
1304 {
1305 Double_t aEta = TMath::Abs(eta);
1306 TGraphAsymmErrors* g = 0;
6b079439 1307 //printf("CMS %dGeV |eta|<%f ->", sNN, aEta);
281a2bf8 1308 switch (sNN) {
1309 case 900:
1310 if (aEta <= 0.51) g = ALICEsqrts900eta05();
1311 else if (aEta <= 1.01) g = ALICEsqrts900eta10();
1312 else if (aEta <= 1.31) g = ALICEsqrts900eta13();
1313 break;
1314 }
1315 return g;
1316 }
1317
1318 /* =================================================================
1319 *
1320 * Measurements from other sources, such as published ALICE, CMS, ...
1321 */
1322 static TGraphAsymmErrors* GetData(UShort_t type, Double_t eta, UShort_t sNN)
1323 {
1324 OtherPNch d;
1325 TGraphAsymmErrors* g = 0;
6b079439 1326 //printf("type=%d (%s), eta=%f, sqrt(s)=%dGeV ->",
1327 // type, (type==0 ? "CMS" : "ALICE"), eta, sNN);
281a2bf8 1328 switch (type) {
1329 case 0: g = d.GetCMS(eta, sNN); break;
1330 case 1: g = d.GetALICE(eta, sNN); break;
1331 }
1332 if (!g) {
1333 // Warning("GetOther", "No other data found for type=%d, eta=%f, sNN=%d",
1334 // type, eta, sNN);
1335 return g;
1336 }
6b079439 1337 //printf("%p (%s)\n", g, (g ? g->GetName() : ""));
281a2bf8 1338 return g;
1339 }
1340};
1341
1342//
1343// EOF
1344//