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

# Content
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