ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/GlobalFunctions.h
Revision: 1.9
Committed: Wed Feb 9 08:57:36 2011 UTC (14 years, 3 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.8: +378 -41 lines
Log Message: 
new RA1 curves

File Contents

# User Rev Content
1 auterman 1.1 #include "SusyScan.h"
2     #include "GeneratorMasses.h"
3    
4 auterman 1.5 #include "TGraph.h"
5     #include "TGraphErrors.h"
6     #include "TSpline.h"
7    
8 auterman 1.1 #include <cmath>
9 auterman 1.8 #include <iostream>
10 auterman 1.1
11    
12 auterman 1.3 double Luminosity = 36.3; //[pb^-1]
13 auterman 1.1 double Mzero(const SusyScan* p){ return p->Mzero; }
14     double Mhalf(const SusyScan* p){ return p->Mhalf; }
15     double MGluino(const SusyScan* p){ return p->MGL; }
16     double MSquarkL(const SusyScan* p){ return p->MUL; }
17     double MSquarkR(const SusyScan* p){ return p->MUR; }
18     double MChi1(const SusyScan* p){ return p->MZ1; }
19     double MChi2(const SusyScan* p){ return p->MZ2; }
20     double MChi3(const SusyScan* p){ return p->MZ3; }
21     double MChi4(const SusyScan* p){ return p->MZ4; }
22     double MCha1(const SusyScan* p){ return p->MW1; }
23     double MCha2(const SusyScan* p){ return p->MW2; }
24     double Xsection(const SusyScan* p){ return p->Xsection; }
25     double ExpXsecLimit(const SusyScan* p){ return p->ExpXsecLimit; }
26     double ObsXsecLimit(const SusyScan* p){ return p->ObsXsecLimit; }
27 auterman 1.5 double Signal(const SusyScan* p){ return p->signal; }
28     double SignalUncertainty(const SusyScan* p){ return p->signal_uncertainty; }
29     double SignalRelUncertainty(const SusyScan* p){ return p->signal_uncertainty/p->signal; }
30 auterman 1.4 double ExpExclusion(const SusyScan* p){ return (ExpXsecLimit(p)<Xsection(p)&&ExpXsecLimit(p)>0.01?1:0.01); }
31     double ObsExclusion(const SusyScan* p){ return (ObsXsecLimit(p)<Xsection(p)&&ObsXsecLimit(p)>0.01?1:0.01); }
32 auterman 1.6 double ExpExclCL(const SusyScan* p){ return (p->CLs_b_xsec<=0.05 ? 1:0.01); }
33 auterman 1.7 double ExpExclCLm2sigma(const SusyScan* p){ return (p->CLs_b_n2_xsec<=0.05 ? 1:0.01); }
34     double ExpExclCLm1sigma(const SusyScan* p){ return (p->CLs_b_n1_xsec<=0.05 ? 1:0.01); }
35     double ExpExclCLp1sigma(const SusyScan* p){ return (p->CLs_b_p1_xsec<=0.05 ? 1:0.01); }
36     double ExpExclCLp2sigma(const SusyScan* p){ return (p->CLs_b_p2_xsec<=0.05 ? 1:0.01); }
37 auterman 1.6 double ObsExclCL(const SusyScan* p){ return (p->CLs_xsec<=0.05 ? 1:0.01); }
38 auterman 1.1 double SoverSqrtB(const SusyScan* p){ return p->signal/(sqrt(p->background)+p->background_uncertainty+p->signal_uncertainty); }
39     double XsecOverObserved(const SusyScan* p){ return (ObsXsecLimit(p)==0 ? 9999. : Xsection(p)/ObsXsecLimit(p)); }
40 auterman 1.3 double XsecOverExpected(const SusyScan* p){ return (ExpXsecLimit(p)==0 ? 9999. : Xsection(p)/ExpXsecLimit(p)); }
41 auterman 1.1 double SignalAcceptance(const SusyScan* p){ return p->signal / (Luminosity*Xsection(p)); }
42 auterman 1.5 double ExpNSignLimit(const SusyScan* p){ return p->ExpNsigLimit; }
43     double ObsNSignLimit(const SusyScan* p){ return p->ObsNsigLimit; }
44 auterman 1.3 double PLExpNSignLimit(const SusyScan* p){ return p->PLExpNsigLimit; }
45 auterman 1.5 double PLObsNSignLimit(const SusyScan* p){ return p->PLObsNsigLimit; }
46 auterman 1.3 double PLExpXsecLimit(const SusyScan* p){ return p->PLExpXsecLimit; }
47     double PLObsXsecLimit(const SusyScan* p){ return p->PLObsXsecLimit; }
48 auterman 1.4 double PLExpExclusion(const SusyScan* p){ return (PLExpXsecLimit(p)<Xsection(p)&&PLExpXsecLimit(p)>0.01?1:0.01); }
49     double PLObsExclusion(const SusyScan* p){ return (PLObsXsecLimit(p)<Xsection(p)&&PLObsXsecLimit(p)>0.01?1:0.01); }
50 auterman 1.5 double FCExpNSignLimit(const SusyScan* p){ return p->FCExpNsigLimit; }
51     double FCObsNSignLimit(const SusyScan* p){ return p->FCObsNsigLimit; }
52     double FCExpXsecLimit(const SusyScan* p){ return p->FCExpXsecLimit; }
53     double FCObsXsecLimit(const SusyScan* p){ return p->FCObsXsecLimit; }
54     double FCExpExclusion(const SusyScan* p){ return (FCExpXsecLimit(p)<Xsection(p)&&FCExpXsecLimit(p)>0.01?1:0.01); }
55     double FCObsExclusion(const SusyScan* p){ return (FCObsXsecLimit(p)<Xsection(p)&&FCObsXsecLimit(p)>0.01?1:0.01); }
56     double MCMCExpNSignLimit(const SusyScan* p){ return p->MCMCExpNsigLimit; }
57     double MCMCObsNSignLimit(const SusyScan* p){ return p->MCMCObsNsigLimit; }
58     double MCMCExpXsecLimit(const SusyScan* p){ return p->MCMCExpXsecLimit; }
59     double MCMCObsXsecLimit(const SusyScan* p){ return p->MCMCObsXsecLimit; }
60     double MCMCExpExclusion(const SusyScan* p){ return (MCMCExpXsecLimit(p)<Xsection(p)&&MCMCExpXsecLimit(p)>0.01?1:0.01); }
61     double MCMCObsExclusion(const SusyScan* p){ return (MCMCObsXsecLimit(p)<Xsection(p)&&MCMCObsXsecLimit(p)>0.01?1:0.01); }
62 auterman 1.8 double SignalContamination(const SusyScan* p){return p->signal_contamination; }
63 auterman 1.1
64     double Mzero(const GeneratorMasses* p){ return p->Mzero; }
65     double Mhalf(const GeneratorMasses* p){ return p->Mhalf; }
66     double MGluino(const GeneratorMasses* p){ return p->MGL; }
67     double MSquarkL(const GeneratorMasses* p){ return p->MUL; }
68     double MSquarkR(const GeneratorMasses* p){ return p->MUR; }
69     double MChi1(const GeneratorMasses* p){ return p->MZ1; }
70     double MChi2(const GeneratorMasses* p){ return p->MZ2; }
71     double MChi3(const GeneratorMasses* p){ return p->MZ3; }
72     double MChi4(const GeneratorMasses* p){ return p->MZ4; }
73     double MCha1(const GeneratorMasses* p){ return p->MW1; }
74     double MCha2(const GeneratorMasses* p){ return p->MW2; }
75 auterman 1.5
76    
77     //RA1 Contours
78     TGraph* set_sneutrino_d0_1(){
79     double sn_m0[9]={0,0,55,90,100,110,100,55,0};
80     double sn_m12[9]={0,140,220,240,240,230,210,150,0};
81    
82     TGraph* sn_d0_gr = new TGraph(9,sn_m0,sn_m12);
83    
84     sn_d0_gr->SetFillColor(kGreen+3);
85     sn_d0_gr->SetFillStyle(1001);
86    
87    
88     return sn_d0_gr;
89     }
90    
91     TGraph* set_sneutrino_d0_2(){
92     double sn_m0[5]={0,50,105,190,0};
93     double sn_m12[5]={0,140,205,140,0};
94    
95     TGraph* sn_d0_gr_2 = new TGraph(5,sn_m0,sn_m12);
96    
97     sn_d0_gr_2->SetFillColor(kGreen+3);
98     sn_d0_gr_2->SetFillStyle(1001);
99    
100    
101     return sn_d0_gr_2;
102     }
103    
104    
105    
106     TGraph* set_lep_ch(){
107    
108     double ch_m0[5];
109     double ch_m12[5];
110    
111    
112     ch_m0[0] = 0;
113     ch_m0[1] = 0;
114     ch_m0[2] = 600;
115     ch_m0[3] = 1000;
116     ch_m0[4] = 1000;
117    
118     ch_m12[0] = 0;
119     ch_m12[1] = 130;
120     ch_m12[2] = 120;
121     ch_m12[3] = 113;
122     ch_m12[4] = 0;
123    
124    
125     TGraph* ch_gr = new TGraph(5,ch_m0,ch_m12);
126    
127     ch_gr->SetFillColor(3);
128     ch_gr->SetFillStyle(1001);
129    
130     return ch_gr;
131    
132     }
133    
134    
135     TGraph* set_lep_sl(){
136    
137    
138    
139     double sl_m0[] = {0,0,30,50,60,75,80,90,100};
140     double sl_m12[] = {0,245,240,220,200,150,100,50,0};
141    
142     TGraph* lep_sl = new TGraph(9,sl_m0,sl_m12);
143    
144     lep_sl->SetFillColor(5);
145     lep_sl->SetFillStyle(1001);
146    
147     return lep_sl;
148     }
149    
150     TGraph* set_tev_sg_cdf(){
151    
152     double sg_m0[] = {0,50,100,150,200,250,300,350,400,450,500,550,600,600};
153     double sg_m12[] = {0,170,160,155,150,122,116,112,110,106,105,100,98,0};
154     TGraph* sg_gr = new TGraph(14,sg_m0,sg_m12);
155    
156     sg_gr->SetFillColor(2);
157     sg_gr->SetFillStyle(1001);
158    
159     return sg_gr;
160    
161     }
162    
163     TGraph* set_tev_sg_d0(){
164     double sgd_m0[] = {0,50,100,150,200,250,300,350,400,450,500,550,600,600};
165     double sgd_m12[] = {0,173,170,168,160,150,140,130,125,120,120,120,120,0};
166     TGraph* sgd_gr = new TGraph(14,sgd_m0,sgd_m12);
167    
168     sgd_gr->SetFillColor(41);
169     sgd_gr->SetFillStyle(1001);
170    
171     return sgd_gr;
172    
173     }
174    
175     TGraph* set_tev_tlp_cdf(){
176     double tlp1_m0[] = {0,20,40,60,70,80,90,80,70,60};
177     double tlp1_m12[] = {170,185,200,215,220,215,210,190,175,160};
178     TGraph* tlp1_gr = new TGraph(10,tlp1_m0,tlp1_m12);
179    
180     tlp1_gr->SetFillColor(4);
181     tlp1_gr->SetFillStyle(1001);
182    
183     return tlp1_gr;
184     }
185    
186     TGraph* set_tev_tlp_d0(){
187     double tlp2_m0[] = {70,80,90,100,105,110,120,130,140};
188     double tlp2_m12[] = {160,172,184,196,205,195,185,173,160};
189     TGraph* tlp2_gr = new TGraph(9,tlp2_m0,tlp2_m12);
190    
191     tlp2_gr->SetFillColor(4);
192     tlp2_gr->SetFillStyle(1001);
193    
194     return tlp2_gr;
195    
196     }
197    
198     TGraph* set_tev_stau(){
199     double st_m0[] = {0,30,200,0,0};
200     double st_m12[] = {230,240,1000,1000,230};
201     TGraph* st_gr = new TGraph(5,st_m0,st_m12);
202    
203     st_gr->SetFillColor(40);
204     st_gr->SetFillStyle(1001);
205    
206    
207     return st_gr;
208    
209     }
210    
211    
212 auterman 1.9 TGraphErrors* getLO_tanBeta3(){
213 auterman 1.5
214    
215    
216     Int_t nl = 9;
217     Double_t xl[9];
218     Double_t yl[9];
219     Double_t exl[9];
220     Double_t eyl[9];
221    
222     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
223    
224 auterman 1.9
225     xl[0] = 0;
226 auterman 1.5 yl[0] = 265;
227     xl[1] = 100;
228     yl[1] = 258;
229     xl[2] = 200;
230     yl[2] = 250;
231     xl[3] = 250;
232     yl[3] = 240;
233     xl[4] = 300;
234     yl[4] = 210;
235     xl[5] = 340;
236     yl[5] = 177;
237     xl[6] = 400;
238     yl[6] = 140;
239 auterman 1.9 xl[7] = 450;
240 auterman 1.5 yl[7] = 120;
241 auterman 1.9 xl[8] = 520;
242 auterman 1.5 yl[8] =100;
243    
244 auterman 1.9
245     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
246     gr1->SetMarkerColor(kGreen+2);
247     gr1->SetMarkerStyle(21);
248    
249    
250     //gr1->Draw("LP");
251    
252     TSpline3 *s = new TSpline3("grs",gr1);
253     s->SetLineColor(kGreen+2);
254     s->SetLineStyle(4);
255     s->SetLineWidth(3);
256    
257    
258     return gr1;
259    
260    
261    
262     }
263    
264     TGraphErrors* getLO_tanBeta10(){
265    
266    
267    
268     Int_t nl = 10;
269     Double_t xl[10];
270     Double_t yl[10];
271     Double_t exl[10];
272     Double_t eyl[10];
273    
274     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
275    
276     xl[0] = 0;
277     yl[0] = 270;
278     xl[1] = 100;
279     yl[1] = 260;
280     xl[2] = 200;
281     yl[2] = 250;
282     xl[3] = 250;
283     yl[3] = 240;
284     xl[4] = 300;
285     yl[4] = 210;
286     xl[5] = 350;
287     yl[5] = 174;
288     xl[6] = 400;
289     yl[6] = 147;
290     xl[7] = 450;
291     yl[7] = 127;
292     xl[8] = 500;
293     yl[8] =115;
294     xl[9] = 520;
295     yl[9] = 112;
296    
297    
298    
299    
300    
301    
302     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
303     gr1->SetMarkerColor(kGreen+2);
304     gr1->SetMarkerStyle(21);
305    
306    
307     //gr1->Draw("LP");
308    
309     TSpline3 *s = new TSpline3("grs",gr1);
310     s->SetLineColor(kGreen+2);
311     s->SetLineStyle(4);
312     s->SetLineWidth(3);
313    
314    
315     return gr1;
316    
317    
318    
319     }
320    
321     TGraphErrors* getLO_tanBeta50(){
322    
323 auterman 1.5
324 auterman 1.9
325     Int_t nl = 10;
326     Double_t xl[10];
327     Double_t yl[10];
328     Double_t exl[10];
329     Double_t eyl[10];
330    
331     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
332    
333    
334     xl[0] = 200;
335     yl[0] = 239;
336     xl[1] = 210;
337     yl[1] = 249;
338     xl[2] = 229;
339     yl[2] = 260;
340     xl[3] = 250;
341     yl[3] = 245;
342     xl[4] = 300;
343     yl[4] = 210;
344     xl[5] = 350;
345     yl[5] = 180;
346     xl[6] = 400;
347     yl[6] = 160;
348     xl[7] = 450;
349     yl[7] = 150;
350     xl[8] = 500;
351     yl[8] =140;
352     xl[9] = 520;
353     yl[9] = 137;
354    
355    
356    
357 auterman 1.5
358    
359    
360     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
361     gr1->SetMarkerColor(kGreen+2);
362     gr1->SetMarkerStyle(21);
363    
364    
365     //gr1->Draw("LP");
366    
367     TSpline3 *s = new TSpline3("grs",gr1);
368     s->SetLineColor(kGreen+2);
369     s->SetLineStyle(4);
370     s->SetLineWidth(3);
371    
372    
373     return gr1;
374    
375    
376    
377     }
378    
379    
380 auterman 1.9
381     TGraphErrors* getExpected_NLO_tanBeta3(){
382 auterman 1.5
383     Int_t nl = 11;
384     Double_t xl[11];
385     Double_t yl[11];
386     Double_t exl[11];
387     Double_t eyl[11];
388    
389 auterman 1.9 xl[0] = 0;
390     yl[0] = 283;
391     xl[1] = 100;
392     yl[1] = 280;
393     xl[2] = 150;
394     yl[2] = 279;
395     xl[3] = 200;
396     yl[3] = 275;
397     xl[4] = 250;
398     yl[4] = 270;
399     xl[5] = 300;
400     yl[5] = 255;
401     xl[6] = 350;
402     yl[6] = 225;
403     xl[7] = 400;
404     yl[7] = 195;
405     xl[8] = 450;
406     yl[8] = 175;
407     xl[9] = 500;
408     yl[9] = 155;
409     xl[10] = 550;
410     yl[10] = 150;
411    
412    
413    
414     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
415     gr1->SetMarkerColor(kBlue);
416     gr1->SetMarkerStyle(21);
417    
418    
419     //gr1->Draw("LP");
420    
421     TSpline3 *s = new TSpline3("grs",gr1);
422     s->SetLineColor(kBlue);
423     s->SetLineStyle(2);
424     s->SetLineWidth(3);
425 auterman 1.5
426    
427 auterman 1.9 return gr1;
428    
429 auterman 1.5
430 auterman 1.9
431    
432    
433     }
434    
435     TGraphErrors* getExpected_NLO_tanBeta10(){
436    
437     Int_t nl = 11;
438     Double_t xl[11];
439     Double_t yl[11];
440     Double_t exl[11];
441     Double_t eyl[11];
442    
443     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
444    
445     xl[0] = 0;
446 auterman 1.5 yl[0] = 283;
447     xl[1] = 100;
448     yl[1] = 280;
449     xl[2] = 150;
450     yl[2] = 279;
451     xl[3] = 200;
452     yl[3] = 275;
453     xl[4] = 250;
454     yl[4] = 270;
455     xl[5] = 300;
456     yl[5] = 255;
457     xl[6] = 350;
458     yl[6] = 225;
459     xl[7] = 400;
460     yl[7] = 195;
461     xl[8] = 450;
462     yl[8] = 175;
463     xl[9] = 500;
464 auterman 1.9 yl[9] = 165;
465 auterman 1.5 xl[10] = 550;
466 auterman 1.9 yl[10] = 150;
467    
468    
469    
470    
471    
472     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
473     gr1->SetMarkerColor(kBlue);
474     gr1->SetMarkerStyle(21);
475    
476    
477     //gr1->Draw("LP");
478    
479     TSpline3 *s = new TSpline3("grs",gr1);
480     s->SetLineColor(kBlue);
481     s->SetLineStyle(2);
482     s->SetLineWidth(3);
483    
484    
485     return gr1;
486    
487    
488    
489    
490    
491     }
492    
493    
494     TGraphErrors* getExpected_NLO_tanBeta50(){
495    
496     Int_t nl = 10;
497     Double_t xl[10];
498     Double_t yl[10];
499     Double_t exl[10];
500     Double_t eyl[10];
501 auterman 1.5
502 auterman 1.9 // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
503    
504     xl[0] = 200;
505     yl[0] = 287;
506     xl[1] = 220;
507     yl[1] = 287;
508     xl[2] = 245;
509     yl[2] = 287;
510     xl[3] = 270;
511     yl[3] = 265;
512     xl[4] = 300;
513     yl[4] = 245;
514     xl[5] = 350;
515     yl[5] = 222;
516     xl[6] = 400;
517     yl[6] = 197;
518     xl[7] = 450;
519     yl[7] = 180;
520     xl[8] = 500;
521     yl[8] = 168;
522     xl[9] = 550;
523     yl[9] = 145;
524    
525    
526    
527 auterman 1.5
528    
529     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
530 auterman 1.9 gr1->SetMarkerColor(kBlue);
531     gr1->SetMarkerStyle(21);
532    
533    
534     //gr1->Draw("LP");
535    
536     TSpline3 *s = new TSpline3("grs",gr1);
537     s->SetLineColor(kBlue);
538     s->SetLineStyle(2);
539     s->SetLineWidth(3);
540    
541 auterman 1.5
542     return gr1;
543    
544    
545    
546    
547    
548     }
549    
550    
551 auterman 1.9 TGraphErrors* getObserved_NLO_tanBeta3(){
552 auterman 1.5
553 auterman 1.9 Int_t nl = 11;
554 auterman 1.5 Double_t xl[11];
555     Double_t yl[11];
556     Double_t exl[11];
557     Double_t eyl[11];
558    
559     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
560    
561    
562 auterman 1.9
563     xl[0] = 0;
564     yl[0] = 274;
565 auterman 1.5 xl[1] = 100;
566     yl[1] = 270;
567     xl[2] = 150;
568     yl[2] = 268;
569     xl[3] = 200;
570     yl[3] = 265;
571     xl[4] = 250;
572     yl[4] = 255;
573     xl[5] = 300;
574 auterman 1.9 yl[5] = 230;
575 auterman 1.5 xl[6] = 350;
576 auterman 1.9 yl[6] = 195;
577 auterman 1.5 xl[7] = 400;
578 auterman 1.9 yl[7] = 160;
579 auterman 1.5 xl[8] = 450;
580 auterman 1.9 yl[8] = 140;
581     xl[9] = 480;
582     yl[9] = 130;
583     xl[10] = 530;
584 auterman 1.5 yl[10] = 120;
585 auterman 1.9
586    
587    
588    
589     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
590     gr1->SetMarkerColor(kBlue);
591     gr1->SetMarkerStyle(21);
592    
593    
594     //gr1->Draw("LP");
595    
596     TSpline3 *s = new TSpline3("grs",gr1);
597     s->SetLineColor(kRed);
598     // s->SetLineStyle(2);
599     s->SetLineWidth(3);
600    
601    
602     return gr1;
603    
604    
605    
606    
607    
608     }
609    
610    
611    
612     TGraphErrors* getObserved_NLO_tanBeta10(){
613    
614     Int_t nl = 11;
615     Double_t xl[11];
616     Double_t yl[11];
617     Double_t exl[11];
618     Double_t eyl[11];
619    
620     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
621    
622     xl[0] = 0;
623     yl[0] = 278;
624     xl[1] = 100;
625     yl[1] = 270;
626     xl[2] = 150;
627     yl[2] = 267;
628     xl[3] = 200;
629     yl[3] = 262;
630     xl[4] = 250;
631     yl[4] = 250;
632     xl[5] = 300;
633     yl[5] = 225;
634     xl[6] = 350;
635     yl[6] = 192;
636     xl[7] = 400;
637     yl[7] = 163;
638     xl[8] = 450;
639     yl[8] = 148;
640     xl[9] = 500;
641     yl[9] = 140;
642     xl[10] = 520;
643     yl[10] = 137;
644    
645    
646    
647    
648    
649     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
650     gr1->SetMarkerColor(kBlue);
651     gr1->SetMarkerStyle(21);
652    
653    
654     //gr1->Draw("LP");
655    
656     TSpline3 *s = new TSpline3("grs",gr1);
657     s->SetLineColor(kRed);
658     // s->SetLineStyle(2);
659     s->SetLineWidth(3);
660    
661    
662     return gr1;
663    
664    
665    
666    
667    
668     }
669    
670    
671     TGraphErrors* getObserved_NLO_tanBeta50(){
672    
673     Int_t nl = 10;
674     Double_t xl[10];
675     Double_t yl[10];
676     Double_t exl[10];
677     Double_t eyl[10];
678    
679     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
680    
681    
682     xl[0] = 200;
683     yl[0] = 243;
684     xl[1] = 220;
685     yl[1] = 264;
686     xl[2] = 235;
687     yl[2] = 278;
688     xl[3] = 250;
689     yl[3] = 267;
690     xl[4] = 300;
691     yl[4] = 230;
692     xl[5] = 350;
693     yl[5] = 205;
694     xl[6] = 400;
695     yl[6] = 184;
696     xl[7] = 450;
697     yl[7] = 168;
698     xl[8] = 500;
699     yl[8] = 156;
700     xl[9] = 520;
701     yl[9] = 148;
702 auterman 1.5
703 auterman 1.9
704    
705 auterman 1.5
706    
707     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
708 auterman 1.9 gr1->SetMarkerColor(kBlue);
709     gr1->SetMarkerStyle(21);
710 auterman 1.5
711    
712     //gr1->Draw("LP");
713    
714     TSpline3 *s = new TSpline3("grs",gr1);
715     s->SetLineColor(kRed);
716     // s->SetLineStyle(2);
717     s->SetLineWidth(3);
718    
719    
720     return gr1;
721    
722    
723    
724    
725    
726     }
727