9c3516c97d4b781645c5263631f7d32311c2a41f
[u/mrichter/AliRoot.git] / PWG2 / SPECTRA / Fit / StarPPSpectra.C
1 TGraphErrors ** StarPPSpectra() {
2
3   TGraphErrors ** gStar = new TGraphErrors*[6];
4   const Int_t kNMaxPoints = 19;
5   Double_t pt[kNMaxPoints];
6   Double_t errpt[kNMaxPoints] = {0};
7   Double_t pim [kNMaxPoints];
8   Double_t pime[kNMaxPoints];
9   Double_t pip [kNMaxPoints];
10   Double_t pipe[kNMaxPoints];
11   Double_t km [kNMaxPoints];
12   Double_t kme[kNMaxPoints];
13   Double_t kp [kNMaxPoints];
14   Double_t kpe[kNMaxPoints];
15   Double_t pm [kNMaxPoints];
16   Double_t pme[kNMaxPoints];
17   Double_t pp [kNMaxPoints];
18   Double_t ppe[kNMaxPoints];
19
20   Int_t ipoint = 0;
21   pt[ipoint++] = 0.225;
22   pt[ipoint++] = 0.275;
23   pt[ipoint++] = 0.325;
24   pt[ipoint++] = 0.375;
25   pt[ipoint++] = 0.425;
26   pt[ipoint++] = 0.475;
27   pt[ipoint++] = 0.525;
28   pt[ipoint++] = 0.575;
29   pt[ipoint++] = 0.625;
30   pt[ipoint++] = 0.675;
31   pt[ipoint++] = 0.725;
32   pt[ipoint++] = 0.775;
33   pt[ipoint++] = 0.825;
34   pt[ipoint++] = 0.875;
35   pt[ipoint++] = 0.925;
36   pt[ipoint++] = 0.975;
37   pt[ipoint++] = 1.025;
38   pt[ipoint++] = 1.075;
39   pt[ipoint++] = 1.125;
40
41
42   // pi-
43   ipoint = 0;
44   pim[ipoint++]  = 2.02;
45   pim[ipoint++]  = 1.52;
46   pim[ipoint++]  = 1.13;
47   pim[ipoint++]  = 8.44*1e-1;
48   pim[ipoint++]  = 6.35*1e-1;
49   pim[ipoint++]  = 4.69*1e-1;
50   pim[ipoint++]  = 3.54*1e-1;
51   pim[ipoint++]  = 2.67*1e-1;
52   pim[ipoint++]  = 2.02*1e-1;
53   pim[ipoint++]  = 1.53*1e-1;
54   pim[ipoint++]  = 1.16*1e-1;
55   pim[ipoint++]  = 8.65*1e-2;
56
57   ipoint = 0;
58   pime[ipoint++] = 0.06;
59   pime[ipoint++] = 0.03;
60   pime[ipoint++] = 0.02;
61   pime[ipoint++] = 0.09*1e-1;
62   pime[ipoint++] = 0.07*1e-1;
63   pime[ipoint++] = 0.05*1e-1;
64   pime[ipoint++] = 0.04*1e-1;
65   pime[ipoint++] = 0.03*1e-1;
66   pime[ipoint++] = 0.04*1e-1;
67   pime[ipoint++] = 0.03*1e-1;
68   pime[ipoint++] = 0.03*1e-1;
69   pime[ipoint++] = 0.30*1e-2;
70
71   gStar[0] = new TGraphErrors(ipoint,pt,pim,errpt,pime);
72
73   // pi+
74   ipoint = 0;
75   pip[ipoint++]  = 2.07;
76   pip[ipoint++]  = 1.54;
77   pip[ipoint++]  = 1.14;
78   pip[ipoint++]  = 8.57*1e-1;
79   pip[ipoint++]  = 6.38*1e-1;
80   pip[ipoint++]  = 4.76*1e-1;
81   pip[ipoint++]  = 3.59*1e-1;
82   pip[ipoint++]  = 2.73*1e-1;
83   pip[ipoint++]  = 2.07*1e-1;
84   pip[ipoint++]  = 1.55*1e-1;
85   pip[ipoint++]  = 1.16*1e-1;
86   pip[ipoint++]  = 8.95*1e-2;
87
88   ipoint = 0;
89   pipe[ipoint++] = 0.06;
90   pipe[ipoint++] = 0.03;
91   pipe[ipoint++] = 0.02;
92   pipe[ipoint++] = 0.09*1e-1;
93   pipe[ipoint++] = 0.07*1e-1;
94   pipe[ipoint++] = 0.05*1e-1;
95   pipe[ipoint++] = 0.04*1e-1;
96   pipe[ipoint++] = 0.03*1e-1;
97   pipe[ipoint++] = 0.04*1e-1;
98   pipe[ipoint++] = 0.03*1e-1;
99   pipe[ipoint++] = 0.03*1e-1;
100   pipe[ipoint++] = 0.31*1e-2;
101
102   gStar[1] = new TGraphErrors(ipoint,pt,pip,errpt,pipe);
103
104   // k-
105   ipoint = 0;
106   km[ipoint++]  = 1.43*1e-1;
107   km[ipoint++]  = 1.26*1e-1;
108   km[ipoint++]  = 1.08*1e-1;
109   km[ipoint++]  = 8.77*1e-2;
110   km[ipoint++]  = 7.34*1e-2;
111   km[ipoint++]  = 6.17*1e-2;
112   km[ipoint++]  = 4.87*1e-2;
113   km[ipoint++]  = 4.11*1e-2;
114   km[ipoint++]  = 3.70*1e-2;
115   km[ipoint++]  = 2.86*1e-2;
116   km[ipoint++]  = 2.42*1e-2;
117
118   ipoint = 0;
119   kme[ipoint++] = 0.11*1e-1;
120   kme[ipoint++] = 0.05*1e-1;
121   kme[ipoint++] = 0.02*1e-1;
122   kme[ipoint++] = 0.24*1e-2;
123   kme[ipoint++] = 0.26*1e-2;
124   kme[ipoint++] = 0.63*1e-2;
125   kme[ipoint++] = 0.50*1e-2;
126   kme[ipoint++] = 0.42*1e-2;
127   kme[ipoint++] = 0.28*1e-2;
128   kme[ipoint++] = 0.22*1e-2;
129   kme[ipoint++] = 0.27*1e-2;
130
131   gStar[2] = new TGraphErrors(ipoint,pt,km,errpt,kme);
132
133
134   // k+
135   ipoint = 0;
136   kp[ipoint++]  = 1.52*1e-1;
137   kp[ipoint++]  = 1.30*1e-1;
138   kp[ipoint++]  = 1.08*1e-1;
139   kp[ipoint++]  = 9.16*1e-2;
140   kp[ipoint++]  = 7.47*1e-2;
141   kp[ipoint++]  = 6.26*1e-2;
142   kp[ipoint++]  = 5.26*1e-2;
143   kp[ipoint++]  = 4.41*1e-2;
144   kp[ipoint++]  = 3.81*1e-2;
145   kp[ipoint++]  = 3.06*1e-2;
146   kp[ipoint++]  = 2.49*1e-2;
147
148   ipoint = 0;
149   kpe[ipoint++] = 0.11*1e-1;
150   kpe[ipoint++] = 0.05*1e-1;
151   kpe[ipoint++] = 0.02*1e-1;
152   kpe[ipoint++] = 0.25*1e-2;
153   kpe[ipoint++] = 0.27*1e-2;
154   kpe[ipoint++] = 0.64*1e-2;
155   kpe[ipoint++] = 0.54*1e-2;
156   kpe[ipoint++] = 0.45*1e-2;
157   kpe[ipoint++] = 0.28*1e-2;
158   kpe[ipoint++] = 0.24*1e-2;
159   kpe[ipoint++] = 0.28*1e-2;
160
161   gStar[3] = new TGraphErrors(ipoint,pt,kp,errpt,kpe);
162
163
164   // pbar
165   ipoint = 0;
166   pm[ipoint++]  = 0;
167   pm[ipoint++]  = 0;
168   pm[ipoint++]  = 0;
169   pm[ipoint++]  = 5.54*1e-2;
170   pm[ipoint++]  = 4.81*1e-2;
171   pm[ipoint++]  = 4.25*1e-2;
172   pm[ipoint++]  = 3.77*1e-2;
173   pm[ipoint++]  = 3.38*1e-2;
174   pm[ipoint++]  = 2.78*1e-2;
175   pm[ipoint++]  = 2.49*1e-2;
176   pm[ipoint++]  = 2.03*1e-2;
177   pm[ipoint++]  = 1.75*1e-2;
178   pm[ipoint++]  = 1.52*1e-2;
179   pm[ipoint++]  = 1.27*1e-2;
180   pm[ipoint++]  = 1.05*1e-2;
181   pm[ipoint++]  = 8.95*1e-3;
182   pm[ipoint++]  = 7.36*1e-3;
183   pm[ipoint++]  = 6.59*1e-3;
184   pm[ipoint++]  = 5.21*1e-3;
185
186   ipoint = 0;
187   pme[ipoint++] = 0;
188   pme[ipoint++] = 0;
189   pme[ipoint++] = 0;
190   pme[ipoint++] = 0.13*1e-2;
191   pme[ipoint++] = 0.11*1e-2;
192   pme[ipoint++] = 0.10*1e-2;
193   pme[ipoint++] = 0.09*1e-2;
194   pme[ipoint++] = 0.08*1e-2;
195   pme[ipoint++] = 0.07*1e-2;
196   pme[ipoint++] = 0.06*1e-2;
197   pme[ipoint++] = 0.06*1e-2;
198   pme[ipoint++] = 0.06*1e-2;
199   pme[ipoint++] = 0.05*1e-2;
200   pme[ipoint++] = 0.04*1e-2;
201   pme[ipoint++] = 0.04*1e-2;
202   pme[ipoint++] = 0.34*1e-3;
203   pme[ipoint++] = 0.34*1e-3;
204   pme[ipoint++] = 0.32*1e-3;
205   pme[ipoint++] = 0.29*1e-3;
206
207   gStar[4] =  new TGraphErrors(ipoint,pt,pm,errpt,pme);
208
209
210   // p
211   ipoint = 0;
212   pp[ipoint++]  = 0;
213   pp[ipoint++]  = 0;
214   pp[ipoint++]  = 0;
215   pp[ipoint++]  = 0;
216   pp[ipoint++]  = 0;
217   pp[ipoint++]  = 5.07*1e-2;
218   pp[ipoint++]  = 4.69*1e-2;
219   pp[ipoint++]  = 4.08*1e-2;
220   pp[ipoint++]  = 3.42*1e-2;
221   pp[ipoint++]  = 2.87*1e-2;
222   pp[ipoint++]  = 2.41*1e-2;
223   pp[ipoint++]  = 2.13*1e-2;
224   pp[ipoint++]  = 1.82*1e-2;
225   pp[ipoint++]  = 1.54*1e-2;
226   pp[ipoint++]  = 1.31*1e-2;
227   pp[ipoint++]  = 1.11*1e-2;
228   pp[ipoint++]  = 9.78*1e-3;
229   pp[ipoint++]  = 8.56*1e-3;
230   pp[ipoint++]  = 7.38*1e-3;
231
232   ipoint = 0;
233   ppe[ipoint++] = 0;
234   ppe[ipoint++] = 0;
235   ppe[ipoint++] = 0;
236   ppe[ipoint++] = 0;
237   ppe[ipoint++] = 0;
238   ppe[ipoint++] = 0.25*1e-2;
239   ppe[ipoint++] = 0.19*1e-2;
240   ppe[ipoint++] = 0.14*1e-2;
241   ppe[ipoint++] = 0.10*1e-2;
242   ppe[ipoint++] = 0.07*1e-2;
243   ppe[ipoint++] = 0.07*1e-2;
244   ppe[ipoint++] = 0.06*1e-2;
245   ppe[ipoint++] = 0.05*1e-2;
246   ppe[ipoint++] = 0.05*1e-2;
247   ppe[ipoint++] = 0.04*1e-2;
248   ppe[ipoint++] = 0.04*1e-2;
249   ppe[ipoint++] = 0.40*1e-3;
250   ppe[ipoint++] = 0.37*1e-3;
251   ppe[ipoint++] = 0.38*1e-3;
252
253
254
255   gStar[5] = new TGraphErrors(ipoint,pt,pp,errpt,ppe);
256
257 //   gStar[0]->SetMarkerStyle(kOpenTriangleUp);
258 //   gStar[1]->SetMarkerStyle(kFullTriangleUp);
259 //   gStar[2]->SetMarkerStyle(kOpenCircle);
260 //   gStar[3]->SetMarkerStyle(kFullCircle);
261 //   gStar[4]->SetMarkerStyle(kOpenSquare);
262 //   gStar[5]->SetMarkerStyle(kFullSquare);
263   gStar[0]->SetMarkerStyle(kOpenStar);
264   gStar[1]->SetMarkerStyle(kFullStar);
265   gStar[2]->SetMarkerStyle(kOpenStar);
266   gStar[3]->SetMarkerStyle(kFullStar);
267   gStar[4]->SetMarkerStyle(kOpenStar);
268   gStar[5]->SetMarkerStyle(kFullStar);
269
270   gStar[0]->SetMarkerColor(kBlack);
271   gStar[1]->SetMarkerColor(kBlack);
272   gStar[2]->SetMarkerColor(kRed);
273   gStar[3]->SetMarkerColor(kRed);
274   gStar[4]->SetMarkerColor(kBlue);
275   gStar[5]->SetMarkerColor(kBlue);
276
277   gStar[0]->SetLineColor (kBlack);
278   gStar[1]->SetLineColor (kBlack);
279   gStar[2]->SetLineColor (kRed);
280   gStar[3]->SetLineColor (kRed);
281   gStar[4]->SetLineColor (kBlue);
282   gStar[5]->SetLineColor (kBlue);
283
284   // Scale all for 2*pi and pt in order to go to dNdpt
285   for(Int_t istar = 0; istar < 6; istar++){
286     Int_t npoint = gStar[istar]->GetN();
287     for(Int_t ipoint = 0; ipoint < npoint; ipoint++){
288       Double_t ptbin = gStar[istar]->GetX()[ipoint];
289       gStar[istar]->SetPoint(ipoint,
290                              gStar[istar]->GetX()[ipoint],
291                              gStar[istar]->GetY()[ipoint]*2*TMath::Pi()*ptbin);
292       gStar[istar]->SetPointError(ipoint,
293                                   gStar[istar]->GetEX()[ipoint],
294                                   gStar[istar]->GetEY()[ipoint]*2*TMath::Pi()*ptbin);
295     }
296     
297   }
298   
299
300
301 //   gStar[0]->Draw("AP");
302 //   gStar[1]->Draw("P");
303 //   gStar[2]->Draw("P");
304 //   gStar[3]->Draw("P");
305 //   gStar[4]->Draw("P");
306 //   gStar[5]->Draw("P");
307
308
309   return gStar;
310
311
312   
313   // table from PRC79  79 (2009) 34909
314 //         pi-                   pi+               k-                 K+               pbar              p
315 //  0.225  2.02+-0.06            2.07+-0.06        (1.43+-0.11)*10-1 (1.52+-0.11)*10-1
316 //  0.275  1.52+-0.03            1.54+-0.03        (1.26+-0.05)*10-1 (1.30+-0.05)*10-1
317 //  0.325  1.13+-0.02            1.14+-0.02        (1.08+-0.02)*10-1 (1.08+-0.02)*10-1
318 //  0.375  (8.44+-0.09)*10-1     (8.57+-0.09)*10-1 (8.77+-0.24)*10-2 (9.16+-0.25)*10-2 (5.54+-0.13)*10-2
319 //  0.425  (6.35+-0.07)*10-1     (6.38+-0.07)*10-1 (7.34+-0.26)*10-2 (7.47+-0.27)*10-2 (4.81+-0.11)*10-2
320 //  0.475  (4.69+-0.05)*10-1     (4.76+-0.05)*10-1 (6.17+-0.63)*10-2 (6.26+-0.64)*10-2 (4.25+-0.10)*10-2 (5.07+-0.25)*10-2
321 //  0.525  (3.54+-0.04)*10-1     (3.59+-0.04)*10-1 (4.87+-0.50)*10-2 (5.26+-0.54)*10-2 (3.77+-0.09)*10-2 (4.69+-0.19)*10-2
322 //  0.575  (2.67+-0.03)*10-1     (2.73+-0.03)*10-1 (4.11+-0.42)*10-2 (4.41+-0.45)*10-2 (3.38+-0.08)*10-2 (4.08+-0.14)*10-2
323 //  0.625  (2.02+-0.04)*10-1     (2.07+-0.04)*10-1 (3.70+-0.28)*10-2 (3.81+-0.28)*10-2 (2.78+-0.07)*10-2 (3.42+-0.10)*10-2
324 //  0.675  (1.53+-0.03)*10-1     (1.55+-0.03)*10-1 (2.86+-0.22)*10-2 (3.06+-0.24)*10-2 (2.49+-0.06)*10-2 (2.87+-0.07)*10-2
325 //  0.725  (1.16+-0.03)*10-1     (1.16+-0.03)*10-1 (2.42+-0.27)*10-2 (2.49+-0.28)*10-2 (2.03+-0.06)*10-2 (2.41+-0.07)*10-2
326 //  0.775  (8.65+-0.30)*10-2     (8.95+-0.31)*10-2                                     (1.75+-0.06)*10-2 (2.13+-0.06)*10-2
327 //  0.825                                                                              (1.52+-0.05)*10-2 (1.82+-0.05)*10-2
328 //  0.875                                                                              (1.27+-0.04)*10-2 (1.54+-0.05)*10-2
329 //  0.925                                                                              (1.05+-0.04)*10-2 (1.31+-0.04)*10-2
330 //  0.975                                                                              (8.95+-0.34)*10-3 (1.11+-0.04)*10-2
331 //  1.025                                                                              (7.36+-0.34)*10-3 (9.78+-0.40)*10-3
332 //  1.075                                                                              (6.59+-0.32)*10-3 (8.56+-0.37)*10-3
333 //  1.125                                                                              (5.21+-0.29)*10-3 (7.38+-0.38)*10-3
334
335
336   }