PWG2/FORWARD/analysis2/OtherData.C

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