/**
* @defgroup pwglf_forward_otherdata External data
*
+ * Collection of external data points for comparisons and the like
+ *
* @ingroup pwglf_forward_scripts
*/
/**
*
* @ingroup pwglf_forward_script_otherdata
*/
-#include <TGraphAsymmErrors.h>
-#include <TMultiGraph.h>
-#include <TStyle.h>
-#include <TMath.h>
-#include <TCanvas.h>
-#include <TLegend.h>
-
-//____________________________________________________________________
-/**
- * Values used
- *
- * @ingroup pwglf_forward_otherdata
- */
-enum {
- UA5,
- CMS,
- ALICE,
- PYTHIA,
- INEL,
- INELGt0,
- NSD
-};
-enum {
- /** Style used for UA5 data */
- UA5Style = 21,
- /** Style used for CMS data */
- CMSStyle = 29,
- /** Style used for ALICE published data */
- ALICEStyle = 27,
- /** Style used for Pythia data */
- PYTHIAStyle = 28,
- /** Color used for UA5 data */
- UA5Color = kBlue+1,
- /** Color used for Pytia data */
- PYTHIAColor = kGray+2,
- /** Color used for CMS data */
- CMSColor = kGreen+1,
- /** Color used for ALICE data */
- ALICEColor = kMagenta+1
-};
-enum {
- /** Marker style INEL data */
- INELStyle = 22,
- /** Marker style INEL>0 data */
- INELGt0Style= 29,
- /** Marker style NSD data */
- NSDStyle = 23,
- /** Color used for UA5 data */
- INELColor = kBlue+1,
- /** Color used for CMS data */
- INELGt0Color = kGreen+1,
- /** Color used for ALICE data */
- NSDColor = kMagenta+1
-};
-enum {
- /** Style offset for mirror data */
- MirrorOff = 4
-};
+#ifndef __CINT__
+# include <TGraphAsymmErrors.h>
+# include <TMultiGraph.h>
+# include <TStyle.h>
+# include <TMath.h>
+# include <TCanvas.h>
+# include <TLegend.h>
+# include <TLegendEntry.h>
+# include <TList.h>
+# include <TAxis.h>
+# include <TH1F.h>
+#else
+class TGraphAsymmErrors;
+class TMultiGraph;
+class TStyle;
+class TCanvas;
+class TLegend;
+class TLegendEntry;
+class TList;
+class TAxis;
+class TH1F;
+class TGraph;
+#endif
-//____________________________________________________________________
-/**
- * Set graph attributes based on trigger type and experiment.
- *
- * @param g Graph
- * @param trig Trigger (INEL, INEL>0, NSD)
- * @param exp Experiment
- * @param mirror True if mirrored data
- * @param name Name of graph
- * @param title Title of graph
- *
- * @ingroup pwglf_forward_otherdata
- */
-void
-SetGraphAttributes(TGraph* g, Int_t trig, Int_t exp, bool mirror,
- const Char_t* name, const Char_t* title)
+struct RefData
{
- Int_t color = 0;
- switch (exp) {
- case UA5: color = UA5Color; break;
- case CMS: color = CMSColor; break;
- case ALICE: color = ALICEColor; break;
- case PYTHIA: color = PYTHIAColor; break;
- }
- Int_t style = 0;
- switch (exp) {
- case UA5: style = UA5Style; break;
- case CMS: style = CMSStyle; break;
- case ALICE: style = ALICEStyle; break;
- case PYTHIA: style = PYTHIAStyle; break;
- }
- Float_t size = g->GetMarkerSize();
- switch (style) {
- case 21:
- case 25: size *= 0.8; break;
- case 27: size *= 1.4; break;
- }
+ //____________________________________________________________________
+ /**
+ * Values used
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ enum {
+ UA5,
+ CMS,
+ ALICE,
+ WIP,
+ PYTHIA,
+ INEL,
+ INELGt0,
+ NSD
+ };
+ enum {
+ /** Style used for UA5 data */
+ UA5Style = 21,
+ /** Style used for CMS data */
+ CMSStyle = 29,
+ /** Style used for ALICE published data */
+ ALICEStyle = 27,
+ /** Color used for ALICE work-in-progress data */
+ WIPStyle = 33,
+ /** Style used for Pythia data */
+ PYTHIAStyle = 28,
+ /** Color used for UA5 data */
+ UA5Color = kBlue+1,
+ /** Color used for Pytia data */
+ PYTHIAColor = kGray+2,
+ /** Color used for CMS data */
+ CMSColor = kGreen+1,
+ /** Color used for ALICE data */
+ ALICEColor = kMagenta+1,
+ /** Color used for ALICE work-in-progress data */
+ WIPColor = kMagenta+3
+ };
+ enum {
+ /** Marker style INEL data */
+ INELStyle = 22,
+ /** Marker style INEL>0 data */
+ INELGt0Style= 29,
+ /** Marker style NSD data */
+ NSDStyle = 23,
+ /** Color used for UA5 data */
+ INELColor = kBlue+1,
+ /** Color used for CMS data */
+ INELGt0Color = kGreen+1,
+ /** Color used for ALICE data */
+ NSDColor = kMagenta+1
+ };
+ enum {
+ /** Style offset for mirror data */
+ MirrorOff = 4
+ };
+
+ //____________________________________________________________________
+ /**
+ * Set graph attributes based on trigger type and experiment.
+ *
+ * @param g Graph
+ * @param exp Experiment
+ * @param mirror True if mirrored data
+ * @param name Name of graph
+ * @param title Title of graph
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static void SetGraphAttributes(TGraph* g,
+ Int_t /*trig*/,
+ Int_t exp,
+ bool mirror,
+ const Char_t* name,
+ const Char_t* title)
+ {
+ Int_t color = 0;
+ switch (exp) {
+ case UA5: color = UA5Color; break;
+ case CMS: color = CMSColor; break;
+ case ALICE: color = ALICEColor; break;
+ case WIP: color = WIPColor; break;
+ case PYTHIA: color = PYTHIAColor; break;
+ }
+ Int_t style = 0;
+ switch (exp) {
+ case UA5: style = UA5Style; break;
+ case CMS: style = CMSStyle; break;
+ case ALICE: style = ALICEStyle; break;
+ case WIP: style = WIPStyle; break;
+ case PYTHIA: style = PYTHIAStyle; break;
+ }
+ Float_t size = g->GetMarkerSize();
+ switch (style) {
+ case 21:
+ case 25: size *= 0.8; break;
+ case 27: size *= 1.4; break;
+ case 33: size *= 1.4; break;
+ }
- if (mirror) style += MirrorOff;
+ if (mirror) style += MirrorOff;
- g->SetName(name);
- g->SetTitle(title);
- g->SetMarkerStyle(style);
- g->SetMarkerSize(size);
- g->SetMarkerColor(color);
- g->SetLineColor(color);
- g->SetFillColor(0);
- g->SetFillStyle(0);
- g->GetHistogram()->SetStats(kFALSE);
- g->GetHistogram()->SetXTitle("#eta");
- g->GetHistogram()->SetYTitle("#frac{1}{N} #frac{dN_{ch}}{#eta}");
-}
+ g->SetName(name);
+ g->SetTitle(title);
+ g->SetMarkerStyle(style);
+ g->SetMarkerSize(size);
+ g->SetMarkerColor(color);
+ g->SetLineColor(color);
+ g->SetFillColor(0);
+ g->SetFillStyle(0);
+ g->GetHistogram()->SetStats(kFALSE);
+ g->GetHistogram()->SetXTitle("#eta");
+ g->GetHistogram()->SetYTitle("#frac{1}{N} #frac{dN_{ch}}{#eta}");
+ }
-//____________________________________________________________________
-TGraphAsymmErrors*
-Pythia900INEL()
-{
+ //____________________________________________________________________
+ /**
+ * Get PYTHIA 900GeV INEL data
+ *
+ *
+ * @return Data graph
+ */
+ static TGraphAsymmErrors* Pythia900INEL()
+ {
- TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
- SetGraphAttributes(gre, INEL, PYTHIA, false, "pythia900Inel",
- "Pythia INEL");
- gre->SetPoint(0,-3.95,1.78199);
- gre->SetPointError(0, 0, 0, 0.0145305, 0.0145305);
- gre->SetPoint(1,-3.85,1.85486);
- gre->SetPointError(1,0,0,0.0148246,0.0148246);
- gre->SetPoint(2,-3.75,1.93886);
- gre->SetPointError(2,0,0,0.0151566,0.0151566);
- gre->SetPoint(3,-3.65,1.96055);
- gre->SetPointError(3,0,0,0.0152411,0.0152411);
- gre->SetPoint(4,-3.55,1.98756);
- gre->SetPointError(4,0,0,0.0153458,0.0153458);
- gre->SetPoint(5,-3.45,2.02844);
- gre->SetPointError(5,0,0,0.0155028,0.0155028);
- gre->SetPoint(6,-3.35,2.09585);
- gre->SetPointError(6,0,0,0.0157583,0.0157583);
- gre->SetPoint(7,-3.25,2.13732);
- gre->SetPointError(7,0,0,0.0159134,0.0159134);
- gre->SetPoint(8,-3.15,2.1686);
- gre->SetPointError(8,0,0,0.0160295,0.0160295);
- gre->SetPoint(9,-3.05,2.25296);
- gre->SetPointError(9,0,0,0.0163383,0.0163383);
- gre->SetPoint(10,-2.95,2.29265);
- gre->SetPointError(10,0,0,0.0164815,0.0164815);
- gre->SetPoint(11,-2.85,2.34799);
- gre->SetPointError(11,0,0,0.0166792,0.0166792);
- gre->SetPoint(12,-2.75,2.35652);
- gre->SetPointError(12,0,0,0.0167095,0.0167095);
- gre->SetPoint(13,-2.65,2.40545);
- gre->SetPointError(13,0,0,0.0168821,0.0168821);
- gre->SetPoint(14,-2.55,2.43934);
- gre->SetPointError(14,0,0,0.0170006,0.0170006);
- gre->SetPoint(15,-2.45,2.45735);
- gre->SetPointError(15,0,0,0.0170633,0.0170633);
- gre->SetPoint(16,-2.35,2.48945);
- gre->SetPointError(16,0,0,0.0171744,0.0171744);
- gre->SetPoint(17,-2.25,2.51635);
- gre->SetPointError(17,0,0,0.0172669,0.0172669);
- gre->SetPoint(18,-2.15,2.55047);
- gre->SetPointError(18,0,0,0.0173836,0.0173836);
- gre->SetPoint(19,-2.05,2.58021);
- gre->SetPointError(19,0,0,0.0174846,0.0174846);
- gre->SetPoint(20,-1.95,2.58732);
- gre->SetPointError(20,0,0,0.0175087,0.0175087);
- gre->SetPoint(21,-1.85,2.60095);
- gre->SetPointError(21,0,0,0.0175547,0.0175547);
- gre->SetPoint(22,-1.75,2.59941);
- gre->SetPointError(22,0,0,0.0175495,0.0175495);
- gre->SetPoint(23,-1.65,2.63021);
- gre->SetPointError(23,0,0,0.0176532,0.0176532);
- gre->SetPoint(24,-1.55,2.61043);
- gre->SetPointError(24,0,0,0.0175867,0.0175867);
- gre->SetPoint(25,-1.45,2.61363);
- gre->SetPointError(25,0,0,0.0175975,0.0175975);
- gre->SetPoint(26,-1.35,2.60829);
- gre->SetPointError(26,0,0,0.0175795,0.0175795);
- gre->SetPoint(27,-1.25,2.61434);
- gre->SetPointError(27,0,0,0.0175999,0.0175999);
- gre->SetPoint(28,-1.15,2.61327);
- gre->SetPointError(28,0,0,0.0175963,0.0175963);
- gre->SetPoint(29,-1.05,2.57145);
- gre->SetPointError(29,0,0,0.0174549,0.0174549);
- gre->SetPoint(30,-0.95,2.55723);
- gre->SetPointError(30,0,0,0.0174066,0.0174066);
- gre->SetPoint(31,-0.85,2.57879);
- gre->SetPointError(31,0,0,0.0174798,0.0174798);
- gre->SetPoint(32,-0.75,2.516);
- gre->SetPointError(32,0,0,0.0172657,0.0172657);
- gre->SetPoint(33,-0.65,2.53709);
- gre->SetPointError(33,0,0,0.0173379,0.0173379);
- gre->SetPoint(34,-0.55,2.51197);
- gre->SetPointError(34,0,0,0.0172519,0.0172519);
- gre->SetPoint(35,-0.45,2.44052);
- gre->SetPointError(35,0,0,0.0170047,0.0170047);
- gre->SetPoint(36,-0.35,2.44882);
- gre->SetPointError(36,0,0,0.0170336,0.0170336);
- gre->SetPoint(37,-0.25,2.45308);
- gre->SetPointError(37,0,0,0.0170484,0.0170484);
- gre->SetPoint(38,-0.15,2.4622);
- gre->SetPointError(38,0,0,0.0170801,0.0170801);
- gre->SetPoint(39,-0.05,2.45735);
- gre->SetPointError(39,0,0,0.0170633,0.0170633);
- gre->SetPoint(40,0.05,2.49254);
- gre->SetPointError(40,0,0,0.017185,0.017185);
- gre->SetPoint(41,0.15,2.49479);
- gre->SetPointError(41,0,0,0.0171928,0.0171928);
- gre->SetPoint(42,0.25,2.49289);
- gre->SetPointError(42,0,0,0.0171862,0.0171862);
- gre->SetPoint(43,0.35,2.4628);
- gre->SetPointError(43,0,0,0.0170822,0.0170822);
- gre->SetPoint(44,0.45,2.51422);
- gre->SetPointError(44,0,0,0.0172596,0.0172596);
- gre->SetPoint(45,0.55,2.51268);
- gre->SetPointError(45,0,0,0.0172543,0.0172543);
- gre->SetPoint(46,0.65,2.51066);
- gre->SetPointError(46,0,0,0.0172474,0.0172474);
- gre->SetPoint(47,0.75,2.53661);
- gre->SetPointError(47,0,0,0.0173363,0.0173363);
- gre->SetPoint(48,0.85,2.54479);
- gre->SetPointError(48,0,0,0.0173642,0.0173642);
- gre->SetPoint(49,0.95,2.55391);
- gre->SetPointError(49,0,0,0.0173953,0.0173953);
- gre->SetPoint(50,1.05,2.5872);
- gre->SetPointError(50,0,0,0.0175083,0.0175083);
- gre->SetPoint(51,1.15,2.60344);
- gre->SetPointError(51,0,0,0.0175631,0.0175631);
- gre->SetPoint(52,1.25,2.60616);
- gre->SetPointError(52,0,0,0.0175723,0.0175723);
- gre->SetPoint(53,1.35,2.62156);
- gre->SetPointError(53,0,0,0.0176242,0.0176242);
- gre->SetPoint(54,1.45,2.61173);
- gre->SetPointError(54,0,0,0.0175911,0.0175911);
- gre->SetPoint(55,1.55,2.60415);
- gre->SetPointError(55,0,0,0.0175655,0.0175655);
- gre->SetPoint(56,1.65,2.60723);
- gre->SetPointError(56,0,0,0.0175759,0.0175759);
- gre->SetPoint(57,1.75,2.60427);
- gre->SetPointError(57,0,0,0.0175659,0.0175659);
- gre->SetPoint(58,1.85,2.56765);
- gre->SetPointError(58,0,0,0.017442,0.017442);
- gre->SetPoint(59,1.95,2.58602);
- gre->SetPointError(59,0,0,0.0175043,0.0175043);
- gre->SetPoint(60,2.05,2.55936);
- gre->SetPointError(60,0,0,0.0174138,0.0174138);
- gre->SetPoint(61,2.15,2.54858);
- gre->SetPointError(61,0,0,0.0173771,0.0173771);
- gre->SetPoint(62,2.25,2.5205);
- gre->SetPointError(62,0,0,0.0172811,0.0172811);
- gre->SetPoint(63,2.35,2.49491);
- gre->SetPointError(63,0,0,0.0171932,0.0171932);
- gre->SetPoint(64,2.45,2.42773);
- gre->SetPointError(64,0,0,0.0169601,0.0169601);
- gre->SetPoint(65,2.55,2.42879);
- gre->SetPointError(65,0,0,0.0169638,0.0169638);
- gre->SetPoint(66,2.65,2.39372);
- gre->SetPointError(66,0,0,0.0168409,0.0168409);
- gre->SetPoint(67,2.75,2.38412);
- gre->SetPointError(67,0,0,0.0168071,0.0168071);
- gre->SetPoint(68,2.85,2.31896);
- gre->SetPointError(68,0,0,0.0165758,0.0165758);
- gre->SetPoint(69,2.95,2.26209);
- gre->SetPointError(69,0,0,0.0163713,0.0163713);
- gre->SetPoint(70,3.05,2.24313);
- gre->SetPointError(70,0,0,0.0163026,0.0163026);
- gre->SetPoint(71,3.15,2.20403);
- gre->SetPointError(71,0,0,0.0161599,0.0161599);
- gre->SetPoint(72,3.25,2.12855);
- gre->SetPointError(72,0,0,0.0158808,0.0158808);
- gre->SetPoint(73,3.35,2.13104);
- gre->SetPointError(73,0,0,0.01589,0.01589);
- gre->SetPoint(74,3.45,2.06339);
- gre->SetPointError(74,0,0,0.0156358,0.0156358);
- gre->SetPoint(75,3.55,1.9846);
- gre->SetPointError(75,0,0,0.0153343,0.0153343);
- gre->SetPoint(76,3.65,1.95391);
- gre->SetPointError(76,0,0,0.0152153,0.0152153);
- gre->SetPoint(77,3.75,1.87998);
- gre->SetPointError(77,0,0,0.0149247,0.0149247);
- gre->SetPoint(78,3.85,1.86256);
- gre->SetPointError(78,0,0,0.0148554,0.0148554);
- gre->SetPoint(79,3.95,1.77239);
- gre->SetPointError(79,0,0,0.0144913,0.0144913);
- gre->SetPoint(80,4.05,1.72855);
- gre->SetPointError(80,0,0,0.014311,0.014311);
- gre->SetPoint(81,4.15,1.69479);
- gre->SetPointError(81,0,0,0.0141705,0.0141705);
- gre->SetPoint(82,4.25,1.64147);
- gre->SetPointError(82,0,0,0.0139459,0.0139459);
- gre->SetPoint(83,4.35,1.58116);
- gre->SetPointError(83,0,0,0.0136873,0.0136873);
- gre->SetPoint(84,4.45,1.55735);
- gre->SetPointError(84,0,0,0.0135838,0.0135838);
- gre->SetPoint(85,4.55,1.48815);
- gre->SetPointError(85,0,0,0.0132786,0.0132786);
- gre->SetPoint(86,4.65,1.40853);
- gre->SetPointError(86,0,0,0.0129185,0.0129185);
- gre->SetPoint(87,4.75,1.36979);
- gre->SetPointError(87,0,0,0.0127396,0.0127396);
- gre->SetPoint(88,4.85,1.32666);
- gre->SetPointError(88,0,0,0.0125374,0.0125374);
- gre->SetPoint(89,4.95,1.29763);
- gre->SetPointError(89,0,0,0.0123995,0.0123995);
- gre->SetPoint(90,5.05,1.25533);
- gre->SetPointError(90,0,0,0.0121957,0.0121957);
- gre->SetPoint(91,5.15,1.20912);
- gre->SetPointError(91,0,0,0.0119692,0.0119692);
- gre->SetPoint(92,5.25,1.18839);
- gre->SetPointError(92,0,0,0.0118661,0.0118661);
- gre->SetPoint(93,5.35,1.15948);
- gre->SetPointError(93,0,0,0.0117209,0.0117209);
- gre->SetPoint(94,5.45,1.1141);
- gre->SetPointError(94,0,0,0.0114892,0.0114892);
- gre->SetPoint(95,5.55,1.06315);
- gre->SetPointError(95,0,0,0.0112235,0.0112235);
- gre->SetPoint(96,5.65,1.05213);
- gre->SetPointError(96,0,0,0.0111651,0.0111651);
- gre->SetPoint(97,5.75,1.02476);
- gre->SetPointError(97,0,0,0.011019,0.011019);
- gre->SetPoint(98,5.85,0.984834);
- gre->SetPointError(98,0,0,0.0108022,0.0108022);
- gre->SetPoint(99,5.95,0.952844);
- gre->SetPointError(99,0,0,0.0106253,0.0106253);
+ TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
+ SetGraphAttributes(gre, INEL, PYTHIA, false, "pythia900Inel",
+ "Pythia");
+ gre->SetPoint(0,-3.95,1.78199);
+ gre->SetPointError(0, 0, 0, 0.0145305, 0.0145305);
+ gre->SetPoint(1,-3.85,1.85486);
+ gre->SetPointError(1,0,0,0.0148246,0.0148246);
+ gre->SetPoint(2,-3.75,1.93886);
+ gre->SetPointError(2,0,0,0.0151566,0.0151566);
+ gre->SetPoint(3,-3.65,1.96055);
+ gre->SetPointError(3,0,0,0.0152411,0.0152411);
+ gre->SetPoint(4,-3.55,1.98756);
+ gre->SetPointError(4,0,0,0.0153458,0.0153458);
+ gre->SetPoint(5,-3.45,2.02844);
+ gre->SetPointError(5,0,0,0.0155028,0.0155028);
+ gre->SetPoint(6,-3.35,2.09585);
+ gre->SetPointError(6,0,0,0.0157583,0.0157583);
+ gre->SetPoint(7,-3.25,2.13732);
+ gre->SetPointError(7,0,0,0.0159134,0.0159134);
+ gre->SetPoint(8,-3.15,2.1686);
+ gre->SetPointError(8,0,0,0.0160295,0.0160295);
+ gre->SetPoint(9,-3.05,2.25296);
+ gre->SetPointError(9,0,0,0.0163383,0.0163383);
+ gre->SetPoint(10,-2.95,2.29265);
+ gre->SetPointError(10,0,0,0.0164815,0.0164815);
+ gre->SetPoint(11,-2.85,2.34799);
+ gre->SetPointError(11,0,0,0.0166792,0.0166792);
+ gre->SetPoint(12,-2.75,2.35652);
+ gre->SetPointError(12,0,0,0.0167095,0.0167095);
+ gre->SetPoint(13,-2.65,2.40545);
+ gre->SetPointError(13,0,0,0.0168821,0.0168821);
+ gre->SetPoint(14,-2.55,2.43934);
+ gre->SetPointError(14,0,0,0.0170006,0.0170006);
+ gre->SetPoint(15,-2.45,2.45735);
+ gre->SetPointError(15,0,0,0.0170633,0.0170633);
+ gre->SetPoint(16,-2.35,2.48945);
+ gre->SetPointError(16,0,0,0.0171744,0.0171744);
+ gre->SetPoint(17,-2.25,2.51635);
+ gre->SetPointError(17,0,0,0.0172669,0.0172669);
+ gre->SetPoint(18,-2.15,2.55047);
+ gre->SetPointError(18,0,0,0.0173836,0.0173836);
+ gre->SetPoint(19,-2.05,2.58021);
+ gre->SetPointError(19,0,0,0.0174846,0.0174846);
+ gre->SetPoint(20,-1.95,2.58732);
+ gre->SetPointError(20,0,0,0.0175087,0.0175087);
+ gre->SetPoint(21,-1.85,2.60095);
+ gre->SetPointError(21,0,0,0.0175547,0.0175547);
+ gre->SetPoint(22,-1.75,2.59941);
+ gre->SetPointError(22,0,0,0.0175495,0.0175495);
+ gre->SetPoint(23,-1.65,2.63021);
+ gre->SetPointError(23,0,0,0.0176532,0.0176532);
+ gre->SetPoint(24,-1.55,2.61043);
+ gre->SetPointError(24,0,0,0.0175867,0.0175867);
+ gre->SetPoint(25,-1.45,2.61363);
+ gre->SetPointError(25,0,0,0.0175975,0.0175975);
+ gre->SetPoint(26,-1.35,2.60829);
+ gre->SetPointError(26,0,0,0.0175795,0.0175795);
+ gre->SetPoint(27,-1.25,2.61434);
+ gre->SetPointError(27,0,0,0.0175999,0.0175999);
+ gre->SetPoint(28,-1.15,2.61327);
+ gre->SetPointError(28,0,0,0.0175963,0.0175963);
+ gre->SetPoint(29,-1.05,2.57145);
+ gre->SetPointError(29,0,0,0.0174549,0.0174549);
+ gre->SetPoint(30,-0.95,2.55723);
+ gre->SetPointError(30,0,0,0.0174066,0.0174066);
+ gre->SetPoint(31,-0.85,2.57879);
+ gre->SetPointError(31,0,0,0.0174798,0.0174798);
+ gre->SetPoint(32,-0.75,2.516);
+ gre->SetPointError(32,0,0,0.0172657,0.0172657);
+ gre->SetPoint(33,-0.65,2.53709);
+ gre->SetPointError(33,0,0,0.0173379,0.0173379);
+ gre->SetPoint(34,-0.55,2.51197);
+ gre->SetPointError(34,0,0,0.0172519,0.0172519);
+ gre->SetPoint(35,-0.45,2.44052);
+ gre->SetPointError(35,0,0,0.0170047,0.0170047);
+ gre->SetPoint(36,-0.35,2.44882);
+ gre->SetPointError(36,0,0,0.0170336,0.0170336);
+ gre->SetPoint(37,-0.25,2.45308);
+ gre->SetPointError(37,0,0,0.0170484,0.0170484);
+ gre->SetPoint(38,-0.15,2.4622);
+ gre->SetPointError(38,0,0,0.0170801,0.0170801);
+ gre->SetPoint(39,-0.05,2.45735);
+ gre->SetPointError(39,0,0,0.0170633,0.0170633);
+ gre->SetPoint(40,0.05,2.49254);
+ gre->SetPointError(40,0,0,0.017185,0.017185);
+ gre->SetPoint(41,0.15,2.49479);
+ gre->SetPointError(41,0,0,0.0171928,0.0171928);
+ gre->SetPoint(42,0.25,2.49289);
+ gre->SetPointError(42,0,0,0.0171862,0.0171862);
+ gre->SetPoint(43,0.35,2.4628);
+ gre->SetPointError(43,0,0,0.0170822,0.0170822);
+ gre->SetPoint(44,0.45,2.51422);
+ gre->SetPointError(44,0,0,0.0172596,0.0172596);
+ gre->SetPoint(45,0.55,2.51268);
+ gre->SetPointError(45,0,0,0.0172543,0.0172543);
+ gre->SetPoint(46,0.65,2.51066);
+ gre->SetPointError(46,0,0,0.0172474,0.0172474);
+ gre->SetPoint(47,0.75,2.53661);
+ gre->SetPointError(47,0,0,0.0173363,0.0173363);
+ gre->SetPoint(48,0.85,2.54479);
+ gre->SetPointError(48,0,0,0.0173642,0.0173642);
+ gre->SetPoint(49,0.95,2.55391);
+ gre->SetPointError(49,0,0,0.0173953,0.0173953);
+ gre->SetPoint(50,1.05,2.5872);
+ gre->SetPointError(50,0,0,0.0175083,0.0175083);
+ gre->SetPoint(51,1.15,2.60344);
+ gre->SetPointError(51,0,0,0.0175631,0.0175631);
+ gre->SetPoint(52,1.25,2.60616);
+ gre->SetPointError(52,0,0,0.0175723,0.0175723);
+ gre->SetPoint(53,1.35,2.62156);
+ gre->SetPointError(53,0,0,0.0176242,0.0176242);
+ gre->SetPoint(54,1.45,2.61173);
+ gre->SetPointError(54,0,0,0.0175911,0.0175911);
+ gre->SetPoint(55,1.55,2.60415);
+ gre->SetPointError(55,0,0,0.0175655,0.0175655);
+ gre->SetPoint(56,1.65,2.60723);
+ gre->SetPointError(56,0,0,0.0175759,0.0175759);
+ gre->SetPoint(57,1.75,2.60427);
+ gre->SetPointError(57,0,0,0.0175659,0.0175659);
+ gre->SetPoint(58,1.85,2.56765);
+ gre->SetPointError(58,0,0,0.017442,0.017442);
+ gre->SetPoint(59,1.95,2.58602);
+ gre->SetPointError(59,0,0,0.0175043,0.0175043);
+ gre->SetPoint(60,2.05,2.55936);
+ gre->SetPointError(60,0,0,0.0174138,0.0174138);
+ gre->SetPoint(61,2.15,2.54858);
+ gre->SetPointError(61,0,0,0.0173771,0.0173771);
+ gre->SetPoint(62,2.25,2.5205);
+ gre->SetPointError(62,0,0,0.0172811,0.0172811);
+ gre->SetPoint(63,2.35,2.49491);
+ gre->SetPointError(63,0,0,0.0171932,0.0171932);
+ gre->SetPoint(64,2.45,2.42773);
+ gre->SetPointError(64,0,0,0.0169601,0.0169601);
+ gre->SetPoint(65,2.55,2.42879);
+ gre->SetPointError(65,0,0,0.0169638,0.0169638);
+ gre->SetPoint(66,2.65,2.39372);
+ gre->SetPointError(66,0,0,0.0168409,0.0168409);
+ gre->SetPoint(67,2.75,2.38412);
+ gre->SetPointError(67,0,0,0.0168071,0.0168071);
+ gre->SetPoint(68,2.85,2.31896);
+ gre->SetPointError(68,0,0,0.0165758,0.0165758);
+ gre->SetPoint(69,2.95,2.26209);
+ gre->SetPointError(69,0,0,0.0163713,0.0163713);
+ gre->SetPoint(70,3.05,2.24313);
+ gre->SetPointError(70,0,0,0.0163026,0.0163026);
+ gre->SetPoint(71,3.15,2.20403);
+ gre->SetPointError(71,0,0,0.0161599,0.0161599);
+ gre->SetPoint(72,3.25,2.12855);
+ gre->SetPointError(72,0,0,0.0158808,0.0158808);
+ gre->SetPoint(73,3.35,2.13104);
+ gre->SetPointError(73,0,0,0.01589,0.01589);
+ gre->SetPoint(74,3.45,2.06339);
+ gre->SetPointError(74,0,0,0.0156358,0.0156358);
+ gre->SetPoint(75,3.55,1.9846);
+ gre->SetPointError(75,0,0,0.0153343,0.0153343);
+ gre->SetPoint(76,3.65,1.95391);
+ gre->SetPointError(76,0,0,0.0152153,0.0152153);
+ gre->SetPoint(77,3.75,1.87998);
+ gre->SetPointError(77,0,0,0.0149247,0.0149247);
+ gre->SetPoint(78,3.85,1.86256);
+ gre->SetPointError(78,0,0,0.0148554,0.0148554);
+ gre->SetPoint(79,3.95,1.77239);
+ gre->SetPointError(79,0,0,0.0144913,0.0144913);
+ gre->SetPoint(80,4.05,1.72855);
+ gre->SetPointError(80,0,0,0.014311,0.014311);
+ gre->SetPoint(81,4.15,1.69479);
+ gre->SetPointError(81,0,0,0.0141705,0.0141705);
+ gre->SetPoint(82,4.25,1.64147);
+ gre->SetPointError(82,0,0,0.0139459,0.0139459);
+ gre->SetPoint(83,4.35,1.58116);
+ gre->SetPointError(83,0,0,0.0136873,0.0136873);
+ gre->SetPoint(84,4.45,1.55735);
+ gre->SetPointError(84,0,0,0.0135838,0.0135838);
+ gre->SetPoint(85,4.55,1.48815);
+ gre->SetPointError(85,0,0,0.0132786,0.0132786);
+ gre->SetPoint(86,4.65,1.40853);
+ gre->SetPointError(86,0,0,0.0129185,0.0129185);
+ gre->SetPoint(87,4.75,1.36979);
+ gre->SetPointError(87,0,0,0.0127396,0.0127396);
+ gre->SetPoint(88,4.85,1.32666);
+ gre->SetPointError(88,0,0,0.0125374,0.0125374);
+ gre->SetPoint(89,4.95,1.29763);
+ gre->SetPointError(89,0,0,0.0123995,0.0123995);
+ gre->SetPoint(90,5.05,1.25533);
+ gre->SetPointError(90,0,0,0.0121957,0.0121957);
+ gre->SetPoint(91,5.15,1.20912);
+ gre->SetPointError(91,0,0,0.0119692,0.0119692);
+ gre->SetPoint(92,5.25,1.18839);
+ gre->SetPointError(92,0,0,0.0118661,0.0118661);
+ gre->SetPoint(93,5.35,1.15948);
+ gre->SetPointError(93,0,0,0.0117209,0.0117209);
+ gre->SetPoint(94,5.45,1.1141);
+ gre->SetPointError(94,0,0,0.0114892,0.0114892);
+ gre->SetPoint(95,5.55,1.06315);
+ gre->SetPointError(95,0,0,0.0112235,0.0112235);
+ gre->SetPoint(96,5.65,1.05213);
+ gre->SetPointError(96,0,0,0.0111651,0.0111651);
+ gre->SetPoint(97,5.75,1.02476);
+ gre->SetPointError(97,0,0,0.011019,0.011019);
+ gre->SetPoint(98,5.85,0.984834);
+ gre->SetPointError(98,0,0,0.0108022,0.0108022);
+ gre->SetPoint(99,5.95,0.952844);
+ gre->SetPointError(99,0,0,0.0106253,0.0106253);
- return gre;
-}
+ return gre;
+ }
-//____________________________________________________________________
-TGraphAsymmErrors*
-Pythia900NSD()
-{
+ //____________________________________________________________________
+ /**
+ * Get PYTHIA 900GeV NSD data
+ *
+ *
+ * @return Data graph
+ */
+ static TGraphAsymmErrors* Pythia900NSD()
+ {
- TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
- SetGraphAttributes(gre, NSD, PYTHIA, false, "pythia900NSD",
- "Pythia NSD");
+ TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
+ SetGraphAttributes(gre, NSD, PYTHIA, false, "pythia900NSD",
+ "Pythia");
- gre->SetPoint(0,-3.95,2.11766);
- gre->SetPointError(0,0,0,0.0179417,0.0179417);
- gre->SetPoint(1,-3.85,2.20415);
- gre->SetPointError(1,0,0,0.0183045,0.0183045);
- gre->SetPoint(2,-3.75,2.30949);
- gre->SetPointError(2,0,0,0.0187368,0.0187368);
- gre->SetPoint(3,-3.65,2.34582);
- gre->SetPointError(3,0,0,0.0188836,0.0188836);
- gre->SetPoint(4,-3.55,2.38322);
- gre->SetPointError(4,0,0,0.0190335,0.0190335);
- gre->SetPoint(5,-3.45,2.43353);
- gre->SetPointError(5,0,0,0.0192334,0.0192334);
- gre->SetPoint(6,-3.35,2.51106);
- gre->SetPointError(6,0,0,0.0195373,0.0195373);
- gre->SetPoint(7,-3.25,2.56578);
- gre->SetPointError(7,0,0,0.0197491,0.0197491);
- gre->SetPoint(8,-3.15,2.60515);
- gre->SetPointError(8,0,0,0.0199,0.0199);
- gre->SetPoint(9,-3.05,2.7105);
- gre->SetPointError(9,0,0,0.0202984,0.0202984);
- gre->SetPoint(10,-2.95,2.77008);
- gre->SetPointError(10,0,0,0.0205203,0.0205203);
- gre->SetPoint(11,-2.85,2.83332);
- gre->SetPointError(11,0,0,0.0207532,0.0207532);
- gre->SetPoint(12,-2.75,2.84715);
- gre->SetPointError(12,0,0,0.0208038,0.0208038);
- gre->SetPoint(13,-2.65,2.91693);
- gre->SetPointError(13,0,0,0.0210571,0.0210571);
- gre->SetPoint(14,-2.55,2.95797);
- gre->SetPointError(14,0,0,0.0212048,0.0212048);
- gre->SetPoint(15,-2.45,2.97499);
- gre->SetPointError(15,0,0,0.0212657,0.0212657);
- gre->SetPoint(16,-2.35,3.01345);
- gre->SetPointError(16,0,0,0.0214027,0.0214027);
- gre->SetPoint(17,-2.25,3.04659);
- gre->SetPointError(17,0,0,0.0215201,0.0215201);
- gre->SetPoint(18,-2.15,3.09341);
- gre->SetPointError(18,0,0,0.0216848,0.0216848);
- gre->SetPoint(19,-2.05,3.13187);
- gre->SetPointError(19,0,0,0.0218192,0.0218192);
- gre->SetPoint(20,-1.95,3.13917);
- gre->SetPointError(20,0,0,0.0218446,0.0218446);
- gre->SetPoint(21,-1.85,3.16911);
- gre->SetPointError(21,0,0,0.0219485,0.0219485);
- gre->SetPoint(22,-1.75,3.15665);
- gre->SetPointError(22,0,0,0.0219053,0.0219053);
- gre->SetPoint(23,-1.65,3.19693);
- gre->SetPointError(23,0,0,0.0220446,0.0220446);
- gre->SetPoint(24,-1.55,3.17002);
- gre->SetPointError(24,0,0,0.0219517,0.0219517);
- gre->SetPoint(25,-1.45,3.18538);
- gre->SetPointError(25,0,0,0.0220048,0.0220048);
- gre->SetPoint(26,-1.35,3.18066);
- gre->SetPointError(26,0,0,0.0219885,0.0219885);
- gre->SetPoint(27,-1.25,3.19754);
- gre->SetPointError(27,0,0,0.0220467,0.0220467);
- gre->SetPoint(28,-1.15,3.18021);
- gre->SetPointError(28,0,0,0.0219869,0.0219869);
- gre->SetPoint(29,-1.05,3.13111);
- gre->SetPointError(29,0,0,0.0218165,0.0218165);
- gre->SetPoint(30,-0.95,3.12153);
- gre->SetPointError(30,0,0,0.0217831,0.0217831);
- gre->SetPoint(31,-0.85,3.14798);
- gre->SetPointError(31,0,0,0.0218752,0.0218752);
- gre->SetPoint(32,-0.75,3.07912);
- gre->SetPointError(32,0,0,0.0216347,0.0216347);
- gre->SetPoint(33,-0.65,3.10207);
- gre->SetPointError(33,0,0,0.0217151,0.0217151);
- gre->SetPoint(34,-0.55,3.06346);
- gre->SetPointError(34,0,0,0.0215796,0.0215796);
- gre->SetPoint(35,-0.45,2.97651);
- gre->SetPointError(35,0,0,0.0212711,0.0212711);
- gre->SetPoint(36,-0.35,2.98715);
- gre->SetPointError(36,0,0,0.0213091,0.0213091);
- gre->SetPoint(37,-0.25,2.98548);
- gre->SetPointError(37,0,0,0.0213032,0.0213032);
- gre->SetPoint(38,-0.15,3.00555);
- gre->SetPointError(38,0,0,0.0213746,0.0213746);
- gre->SetPoint(39,-0.05,3.01193);
- gre->SetPointError(39,0,0,0.0213973,0.0213973);
- gre->SetPoint(40,0.05,3.04385);
- gre->SetPointError(40,0,0,0.0215104,0.0215104);
- gre->SetPoint(41,0.15,3.04933);
- gre->SetPointError(41,0,0,0.0215297,0.0215297);
- gre->SetPoint(42,0.25,3.04659);
- gre->SetPointError(42,0,0,0.0215201,0.0215201);
- gre->SetPoint(43,0.35,3.00813);
- gre->SetPointError(43,0,0,0.0213838,0.0213838);
- gre->SetPoint(44,0.45,3.06666);
- gre->SetPointError(44,0,0,0.0215908,0.0215908);
- gre->SetPoint(45,0.55,3.07167);
- gre->SetPointError(45,0,0,0.0216085,0.0216085);
- gre->SetPoint(46,0.65,3.0659);
- gre->SetPointError(46,0,0,0.0215881,0.0215881);
- gre->SetPoint(47,0.75,3.09159);
- gre->SetPointError(47,0,0,0.0216784,0.0216784);
- gre->SetPoint(48,0.85,3.10846);
- gre->SetPointError(48,0,0,0.0217375,0.0217375);
- gre->SetPoint(49,0.95,3.11925);
- gre->SetPointError(49,0,0,0.0217752,0.0217752);
- gre->SetPoint(50,1.05,3.15558);
- gre->SetPointError(50,0,0,0.0219016,0.0219016);
- gre->SetPoint(51,1.15,3.16911);
- gre->SetPointError(51,0,0,0.0219485,0.0219485);
- gre->SetPoint(52,1.25,3.17246);
- gre->SetPointError(52,0,0,0.0219601,0.0219601);
- gre->SetPoint(53,1.35,3.19146);
- gre->SetPointError(53,0,0,0.0220258,0.0220258);
- gre->SetPoint(54,1.45,3.17458);
- gre->SetPointError(54,0,0,0.0219675,0.0219675);
- gre->SetPoint(55,1.55,3.16866);
- gre->SetPointError(55,0,0,0.0219469,0.0219469);
- gre->SetPoint(56,1.65,3.16592);
- gre->SetPointError(56,0,0,0.0219375,0.0219375);
- gre->SetPoint(57,1.75,3.16394);
- gre->SetPointError(57,0,0,0.0219306,0.0219306);
- gre->SetPoint(58,1.85,3.11956);
- gre->SetPointError(58,0,0,0.0217762,0.0217762);
- gre->SetPoint(59,1.95,3.14646);
- gre->SetPointError(59,0,0,0.02187,0.02187);
- gre->SetPoint(60,2.05,3.10147);
- gre->SetPointError(60,0,0,0.021713,0.021713);
- gre->SetPoint(61,2.15,3.09356);
- gre->SetPointError(61,0,0,0.0216853,0.0216853);
- gre->SetPoint(62,2.25,3.05328);
- gre->SetPointError(62,0,0,0.0215437,0.0215437);
- gre->SetPoint(63,2.35,3.01953);
- gre->SetPointError(63,0,0,0.0214243,0.0214243);
- gre->SetPoint(64,2.45,2.9373);
- gre->SetPointError(64,0,0,0.0211305,0.0211305);
- gre->SetPoint(65,2.55,2.92772);
- gre->SetPointError(65,0,0,0.0210961,0.0210961);
- gre->SetPoint(66,2.65,2.89154);
- gre->SetPointError(66,0,0,0.0209653,0.0209653);
- gre->SetPoint(67,2.75,2.87619);
- gre->SetPointError(67,0,0,0.0209096,0.0209096);
- gre->SetPoint(68,2.85,2.78924);
- gre->SetPointError(68,0,0,0.0205911,0.0205911);
- gre->SetPoint(69,2.95,2.72159);
- gre->SetPointError(69,0,0,0.0203399,0.0203399);
- gre->SetPoint(70,3.05,2.69089);
- gre->SetPointError(70,0,0,0.0202248,0.0202248);
- gre->SetPoint(71,3.15,2.64939);
- gre->SetPointError(71,0,0,0.0200682,0.0200682);
- gre->SetPoint(72,3.25,2.55545);
- gre->SetPointError(72,0,0,0.0197092,0.0197092);
- gre->SetPoint(73,3.35,2.56745);
- gre->SetPointError(73,0,0,0.0197555,0.0197555);
- gre->SetPoint(74,3.45,2.47503);
- gre->SetPointError(74,0,0,0.0193967,0.0193967);
- gre->SetPoint(75,3.55,2.36741);
- gre->SetPointError(75,0,0,0.0189703,0.0189703);
- gre->SetPoint(76,3.65,2.33412);
- gre->SetPointError(76,0,0,0.0188364,0.0188364);
- gre->SetPoint(77,3.75,2.2385);
- gre->SetPointError(77,0,0,0.0184466,0.0184466);
- gre->SetPoint(78,3.85,2.21768);
- gre->SetPointError(78,0,0,0.0183606,0.0183606);
- gre->SetPoint(79,3.95,2.1055);
- gre->SetPointError(79,0,0,0.0178901,0.0178901);
- gre->SetPoint(80,4.05,2.05047);
- gre->SetPointError(80,0,0,0.0176548,0.0176548);
- gre->SetPoint(81,4.15,2.00486);
- gre->SetPointError(81,0,0,0.0174574,0.0174574);
- gre->SetPoint(82,4.25,1.94573);
- gre->SetPointError(82,0,0,0.017198,0.017198);
- gre->SetPoint(83,4.35,1.87064);
- gre->SetPointError(83,0,0,0.0168629,0.0168629);
- gre->SetPoint(84,4.45,1.83735);
- gre->SetPointError(84,0,0,0.0167122,0.0167122);
- gre->SetPoint(85,4.55,1.75314);
- gre->SetPointError(85,0,0,0.0163247,0.0163247);
- gre->SetPoint(86,4.65,1.65828);
- gre->SetPointError(86,0,0,0.0158769,0.0158769);
- gre->SetPoint(87,4.75,1.60751);
- gre->SetPointError(87,0,0,0.015632,0.015632);
- gre->SetPoint(88,4.85,1.56312);
- gre->SetPointError(88,0,0,0.0154146,0.0154146);
- gre->SetPoint(89,4.95,1.52117);
- gre->SetPointError(89,0,0,0.0152064,0.0152064);
- gre->SetPoint(90,5.05,1.46553);
- gre->SetPointError(90,0,0,0.0149257,0.0149257);
- gre->SetPoint(91,5.15,1.42038);
- gre->SetPointError(91,0,0,0.014694,0.014694);
- gre->SetPoint(92,5.25,1.38816);
- gre->SetPointError(92,0,0,0.0145263,0.0145263);
- gre->SetPoint(93,5.35,1.35046);
- gre->SetPointError(93,0,0,0.0143277,0.0143277);
- gre->SetPoint(94,5.45,1.30075);
- gre->SetPointError(94,0,0,0.0140616,0.0140616);
- gre->SetPoint(95,5.55,1.24025);
- gre->SetPointError(95,0,0,0.0137307,0.0137307);
- gre->SetPoint(96,5.65,1.21806);
- gre->SetPointError(96,0,0,0.0136073,0.0136073);
- gre->SetPoint(97,5.75,1.19435);
- gre->SetPointError(97,0,0,0.0134742,0.0134742);
- gre->SetPoint(98,5.85,1.14175);
- gre->SetPointError(98,0,0,0.0131741,0.0131741);
- gre->SetPoint(99,5.95,1.09235);
- gre->SetPointError(99,0,0,0.012886,0.012886);
+ gre->SetPoint(0,-3.95,2.11766);
+ gre->SetPointError(0,0,0,0.0179417,0.0179417);
+ gre->SetPoint(1,-3.85,2.20415);
+ gre->SetPointError(1,0,0,0.0183045,0.0183045);
+ gre->SetPoint(2,-3.75,2.30949);
+ gre->SetPointError(2,0,0,0.0187368,0.0187368);
+ gre->SetPoint(3,-3.65,2.34582);
+ gre->SetPointError(3,0,0,0.0188836,0.0188836);
+ gre->SetPoint(4,-3.55,2.38322);
+ gre->SetPointError(4,0,0,0.0190335,0.0190335);
+ gre->SetPoint(5,-3.45,2.43353);
+ gre->SetPointError(5,0,0,0.0192334,0.0192334);
+ gre->SetPoint(6,-3.35,2.51106);
+ gre->SetPointError(6,0,0,0.0195373,0.0195373);
+ gre->SetPoint(7,-3.25,2.56578);
+ gre->SetPointError(7,0,0,0.0197491,0.0197491);
+ gre->SetPoint(8,-3.15,2.60515);
+ gre->SetPointError(8,0,0,0.0199,0.0199);
+ gre->SetPoint(9,-3.05,2.7105);
+ gre->SetPointError(9,0,0,0.0202984,0.0202984);
+ gre->SetPoint(10,-2.95,2.77008);
+ gre->SetPointError(10,0,0,0.0205203,0.0205203);
+ gre->SetPoint(11,-2.85,2.83332);
+ gre->SetPointError(11,0,0,0.0207532,0.0207532);
+ gre->SetPoint(12,-2.75,2.84715);
+ gre->SetPointError(12,0,0,0.0208038,0.0208038);
+ gre->SetPoint(13,-2.65,2.91693);
+ gre->SetPointError(13,0,0,0.0210571,0.0210571);
+ gre->SetPoint(14,-2.55,2.95797);
+ gre->SetPointError(14,0,0,0.0212048,0.0212048);
+ gre->SetPoint(15,-2.45,2.97499);
+ gre->SetPointError(15,0,0,0.0212657,0.0212657);
+ gre->SetPoint(16,-2.35,3.01345);
+ gre->SetPointError(16,0,0,0.0214027,0.0214027);
+ gre->SetPoint(17,-2.25,3.04659);
+ gre->SetPointError(17,0,0,0.0215201,0.0215201);
+ gre->SetPoint(18,-2.15,3.09341);
+ gre->SetPointError(18,0,0,0.0216848,0.0216848);
+ gre->SetPoint(19,-2.05,3.13187);
+ gre->SetPointError(19,0,0,0.0218192,0.0218192);
+ gre->SetPoint(20,-1.95,3.13917);
+ gre->SetPointError(20,0,0,0.0218446,0.0218446);
+ gre->SetPoint(21,-1.85,3.16911);
+ gre->SetPointError(21,0,0,0.0219485,0.0219485);
+ gre->SetPoint(22,-1.75,3.15665);
+ gre->SetPointError(22,0,0,0.0219053,0.0219053);
+ gre->SetPoint(23,-1.65,3.19693);
+ gre->SetPointError(23,0,0,0.0220446,0.0220446);
+ gre->SetPoint(24,-1.55,3.17002);
+ gre->SetPointError(24,0,0,0.0219517,0.0219517);
+ gre->SetPoint(25,-1.45,3.18538);
+ gre->SetPointError(25,0,0,0.0220048,0.0220048);
+ gre->SetPoint(26,-1.35,3.18066);
+ gre->SetPointError(26,0,0,0.0219885,0.0219885);
+ gre->SetPoint(27,-1.25,3.19754);
+ gre->SetPointError(27,0,0,0.0220467,0.0220467);
+ gre->SetPoint(28,-1.15,3.18021);
+ gre->SetPointError(28,0,0,0.0219869,0.0219869);
+ gre->SetPoint(29,-1.05,3.13111);
+ gre->SetPointError(29,0,0,0.0218165,0.0218165);
+ gre->SetPoint(30,-0.95,3.12153);
+ gre->SetPointError(30,0,0,0.0217831,0.0217831);
+ gre->SetPoint(31,-0.85,3.14798);
+ gre->SetPointError(31,0,0,0.0218752,0.0218752);
+ gre->SetPoint(32,-0.75,3.07912);
+ gre->SetPointError(32,0,0,0.0216347,0.0216347);
+ gre->SetPoint(33,-0.65,3.10207);
+ gre->SetPointError(33,0,0,0.0217151,0.0217151);
+ gre->SetPoint(34,-0.55,3.06346);
+ gre->SetPointError(34,0,0,0.0215796,0.0215796);
+ gre->SetPoint(35,-0.45,2.97651);
+ gre->SetPointError(35,0,0,0.0212711,0.0212711);
+ gre->SetPoint(36,-0.35,2.98715);
+ gre->SetPointError(36,0,0,0.0213091,0.0213091);
+ gre->SetPoint(37,-0.25,2.98548);
+ gre->SetPointError(37,0,0,0.0213032,0.0213032);
+ gre->SetPoint(38,-0.15,3.00555);
+ gre->SetPointError(38,0,0,0.0213746,0.0213746);
+ gre->SetPoint(39,-0.05,3.01193);
+ gre->SetPointError(39,0,0,0.0213973,0.0213973);
+ gre->SetPoint(40,0.05,3.04385);
+ gre->SetPointError(40,0,0,0.0215104,0.0215104);
+ gre->SetPoint(41,0.15,3.04933);
+ gre->SetPointError(41,0,0,0.0215297,0.0215297);
+ gre->SetPoint(42,0.25,3.04659);
+ gre->SetPointError(42,0,0,0.0215201,0.0215201);
+ gre->SetPoint(43,0.35,3.00813);
+ gre->SetPointError(43,0,0,0.0213838,0.0213838);
+ gre->SetPoint(44,0.45,3.06666);
+ gre->SetPointError(44,0,0,0.0215908,0.0215908);
+ gre->SetPoint(45,0.55,3.07167);
+ gre->SetPointError(45,0,0,0.0216085,0.0216085);
+ gre->SetPoint(46,0.65,3.0659);
+ gre->SetPointError(46,0,0,0.0215881,0.0215881);
+ gre->SetPoint(47,0.75,3.09159);
+ gre->SetPointError(47,0,0,0.0216784,0.0216784);
+ gre->SetPoint(48,0.85,3.10846);
+ gre->SetPointError(48,0,0,0.0217375,0.0217375);
+ gre->SetPoint(49,0.95,3.11925);
+ gre->SetPointError(49,0,0,0.0217752,0.0217752);
+ gre->SetPoint(50,1.05,3.15558);
+ gre->SetPointError(50,0,0,0.0219016,0.0219016);
+ gre->SetPoint(51,1.15,3.16911);
+ gre->SetPointError(51,0,0,0.0219485,0.0219485);
+ gre->SetPoint(52,1.25,3.17246);
+ gre->SetPointError(52,0,0,0.0219601,0.0219601);
+ gre->SetPoint(53,1.35,3.19146);
+ gre->SetPointError(53,0,0,0.0220258,0.0220258);
+ gre->SetPoint(54,1.45,3.17458);
+ gre->SetPointError(54,0,0,0.0219675,0.0219675);
+ gre->SetPoint(55,1.55,3.16866);
+ gre->SetPointError(55,0,0,0.0219469,0.0219469);
+ gre->SetPoint(56,1.65,3.16592);
+ gre->SetPointError(56,0,0,0.0219375,0.0219375);
+ gre->SetPoint(57,1.75,3.16394);
+ gre->SetPointError(57,0,0,0.0219306,0.0219306);
+ gre->SetPoint(58,1.85,3.11956);
+ gre->SetPointError(58,0,0,0.0217762,0.0217762);
+ gre->SetPoint(59,1.95,3.14646);
+ gre->SetPointError(59,0,0,0.02187,0.02187);
+ gre->SetPoint(60,2.05,3.10147);
+ gre->SetPointError(60,0,0,0.021713,0.021713);
+ gre->SetPoint(61,2.15,3.09356);
+ gre->SetPointError(61,0,0,0.0216853,0.0216853);
+ gre->SetPoint(62,2.25,3.05328);
+ gre->SetPointError(62,0,0,0.0215437,0.0215437);
+ gre->SetPoint(63,2.35,3.01953);
+ gre->SetPointError(63,0,0,0.0214243,0.0214243);
+ gre->SetPoint(64,2.45,2.9373);
+ gre->SetPointError(64,0,0,0.0211305,0.0211305);
+ gre->SetPoint(65,2.55,2.92772);
+ gre->SetPointError(65,0,0,0.0210961,0.0210961);
+ gre->SetPoint(66,2.65,2.89154);
+ gre->SetPointError(66,0,0,0.0209653,0.0209653);
+ gre->SetPoint(67,2.75,2.87619);
+ gre->SetPointError(67,0,0,0.0209096,0.0209096);
+ gre->SetPoint(68,2.85,2.78924);
+ gre->SetPointError(68,0,0,0.0205911,0.0205911);
+ gre->SetPoint(69,2.95,2.72159);
+ gre->SetPointError(69,0,0,0.0203399,0.0203399);
+ gre->SetPoint(70,3.05,2.69089);
+ gre->SetPointError(70,0,0,0.0202248,0.0202248);
+ gre->SetPoint(71,3.15,2.64939);
+ gre->SetPointError(71,0,0,0.0200682,0.0200682);
+ gre->SetPoint(72,3.25,2.55545);
+ gre->SetPointError(72,0,0,0.0197092,0.0197092);
+ gre->SetPoint(73,3.35,2.56745);
+ gre->SetPointError(73,0,0,0.0197555,0.0197555);
+ gre->SetPoint(74,3.45,2.47503);
+ gre->SetPointError(74,0,0,0.0193967,0.0193967);
+ gre->SetPoint(75,3.55,2.36741);
+ gre->SetPointError(75,0,0,0.0189703,0.0189703);
+ gre->SetPoint(76,3.65,2.33412);
+ gre->SetPointError(76,0,0,0.0188364,0.0188364);
+ gre->SetPoint(77,3.75,2.2385);
+ gre->SetPointError(77,0,0,0.0184466,0.0184466);
+ gre->SetPoint(78,3.85,2.21768);
+ gre->SetPointError(78,0,0,0.0183606,0.0183606);
+ gre->SetPoint(79,3.95,2.1055);
+ gre->SetPointError(79,0,0,0.0178901,0.0178901);
+ gre->SetPoint(80,4.05,2.05047);
+ gre->SetPointError(80,0,0,0.0176548,0.0176548);
+ gre->SetPoint(81,4.15,2.00486);
+ gre->SetPointError(81,0,0,0.0174574,0.0174574);
+ gre->SetPoint(82,4.25,1.94573);
+ gre->SetPointError(82,0,0,0.017198,0.017198);
+ gre->SetPoint(83,4.35,1.87064);
+ gre->SetPointError(83,0,0,0.0168629,0.0168629);
+ gre->SetPoint(84,4.45,1.83735);
+ gre->SetPointError(84,0,0,0.0167122,0.0167122);
+ gre->SetPoint(85,4.55,1.75314);
+ gre->SetPointError(85,0,0,0.0163247,0.0163247);
+ gre->SetPoint(86,4.65,1.65828);
+ gre->SetPointError(86,0,0,0.0158769,0.0158769);
+ gre->SetPoint(87,4.75,1.60751);
+ gre->SetPointError(87,0,0,0.015632,0.015632);
+ gre->SetPoint(88,4.85,1.56312);
+ gre->SetPointError(88,0,0,0.0154146,0.0154146);
+ gre->SetPoint(89,4.95,1.52117);
+ gre->SetPointError(89,0,0,0.0152064,0.0152064);
+ gre->SetPoint(90,5.05,1.46553);
+ gre->SetPointError(90,0,0,0.0149257,0.0149257);
+ gre->SetPoint(91,5.15,1.42038);
+ gre->SetPointError(91,0,0,0.014694,0.014694);
+ gre->SetPoint(92,5.25,1.38816);
+ gre->SetPointError(92,0,0,0.0145263,0.0145263);
+ gre->SetPoint(93,5.35,1.35046);
+ gre->SetPointError(93,0,0,0.0143277,0.0143277);
+ gre->SetPoint(94,5.45,1.30075);
+ gre->SetPointError(94,0,0,0.0140616,0.0140616);
+ gre->SetPoint(95,5.55,1.24025);
+ gre->SetPointError(95,0,0,0.0137307,0.0137307);
+ gre->SetPoint(96,5.65,1.21806);
+ gre->SetPointError(96,0,0,0.0136073,0.0136073);
+ gre->SetPoint(97,5.75,1.19435);
+ gre->SetPointError(97,0,0,0.0134742,0.0134742);
+ gre->SetPoint(98,5.85,1.14175);
+ gre->SetPointError(98,0,0,0.0131741,0.0131741);
+ gre->SetPoint(99,5.95,1.09235);
+ gre->SetPointError(99,0,0,0.012886,0.012886);
- return gre;
-}
+ return gre;
+ }
-//____________________________________________________________________
-/**
- * Get the UA5 NSD data for pp at @f$ \sqrt{s} = 900GeV@f$
- * p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
- *
- * @param mirrored Wether to produce the mirrored plot
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* UA5Nsd(Bool_t mirrored=false)
-{
- //UA5 data NSD - p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
- double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
- 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
- 4.625 };
- double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
- 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
- double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
- 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
- double y[] = { 3.48, 3.38, 3.52, 3.68, 3.71, 3.86, 3.76, 3.66, 3.72,
- 3.69, 3.56, 3.41, 3.15, 3.09, 2.74, 2.73, 2.32, 1.99, 1.69 };
- double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
- 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
- double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
- 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
- const int np = 19;
- double xm[np];
- double exmm[np];
- double expm[np];
- double ym[np];
- double eymm[np];
- double eypm[np];
- for (int i = 0; i < np; i++) {
- int j = np-1-i;
- xm[i] = -x[j];
- exmm[i] = exm[j];
- expm[i] = exp[j];
- ym[i] = y[j];
- eymm[i] = eym[j];
- eypm[i] = eyp[j];
- }
-
- TGraphAsymmErrors* g = new TGraphAsymmErrors(19,x, y, exm, exp, eym, eyp);
- TGraphAsymmErrors* gm = new TGraphAsymmErrors(19,xm,ym,exmm,expm,eymm,eypm);
- SetGraphAttributes(g, NSD, UA5, false,"ua5_nsd", "UA5 NSD");
- SetGraphAttributes(gm, NSD, UA5, true,"ua5_nsd_mirrored",
- "UA5 NSD (mirrored)");
+ //____________________________________________________________________
+ /**
+ * Get the UA5 NSD data for pp at @f$ \sqrt{s} = 900GeV@f$
+ * p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
+ *
+ * @param mirrored Wether to produce the mirrored plot
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* UA5Nsd(Bool_t mirrored=false)
+ {
+ //UA5 data NSD - p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
+ double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875,2.125,
+ 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125,4.375,
+ 4.625 };
+ double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
+ double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
+ double y[] = { 3.48, 3.38, 3.52, 3.68, 3.71, 3.86, 3.76, 3.66, 3.72,
+ 3.69, 3.56, 3.41, 3.15, 3.09, 2.74, 2.73, 2.32, 1.99, 1.69 };
+ double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
+ 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1,
+ 0.13 };
+ double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
+ 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1,
+ 0.13 };
+ const int np = 19;
+ double xm[np];
+ double exmm[np];
+ double expm[np];
+ double ym[np];
+ double eymm[np];
+ double eypm[np];
+ for (int i = 0; i < np; i++) {
+ int j = np-1-i;
+ xm[i] = -x[j];
+ exmm[i] = exm[j];
+ expm[i] = exp[j];
+ ym[i] = y[j];
+ eymm[i] = eym[j];
+ eypm[i] = eyp[j];
+ }
- return (mirrored ? gm : g);
-}
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(19,x, y, exm, exp, eym, eyp);
+ TGraphAsymmErrors* gm = new TGraphAsymmErrors(19,xm,ym,exmm,expm,eymm,eypm);
+ SetGraphAttributes(g, NSD, UA5, false,"ua5_nsd", "UA5");
+ SetGraphAttributes(gm, NSD, UA5, true,"ua5_nsd_mirrored",
+ "UA5 (mirrored)");
-//____________________________________________________________________
-/**
- * Get the UA5 INEL data for pp at @f$ \sqrt{s} = 900GeV@f$
- * p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
- *
- * @param mirrored Wether to produce the mirrored plot
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* UA5Inel(Bool_t mirrored=false)
-{
- //UA5 data INEL - p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
- double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
- 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
- 4.625 };
- double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
- 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
- double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
- 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
- double y[] = { 3.14, 3.04, 3.17, 3.33, 3.33, 3.53, 3.46, 3.41, 3.45,
- 3.39, 3.07, 3.07, 2.93, 2.93, 2.55, 2.48, 2.18, 1.91, 1.52 };
- double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
- 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
- double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
- 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
- const int np = 19;
- double xm[np];
- double exmm[np];
- double expm[np];
- double ym[np];
- double eymm[np];
- double eypm[np];
- for (int i = 0; i < np; i++) {
- int j = np-1-i;
- xm[i] = -x[j];
- exmm[i] = exm[j];
- expm[i] = exp[j];
- ym[i] = y[j];
- eymm[i] = eym[j];
- eypm[i] = eyp[j];
+ return (mirrored ? gm : g);
}
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np,x, y, exm, exp, eym, eyp);
- TGraphAsymmErrors* gm = new TGraphAsymmErrors(np,xm,ym,exmm,expm,eymm,eypm);
- SetGraphAttributes(g, INEL, UA5, false, "ua5_inel", "UA5 INEL");
- SetGraphAttributes(gm, INEL, UA5, true, "ua5_inel_mirrored",
- "UA5 INEL (mirrored)");
-
- return (mirrored ? gm : g);
-}
+ //____________________________________________________________________
+ /**
+ * Get the UA5 INEL data for pp at @f$ \sqrt{s} = 900GeV@f$
+ * p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
+ *
+ * @param mirrored Wether to produce the mirrored plot
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* UA5Inel(Bool_t mirrored=false)
+ {
+ //UA5 data INEL - p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
+ double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875,2.125,
+ 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125,4.375,
+ 4.625 };
+ double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
+ double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+ 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
+ double y[] = { 3.14, 3.04, 3.17, 3.33, 3.33, 3.53, 3.46, 3.41, 3.45,
+ 3.39, 3.07, 3.07, 2.93, 2.93, 2.55, 2.48, 2.18, 1.91, 1.52 };
+ double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
+ 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1,
+ 0.13 };
+ double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
+ 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1,
+ 0.13 };
+ const int np = 19;
+ double xm[np];
+ double exmm[np];
+ double expm[np];
+ double ym[np];
+ double eymm[np];
+ double eypm[np];
+ for (int i = 0; i < np; i++) {
+ int j = np-1-i;
+ xm[i] = -x[j];
+ exmm[i] = exm[j];
+ expm[i] = exp[j];
+ ym[i] = y[j];
+ eymm[i] = eym[j];
+ eypm[i] = eyp[j];
+ }
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np,x, y, exm, exp, eym, eyp);
+ TGraphAsymmErrors* gm = new TGraphAsymmErrors(np,xm,ym,exmm,expm,eymm,eypm);
-//____________________________________________________________________
-/**
- * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$ \sqrt{s}
- * = 900GeV@f$
- * p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralInel900()
-{
- // INEL - p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
- TGraphAsymmErrors* g = 0;
- double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3,
- -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 3.28, 3.28, 3.22, 3.12, 3.06, 3.02, 2.98, 3.02, 3.02,
- 3.05, 3.15, 3.21, 3.26, 3.33 };
- double eym[] = { 0.06324555320336758, 0.06324555320336758,
- 0.06324555320336758, 0.06324555320336758,
- 0.06324555320336758, 0.05385164807134505,
- 0.05385164807134505, 0.05385164807134505,
- 0.05385164807134505, 0.06324555320336758,
- 0.06324555320336758, 0.06324555320336758,
- 0.06324555320336758, 0.06324555320336758 };
- double eyp[] = { 0.08246211251235322, 0.08246211251235322,
- 0.08246211251235322, 0.08246211251235322,
- 0.08246211251235322, 0.08246211251235322,
- 0.07280109889280519, 0.08246211251235322,
- 0.08246211251235322, 0.08246211251235322,
- 0.08246211251235322, 0.08246211251235322,
- 0.08246211251235322, 0.08246211251235322 };
- const int np = 14;
- for (int i = 0; i < np; i++) {
- eym[i] += 0.02;
- eyp[i] += 0.02;
+ SetGraphAttributes(g, INEL, UA5, false, "ua5_inel", "UA5");
+ SetGraphAttributes(gm, INEL, UA5, true, "ua5_inel_mirrored",
+ "UA5 (mirrored)");
+
+ return (mirrored ? gm : g);
}
- g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, INEL, ALICE, false, "alice_inel",
- "ALICE INEL (publ.)");
- return g;
-}
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$ \sqrt{s}
+ * = 900GeV@f$
+ * p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel900()
+ {
+ // INEL - p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
+ TGraphAsymmErrors* g = 0;
+ double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3,
+ -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 3.28, 3.28, 3.22, 3.12, 3.06, 3.02, 2.98, 3.02, 3.02,
+ 3.05, 3.15, 3.21, 3.26, 3.33 };
+ double eym[] = { 0.06324555320336758, 0.06324555320336758,
+ 0.06324555320336758, 0.06324555320336758,
+ 0.06324555320336758, 0.05385164807134505,
+ 0.05385164807134505, 0.05385164807134505,
+ 0.05385164807134505, 0.06324555320336758,
+ 0.06324555320336758, 0.06324555320336758,
+ 0.06324555320336758, 0.06324555320336758 };
+ double eyp[] = { 0.08246211251235322, 0.08246211251235322,
+ 0.08246211251235322, 0.08246211251235322,
+ 0.08246211251235322, 0.08246211251235322,
+ 0.07280109889280519, 0.08246211251235322,
+ 0.08246211251235322, 0.08246211251235322,
+ 0.08246211251235322, 0.08246211251235322,
+ 0.08246211251235322, 0.08246211251235322 };
+ const int np = 14;
+ for (int i = 0; i < np; i++) {
+ eym[i] += 0.02;
+ eyp[i] += 0.02;
+ }
+ g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, INEL, ALICE, false, "alice_inel",
+ "#it{Eur.Phys.J}.#bf{C68}:89-108,2010"
+ /* "ALICE INEL (publ.)" */);
-//____________________________________________________________________
-/**
- * Get the ALICE INEL>0 data in @f$ |\eta|<1.3@f$ for pp at @f$
- * \sqrt{s} = 900GeV@f$
- *
- * p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralInelGt900()
-{
-#if 0
- // INEL>0 - p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
- double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7,
- 0.9 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1 };
- double y[] = { 4.0, 3.87, 3.8, 3.7, 3.67, 3.73, 3.72, 3.77, 3.92,
- 4.01 };
- double eym[] = { 0.07615773105863909, 0.07615773105863909,
- 0.07615773105863909, 0.06324555320336758,
- 0.06324555320336758, 0.06324555320336758,
- 0.0670820393249937, 0.07615773105863909,
- 0.07615773105863909, 0.07615773105863909 };
- double eyp[] = { 0.08544003745317531, 0.07615773105863909,
- 0.07615773105863909, 0.07280109889280519,
- 0.07280109889280519, 0.07280109889280519,
- 0.07615773105863909, 0.07615773105863909,
- 0.08544003745317531, 0.08544003745317531 };
- const int np = 10;
- for (int i = 0; i < np; i++) {
- double stat = (i >= 3 && i<=5) ? 0.02 : 0.03;
- eym[i] += stat;
- eyp[i] += stat;
+ return g;
}
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL>0 data in @f$ |\eta|<1.3@f$ for pp at @f$
+ * \sqrt{s} = 900GeV@f$
+ *
+ * p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt900()
+ {
+#if 0
+ // INEL>0 - p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
+ double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7,
+ 0.9 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1 };
+ double y[] = { 4.0, 3.87, 3.8, 3.7, 3.67, 3.73, 3.72, 3.77, 3.92,
+ 4.01 };
+ double eym[] = { 0.07615773105863909, 0.07615773105863909,
+ 0.07615773105863909, 0.06324555320336758,
+ 0.06324555320336758, 0.06324555320336758,
+ 0.0670820393249937, 0.07615773105863909,
+ 0.07615773105863909, 0.07615773105863909 };
+ double eyp[] = { 0.08544003745317531, 0.07615773105863909,
+ 0.07615773105863909, 0.07280109889280519,
+ 0.07280109889280519, 0.07280109889280519,
+ 0.07615773105863909, 0.07615773105863909,
+ 0.08544003745317531, 0.08544003745317531 };
+ const int np = 10;
+ for (int i = 0; i < np; i++) {
+ double stat = (i >= 3 && i<=5) ? 0.02 : 0.03;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
+
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
#else
- // These are from JFGO
- TGraphAsymmErrors *g = new TGraphAsymmErrors(15);
- g->SetPoint(0,-1.5,4.12575);
- g->SetPointError(0,0.1,0.1,0.0742967,0.0814571);
- g->SetPoint(1,-1.3,3.91209);
- g->SetPointError(1,0.1,0.1,0.0697701,0.0766199);
- g->SetPoint(2,-1.1,3.98377);
- g->SetPointError(2,0.1,0.1,0.0704503,0.0774795);
- g->SetPoint(3,-0.9,4.00035);
- g->SetPointError(3,0.1,0.1,0.0702388,0.0773433);
- g->SetPoint(4,-0.7,3.87228);
- g->SetPointError(4,0.1,0.1,0.067597,0.0745103);
- g->SetPoint(5,-0.5,3.79613);
- g->SetPointError(5,0.1,0.1,0.0659771,0.0727816);
- g->SetPoint(6,-0.3,3.70489);
- g->SetPointError(6,0.1,0.1,0.0642016,0.0708603);
- g->SetPoint(7,-0.1,3.67423);
- g->SetPointError(7,0.1,0.1,0.0635759,0.0701884);
- g->SetPoint(8,0.1,3.72765);
- g->SetPointError(8,0.1,0.1,0.0645004,0.071209);
- g->SetPoint(9,0.3,3.72171);
- g->SetPointError(9,0.1,0.1,0.064493,0.071182);
- g->SetPoint(10,0.5,3.77428);
- g->SetPointError(10,0.1,0.1,0.0655974,0.0723627);
- g->SetPoint(11,0.7,3.91704);
- g->SetPointError(11,0.1,0.1,0.0683783,0.0753716);
- g->SetPoint(12,0.9,4.00674);
- g->SetPointError(12,0.1,0.1,0.0703511,0.0774669);
- g->SetPoint(13,1.1,3.97948);
- g->SetPointError(13,0.1,0.1,0.0703744,0.077396);
- g->SetPoint(14,1.3,3.99165);
- g->SetPointError(14,0.1,0.1,0.0711888,0.078178);
+ // These are from JFGO
+ TGraphAsymmErrors *g = new TGraphAsymmErrors(15);
+ g->SetPoint(0,-1.5,4.12575);
+ g->SetPointError(0,0.1,0.1,0.0742967,0.0814571);
+ g->SetPoint(1,-1.3,3.91209);
+ g->SetPointError(1,0.1,0.1,0.0697701,0.0766199);
+ g->SetPoint(2,-1.1,3.98377);
+ g->SetPointError(2,0.1,0.1,0.0704503,0.0774795);
+ g->SetPoint(3,-0.9,4.00035);
+ g->SetPointError(3,0.1,0.1,0.0702388,0.0773433);
+ g->SetPoint(4,-0.7,3.87228);
+ g->SetPointError(4,0.1,0.1,0.067597,0.0745103);
+ g->SetPoint(5,-0.5,3.79613);
+ g->SetPointError(5,0.1,0.1,0.0659771,0.0727816);
+ g->SetPoint(6,-0.3,3.70489);
+ g->SetPointError(6,0.1,0.1,0.0642016,0.0708603);
+ g->SetPoint(7,-0.1,3.67423);
+ g->SetPointError(7,0.1,0.1,0.0635759,0.0701884);
+ g->SetPoint(8,0.1,3.72765);
+ g->SetPointError(8,0.1,0.1,0.0645004,0.071209);
+ g->SetPoint(9,0.3,3.72171);
+ g->SetPointError(9,0.1,0.1,0.064493,0.071182);
+ g->SetPoint(10,0.5,3.77428);
+ g->SetPointError(10,0.1,0.1,0.0655974,0.0723627);
+ g->SetPoint(11,0.7,3.91704);
+ g->SetPointError(11,0.1,0.1,0.0683783,0.0753716);
+ g->SetPoint(12,0.9,4.00674);
+ g->SetPointError(12,0.1,0.1,0.0703511,0.0774669);
+ g->SetPoint(13,1.1,3.97948);
+ g->SetPointError(13,0.1,0.1,0.0703744,0.077396);
+ g->SetPoint(14,1.3,3.99165);
+ g->SetPointError(14,0.1,0.1,0.0711888,0.078178);
#endif
- SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt900",
- "ALICE INEL>0 (publ.)");
- return g;
-}
+ SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt900",
+ // "ALICE INEL>0 (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:345-354,2010");
+ return g;
+ }
-//____________________________________________________________________
-/**
- * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
- * \sqrt{s} = 2.36TeV@f$
- *
- * p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralInelGt2360()
-{
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
+ * \sqrt{s} = 2.36TeV@f$
+ *
+ * p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt2360()
+ {
#if 0
- // INEL>0 - p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
- double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 4.91, 4.76, 4.63, 4.64, 4.55, 4.55, 4.64, 4.66, 4.82, 4.88 };
- double eym[] = { 0.08544003745317531, 0.08544003745317531,
- 0.08544003745317531, 0.08544003745317531,
- 0.08544003745317531, 0.08544003745317531,
- 0.08544003745317531, 0.08544003745317531,
- 0.08544003745317531, 0.08544003745317531 };
- double eyp[] = { 0.11401754250991379, 0.11401754250991379,
- 0.1044030650891055, 0.1044030650891055,
- 0.1044030650891055, 0.1044030650891055,
- 0.1044030650891055, 0.1044030650891055,
- 0.11401754250991379, 0.11401754250991379 };
- const int np = 10;
- for (int i = 0; i < np; i++) {
- double stat = 0.3;
- eym[i] += stat;
- eyp[i] += stat;
- }
+ // INEL>0 - p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
+ double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 4.91, 4.76, 4.63, 4.64, 4.55, 4.55, 4.64, 4.66, 4.82, 4.88 };
+ double eym[] = { 0.08544003745317531, 0.08544003745317531,
+ 0.08544003745317531, 0.08544003745317531,
+ 0.08544003745317531, 0.08544003745317531,
+ 0.08544003745317531, 0.08544003745317531,
+ 0.08544003745317531, 0.08544003745317531 };
+ double eyp[] = { 0.11401754250991379, 0.11401754250991379,
+ 0.1044030650891055, 0.1044030650891055,
+ 0.1044030650891055, 0.1044030650891055,
+ 0.1044030650891055, 0.1044030650891055,
+ 0.11401754250991379, 0.11401754250991379 };
+ const int np = 10;
+ for (int i = 0; i < np; i++) {
+ double stat = 0.3;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
#else
- // These are from JFGO
- TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
- g->SetPoint(0,-1.5,4.79047);
- g->SetPointError(0,0.1,0.1,0.0844278,0.109947);
- g->SetPoint(1,-1.3,4.91068);
- g->SetPointError(1,0.1,0.1,0.0856751,0.112038);
- g->SetPoint(2,-1.1,4.87386);
- g->SetPointError(2,0.1,0.1,0.0842846,0.110628);
- g->SetPoint(3,-0.9,4.91365);
- g->SetPointError(3,0.1,0.1,0.084339,0.111049);
- g->SetPoint(4,-0.7,4.7601);
- g->SetPointError(4,0.1,0.1,0.0812087,0.107203);
- g->SetPoint(5,-0.5,4.63355);
- g->SetPointError(5,0.1,0.1,0.078687,0.104079);
- g->SetPoint(6,-0.3,4.63885);
- g->SetPointError(6,0.1,0.1,0.0785337,0.104014);
- g->SetPoint(7,-0.1,4.55439);
- g->SetPointError(7,0.1,0.1,0.0769842,0.10203);
- g->SetPoint(8,0.1,4.55087);
- g->SetPointError(8,0.1,0.1,0.0769246,0.101951);
- g->SetPoint(9,0.3,4.64118);
- g->SetPointError(9,0.1,0.1,0.0785732,0.104066);
- g->SetPoint(10,0.5,4.66172);
- g->SetPointError(10,0.1,0.1,0.0791652,0.104711);
- g->SetPoint(11,0.7,4.81871);
- g->SetPointError(11,0.1,0.1,0.0822086,0.108523);
- g->SetPoint(12,0.9,4.88193);
- g->SetPointError(12,0.1,0.1,0.0837944,0.110332);
- g->SetPoint(13,1.1,4.89068);
- g->SetPointError(13,0.1,0.1,0.0845754,0.111009);
- g->SetPoint(14,1.3,5.05663);
- g->SetPointError(14,0.1,0.1,0.0882216,0.115368);
+ // These are from JFGO
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
+ g->SetPoint(0,-1.5,4.79047);
+ g->SetPointError(0,0.1,0.1,0.0844278,0.109947);
+ g->SetPoint(1,-1.3,4.91068);
+ g->SetPointError(1,0.1,0.1,0.0856751,0.112038);
+ g->SetPoint(2,-1.1,4.87386);
+ g->SetPointError(2,0.1,0.1,0.0842846,0.110628);
+ g->SetPoint(3,-0.9,4.91365);
+ g->SetPointError(3,0.1,0.1,0.084339,0.111049);
+ g->SetPoint(4,-0.7,4.7601);
+ g->SetPointError(4,0.1,0.1,0.0812087,0.107203);
+ g->SetPoint(5,-0.5,4.63355);
+ g->SetPointError(5,0.1,0.1,0.078687,0.104079);
+ g->SetPoint(6,-0.3,4.63885);
+ g->SetPointError(6,0.1,0.1,0.0785337,0.104014);
+ g->SetPoint(7,-0.1,4.55439);
+ g->SetPointError(7,0.1,0.1,0.0769842,0.10203);
+ g->SetPoint(8,0.1,4.55087);
+ g->SetPointError(8,0.1,0.1,0.0769246,0.101951);
+ g->SetPoint(9,0.3,4.64118);
+ g->SetPointError(9,0.1,0.1,0.0785732,0.104066);
+ g->SetPoint(10,0.5,4.66172);
+ g->SetPointError(10,0.1,0.1,0.0791652,0.104711);
+ g->SetPoint(11,0.7,4.81871);
+ g->SetPointError(11,0.1,0.1,0.0822086,0.108523);
+ g->SetPoint(12,0.9,4.88193);
+ g->SetPointError(12,0.1,0.1,0.0837944,0.110332);
+ g->SetPoint(13,1.1,4.89068);
+ g->SetPointError(13,0.1,0.1,0.0845754,0.111009);
+ g->SetPoint(14,1.3,5.05663);
+ g->SetPointError(14,0.1,0.1,0.0882216,0.115368);
#endif
- SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt2360",
- "ALICE INEL>0 (publ.)");
- return g;
-}
+ SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt2360",
+ // "ALICE INEL>0 (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:345-354,2010");
+ return g;
+ }
-//____________________________________________________________________
-/**
- * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
- * \sqrt{s} = 7TeV@f$
- *
- * p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralInelGt7000()
-{
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
+ * \sqrt{s} = 7TeV@f$
+ *
+ * p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt7000()
+ {
#if 0
- // INEL>0 - p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
-// Plot: p7741_d6x1y1
- double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 6.22, 6.07, 6.01, 5.84, 5.85, 5.85, 5.91, 6.01, 6.17, 6.26 };
- double eym[] = { 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.11180339887498948,
- 0.11180339887498948, 0.11180339887498948,
- 0.11180339887498948, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644 };
- double eyp[] = { 0.21095023109728983, 0.21095023109728983,
- 0.2009975124224178, 0.2009975124224178,
- 0.2009975124224178, 0.2009975124224178,
- 0.2009975124224178, 0.2009975124224178,
- 0.21095023109728983, 0.21095023109728983 };
- const int np = 10;
- for (int i = 0; i < np; i++) {
- double stat = 0.2;
- eym[i] += stat;
- eyp[i] += stat;
- }
+ // INEL>0 - p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
+ // Plot: p7741_d6x1y1
+ double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 6.22, 6.07, 6.01, 5.84, 5.85, 5.85, 5.91, 6.01, 6.17, 6.26 };
+ double eym[] = { 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.11180339887498948,
+ 0.11180339887498948, 0.11180339887498948,
+ 0.11180339887498948, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644 };
+ double eyp[] = { 0.21095023109728983, 0.21095023109728983,
+ 0.2009975124224178, 0.2009975124224178,
+ 0.2009975124224178, 0.2009975124224178,
+ 0.2009975124224178, 0.2009975124224178,
+ 0.21095023109728983, 0.21095023109728983 };
+ const int np = 10;
+ for (int i = 0; i < np; i++) {
+ double stat = 0.2;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
#else
- // These are from JFGO
- TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
- g->SetPoint(0,-1.5,6.28573);
- g->SetPointError(0,0.1,0.1,0.125928,0.215392);
- g->SetPoint(1,-1.3,6.25573);
- g->SetPointError(1,0.1,0.1,0.124352,0.213795);
- g->SetPoint(2,-1.1,6.28779);
- g->SetPointError(2,0.1,0.1,0.124143,0.214399);
- g->SetPoint(3,-0.9,6.21881);
- g->SetPointError(3,0.1,0.1,0.122079,0.211642);
- g->SetPoint(4,-0.7,6.0728);
- g->SetPointError(4,0.1,0.1,0.118661,0.206355);
- g->SetPoint(5,-0.5,6.011);
- g->SetPointError(5,0.1,0.1,0.117043,0.204019);
- g->SetPoint(6,-0.3,5.84071);
- g->SetPointError(6,0.1,0.1,0.11346,0.198086);
- g->SetPoint(7,-0.1,5.8532);
- g->SetPointError(7,0.1,0.1,0.113569,0.198433);
- g->SetPoint(8,0.1,5.84811);
- g->SetPointError(8,0.1,0.1,0.11347,0.198261);
- g->SetPoint(9,0.3,5.91022);
- g->SetPointError(9,0.1,0.1,0.11481,0.200444);
- g->SetPoint(10,0.5,6.00649);
- g->SetPointError(10,0.1,0.1,0.116955,0.203866);
- g->SetPoint(11,0.7,6.17115);
- g->SetPointError(11,0.1,0.1,0.120583,0.209697);
- g->SetPoint(12,0.9,6.2645);
- g->SetPointError(12,0.1,0.1,0.122976,0.213197);
- g->SetPoint(13,1.1,6.36448);
- g->SetPointError(13,0.1,0.1,0.125657,0.217014);
- g->SetPoint(14,1.3,6.39489);
- g->SetPointError(14,0.1,0.1,0.127118,0.218551);
+ // These are from JFGO
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
+ g->SetPoint(0,-1.5,6.28573);
+ g->SetPointError(0,0.1,0.1,0.125928,0.215392);
+ g->SetPoint(1,-1.3,6.25573);
+ g->SetPointError(1,0.1,0.1,0.124352,0.213795);
+ g->SetPoint(2,-1.1,6.28779);
+ g->SetPointError(2,0.1,0.1,0.124143,0.214399);
+ g->SetPoint(3,-0.9,6.21881);
+ g->SetPointError(3,0.1,0.1,0.122079,0.211642);
+ g->SetPoint(4,-0.7,6.0728);
+ g->SetPointError(4,0.1,0.1,0.118661,0.206355);
+ g->SetPoint(5,-0.5,6.011);
+ g->SetPointError(5,0.1,0.1,0.117043,0.204019);
+ g->SetPoint(6,-0.3,5.84071);
+ g->SetPointError(6,0.1,0.1,0.11346,0.198086);
+ g->SetPoint(7,-0.1,5.8532);
+ g->SetPointError(7,0.1,0.1,0.113569,0.198433);
+ g->SetPoint(8,0.1,5.84811);
+ g->SetPointError(8,0.1,0.1,0.11347,0.198261);
+ g->SetPoint(9,0.3,5.91022);
+ g->SetPointError(9,0.1,0.1,0.11481,0.200444);
+ g->SetPoint(10,0.5,6.00649);
+ g->SetPointError(10,0.1,0.1,0.116955,0.203866);
+ g->SetPoint(11,0.7,6.17115);
+ g->SetPointError(11,0.1,0.1,0.120583,0.209697);
+ g->SetPoint(12,0.9,6.2645);
+ g->SetPointError(12,0.1,0.1,0.122976,0.213197);
+ g->SetPoint(13,1.1,6.36448);
+ g->SetPointError(13,0.1,0.1,0.125657,0.217014);
+ g->SetPoint(14,1.3,6.39489);
+ g->SetPointError(14,0.1,0.1,0.127118,0.218551);
#endif
- SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt7000",
- "ALICE INEL>0 (publ.)");
- return g;
-}
+ SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt7000",
+ // "ALICE INEL>0 (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:345-354,2010");
+ return g;
+ }
-//____________________________________________________________________
-/**
- * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
- * \sqrt{s} = 900GeV@f$
- *
- * p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralNsd900()
-{
- //NSD - p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
- double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
- 0.5, 0.7, 0.9, 1.1, 1.3 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 3.9, 3.89, 3.81, 3.7, 3.64, 3.59, 3.53, 3.58, 3.59,
- 3.61, 3.74, 3.8, 3.87, 3.95 };
- double eym[] = { 0.13341664064126335, 0.13152946437965907,
- 0.13152946437965907, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.13152946437965907,
- 0.13152946437965907, 0.13341664064126335 };
- double eyp[] = { 0.13341664064126335, 0.13152946437965907,
- 0.13152946437965907, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.13152946437965907,
- 0.13152946437965907, 0.13341664064126335 };
- const int np = 14;
- for (int i = 0; i < np; i++) {
- double stat = (i == 0 || i == np-1) ? 0.03 : 0.02;
- eym[i] += stat;
- eyp[i] += stat;
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
+ * \sqrt{s} = 900GeV@f$
+ *
+ * p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd900()
+ {
+ //NSD - p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
+ double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
+ 0.5, 0.7, 0.9, 1.1, 1.3 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 3.9, 3.89, 3.81, 3.7, 3.64, 3.59, 3.53, 3.58, 3.59,
+ 3.61, 3.74, 3.8, 3.87, 3.95 };
+ double eym[] = { 0.13341664064126335, 0.13152946437965907,
+ 0.13152946437965907, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.13152946437965907,
+ 0.13152946437965907, 0.13341664064126335 };
+ double eyp[] = { 0.13341664064126335, 0.13152946437965907,
+ 0.13152946437965907, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.13152946437965907,
+ 0.13152946437965907, 0.13341664064126335 };
+ const int np = 14;
+ for (int i = 0; i < np; i++) {
+ double stat = (i == 0 || i == np-1) ? 0.03 : 0.02;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
+
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd",
+ //"ALICE NSD (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:89-108,2010");
+
+ return g;
}
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd", "ALICE NSD (publ.)");
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$
+ * \sqrt{s} = 2.36TeV@f$
+ *
+ * p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel2360()
+ {
+ // INEL - p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
+ double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
+ 0.5, 0.7, 0.9, 1.1, 1.3 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 4.08, 4.01, 4.0, 3.88, 3.77, 3.8, 3.73, 3.71, 3.79,
+ 3.82, 3.94, 4.02, 4.05, 4.16 };
+ double eym[] = { 0.13341664064126335, 0.12369316876852982,
+ 0.12369316876852982, 0.1216552506059644,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.1216552506059644, 0.11180339887498948,
+ 0.1216552506059644, 0.1216552506059644,
+ 0.12369316876852982, 0.12369316876852982,
+ 0.13341664064126335, 0.13341664064126335 };
+ double eyp[] = { 0.2716615541441225, 0.2716615541441225,
+ 0.2716615541441225, 0.260768096208106,
+ 0.25079872407968906, 0.25079872407968906,
+ 0.25079872407968906, 0.25079872407968906,
+ 0.25079872407968906, 0.260768096208106,
+ 0.261725046566048, 0.2716615541441225,
+ 0.2716615541441225, 0.2816025568065745 };
+ const int np = 14;
+ for (int i = 0; i < np; i++) {
+ double stat = (i < 3 || i > np-1-4) ? 0.03 : 0.02;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
- return g;
-}
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, ALICE, false, "alice_inel2360",
+ // "ALICE INEL (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:89-108,2010");
+ return g;
+ }
-//____________________________________________________________________
-/**
- * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$
- * \sqrt{s} = 2.36TeV@f$
- *
- * p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralInel2360()
-{
- // INEL - p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
- double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
- 0.5, 0.7, 0.9, 1.1, 1.3 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 4.08, 4.01, 4.0, 3.88, 3.77, 3.8, 3.73, 3.71, 3.79,
- 3.82, 3.94, 4.02, 4.05, 4.16 };
- double eym[] = { 0.13341664064126335, 0.12369316876852982,
- 0.12369316876852982, 0.1216552506059644,
- 0.1216552506059644, 0.1216552506059644,
- 0.1216552506059644, 0.11180339887498948,
- 0.1216552506059644, 0.1216552506059644,
- 0.12369316876852982, 0.12369316876852982,
- 0.13341664064126335, 0.13341664064126335 };
- double eyp[] = { 0.2716615541441225, 0.2716615541441225,
- 0.2716615541441225, 0.260768096208106,
- 0.25079872407968906, 0.25079872407968906,
- 0.25079872407968906, 0.25079872407968906,
- 0.25079872407968906, 0.260768096208106,
- 0.261725046566048, 0.2716615541441225,
- 0.2716615541441225, 0.2816025568065745 };
- const int np = 14;
- for (int i = 0; i < np; i++) {
- double stat = (i < 3 || i > np-1-4) ? 0.03 : 0.02;
- eym[i] += stat;
- eyp[i] += stat;
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
+ * \sqrt{s} = 2.36TeV@f$
+ *
+ * p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd2360()
+ {
+ // NSD - p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
+ double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
+ 0.5, 0.7, 0.9, 1.1, 1.3 };
+ double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
+ 0.1, 0.1, 0.1, 0.1, 0.1 };
+ double y[] = { 4.79, 4.72, 4.7, 4.56, 4.44, 4.47, 4.39, 4.37, 4.47,
+ 4.5, 4.64, 4.73, 4.76, 4.9 };
+ double eym[] = { 0.13601470508735444, 0.13341664064126335,
+ 0.13341664064126335, 0.12369316876852982,
+ 0.12369316876852982, 0.12369316876852982,
+ 0.12369316876852982, 0.12369316876852982,
+ 0.12369316876852982, 0.12369316876852982,
+ 0.12369316876852982, 0.13341664064126335,
+ 0.13341664064126335, 0.13341664064126335 };
+ double eyp[] = { 0.18439088914585774, 0.18248287590894657,
+ 0.18248287590894657, 0.1726267650163207,
+ 0.1726267650163207, 0.1726267650163207,
+ 0.16278820596099708, 0.16278820596099708,
+ 0.1726267650163207, 0.1726267650163207,
+ 0.1726267650163207, 0.18248287590894657,
+ 0.18248287590894657, 0.18248287590894657 };
+ const int np = 14;
+
+ for (int i = 0; i < np; i++) {
+ double stat = (i < 1) ? 0.03 : 0.02;
+ eym[i] += stat;
+ eyp[i] += stat;
+ }
+
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd2360",
+ // "ALICE NSD (publ.)"
+ "#it{Eur.Phys.J}.#bf{C68}:89-108,2010");
+ return g;
}
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, ALICE, false, "alice_inel2360",
- "ALICE INEL (publ.)");
- return g;
-}
+ //____________________________________________________________________
+ /**
+ * Get the ALICE MB (CINT5, V0AND) data in @f$ |\eta|<2@f$ for pPb at
+ * @f$ \sqrt{s} = 5.023TeV@f$
+ *
+ * arXiv:1210.3615 [nucl-ex]
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralpPb5023()
+ {
+ TGraphAsymmErrors *g = new TGraphAsymmErrors(43);
+ g->SetName("");
+ g->SetTitle("");
+ g->SetFillColor(1);
+ g->SetMarkerStyle(24);
+ g->SetMarkerSize(1.3);
+ g->SetPoint(0,-2.05,16.1757);
+ g->SetPointError(0,0.05,0.05,0.360241,0.360241);
+ g->SetPoint(1,-1.95,16.8454);
+ g->SetPointError(1,0.05,0.05,0.184436,0.184436);
+ g->SetPoint(2,-1.85,16.9748);
+ g->SetPointError(2,0.05,0.05,0.112833,0.112833);
+ g->SetPoint(3,-1.75,17.2504);
+ g->SetPointError(3,0.05,0.05,0.0851792,0.0851792);
+ g->SetPoint(4,-1.65,17.336);
+ g->SetPointError(4,0.05,0.05,0.0674453,0.0674453);
+ g->SetPoint(5,-1.55,17.4649);
+ g->SetPointError(5,0.05,0.05,0.0558537,0.0558537);
+ g->SetPoint(6,-1.45,17.5884);
+ g->SetPointError(6,0.05,0.05,0.0481065,0.0481065);
+ g->SetPoint(7,-1.35,17.7509);
+ g->SetPointError(7,0.05,0.05,0.0448122,0.0448122);
+ g->SetPoint(8,-1.25,17.8052);
+ g->SetPointError(8,0.05,0.05,0.0403298,0.0403298);
+ g->SetPoint(9,-1.15,17.8357);
+ g->SetPointError(9,0.05,0.05,0.0383174,0.0383174);
+ g->SetPoint(10,-1.05,17.7547);
+ g->SetPointError(10,0.05,0.05,0.0356689,0.0356689);
+ g->SetPoint(11,-0.95,17.6859);
+ g->SetPointError(11,0.05,0.05,0.034326,0.034326);
+ g->SetPoint(12,-0.85,17.6665);
+ g->SetPointError(12,0.05,0.05,0.0333449,0.0333449);
+ g->SetPoint(13,-0.75,17.6044);
+ g->SetPointError(13,0.05,0.05,0.0325356,0.0325356);
+ g->SetPoint(14,-0.65,17.4815);
+ g->SetPointError(14,0.05,0.05,0.0318004,0.0318004);
+ g->SetPoint(15,-0.55,17.4);
+ g->SetPointError(15,0.05,0.05,0.0312675,0.0312675);
+ g->SetPoint(16,-0.45,17.3425);
+ g->SetPointError(16,0.05,0.05,0.0310344,0.0310344);
+ g->SetPoint(17,-0.35,17.2885);
+ g->SetPointError(17,0.05,0.05,0.0306043,0.0306043);
+ g->SetPoint(18,-0.25,17.2646);
+ g->SetPointError(18,0.05,0.05,0.0303226,0.0303226);
+ g->SetPoint(19,-0.15,17.316);
+ g->SetPointError(19,0.05,0.05,0.0302368,0.0302368);
+ g->SetPoint(20,-0.05,17.312);
+ g->SetPointError(20,0.05,0.05,0.0301444,0.0301444);
+ g->SetPoint(21,0.05,17.4418);
+ g->SetPointError(21,0.05,0.05,0.0301526,0.0301526);
+ g->SetPoint(22,0.15,17.4944);
+ g->SetPointError(22,0.05,0.05,0.0303199,0.0303199);
+ g->SetPoint(23,0.25,17.642);
+ g->SetPointError(23,0.05,0.05,0.0303867,0.0303867);
+ g->SetPoint(24,0.35,17.8153);
+ g->SetPointError(24,0.05,0.05,0.0306752,0.0306752);
+ g->SetPoint(25,0.45,18.0244);
+ g->SetPointError(25,0.05,0.05,0.0310274,0.0310274);
+ g->SetPoint(26,0.55,18.1993);
+ g->SetPointError(26,0.05,0.05,0.0314353,0.0314353);
+ g->SetPoint(27,0.65,18.349);
+ g->SetPointError(27,0.05,0.05,0.0316803,0.0316803);
+ g->SetPoint(28,0.75,18.5976);
+ g->SetPointError(28,0.05,0.05,0.0322819,0.0322819);
+ g->SetPoint(29,0.85,18.8045);
+ g->SetPointError(29,0.05,0.05,0.0329447,0.0329447);
+ g->SetPoint(30,0.95,18.9865);
+ g->SetPointError(30,0.05,0.05,0.0337513,0.0337513);
+ g->SetPoint(31,1.05,19.2313);
+ g->SetPointError(31,0.05,0.05,0.0354009,0.0354009);
+ g->SetPoint(32,1.15,19.4055);
+ g->SetPointError(32,0.05,0.05,0.0367366,0.0367366);
+ g->SetPoint(33,1.25,19.5893);
+ g->SetPointError(33,0.05,0.05,0.0385048,0.0385048);
+ g->SetPoint(34,1.35,19.8196);
+ g->SetPointError(34,0.05,0.05,0.0421699,0.0421699);
+ g->SetPoint(35,1.45,19.9476);
+ g->SetPointError(35,0.05,0.05,0.0451541,0.0451541);
+ g->SetPoint(36,1.55,20.1012);
+ g->SetPointError(36,0.05,0.05,0.0513641,0.0513641);
+ g->SetPoint(37,1.65,20.1082);
+ g->SetPointError(37,0.05,0.05,0.060302,0.060302);
+ g->SetPoint(38,1.75,20.1732);
+ g->SetPointError(38,0.05,0.05,0.0739969,0.0739969);
+ g->SetPoint(39,1.85,20.1964);
+ g->SetPointError(39,0.05,0.05,0.0953757,0.0953757);
+ g->SetPoint(40,1.95,20.0509);
+ g->SetPointError(40,0.05,0.05,0.147212,0.147212);
+ g->SetPoint(41,2.05,20.3151);
+ g->SetPointError(41,0.05,0.05,0.272151,0.272151);
+ g->SetPoint(42,2.15,20.1319);
+ g->SetPointError(42,0.05,0.05,0.802706,0.802706);
-//____________________________________________________________________
-/**
- * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
- * \sqrt{s} = 2.36TeV@f$
- *
- * p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* AliceCentralNsd2360()
-{
- // NSD - p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
- double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
- 0.5, 0.7, 0.9, 1.1, 1.3 };
- double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
- 0.1, 0.1, 0.1, 0.1, 0.1 };
- double y[] = { 4.79, 4.72, 4.7, 4.56, 4.44, 4.47, 4.39, 4.37, 4.47,
- 4.5, 4.64, 4.73, 4.76, 4.9 };
- double eym[] = { 0.13601470508735444, 0.13341664064126335,
- 0.13341664064126335, 0.12369316876852982,
- 0.12369316876852982, 0.12369316876852982,
- 0.12369316876852982, 0.12369316876852982,
- 0.12369316876852982, 0.12369316876852982,
- 0.12369316876852982, 0.13341664064126335,
- 0.13341664064126335, 0.13341664064126335 };
- double eyp[] = { 0.18439088914585774, 0.18248287590894657,
- 0.18248287590894657, 0.1726267650163207,
- 0.1726267650163207, 0.1726267650163207,
- 0.16278820596099708, 0.16278820596099708,
- 0.1726267650163207, 0.1726267650163207,
- 0.1726267650163207, 0.18248287590894657,
- 0.18248287590894657, 0.18248287590894657 };
- const int np = 14;
+
+ SetGraphAttributes(g, NSD, ALICE, false, "alice_ppb50230",
+ "arXiv:1210.3615");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
+ * = 900GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel900Work()
+ {
+ TGraphAsymmErrors *g = new TGraphAsymmErrors(18);
+ // g->SetPoint(0,-1.9,0); g->SetPointError(0,0.1,0.1,0,0);
+ g->SetPoint(0,-1.7,3.13935);g->SetPointError(0,0.1,0.1,0.0726186,0.0525276);
+ g->SetPoint(1,-1.5,3.15634);g->SetPointError(1,0.1,0.1,0.0338547,0.0380273);
+ g->SetPoint(2,-1.3,3.13683);g->SetPointError(2,0.1,0.1,0.0295176,0.0295638);
+ g->SetPoint(3,-1.1,3.10618);g->SetPointError(3,0.1,0.1,0.0306925,0.0329387);
+ g->SetPoint(4,-0.9,3.05921);g->SetPointError(4,0.1,0.1,0.0224684,0.025408);
+ g->SetPoint(5,-0.7,3.00303);g->SetPointError(5,0.1,0.1,0.0389278,0.0238328);
+ g->SetPoint(6,-0.5,2.94604);g->SetPointError(6,0.1,0.1,0.0211986,0.0322219);
+ g->SetPoint(7,-0.3,2.91507);g->SetPointError(7,0.1,0.1,0.030029,0.0209573);
+ g->SetPoint(8,-0.1,2.88965);g->SetPointError(8,0.1,0.1,0.0286516,0.0253694);
+ g->SetPoint(9,0.1,2.89731); g->SetPointError(9,0.1,0.1,0.0334615,0.0192116);
+ g->SetPoint(10,0.3,2.91188);g->SetPointError(10,0.1,0.1,0.0503868,0.024911);
+ g->SetPoint(11,0.5,2.96295);g->SetPointError(11,0.1,0.1,0.030009,0.0284692);
+ g->SetPoint(12,0.7,3.0089); g->SetPointError(12,0.1,0.1,0.0189095,0.026319);
+ g->SetPoint(13,0.9,3.07028);g->SetPointError(13,0.1,0.1,0.0449128,0.030738);
+ g->SetPoint(14,1.1,3.10215);g->SetPointError(14,0.1,0.1,0.0288688,0.026301);
+ g->SetPoint(15,1.3,3.12946);g->SetPointError(15,0.1,0.1,0.0431495,0.026355);
+ g->SetPoint(16,1.5,3.14549);g->SetPointError(16,0.1,0.1,0.0322482,0.033611);
+ g->SetPoint(17,1.7,3.15729);g->SetPointError(17,0.1,0.1,0.105509,0.0523796);
+ // g->SetPoint(19,1.9,0); g->SetPointError(19,0.1,0.1,0,0);
+
+ SetGraphAttributes(g, INEL, WIP, false, "alice_pp900work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
+ * = 900GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd900Work()
+ {
+ TGraphAsymmErrors *g = new TGraphAsymmErrors(18);
- for (int i = 0; i < np; i++) {
- double stat = (i < 1) ? 0.03 : 0.02;
- eym[i] += stat;
- eyp[i] += stat;
+ g->SetPoint(0,-1.7,3.84726);g->SetPointError(0,0.1,0.1,0.114853,0.118974);
+ g->SetPoint(1,-1.5,3.87094);g->SetPointError(1,0.1,0.1,0.10574,0.108613);
+ g->SetPoint(2,-1.3,3.84769);g->SetPointError(2,0.1,0.1,0.105942,0.107644);
+ g->SetPoint(3,-1.1,3.8122); g->SetPointError(3,0.1,0.1,0.100838,0.101818);
+ g->SetPoint(4,-0.9,3.75388);g->SetPointError(4,0.1,0.1,0.0967073,0.0972099);
+ g->SetPoint(5,-0.7,3.68733);g->SetPointError(5,0.1,0.1,0.0923424,0.0925662);
+ g->SetPoint(6,-0.5,3.61874);g->SetPointError(6,0.1,0.1,0.0904027,0.090477);
+ g->SetPoint(7,-0.3,3.58091);g->SetPointError(7,0.1,0.1,0.0875915,0.0875992);
+ g->SetPoint(8,-0.1,3.54905);g->SetPointError(8,0.1,0.1,0.086046,0.0860293);
+ g->SetPoint(9,0.1,3.55968); g->SetPointError(9,0.1,0.1,0.0884776,0.0884612);
+ g->SetPoint(10,0.3,3.57729);g->SetPointError(10,0.1,0.1,0.0857614,0.085769);
+ g->SetPoint(11,0.5,3.63879);g->SetPointError(11,0.1,0.1,0.0879787,0.088056);
+ g->SetPoint(12,0.7,3.69422);g->SetPointError(12,0.1,0.1,0.0931736,0.093396);
+ g->SetPoint(13,0.9,3.76835);g->SetPointError(13,0.1,0.1,0.0928833,0.093410);
+ g->SetPoint(14,1.1,3.80647);g->SetPointError(14,0.1,0.1,0.0998124,0.100799);
+ g->SetPoint(15,1.3,3.83824);g->SetPointError(15,0.1,0.1,0.102549,0.104298);
+ g->SetPoint(16,1.5,3.85778);g->SetPointError(16,0.1,0.1,0.10267,0.105607);
+ g->SetPoint(17,1.7,3.87075);g->SetPointError(17,0.1,0.1,0.114093,0.11829);
+
+ SetGraphAttributes(g, NSD, WIP, false, "alice_pp900nsdwork",
+ "PWG-UD/MULT - work in progress");
+ return g;
}
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd2360",
- "ALICE NSD (publ.)");
- return g;
-}
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL>0 data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
+ * = 900GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt900Work()
+ {
+ TGraphAsymmErrors *g = new TGraphAsymmErrors(10);
+
+ g->SetPoint(0,-0.9,3.90755);g->SetPointError(0,0.1,0.1,0.0377085,0.0355101);
+ g->SetPoint(1,-0.7,3.8357); g->SetPointError(1,0.1,0.1,0.0477674,0.0349303);
+ g->SetPoint(2,-0.5,3.76291);g->SetPointError(2,0.1,0.1,0.0277709,0.040401);
+ g->SetPoint(3,-0.3,3.72336);g->SetPointError(3,0.1,0.1,0.0343553,0.0250805);
+ g->SetPoint(4,-0.1,3.69098);g->SetPointError(4,0.1,0.1,0.0324842,0.0324248);
+ g->SetPoint(5,0.1,3.70076); g->SetPointError(5,0.1,0.1,0.0390932,0.0246738);
+ g->SetPoint(6,0.3,3.71924); g->SetPointError(6,0.1,0.1,0.0576054,0.0287106);
+ g->SetPoint(7,0.5,3.7844); g->SetPointError(7,0.1,0.1,0.0316759,0.0295124);
+ g->SetPoint(8,0.7,3.84319); g->SetPointError(8,0.1,0.1,0.0293134,0.0332125);
+ g->SetPoint(9,0.9,3.92163); g->SetPointError(9,0.1,0.1,0.0558339,0.0394925);
+ SetGraphAttributes(g, INELGt0, WIP, false, "alice_pp900inelgtwork",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 2760GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel2760Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,4.0313); g->SetPointError(0,0.1,0.1,0.100951,0.104678);
+ g->SetPoint(1,-1.5,4.0431); g->SetPointError(1,0.1,0.1,0.10129,0.103883);
+ g->SetPoint(2,-1.3,4.01251); g->SetPointError(2,0.1,0.1,0.10847,0.110089);
+ g->SetPoint(3,-1.1,3.96799); g->SetPointError(3,0.1,0.1,0.105543,0.106623);
+ g->SetPoint(4,-0.9,3.89669); g->SetPointError(4,0.1,0.1,0.110974,0.111625);
+ g->SetPoint(5,-0.7,3.81051); g->SetPointError(5,0.1,0.1,0.108463,0.108882);
+ g->SetPoint(6,-0.5,3.76537); g->SetPointError(6,0.1,0.1,0.105488,0.105773);
+ g->SetPoint(7,-0.3,3.69733); g->SetPointError(7,0.1,0.1,0.110156,0.11035);
+ g->SetPoint(8,-0.1,3.68148); g->SetPointError(8,0.1,0.1,0.105564,0.105733);
+ g->SetPoint(9,0.1,3.67386); g->SetPointError(9,0.1,0.1,0.1058,0.105968);
+ g->SetPoint(10,0.3,3.69873); g->SetPointError(10,0.1,0.1,0.107167,0.107367);
+ g->SetPoint(11,0.5,3.76377); g->SetPointError(11,0.1,0.1,0.111177,0.111448);
+ g->SetPoint(12,0.7,3.81956); g->SetPointError(12,0.1,0.1,0.107198,0.107623);
+ g->SetPoint(13,0.9,3.89506); g->SetPointError(13,0.1,0.1,0.105617,0.1063);
+ g->SetPoint(14,1.1,3.95888); g->SetPointError(14,0.1,0.1,0.111316,0.112336);
+ g->SetPoint(15,1.3,4.00176); g->SetPointError(15,0.1,0.1,0.111751,0.113315);
+ g->SetPoint(16,1.5,4.03247); g->SetPointError(16,0.1,0.1,0.114383,0.116674);
+ g->SetPoint(17,1.7,4.061); g->SetPointError(17,0.1,0.1,0.107094,0.110665);
-//____________________________________________________________________
-/**
- * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
- * \sqrt{s} = 900GeV@f$
- *
- * p7743_d8x1y1
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* CMSNsd900()
-{
- // CMS published NSD data - p7743_d8x1y1
- double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75, 1.25, 1.75,
- 2.25 };
- double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double y[] = { 3.6, 3.73, 3.62, 3.54, 3.48, 3.48, 3.54, 3.62, 3.73, 3.6 };
- double eym[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14,0.13 };
- double eyp[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14, 0.13 };
- const int np = 10;
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, CMS, false, "cms_nsd900", "CMS NSD");
+ SetGraphAttributes(g, INEL, WIP, false,
+ "alice_ppInel2760Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 2760GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd2760Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,4.8704); g->SetPointError(0,0.1,0.1,0.221293,0.224755);
+ g->SetPoint(1,-1.5,4.88859); g->SetPointError(1,0.1,0.1,0.221269,0.223478);
+ g->SetPoint(2,-1.3,4.85326); g->SetPointError(2,0.1,0.1,0.21109,0.212455);
+ g->SetPoint(3,-1.1,4.80085); g->SetPointError(3,0.1,0.1,0.21041,0.211157);
+ g->SetPoint(4,-0.9,4.71513); g->SetPointError(4,0.1,0.1,0.198361,0.198749);
+ g->SetPoint(5,-0.7,4.61153); g->SetPointError(5,0.1,0.1,0.194009,0.194176);
+ g->SetPoint(6,-0.5,4.55715); g->SetPointError(6,0.1,0.1,0.193226,0.193281);
+ g->SetPoint(7,-0.3,4.47508); g->SetPointError(7,0.1,0.1,0.182433,0.182439);
+ g->SetPoint(8,-0.1,4.45709); g->SetPointError(8,0.1,0.1,0.186518,0.186506);
+ g->SetPoint(9,0.1,4.44707); g->SetPointError(9,0.1,0.1,0.185747,0.185735);
+ g->SetPoint(10,0.3,4.47734); g->SetPointError(10,0.1,0.1,0.185835,0.185841);
+ g->SetPoint(11,0.5,4.55477); g->SetPointError(11,0.1,0.1,0.186934,0.186991);
+ g->SetPoint(12,0.7,4.62236); g->SetPointError(12,0.1,0.1,0.196631,0.196796);
+ g->SetPoint(13,0.9,4.71277); g->SetPointError(13,0.1,0.1,0.204034,0.20441);
+ g->SetPoint(14,1.1,4.78902); g->SetPointError(14,0.1,0.1,0.20317,0.20394);
+ g->SetPoint(15,1.3,4.84008); g->SetPointError(15,0.1,0.1,0.205573,0.206967);
+ g->SetPoint(16,1.5,4.87453); g->SetPointError(16,0.1,0.1,0.206314,0.208667);
+ g->SetPoint(17,1.7,4.90614); g->SetPointError(17,0.1,0.1,0.218996,0.222545);
+
+ SetGraphAttributes(g, NSD, WIP, false,
+ "alice_ppNsd2760Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 2760GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt2760Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
+ g->SetPoint(0,-0.9,4.96315); g->SetPointError(0,0.1,0.1,0.0439746,0.0440108);
+ g->SetPoint(1,-0.7,4.8532); g->SetPointError(1,0.1,0.1,0.0426373,0.0600727);
+ g->SetPoint(2,-0.5,4.79582);g->SetPointError(2,0.1,0.1,0.0475367,0.0466255);
+ g->SetPoint(3,-0.3,4.70907);g->SetPointError(3,0.1,0.1,0.0313084,0.0468084);
+ g->SetPoint(4,-0.1,4.68906);g->SetPointError(4,0.1,0.1,0.0413149,0.0397909);
+ g->SetPoint(5,0.1,4.67937); g->SetPointError(5,0.1,0.1,0.0346151,0.0450248);
+ g->SetPoint(6,0.3,4.7109); g->SetPointError(6,0.1,0.1,0.0408403,0.0839992);
+ g->SetPoint(7,0.5,4.79359); g->SetPointError(7,0.1,0.1,0.0324516,0.0357053);
+ g->SetPoint(8,0.7,4.86469); g->SetPointError(8,0.1,0.1,0.0452175,0.0477304);
+ g->SetPoint(9,0.9,4.96078); g->SetPointError(9,0.1,0.1,0.0566798,0.0804077);
+
+ SetGraphAttributes(g, INELGt0, WIP, false,
+ "alice_ppInelGt2760Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
- return g;
-}
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 7000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel7000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,4.97541); g->SetPointError(0,0.1,0.1,0.187526,0.162049);
+ g->SetPoint(1,-1.5,4.98161); g->SetPointError(1,0.1,0.1,0.128353,0.149085);
+ g->SetPoint(2,-1.3,4.94853); g->SetPointError(2,0.1,0.1,0.129841,0.144762);
+ g->SetPoint(3,-1.1,4.88924); g->SetPointError(3,0.1,0.1,0.137866,0.157862);
+ g->SetPoint(4,-0.9,4.79998); g->SetPointError(4,0.1,0.1,0.144492,0.158783);
+ g->SetPoint(5,-0.7,4.71399); g->SetPointError(5,0.1,0.1,0.132703,0.156135);
+ g->SetPoint(6,-0.5,4.63098); g->SetPointError(6,0.1,0.1,0.129938,0.147085);
+ g->SetPoint(7,-0.3,4.56815); g->SetPointError(7,0.1,0.1,0.129424,0.145485);
+ g->SetPoint(8,-0.1,4.52372); g->SetPointError(8,0.1,0.1,0.129049,0.145285);
+ g->SetPoint(9,0.1,4.52946); g->SetPointError(9,0.1,0.1,0.131266,0.144285);
+ g->SetPoint(10,0.3,4.56411); g->SetPointError(10,0.1,0.1,0.130652,0.149019);
+ g->SetPoint(11,0.5,4.63554); g->SetPointError(11,0.1,0.1,0.133415,0.144298);
+ g->SetPoint(12,0.7,4.71592); g->SetPointError(12,0.1,0.1,0.136436,0.151768);
+ g->SetPoint(13,0.9,4.8059); g->SetPointError(13,0.1,0.1,0.136996,0.142551);
+ g->SetPoint(14,1.1,4.88457); g->SetPointError(14,0.1,0.1,0.134237,0.142764);
+ g->SetPoint(15,1.3,4.92903); g->SetPointError(15,0.1,0.1,0.131933,0.152767);
+ g->SetPoint(16,1.5,4.96487); g->SetPointError(16,0.1,0.1,0.140214,0.147354);
+ g->SetPoint(17,1.7,4.95502); g->SetPointError(17,0.1,0.1,0.156906,0.14759);
+
+ SetGraphAttributes(g, INEL, WIP, false,
+ "alice_ppInel7000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 7000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd7000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,6.12747);g->SetPointError(0,0.1,0.1,0.157334,0.164918);
+ g->SetPoint(1,-1.5,6.1353); g->SetPointError(1,0.1,0.1,0.146834,0.152011);
+ g->SetPoint(2,-1.3,6.09648);g->SetPointError(2,0.1,0.1,0.140067,0.143287);
+ g->SetPoint(3,-1.1,6.02552);g->SetPointError(3,0.1,0.1,0.133435,0.135281);
+ g->SetPoint(4,-0.9,5.91705);g->SetPointError(4,0.1,0.1,0.129449,0.130381);
+ g->SetPoint(5,-0.7,5.81246);g->SetPointError(5,0.1,0.1,0.126477,0.126883);
+ g->SetPoint(6,-0.5,5.71104);g->SetPointError(6,0.1,0.1,0.124521,0.124655);
+ g->SetPoint(7,-0.3,5.63422);g->SetPointError(7,0.1,0.1,0.120116,0.12013);
+ g->SetPoint(8,-0.1,5.57977);g->SetPointError(8,0.1,0.1,0.119286,0.119256);
+ g->SetPoint(9,0.1,5.58662); g->SetPointError(9,0.1,0.1,0.119331,0.119301);
+ g->SetPoint(10,0.3,5.6291); g->SetPointError(10,0.1,0.1,0.120683,0.120697);
+ g->SetPoint(11,0.5,5.7166); g->SetPointError(11,0.1,0.1,0.122787,0.122923);
+ g->SetPoint(12,0.7,5.81463);g->SetPointError(12,0.1,0.1,0.126293,0.1267);
+ g->SetPoint(13,0.9,5.92404);g->SetPointError(13,0.1,0.1,0.129522,0.130456);
+ g->SetPoint(14,1.1,6.01958);g->SetPointError(14,0.1,0.1,0.134505,0.136333);
+ g->SetPoint(15,1.3,6.07232);g->SetPointError(15,0.1,0.1,0.140728,0.143909);
+ g->SetPoint(16,1.5,6.11596);g->SetPointError(16,0.1,0.1,0.14756,0.15268);
+ g->SetPoint(17,1.7,6.10155);g->SetPointError(17,0.1,0.1,0.151731,0.159518);
+
+ SetGraphAttributes(g, NSD, WIP, false,
+ "alice_ppNsd7000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 7000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt7000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
+ g->SetPoint(0,-0.9,6.22689);g->SetPointError(0,0.1,0.1,0.0959094,0.10395);
+ g->SetPoint(1,-0.7,6.11603);g->SetPointError(1,0.1,0.1,0.0609242,0.0984269);
+ g->SetPoint(2,-0.5,6.00881);g->SetPointError(2,0.1,0.1,0.0595691,0.0842045);
+ g->SetPoint(3,-0.3,5.9274); g->SetPointError(3,0.1,0.1,0.0560837,0.0780806);
+ g->SetPoint(4,-0.1,5.86988);g->SetPointError(4,0.1,0.1,0.0552611,0.0798584);
+ g->SetPoint(5,0.1,5.8773); g->SetPointError(5,0.1,0.1,0.062512,0.077947);
+ g->SetPoint(6,0.3,5.92215); g->SetPointError(6,0.1,0.1,0.0535152,0.0863595);
+ g->SetPoint(7,0.5,6.01458); g->SetPointError(7,0.1,0.1,0.0578218,0.0745799);
+ g->SetPoint(8,0.7,6.1186); g->SetPointError(8,0.1,0.1,0.0767397,0.0899574);
+ g->SetPoint(9,0.9,6.23468); g->SetPointError(9,0.1,0.1,0.0786932,0.073295);
+
+ SetGraphAttributes(g, INELGt0, WIP, false,
+ "alice_ppInelGt7000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 8000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInel8000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,5.1333); g->SetPointError(0,0.1,0.1,0.0613865,0.0707879);
+ g->SetPoint(1,-1.5,5.14341);g->SetPointError(1,0.1,0.1,0.117829,0.121428);
+ g->SetPoint(2,-1.3,5.13589);g->SetPointError(2,0.1,0.1,0.110807,0.113393);
+ g->SetPoint(3,-1.1,5.06167);g->SetPointError(3,0.1,0.1,0.120091,0.121635);
+ g->SetPoint(4,-0.9,4.97796);g->SetPointError(4,0.1,0.1,0.116141,0.117154);
+ g->SetPoint(5,-0.7,4.88431);g->SetPointError(5,0.1,0.1,0.11944,0.120064);
+ g->SetPoint(6,-0.5,4.81236);g->SetPointError(6,0.1,0.1,0.11049,0.110936);
+ g->SetPoint(7,-0.3,4.72239);g->SetPointError(7,0.1,0.1,0.110969,0.111284);
+ g->SetPoint(8,-0.1,4.66962);g->SetPointError(8,0.1,0.1,0.125108,0.125337);
+ g->SetPoint(9,0.1,4.69441); g->SetPointError(9,0.1,0.1,0.113766,0.114021);
+ g->SetPoint(10,0.3,4.7335); g->SetPointError(10,0.1,0.1,0.104531,0.104866);
+ g->SetPoint(11,0.5,4.79917);g->SetPointError(11,0.1,0.1,0.107076,0.107534);
+ g->SetPoint(12,0.7,4.88713);g->SetPointError(12,0.1,0.1,0.106124,0.106827);
+ g->SetPoint(13,0.9,4.98035);g->SetPointError(13,0.1,0.1,0.120107,0.121087);
+ g->SetPoint(14,1.1,5.05366);g->SetPointError(14,0.1,0.1,0.115795,0.11739);
+ g->SetPoint(15,1.3,5.11276);g->SetPointError(15,0.1,0.1,0.123574,0.125877);
+ g->SetPoint(16,1.5,5.16105);g->SetPointError(16,0.1,0.1,0.0979751,0.102305);
+ g->SetPoint(17,1.7,5.16477);g->SetPointError(17,0.1,0.1,0.116096,0.121392);
+
+ SetGraphAttributes(g, INEL, WIP, false,
+ "alice_ppInel8000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 8000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralNsd8000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
+ g->SetPoint(0,-1.7,6.28722);g->SetPointError(0,0.1,0.1,0.220541,0.2263);
+ g->SetPoint(1,-1.5,6.29211);g->SetPointError(1,0.1,0.1,0.126653,0.132922);
+ g->SetPoint(2,-1.3,6.29313);g->SetPointError(2,0.1,0.1,0.142432,0.145805);
+ g->SetPoint(3,-1.1,6.1944); g->SetPointError(3,0.1,0.1,0.116871,0.119093);
+ g->SetPoint(4,-0.9,6.09529);g->SetPointError(4,0.1,0.1,0.121064,0.122121);
+ g->SetPoint(5,-0.7,5.97811);g->SetPointError(5,0.1,0.1,0.111525,0.112012);
+ g->SetPoint(6,-0.5,5.88992);g->SetPointError(6,0.1,0.1,0.119488,0.119637);
+ g->SetPoint(7,-0.3,5.78296);g->SetPointError(7,0.1,0.1,0.114947,0.114962);
+ g->SetPoint(8,-0.1,5.71633);g->SetPointError(8,0.1,0.1,0.0933,0.09326);
+ g->SetPoint(9,0.1,5.74663); g->SetPointError(9,0.1,0.1,0.109892,0.109857);
+ g->SetPoint(10,0.3,5.79472);g->SetPointError(10,0.1,0.1,0.123704,0.123718);
+ g->SetPoint(11,0.5,5.87545);g->SetPointError(11,0.1,0.1,0.122522,0.122667);
+ g->SetPoint(12,0.7,5.98273);g->SetPointError(12,0.1,0.1,0.128316,0.128739);
+ g->SetPoint(13,0.9,6.09037);g->SetPointError(13,0.1,0.1,0.114321,0.115437);
+ g->SetPoint(14,1.1,6.18105);g->SetPointError(14,0.1,0.1,0.125412,0.127476);
+ g->SetPoint(15,1.3,6.24275);g->SetPointError(15,0.1,0.1,0.118631,0.122597);
+ g->SetPoint(16,1.5,6.28916);g->SetPointError(16,0.1,0.1,0.144205,0.149736);
+ g->SetPoint(17,1.7,6.28878);g->SetPointError(17,0.1,0.1,0.134438,0.143695);
+
+ SetGraphAttributes(g, NSD, WIP, false,
+ "alice_ppNsd8000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+ //____________________________________________________________________
+ /**
+ * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
+ * at @f$ \sqrt{s} = 8000GeV@f$
+ * Work in progress
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* AliceCentralInelGt8000Work()
+ {
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
+ g->SetPoint(0,-0.9,6.38567);g->SetPointError(0,0.1,0.1,0.0436571,0.0436571);
+ g->SetPoint(1,-0.7,6.26363);g->SetPointError(1,0.1,0.1,0.0312036,0.0312036);
+ g->SetPoint(2,-0.5,6.17205);g->SetPointError(2,0.1,0.1,0.0351509,0.0351509);
+ g->SetPoint(3,-0.3,6.05629);g->SetPointError(3,0.1,0.1,0.0302028,0.0302028);
+ g->SetPoint(4,-0.1,5.98823);g->SetPointError(4,0.1,0.1,0.0141541,0.0141541);
+ g->SetPoint(5,0.1,6.02043); g->SetPointError(5,0.1,0.1,0.0256893,0.0256893);
+ g->SetPoint(6,0.3,6.07111); g->SetPointError(6,0.1,0.1,0.0380304,0.0380304);
+ g->SetPoint(7,0.5,6.15492); g->SetPointError(7,0.1,0.1,0.0384435,0.0384435);
+ g->SetPoint(8,0.7,6.26781); g->SetPointError(8,0.1,0.1,0.0450579,0.0450579);
+ g->SetPoint(9,0.9,6.38491); g->SetPointError(9,0.1,0.1,0.0396431,0.0396431);
+
+ SetGraphAttributes(g, INELGt0, WIP, false,
+ "alice_ppInelGt8000Work",
+ "PWG-UD/MULT - work in progress");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
+ * \sqrt{s} = 900GeV@f$
+ *
+ * p7743_d8x1y1
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* CMSNsd900()
+ {
+ // CMS published NSD data - p7743_d8x1y1
+ double x[] ={ -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25};
+ double exm[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double exp[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double y[] = { 3.6, 3.73, 3.62, 3.54, 3.48, 3.48, 3.54, 3.62, 3.73, 3.6 };
+ double eym[] ={ 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14,0.13 };
+ double eyp[] ={ 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14,0.13 };
+ const int np = 10;
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, CMS, false, "cms_nsd900", "CMS");
-//____________________________________________________________________
-/**
- * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
- * \sqrt{s} = 2.36GeV@f$
- *
- * p7743_d8x1y2
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* CMSNsd2360()
-{
- // CMS NSD 2360 - p7743_d8x1y2
- double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
- double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double y[] = { 4.78, 4.81, 4.66, 4.61, 4.47, 4.47, 4.61, 4.66, 4.81, 4.78 };
- double eym[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
- double eyp[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
- const int np = 10;
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, CMS, false, "cms_nsd2360", "CMS NSD");
- return g;
-}
+ return g;
+ }
-//____________________________________________________________________
-/**
- * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
- * \sqrt{s} = 7TeV@f$
- *
- * p7838_d5x1y1
- *
- * @return graph of data
- *
- * @ingroup pwglf_forward_otherdata
- */
-TGraphAsymmErrors* CMSNsd7000()
-{
- // CMS NSD 7000 - Plot: p7838_d5x1y1
- double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
- double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
- double y[] = { 6.18, 6.26, 6.14, 6.01, 5.78, 5.78, 6.01, 6.14, 6.26, 6.18 };
- double eym[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
- double eyp[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
- const int np = 10;
- TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
- SetGraphAttributes(g, NSD, CMS, false, "cms_nsd7000", "CMS NSD");
- return g;
-}
+ //____________________________________________________________________
+ /**
+ * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
+ * \sqrt{s} = 2.36GeV@f$
+ *
+ * p7743_d8x1y2
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* CMSNsd2360()
+ {
+ // CMS NSD 2360 - p7743_d8x1y2
+ double x[] ={ -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25};
+ double exm[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double exp[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double y[] = { 4.78, 4.81, 4.66, 4.61, 4.47, 4.47, 4.61, 4.66, 4.81, 4.78};
+ double eym[] ={ 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17};
+ double eyp[] ={ 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17};
+ const int np = 10;
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, CMS, false, "cms_nsd2360", "CMS");
+ return g;
+ }
-//____________________________________________________________________
-TGraphAsymmErrors*
-GetSingle(UShort_t which,
- UShort_t sys,
- UShort_t energy,
- UShort_t type=0x1,
- UShort_t centLow=0,
- UShort_t centHigh=0)
-{
- TGraphAsymmErrors* ret = 0;
- if (sys == 1) {
- if (TMath::Abs(energy-900) < 10) {
- switch (type) {
- case 1: // INEL
- switch (which) {
- case PYTHIA: ret = Pythia900INEL(); break;
- case UA5: ret = UA5Inel(false); break;
- case UA5+10: ret = UA5Inel(true); break;
- case ALICE: ret = AliceCentralInel900(); break;
- }
- break;
- case 2: // INEL>0
- switch (which) {
- case ALICE: ret = AliceCentralInelGt900(); break;
+
+ //____________________________________________________________________
+ /**
+ * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
+ * \sqrt{s} = 7TeV@f$
+ *
+ * p7838_d5x1y1
+ *
+ * @return graph of data
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TGraphAsymmErrors* CMSNsd7000()
+ {
+ // CMS NSD 7000 - Plot: p7838_d5x1y1
+ double x[] ={ -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25};
+ double exm[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double exp[] ={ 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25};
+ double y[] ={ 6.18, 6.26, 6.14, 6.01, 5.78, 5.78, 6.01, 6.14, 6.26, 6.18 };
+ double eym[] ={ 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25};
+ double eyp[] ={ 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25};
+ const int np = 10;
+ TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
+ SetGraphAttributes(g, NSD, CMS, false, "cms_nsd7000", "CMS");
+ return g;
+ }
+
+ //____________________________________________________________________
+ /**
+ * Get a single data graph
+ *
+ * @param which Which type
+ * @param sys Collisition system
+ * @param energy Collision energy
+ * @param type Trigger type
+ * @param centLow Low cut on centrality
+ * @param centHigh Up cut on centraltiy
+ *
+ * @return Data graph
+ */
+ static TGraphAsymmErrors* GetSingle(UShort_t which,
+ UShort_t sys,
+ UShort_t energy,
+ UShort_t type=0x1,
+ UShort_t centLow=0,
+ UShort_t centHigh=0)
+ {
+ TGraphAsymmErrors* ret = 0;
+ if (sys == 1) {
+ if (TMath::Abs(energy-900) < 10) {
+ switch (type) {
+ case 1: // INEL
+ switch (which) {
+ case PYTHIA: ret = Pythia900INEL(); break;
+ case UA5: ret = UA5Inel(false); break;
+ case UA5+10: ret = UA5Inel(true); break;
+ case ALICE: ret = AliceCentralInel900(); break;
+ case WIP: ret = AliceCentralInel900Work(); break;
+ }
+ break;
+ case 2: // INEL>0
+ switch (which) {
+ case ALICE: ret = AliceCentralInelGt900(); break;
+ case WIP: ret = AliceCentralInelGt900Work(); break;
+ }
+ break;
+ case 4: // NSD
+ switch (which) {
+ case PYTHIA: ret = Pythia900NSD(); break;
+ case UA5: ret = UA5Nsd(false); break;
+ case UA5+10: ret = UA5Nsd(true); break;
+ case ALICE: ret = AliceCentralNsd900(); break;
+ case WIP: ret = AliceCentralNsd900Work(); break;
+ case CMS: ret = CMSNsd900(); break;
+ }
+ break;
+ } // type
+ }
+ else if (TMath::Abs(energy-2360) < 10) {
+ switch (type) {
+ case 1: // INEL
+ switch (which) {
+ case ALICE: ret = AliceCentralInel2360(); break;
+ case WIP: ret = AliceCentralInel2760Work(); break;
+ }
+ break;
+ case 2: // INEL > 0
+ switch (which) {
+ case ALICE: ret = AliceCentralInelGt2360(); break;
+ case WIP: ret = AliceCentralInelGt2760Work(); break;
+ }
+ break;
+ case 4: // NSD
+ switch (which) {
+ case ALICE: ret = AliceCentralNsd2360(); break;
+ case CMS: ret = CMSNsd2360(); break;
+ case WIP: ret = AliceCentralNsd2760Work(); break;
+ }
+ break;
}
- break;
- case 4: // NSD
- switch (which) {
- case PYTHIA: ret = Pythia900NSD(); break;
- case UA5: ret = UA5Nsd(false); break;
- case UA5+10: ret = UA5Nsd(true); break;
- case ALICE: ret = AliceCentralNsd900(); break;
- case CMS: ret = CMSNsd900(); break;
+ }
+ else if (TMath::Abs(energy-2760) < 10) {
+ switch (type) {
+ case 1: // INEL
+ switch (which) {
+ case WIP: ret = AliceCentralInel2760Work(); break;
+ }
+ break;
+ case 2: // INEL > 0
+ switch (which) {
+ case WIP: ret = AliceCentralInelGt2760Work(); break;
+ }
+ break;
+ case 4: // NSD
+ switch (which) {
+ case WIP: ret = AliceCentralNsd2760Work(); break;
+ }
+ break;
}
- break;
- } // type
- }
- else if (TMath::Abs(energy-2360) < 10) {
- switch (type) {
- case 1: // INEL
- switch (which) {
- case ALICE: ret = AliceCentralInel2360(); break;
+ }
+ else if (TMath::Abs(energy-7000) < 10) {
+ switch (type) {
+ case 1:
+ switch (which) {
+ case WIP: ret = AliceCentralInel7000Work(); break;
+ }
+ break;
+ case 2: // INEL > 0
+ switch (which) {
+ case ALICE: ret = AliceCentralInelGt7000(); break;
+ case WIP: ret = AliceCentralInelGt7000Work(); break;
+ }
+ break;
+ case 4: // NSD
+ switch (which) {
+ case CMS: ret = CMSNsd7000(); break;
+ case WIP: ret = AliceCentralNsd7000Work(); break;
+ }
+ break;
}
- break;
- case 2: // INEL > 0
- switch (which) {
- case ALICE: ret = AliceCentralInelGt2360(); break;
+ }
+ else if (TMath::Abs(energy-8000) < 10) {
+ switch (type) {
+ case 1:
+ switch (which) {
+ case WIP: ret = AliceCentralInel8000Work(); break;
+ }
+ break;
+ case 2: // INEL > 0
+ switch (which) {
+ case WIP: ret = AliceCentralInelGt8000Work(); break;
+ }
+ break;
+ case 4: // NSD
+ switch (which) {
+ case WIP: ret = AliceCentralNsd8000Work(); break;
+ }
+ break;
}
- break;
- case 4: // NSD
+ }
+ }
+ else if (sys == 3) { // pPb
+ if (TMath::Abs(energy - 5023) < 10 ||
+ TMath::Abs(energy - 8000) < 10) {
switch (which) {
- case ALICE: ret = AliceCentralNsd2360(); break;
- case CMS: ret = CMSNsd2360(); break;
+ case ALICE: ret = AliceCentralpPb5023(); break;
}
- break;
}
}
- else if (TMath::Abs(energy-7000) < 10) {
+ if (!ret) {
+ TString w;
+ switch (which) {
+ case UA5: w = "UA5"; break;
+ case CMS: w = "CMS"; break;
+ case ALICE: w = "ALICE"; break;
+ case WIP: w = "WIP"; break;
+ case PYTHIA:w = "Pyhthia"; break;
+ default: w = Form("unknown(%d)", which);
+ }
+ TString sy;
+ switch (sys) {
+ case 1: sy = "pp"; break;
+ case 2: sy = "PbPb"; break;
+ case 3: sy = "pPb"; break;
+ default: sy = Form("unknown(%d)", sys);
+ }
+ TString tr;
switch (type) {
- case 1: ret = 0; break;
- case 2: // INEL > 0
- switch (which) {
- case ALICE: ret = AliceCentralInelGt7000(); break;
- }
- break;
- case 4: // NSD
- switch (which) {
- case CMS: ret = CMSNsd7000(); break;
- }
- break;
+ case 1: tr = "INEL"; break;
+ case 2: tr = "INEL>0"; break;
+ case 4: tr = "NSD"; break;
+ default: tr = Form("unknown(%d)", sys);
}
+
+ Warning("GetSingle", "Nothing to get for "
+ "which=%s, sys=%s, energy=%dGeV, type=%s, "
+ "centLow=%d, centHigh=%d",
+ w.Data(), sy.Data(), energy, tr.Data(), centLow, centHigh);
}
- }
#if 0
- if (ret) {
- if (!wn.IsNull()) wn.Append(",");
- switch (which) {
- case PYTHIA: wn.Append("Pythia"); break;
- case UA5: wn.Append("UA5"); break;
- case CMS: wn.Append("CMS"); break;
- case ALICE: wn.Append("ALICE"); break;
- }
- }
+ if (ret) {
+ if (!wn.IsNull()) wn.Append(",");
+ switch (which) {
+ case PYTHIA: wn.Append("Pythia"); break;
+ case UA5: wn.Append("UA5"); break;
+ case CMS: wn.Append("CMS"); break;
+ case ALICE: wn.Append("ALICE"); break;
+ }
+ }
#endif
- return ret;
-}
-
-//____________________________________________________________________
-TString&
-AppendItem(TString& s, char delim, const char* what)
-{
- if (!s.IsNull()) s.Append(Form("%c", delim));
- s.Append(what);
- return s;
-}
-
-//____________________________________________________________________
-/**
- * Get a multi graph of data for a given energy and trigger type
- *
- * @param sys Collision system (1: pp, 2: PbPb)
- * @param energy Energy in GeV (900, 2360, 7000)
- * @param type Bit pattern of trigger type
- * - 0x1 INEL
- * - 0x2 INEL>0
- * - 0x4 NSD
- * @param centLow Low centrality cut (only for PbPB)
- * @param centHigh High centrality cut (only for PbPB)
- * @param which What to get
- *
- * @return A multi graph with the selected data.
- *
- * @ingroup pwglf_forward_otherdata
- */
-TMultiGraph*
-GetData(UShort_t sys,
- UShort_t energy,
- UShort_t type=0x1,
- UShort_t centLow=0,
- UShort_t centHigh=0,
- UShort_t which=0x7)
-{
- TMultiGraph* mp = new TMultiGraph(Form("dndeta_%dGeV_%d_%03d_%03d",
- energy, type, centLow, centHigh),"");
- TString tn;
- TString en;
- TString sn;
- TString cn;
- bool ua5 = (which & (1 << UA5)); // 0x1
- bool cms = (which & (1 << CMS)); // 0x2
- bool alice = (which & (1 << ALICE)); // 0x4
- bool pythia = (which & (1 << PYTHIA)); // 0x8
-
- en.Append(Form(", #sqrt{s%s}=", sys == 1 ? "" : "_{NN}"));
- if (energy < 1000)
- en.Append(Form("%dGeV", energy));
- else {
- if (energy % 1000 == 0)
- en.Append(Form("%dTeV", energy/1000));
- else
- en.Append(Form("%4.2fTeV", float(energy)/1000));
+ return ret;
}
- if (sys == 1) {
- if (!(type & 0x7))
- Warning("GetData", "Unknown trigger mask 0x%x", type);
+ //____________________________________________________________________
+ /**
+ * Append an item to a list
+ *
+ * @param s List to append to
+ * @param delim Delimiter
+ * @param what What to append
+ *
+ * @return New string value
+ */
+ static TString& AppendItem(TString& s, char delim, const char* what)
+ {
+ if (!s.IsNull()) s.Append(Form("%c", delim));
+ s.Append(what);
+ return s;
+ }
- if (TMath::Abs(energy-2750) < 11) {
- Warning("GetData", "Using 2360GeV data for %dGeV comparison", energy);
- energy = 2360;
- }
- if (!(TMath::Abs(energy-900) < 10 ||
- TMath::Abs(energy-2360) < 10 ||
- TMath::Abs(energy-7000) < 10)) {
- Warning("GetData", "No other results for sys=%d, energy=%d",
- sys, energy);
- return 0;
+ //____________________________________________________________________
+ /**
+ * Get a multi graph of data for a given energy and trigger type
+ *
+ * @param sys Collision system (1: pp, 2: PbPb)
+ * @param energy Energy in GeV (900, 2360, 7000)
+ * @param type Bit pattern of trigger type
+ * - 0x1 INEL
+ * - 0x2 INEL>0
+ * - 0x4 NSD
+ * @param centLow Low centrality cut (only for PbPB)
+ * @param centHigh High centrality cut (only for PbPB)
+ * @param which What to get
+ *
+ * @return A multi graph with the selected data.
+ *
+ * @ingroup pwglf_forward_otherdata
+ */
+ static TMultiGraph* GetData(UShort_t sys,
+ UShort_t energy,
+ UShort_t type=0x1,
+ UShort_t centLow=0,
+ UShort_t centHigh=0,
+ UShort_t which=0x7)
+ {
+ TMultiGraph* mp = new TMultiGraph(Form("dndeta_%dGeV_%d_%03d_%03d",
+ energy, type, centLow, centHigh),"");
+ TString tn;
+ TString en;
+ TString sn;
+ TString cn;
+ bool ua5 = (which & (1 << UA5)); // 0x1
+ bool cms = (which & (1 << CMS)); // 0x2
+ bool alice = (which & (1 << ALICE)); // 0x4
+ bool work = (which & (1 << WIP)); // 0x8
+ bool pythia = (which & (1 << PYTHIA)); // 0x10
+ en.Append(Form(", #sqrt{s%s}=", sys == 1 ? "" : "_{NN}"));
+ if (energy < 1000)
+ en.Append(Form("%dGeV", energy));
+ else {
+ if (energy % 1000 == 0)
+ en.Append(Form("%dTeV", energy/1000));
+ else
+ en.Append(Form("%4.2fTeV", float(energy)/1000));
}
+
+ if (sys == 1) {
+ if (!(type & 0x7))
+ Warning("GetData", "Unknown trigger mask 0x%x", type);
+
+ if (TMath::Abs(energy-2750) < 11) {
+ Warning("GetData", "Using 2360GeV data for %dGeV comparison", energy);
+ energy = 2360;
+ }
+ if (!(TMath::Abs(energy-900) < 10 ||
+ TMath::Abs(energy-2360) < 10 ||
+ TMath::Abs(energy-7000) < 10 ||
+ TMath::Abs(energy-8000) < 10)) {
+ Warning("GetData", "No other results for sys=%d, energy=%d",
+ sys, energy);
+ return 0;
+ }
- sn = "pp";
+ sn = "pp";
- if (type & 0x1) AppendItem(tn, '|', "INEL");
- if (type & 0x2) AppendItem(tn, '|', "INEL>0");
- if (type & 0x4) AppendItem(tn, '|', "NSD");
+ // Substitute NSD for V0-AND
+ if (type == 0x2000) type = 0x4;
+ if (type & 0x1) AppendItem(tn, '|', "INEL");
+ if (type & 0x2) AppendItem(tn, '|', "INEL>0");
+ if (type & 0x4) AppendItem(tn, '|', "NSD");
- Bool_t seenUA5 = false;
- for (Int_t i = 0; i < 3; i++) {
- UShort_t mask = (1 << i);
- if ((type & mask) == 0) continue;
- TGraphAsymmErrors* gUAp =(ua5 ?GetSingle(UA5, sys,energy,mask):0);
- TGraphAsymmErrors* gUAn =(ua5 ?GetSingle(UA5+10,sys,energy,mask):0);
- TGraphAsymmErrors* gCMS =(cms ?GetSingle(CMS, sys,energy,mask):0);
- TGraphAsymmErrors* gALI =(alice ?GetSingle(ALICE, sys,energy,mask):0);
- TGraphAsymmErrors* gPYT =(pythia?GetSingle(PYTHIA,sys,energy,mask):0);
- if (gUAp) mp->Add(gUAp);
- if (gUAn) mp->Add(gUAn);
- if (gCMS) mp->Add(gCMS);
+ Bool_t seenUA5 = false;
+ for (Int_t i = 0; i < 3; i++) {
+ UShort_t mask = (1 << i);
+ if ((type & mask) == 0) continue;
+ TGraphAsymmErrors* gUAp =(ua5 ?GetSingle(UA5, sys,energy,mask):0);
+ TGraphAsymmErrors* gUAn =(ua5 ?GetSingle(UA5+10,sys,energy,mask):0);
+ TGraphAsymmErrors* gCMS =(cms ?GetSingle(CMS, sys,energy,mask):0);
+ TGraphAsymmErrors* gALI =(alice ?GetSingle(ALICE, sys,energy,mask):0);
+ TGraphAsymmErrors* gPYT =(pythia?GetSingle(PYTHIA,sys,energy,mask):0);
+ TGraphAsymmErrors* gWRK =(work ?GetSingle(WIP, sys,energy,mask):0);
+ if (gUAp) mp->Add(gUAp);
+ if (gUAn) mp->Add(gUAn);
+ if (gCMS) mp->Add(gCMS);
+ if (gALI) mp->Add(gALI);
+ if (gPYT) mp->Add(gPYT);
+ if (gWRK) mp->Add(gWRK);
+ if (gUAp || gUAn) seenUA5 = true;
+ }
+ if (seenUA5) sn.Append("(p#bar{p})");
+ }
+ else if (sys == 2) {
+ // Nothing for PbPb so far
+ cn = Form(", %d%%-%d%% central", centLow, centHigh);
+ sn = "PbPb";
+ // Warning("GetData", "No other data for PbPb yet");
+ }
+ else if (sys == 3) {
+ if (!(TMath::Abs(energy-5023) < 10 ||
+ TMath::Abs(energy-8000) < 10) ) {
+ Warning("GetData", "No other results for sys=%d, energy=%d",
+ sys, energy);
+ return 0;
+ }
+
+ // Info("GetData", "Getting ALICE pPb data");
+ sn = "pPb";
+ TGraphAsymmErrors* gALI =(alice ?GetSingle(ALICE, sys,energy, 0):0);
if (gALI) mp->Add(gALI);
- if (gPYT) mp->Add(gPYT);
- if (gUAp || gUAn) seenUA5 = true;
+ // Info("GetData", "Got %p", gALI);
+ // Warning("GetData", "Unknown system %d", sys);
}
- if (seenUA5) sn.Append("(p#bar{p})");
- }
- else if (sys == 2) {
- // Nothing for PbPb so far
- cn = Form(", %d%%-%d%% central", centLow, centHigh);
- sn = "PbPb";
- // Warning("GetData", "No other data for PbPb yet");
- }
- else
- Warning("GetData", "Unknown system %d", sys);
- if (!mp->GetListOfGraphs() || mp->GetListOfGraphs()->GetEntries() <= 0) {
- delete mp;
- mp = 0;
- return 0;
+ if (!mp->GetListOfGraphs() || mp->GetListOfGraphs()->GetEntries() <= 0) {
+ delete mp;
+ mp = 0;
+ return 0;
+ }
+ TString tit(Form("%s%s, %s%s",
+ sn.Data(), en.Data(), tn.Data(), cn.Data()));
+ mp->SetTitle(tit.Data());
+ return mp;
}
- TString tit(Form("%s%s, %s%s",
- sn.Data(), en.Data(), tn.Data(), cn.Data()));
- mp->SetTitle(tit.Data());
- return mp;
-}
-
+};
//____________________________________________________________________
/**
* Plot external data for a given selection of energy and trigger type
UShort_t type=0x1,
UShort_t centLow=0,
UShort_t centHigh=5,
- UShort_t which=0x7)
+ UShort_t which=0xf)
{
- TMultiGraph* mp = GetData(sys, energy, type, centLow, centHigh, which);
+ TMultiGraph* mp =
+ RefData::GetData(sys, energy, type, centLow, centHigh, which);
if (!mp) return;
gStyle->SetTitleX(0.1);
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff