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