PWGLF/FORWARD/analysis2/scripts/OtherPNchData.C

   1 #include <TGraphAsymmErrors.h>
   2
   3 /**
   4  * Class that holds information on other @f$ P(N_{ch})@f$ results
   5  *
   6  * @ingroup pwglf_forward_multdist
   7  */
   8 struct OtherPNch
   9 {
  10   /* ===================================================================
  11    *
  12    * CMS data
  13    */
  14   /**
  15    * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le0.5@f$
  16    * http://hepdata.cedar.ac.uk/view/p8068/d2
  17    *
  18    * @return p8068_d2x1y1
  19    */
  20   TGraphAsymmErrors* CMSsqrts900eta05()
  21   {
  22     //printf("CMS sqrt(s)=900GeV |eta|<05");
  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   /**
 108    * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.0@f$
 109    * http://hepdata.cedar.ac.uk/view/p8068/d3
 110    *
 111    * @return p8068_d3x1y1
 112    */
 113   TGraphAsymmErrors* CMSsqrts900eta10()
 114   {
 115     //printf("CMS sqrt(s)=900GeV |eta|<10");
 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   /**
 250    * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.5@f$
 251    * http://hepdata.cedar.ac.uk/view/p8068/d4
 252    *
 253    * @return p8068_d4x1y1
 254    */
 255   TGraphAsymmErrors* CMSsqrts900eta15()
 256   {
 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
 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   }
 426   /**
 427    * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le2.0@f$
 428    * http://hepdata.cedar.ac.uk/view/p8068/d6
 429    *
 430    * @return p8068_d6x1y1
 431    */
 432   TGraphAsymmErrors* CMSsqrts900eta20()
 433   {
 434     //printf("CMS sqrt(s)=900GeV |eta|<20");
 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   }
 651   /**
 652    * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le2.4@f$
 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   }
 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;
 870     //printf("CMS %dGeV |eta|<%f ->", sNN, aEta);
 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   /**
 890    * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le0.5@f$
 891    * http://hepdata.cedar.ac.uk/view/p7742/d8
 892    *
 893    * @return p7742_d8x1y1
 894    */
 895   TGraphAsymmErrors* ALICEsqrts900eta05()
 896   {
 897     //printf("ALICE sqrt(s)=900GeV |eta|<05");
 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   /**
 983    * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.0@f$
 984    * http://hepdata.cedar.ac.uk/view/p7742/d9
 985    *
 986    * @return p7742_d9x1y1
 987    */
 988   TGraphAsymmErrors* ALICEsqrts900eta10()
 989   {
 990     //printf("ALICE sqrt(s)=900GeV |eta|<10");
 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   /**
1131    * ALICE @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\le1.3@f$
1132    * http://hepdata.cedar.ac.uk/view/p7742/d10
1133    *
1134    * @return p7742_d10x1y1
1135    */
1136   TGraphAsymmErrors* ALICEsqrts900eta13()
1137   {
1138     //printf("ALICE sqrt(s)=900GeV |eta|<13");
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;
1307     //printf("CMS %dGeV |eta|<%f ->", sNN, aEta);
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;
1326     //printf("type=%d (%s), eta=%f, sqrt(s)=%dGeV ->",
1327     // type, (type==0 ? "CMS" : "ALICE"), eta, sNN);
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     }
1337     //printf("%p (%s)\n", g, (g ? g->GetName() : ""));
1338     return g;
1339   }
1340 };
1341
1342 //
1343 // EOF
1344 //