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