ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/ExclusionPlot.h
Revision: 1.1
Committed: Mon May 23 16:48:15 2011 UTC (13 years, 11 months ago) by auterman
Content type: text/plain
Branch: MAIN
CVS Tags: JHEP2010, HEAD
Log Message: 
bla

File Contents

# User Rev Content
1 auterman 1.1 #ifndef MAKE_EXCLUSIONPLOT_HH
2     #define MAKE_EXCLUSIONPLOT_HH
3    
4     #include "TH1F.h"
5     #include "TH2F.h"
6     #include "TGraph.h"
7     #include "TF1.h"
8     #include "TLegend.h"
9     #include "TLatex.h"
10     #include "TFile.h"
11     #include "TSpline.h"
12     #include "TGraphErrors.h"
13    
14    
15     #include <vector>
16    
17    
18    
19    
20    
21     void ExclusionPlot();
22     void CommandMSUGRA(TString plotName,Int_t tanBeta);
23     void setPlottingStyle(TH1F& hsig);
24    
25     TGraph* set_lep_ch(Int_t tanBeta);
26     TGraph* set_lep_ch_tanBeta3();
27     TGraph* set_lep_ch_tanBeta10();
28     TGraph* set_lep_ch_tanBeta50();
29     TGraph* set_lep_sl(Int_t tanBeta);//slepton curve
30     TGraph* set_tev_sg_cdf(Int_t tanBeta);//squark gluino cdf
31     TGraph* set_tev_sg_d0(Int_t tanBeta);//squark gluino d0
32     TGraph* set_tev_stau(Int_t tanBeta);//stau
33     TGraph* set_sneutrino_d0_1(Int_t tanBeta);
34     TGraph* set_sneutrino_d0_2(Int_t tanBeta);
35     TF1* constant_squark(int tanBeta,int i);
36     TF1* constant_gluino(int tanBeta,int i);
37     TLatex* constant_squark_text(Int_t it,TF1& lnsq,Int_t tanBeta);
38     TLatex* constant_gluino_text(Int_t it,TF1& lngl);
39     TLegend* makeStauLegend(Double_t txtsz,Int_t tanBeta_);
40     TLegend* makeExpLegend(TGraph& sg_gr, TGraph& sgd_gr,TGraph& ch_gr,TGraph& sl_gr,TGraph&,Double_t txtsz,Int_t tanbeta);
41    
42    
43    
44    
45    
46    
47    
48     TGraphErrors* getLO_tanBeta3(){
49    
50    
51    
52     Int_t nl = 9;
53     Double_t xl[9];
54     Double_t yl[9];
55     Double_t exl[9];
56     Double_t eyl[9];
57    
58     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
59    
60    
61     xl[0] = 0;
62     yl[0] = 265;
63     xl[1] = 100;
64     yl[1] = 258;
65     xl[2] = 200;
66     yl[2] = 250;
67     xl[3] = 250;
68     yl[3] = 240;
69     xl[4] = 300;
70     yl[4] = 210;
71     xl[5] = 340;
72     yl[5] = 177;
73     xl[6] = 400;
74     yl[6] = 140;
75     xl[7] = 450;
76     yl[7] = 120;
77     xl[8] = 520;
78     yl[8] =100;
79    
80    
81     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
82     gr1->SetMarkerColor(kGreen+2);
83     gr1->SetMarkerStyle(21);
84    
85    
86     //gr1->Draw("LP");
87    
88     TSpline3 *s = new TSpline3("grs",gr1);
89     s->SetLineColor(kGreen+2);
90     s->SetLineStyle(4);
91     s->SetLineWidth(3);
92    
93    
94     return gr1;
95    
96    
97    
98     }
99    
100     TGraphErrors* getLO_tanBeta10(){
101    
102    
103    
104     Int_t nl = 10;
105     Double_t xl[10];
106     Double_t yl[10];
107     Double_t exl[10];
108     Double_t eyl[10];
109    
110     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
111    
112     xl[0] = 0;
113     yl[0] = 270;
114     xl[1] = 100;
115     yl[1] = 260;
116     xl[2] = 200;
117     yl[2] = 250;
118     xl[3] = 250;
119     yl[3] = 240;
120     xl[4] = 300;
121     yl[4] = 210;
122     xl[5] = 350;
123     yl[5] = 174;
124     xl[6] = 400;
125     yl[6] = 147;
126     xl[7] = 450;
127     yl[7] = 127;
128     xl[8] = 500;
129     yl[8] =115;
130     xl[9] = 520;
131     yl[9] = 112;
132    
133    
134    
135    
136    
137    
138     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
139     gr1->SetMarkerColor(kGreen+2);
140     gr1->SetMarkerStyle(21);
141    
142    
143     //gr1->Draw("LP");
144    
145     TSpline3 *s = new TSpline3("grs",gr1);
146     s->SetLineColor(kGreen+2);
147     s->SetLineStyle(4);
148     s->SetLineWidth(3);
149    
150    
151     return gr1;
152    
153    
154    
155     }
156    
157     TGraphErrors* getLO_tanBeta50(){
158    
159    
160    
161     Int_t nl = 10;
162     Double_t xl[10];
163     Double_t yl[10];
164     Double_t exl[10];
165     Double_t eyl[10];
166    
167     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
168    
169    
170     xl[0] = 200;
171     yl[0] = 239;
172     xl[1] = 210;
173     yl[1] = 249;
174     xl[2] = 229;
175     yl[2] = 260;
176     xl[3] = 250;
177     yl[3] = 245;
178     xl[4] = 300;
179     yl[4] = 210;
180     xl[5] = 350;
181     yl[5] = 180;
182     xl[6] = 400;
183     yl[6] = 160;
184     xl[7] = 450;
185     yl[7] = 150;
186     xl[8] = 500;
187     yl[8] =140;
188     xl[9] = 520;
189     yl[9] = 137;
190    
191    
192    
193    
194    
195    
196     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
197     gr1->SetMarkerColor(kGreen+2);
198     gr1->SetMarkerStyle(21);
199    
200    
201     //gr1->Draw("LP");
202    
203     TSpline3 *s = new TSpline3("grs",gr1);
204     s->SetLineColor(kGreen+2);
205     s->SetLineStyle(4);
206     s->SetLineWidth(3);
207    
208    
209     return gr1;
210    
211    
212    
213     }
214    
215    
216    
217     TGraphErrors* getExpected_NLO_tanBeta3(){
218    
219     Int_t nl = 11;
220     Double_t xl[11];
221     Double_t yl[11];
222     Double_t exl[11];
223     Double_t eyl[11];
224    
225     xl[0] = 0;
226     yl[0] = 283;
227     xl[1] = 100;
228     yl[1] = 280;
229     xl[2] = 150;
230     yl[2] = 279;
231     xl[3] = 200;
232     yl[3] = 275;
233     xl[4] = 250;
234     yl[4] = 270;
235     xl[5] = 300;
236     yl[5] = 255;
237     xl[6] = 350;
238     yl[6] = 225;
239     xl[7] = 400;
240     yl[7] = 195;
241     xl[8] = 450;
242     yl[8] = 175;
243     xl[9] = 500;
244     yl[9] = 155;
245     xl[10] = 550;
246     yl[10] = 150;
247    
248    
249    
250     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
251     gr1->SetMarkerColor(kBlue);
252     gr1->SetMarkerStyle(21);
253    
254    
255     //gr1->Draw("LP");
256    
257     TSpline3 *s = new TSpline3("grs",gr1);
258     s->SetLineColor(kBlue);
259     s->SetLineStyle(2);
260     s->SetLineWidth(3);
261    
262    
263     return gr1;
264    
265    
266    
267    
268    
269     }
270    
271     TGraphErrors* getExpected_NLO_tanBeta10(){
272    
273     Int_t nl = 11;
274     Double_t xl[11];
275     Double_t yl[11];
276     Double_t exl[11];
277     Double_t eyl[11];
278    
279     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
280    
281     xl[0] = 0;
282     yl[0] = 283;
283     xl[1] = 100;
284     yl[1] = 280;
285     xl[2] = 150;
286     yl[2] = 279;
287     xl[3] = 200;
288     yl[3] = 275;
289     xl[4] = 250;
290     yl[4] = 270;
291     xl[5] = 300;
292     yl[5] = 255;
293     xl[6] = 350;
294     yl[6] = 225;
295     xl[7] = 400;
296     yl[7] = 195;
297     xl[8] = 450;
298     yl[8] = 175;
299     xl[9] = 500;
300     yl[9] = 165;
301     xl[10] = 550;
302     yl[10] = 150;
303    
304    
305    
306    
307    
308     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
309     gr1->SetMarkerColor(kBlue);
310     gr1->SetMarkerStyle(21);
311    
312    
313     //gr1->Draw("LP");
314    
315     TSpline3 *s = new TSpline3("grs",gr1);
316     s->SetLineColor(kBlue);
317     s->SetLineStyle(2);
318     s->SetLineWidth(3);
319    
320    
321     return gr1;
322    
323    
324    
325    
326    
327     }
328    
329    
330     TGraphErrors* getExpected_NLO_tanBeta10_2(){
331    
332     Int_t nl = 11;
333     Double_t xl[11];
334     Double_t yl[11];
335     Double_t exl[11];
336     Double_t eyl[11];
337    
338     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
339    
340     xl[0] = 0;
341     yl[0] = 300;
342     xl[1] = 100;
343     yl[1] = 288;
344     xl[2] = 150;
345     yl[2] = 283;
346     xl[3] = 200;
347     yl[3] = 278;
348     xl[4] = 250;
349     yl[4] = 268;
350     xl[5] = 300;
351     yl[5] = 255;
352     xl[6] = 350;
353     yl[6] = 225;
354     xl[7] = 400;
355     yl[7] = 197;
356     xl[8] = 450;
357     yl[8] = 175;
358     xl[9] = 500;
359     yl[9] = 160;
360     xl[10] = 650;
361     yl[10] = 145;
362    
363    
364    
365    
366     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
367     gr1->SetMarkerColor(kBlue);
368     gr1->SetMarkerStyle(21);
369    
370    
371     //gr1->Draw("LP");
372    
373     TSpline3 *s = new TSpline3("grs",gr1);
374     s->SetLineColor(kBlue);
375     s->SetLineStyle(2);
376     s->SetLineWidth(3);
377    
378    
379     return gr1;
380    
381    
382    
383    
384    
385     }
386    
387    
388     TGraphErrors* getExpected_NLO_tanBeta3_2(){
389    
390     Int_t nl = 11;
391     Double_t xl[11];
392     Double_t yl[11];
393     Double_t exl[11];
394     Double_t eyl[11];
395    
396     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
397    
398     xl[0] = 0;
399     yl[0] = 295;
400     xl[1] = 100;
401     yl[1] = 290;
402     xl[2] = 150;
403     yl[2] = 283;
404     xl[3] = 200;
405     yl[3] = 275;
406     xl[4] = 250;
407     yl[4] = 270;
408     xl[5] = 300;
409     yl[5] = 255;
410     xl[6] = 350;
411     yl[6] = 225;
412     xl[7] = 400;
413     yl[7] = 195;
414     xl[8] = 450;
415     yl[8] = 175;
416     xl[9] = 500;
417     yl[9] = 160;
418     xl[10] = 600;
419     yl[10] = 140;
420    
421    
422    
423    
424     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
425     gr1->SetMarkerColor(kBlue);
426     gr1->SetMarkerStyle(21);
427    
428    
429     //gr1->Draw("LP");
430    
431     TSpline3 *s = new TSpline3("grs",gr1);
432     s->SetLineColor(kBlue);
433     s->SetLineStyle(2);
434     s->SetLineWidth(3);
435    
436    
437     return gr1;
438    
439    
440    
441    
442    
443     }
444    
445    
446     TGraphErrors* getExpected_NLO_tanBeta10_up(){
447    
448     Int_t nl = 11;
449     Double_t xl[11];
450     Double_t yl[11];
451     Double_t exl[11];
452     Double_t eyl[11];
453    
454     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
455    
456     xl[0] = 0;
457     yl[0] = 318;
458     xl[1] = 100;
459     yl[1] = 318;
460     xl[2] = 150;
461     yl[2] = 312;
462     xl[3] = 200;
463     yl[3] = 308;
464     xl[4] = 250;
465     yl[4] = 300;
466     xl[5] = 300;
467     yl[5] = 290;
468     xl[6] = 350;
469     yl[6] = 265;
470     xl[7] = 400;
471     yl[7] = 239;
472     xl[8] = 450;
473     yl[8] = 208;
474     xl[9] = 500;
475     yl[9] = 190;
476     xl[10] = 620;
477     yl[10] = 180;
478    
479    
480    
481    
482     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
483     gr1->SetMarkerColor(kBlue);
484     gr1->SetMarkerStyle(21);
485    
486    
487     //gr1->Draw("LP");
488    
489     TSpline3 *s = new TSpline3("grs",gr1);
490     s->SetLineColor(kBlue);
491     s->SetLineStyle(2);
492     s->SetLineWidth(3);
493    
494    
495     return gr1;
496    
497    
498    
499    
500    
501     }
502    
503    
504    
505     TGraphErrors* getExpected_NLO_tanBeta3_up(){
506    
507     Int_t nl = 11;
508     Double_t xl[11];
509     Double_t yl[11];
510     Double_t exl[11];
511     Double_t eyl[11];
512    
513     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
514    
515     xl[0] = 0;
516     yl[0] = 320;
517     xl[1] = 100;
518     yl[1] = 313;
519     xl[2] = 150;
520     yl[2] = 310;
521     xl[3] = 200;
522     yl[3] = 307;
523     xl[4] = 250;
524     yl[4] = 300;
525     xl[5] = 300;
526     yl[5] = 290;
527     xl[6] = 350;
528     yl[6] = 258;
529     xl[7] = 400;
530     yl[7] = 230;
531     xl[8] = 450;
532     yl[8] = 200;
533     xl[9] = 500;
534     yl[9] = 188;
535     xl[10] = 650;
536     yl[10] = 175;
537    
538    
539    
540    
541     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
542     gr1->SetMarkerColor(kBlue);
543     gr1->SetMarkerStyle(21);
544    
545    
546     //gr1->Draw("LP");
547    
548     TSpline3 *s = new TSpline3("grs",gr1);
549     s->SetLineColor(kBlue);
550     s->SetLineStyle(2);
551     s->SetLineWidth(3);
552    
553    
554     return gr1;
555    
556    
557    
558    
559    
560     }
561    
562    
563    
564     TGraphErrors* getExpected_NLO_tanBeta10_low(){
565    
566     Int_t nl = 11;
567     Double_t xl[11];
568     Double_t yl[11];
569     Double_t exl[11];
570     Double_t eyl[11];
571    
572     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
573    
574     xl[0] = 0;
575     yl[0] = 275;
576     xl[1] = 100;
577     yl[1] = 268;
578     xl[2] = 150;
579     yl[2] = 262;
580     xl[3] = 200;
581     yl[3] = 255;
582     xl[4] = 250;
583     yl[4] = 250;
584     xl[5] = 300;
585     yl[5] = 230;
586     xl[6] = 350;
587     yl[6] = 190;
588     xl[7] = 400;
589     yl[7] = 160;
590     xl[8] = 450;
591     yl[8] = 148;
592     xl[9] = 500;
593     yl[9] = 140;
594     xl[10] = 650;
595     yl[10] = 120;
596    
597    
598    
599    
600     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
601     gr1->SetMarkerColor(kBlue);
602     gr1->SetMarkerStyle(21);
603    
604    
605     //gr1->Draw("LP");
606    
607     TSpline3 *s = new TSpline3("grs",gr1);
608     s->SetLineColor(kBlue);
609     s->SetLineStyle(2);
610     s->SetLineWidth(3);
611    
612    
613     return gr1;
614    
615    
616    
617    
618    
619     }
620    
621     TGraphErrors* getExpected_NLO_tanBeta3_low(){
622    
623     Int_t nl = 11;
624     Double_t xl[11];
625     Double_t yl[11];
626     Double_t exl[11];
627     Double_t eyl[11];
628    
629     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
630    
631     xl[0] = 0;
632     yl[0] = 265;
633     xl[1] = 100;
634     yl[1] = 260;
635     xl[2] = 150;
636     yl[2] = 258;
637     xl[3] = 200;
638     yl[3] = 255;
639     xl[4] = 250;
640     yl[4] = 250;
641     xl[5] = 300;
642     yl[5] = 230;
643     xl[6] = 350;
644     yl[6] = 185;
645     xl[7] = 400;
646     yl[7] = 155;
647     xl[8] = 450;
648     yl[8] = 138;
649     xl[9] = 480;
650     yl[9] = 128;
651     xl[10] = 530;
652     yl[10] = 118;
653    
654     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
655     gr1->SetMarkerColor(kBlue);
656     gr1->SetMarkerStyle(21);
657    
658    
659     //gr1->Draw("LP");
660    
661     TSpline3 *s = new TSpline3("grs",gr1);
662     s->SetLineColor(kBlue);
663     s->SetLineStyle(2);
664     s->SetLineWidth(3);
665    
666     return gr1;
667     }
668    
669    
670    
671    
672     TGraphErrors* getExpected_NLO_tanBeta50(){
673    
674     Int_t nl = 10;
675     Double_t xl[10];
676     Double_t yl[10];
677     Double_t exl[10];
678     Double_t eyl[10];
679    
680     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
681    
682     xl[0] = 200;
683     yl[0] = 287;
684     xl[1] = 220;
685     yl[1] = 287;
686     xl[2] = 245;
687     yl[2] = 287;
688     xl[3] = 270;
689     yl[3] = 265;
690     xl[4] = 300;
691     yl[4] = 245;
692     xl[5] = 350;
693     yl[5] = 222;
694     xl[6] = 400;
695     yl[6] = 197;
696     xl[7] = 450;
697     yl[7] = 180;
698     xl[8] = 500;
699     yl[8] = 168;
700     xl[9] = 550;
701     yl[9] = 145;
702    
703    
704    
705    
706    
707     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
708     gr1->SetMarkerColor(kBlue);
709     gr1->SetMarkerStyle(21);
710    
711    
712     //gr1->Draw("LP");
713    
714     TSpline3 *s = new TSpline3("grs",gr1);
715     s->SetLineColor(kBlue);
716     s->SetLineStyle(2);
717     s->SetLineWidth(3);
718    
719    
720     return gr1;
721    
722    
723    
724    
725    
726     }
727    
728    
729    
730    
731     TGraphErrors* getExpected_NLO_tanBeta50_2(){
732    
733     Int_t nl = 11;
734     Double_t xl[11];
735     Double_t yl[11];
736     Double_t exl[11];
737     Double_t eyl[11];
738    
739     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
740    
741     xl[0] = 200;
742     yl[0] = 287;
743     xl[1] = 220;
744     yl[1] = 287;
745     xl[2] = 245;
746     yl[2] = 287;
747     xl[3] = 270;
748     yl[3] = 270;
749     xl[4] = 300;
750     yl[4] = 255;
751     xl[5] = 350;
752     yl[5] = 225;
753     xl[6] = 400;
754     yl[6] = 203;
755     xl[7] = 450;
756     yl[7] = 180;
757     xl[8] = 500;
758     yl[8] = 158;
759     xl[9] = 550;
760     yl[9] = 150;
761     xl[10] = 650;
762     yl[10] = 140;
763    
764    
765    
766    
767    
768     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
769     gr1->SetMarkerColor(kBlue);
770     gr1->SetMarkerStyle(21);
771    
772    
773     //gr1->Draw("LP");
774    
775     TSpline3 *s = new TSpline3("grs",gr1);
776     s->SetLineColor(kBlue);
777     s->SetLineStyle(2);
778     s->SetLineWidth(3);
779    
780    
781     return gr1;
782    
783    
784    
785    
786    
787     }
788    
789    
790     TGraphErrors* getExpected_NLO_tanBeta50_up(){
791    
792     Int_t nl = 11;
793     Double_t xl[11];
794     Double_t yl[11];
795     Double_t exl[11];
796     Double_t eyl[11];
797    
798     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
799    
800     xl[0] = 200;
801     yl[0] = 340;
802     xl[1] = 220;
803     yl[1] = 320;
804     xl[2] = 245;
805     yl[2] = 320;
806     xl[3] = 270;
807     yl[3] = 310;
808     xl[4] = 300;
809     yl[4] = 290;
810     xl[5] = 350;
811     yl[5] = 270;
812     xl[6] = 400;
813     yl[6] = 245;
814     xl[7] = 450;
815     yl[7] = 220;
816     xl[8] = 500;
817     yl[8] = 200;
818     xl[9] = 550;
819     yl[9] = 180;
820     xl[10] = 650;
821     yl[10] = 170;
822    
823    
824    
825    
826     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
827     gr1->SetMarkerColor(kBlue);
828     gr1->SetMarkerStyle(21);
829    
830    
831     //gr1->Draw("LP");
832    
833     TSpline3 *s = new TSpline3("grs",gr1);
834     s->SetLineColor(kBlue);
835     s->SetLineStyle(2);
836     s->SetLineWidth(3);
837    
838    
839     return gr1;
840    
841    
842    
843    
844    
845     }
846    
847     TGraphErrors* getExpected_NLO_tanBeta50_low(){
848    
849     Int_t nl = 8;
850     Double_t xl[8];
851     Double_t yl[8];
852     Double_t exl[8];
853     Double_t eyl[8];
854    
855     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
856    
857     xl[0] = 200;
858     yl[0] = 280;
859     xl[1] = 220;
860     yl[1] = 270;
861     xl[2] = 245;
862     yl[2] = 262;
863     xl[3] = 270;
864     yl[3] = 242;
865     xl[4] = 300;
866     yl[4] = 227;
867     xl[5] = 350;
868     yl[5] = 197;
869     xl[6] = 400;
870     yl[6] = 180;
871     xl[7] = 450;
872     yl[7] = 150;
873     // xl[8] = 500;
874     //yl[8] = 168;
875     //xl[9] = 550;
876     //yl[9] = 145;
877    
878    
879    
880    
881    
882     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
883     gr1->SetMarkerColor(kBlue);
884     gr1->SetMarkerStyle(21);
885    
886    
887     //gr1->Draw("LP");
888    
889     TSpline3 *s = new TSpline3("grs",gr1);
890     s->SetLineColor(kBlue);
891     s->SetLineStyle(2);
892     s->SetLineWidth(3);
893    
894    
895     return gr1;
896    
897    
898    
899    
900    
901     }
902    
903    
904    
905    
906    
907    
908     TGraphErrors* getObserved_NLO_tanBeta3(){
909    
910     Int_t nl = 11;
911     Double_t xl[11];
912     Double_t yl[11];
913     Double_t exl[11];
914     Double_t eyl[11];
915    
916     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
917    
918    
919    
920     xl[0] = 0;
921     yl[0] = 274;
922     xl[1] = 100;
923     yl[1] = 270;
924     xl[2] = 150;
925     yl[2] = 268;
926     xl[3] = 200;
927     yl[3] = 265;
928     xl[4] = 250;
929     yl[4] = 255;
930     xl[5] = 300;
931     yl[5] = 230;
932     xl[6] = 350;
933     yl[6] = 195;
934     xl[7] = 400;
935     yl[7] = 160;
936     xl[8] = 450;
937     yl[8] = 140;
938     xl[9] = 480;
939     yl[9] = 130;
940     xl[10] = 530;
941     yl[10] = 120;
942    
943    
944    
945    
946     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
947     gr1->SetMarkerColor(kBlue);
948     gr1->SetMarkerStyle(21);
949    
950    
951     //gr1->Draw("LP");
952    
953     TSpline3 *s = new TSpline3("grs",gr1);
954     s->SetLineColor(kRed);
955     // s->SetLineStyle(2);
956     s->SetLineWidth(3);
957    
958    
959     return gr1;
960    
961    
962    
963    
964    
965     }
966    
967     TGraphErrors* getObserved_NLO_tanBeta3_2(){
968    
969     Int_t nl = 11;
970     Double_t xl[11];
971     Double_t yl[11];
972     Double_t exl[11];
973     Double_t eyl[11];
974    
975     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
976    
977    
978    
979     xl[0] = 0;
980     yl[0] = 277;
981     xl[1] = 100;
982     yl[1] = 272;
983     xl[2] = 150;
984     yl[2] = 268;
985     xl[3] = 200;
986     yl[3] = 263;
987     xl[4] = 250;
988     yl[4] = 258;
989     xl[5] = 300;
990     yl[5] = 237;
991     xl[6] = 350;
992     yl[6] = 198;
993     xl[7] = 400;
994     yl[7] = 168;
995     xl[8] = 450;
996     yl[8] = 152;
997     xl[9] = 480;
998     yl[9] = 145;
999     xl[10] = 530;
1000     yl[10] = 133;
1001    
1002    
1003    
1004    
1005     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
1006     gr1->SetMarkerColor(kBlue);
1007     gr1->SetMarkerStyle(21);
1008    
1009    
1010     //gr1->Draw("LP");
1011    
1012     TSpline3 *s = new TSpline3("grs",gr1);
1013     s->SetLineColor(kRed);
1014     // s->SetLineStyle(2);
1015     s->SetLineWidth(3);
1016    
1017    
1018     return gr1;
1019    
1020    
1021    
1022    
1023    
1024     }
1025    
1026    
1027    
1028     TGraphErrors* getObserved_NLO_tanBeta10(){
1029    
1030     Int_t nl = 11;
1031     Double_t xl[11];
1032     Double_t yl[11];
1033     Double_t exl[11];
1034     Double_t eyl[11];
1035    
1036     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
1037    
1038     xl[0] = 0;
1039     yl[0] = 278;
1040     xl[1] = 100;
1041     yl[1] = 270;
1042     xl[2] = 150;
1043     yl[2] = 267;
1044     xl[3] = 200;
1045     yl[3] = 262;
1046     xl[4] = 250;
1047     yl[4] = 250;
1048     xl[5] = 300;
1049     yl[5] = 225;
1050     xl[6] = 350;
1051     yl[6] = 192;
1052     xl[7] = 400;
1053     yl[7] = 163;
1054     xl[8] = 450;
1055     yl[8] = 148;
1056     xl[9] = 500;
1057     yl[9] = 140;
1058     xl[10] = 520;
1059     yl[10] = 137;
1060    
1061    
1062    
1063    
1064    
1065     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
1066     gr1->SetMarkerColor(kBlue);
1067     gr1->SetMarkerStyle(21);
1068    
1069    
1070     //gr1->Draw("LP");
1071    
1072     TSpline3 *s = new TSpline3("grs",gr1);
1073     s->SetLineColor(kRed);
1074     // s->SetLineStyle(2);
1075     s->SetLineWidth(3);
1076    
1077    
1078     return gr1;
1079    
1080    
1081    
1082    
1083    
1084     }
1085    
1086    
1087     TGraphErrors* getObserved_NLO_tanBeta10_2(){
1088    
1089     Int_t nl = 11;
1090     Double_t xl[11];
1091     Double_t yl[11];
1092     Double_t exl[11];
1093     Double_t eyl[11];
1094    
1095     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
1096    
1097     xl[0] = 0;
1098     yl[0] = 288;
1099     xl[1] = 100;
1100     yl[1] = 280;
1101     xl[2] = 150;
1102     yl[2] = 275;
1103     xl[3] = 200;
1104     yl[3] = 268;
1105     xl[4] = 250;
1106     yl[4] = 260;
1107     xl[5] = 300;
1108     yl[5] = 237;
1109     xl[6] = 350;
1110     yl[6] = 203;
1111     xl[7] = 400;
1112     yl[7] = 172;
1113     xl[8] = 450;
1114     yl[8] = 156;
1115     xl[9] = 500;
1116     yl[9] = 147;
1117     xl[10] = 650;
1118     yl[10] = 131;
1119    
1120    
1121    
1122    
1123    
1124     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
1125     gr1->SetMarkerColor(kBlue);
1126     gr1->SetMarkerStyle(21);
1127    
1128    
1129     //gr1->Draw("LP");
1130    
1131     TSpline3 *s = new TSpline3("grs",gr1);
1132     s->SetLineColor(kRed);
1133     // s->SetLineStyle(2);
1134     s->SetLineWidth(3);
1135    
1136    
1137     return gr1;
1138    
1139    
1140    
1141    
1142    
1143     }
1144    
1145    
1146     TGraphErrors* getObserved_NLO_tanBeta50(){
1147    
1148     Int_t nl = 10;
1149     Double_t xl[10];
1150     Double_t yl[10];
1151     Double_t exl[10];
1152     Double_t eyl[10];
1153    
1154     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
1155    
1156    
1157     xl[0] = 200;
1158     yl[0] = 243;
1159     xl[1] = 220;
1160     yl[1] = 266;
1161     xl[2] = 235;
1162     yl[2] = 278;
1163     xl[3] = 250;
1164     yl[3] = 267;
1165     xl[4] = 300;
1166     yl[4] = 230;
1167     xl[5] = 350;
1168     yl[5] = 205;
1169     xl[6] = 400;
1170     yl[6] = 184;
1171     xl[7] = 450;
1172     yl[7] = 168;
1173     xl[8] = 500;
1174     yl[8] = 156;
1175     xl[9] = 520;
1176     yl[9] = 148;
1177    
1178    
1179    
1180    
1181    
1182     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
1183     gr1->SetMarkerColor(kBlue);
1184     gr1->SetMarkerStyle(21);
1185    
1186    
1187     //gr1->Draw("LP");
1188    
1189     TSpline3 *s = new TSpline3("grs",gr1);
1190     s->SetLineColor(kRed);
1191     // s->SetLineStyle(2);
1192     s->SetLineWidth(3);
1193    
1194    
1195     return gr1;
1196    
1197    
1198    
1199    
1200    
1201     }
1202    
1203    
1204     TGraphErrors* getObserved_NLO_tanBeta50_2(){
1205    
1206     Int_t nl = 9;
1207     Double_t xl[9];
1208     Double_t yl[9];
1209     Double_t exl[9];
1210     Double_t eyl[9];
1211    
1212     // cout << " n " << hist->GetXaxis()->GetNbins() << endl;
1213    
1214    
1215     xl[0] = 200;
1216     yl[0] = 240;
1217     xl[1] = 220;
1218     yl[1] = 266;
1219     xl[2] = 235;
1220     yl[2] = 275;
1221     xl[3] = 250;
1222     yl[3] = 270;
1223     xl[4] = 300;
1224     yl[4] = 240;
1225     xl[5] = 350;
1226     yl[5] = 210;
1227     xl[6] = 400;
1228     yl[6] = 190;
1229     xl[7] = 450;
1230     yl[7] = 168;
1231     xl[8] = 500;
1232     yl[8] = 150;
1233    
1234    
1235    
1236    
1237    
1238    
1239     TGraphErrors* gr1 = new TGraphErrors(nl,xl,yl,exl,eyl);
1240     gr1->SetMarkerColor(kBlue);
1241     gr1->SetMarkerStyle(21);
1242    
1243    
1244     //gr1->Draw("LP");
1245    
1246     TSpline3 *s = new TSpline3("grs",gr1);
1247     s->SetLineColor(kRed);
1248     // s->SetLineStyle(2);
1249     s->SetLineWidth(3);
1250    
1251    
1252     return gr1;
1253    
1254    
1255    
1256    
1257    
1258     }
1259    
1260    
1261    
1262    
1263     TH1F* getHisto_1d( TString path,
1264     TString nameHist,
1265     TString nameFile ) {
1266     TString name = path + nameFile;
1267     TFile* file = new TFile(name);
1268     // TDirectory* dir = (TDirectory*)file->Get(Dirname);
1269     TH2F* hist = (TH2F*)file->Get(nameHist);
1270     if (!hist) {
1271     std::cout << " name: " << nameHist
1272     << " file: " << nameFile
1273     << std::endl;
1274     abort();
1275    
1276     }
1277    
1278    
1279    
1280     TH1F* Onedhist = new TH1F(nameHist,nameHist,int(hist->GetXaxis()->GetNbins()+0.5),hist->GetXaxis()->GetXmin(),hist->GetXaxis()->GetXmax());
1281    
1282     for(int x = 0;x < hist->GetXaxis()->GetNbins();x++){
1283    
1284     bool firsthit = false;
1285     for(int y = hist->GetYaxis()->GetNbins(); y>0; y--){
1286    
1287     double y_height = hist->GetYaxis()->GetXmin()+hist->GetYaxis()->GetBinWidth(0)*y;
1288    
1289     if(firsthit == false && hist->GetBinContent(x+1,y) > 0){
1290     Onedhist->SetBinContent(x+1,y_height);
1291     firsthit = true;
1292     }
1293    
1294    
1295     }
1296    
1297     }
1298    
1299    
1300     Onedhist->SetLineWidth(1);
1301    
1302     Onedhist->GetXaxis()->SetTitleSize(0.055);
1303     Onedhist->GetYaxis()->SetTitleSize(0.055);
1304     Onedhist->GetXaxis()->SetLabelSize(0.05);
1305     Onedhist->GetYaxis()->SetLabelSize(0.05);
1306     Onedhist->SetStats(kFALSE);
1307    
1308    
1309    
1310    
1311     return Onedhist;
1312     }
1313    
1314    
1315     #endif