AliPhysics  64f4410 (64f4410)
ATLASpPb.C
Go to the documentation of this file.
1 //====================================================================
2 TGraph*
3 Get6090(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
4 {
5  TGraph* graph = new TGraph(54);
6  graph->SetName("g6090");
7  graph->SetTitle("60%-90%");
8  graph->SetFillColor(1);
9  graph->SetMarkerStyle(mStyle);
10  graph->SetMarkerSize(mSize);
11  graph->SetMarkerColor(mColor);;
12  graph->SetPoint(0,-2.65,9.322033898);
13  graph->SetPoint(1,-2.55,9.322033898);
14  graph->SetPoint(2,-2.45,9.322033898);
15  graph->SetPoint(3,-2.35,9.491525424);
16  graph->SetPoint(4,-2.25,9.491525424);
17  graph->SetPoint(5,-2.15,9.322033898);
18  graph->SetPoint(6,-2.05,9.491525424);
19  graph->SetPoint(7,-1.95,9.491525424);
20  graph->SetPoint(8,-1.85,9.491525424);
21  graph->SetPoint(9,-1.75,9.491525424);
22  graph->SetPoint(10,-1.65,9.491525424);
23  graph->SetPoint(11,-1.55,9.322033898);
24  graph->SetPoint(12,-1.45,9.152542373);
25  graph->SetPoint(13,-1.35,9.322033898);
26  graph->SetPoint(14,-1.25,9.152542373);
27  graph->SetPoint(15,-1.15,8.983050847);
28  graph->SetPoint(16,-1.05,8.983050847);
29  graph->SetPoint(17,-0.95,8.813559322);
30  graph->SetPoint(18,-0.85,8.813559322);
31  graph->SetPoint(19,-0.75,8.644067797);
32  graph->SetPoint(20,-0.65,8.644067797);
33  graph->SetPoint(21,-0.55,8.644067797);
34  graph->SetPoint(22,-0.45,8.474576271);
35  graph->SetPoint(23,-0.35,8.474576271);
36  graph->SetPoint(24,-0.25,8.305084746);
37  graph->SetPoint(25,-0.15,8.305084746);
38  graph->SetPoint(26,-0.05,8.13559322);
39  graph->SetPoint(27,0.05,8.305084746);
40  graph->SetPoint(28,0.15,8.305084746);
41  graph->SetPoint(29,0.25,8.305084746);
42  graph->SetPoint(30,0.35,8.305084746);
43  graph->SetPoint(31,0.45,8.305084746);
44  graph->SetPoint(32,0.55,8.644067797);
45  graph->SetPoint(33,0.65,8.474576271);
46  graph->SetPoint(34,0.75,8.644067797);
47  graph->SetPoint(35,0.85,8.644067797);
48  graph->SetPoint(36,0.95,8.644067797);
49  graph->SetPoint(37,1.05,8.813559322);
50  graph->SetPoint(38,1.15,8.813559322);
51  graph->SetPoint(39,1.25,8.813559322);
52  graph->SetPoint(40,1.35,8.983050847);
53  graph->SetPoint(41,1.45,9.152542373);
54  graph->SetPoint(42,1.55,8.983050847);
55  graph->SetPoint(43,1.65,9.152542373);
56  graph->SetPoint(44,1.75,9.152542373);
57  graph->SetPoint(45,1.85,8.983050847);
58  graph->SetPoint(46,1.95,9.152542373);
59  graph->SetPoint(47,2.05,8.983050847);
60  graph->SetPoint(48,2.15,8.983050847);
61  graph->SetPoint(49,2.25,8.983050847);
62  graph->SetPoint(50,2.35,8.983050847);
63  graph->SetPoint(51,2.45,8.813559322);
64  graph->SetPoint(52,2.55,8.644067797);
65  graph->SetPoint(53,2.65,8.813559322);
66 
67  return graph;
68 }
69 
70 //--------------------------------------------------------------------
71 TGraph*
72 Get4060(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
73 {
74  TGraph* graph = new TGraph(54);
75  graph->SetName("g4060");
76  graph->SetTitle("40%-60%");
77  graph->SetFillColor(1);
78  graph->SetMarkerStyle(mStyle);
79  graph->SetMarkerSize(mSize);
80  graph->SetMarkerColor(mColor);;
81  graph->SetPoint(0,-2.65,16.10169492);
82  graph->SetPoint(1,-2.55,16.61016949);
83  graph->SetPoint(2,-2.45,16.77966102);
84  graph->SetPoint(3,-2.35,16.94915254);
85  graph->SetPoint(4,-2.25,17.11864407);
86  graph->SetPoint(5,-2.15,17.11864407);
87  graph->SetPoint(6,-2.05,17.28813559);
88  graph->SetPoint(7,-1.95,17.11864407);
89  graph->SetPoint(8,-1.85,17.28813559);
90  graph->SetPoint(9,-1.75,17.28813559);
91  graph->SetPoint(10,-1.65,17.28813559);
92  graph->SetPoint(11,-1.55,17.28813559);
93  graph->SetPoint(12,-1.45,17.28813559);
94  graph->SetPoint(13,-1.35,17.28813559);
95  graph->SetPoint(14,-1.25,17.28813559);
96  graph->SetPoint(15,-1.15,17.11864407);
97  graph->SetPoint(16,-1.05,17.11864407);
98  graph->SetPoint(17,-0.95,16.94915254);
99  graph->SetPoint(18,-0.85,16.77966102);
100  graph->SetPoint(19,-0.75,16.77966102);
101  graph->SetPoint(20,-0.65,16.77966102);
102  graph->SetPoint(21,-0.55,16.44067797);
103  graph->SetPoint(22,-0.45,16.44067797);
104  graph->SetPoint(23,-0.35,16.27118644);
105  graph->SetPoint(24,-0.25,16.27118644);
106  graph->SetPoint(25,-0.15,16.10169492);
107  graph->SetPoint(26,-0.05,16.10169492);
108  graph->SetPoint(27,0.05,16.10169492);
109  graph->SetPoint(28,0.15,16.27118644);
110  graph->SetPoint(29,0.25,16.27118644);
111  graph->SetPoint(30,0.35,16.44067797);
112  graph->SetPoint(31,0.45,16.61016949);
113  graph->SetPoint(32,0.55,16.94915254);
114  graph->SetPoint(33,0.65,16.94915254);
115  graph->SetPoint(34,0.75,17.11864407);
116  graph->SetPoint(35,0.85,17.45762712);
117  graph->SetPoint(36,0.95,17.62711864);
118  graph->SetPoint(37,1.05,17.79661017);
119  graph->SetPoint(38,1.15,17.79661017);
120  graph->SetPoint(39,1.25,18.13559322);
121  graph->SetPoint(40,1.35,18.30508475);
122  graph->SetPoint(41,1.45,18.30508475);
123  graph->SetPoint(42,1.55,18.30508475);
124  graph->SetPoint(43,1.65,18.47457627);
125  graph->SetPoint(44,1.75,18.6440678);
126  graph->SetPoint(45,1.85,18.47457627);
127  graph->SetPoint(46,1.95,18.6440678);
128  graph->SetPoint(47,2.05,18.6440678);
129  graph->SetPoint(48,2.15,18.47457627);
130  graph->SetPoint(49,2.25,18.47457627);
131  graph->SetPoint(50,2.35,18.47457627);
132  graph->SetPoint(51,2.45,18.30508475);
133  graph->SetPoint(52,2.55,17.79661017);
134  graph->SetPoint(53,2.65,17.96610169);
135 
136  return graph;
137 }
138 
139 //--------------------------------------------------------------------
140 TGraph*
141 Get3040(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
142 {
143  TGraph* graph = new TGraph(54);
144  graph->SetName("g3040");
145  graph->SetTitle("30%-40%");
146  graph->SetFillColor(1);
147  graph->SetMarkerStyle(mStyle);
148  graph->SetMarkerSize(mSize);
149  graph->SetMarkerColor(mColor);;
150  graph->SetPoint(0,-2.65,20.6779661);
151  graph->SetPoint(1,-2.55,20.84745763);
152  graph->SetPoint(2,-2.45,21.18644068);
153  graph->SetPoint(3,-2.35,21.18644068);
154  graph->SetPoint(4,-2.25,21.52542373);
155  graph->SetPoint(5,-2.15,21.69491525);
156  graph->SetPoint(6,-2.05,21.86440678);
157  graph->SetPoint(7,-1.95,22.03389831);
158  graph->SetPoint(8,-1.85,22.20338983);
159  graph->SetPoint(9,-1.75,22.20338983);
160  graph->SetPoint(10,-1.65,22.37288136);
161  graph->SetPoint(11,-1.55,22.37288136);
162  graph->SetPoint(12,-1.45,22.54237288);
163  graph->SetPoint(13,-1.35,22.37288136);
164  graph->SetPoint(14,-1.25,22.37288136);
165  graph->SetPoint(15,-1.15,22.37288136);
166  graph->SetPoint(16,-1.05,22.37288136);
167  graph->SetPoint(17,-0.95,22.20338983);
168  graph->SetPoint(18,-0.85,22.03389831);
169  graph->SetPoint(19,-0.75,22.03389831);
170  graph->SetPoint(20,-0.65,21.69491525);
171  graph->SetPoint(21,-0.55,21.86440678);
172  graph->SetPoint(22,-0.45,21.69491525);
173  graph->SetPoint(23,-0.35,21.52542373);
174  graph->SetPoint(24,-0.25,21.52542373);
175  graph->SetPoint(25,-0.15,21.52542373);
176  graph->SetPoint(26,-0.05,21.3559322);
177  graph->SetPoint(27,0.05,21.69491525);
178  graph->SetPoint(28,0.15,21.69491525);
179  graph->SetPoint(29,0.25,21.86440678);
180  graph->SetPoint(30,0.35,22.20338983);
181  graph->SetPoint(31,0.45,22.20338983);
182  graph->SetPoint(32,0.55,22.71186441);
183  graph->SetPoint(33,0.65,22.71186441);
184  graph->SetPoint(34,0.75,23.22033898);
185  graph->SetPoint(35,0.85,23.38983051);
186  graph->SetPoint(36,0.95,23.55932203);
187  graph->SetPoint(37,1.05,23.89830508);
188  graph->SetPoint(38,1.15,24.06779661);
189  graph->SetPoint(39,1.25,24.40677966);
190  graph->SetPoint(40,1.35,24.57627119);
191  graph->SetPoint(41,1.45,24.74576271);
192  graph->SetPoint(42,1.55,24.91525424);
193  graph->SetPoint(43,1.65,25.08474576);
194  graph->SetPoint(44,1.75,25.25423729);
195  graph->SetPoint(45,1.85,25.25423729);
196  graph->SetPoint(46,1.95,25.25423729);
197  graph->SetPoint(47,2.05,25.08474576);
198  graph->SetPoint(48,2.15,25.25423729);
199  graph->SetPoint(49,2.25,25.25423729);
200  graph->SetPoint(50,2.35,25.08474576);
201  graph->SetPoint(51,2.45,24.74576271);
202  graph->SetPoint(52,2.55,24.23728814);
203  graph->SetPoint(53,2.65,24.57627119);
204 
205  return graph;
206 }
207 
208 //--------------------------------------------------------------------
209 TGraph*
210 Get2030(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
211 {
212  TGraph* graph = new TGraph(54);
213  graph->SetName("g2030");
214  graph->SetTitle("20%-30%");
215  graph->SetFillColor(1);
216  graph->SetMarkerStyle(mStyle);
217  graph->SetMarkerSize(mSize);
218  graph->SetMarkerColor(mColor);;
219 
220  graph->SetPoint(0,-2.65,23.38983051);
221  graph->SetPoint(1,-2.55,23.55932203);
222  graph->SetPoint(2,-2.45,23.89830508);
223  graph->SetPoint(3,-2.35,24.06779661);
224  graph->SetPoint(4,-2.25,24.40677966);
225  graph->SetPoint(5,-2.15,24.74576271);
226  graph->SetPoint(6,-2.05,24.91525424);
227  graph->SetPoint(7,-1.95,25.08474576);
228  graph->SetPoint(8,-1.85,25.42372881);
229  graph->SetPoint(9,-1.75,25.59322034);
230  graph->SetPoint(10,-1.65,25.93220339);
231  graph->SetPoint(11,-1.55,25.93220339);
232  graph->SetPoint(12,-1.45,26.10169492);
233  graph->SetPoint(13,-1.35,26.27118644);
234  graph->SetPoint(14,-1.25,26.10169492);
235  graph->SetPoint(15,-1.15,26.44067797);
236  graph->SetPoint(16,-1.05,26.27118644);
237  graph->SetPoint(17,-0.95,26.10169492);
238  graph->SetPoint(18,-0.85,26.10169492);
239  graph->SetPoint(19,-0.75,25.93220339);
240  graph->SetPoint(20,-0.65,25.76271186);
241  graph->SetPoint(21,-0.55,25.76271186);
242  graph->SetPoint(22,-0.45,25.59322034);
243  graph->SetPoint(23,-0.35,25.59322034);
244  graph->SetPoint(24,-0.25,25.42372881);
245  graph->SetPoint(25,-0.15,25.59322034);
246  graph->SetPoint(26,-0.05,25.59322034);
247  graph->SetPoint(27,0.05,25.76271186);
248  graph->SetPoint(28,0.15,25.93220339);
249  graph->SetPoint(29,0.25,26.10169492);
250  graph->SetPoint(30,0.35,26.27118644);
251  graph->SetPoint(31,0.45,26.77966102);
252  graph->SetPoint(32,0.55,27.11864407);
253  graph->SetPoint(33,0.65,27.28813559);
254  graph->SetPoint(34,0.75,27.62711864);
255  graph->SetPoint(35,0.85,28.13559322);
256  graph->SetPoint(36,0.95,28.47457627);
257  graph->SetPoint(37,1.05,28.6440678);
258  graph->SetPoint(38,1.15,29.3220339);
259  graph->SetPoint(39,1.25,29.15254237);
260  graph->SetPoint(40,1.35,29.66101695);
261  graph->SetPoint(41,1.45,29.83050847);
262  graph->SetPoint(42,1.55,30.16949153);
263  graph->SetPoint(43,1.65,30.33898305);
264  graph->SetPoint(44,1.75,30.6779661);
265  graph->SetPoint(45,1.85,30.33898305);
266  graph->SetPoint(46,1.95,30.50847458);
267  graph->SetPoint(47,2.05,30.50847458);
268  graph->SetPoint(48,2.15,30.6779661);
269  graph->SetPoint(49,2.25,30.50847458);
270  graph->SetPoint(50,2.35,30.33898305);
271  graph->SetPoint(51,2.45,30.33898305);
272  graph->SetPoint(52,2.55,29.66101695);
273  graph->SetPoint(53,2.65,30);
274 
275  return graph;
276 }
277 
278 
279 //--------------------------------------------------------------------
280 TGraph*
281 Get1020(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
282 {
283  TGraph* graph = new TGraph(54);
284  graph->SetName("g1020");
285  graph->SetTitle("10%-20%");
286  graph->SetFillColor(1);
287  graph->SetMarkerStyle(mStyle);
288  graph->SetMarkerSize(mSize);
289  graph->SetMarkerColor(mColor);;
290 
291  graph->SetPoint(0,-2.65,26.94915254);
292  graph->SetPoint(1,-2.55,27.28813559);
293  graph->SetPoint(2,-2.45,27.62711864);
294  graph->SetPoint(3,-2.35,28.30508475);
295  graph->SetPoint(4,-2.25,28.47457627);
296  graph->SetPoint(5,-2.15,28.81355932);
297  graph->SetPoint(6,-2.05,29.15254237);
298  graph->SetPoint(7,-1.95,29.66101695);
299  graph->SetPoint(8,-1.85,29.83050847);
300  graph->SetPoint(9,-1.75,30.33898305);
301  graph->SetPoint(10,-1.65,30.50847458);
302  graph->SetPoint(11,-1.55,30.84745763);
303  graph->SetPoint(12,-1.45,30.84745763);
304  graph->SetPoint(13,-1.35,31.18644068);
305  graph->SetPoint(14,-1.25,31.18644068);
306  graph->SetPoint(15,-1.15,31.18644068);
307  graph->SetPoint(16,-1.05,31.18644068);
308  graph->SetPoint(17,-0.95,31.18644068);
309  graph->SetPoint(18,-0.85,31.01694915);
310  graph->SetPoint(19,-0.75,31.01694915);
311  graph->SetPoint(20,-0.65,31.01694915);
312  graph->SetPoint(21,-0.55,30.84745763);
313  graph->SetPoint(22,-0.45,30.84745763);
314  graph->SetPoint(23,-0.35,30.84745763);
315  graph->SetPoint(24,-0.25,30.84745763);
316  graph->SetPoint(25,-0.15,30.84745763);
317  graph->SetPoint(26,-0.05,31.01694915);
318  graph->SetPoint(27,0.05,31.01694915);
319  graph->SetPoint(28,0.15,31.3559322);
320  graph->SetPoint(29,0.25,31.52542373);
321  graph->SetPoint(30,0.35,32.03389831);
322  graph->SetPoint(31,0.45,32.37288136);
323  graph->SetPoint(32,0.55,32.88135593);
324  graph->SetPoint(33,0.65,33.22033898);
325  graph->SetPoint(34,0.75,33.55932203);
326  graph->SetPoint(35,0.85,34.23728814);
327  graph->SetPoint(36,0.95,34.57627119);
328  graph->SetPoint(37,1.05,35.25423729);
329  graph->SetPoint(38,1.15,35.42372881);
330  graph->SetPoint(39,1.25,35.93220339);
331  graph->SetPoint(40,1.35,36.27118644);
332  graph->SetPoint(41,1.45,36.61016949);
333  graph->SetPoint(42,1.55,37.11864407);
334  graph->SetPoint(43,1.65,36.94915254);
335  graph->SetPoint(44,1.75,37.28813559);
336  graph->SetPoint(45,1.85,37.28813559);
337  graph->SetPoint(46,1.95,37.45762712);
338  graph->SetPoint(47,2.05,37.62711864);
339  graph->SetPoint(48,2.15,37.45762712);
340  graph->SetPoint(49,2.25,37.45762712);
341  graph->SetPoint(50,2.35,37.28813559);
342  graph->SetPoint(51,2.45,37.11864407);
343  graph->SetPoint(52,2.55,36.27118644);
344  graph->SetPoint(53,2.65,36.94915254);
345 
346  return graph;
347 }
348 
349 //--------------------------------------------------------------------
350 TGraph*
351 Get0510(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
352 {
353  TGraph* graph = new TGraph(54);
354  graph->SetName("g0510");
355  graph->SetTitle("5%-10%");
356  graph->SetFillColor(1);
357  graph->SetMarkerStyle(mStyle);
358  graph->SetMarkerSize(mSize);
359  graph->SetMarkerColor(mColor);;
360 
361  graph->SetPoint(0,-2.65,30.33898305);
362  graph->SetPoint(1,-2.55,31.3559322);
363  graph->SetPoint(2,-2.45,31.69491525);
364  graph->SetPoint(3,-2.35,31.86440678);
365  graph->SetPoint(4,-2.25,32.54237288);
366  graph->SetPoint(5,-2.15,33.22033898);
367  graph->SetPoint(6,-2.05,33.72881356);
368  graph->SetPoint(7,-1.95,34.06779661);
369  graph->SetPoint(8,-1.85,34.40677966);
370  graph->SetPoint(9,-1.75,34.74576271);
371  graph->SetPoint(10,-1.65,35.42372881);
372  graph->SetPoint(11,-1.55,35.59322034);
373  graph->SetPoint(12,-1.45,35.93220339);
374  graph->SetPoint(13,-1.35,36.10169492);
375  graph->SetPoint(14,-1.25,36.27118644);
376  graph->SetPoint(15,-1.15,36.44067797);
377  graph->SetPoint(16,-1.05,36.61016949);
378  graph->SetPoint(17,-0.95,36.44067797);
379  graph->SetPoint(18,-0.85,36.61016949);
380  graph->SetPoint(19,-0.75,36.44067797);
381  graph->SetPoint(20,-0.65,36.61016949);
382  graph->SetPoint(21,-0.55,36.44067797);
383  graph->SetPoint(22,-0.45,36.27118644);
384  graph->SetPoint(23,-0.35,36.61016949);
385  graph->SetPoint(24,-0.25,36.61016949);
386  graph->SetPoint(25,-0.15,36.61016949);
387  graph->SetPoint(26,-0.05,36.94915254);
388  graph->SetPoint(27,0.05,36.77966102);
389  graph->SetPoint(28,0.15,37.28813559);
390  graph->SetPoint(29,0.25,37.79661017);
391  graph->SetPoint(30,0.35,38.30508475);
392  graph->SetPoint(31,0.45,38.98305085);
393  graph->SetPoint(32,0.55,39.3220339);
394  graph->SetPoint(33,0.65,39.83050847);
395  graph->SetPoint(34,0.75,40.6779661);
396  graph->SetPoint(35,0.85,41.18644068);
397  graph->SetPoint(36,0.95,41.86440678);
398  graph->SetPoint(37,1.05,42.20338983);
399  graph->SetPoint(38,1.15,42.88135593);
400  graph->SetPoint(39,1.25,43.55932203);
401  graph->SetPoint(40,1.35,43.72881356);
402  graph->SetPoint(41,1.45,44.23728814);
403  graph->SetPoint(42,1.55,44.74576271);
404  graph->SetPoint(43,1.65,45.25423729);
405  graph->SetPoint(44,1.75,45.25423729);
406  graph->SetPoint(45,1.85,45.42372881);
407  graph->SetPoint(46,1.95,45.76271186);
408  graph->SetPoint(47,2.05,45.76271186);
409  graph->SetPoint(48,2.15,45.76271186);
410  graph->SetPoint(49,2.25,45.42372881);
411  graph->SetPoint(50,2.35,45.42372881);
412  graph->SetPoint(51,2.45,45.42372881);
413  graph->SetPoint(52,2.55,44.74576271);
414  graph->SetPoint(53,2.65,45.08474576);
415 
416  return graph;
417 }
418 
419 //--------------------------------------------------------------------
420 TGraph*
421 Get0105(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
422 {
423  TGraph* graph = new TGraph(54);
424  graph->SetName("g0105");
425  graph->SetTitle("1%-5%");
426  graph->SetFillColor(1);
427  graph->SetMarkerStyle(mStyle);
428  graph->SetMarkerSize(mSize);
429  graph->SetMarkerColor(mColor);;
430 
431  graph->SetPoint(0,-2.65,33.72881356);
432  graph->SetPoint(1,-2.55,34.91525424);
433  graph->SetPoint(2,-2.45,35.59322034);
434  graph->SetPoint(3,-2.35,36.44067797);
435  graph->SetPoint(4,-2.25,37.28813559);
436  graph->SetPoint(5,-2.15,37.79661017);
437  graph->SetPoint(6,-2.05,38.47457627);
438  graph->SetPoint(7,-1.95,39.15254237);
439  graph->SetPoint(8,-1.85,39.83050847);
440  graph->SetPoint(9,-1.75,40.50847458);
441  graph->SetPoint(10,-1.65,41.01694915);
442  graph->SetPoint(11,-1.55,41.3559322);
443  graph->SetPoint(12,-1.45,41.69491525);
444  graph->SetPoint(13,-1.35,42.20338983);
445  graph->SetPoint(14,-1.25,42.37288136);
446  graph->SetPoint(15,-1.15,42.71186441);
447  graph->SetPoint(16,-1.05,42.71186441);
448  graph->SetPoint(17,-0.95,42.88135593);
449  graph->SetPoint(18,-0.85,42.71186441);
450  graph->SetPoint(19,-0.75,42.88135593);
451  graph->SetPoint(20,-0.65,42.88135593);
452  graph->SetPoint(21,-0.55,42.88135593);
453  graph->SetPoint(22,-0.45,43.05084746);
454  graph->SetPoint(23,-0.35,43.22033898);
455  graph->SetPoint(24,-0.25,43.22033898);
456  graph->SetPoint(25,-0.15,43.38983051);
457  graph->SetPoint(26,-0.05,43.72881356);
458  graph->SetPoint(27,0.05,44.06779661);
459  graph->SetPoint(28,0.15,44.57627119);
460  graph->SetPoint(29,0.25,45.08474576);
461  graph->SetPoint(30,0.35,45.76271186);
462  graph->SetPoint(31,0.45,46.61016949);
463  graph->SetPoint(32,0.55,46.77966102);
464  graph->SetPoint(33,0.65,48.13559322);
465  graph->SetPoint(34,0.75,48.81355932);
466  graph->SetPoint(35,0.85,49.49152542);
467  graph->SetPoint(36,0.95,50.16949153);
468  graph->SetPoint(37,1.05,51.01694915);
469  graph->SetPoint(38,1.15,51.86440678);
470  graph->SetPoint(39,1.25,52.20338983);
471  graph->SetPoint(40,1.35,53.05084746);
472  graph->SetPoint(41,1.45,53.22033898);
473  graph->SetPoint(42,1.55,54.06779661);
474  graph->SetPoint(43,1.65,54.57627119);
475  graph->SetPoint(44,1.75,54.57627119);
476  graph->SetPoint(45,1.85,55.08474576);
477  graph->SetPoint(46,1.95,55.42372881);
478  graph->SetPoint(47,2.05,55.25423729);
479  graph->SetPoint(48,2.15,55.25423729);
480  graph->SetPoint(49,2.25,55.42372881);
481  graph->SetPoint(50,2.35,55.25423729);
482  graph->SetPoint(51,2.45,55.25423729);
483  graph->SetPoint(52,2.55,54.57627119);
484  graph->SetPoint(53,2.65,54.74576271);
485 
486  return graph;
487 }
488 
489 //--------------------------------------------------------------------
490 TGraph*
491 Get0001(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
492 {
493  TGraph* graph = new TGraph(54);
494  graph->SetName("g0001");
495  graph->SetTitle("0%-1%");
496  graph->SetFillColor(1);
497  graph->SetMarkerStyle(mStyle);
498  graph->SetMarkerSize(mSize);
499  graph->SetMarkerColor(mColor);;
500 
501  graph->SetPoint(0,-2.65,40.33898305);
502  graph->SetPoint(1,-2.55,41.52542373);
503  graph->SetPoint(2,-2.45,43.05084746);
504  graph->SetPoint(3,-2.35,43.89830508);
505  graph->SetPoint(4,-2.25,44.74576271);
506  graph->SetPoint(5,-2.15,45.76271186);
507  graph->SetPoint(6,-2.05,46.27118644);
508  graph->SetPoint(7,-1.95,47.28813559);
509  graph->SetPoint(8,-1.85,48.81355932);
510  graph->SetPoint(9,-1.75,49.66101695);
511  graph->SetPoint(10,-1.65,50.33898305);
512  graph->SetPoint(11,-1.55,50.84745763);
513  graph->SetPoint(12,-1.45,51.52542373);
514  graph->SetPoint(13,-1.35,51.86440678);
515  graph->SetPoint(14,-1.25,52.71186441);
516  graph->SetPoint(15,-1.15,52.71186441);
517  graph->SetPoint(16,-1.05,53.05084746);
518  graph->SetPoint(17,-0.95,53.55932203);
519  graph->SetPoint(18,-0.85,53.55932203);
520  graph->SetPoint(19,-0.75,53.55932203);
521  graph->SetPoint(20,-0.65,53.89830508);
522  graph->SetPoint(21,-0.55,54.23728814);
523  graph->SetPoint(22,-0.45,54.57627119);
524  graph->SetPoint(23,-0.35,54.57627119);
525  graph->SetPoint(24,-0.25,54.74576271);
526  graph->SetPoint(25,-0.15,54.74576271);
527  graph->SetPoint(26,-0.05,55.08474576);
528  graph->SetPoint(27,0.05,55.93220339);
529  graph->SetPoint(28,0.15,56.77966102);
530  graph->SetPoint(29,0.25,57.28813559);
531  graph->SetPoint(30,0.35,58.6440678);
532  graph->SetPoint(31,0.45,58.98305085);
533  graph->SetPoint(32,0.55,60.50847458);
534  graph->SetPoint(33,0.65,61.69491525);
535  graph->SetPoint(34,0.75,62.88135593);
536  graph->SetPoint(35,0.85,63.05084746);
537  graph->SetPoint(36,0.95,64.40677966);
538  graph->SetPoint(37,1.05,64.91525424);
539  graph->SetPoint(38,1.15,66.77966102);
540  graph->SetPoint(39,1.25,67.28813559);
541  graph->SetPoint(40,1.35,68.30508475);
542  graph->SetPoint(41,1.45,70.16949153);
543  graph->SetPoint(42,1.55,69.66101695);
544  graph->SetPoint(43,1.65,70.84745763);
545  graph->SetPoint(44,1.75,71.69491525);
546  graph->SetPoint(45,1.85,72.54237288);
547  graph->SetPoint(46,1.95,71.86440678);
548  graph->SetPoint(47,2.05,72.20338983);
549  graph->SetPoint(48,2.15,72.54237288);
550  graph->SetPoint(49,2.25,72.88135593);
551  graph->SetPoint(50,2.35,72.20338983);
552  graph->SetPoint(51,2.45,72.37288136);
553  graph->SetPoint(52,2.55,72.20338983);
554  graph->SetPoint(53,2.65,72.37288136);
555 
556  return graph;
557 }
558 
559 
560 //====================================================================
562 {
563  TGraphErrors* ge = new TGraphErrors(g->GetN());
564  ge->SetName(g->GetName());
565  ge->SetTitle(g->GetTitle());
566  ge->SetMarkerSize(g->GetMarkerSize());
567  ge->SetMarkerColor(g->GetMarkerColor());
568  ge->SetMarkerStyle(g->GetMarkerStyle());
569  ge->SetFillColor(kBlue-10);
570  ge->SetFillStyle(3002);
571 
572  for (Int_t i = 0; i < g->GetN(); i++) {
573  ge->SetPoint(i, g->GetX()[i], g->GetY()[i]);
574  ge->SetPointError(i, 0.05, e * g->GetY()[i]);
575  }
576  return ge;
577 }
578 
579 //====================================================================
580 TMultiGraph*
581 ATLASpPb(Bool_t draw=false, TLegend* l=0)
582 {
583  TCanvas* c = 0;
584  if (draw) {
585  c = new TCanvas("C","C");
586  c->SetTopMargin(0.02);
587  c->SetRightMargin(0.02);
588  if (!l) {
589  l = new TLegend(0.15, 0.7, 0.4, 0.95);
590  l->SetFillColor(0);
591  l->SetFillStyle(0);
592  l->SetBorderSize(0);
593  }
594  }
595 
596  Double_t mSize = 1.6;
597  Int_t mStyle = 22;
598 
599  TMultiGraph* mg = new TMultiGraph("atlas", "ATLAS");
600  TGraph* g0001 = GetGE(Get0001(mSize,mStyle,kViolet)); // 0-1
601  TGraph* g0105 = GetGE(Get0105(mSize,mStyle,52)); // 1-5
602  TGraph* g0510 = GetGE(Get0510(mSize,mStyle,55)); // 5-10
603  TGraph* g1020 = GetGE(Get1020(mSize,mStyle,59)); // 10-20
604  TGraph* g2030 = GetGE(Get2030(mSize,mStyle,64)); // 20-30
605  TGraph* g3040 = GetGE(Get3040(mSize,mStyle,kTeal-5)); // 30-40
606  TGraph* g4060 = GetGE(Get4060(mSize,mStyle,76), 0.06); // 40-60
607  TGraph* g6090 = GetGE(Get6090(mSize,mStyle,kYellow-3), 0.1); // 60-90
608 
609 
610  mg->Add(g0001); if (l) l->AddEntry(g0001,g0001->GetTitle(),"pl");
611  mg->Add(g0105); if (l) l->AddEntry(g0105,g0105->GetTitle(),"pl");
612  mg->Add(g0510); if (l) l->AddEntry(g0510,g0510->GetTitle(),"pl");
613  mg->Add(g1020); if (l) l->AddEntry(g1020,g1020->GetTitle(),"pl");
614  mg->Add(g2030); if (l) l->AddEntry(g2030,g2030->GetTitle(),"pl");
615  mg->Add(g3040); if (l) l->AddEntry(g3040,g3040->GetTitle(),"pl");
616  mg->Add(g4060); if (l) l->AddEntry(g4060,g4060->GetTitle(),"pl");
617  mg->Add(g6090); if (l) l->AddEntry(g6090,g6090->GetTitle(),"pl");
618 
619  if (draw) mg->Draw("ap");
620  if (l) l->Draw();
621 
622  return mg;
623 }
624 
625 //====================================================================
626 TGraphErrors*
627 Average(const char* name,
628  const char* title,
629  TGraphErrors* g1, Double_t w1,
630  TGraphErrors* g2, Double_t w2)
631 {
632  TGraphErrors* ret = new TGraphErrors(g1->GetN());
633  ret->SetName(name);
634  ret->SetTitle(title);
635  ret->SetMarkerColor(g1->GetMarkerColor());
636  ret->SetLineColor(kBlack);
637  ret->SetMarkerStyle(g1->GetMarkerStyle());
638  ret->SetMarkerSize(g1->GetMarkerSize());
639  ret->SetFillStyle(0);
640  for (Int_t i = 0; i < g1->GetN(); i++) {
641  Double_t y1 = w1 * g1->GetY()[i];
642  Double_t y2 = w2 * g2->GetY()[i];
643  Double_t e1 = w1 * g1->GetEY()[i];
644  Double_t e2 = w2 * g2->GetEY()[i];
645  Double_t y = (y1 + y2) / (w1+w2);
646  Double_t ee = TMath::Sqrt(e1*e1+e2*e2);
647  ret->SetPoint(i, g1->GetX()[i], y);
648  ret->SetPointError(i, g1->GetEX()[i], ee);
649  }
650  return ret;
651 }
652 //____________________________________________________________________
653 TH1*
654 Average(const char* name,
655  const char* title,
656  TH1* g1, Double_t w1,
657  TH1* g2, Double_t w2)
658 {
659  TH1* ret = static_cast<TH1*>(g1->Clone(name));
660  ret->SetTitle(title);
661  ret->SetLineColor(kBlack);
662  ret->SetDirectory(0);
663  ret->Reset();
664  for (Int_t i = 0; i < g1->GetNbinsX(); i++) {
665  Double_t y1 = w1 * g1->GetBinContent(i);
666  Double_t y2 = w2 * g2->GetBinContent(i);
667  Double_t e1 = w1 * g1->GetBinError(i);
668  Double_t e2 = w2 * g2->GetBinError(i);
669  Double_t y = (y1 + y2) / (w1+w2);
670  Double_t ee = TMath::Sqrt(e1*e1+e2*e2);
671  ret->SetBinContent(i, y);
672  ret->SetBinError(i, ee);
673  }
674  return ret;
675 }
676 //____________________________________________________________________
678 Average(const char* title,
679  TList* gl, Int_t idx1, Double_t w1, Int_t idx2, Double_t w2)
680 {
681  TGraphErrors* g1 = static_cast<TGraphErrors*>(gl->At(idx1));
682  TGraphErrors* g2 = static_cast<TGraphErrors*>(gl->At(idx2));
683 
684  return Average(g1->GetName(), title, g1, w1, g2, w2);
685 }
686 
687 //====================================================================
688 TMultiGraph*
689 ATLASpPbResum(Bool_t draw=false, TLegend* l=0)
690 {
691  TCanvas* c = 0;
692  if (draw) {
693  c = new TCanvas("C","C");
694  c->SetTopMargin(0.02);
695  c->SetRightMargin(0.02);
696  if (!l) {
697  l = new TLegend(0.15, 0.7, 0.4, 0.95);
698  l->SetFillColor(0);
699  l->SetFillStyle(0);
700  l->SetBorderSize(0);
701  }
702  }
703 
704  Option_t* op = "p 5 z";
705  TMultiGraph* mg = ATLASpPb(false);
706  TList* gl = mg->GetListOfGraphs();
707  TMultiGraph* rm = new TMultiGraph;
708 
709  TGraphErrors* ge = 0;
710  // Add 0-1% to 1-5%
711  ge = Average("0-5%", gl, 0, .2, 1, .8);
712  rm->Add(ge, op);
713  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
714  // 5-10%
715  ge = static_cast<TGraphErrors*>(gl->At(2));
716  ge->SetFillStyle(0);
717  rm->Add(ge, op);
718  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
719  // 10-20%
720  ge = static_cast<TGraphErrors*>(gl->At(3));
721  ge->SetFillStyle(0);
722  rm->Add(ge, op);
723  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
724  // Add 20-30% to 30-40%
725  ge = Average("20-40%", gl, 4, .5, 5, .5);
726  rm->Add(ge, op);
727  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
728  // 40-60%
729  ge = static_cast<TGraphErrors*>(gl->At(6));
730  ge->SetFillStyle(0);
731  rm->Add(ge, op);
732  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
733  // 60-90%
734  ge = static_cast<TGraphErrors*>(gl->At(7));
735  ge->SetFillStyle(0);
736  rm->Add(ge, op);
737  if (l) l->AddEntry(ge, ge->GetTitle(), "p");
738 
739  if (draw) rm->Draw("ap");
740  if (l) l->Draw();
741 
742  return rm;
743 }
744 
745 
746 // EOF
TGraph * Get6090(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:3
TGraph * Get0001(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:491
TMultiGraph * ATLASpPb(Bool_t draw=false, TLegend *l=0)
Definition: ATLASpPb.C:581
double Double_t
Definition: External.C:58
TGraph * Get2030(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:210
const char * title
Definition: MakeQAPdf.C:27
TGraph * Get0510(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:351
TGraph * Get4060(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:72
TCanvas * c
Definition: TestFitELoss.C:172
TMultiGraph * ATLASpPbResum(Bool_t draw=false, TLegend *l=0)
Definition: ATLASpPb.C:689
TGraph * Get0105(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:421
TGraph * Get3040(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:141
int Int_t
Definition: External.C:63
TGraphErrors * Average(const char *name, const char *title, TGraphErrors *g1, Double_t w1, TGraphErrors *g2, Double_t w2)
Definition: ATLASpPb.C:627
Bool_t draw[nPtBins]
const char Option_t
Definition: External.C:48
TGraph * Get1020(Double_t mSize=1.2, Int_t mStyle=22, Color_t mColor=kBlack)
Definition: ATLASpPb.C:281
bool Bool_t
Definition: External.C:53
Definition: External.C:196
TGraph * GetGE(TGraph *g, Double_t e=0.045)
Definition: ATLASpPb.C:561