]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/FORWARD/analysis2/OtherData.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
First version of the new GUI in development. You can enable this UI with "alieve...
[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / OtherData.C
CommitLineData
ffca499d 1//____________________________________________________________________
2/**
bd6f5206 3 * @defgroup pwglf_forward_otherdata External data
ffca499d 4 *
290052e7 5 * Collection of external data points for comparisons and the like
6 *
bd6f5206 7 * @ingroup pwglf_forward_scripts
ffca499d 8 */
7c1a1f1d 9/**
10 * @file
11 *
bd6f5206 12 * @ingroup pwglf_forward_script_otherdata
7c1a1f1d 13 */
05cc5a69 14#ifndef __CINT__
15# include <TGraphAsymmErrors.h>
16# include <TMultiGraph.h>
17# include <TStyle.h>
18# include <TMath.h>
19# include <TCanvas.h>
20# include <TLegend.h>
21# include <TLegendEntry.h>
22# include <TList.h>
23# include <TAxis.h>
24# include <TH1F.h>
25#else
26class TGraphAsymmErrors;
27class TMultiGraph;
28class TStyle;
29class TCanvas;
30class TLegend;
31class TLegendEntry;
32class TList;
33class TAxis;
34class TH1F;
35class TGraph;
36#endif
7e4038b5 37
7e4038b5 38//____________________________________________________________________
39/**
40 * Values used
41 *
bd6f5206 42 * @ingroup pwglf_forward_otherdata
7e4038b5 43 */
44enum {
ffca499d 45 UA5,
46 CMS,
47 ALICE,
bfab35d9 48 WIP,
b6b35c77 49 PYTHIA,
ffca499d 50 INEL,
51 INELGt0,
52 NSD
53};
54enum {
55 /** Style used for UA5 data */
56 UA5Style = 21,
57 /** Style used for CMS data */
58 CMSStyle = 29,
59 /** Style used for ALICE published data */
60 ALICEStyle = 27,
bfab35d9 61 /** Color used for ALICE work-in-progress data */
62 WIPStyle = 33,
b6b35c77 63 /** Style used for Pythia data */
64 PYTHIAStyle = 28,
7e4038b5 65 /** Color used for UA5 data */
66 UA5Color = kBlue+1,
b6b35c77 67 /** Color used for Pytia data */
68 PYTHIAColor = kGray+2,
7e4038b5 69 /** Color used for CMS data */
70 CMSColor = kGreen+1,
71 /** Color used for ALICE data */
bfab35d9 72 ALICEColor = kMagenta+1,
73 /** Color used for ALICE work-in-progress data */
74 WIPColor = kMagenta+3
ffca499d 75};
76enum {
7e4038b5 77 /** Marker style INEL data */
ffca499d 78 INELStyle = 22,
7e4038b5 79 /** Marker style INEL>0 data */
ffca499d 80 INELGt0Style= 29,
7e4038b5 81 /** Marker style NSD data */
ffca499d 82 NSDStyle = 23,
83 /** Color used for UA5 data */
84 INELColor = kBlue+1,
85 /** Color used for CMS data */
86 INELGt0Color = kGreen+1,
87 /** Color used for ALICE data */
88 NSDColor = kMagenta+1
89};
90enum {
91 /** Style offset for mirror data */
7e4038b5 92 MirrorOff = 4
93};
94
95//____________________________________________________________________
96/**
ffca499d 97 * Set graph attributes based on trigger type and experiment.
7e4038b5 98 *
99 * @param g Graph
ffca499d 100 * @param exp Experiment
101 * @param mirror True if mirrored data
7e4038b5 102 * @param name Name of graph
103 * @param title Title of graph
104 *
bd6f5206 105 * @ingroup pwglf_forward_otherdata
7e4038b5 106 */
107void
bfab35d9 108SetGraphAttributes(TGraph* g, Int_t trig, Int_t exp, bool mirror,
7e4038b5 109 const Char_t* name, const Char_t* title)
110{
ffca499d 111 Int_t color = 0;
112 switch (exp) {
bfab35d9 113 case UA5: color = UA5Color; break;
114 case CMS: color = CMSColor; break;
115 case ALICE: color = ALICEColor; break;
116 case WIP: color = WIPColor; break;
117 case PYTHIA: color = PYTHIAColor; break;
ffca499d 118 }
119 Int_t style = 0;
120 switch (exp) {
bfab35d9 121 case UA5: style = UA5Style; break;
122 case CMS: style = CMSStyle; break;
123 case ALICE: style = ALICEStyle; break;
124 case WIP: style = WIPStyle; break;
125 case PYTHIA: style = PYTHIAStyle; break;
ffca499d 126 }
127 Float_t size = g->GetMarkerSize();
128 switch (style) {
129 case 21:
130 case 25: size *= 0.8; break;
131 case 27: size *= 1.4; break;
bfab35d9 132 case 33: size *= 1.4; break;
ffca499d 133 }
134
135 if (mirror) style += MirrorOff;
136
7e4038b5 137 g->SetName(name);
138 g->SetTitle(title);
ffca499d 139 g->SetMarkerStyle(style);
140 g->SetMarkerSize(size);
7e4038b5 141 g->SetMarkerColor(color);
142 g->SetLineColor(color);
143 g->SetFillColor(0);
144 g->SetFillStyle(0);
145 g->GetHistogram()->SetStats(kFALSE);
146 g->GetHistogram()->SetXTitle("#eta");
147 g->GetHistogram()->SetYTitle("#frac{1}{N} #frac{dN_{ch}}{#eta}");
148}
149
b6b35c77 150//____________________________________________________________________
290052e7 151/**
152 * Get PYTHIA 900GeV INEL data
153 *
154 *
155 * @return Data graph
156 */
b6b35c77 157TGraphAsymmErrors*
158Pythia900INEL()
159{
160
161 TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
162 SetGraphAttributes(gre, INEL, PYTHIA, false, "pythia900Inel",
163 "Pythia INEL");
164 gre->SetPoint(0,-3.95,1.78199);
165 gre->SetPointError(0, 0, 0, 0.0145305, 0.0145305);
166 gre->SetPoint(1,-3.85,1.85486);
167 gre->SetPointError(1,0,0,0.0148246,0.0148246);
168 gre->SetPoint(2,-3.75,1.93886);
169 gre->SetPointError(2,0,0,0.0151566,0.0151566);
170 gre->SetPoint(3,-3.65,1.96055);
171 gre->SetPointError(3,0,0,0.0152411,0.0152411);
172 gre->SetPoint(4,-3.55,1.98756);
173 gre->SetPointError(4,0,0,0.0153458,0.0153458);
174 gre->SetPoint(5,-3.45,2.02844);
175 gre->SetPointError(5,0,0,0.0155028,0.0155028);
176 gre->SetPoint(6,-3.35,2.09585);
177 gre->SetPointError(6,0,0,0.0157583,0.0157583);
178 gre->SetPoint(7,-3.25,2.13732);
179 gre->SetPointError(7,0,0,0.0159134,0.0159134);
180 gre->SetPoint(8,-3.15,2.1686);
181 gre->SetPointError(8,0,0,0.0160295,0.0160295);
182 gre->SetPoint(9,-3.05,2.25296);
183 gre->SetPointError(9,0,0,0.0163383,0.0163383);
184 gre->SetPoint(10,-2.95,2.29265);
185 gre->SetPointError(10,0,0,0.0164815,0.0164815);
186 gre->SetPoint(11,-2.85,2.34799);
187 gre->SetPointError(11,0,0,0.0166792,0.0166792);
188 gre->SetPoint(12,-2.75,2.35652);
189 gre->SetPointError(12,0,0,0.0167095,0.0167095);
190 gre->SetPoint(13,-2.65,2.40545);
191 gre->SetPointError(13,0,0,0.0168821,0.0168821);
192 gre->SetPoint(14,-2.55,2.43934);
193 gre->SetPointError(14,0,0,0.0170006,0.0170006);
194 gre->SetPoint(15,-2.45,2.45735);
195 gre->SetPointError(15,0,0,0.0170633,0.0170633);
196 gre->SetPoint(16,-2.35,2.48945);
197 gre->SetPointError(16,0,0,0.0171744,0.0171744);
198 gre->SetPoint(17,-2.25,2.51635);
199 gre->SetPointError(17,0,0,0.0172669,0.0172669);
200 gre->SetPoint(18,-2.15,2.55047);
201 gre->SetPointError(18,0,0,0.0173836,0.0173836);
202 gre->SetPoint(19,-2.05,2.58021);
203 gre->SetPointError(19,0,0,0.0174846,0.0174846);
204 gre->SetPoint(20,-1.95,2.58732);
205 gre->SetPointError(20,0,0,0.0175087,0.0175087);
206 gre->SetPoint(21,-1.85,2.60095);
207 gre->SetPointError(21,0,0,0.0175547,0.0175547);
208 gre->SetPoint(22,-1.75,2.59941);
209 gre->SetPointError(22,0,0,0.0175495,0.0175495);
210 gre->SetPoint(23,-1.65,2.63021);
211 gre->SetPointError(23,0,0,0.0176532,0.0176532);
212 gre->SetPoint(24,-1.55,2.61043);
213 gre->SetPointError(24,0,0,0.0175867,0.0175867);
214 gre->SetPoint(25,-1.45,2.61363);
215 gre->SetPointError(25,0,0,0.0175975,0.0175975);
216 gre->SetPoint(26,-1.35,2.60829);
217 gre->SetPointError(26,0,0,0.0175795,0.0175795);
218 gre->SetPoint(27,-1.25,2.61434);
219 gre->SetPointError(27,0,0,0.0175999,0.0175999);
220 gre->SetPoint(28,-1.15,2.61327);
221 gre->SetPointError(28,0,0,0.0175963,0.0175963);
222 gre->SetPoint(29,-1.05,2.57145);
223 gre->SetPointError(29,0,0,0.0174549,0.0174549);
224 gre->SetPoint(30,-0.95,2.55723);
225 gre->SetPointError(30,0,0,0.0174066,0.0174066);
226 gre->SetPoint(31,-0.85,2.57879);
227 gre->SetPointError(31,0,0,0.0174798,0.0174798);
228 gre->SetPoint(32,-0.75,2.516);
229 gre->SetPointError(32,0,0,0.0172657,0.0172657);
230 gre->SetPoint(33,-0.65,2.53709);
231 gre->SetPointError(33,0,0,0.0173379,0.0173379);
232 gre->SetPoint(34,-0.55,2.51197);
233 gre->SetPointError(34,0,0,0.0172519,0.0172519);
234 gre->SetPoint(35,-0.45,2.44052);
235 gre->SetPointError(35,0,0,0.0170047,0.0170047);
236 gre->SetPoint(36,-0.35,2.44882);
237 gre->SetPointError(36,0,0,0.0170336,0.0170336);
238 gre->SetPoint(37,-0.25,2.45308);
239 gre->SetPointError(37,0,0,0.0170484,0.0170484);
240 gre->SetPoint(38,-0.15,2.4622);
241 gre->SetPointError(38,0,0,0.0170801,0.0170801);
242 gre->SetPoint(39,-0.05,2.45735);
243 gre->SetPointError(39,0,0,0.0170633,0.0170633);
244 gre->SetPoint(40,0.05,2.49254);
245 gre->SetPointError(40,0,0,0.017185,0.017185);
246 gre->SetPoint(41,0.15,2.49479);
247 gre->SetPointError(41,0,0,0.0171928,0.0171928);
248 gre->SetPoint(42,0.25,2.49289);
249 gre->SetPointError(42,0,0,0.0171862,0.0171862);
250 gre->SetPoint(43,0.35,2.4628);
251 gre->SetPointError(43,0,0,0.0170822,0.0170822);
252 gre->SetPoint(44,0.45,2.51422);
253 gre->SetPointError(44,0,0,0.0172596,0.0172596);
254 gre->SetPoint(45,0.55,2.51268);
255 gre->SetPointError(45,0,0,0.0172543,0.0172543);
256 gre->SetPoint(46,0.65,2.51066);
257 gre->SetPointError(46,0,0,0.0172474,0.0172474);
258 gre->SetPoint(47,0.75,2.53661);
259 gre->SetPointError(47,0,0,0.0173363,0.0173363);
260 gre->SetPoint(48,0.85,2.54479);
261 gre->SetPointError(48,0,0,0.0173642,0.0173642);
262 gre->SetPoint(49,0.95,2.55391);
263 gre->SetPointError(49,0,0,0.0173953,0.0173953);
264 gre->SetPoint(50,1.05,2.5872);
265 gre->SetPointError(50,0,0,0.0175083,0.0175083);
266 gre->SetPoint(51,1.15,2.60344);
267 gre->SetPointError(51,0,0,0.0175631,0.0175631);
268 gre->SetPoint(52,1.25,2.60616);
269 gre->SetPointError(52,0,0,0.0175723,0.0175723);
270 gre->SetPoint(53,1.35,2.62156);
271 gre->SetPointError(53,0,0,0.0176242,0.0176242);
272 gre->SetPoint(54,1.45,2.61173);
273 gre->SetPointError(54,0,0,0.0175911,0.0175911);
274 gre->SetPoint(55,1.55,2.60415);
275 gre->SetPointError(55,0,0,0.0175655,0.0175655);
276 gre->SetPoint(56,1.65,2.60723);
277 gre->SetPointError(56,0,0,0.0175759,0.0175759);
278 gre->SetPoint(57,1.75,2.60427);
279 gre->SetPointError(57,0,0,0.0175659,0.0175659);
280 gre->SetPoint(58,1.85,2.56765);
281 gre->SetPointError(58,0,0,0.017442,0.017442);
282 gre->SetPoint(59,1.95,2.58602);
283 gre->SetPointError(59,0,0,0.0175043,0.0175043);
284 gre->SetPoint(60,2.05,2.55936);
285 gre->SetPointError(60,0,0,0.0174138,0.0174138);
286 gre->SetPoint(61,2.15,2.54858);
287 gre->SetPointError(61,0,0,0.0173771,0.0173771);
288 gre->SetPoint(62,2.25,2.5205);
289 gre->SetPointError(62,0,0,0.0172811,0.0172811);
290 gre->SetPoint(63,2.35,2.49491);
291 gre->SetPointError(63,0,0,0.0171932,0.0171932);
292 gre->SetPoint(64,2.45,2.42773);
293 gre->SetPointError(64,0,0,0.0169601,0.0169601);
294 gre->SetPoint(65,2.55,2.42879);
295 gre->SetPointError(65,0,0,0.0169638,0.0169638);
296 gre->SetPoint(66,2.65,2.39372);
297 gre->SetPointError(66,0,0,0.0168409,0.0168409);
298 gre->SetPoint(67,2.75,2.38412);
299 gre->SetPointError(67,0,0,0.0168071,0.0168071);
300 gre->SetPoint(68,2.85,2.31896);
301 gre->SetPointError(68,0,0,0.0165758,0.0165758);
302 gre->SetPoint(69,2.95,2.26209);
303 gre->SetPointError(69,0,0,0.0163713,0.0163713);
304 gre->SetPoint(70,3.05,2.24313);
305 gre->SetPointError(70,0,0,0.0163026,0.0163026);
306 gre->SetPoint(71,3.15,2.20403);
307 gre->SetPointError(71,0,0,0.0161599,0.0161599);
308 gre->SetPoint(72,3.25,2.12855);
309 gre->SetPointError(72,0,0,0.0158808,0.0158808);
310 gre->SetPoint(73,3.35,2.13104);
311 gre->SetPointError(73,0,0,0.01589,0.01589);
312 gre->SetPoint(74,3.45,2.06339);
313 gre->SetPointError(74,0,0,0.0156358,0.0156358);
314 gre->SetPoint(75,3.55,1.9846);
315 gre->SetPointError(75,0,0,0.0153343,0.0153343);
316 gre->SetPoint(76,3.65,1.95391);
317 gre->SetPointError(76,0,0,0.0152153,0.0152153);
318 gre->SetPoint(77,3.75,1.87998);
319 gre->SetPointError(77,0,0,0.0149247,0.0149247);
320 gre->SetPoint(78,3.85,1.86256);
321 gre->SetPointError(78,0,0,0.0148554,0.0148554);
322 gre->SetPoint(79,3.95,1.77239);
323 gre->SetPointError(79,0,0,0.0144913,0.0144913);
324 gre->SetPoint(80,4.05,1.72855);
325 gre->SetPointError(80,0,0,0.014311,0.014311);
326 gre->SetPoint(81,4.15,1.69479);
327 gre->SetPointError(81,0,0,0.0141705,0.0141705);
328 gre->SetPoint(82,4.25,1.64147);
329 gre->SetPointError(82,0,0,0.0139459,0.0139459);
330 gre->SetPoint(83,4.35,1.58116);
331 gre->SetPointError(83,0,0,0.0136873,0.0136873);
332 gre->SetPoint(84,4.45,1.55735);
333 gre->SetPointError(84,0,0,0.0135838,0.0135838);
334 gre->SetPoint(85,4.55,1.48815);
335 gre->SetPointError(85,0,0,0.0132786,0.0132786);
336 gre->SetPoint(86,4.65,1.40853);
337 gre->SetPointError(86,0,0,0.0129185,0.0129185);
338 gre->SetPoint(87,4.75,1.36979);
339 gre->SetPointError(87,0,0,0.0127396,0.0127396);
340 gre->SetPoint(88,4.85,1.32666);
341 gre->SetPointError(88,0,0,0.0125374,0.0125374);
342 gre->SetPoint(89,4.95,1.29763);
343 gre->SetPointError(89,0,0,0.0123995,0.0123995);
344 gre->SetPoint(90,5.05,1.25533);
345 gre->SetPointError(90,0,0,0.0121957,0.0121957);
346 gre->SetPoint(91,5.15,1.20912);
347 gre->SetPointError(91,0,0,0.0119692,0.0119692);
348 gre->SetPoint(92,5.25,1.18839);
349 gre->SetPointError(92,0,0,0.0118661,0.0118661);
350 gre->SetPoint(93,5.35,1.15948);
351 gre->SetPointError(93,0,0,0.0117209,0.0117209);
352 gre->SetPoint(94,5.45,1.1141);
353 gre->SetPointError(94,0,0,0.0114892,0.0114892);
354 gre->SetPoint(95,5.55,1.06315);
355 gre->SetPointError(95,0,0,0.0112235,0.0112235);
356 gre->SetPoint(96,5.65,1.05213);
357 gre->SetPointError(96,0,0,0.0111651,0.0111651);
358 gre->SetPoint(97,5.75,1.02476);
359 gre->SetPointError(97,0,0,0.011019,0.011019);
360 gre->SetPoint(98,5.85,0.984834);
361 gre->SetPointError(98,0,0,0.0108022,0.0108022);
362 gre->SetPoint(99,5.95,0.952844);
363 gre->SetPointError(99,0,0,0.0106253,0.0106253);
364
365 return gre;
366}
367
368//____________________________________________________________________
290052e7 369/**
370 * Get PYTHIA 900GeV NSD data
371 *
372 *
373 * @return Data graph
374 */
b6b35c77 375TGraphAsymmErrors*
376Pythia900NSD()
377{
378
379 TGraphAsymmErrors *gre = new TGraphAsymmErrors(100);
380 SetGraphAttributes(gre, NSD, PYTHIA, false, "pythia900NSD",
381 "Pythia NSD");
382
383 gre->SetPoint(0,-3.95,2.11766);
384 gre->SetPointError(0,0,0,0.0179417,0.0179417);
385 gre->SetPoint(1,-3.85,2.20415);
386 gre->SetPointError(1,0,0,0.0183045,0.0183045);
387 gre->SetPoint(2,-3.75,2.30949);
388 gre->SetPointError(2,0,0,0.0187368,0.0187368);
389 gre->SetPoint(3,-3.65,2.34582);
390 gre->SetPointError(3,0,0,0.0188836,0.0188836);
391 gre->SetPoint(4,-3.55,2.38322);
392 gre->SetPointError(4,0,0,0.0190335,0.0190335);
393 gre->SetPoint(5,-3.45,2.43353);
394 gre->SetPointError(5,0,0,0.0192334,0.0192334);
395 gre->SetPoint(6,-3.35,2.51106);
396 gre->SetPointError(6,0,0,0.0195373,0.0195373);
397 gre->SetPoint(7,-3.25,2.56578);
398 gre->SetPointError(7,0,0,0.0197491,0.0197491);
399 gre->SetPoint(8,-3.15,2.60515);
400 gre->SetPointError(8,0,0,0.0199,0.0199);
401 gre->SetPoint(9,-3.05,2.7105);
402 gre->SetPointError(9,0,0,0.0202984,0.0202984);
403 gre->SetPoint(10,-2.95,2.77008);
404 gre->SetPointError(10,0,0,0.0205203,0.0205203);
405 gre->SetPoint(11,-2.85,2.83332);
406 gre->SetPointError(11,0,0,0.0207532,0.0207532);
407 gre->SetPoint(12,-2.75,2.84715);
408 gre->SetPointError(12,0,0,0.0208038,0.0208038);
409 gre->SetPoint(13,-2.65,2.91693);
410 gre->SetPointError(13,0,0,0.0210571,0.0210571);
411 gre->SetPoint(14,-2.55,2.95797);
412 gre->SetPointError(14,0,0,0.0212048,0.0212048);
413 gre->SetPoint(15,-2.45,2.97499);
414 gre->SetPointError(15,0,0,0.0212657,0.0212657);
415 gre->SetPoint(16,-2.35,3.01345);
416 gre->SetPointError(16,0,0,0.0214027,0.0214027);
417 gre->SetPoint(17,-2.25,3.04659);
418 gre->SetPointError(17,0,0,0.0215201,0.0215201);
419 gre->SetPoint(18,-2.15,3.09341);
420 gre->SetPointError(18,0,0,0.0216848,0.0216848);
421 gre->SetPoint(19,-2.05,3.13187);
422 gre->SetPointError(19,0,0,0.0218192,0.0218192);
423 gre->SetPoint(20,-1.95,3.13917);
424 gre->SetPointError(20,0,0,0.0218446,0.0218446);
425 gre->SetPoint(21,-1.85,3.16911);
426 gre->SetPointError(21,0,0,0.0219485,0.0219485);
427 gre->SetPoint(22,-1.75,3.15665);
428 gre->SetPointError(22,0,0,0.0219053,0.0219053);
429 gre->SetPoint(23,-1.65,3.19693);
430 gre->SetPointError(23,0,0,0.0220446,0.0220446);
431 gre->SetPoint(24,-1.55,3.17002);
432 gre->SetPointError(24,0,0,0.0219517,0.0219517);
433 gre->SetPoint(25,-1.45,3.18538);
434 gre->SetPointError(25,0,0,0.0220048,0.0220048);
435 gre->SetPoint(26,-1.35,3.18066);
436 gre->SetPointError(26,0,0,0.0219885,0.0219885);
437 gre->SetPoint(27,-1.25,3.19754);
438 gre->SetPointError(27,0,0,0.0220467,0.0220467);
439 gre->SetPoint(28,-1.15,3.18021);
440 gre->SetPointError(28,0,0,0.0219869,0.0219869);
441 gre->SetPoint(29,-1.05,3.13111);
442 gre->SetPointError(29,0,0,0.0218165,0.0218165);
443 gre->SetPoint(30,-0.95,3.12153);
444 gre->SetPointError(30,0,0,0.0217831,0.0217831);
445 gre->SetPoint(31,-0.85,3.14798);
446 gre->SetPointError(31,0,0,0.0218752,0.0218752);
447 gre->SetPoint(32,-0.75,3.07912);
448 gre->SetPointError(32,0,0,0.0216347,0.0216347);
449 gre->SetPoint(33,-0.65,3.10207);
450 gre->SetPointError(33,0,0,0.0217151,0.0217151);
451 gre->SetPoint(34,-0.55,3.06346);
452 gre->SetPointError(34,0,0,0.0215796,0.0215796);
453 gre->SetPoint(35,-0.45,2.97651);
454 gre->SetPointError(35,0,0,0.0212711,0.0212711);
455 gre->SetPoint(36,-0.35,2.98715);
456 gre->SetPointError(36,0,0,0.0213091,0.0213091);
457 gre->SetPoint(37,-0.25,2.98548);
458 gre->SetPointError(37,0,0,0.0213032,0.0213032);
459 gre->SetPoint(38,-0.15,3.00555);
460 gre->SetPointError(38,0,0,0.0213746,0.0213746);
461 gre->SetPoint(39,-0.05,3.01193);
462 gre->SetPointError(39,0,0,0.0213973,0.0213973);
463 gre->SetPoint(40,0.05,3.04385);
464 gre->SetPointError(40,0,0,0.0215104,0.0215104);
465 gre->SetPoint(41,0.15,3.04933);
466 gre->SetPointError(41,0,0,0.0215297,0.0215297);
467 gre->SetPoint(42,0.25,3.04659);
468 gre->SetPointError(42,0,0,0.0215201,0.0215201);
469 gre->SetPoint(43,0.35,3.00813);
470 gre->SetPointError(43,0,0,0.0213838,0.0213838);
471 gre->SetPoint(44,0.45,3.06666);
472 gre->SetPointError(44,0,0,0.0215908,0.0215908);
473 gre->SetPoint(45,0.55,3.07167);
474 gre->SetPointError(45,0,0,0.0216085,0.0216085);
475 gre->SetPoint(46,0.65,3.0659);
476 gre->SetPointError(46,0,0,0.0215881,0.0215881);
477 gre->SetPoint(47,0.75,3.09159);
478 gre->SetPointError(47,0,0,0.0216784,0.0216784);
479 gre->SetPoint(48,0.85,3.10846);
480 gre->SetPointError(48,0,0,0.0217375,0.0217375);
481 gre->SetPoint(49,0.95,3.11925);
482 gre->SetPointError(49,0,0,0.0217752,0.0217752);
483 gre->SetPoint(50,1.05,3.15558);
484 gre->SetPointError(50,0,0,0.0219016,0.0219016);
485 gre->SetPoint(51,1.15,3.16911);
486 gre->SetPointError(51,0,0,0.0219485,0.0219485);
487 gre->SetPoint(52,1.25,3.17246);
488 gre->SetPointError(52,0,0,0.0219601,0.0219601);
489 gre->SetPoint(53,1.35,3.19146);
490 gre->SetPointError(53,0,0,0.0220258,0.0220258);
491 gre->SetPoint(54,1.45,3.17458);
492 gre->SetPointError(54,0,0,0.0219675,0.0219675);
493 gre->SetPoint(55,1.55,3.16866);
494 gre->SetPointError(55,0,0,0.0219469,0.0219469);
495 gre->SetPoint(56,1.65,3.16592);
496 gre->SetPointError(56,0,0,0.0219375,0.0219375);
497 gre->SetPoint(57,1.75,3.16394);
498 gre->SetPointError(57,0,0,0.0219306,0.0219306);
499 gre->SetPoint(58,1.85,3.11956);
500 gre->SetPointError(58,0,0,0.0217762,0.0217762);
501 gre->SetPoint(59,1.95,3.14646);
502 gre->SetPointError(59,0,0,0.02187,0.02187);
503 gre->SetPoint(60,2.05,3.10147);
504 gre->SetPointError(60,0,0,0.021713,0.021713);
505 gre->SetPoint(61,2.15,3.09356);
506 gre->SetPointError(61,0,0,0.0216853,0.0216853);
507 gre->SetPoint(62,2.25,3.05328);
508 gre->SetPointError(62,0,0,0.0215437,0.0215437);
509 gre->SetPoint(63,2.35,3.01953);
510 gre->SetPointError(63,0,0,0.0214243,0.0214243);
511 gre->SetPoint(64,2.45,2.9373);
512 gre->SetPointError(64,0,0,0.0211305,0.0211305);
513 gre->SetPoint(65,2.55,2.92772);
514 gre->SetPointError(65,0,0,0.0210961,0.0210961);
515 gre->SetPoint(66,2.65,2.89154);
516 gre->SetPointError(66,0,0,0.0209653,0.0209653);
517 gre->SetPoint(67,2.75,2.87619);
518 gre->SetPointError(67,0,0,0.0209096,0.0209096);
519 gre->SetPoint(68,2.85,2.78924);
520 gre->SetPointError(68,0,0,0.0205911,0.0205911);
521 gre->SetPoint(69,2.95,2.72159);
522 gre->SetPointError(69,0,0,0.0203399,0.0203399);
523 gre->SetPoint(70,3.05,2.69089);
524 gre->SetPointError(70,0,0,0.0202248,0.0202248);
525 gre->SetPoint(71,3.15,2.64939);
526 gre->SetPointError(71,0,0,0.0200682,0.0200682);
527 gre->SetPoint(72,3.25,2.55545);
528 gre->SetPointError(72,0,0,0.0197092,0.0197092);
529 gre->SetPoint(73,3.35,2.56745);
530 gre->SetPointError(73,0,0,0.0197555,0.0197555);
531 gre->SetPoint(74,3.45,2.47503);
532 gre->SetPointError(74,0,0,0.0193967,0.0193967);
533 gre->SetPoint(75,3.55,2.36741);
534 gre->SetPointError(75,0,0,0.0189703,0.0189703);
535 gre->SetPoint(76,3.65,2.33412);
536 gre->SetPointError(76,0,0,0.0188364,0.0188364);
537 gre->SetPoint(77,3.75,2.2385);
538 gre->SetPointError(77,0,0,0.0184466,0.0184466);
539 gre->SetPoint(78,3.85,2.21768);
540 gre->SetPointError(78,0,0,0.0183606,0.0183606);
541 gre->SetPoint(79,3.95,2.1055);
542 gre->SetPointError(79,0,0,0.0178901,0.0178901);
543 gre->SetPoint(80,4.05,2.05047);
544 gre->SetPointError(80,0,0,0.0176548,0.0176548);
545 gre->SetPoint(81,4.15,2.00486);
546 gre->SetPointError(81,0,0,0.0174574,0.0174574);
547 gre->SetPoint(82,4.25,1.94573);
548 gre->SetPointError(82,0,0,0.017198,0.017198);
549 gre->SetPoint(83,4.35,1.87064);
550 gre->SetPointError(83,0,0,0.0168629,0.0168629);
551 gre->SetPoint(84,4.45,1.83735);
552 gre->SetPointError(84,0,0,0.0167122,0.0167122);
553 gre->SetPoint(85,4.55,1.75314);
554 gre->SetPointError(85,0,0,0.0163247,0.0163247);
555 gre->SetPoint(86,4.65,1.65828);
556 gre->SetPointError(86,0,0,0.0158769,0.0158769);
557 gre->SetPoint(87,4.75,1.60751);
558 gre->SetPointError(87,0,0,0.015632,0.015632);
559 gre->SetPoint(88,4.85,1.56312);
560 gre->SetPointError(88,0,0,0.0154146,0.0154146);
561 gre->SetPoint(89,4.95,1.52117);
562 gre->SetPointError(89,0,0,0.0152064,0.0152064);
563 gre->SetPoint(90,5.05,1.46553);
564 gre->SetPointError(90,0,0,0.0149257,0.0149257);
565 gre->SetPoint(91,5.15,1.42038);
566 gre->SetPointError(91,0,0,0.014694,0.014694);
567 gre->SetPoint(92,5.25,1.38816);
568 gre->SetPointError(92,0,0,0.0145263,0.0145263);
569 gre->SetPoint(93,5.35,1.35046);
570 gre->SetPointError(93,0,0,0.0143277,0.0143277);
571 gre->SetPoint(94,5.45,1.30075);
572 gre->SetPointError(94,0,0,0.0140616,0.0140616);
573 gre->SetPoint(95,5.55,1.24025);
574 gre->SetPointError(95,0,0,0.0137307,0.0137307);
575 gre->SetPoint(96,5.65,1.21806);
576 gre->SetPointError(96,0,0,0.0136073,0.0136073);
577 gre->SetPoint(97,5.75,1.19435);
578 gre->SetPointError(97,0,0,0.0134742,0.0134742);
579 gre->SetPoint(98,5.85,1.14175);
580 gre->SetPointError(98,0,0,0.0131741,0.0131741);
581 gre->SetPoint(99,5.95,1.09235);
582 gre->SetPointError(99,0,0,0.012886,0.012886);
583
584 return gre;
585}
586
7e4038b5 587//____________________________________________________________________
588/**
589 * Get the UA5 NSD data for pp at @f$ \sqrt{s} = 900GeV@f$
590 * p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
591 *
592 * @param mirrored Wether to produce the mirrored plot
593 *
594 * @return graph of data
595 *
bd6f5206 596 * @ingroup pwglf_forward_otherdata
7e4038b5 597 */
598TGraphAsymmErrors* UA5Nsd(Bool_t mirrored=false)
599{
600 //UA5 data NSD - p7886_d1x1y4 - Z.Phys.C33:1-6,1986.
601 double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
602 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
603 4.625 };
604 double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
605 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
606 double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
607 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
608 double y[] = { 3.48, 3.38, 3.52, 3.68, 3.71, 3.86, 3.76, 3.66, 3.72,
609 3.69, 3.56, 3.41, 3.15, 3.09, 2.74, 2.73, 2.32, 1.99, 1.69 };
610 double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
611 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
612 double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
613 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
860e603e 614 const int np = 19;
7e4038b5 615 double xm[np];
b2e7f2d6 616 double exmm[np];
617 double expm[np];
618 double ym[np];
619 double eymm[np];
620 double eypm[np];
621 for (int i = 0; i < np; i++) {
622 int j = np-1-i;
623 xm[i] = -x[j];
624 exmm[i] = exm[j];
625 expm[i] = exp[j];
626 ym[i] = y[j];
627 eymm[i] = eym[j];
628 eypm[i] = eyp[j];
629 }
7e4038b5 630
b2e7f2d6 631 TGraphAsymmErrors* g = new TGraphAsymmErrors(19,x, y, exm, exp, eym, eyp);
632 TGraphAsymmErrors* gm = new TGraphAsymmErrors(19,xm,ym,exmm,expm,eymm,eypm);
ffca499d 633 SetGraphAttributes(g, NSD, UA5, false,"ua5_nsd", "UA5 NSD");
634 SetGraphAttributes(gm, NSD, UA5, true,"ua5_nsd_mirrored",
7e4038b5 635 "UA5 NSD (mirrored)");
636
637 return (mirrored ? gm : g);
638}
639
640//____________________________________________________________________
641/**
642 * Get the UA5 INEL data for pp at @f$ \sqrt{s} = 900GeV@f$
643 * p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
644 *
645 * @param mirrored Wether to produce the mirrored plot
646 *
647 * @return graph of data
648 *
bd6f5206 649 * @ingroup pwglf_forward_otherdata
7e4038b5 650 */
651TGraphAsymmErrors* UA5Inel(Bool_t mirrored=false)
652{
653 //UA5 data INEL - p7886_d2x1y2 - Z.Phys.C33:1-6,1986.
654 double x[] = { 0.125, 0.375, 0.625, 0.875, 1.125, 1.375, 1.625, 1.875, 2.125,
655 2.375, 2.625, 2.875, 3.125, 3.375, 3.625, 3.875, 4.125, 4.375,
656 4.625 };
657 double exm[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
658 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
659 double exp[] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
660 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 };
661 double y[] = { 3.14, 3.04, 3.17, 3.33, 3.33, 3.53, 3.46, 3.41, 3.45,
662 3.39, 3.07, 3.07, 2.93, 2.93, 2.55, 2.48, 2.18, 1.91, 1.52 };
663 double eym[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
664 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
665 double eyp[] = { 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07, 0.07,
666 0.07, 0.07, 0.07, 0.07, 0.08, 0.08, 0.09, 0.09, 0.1, 0.13 };
667 const int np = 19;
668 double xm[np];
b2e7f2d6 669 double exmm[np];
670 double expm[np];
671 double ym[np];
672 double eymm[np];
673 double eypm[np];
674 for (int i = 0; i < np; i++) {
675 int j = np-1-i;
676 xm[i] = -x[j];
677 exmm[i] = exm[j];
678 expm[i] = exp[j];
679 ym[i] = y[j];
680 eymm[i] = eym[j];
681 eypm[i] = eyp[j];
682 }
683 TGraphAsymmErrors* g = new TGraphAsymmErrors(np,x, y, exm, exp, eym, eyp);
684 TGraphAsymmErrors* gm = new TGraphAsymmErrors(np,xm,ym,exmm,expm,eymm,eypm);
7e4038b5 685
ffca499d 686 SetGraphAttributes(g, INEL, UA5, false, "ua5_inel", "UA5 INEL");
687 SetGraphAttributes(gm, INEL, UA5, true, "ua5_inel_mirrored",
7e4038b5 688 "UA5 INEL (mirrored)");
689
690 return (mirrored ? gm : g);
691}
692
693//____________________________________________________________________
694/**
695 * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$ \sqrt{s}
696 * = 900GeV@f$
697 * p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
698 *
699 * @return graph of data
700 *
bd6f5206 701 * @ingroup pwglf_forward_otherdata
7e4038b5 702 */
703TGraphAsymmErrors* AliceCentralInel900()
704{
705 // INEL - p7742_d1x1y1 - Eur.Phys.J.C68:89-108,2010.
706 TGraphAsymmErrors* g = 0;
707 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3,
708 -0.1, 0.1, 0.3, 0.5, 0.7, 0.9, 1.1, 1.3 };
709 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
710 0.1, 0.1, 0.1, 0.1, 0.1 };
711 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
712 0.1, 0.1, 0.1, 0.1, 0.1 };
713 double y[] = { 3.28, 3.28, 3.22, 3.12, 3.06, 3.02, 2.98, 3.02, 3.02,
714 3.05, 3.15, 3.21, 3.26, 3.33 };
715 double eym[] = { 0.06324555320336758, 0.06324555320336758,
716 0.06324555320336758, 0.06324555320336758,
717 0.06324555320336758, 0.05385164807134505,
718 0.05385164807134505, 0.05385164807134505,
719 0.05385164807134505, 0.06324555320336758,
720 0.06324555320336758, 0.06324555320336758,
721 0.06324555320336758, 0.06324555320336758 };
722 double eyp[] = { 0.08246211251235322, 0.08246211251235322,
723 0.08246211251235322, 0.08246211251235322,
724 0.08246211251235322, 0.08246211251235322,
725 0.07280109889280519, 0.08246211251235322,
726 0.08246211251235322, 0.08246211251235322,
727 0.08246211251235322, 0.08246211251235322,
728 0.08246211251235322, 0.08246211251235322 };
729 const int np = 14;
730 for (int i = 0; i < np; i++) {
731 eym[i] += 0.02;
732 eyp[i] += 0.02;
733 }
734 g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 735 SetGraphAttributes(g, INEL, ALICE, false, "alice_inel",
bfab35d9 736 "INEL - Eur.Phys.J.C68:89-108,2010"
737 /* "ALICE INEL (publ.)" */);
7e4038b5 738
739 return g;
740}
741
742//____________________________________________________________________
743/**
744 * Get the ALICE INEL>0 data in @f$ |\eta|<1.3@f$ for pp at @f$
745 * \sqrt{s} = 900GeV@f$
746 *
747 * p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
748 *
749 * @return graph of data
750 *
bd6f5206 751 * @ingroup pwglf_forward_otherdata
7e4038b5 752 */
753TGraphAsymmErrors* AliceCentralInelGt900()
754{
797161e8 755#if 0
7e4038b5 756 // INEL>0 - p7741_d4x1y1 - Eur.Phys.J.C68:345-354,2010.
757 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7,
758 0.9 };
759 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
760 0.1 };
761 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
762 0.1 };
763 double y[] = { 4.0, 3.87, 3.8, 3.7, 3.67, 3.73, 3.72, 3.77, 3.92,
764 4.01 };
765 double eym[] = { 0.07615773105863909, 0.07615773105863909,
766 0.07615773105863909, 0.06324555320336758,
767 0.06324555320336758, 0.06324555320336758,
768 0.0670820393249937, 0.07615773105863909,
769 0.07615773105863909, 0.07615773105863909 };
770 double eyp[] = { 0.08544003745317531, 0.07615773105863909,
771 0.07615773105863909, 0.07280109889280519,
772 0.07280109889280519, 0.07280109889280519,
773 0.07615773105863909, 0.07615773105863909,
774 0.08544003745317531, 0.08544003745317531 };
775 const int np = 10;
776 for (int i = 0; i < np; i++) {
777 double stat = (i >= 3 && i<=5) ? 0.02 : 0.03;
778 eym[i] += stat;
779 eyp[i] += stat;
780 }
781
782 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
797161e8 783#else
784 // These are from JFGO
785 TGraphAsymmErrors *g = new TGraphAsymmErrors(15);
786 g->SetPoint(0,-1.5,4.12575);
787 g->SetPointError(0,0.1,0.1,0.0742967,0.0814571);
788 g->SetPoint(1,-1.3,3.91209);
789 g->SetPointError(1,0.1,0.1,0.0697701,0.0766199);
790 g->SetPoint(2,-1.1,3.98377);
791 g->SetPointError(2,0.1,0.1,0.0704503,0.0774795);
792 g->SetPoint(3,-0.9,4.00035);
793 g->SetPointError(3,0.1,0.1,0.0702388,0.0773433);
794 g->SetPoint(4,-0.7,3.87228);
795 g->SetPointError(4,0.1,0.1,0.067597,0.0745103);
796 g->SetPoint(5,-0.5,3.79613);
797 g->SetPointError(5,0.1,0.1,0.0659771,0.0727816);
798 g->SetPoint(6,-0.3,3.70489);
799 g->SetPointError(6,0.1,0.1,0.0642016,0.0708603);
800 g->SetPoint(7,-0.1,3.67423);
801 g->SetPointError(7,0.1,0.1,0.0635759,0.0701884);
802 g->SetPoint(8,0.1,3.72765);
803 g->SetPointError(8,0.1,0.1,0.0645004,0.071209);
804 g->SetPoint(9,0.3,3.72171);
805 g->SetPointError(9,0.1,0.1,0.064493,0.071182);
806 g->SetPoint(10,0.5,3.77428);
807 g->SetPointError(10,0.1,0.1,0.0655974,0.0723627);
808 g->SetPoint(11,0.7,3.91704);
809 g->SetPointError(11,0.1,0.1,0.0683783,0.0753716);
810 g->SetPoint(12,0.9,4.00674);
811 g->SetPointError(12,0.1,0.1,0.0703511,0.0774669);
812 g->SetPoint(13,1.1,3.97948);
813 g->SetPointError(13,0.1,0.1,0.0703744,0.077396);
814 g->SetPoint(14,1.3,3.99165);
815 g->SetPointError(14,0.1,0.1,0.0711888,0.078178);
816#endif
ffca499d 817 SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt900",
bfab35d9 818 // "ALICE INEL>0 (publ.)"
819 "INEL>0 - Eur.Phys.J.C68:345-354,2010");
7e4038b5 820 return g;
7e4038b5 821}
822
823//____________________________________________________________________
824/**
825 * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
826 * \sqrt{s} = 2.36TeV@f$
827 *
828 * p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
829 *
830 * @return graph of data
831 *
bd6f5206 832 * @ingroup pwglf_forward_otherdata
7e4038b5 833 */
834TGraphAsymmErrors* AliceCentralInelGt2360()
835{
797161e8 836#if 0
7e4038b5 837 // INEL>0 - p7741_d5x1y1 - Eur.Phys.J.C68:345-354,2010.
838 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
839 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
840 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
841 double y[] = { 4.91, 4.76, 4.63, 4.64, 4.55, 4.55, 4.64, 4.66, 4.82, 4.88 };
842 double eym[] = { 0.08544003745317531, 0.08544003745317531,
843 0.08544003745317531, 0.08544003745317531,
844 0.08544003745317531, 0.08544003745317531,
845 0.08544003745317531, 0.08544003745317531,
846 0.08544003745317531, 0.08544003745317531 };
847 double eyp[] = { 0.11401754250991379, 0.11401754250991379,
848 0.1044030650891055, 0.1044030650891055,
849 0.1044030650891055, 0.1044030650891055,
850 0.1044030650891055, 0.1044030650891055,
851 0.11401754250991379, 0.11401754250991379 };
852 const int np = 10;
853 for (int i = 0; i < np; i++) {
854 double stat = 0.3;
855 eym[i] += stat;
856 eyp[i] += stat;
857 }
858
859 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
797161e8 860#else
861 // These are from JFGO
862 TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
863 g->SetPoint(0,-1.5,4.79047);
864 g->SetPointError(0,0.1,0.1,0.0844278,0.109947);
865 g->SetPoint(1,-1.3,4.91068);
866 g->SetPointError(1,0.1,0.1,0.0856751,0.112038);
867 g->SetPoint(2,-1.1,4.87386);
868 g->SetPointError(2,0.1,0.1,0.0842846,0.110628);
869 g->SetPoint(3,-0.9,4.91365);
870 g->SetPointError(3,0.1,0.1,0.084339,0.111049);
871 g->SetPoint(4,-0.7,4.7601);
872 g->SetPointError(4,0.1,0.1,0.0812087,0.107203);
873 g->SetPoint(5,-0.5,4.63355);
874 g->SetPointError(5,0.1,0.1,0.078687,0.104079);
875 g->SetPoint(6,-0.3,4.63885);
876 g->SetPointError(6,0.1,0.1,0.0785337,0.104014);
877 g->SetPoint(7,-0.1,4.55439);
878 g->SetPointError(7,0.1,0.1,0.0769842,0.10203);
879 g->SetPoint(8,0.1,4.55087);
880 g->SetPointError(8,0.1,0.1,0.0769246,0.101951);
881 g->SetPoint(9,0.3,4.64118);
882 g->SetPointError(9,0.1,0.1,0.0785732,0.104066);
883 g->SetPoint(10,0.5,4.66172);
884 g->SetPointError(10,0.1,0.1,0.0791652,0.104711);
885 g->SetPoint(11,0.7,4.81871);
886 g->SetPointError(11,0.1,0.1,0.0822086,0.108523);
887 g->SetPoint(12,0.9,4.88193);
888 g->SetPointError(12,0.1,0.1,0.0837944,0.110332);
889 g->SetPoint(13,1.1,4.89068);
890 g->SetPointError(13,0.1,0.1,0.0845754,0.111009);
891 g->SetPoint(14,1.3,5.05663);
892 g->SetPointError(14,0.1,0.1,0.0882216,0.115368);
893#endif
894
ffca499d 895 SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt2360",
bfab35d9 896 // "ALICE INEL>0 (publ.)"
897 "INEL>0 - Eur.Phys.J.C68:345-354,2010");
7e4038b5 898 return g;
899}
900
901//____________________________________________________________________
902/**
903 * Get the ALICE INEL>0 data in @f$ |\eta|<0.9@f$ for pp at @f$
904 * \sqrt{s} = 7TeV@f$
905 *
906 * p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
907 *
908 * @return graph of data
909 *
bd6f5206 910 * @ingroup pwglf_forward_otherdata
7e4038b5 911 */
912TGraphAsymmErrors* AliceCentralInelGt7000()
913{
797161e8 914#if 0
7e4038b5 915 // INEL>0 - p7741_d6x1y1 - Eur.Phys.J.C68:345-354,2010.
916// Plot: p7741_d6x1y1
917 double x[] = { -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3, 0.5, 0.7, 0.9 };
918 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
919 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1 };
920 double y[] = { 6.22, 6.07, 6.01, 5.84, 5.85, 5.85, 5.91, 6.01, 6.17, 6.26 };
921 double eym[] = { 0.1216552506059644, 0.1216552506059644,
922 0.1216552506059644, 0.11180339887498948,
923 0.11180339887498948, 0.11180339887498948,
924 0.11180339887498948, 0.1216552506059644,
925 0.1216552506059644, 0.1216552506059644 };
926 double eyp[] = { 0.21095023109728983, 0.21095023109728983,
927 0.2009975124224178, 0.2009975124224178,
928 0.2009975124224178, 0.2009975124224178,
929 0.2009975124224178, 0.2009975124224178,
930 0.21095023109728983, 0.21095023109728983 };
931 const int np = 10;
932 for (int i = 0; i < np; i++) {
933 double stat = 0.2;
934 eym[i] += stat;
935 eyp[i] += stat;
936 }
937
938 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
797161e8 939#else
940 // These are from JFGO
941 TGraphAsymmErrors* g = new TGraphAsymmErrors(15);
942 g->SetPoint(0,-1.5,6.28573);
943 g->SetPointError(0,0.1,0.1,0.125928,0.215392);
944 g->SetPoint(1,-1.3,6.25573);
945 g->SetPointError(1,0.1,0.1,0.124352,0.213795);
946 g->SetPoint(2,-1.1,6.28779);
947 g->SetPointError(2,0.1,0.1,0.124143,0.214399);
948 g->SetPoint(3,-0.9,6.21881);
949 g->SetPointError(3,0.1,0.1,0.122079,0.211642);
950 g->SetPoint(4,-0.7,6.0728);
951 g->SetPointError(4,0.1,0.1,0.118661,0.206355);
952 g->SetPoint(5,-0.5,6.011);
953 g->SetPointError(5,0.1,0.1,0.117043,0.204019);
954 g->SetPoint(6,-0.3,5.84071);
955 g->SetPointError(6,0.1,0.1,0.11346,0.198086);
956 g->SetPoint(7,-0.1,5.8532);
957 g->SetPointError(7,0.1,0.1,0.113569,0.198433);
958 g->SetPoint(8,0.1,5.84811);
959 g->SetPointError(8,0.1,0.1,0.11347,0.198261);
960 g->SetPoint(9,0.3,5.91022);
961 g->SetPointError(9,0.1,0.1,0.11481,0.200444);
962 g->SetPoint(10,0.5,6.00649);
963 g->SetPointError(10,0.1,0.1,0.116955,0.203866);
964 g->SetPoint(11,0.7,6.17115);
965 g->SetPointError(11,0.1,0.1,0.120583,0.209697);
966 g->SetPoint(12,0.9,6.2645);
967 g->SetPointError(12,0.1,0.1,0.122976,0.213197);
968 g->SetPoint(13,1.1,6.36448);
969 g->SetPointError(13,0.1,0.1,0.125657,0.217014);
970 g->SetPoint(14,1.3,6.39489);
971 g->SetPointError(14,0.1,0.1,0.127118,0.218551);
972#endif
ffca499d 973 SetGraphAttributes(g, INELGt0, ALICE, false, "alice_inelgt7000",
bfab35d9 974 // "ALICE INEL>0 (publ.)"
975 "INEL > 0 - Eur.Phys.J.C68:345-354,2010");
7e4038b5 976 return g;
977}
978
979//____________________________________________________________________
980/**
981 * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
982 * \sqrt{s} = 900GeV@f$
983 *
984 * p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
985 *
986 * @return graph of data
987 *
bd6f5206 988 * @ingroup pwglf_forward_otherdata
7e4038b5 989 */
990TGraphAsymmErrors* AliceCentralNsd900()
991{
992 //NSD - p7742_d2x1y1 - Eur.Phys.J.C68:89-108,2010.
993 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
994 0.5, 0.7, 0.9, 1.1, 1.3 };
995 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
996 0.1, 0.1, 0.1, 0.1, 0.1 };
997 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
998 0.1, 0.1, 0.1, 0.1, 0.1 };
999 double y[] = { 3.9, 3.89, 3.81, 3.7, 3.64, 3.59, 3.53, 3.58, 3.59,
1000 3.61, 3.74, 3.8, 3.87, 3.95 };
1001 double eym[] = { 0.13341664064126335, 0.13152946437965907,
1002 0.13152946437965907, 0.1216552506059644,
1003 0.1216552506059644, 0.1216552506059644,
1004 0.1216552506059644, 0.1216552506059644,
1005 0.1216552506059644, 0.1216552506059644,
1006 0.1216552506059644, 0.13152946437965907,
1007 0.13152946437965907, 0.13341664064126335 };
1008 double eyp[] = { 0.13341664064126335, 0.13152946437965907,
1009 0.13152946437965907, 0.1216552506059644,
1010 0.1216552506059644, 0.1216552506059644,
1011 0.1216552506059644, 0.1216552506059644,
1012 0.1216552506059644, 0.1216552506059644,
1013 0.1216552506059644, 0.13152946437965907,
1014 0.13152946437965907, 0.13341664064126335 };
1015 const int np = 14;
1016 for (int i = 0; i < np; i++) {
1017 double stat = (i == 0 || i == np-1) ? 0.03 : 0.02;
1018 eym[i] += stat;
1019 eyp[i] += stat;
1020 }
1021
1022 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
bfab35d9 1023 SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd",
1024 //"ALICE NSD (publ.)"
1025 "NSD - Eur.Phys.J.C68:89-108,2010");
7e4038b5 1026
1027 return g;
1028}
1029
1030//____________________________________________________________________
1031/**
1032 * Get the ALICE INEL data in @f$ |\eta|<1.3@f$ for pp at @f$
1033 * \sqrt{s} = 2.36TeV@f$
1034 *
1035 * p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
1036 *
1037 * @return graph of data
1038 *
bd6f5206 1039 * @ingroup pwglf_forward_otherdata
7e4038b5 1040 */
1041TGraphAsymmErrors* AliceCentralInel2360()
1042{
1043 // INEL - p7742_d3x1y1 - Eur.Phys.J.C68:89-108,2010.
1044 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
1045 0.5, 0.7, 0.9, 1.1, 1.3 };
1046 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
1047 0.1, 0.1, 0.1, 0.1, 0.1 };
1048 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
1049 0.1, 0.1, 0.1, 0.1, 0.1 };
1050 double y[] = { 4.08, 4.01, 4.0, 3.88, 3.77, 3.8, 3.73, 3.71, 3.79,
1051 3.82, 3.94, 4.02, 4.05, 4.16 };
1052 double eym[] = { 0.13341664064126335, 0.12369316876852982,
1053 0.12369316876852982, 0.1216552506059644,
1054 0.1216552506059644, 0.1216552506059644,
1055 0.1216552506059644, 0.11180339887498948,
1056 0.1216552506059644, 0.1216552506059644,
1057 0.12369316876852982, 0.12369316876852982,
1058 0.13341664064126335, 0.13341664064126335 };
1059 double eyp[] = { 0.2716615541441225, 0.2716615541441225,
1060 0.2716615541441225, 0.260768096208106,
1061 0.25079872407968906, 0.25079872407968906,
1062 0.25079872407968906, 0.25079872407968906,
1063 0.25079872407968906, 0.260768096208106,
1064 0.261725046566048, 0.2716615541441225,
1065 0.2716615541441225, 0.2816025568065745 };
1066 const int np = 14;
1067 for (int i = 0; i < np; i++) {
1068 double stat = (i < 3 || i > np-1-4) ? 0.03 : 0.02;
1069 eym[i] += stat;
1070 eyp[i] += stat;
1071 }
1072
1073 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 1074 SetGraphAttributes(g, NSD, ALICE, false, "alice_inel2360",
bfab35d9 1075 // "ALICE INEL (publ.)"
1076 "INEL - Eur.Phys.J.C68:89-108,2010");
7e4038b5 1077 return g;
1078}
1079
1080//____________________________________________________________________
1081/**
1082 * Get the ALICE NSD data in @f$ |\eta|<1.3@f$ for pp at @f$
1083 * \sqrt{s} = 2.36TeV@f$
1084 *
1085 * p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
1086 *
1087 * @return graph of data
1088 *
bd6f5206 1089 * @ingroup pwglf_forward_otherdata
7e4038b5 1090 */
1091TGraphAsymmErrors* AliceCentralNsd2360()
1092{
1093 // NSD - p7742_d4x1y1 - Eur.Phys.J.C68:89-108,2010.
1094 double x[] = { -1.3, -1.1, -0.9, -0.7, -0.5, -0.3, -0.1, 0.1, 0.3,
1095 0.5, 0.7, 0.9, 1.1, 1.3 };
1096 double exm[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
1097 0.1, 0.1, 0.1, 0.1, 0.1 };
1098 double exp[] = { 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1,
1099 0.1, 0.1, 0.1, 0.1, 0.1 };
1100 double y[] = { 4.79, 4.72, 4.7, 4.56, 4.44, 4.47, 4.39, 4.37, 4.47,
1101 4.5, 4.64, 4.73, 4.76, 4.9 };
1102 double eym[] = { 0.13601470508735444, 0.13341664064126335,
1103 0.13341664064126335, 0.12369316876852982,
1104 0.12369316876852982, 0.12369316876852982,
1105 0.12369316876852982, 0.12369316876852982,
1106 0.12369316876852982, 0.12369316876852982,
1107 0.12369316876852982, 0.13341664064126335,
1108 0.13341664064126335, 0.13341664064126335 };
1109 double eyp[] = { 0.18439088914585774, 0.18248287590894657,
1110 0.18248287590894657, 0.1726267650163207,
1111 0.1726267650163207, 0.1726267650163207,
1112 0.16278820596099708, 0.16278820596099708,
1113 0.1726267650163207, 0.1726267650163207,
1114 0.1726267650163207, 0.18248287590894657,
1115 0.18248287590894657, 0.18248287590894657 };
1116 const int np = 14;
1117
1118 for (int i = 0; i < np; i++) {
1119 double stat = (i < 1) ? 0.03 : 0.02;
1120 eym[i] += stat;
1121 eyp[i] += stat;
1122 }
1123
1124 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 1125 SetGraphAttributes(g, NSD, ALICE, false, "alice_nsd2360",
bfab35d9 1126 // "ALICE NSD (publ.)"
1127 "NSD - Eur.Phys.J.C68:89-108,2010");
7e4038b5 1128 return g;
1129}
1130
6d233add 1131//____________________________________________________________________
1132/**
1133 * Get the ALICE MB (CINT5, V0AND) data in @f$ |\eta|<2@f$ for pPb at
1134 * @f$ \sqrt{s} = 5.023TeV@f$
1135 *
1136 * arXiv:1210.3615 [nucl-ex]
1137 *
1138 * @return graph of data
1139 *
1140 * @ingroup pwglf_forward_otherdata
1141 */
1142TGraphAsymmErrors* AliceCentralpPb5023()
1143{
1144 TGraphAsymmErrors *g = new TGraphAsymmErrors(43);
1145 g->SetName("");
1146 g->SetTitle("");
1147 g->SetFillColor(1);
1148 g->SetMarkerStyle(24);
1149 g->SetMarkerSize(1.3);
1150 g->SetPoint(0,-2.05,16.1757);
1151 g->SetPointError(0,0.05,0.05,0.360241,0.360241);
1152 g->SetPoint(1,-1.95,16.8454);
1153 g->SetPointError(1,0.05,0.05,0.184436,0.184436);
1154 g->SetPoint(2,-1.85,16.9748);
1155 g->SetPointError(2,0.05,0.05,0.112833,0.112833);
1156 g->SetPoint(3,-1.75,17.2504);
1157 g->SetPointError(3,0.05,0.05,0.0851792,0.0851792);
1158 g->SetPoint(4,-1.65,17.336);
1159 g->SetPointError(4,0.05,0.05,0.0674453,0.0674453);
1160 g->SetPoint(5,-1.55,17.4649);
1161 g->SetPointError(5,0.05,0.05,0.0558537,0.0558537);
1162 g->SetPoint(6,-1.45,17.5884);
1163 g->SetPointError(6,0.05,0.05,0.0481065,0.0481065);
1164 g->SetPoint(7,-1.35,17.7509);
1165 g->SetPointError(7,0.05,0.05,0.0448122,0.0448122);
1166 g->SetPoint(8,-1.25,17.8052);
1167 g->SetPointError(8,0.05,0.05,0.0403298,0.0403298);
1168 g->SetPoint(9,-1.15,17.8357);
1169 g->SetPointError(9,0.05,0.05,0.0383174,0.0383174);
1170 g->SetPoint(10,-1.05,17.7547);
1171 g->SetPointError(10,0.05,0.05,0.0356689,0.0356689);
1172 g->SetPoint(11,-0.95,17.6859);
1173 g->SetPointError(11,0.05,0.05,0.034326,0.034326);
1174 g->SetPoint(12,-0.85,17.6665);
1175 g->SetPointError(12,0.05,0.05,0.0333449,0.0333449);
1176 g->SetPoint(13,-0.75,17.6044);
1177 g->SetPointError(13,0.05,0.05,0.0325356,0.0325356);
1178 g->SetPoint(14,-0.65,17.4815);
1179 g->SetPointError(14,0.05,0.05,0.0318004,0.0318004);
1180 g->SetPoint(15,-0.55,17.4);
1181 g->SetPointError(15,0.05,0.05,0.0312675,0.0312675);
1182 g->SetPoint(16,-0.45,17.3425);
1183 g->SetPointError(16,0.05,0.05,0.0310344,0.0310344);
1184 g->SetPoint(17,-0.35,17.2885);
1185 g->SetPointError(17,0.05,0.05,0.0306043,0.0306043);
1186 g->SetPoint(18,-0.25,17.2646);
1187 g->SetPointError(18,0.05,0.05,0.0303226,0.0303226);
1188 g->SetPoint(19,-0.15,17.316);
1189 g->SetPointError(19,0.05,0.05,0.0302368,0.0302368);
1190 g->SetPoint(20,-0.05,17.312);
1191 g->SetPointError(20,0.05,0.05,0.0301444,0.0301444);
1192 g->SetPoint(21,0.05,17.4418);
1193 g->SetPointError(21,0.05,0.05,0.0301526,0.0301526);
1194 g->SetPoint(22,0.15,17.4944);
1195 g->SetPointError(22,0.05,0.05,0.0303199,0.0303199);
1196 g->SetPoint(23,0.25,17.642);
1197 g->SetPointError(23,0.05,0.05,0.0303867,0.0303867);
1198 g->SetPoint(24,0.35,17.8153);
1199 g->SetPointError(24,0.05,0.05,0.0306752,0.0306752);
1200 g->SetPoint(25,0.45,18.0244);
1201 g->SetPointError(25,0.05,0.05,0.0310274,0.0310274);
1202 g->SetPoint(26,0.55,18.1993);
1203 g->SetPointError(26,0.05,0.05,0.0314353,0.0314353);
1204 g->SetPoint(27,0.65,18.349);
1205 g->SetPointError(27,0.05,0.05,0.0316803,0.0316803);
1206 g->SetPoint(28,0.75,18.5976);
1207 g->SetPointError(28,0.05,0.05,0.0322819,0.0322819);
1208 g->SetPoint(29,0.85,18.8045);
1209 g->SetPointError(29,0.05,0.05,0.0329447,0.0329447);
1210 g->SetPoint(30,0.95,18.9865);
1211 g->SetPointError(30,0.05,0.05,0.0337513,0.0337513);
1212 g->SetPoint(31,1.05,19.2313);
1213 g->SetPointError(31,0.05,0.05,0.0354009,0.0354009);
1214 g->SetPoint(32,1.15,19.4055);
1215 g->SetPointError(32,0.05,0.05,0.0367366,0.0367366);
1216 g->SetPoint(33,1.25,19.5893);
1217 g->SetPointError(33,0.05,0.05,0.0385048,0.0385048);
1218 g->SetPoint(34,1.35,19.8196);
1219 g->SetPointError(34,0.05,0.05,0.0421699,0.0421699);
1220 g->SetPoint(35,1.45,19.9476);
1221 g->SetPointError(35,0.05,0.05,0.0451541,0.0451541);
1222 g->SetPoint(36,1.55,20.1012);
1223 g->SetPointError(36,0.05,0.05,0.0513641,0.0513641);
1224 g->SetPoint(37,1.65,20.1082);
1225 g->SetPointError(37,0.05,0.05,0.060302,0.060302);
1226 g->SetPoint(38,1.75,20.1732);
1227 g->SetPointError(38,0.05,0.05,0.0739969,0.0739969);
1228 g->SetPoint(39,1.85,20.1964);
1229 g->SetPointError(39,0.05,0.05,0.0953757,0.0953757);
1230 g->SetPoint(40,1.95,20.0509);
1231 g->SetPointError(40,0.05,0.05,0.147212,0.147212);
1232 g->SetPoint(41,2.05,20.3151);
1233 g->SetPointError(41,0.05,0.05,0.272151,0.272151);
1234 g->SetPoint(42,2.15,20.1319);
1235 g->SetPointError(42,0.05,0.05,0.802706,0.802706);
1236
1237
bfab35d9 1238 SetGraphAttributes(g, NSD, ALICE, false, "alice_ppb50230",
1239 "NSD arXiv:1210.3615");
1240 return g;
1241}
1242
1243//____________________________________________________________________
1244/**
1245 * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
1246 * = 900GeV@f$
1247 * Work in progress
1248 *
1249 * @return graph of data
1250 *
1251 * @ingroup pwglf_forward_otherdata
1252 */
1253TGraphAsymmErrors* AliceCentralInel900Work()
1254{
1255 TGraphAsymmErrors *g = new TGraphAsymmErrors(18);
1256 // g->SetPoint(0,-1.9,0); g->SetPointError(0,0.1,0.1,0,0);
1257 g->SetPoint(0,-1.7,3.13935); g->SetPointError(0,0.1,0.1,0.0726186,0.0525276);
1258 g->SetPoint(1,-1.5,3.15634); g->SetPointError(1,0.1,0.1,0.0338547,0.0380273);
1259 g->SetPoint(2,-1.3,3.13683); g->SetPointError(2,0.1,0.1,0.0295176,0.0295638);
1260 g->SetPoint(3,-1.1,3.10618); g->SetPointError(3,0.1,0.1,0.0306925,0.0329387);
1261 g->SetPoint(4,-0.9,3.05921); g->SetPointError(4,0.1,0.1,0.0224684,0.025408);
1262 g->SetPoint(5,-0.7,3.00303); g->SetPointError(5,0.1,0.1,0.0389278,0.0238328);
1263 g->SetPoint(6,-0.5,2.94604); g->SetPointError(6,0.1,0.1,0.0211986,0.0322219);
1264 g->SetPoint(7,-0.3,2.91507); g->SetPointError(7,0.1,0.1,0.030029,0.0209573);
1265 g->SetPoint(8,-0.1,2.88965); g->SetPointError(8,0.1,0.1,0.0286516,0.0253694);
1266 g->SetPoint(9,0.1,2.89731); g->SetPointError(9,0.1,0.1,0.0334615,0.0192116);
1267 g->SetPoint(10,0.3,2.91188); g->SetPointError(10,0.1,0.1,0.0503868,0.0249106);
1268 g->SetPoint(11,0.5,2.96295); g->SetPointError(11,0.1,0.1,0.030009,0.0284692);
1269 g->SetPoint(12,0.7,3.0089); g->SetPointError(12,0.1,0.1,0.0189095,0.026319);
1270 g->SetPoint(13,0.9,3.07028); g->SetPointError(13,0.1,0.1,0.0449128,0.0307375);
1271 g->SetPoint(14,1.1,3.10215); g->SetPointError(14,0.1,0.1,0.0288688,0.0263013);
1272 g->SetPoint(15,1.3,3.12946); g->SetPointError(15,0.1,0.1,0.0431495,0.026355);
1273 g->SetPoint(16,1.5,3.14549); g->SetPointError(16,0.1,0.1,0.0322482,0.0336111);
1274 g->SetPoint(17,1.7,3.15729); g->SetPointError(17,0.1,0.1,0.105509,0.0523796);
1275 // g->SetPoint(19,1.9,0); g->SetPointError(19,0.1,0.1,0,0);
1276
1277 SetGraphAttributes(g, INEL, WIP, false, "alice_pp900work",
1278 "ALICE INEL - work in progress");
1279 return g;
1280}
1281
1282//____________________________________________________________________
1283/**
1284 * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
1285 * = 900GeV@f$
1286 * Work in progress
1287 *
1288 * @return graph of data
1289 *
1290 * @ingroup pwglf_forward_otherdata
1291 */
1292TGraphAsymmErrors* AliceCentralNsd900Work()
1293{
1294 TGraphAsymmErrors *g = new TGraphAsymmErrors(18);
1295
1296 g->SetPoint(0,-1.7,3.84726); g->SetPointError(0,0.1,0.1,0.114853,0.118974);
1297 g->SetPoint(1,-1.5,3.87094); g->SetPointError(1,0.1,0.1,0.10574,0.108613);
1298 g->SetPoint(2,-1.3,3.84769); g->SetPointError(2,0.1,0.1,0.105942,0.107644);
1299 g->SetPoint(3,-1.1,3.8122); g->SetPointError(3,0.1,0.1,0.100838,0.101818);
1300 g->SetPoint(4,-0.9,3.75388); g->SetPointError(4,0.1,0.1,0.0967073,0.0972099);
1301 g->SetPoint(5,-0.7,3.68733); g->SetPointError(5,0.1,0.1,0.0923424,0.0925662);
1302 g->SetPoint(6,-0.5,3.61874); g->SetPointError(6,0.1,0.1,0.0904027,0.090477);
1303 g->SetPoint(7,-0.3,3.58091); g->SetPointError(7,0.1,0.1,0.0875915,0.0875992);
1304 g->SetPoint(8,-0.1,3.54905); g->SetPointError(8,0.1,0.1,0.086046,0.0860293);
1305 g->SetPoint(9,0.1,3.55968); g->SetPointError(9,0.1,0.1,0.0884776,0.0884612);
1306 g->SetPoint(10,0.3,3.57729); g->SetPointError(10,0.1,0.1,0.0857614,0.0857693);
1307 g->SetPoint(11,0.5,3.63879); g->SetPointError(11,0.1,0.1,0.0879787,0.0880559);
1308 g->SetPoint(12,0.7,3.69422); g->SetPointError(12,0.1,0.1,0.0931736,0.0933962);
1309 g->SetPoint(13,0.9,3.76835); g->SetPointError(13,0.1,0.1,0.0928833,0.0934106);
1310 g->SetPoint(14,1.1,3.80647); g->SetPointError(14,0.1,0.1,0.0998124,0.100799);
1311 g->SetPoint(15,1.3,3.83824); g->SetPointError(15,0.1,0.1,0.102549,0.104298);
1312 g->SetPoint(16,1.5,3.85778); g->SetPointError(16,0.1,0.1,0.10267,0.105607);
1313 g->SetPoint(17,1.7,3.87075); g->SetPointError(17,0.1,0.1,0.114093,0.11829);
1314
1315 SetGraphAttributes(g, NSD, WIP, false, "alice_pp900nsdwork",
1316 "ALICE NSD - work in progress");
1317 return g;
1318}
1319
1320//____________________________________________________________________
1321/**
1322 * Get the ALICE INEL>0 data in @f$ |\eta|<1.8@f$ for pp at @f$ \sqrt{s}
1323 * = 900GeV@f$
1324 * Work in progress
1325 *
1326 * @return graph of data
1327 *
1328 * @ingroup pwglf_forward_otherdata
1329 */
1330TGraphAsymmErrors* AliceCentralInelGt900Work()
1331{
1332 TGraphAsymmErrors *g = new TGraphAsymmErrors(10);
1333
1334 g->SetPoint(0,-0.9,3.90755); g->SetPointError(0,0.1,0.1,0.0377085,0.0355101);
1335 g->SetPoint(1,-0.7,3.8357); g->SetPointError(1,0.1,0.1,0.0477674,0.0349303);
1336 g->SetPoint(2,-0.5,3.76291); g->SetPointError(2,0.1,0.1,0.0277709,0.040401);
1337 g->SetPoint(3,-0.3,3.72336); g->SetPointError(3,0.1,0.1,0.0343553,0.0250805);
1338 g->SetPoint(4,-0.1,3.69098); g->SetPointError(4,0.1,0.1,0.0324842,0.0324248);
1339 g->SetPoint(5,0.1,3.70076); g->SetPointError(5,0.1,0.1,0.0390932,0.0246738);
1340 g->SetPoint(6,0.3,3.71924); g->SetPointError(6,0.1,0.1,0.0576054,0.0287106);
1341 g->SetPoint(7,0.5,3.7844); g->SetPointError(7,0.1,0.1,0.0316759,0.0295124);
1342 g->SetPoint(8,0.7,3.84319); g->SetPointError(8,0.1,0.1,0.0293134,0.0332125);
1343 g->SetPoint(9,0.9,3.92163); g->SetPointError(9,0.1,0.1,0.0558339,0.0394925);
1344
1345 SetGraphAttributes(g, INELGt0, WIP, false, "alice_pp900inelgtwork",
1346 "ALICE INEL>0 - work in progress");
1347 return g;
1348}
1349
1350//____________________________________________________________________
1351/**
1352 * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
1353 * at @f$ \sqrt{s} = 2760GeV@f$
1354 * Work in progress
1355 *
1356 * @return graph of data
1357 *
1358 * @ingroup pwglf_forward_otherdata
1359 */
1360TGraphAsymmErrors* AliceCentralInel2760Work()
1361{
1362 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1363 g->SetPoint(0,-1.7,4.0313); g->SetPointError(0,0.1,0.1,0.100951,0.104678);
1364 g->SetPoint(1,-1.5,4.0431); g->SetPointError(1,0.1,0.1,0.10129,0.103883);
1365 g->SetPoint(2,-1.3,4.01251); g->SetPointError(2,0.1,0.1,0.10847,0.110089);
1366 g->SetPoint(3,-1.1,3.96799); g->SetPointError(3,0.1,0.1,0.105543,0.106623);
1367 g->SetPoint(4,-0.9,3.89669); g->SetPointError(4,0.1,0.1,0.110974,0.111625);
1368 g->SetPoint(5,-0.7,3.81051); g->SetPointError(5,0.1,0.1,0.108463,0.108882);
1369 g->SetPoint(6,-0.5,3.76537); g->SetPointError(6,0.1,0.1,0.105488,0.105773);
1370 g->SetPoint(7,-0.3,3.69733); g->SetPointError(7,0.1,0.1,0.110156,0.11035);
1371 g->SetPoint(8,-0.1,3.68148); g->SetPointError(8,0.1,0.1,0.105564,0.105733);
1372 g->SetPoint(9,0.1,3.67386); g->SetPointError(9,0.1,0.1,0.1058,0.105968);
1373 g->SetPoint(10,0.3,3.69873); g->SetPointError(10,0.1,0.1,0.107167,0.107367);
1374 g->SetPoint(11,0.5,3.76377); g->SetPointError(11,0.1,0.1,0.111177,0.111448);
1375 g->SetPoint(12,0.7,3.81956); g->SetPointError(12,0.1,0.1,0.107198,0.107623);
1376 g->SetPoint(13,0.9,3.89506); g->SetPointError(13,0.1,0.1,0.105617,0.1063);
1377 g->SetPoint(14,1.1,3.95888); g->SetPointError(14,0.1,0.1,0.111316,0.112336);
1378 g->SetPoint(15,1.3,4.00176); g->SetPointError(15,0.1,0.1,0.111751,0.113315);
1379 g->SetPoint(16,1.5,4.03247); g->SetPointError(16,0.1,0.1,0.114383,0.116674);
1380 g->SetPoint(17,1.7,4.061); g->SetPointError(17,0.1,0.1,0.107094,0.110665);
1381
1382 SetGraphAttributes(g, INEL, WIP, false,
1383 "alice_ppInel2760Work",
1384 "ALICE INEL - work in progress");
1385 return g;
1386}
1387//____________________________________________________________________
1388/**
1389 * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
1390 * at @f$ \sqrt{s} = 2760GeV@f$
1391 * Work in progress
1392 *
1393 * @return graph of data
1394 *
1395 * @ingroup pwglf_forward_otherdata
1396 */
1397TGraphAsymmErrors* AliceCentralNsd2760Work()
1398{
1399 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1400 g->SetPoint(0,-1.7,4.8704); g->SetPointError(0,0.1,0.1,0.221293,0.224755);
1401 g->SetPoint(1,-1.5,4.88859); g->SetPointError(1,0.1,0.1,0.221269,0.223478);
1402 g->SetPoint(2,-1.3,4.85326); g->SetPointError(2,0.1,0.1,0.21109,0.212455);
1403 g->SetPoint(3,-1.1,4.80085); g->SetPointError(3,0.1,0.1,0.21041,0.211157);
1404 g->SetPoint(4,-0.9,4.71513); g->SetPointError(4,0.1,0.1,0.198361,0.198749);
1405 g->SetPoint(5,-0.7,4.61153); g->SetPointError(5,0.1,0.1,0.194009,0.194176);
1406 g->SetPoint(6,-0.5,4.55715); g->SetPointError(6,0.1,0.1,0.193226,0.193281);
1407 g->SetPoint(7,-0.3,4.47508); g->SetPointError(7,0.1,0.1,0.182433,0.182439);
1408 g->SetPoint(8,-0.1,4.45709); g->SetPointError(8,0.1,0.1,0.186518,0.186506);
1409 g->SetPoint(9,0.1,4.44707); g->SetPointError(9,0.1,0.1,0.185747,0.185735);
1410 g->SetPoint(10,0.3,4.47734); g->SetPointError(10,0.1,0.1,0.185835,0.185841);
1411 g->SetPoint(11,0.5,4.55477); g->SetPointError(11,0.1,0.1,0.186934,0.186991);
1412 g->SetPoint(12,0.7,4.62236); g->SetPointError(12,0.1,0.1,0.196631,0.196796);
1413 g->SetPoint(13,0.9,4.71277); g->SetPointError(13,0.1,0.1,0.204034,0.20441);
1414 g->SetPoint(14,1.1,4.78902); g->SetPointError(14,0.1,0.1,0.20317,0.20394);
1415 g->SetPoint(15,1.3,4.84008); g->SetPointError(15,0.1,0.1,0.205573,0.206967);
1416 g->SetPoint(16,1.5,4.87453); g->SetPointError(16,0.1,0.1,0.206314,0.208667);
1417 g->SetPoint(17,1.7,4.90614); g->SetPointError(17,0.1,0.1,0.218996,0.222545);
1418
1419 SetGraphAttributes(g, NSD, WIP, false,
1420 "alice_ppNsd2760Work",
1421 "ALICE NSD - work in progress");
1422 return g;
1423}
1424//____________________________________________________________________
1425/**
1426 * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
1427 * at @f$ \sqrt{s} = 2760GeV@f$
1428 * Work in progress
1429 *
1430 * @return graph of data
1431 *
1432 * @ingroup pwglf_forward_otherdata
1433 */
1434TGraphAsymmErrors* AliceCentralInelGt2760Work()
1435{
1436 TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
1437 g->SetPoint(0,-0.9,4.96315); g->SetPointError(0,0.1,0.1,0.0439746,0.0440108);
1438 g->SetPoint(1,-0.7,4.8532); g->SetPointError(1,0.1,0.1,0.0426373,0.0600727);
1439 g->SetPoint(2,-0.5,4.79582); g->SetPointError(2,0.1,0.1,0.0475367,0.0466255);
1440 g->SetPoint(3,-0.3,4.70907); g->SetPointError(3,0.1,0.1,0.0313084,0.0468084);
1441 g->SetPoint(4,-0.1,4.68906); g->SetPointError(4,0.1,0.1,0.0413149,0.0397909);
1442 g->SetPoint(5,0.1,4.67937); g->SetPointError(5,0.1,0.1,0.0346151,0.0450248);
1443 g->SetPoint(6,0.3,4.7109); g->SetPointError(6,0.1,0.1,0.0408403,0.0839992);
1444 g->SetPoint(7,0.5,4.79359); g->SetPointError(7,0.1,0.1,0.0324516,0.0357053);
1445 g->SetPoint(8,0.7,4.86469); g->SetPointError(8,0.1,0.1,0.0452175,0.0477304);
1446 g->SetPoint(9,0.9,4.96078); g->SetPointError(9,0.1,0.1,0.0566798,0.0804077);
1447
1448 SetGraphAttributes(g, INELGt0, WIP, false,
1449 "alice_ppInelGt2760Work",
1450 "ALICE INELGt0 - work in progress");
1451 return g;
1452}
1453
1454//____________________________________________________________________
1455/**
1456 * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
1457 * at @f$ \sqrt{s} = 7000GeV@f$
1458 * Work in progress
1459 *
1460 * @return graph of data
1461 *
1462 * @ingroup pwglf_forward_otherdata
1463 */
1464TGraphAsymmErrors* AliceCentralInel7000Work()
1465{
1466 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1467 g->SetPoint(0,-1.7,4.97541); g->SetPointError(0,0.1,0.1,0.187526,0.162049);
1468 g->SetPoint(1,-1.5,4.98161); g->SetPointError(1,0.1,0.1,0.128353,0.149085);
1469 g->SetPoint(2,-1.3,4.94853); g->SetPointError(2,0.1,0.1,0.129841,0.144762);
1470 g->SetPoint(3,-1.1,4.88924); g->SetPointError(3,0.1,0.1,0.137866,0.157862);
1471 g->SetPoint(4,-0.9,4.79998); g->SetPointError(4,0.1,0.1,0.144492,0.158783);
1472 g->SetPoint(5,-0.7,4.71399); g->SetPointError(5,0.1,0.1,0.132703,0.156135);
1473 g->SetPoint(6,-0.5,4.63098); g->SetPointError(6,0.1,0.1,0.129938,0.147085);
1474 g->SetPoint(7,-0.3,4.56815); g->SetPointError(7,0.1,0.1,0.129424,0.145485);
1475 g->SetPoint(8,-0.1,4.52372); g->SetPointError(8,0.1,0.1,0.129049,0.145285);
1476 g->SetPoint(9,0.1,4.52946); g->SetPointError(9,0.1,0.1,0.131266,0.144285);
1477 g->SetPoint(10,0.3,4.56411); g->SetPointError(10,0.1,0.1,0.130652,0.149019);
1478 g->SetPoint(11,0.5,4.63554); g->SetPointError(11,0.1,0.1,0.133415,0.144298);
1479 g->SetPoint(12,0.7,4.71592); g->SetPointError(12,0.1,0.1,0.136436,0.151768);
1480 g->SetPoint(13,0.9,4.8059); g->SetPointError(13,0.1,0.1,0.136996,0.142551);
1481 g->SetPoint(14,1.1,4.88457); g->SetPointError(14,0.1,0.1,0.134237,0.142764);
1482 g->SetPoint(15,1.3,4.92903); g->SetPointError(15,0.1,0.1,0.131933,0.152767);
1483 g->SetPoint(16,1.5,4.96487); g->SetPointError(16,0.1,0.1,0.140214,0.147354);
1484 g->SetPoint(17,1.7,4.95502); g->SetPointError(17,0.1,0.1,0.156906,0.14759);
1485
1486 SetGraphAttributes(g, INEL, WIP, false,
1487 "alice_ppInel7000Work",
1488 "ALICE INEL - work in progress");
1489 return g;
1490}
1491//____________________________________________________________________
1492/**
1493 * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
1494 * at @f$ \sqrt{s} = 7000GeV@f$
1495 * Work in progress
1496 *
1497 * @return graph of data
1498 *
1499 * @ingroup pwglf_forward_otherdata
1500 */
1501TGraphAsymmErrors* AliceCentralNsd7000Work()
1502{
1503 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1504 g->SetPoint(0,-1.7,6.12747); g->SetPointError(0,0.1,0.1,0.157334,0.164918);
1505 g->SetPoint(1,-1.5,6.1353); g->SetPointError(1,0.1,0.1,0.146834,0.152011);
1506 g->SetPoint(2,-1.3,6.09648); g->SetPointError(2,0.1,0.1,0.140067,0.143287);
1507 g->SetPoint(3,-1.1,6.02552); g->SetPointError(3,0.1,0.1,0.133435,0.135281);
1508 g->SetPoint(4,-0.9,5.91705); g->SetPointError(4,0.1,0.1,0.129449,0.130381);
1509 g->SetPoint(5,-0.7,5.81246); g->SetPointError(5,0.1,0.1,0.126477,0.126883);
1510 g->SetPoint(6,-0.5,5.71104); g->SetPointError(6,0.1,0.1,0.124521,0.124655);
1511 g->SetPoint(7,-0.3,5.63422); g->SetPointError(7,0.1,0.1,0.120116,0.12013);
1512 g->SetPoint(8,-0.1,5.57977); g->SetPointError(8,0.1,0.1,0.119286,0.119256);
1513 g->SetPoint(9,0.1,5.58662); g->SetPointError(9,0.1,0.1,0.119331,0.119301);
1514 g->SetPoint(10,0.3,5.6291); g->SetPointError(10,0.1,0.1,0.120683,0.120697);
1515 g->SetPoint(11,0.5,5.7166); g->SetPointError(11,0.1,0.1,0.122787,0.122923);
1516 g->SetPoint(12,0.7,5.81463); g->SetPointError(12,0.1,0.1,0.126293,0.1267);
1517 g->SetPoint(13,0.9,5.92404); g->SetPointError(13,0.1,0.1,0.129522,0.130456);
1518 g->SetPoint(14,1.1,6.01958); g->SetPointError(14,0.1,0.1,0.134505,0.136333);
1519 g->SetPoint(15,1.3,6.07232); g->SetPointError(15,0.1,0.1,0.140728,0.143909);
1520 g->SetPoint(16,1.5,6.11596); g->SetPointError(16,0.1,0.1,0.14756,0.15268);
1521 g->SetPoint(17,1.7,6.10155); g->SetPointError(17,0.1,0.1,0.151731,0.159518);
1522
1523 SetGraphAttributes(g, NSD, WIP, false,
1524 "alice_ppNsd7000Work",
1525 "ALICE NSD - work in progress");
1526 return g;
1527}
1528
1529//____________________________________________________________________
1530/**
1531 * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
1532 * at @f$ \sqrt{s} = 7000GeV@f$
1533 * Work in progress
1534 *
1535 * @return graph of data
1536 *
1537 * @ingroup pwglf_forward_otherdata
1538 */
1539TGraphAsymmErrors* AliceCentralInelGt7000Work()
1540{
1541 TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
1542 g->SetPoint(0,-0.9,6.22689); g->SetPointError(0,0.1,0.1,0.0959094,0.10395);
1543 g->SetPoint(1,-0.7,6.11603); g->SetPointError(1,0.1,0.1,0.0609242,0.0984269);
1544 g->SetPoint(2,-0.5,6.00881); g->SetPointError(2,0.1,0.1,0.0595691,0.0842045);
1545 g->SetPoint(3,-0.3,5.9274); g->SetPointError(3,0.1,0.1,0.0560837,0.0780806);
1546 g->SetPoint(4,-0.1,5.86988); g->SetPointError(4,0.1,0.1,0.0552611,0.0798584);
1547 g->SetPoint(5,0.1,5.8773); g->SetPointError(5,0.1,0.1,0.062512,0.077947);
1548 g->SetPoint(6,0.3,5.92215); g->SetPointError(6,0.1,0.1,0.0535152,0.0863595);
1549 g->SetPoint(7,0.5,6.01458); g->SetPointError(7,0.1,0.1,0.0578218,0.0745799);
1550 g->SetPoint(8,0.7,6.1186); g->SetPointError(8,0.1,0.1,0.0767397,0.0899574);
1551 g->SetPoint(9,0.9,6.23468); g->SetPointError(9,0.1,0.1,0.0786932,0.073295);
1552
1553 SetGraphAttributes(g, INELGt0, WIP, false,
1554 "alice_ppInelGt7000Work",
1555 "ALICE INELGt0 - work in progress");
1556 return g;
1557}
1558//____________________________________________________________________
1559/**
1560 * Get the ALICE INEL data in @f$ |\eta|<1.8@f$ for pp
1561 * at @f$ \sqrt{s} = 8000GeV@f$
1562 * Work in progress
1563 *
1564 * @return graph of data
1565 *
1566 * @ingroup pwglf_forward_otherdata
1567 */
1568TGraphAsymmErrors* AliceCentralInel8000Work()
1569{
1570 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1571 g->SetPoint(0,-1.7,5.1333); g->SetPointError(0,0.1,0.1,0.0613865,0.0707879);
1572 g->SetPoint(1,-1.5,5.14341); g->SetPointError(1,0.1,0.1,0.117829,0.121428);
1573 g->SetPoint(2,-1.3,5.13589); g->SetPointError(2,0.1,0.1,0.110807,0.113393);
1574 g->SetPoint(3,-1.1,5.06167); g->SetPointError(3,0.1,0.1,0.120091,0.121635);
1575 g->SetPoint(4,-0.9,4.97796); g->SetPointError(4,0.1,0.1,0.116141,0.117154);
1576 g->SetPoint(5,-0.7,4.88431); g->SetPointError(5,0.1,0.1,0.11944,0.120064);
1577 g->SetPoint(6,-0.5,4.81236); g->SetPointError(6,0.1,0.1,0.11049,0.110936);
1578 g->SetPoint(7,-0.3,4.72239); g->SetPointError(7,0.1,0.1,0.110969,0.111284);
1579 g->SetPoint(8,-0.1,4.66962); g->SetPointError(8,0.1,0.1,0.125108,0.125337);
1580 g->SetPoint(9,0.1,4.69441); g->SetPointError(9,0.1,0.1,0.113766,0.114021);
1581 g->SetPoint(10,0.3,4.7335); g->SetPointError(10,0.1,0.1,0.104531,0.104866);
1582 g->SetPoint(11,0.5,4.79917); g->SetPointError(11,0.1,0.1,0.107076,0.107534);
1583 g->SetPoint(12,0.7,4.88713); g->SetPointError(12,0.1,0.1,0.106124,0.106827);
1584 g->SetPoint(13,0.9,4.98035); g->SetPointError(13,0.1,0.1,0.120107,0.121087);
1585 g->SetPoint(14,1.1,5.05366); g->SetPointError(14,0.1,0.1,0.115795,0.11739);
1586 g->SetPoint(15,1.3,5.11276); g->SetPointError(15,0.1,0.1,0.123574,0.125877);
1587 g->SetPoint(16,1.5,5.16105); g->SetPointError(16,0.1,0.1,0.0979751,0.102305);
1588 g->SetPoint(17,1.7,5.16477); g->SetPointError(17,0.1,0.1,0.116096,0.121392);
1589
1590 SetGraphAttributes(g, INEL, WIP, false,
1591 "alice_ppInel8000Work",
1592 "ALICE INEL - work in progress");
1593 return g;
1594}
1595
1596//____________________________________________________________________
1597/**
1598 * Get the ALICE NSD data in @f$ |\eta|<1.8@f$ for pp
1599 * at @f$ \sqrt{s} = 8000GeV@f$
1600 * Work in progress
1601 *
1602 * @return graph of data
1603 *
1604 * @ingroup pwglf_forward_otherdata
1605 */
1606TGraphAsymmErrors* AliceCentralNsd8000Work()
1607{
1608 TGraphAsymmErrors* g = new TGraphAsymmErrors(18);
1609 g->SetPoint(0,-1.7,6.28722); g->SetPointError(0,0.1,0.1,0.220541,0.2263);
1610 g->SetPoint(1,-1.5,6.29211); g->SetPointError(1,0.1,0.1,0.126653,0.132922);
1611 g->SetPoint(2,-1.3,6.29313); g->SetPointError(2,0.1,0.1,0.142432,0.145805);
1612 g->SetPoint(3,-1.1,6.1944); g->SetPointError(3,0.1,0.1,0.116871,0.119093);
1613 g->SetPoint(4,-0.9,6.09529); g->SetPointError(4,0.1,0.1,0.121064,0.122121);
1614 g->SetPoint(5,-0.7,5.97811); g->SetPointError(5,0.1,0.1,0.111525,0.112012);
1615 g->SetPoint(6,-0.5,5.88992); g->SetPointError(6,0.1,0.1,0.119488,0.119637);
1616 g->SetPoint(7,-0.3,5.78296); g->SetPointError(7,0.1,0.1,0.114947,0.114962);
1617 g->SetPoint(8,-0.1,5.71633); g->SetPointError(8,0.1,0.1,0.0933,0.09326);
1618 g->SetPoint(9,0.1,5.74663); g->SetPointError(9,0.1,0.1,0.109892,0.109857);
1619 g->SetPoint(10,0.3,5.79472); g->SetPointError(10,0.1,0.1,0.123704,0.123718);
1620 g->SetPoint(11,0.5,5.87545); g->SetPointError(11,0.1,0.1,0.122522,0.122667);
1621 g->SetPoint(12,0.7,5.98273); g->SetPointError(12,0.1,0.1,0.128316,0.128739);
1622 g->SetPoint(13,0.9,6.09037); g->SetPointError(13,0.1,0.1,0.114321,0.115437);
1623 g->SetPoint(14,1.1,6.18105); g->SetPointError(14,0.1,0.1,0.125412,0.127476);
1624 g->SetPoint(15,1.3,6.24275); g->SetPointError(15,0.1,0.1,0.118631,0.122597);
1625 g->SetPoint(16,1.5,6.28916); g->SetPointError(16,0.1,0.1,0.144205,0.149736);
1626 g->SetPoint(17,1.7,6.28878); g->SetPointError(17,0.1,0.1,0.134438,0.143695);
1627
1628 SetGraphAttributes(g, NSD, WIP, false,
1629 "alice_ppNsd8000Work",
1630 "ALICE NSD - work in progress");
1631 return g;
1632}
1633//____________________________________________________________________
1634/**
1635 * Get the ALICE INELGt0 data in @f$ |\eta|<1.8@f$ for pp
1636 * at @f$ \sqrt{s} = 8000GeV@f$
1637 * Work in progress
1638 *
1639 * @return graph of data
1640 *
1641 * @ingroup pwglf_forward_otherdata
1642 */
1643TGraphAsymmErrors* AliceCentralInelGt8000Work()
1644{
1645 TGraphAsymmErrors* g = new TGraphAsymmErrors(10);
1646 g->SetPoint(0,-0.9,6.38567); g->SetPointError(0,0.1,0.1,0.0436571,0.0436571);
1647 g->SetPoint(1,-0.7,6.26363); g->SetPointError(1,0.1,0.1,0.0312036,0.0312036);
1648 g->SetPoint(2,-0.5,6.17205); g->SetPointError(2,0.1,0.1,0.0351509,0.0351509);
1649 g->SetPoint(3,-0.3,6.05629); g->SetPointError(3,0.1,0.1,0.0302028,0.0302028);
1650 g->SetPoint(4,-0.1,5.98823); g->SetPointError(4,0.1,0.1,0.0141541,0.0141541);
1651 g->SetPoint(5,0.1,6.02043); g->SetPointError(5,0.1,0.1,0.0256893,0.0256893);
1652 g->SetPoint(6,0.3,6.07111); g->SetPointError(6,0.1,0.1,0.0380304,0.0380304);
1653 g->SetPoint(7,0.5,6.15492); g->SetPointError(7,0.1,0.1,0.0384435,0.0384435);
1654 g->SetPoint(8,0.7,6.26781); g->SetPointError(8,0.1,0.1,0.0450579,0.0450579);
1655 g->SetPoint(9,0.9,6.38491); g->SetPointError(9,0.1,0.1,0.0396431,0.0396431);
1656
1657 SetGraphAttributes(g, INELGt0, WIP, false,
1658 "alice_ppInelGt8000Work",
1659 "ALICE INELGt0 - work in progress");
6d233add 1660 return g;
1661}
7e4038b5 1662
1663//____________________________________________________________________
1664/**
1665 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
1666 * \sqrt{s} = 900GeV@f$
1667 *
1668 * p7743_d8x1y1
1669 *
1670 * @return graph of data
1671 *
bd6f5206 1672 * @ingroup pwglf_forward_otherdata
7e4038b5 1673 */
1674TGraphAsymmErrors* CMSNsd900()
1675{
1676 // CMS published NSD data - p7743_d8x1y1
1677 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75, 1.25, 1.75,
1678 2.25 };
1679 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1680 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1681 double y[] = { 3.6, 3.73, 3.62, 3.54, 3.48, 3.48, 3.54, 3.62, 3.73, 3.6 };
1682 double eym[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14,0.13 };
1683 double eyp[] = { 0.13, 0.14, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.14, 0.13 };
1684 const int np = 10;
1685 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 1686 SetGraphAttributes(g, NSD, CMS, false, "cms_nsd900", "CMS NSD");
7e4038b5 1687
1688 return g;
1689}
1690
1691
1692//____________________________________________________________________
1693/**
1694 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
1695 * \sqrt{s} = 2.36GeV@f$
1696 *
1697 * p7743_d8x1y2
1698 *
1699 * @return graph of data
1700 *
bd6f5206 1701 * @ingroup pwglf_forward_otherdata
7e4038b5 1702 */
1703TGraphAsymmErrors* CMSNsd2360()
1704{
1705 // CMS NSD 2360 - p7743_d8x1y2
1706 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
1707 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1708 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1709 double y[] = { 4.78, 4.81, 4.66, 4.61, 4.47, 4.47, 4.61, 4.66, 4.81, 4.78 };
1710 double eym[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
1711 double eyp[] = { 0.17, 0.18, 0.17, 0.17, 0.16, 0.16, 0.17, 0.17, 0.18, 0.17 };
1712 const int np = 10;
1713 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 1714 SetGraphAttributes(g, NSD, CMS, false, "cms_nsd2360", "CMS NSD");
7e4038b5 1715 return g;
1716}
1717
1718
1719//____________________________________________________________________
1720/**
1721 * Get the CMS NSD data in @f$ |\eta|<2.25@f$ for pp at @f$
1722 * \sqrt{s} = 7TeV@f$
1723 *
1724 * p7838_d5x1y1
1725 *
1726 * @return graph of data
1727 *
bd6f5206 1728 * @ingroup pwglf_forward_otherdata
7e4038b5 1729 */
1730TGraphAsymmErrors* CMSNsd7000()
1731{
1732 // CMS NSD 7000 - Plot: p7838_d5x1y1
1733 double x[] = { -2.25, -1.75, -1.25, -0.75, -0.25, 0.25, 0.75,1.25,1.75,2.25 };
1734 double exm[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1735 double exp[] = { 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25 };
1736 double y[] = { 6.18, 6.26, 6.14, 6.01, 5.78, 5.78, 6.01, 6.14, 6.26, 6.18 };
1737 double eym[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
1738 double eyp[] = { 0.25, 0.25, 0.24, 0.24, 0.23, 0.23, 0.24, 0.24, 0.25, 0.25 };
1739 const int np = 10;
1740 TGraphAsymmErrors* g = new TGraphAsymmErrors(np, x, y, exm, exp, eym, eyp);
ffca499d 1741 SetGraphAttributes(g, NSD, CMS, false, "cms_nsd7000", "CMS NSD");
7e4038b5 1742 return g;
1743}
1744
797161e8 1745//____________________________________________________________________
290052e7 1746/**
1747 * Get a single data graph
1748 *
1749 * @param which Which type
1750 * @param sys Collisition system
1751 * @param energy Collision energy
1752 * @param type Trigger type
1753 * @param centLow Low cut on centrality
1754 * @param centHigh Up cut on centraltiy
1755 *
1756 * @return Data graph
1757 */
797161e8 1758TGraphAsymmErrors*
1759GetSingle(UShort_t which,
1760 UShort_t sys,
1761 UShort_t energy,
1762 UShort_t type=0x1,
05cc5a69 1763 UShort_t centLow=0,
797161e8 1764 UShort_t centHigh=0)
1765{
1766 TGraphAsymmErrors* ret = 0;
1767 if (sys == 1) {
1768 if (TMath::Abs(energy-900) < 10) {
1769 switch (type) {
1770 case 1: // INEL
1771 switch (which) {
bfab35d9 1772 case PYTHIA: ret = Pythia900INEL(); break;
1773 case UA5: ret = UA5Inel(false); break;
1774 case UA5+10: ret = UA5Inel(true); break;
1775 case ALICE: ret = AliceCentralInel900(); break;
1776 case WIP: ret = AliceCentralInel900Work(); break;
797161e8 1777 }
1778 break;
1779 case 2: // INEL>0
1780 switch (which) {
1781 case ALICE: ret = AliceCentralInelGt900(); break;
bfab35d9 1782 case WIP: ret = AliceCentralInelGt900Work(); break;
797161e8 1783 }
1784 break;
1785 case 4: // NSD
1786 switch (which) {
1787 case PYTHIA: ret = Pythia900NSD(); break;
1788 case UA5: ret = UA5Nsd(false); break;
1789 case UA5+10: ret = UA5Nsd(true); break;
1790 case ALICE: ret = AliceCentralNsd900(); break;
bfab35d9 1791 case WIP: ret = AliceCentralNsd900Work(); break;
797161e8 1792 case CMS: ret = CMSNsd900(); break;
1793 }
1794 break;
1795 } // type
1796 }
1797 else if (TMath::Abs(energy-2360) < 10) {
1798 switch (type) {
1799 case 1: // INEL
1800 switch (which) {
1801 case ALICE: ret = AliceCentralInel2360(); break;
bfab35d9 1802 case WIP: ret = AliceCentralInel2760Work(); break;
797161e8 1803 }
1804 break;
1805 case 2: // INEL > 0
1806 switch (which) {
1807 case ALICE: ret = AliceCentralInelGt2360(); break;
bfab35d9 1808 case WIP: ret = AliceCentralInelGt2760Work(); break;
797161e8 1809 }
1810 break;
1811 case 4: // NSD
1812 switch (which) {
1813 case ALICE: ret = AliceCentralNsd2360(); break;
1814 case CMS: ret = CMSNsd2360(); break;
bfab35d9 1815 case WIP: ret = AliceCentralNsd2760Work(); break;
1816 }
1817 break;
1818 }
1819 }
1820 else if (TMath::Abs(energy-2760) < 10) {
1821 switch (type) {
1822 case 1: // INEL
1823 switch (which) {
1824 case WIP: ret = AliceCentralInel2760Work(); break;
1825 }
1826 break;
1827 case 2: // INEL > 0
1828 switch (which) {
1829 case WIP: ret = AliceCentralInelGt2760Work(); break;
1830 }
1831 break;
1832 case 4: // NSD
1833 switch (which) {
1834 case WIP: ret = AliceCentralNsd2760Work(); break;
797161e8 1835 }
1836 break;
1837 }
1838 }
1839 else if (TMath::Abs(energy-7000) < 10) {
1840 switch (type) {
bfab35d9 1841 case 1:
1842 switch (which) {
1843 case WIP: ret = AliceCentralInel7000Work(); break;
1844 }
1845 break;
797161e8 1846 case 2: // INEL > 0
1847 switch (which) {
1848 case ALICE: ret = AliceCentralInelGt7000(); break;
bfab35d9 1849 case WIP: ret = AliceCentralInelGt7000Work(); break;
797161e8 1850 }
1851 break;
1852 case 4: // NSD
1853 switch (which) {
1854 case CMS: ret = CMSNsd7000(); break;
bfab35d9 1855 case WIP: ret = AliceCentralNsd7000Work(); break;
1856 }
1857 break;
1858 }
1859 }
1860 else if (TMath::Abs(energy-8000) < 10) {
1861 switch (type) {
1862 case 1:
1863 switch (which) {
1864 case WIP: ret = AliceCentralInel8000Work(); break;
1865 }
1866 break;
1867 case 2: // INEL > 0
1868 switch (which) {
1869 case WIP: ret = AliceCentralInelGt8000Work(); break;
1870 }
1871 break;
1872 case 4: // NSD
1873 switch (which) {
1874 case WIP: ret = AliceCentralNsd8000Work(); break;
797161e8 1875 }
1876 break;
1877 }
1878 }
1879 }
6d233add 1880 else if (sys == 3) { // pPb
05cc5a69 1881 if (TMath::Abs(energy - 5023) < 10 ||
1882 TMath::Abs(energy - 8000) < 10) {
6d233add 1883 switch (which) {
1884 case ALICE: ret = AliceCentralpPb5023(); break;
1885 }
1886 }
1887 }
05cc5a69 1888 if (!ret) {
bfab35d9 1889 TString w;
1890 switch (which) {
1891 case UA5: w = "UA5"; break;
1892 case CMS: w = "CMS"; break;
1893 case ALICE: w = "ALICE"; break;
1894 case WIP: w = "WIP"; break;
1895 case PYTHIA:w = "Pyhthia"; break;
1896 default: w = Form("unknown(%d)", which);
1897 }
1898 TString sy;
1899 switch (sys) {
1900 case 1: sy = "pp"; break;
1901 case 2: sy = "PbPb"; break;
1902 case 3: sy = "pPb"; break;
1903 default: sy = Form("unknown(%d)", sys);
1904 }
1905 TString tr;
1906 switch (type) {
1907 case 1: tr = "INEL"; break;
1908 case 2: tr = "INEL>0"; break;
1909 case 4: tr = "NSD"; break;
1910 default: tr = Form("unknown(%d)", sys);
1911 }
1912
1913 Warning("GetSingle", "Nothing to get for "
1914 "which=%s, sys=%s, energy=%dGeV, type=%s, centLow=%d, centHigh=%d",
1915 w.Data(), sy.Data(), energy, tr.Data(), centLow, centHigh);
05cc5a69 1916 }
f7cfc454 1917#if 0
1918 if (ret) {
1919 if (!wn.IsNull()) wn.Append(",");
1920 switch (which) {
1921 case PYTHIA: wn.Append("Pythia"); break;
1922 case UA5: wn.Append("UA5"); break;
1923 case CMS: wn.Append("CMS"); break;
1924 case ALICE: wn.Append("ALICE"); break;
1925 }
1926 }
1927#endif
797161e8 1928 return ret;
1929}
1930
f7cfc454 1931//____________________________________________________________________
290052e7 1932/**
1933 * Append an item to a list
1934 *
1935 * @param s List to append to
1936 * @param delim Delimiter
1937 * @param what What to append
1938 *
1939 * @return New string value
1940 */
f7cfc454 1941TString&
1942AppendItem(TString& s, char delim, const char* what)
1943{
1944 if (!s.IsNull()) s.Append(Form("%c", delim));
1945 s.Append(what);
1946 return s;
1947}
1948
7e4038b5 1949//____________________________________________________________________
1950/**
1951 * Get a multi graph of data for a given energy and trigger type
1952 *
e308a636 1953 * @param sys Collision system (1: pp, 2: PbPb)
7e4038b5 1954 * @param energy Energy in GeV (900, 2360, 7000)
1955 * @param type Bit pattern of trigger type
1956 * - 0x1 INEL
1957 * - 0x2 INEL>0
1958 * - 0x4 NSD
ffca499d 1959 * @param centLow Low centrality cut (only for PbPB)
1960 * @param centHigh High centrality cut (only for PbPB)
c6115ede 1961 * @param which What to get
7e4038b5 1962 *
1963 * @return A multi graph with the selected data.
1964 *
bd6f5206 1965 * @ingroup pwglf_forward_otherdata
7e4038b5 1966 */
1967TMultiGraph*
e308a636 1968GetData(UShort_t sys,
1969 UShort_t energy,
1970 UShort_t type=0x1,
1971 UShort_t centLow=0,
1972 UShort_t centHigh=0,
797161e8 1973 UShort_t which=0x7)
7e4038b5 1974{
e308a636 1975 TMultiGraph* mp = new TMultiGraph(Form("dndeta_%dGeV_%d_%03d_%03d",
1976 energy, type, centLow, centHigh),"");
7e4038b5 1977 TString tn;
1978 TString en;
e308a636 1979 TString sn;
1980 TString cn;
bfab35d9 1981 bool ua5 = (which & (1 << UA5)); // 0x1
1982 bool cms = (which & (1 << CMS)); // 0x2
1983 bool alice = (which & (1 << ALICE)); // 0x4
1984 bool work = (which & (1 << WIP)); // 0x8
1985 bool pythia = (which & (1 << PYTHIA)); // 0x10
797161e8 1986 en.Append(Form(", #sqrt{s%s}=", sys == 1 ? "" : "_{NN}"));
1987 if (energy < 1000)
1988 en.Append(Form("%dGeV", energy));
1989 else {
1990 if (energy % 1000 == 0)
1991 en.Append(Form("%dTeV", energy/1000));
e308a636 1992 else
797161e8 1993 en.Append(Form("%4.2fTeV", float(energy)/1000));
1994 }
1995
1996 if (sys == 1) {
e308a636 1997 if (!(type & 0x7))
1998 Warning("GetData", "Unknown trigger mask 0x%x", type);
1999
f7cfc454 2000 if (TMath::Abs(energy-2750) < 11) {
2001 Warning("GetData", "Using 2360GeV data for %dGeV comparison", energy);
2002 energy = 2360;
2003 }
797161e8 2004 if (!(TMath::Abs(energy-900) < 10 ||
2005 TMath::Abs(energy-2360) < 10 ||
bfab35d9 2006 TMath::Abs(energy-7000) < 10 ||
2007 TMath::Abs(energy-8000) < 10)) {
e308a636 2008 Warning("GetData", "No other results for sys=%d, energy=%d",
2009 sys, energy);
797161e8 2010 return 0;
2011 }
2012
2013 sn = "pp";
2014
bfab35d9 2015 // Substitute NSD for V0-AND
2016 if (type == 0x2000) type = 0x4;
f7cfc454 2017 if (type & 0x1) AppendItem(tn, '|', "INEL");
2018 if (type & 0x2) AppendItem(tn, '|', "INEL>0");
2019 if (type & 0x4) AppendItem(tn, '|', "NSD");
797161e8 2020
2021 Bool_t seenUA5 = false;
2022 for (Int_t i = 0; i < 3; i++) {
2023 UShort_t mask = (1 << i);
2024 if ((type & mask) == 0) continue;
f7cfc454 2025 TGraphAsymmErrors* gUAp =(ua5 ?GetSingle(UA5, sys,energy,mask):0);
2026 TGraphAsymmErrors* gUAn =(ua5 ?GetSingle(UA5+10,sys,energy,mask):0);
2027 TGraphAsymmErrors* gCMS =(cms ?GetSingle(CMS, sys,energy,mask):0);
2028 TGraphAsymmErrors* gALI =(alice ?GetSingle(ALICE, sys,energy,mask):0);
2029 TGraphAsymmErrors* gPYT =(pythia?GetSingle(PYTHIA,sys,energy,mask):0);
bfab35d9 2030 TGraphAsymmErrors* gWRK =(work ?GetSingle(WIP, sys,energy,mask):0);
797161e8 2031 if (gUAp) mp->Add(gUAp);
2032 if (gUAn) mp->Add(gUAn);
2033 if (gCMS) mp->Add(gCMS);
2034 if (gALI) mp->Add(gALI);
2035 if (gPYT) mp->Add(gPYT);
bfab35d9 2036 if (gWRK) mp->Add(gWRK);
797161e8 2037 if (gUAp || gUAn) seenUA5 = true;
2038 }
2039 if (seenUA5) sn.Append("(p#bar{p})");
7e4038b5 2040 }
e308a636 2041 else if (sys == 2) {
2042 // Nothing for PbPb so far
2043 cn = Form(", %d%%-%d%% central", centLow, centHigh);
797161e8 2044 sn = "PbPb";
5bb5d1f6 2045 // Warning("GetData", "No other data for PbPb yet");
7e4038b5 2046 }
6d233add 2047 else if (sys == 3) {
05cc5a69 2048 if (!(TMath::Abs(energy-5023) < 10 ||
2049 TMath::Abs(energy-8000) < 10) ) {
6d233add 2050 Warning("GetData", "No other results for sys=%d, energy=%d",
2051 sys, energy);
2052 return 0;
2053 }
2054
05cc5a69 2055 // Info("GetData", "Getting ALICE pPb data");
6d233add 2056 sn = "pPb";
2057 TGraphAsymmErrors* gALI =(alice ?GetSingle(ALICE, sys,energy, 0):0);
2058 if (gALI) mp->Add(gALI);
05cc5a69 2059 // Info("GetData", "Got %p", gALI);
6d233add 2060 // Warning("GetData", "Unknown system %d", sys);
2061 }
797161e8 2062
82fc512a 2063 if (!mp->GetListOfGraphs() || mp->GetListOfGraphs()->GetEntries() <= 0) {
2064 delete mp;
2065 mp = 0;
797161e8 2066 return 0;
82fc512a 2067 }
797161e8 2068 TString tit(Form("%s%s, %s%s",
2069 sn.Data(), en.Data(), tn.Data(), cn.Data()));
2070 mp->SetTitle(tit.Data());
7e4038b5 2071 return mp;
2072}
2073
2074//____________________________________________________________________
2075/**
2076 * Plot external data for a given selection of energy and trigger type
2077 * (see GetData)
2078 *
ffca499d 2079 * @param sys Collision system (1: pp, 2: PbPb)
2080 * @param energy Energy in GeV (900, 2360, 7000)
2081 * @param type Bit pattern of trigger type
2082 * - 0x1 INEL
2083 * - 0x2 INEL>0
2084 * - 0x4 NSD
2085 * @param centLow Low centrality cut (only for PbPB)
2086 * @param centHigh High centrality cut (only for PbPB)
c6115ede 2087 * @param which Which data to show
7e4038b5 2088 *
bd6f5206 2089 * @ingroup pwglf_forward_otherdata
7e4038b5 2090 */
2091void
e308a636 2092OtherData(UShort_t sys=1,
797161e8 2093 UShort_t energy=900,
2094 UShort_t type=0x1,
2095 UShort_t centLow=0,
2096 UShort_t centHigh=5,
bfab35d9 2097 UShort_t which=0xf)
7e4038b5 2098{
797161e8 2099 TMultiGraph* mp = GetData(sys, energy, type, centLow, centHigh, which);
82fc512a 2100 if (!mp) return;
2101
7e4038b5 2102 gStyle->SetTitleX(0.1);
2103 gStyle->SetTitleY(1.0);
2104 gStyle->SetTitleW(0.85);
2105 gStyle->SetTitleH(0.05);
2106 gStyle->SetTitleBorderSize(0);
2107 gStyle->SetTitleTextColor(kWhite);
2108 gStyle->SetTitleFillColor(kBlack);
2109 gStyle->SetTitleFontSize(0.02);
2110
797161e8 2111 gStyle->SetOptTitle(0);
7e4038b5 2112 gStyle->SetOptStat(0);
2113
2114 TCanvas* c = new TCanvas("c", "dN/deta", 800, 600);
2115 c->SetFillColor(0);
2116 c->SetBorderSize(0);
2117 c->SetBorderMode(0);
797161e8 2118 c->SetRightMargin(0.02);
2119 c->SetTopMargin(0.02);
7e4038b5 2120
82fc512a 2121
7e4038b5 2122 mp->SetMinimum(0);
2123 mp->Draw("ap");
797161e8 2124 if (mp->GetXaxis()) {
7e4038b5 2125 mp->GetXaxis()->SetTitle("#eta");
797161e8 2126 mp->GetXaxis()->SetTitleFont(12);
2127 mp->GetXaxis()->SetLabelFont(132);
2128 }
2129 if (mp->GetYaxis()) {
2130 mp->GetYaxis()->SetTitle("#frac{1}{N} #frac{dN_{#font[132]{ch}}}{d#eta}");
2131 mp->GetYaxis()->SetTitleFont(12);
2132 mp->GetYaxis()->SetLabelFont(132);
2133 }
2134 TLegend* l = c->BuildLegend(0.3, 0.15, 0.7, 0.5,
2135 mp->GetTitle());
7e4038b5 2136 l->SetFillColor(0);
797161e8 2137 l->SetTextFont(132);
7e4038b5 2138 l->SetBorderSize(0);
797161e8 2139 TLegendEntry* h = static_cast<TLegendEntry*>(l->GetListOfPrimitives()->At(0));
2140 if (h) {
2141 h->SetTextFont(22);
2142 }
7e4038b5 2143 c->cd();
2144}
2145
2146//____________________________________________________________________
2147//
2148// EOF
2149//