ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/plot.cc
Revision: 1.20
Committed: Thu Mar 10 20:59:52 2011 UTC (14 years, 2 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.19: +358 -52 lines
Log Message: 
commiting style changes

File Contents

# User Rev Content
1 auterman 1.1 #include "plot.h"
2     #include "SusyScan.h"
3 auterman 1.2 #include "GeneratorMasses.h"
4 auterman 1.4 #include "PlotTools.h"
5 auterman 1.2 #include "TheLimits.h"
6 auterman 1.6 #include "GlobalFunctions.h"
7 auterman 1.19 #include "StyleSettings.h"
8 auterman 1.1
9     #include "TRint.h"
10     #include "TROOT.h"
11     #include "TObjArray.h"
12     #include "TStyle.h"
13    
14     #include "TChain.h"
15     #include "TFile.h"
16 auterman 1.2 #include "TGraph.h"
17 auterman 1.1 #include "TH1.h"
18     #include "TH2.h"
19     #include "TH2F.h"
20     #include "TTree.h"
21     #include "TKey.h"
22     #include "Riostream.h"
23     #include "TCanvas.h"
24 auterman 1.8 #include "TLegend.h"
25 auterman 1.9 #include "TLatex.h"
26 auterman 1.10 #include "TMarker.h"
27 auterman 1.20 #include "TPaveText.h"
28 auterman 1.1
29     #include <string>
30 auterman 1.2 #include <cmath>
31 auterman 1.9 #include <stdio.h>
32 auterman 1.1
33     int plot(int argc, char** argv)
34     {
35 auterman 1.6 //interactive root session
36 auterman 1.7 //TApplication theApp("App", 0, 0);
37 auterman 1.2 if (gROOT->IsBatch()) {
38     fprintf(stderr, "%s: cannot run in batch mode\n", argv[0]);
39     return 1;
40     }
41    
42 auterman 1.6 //Style stuff
43 auterman 1.19 /*
44 auterman 1.1 gStyle->SetHistFillColor(0);
45     gStyle->SetPalette(1);
46     gStyle->SetCanvasColor(0);
47     gStyle->SetCanvasBorderMode(0);
48     gStyle->SetPadColor(0);
49     gStyle->SetPadBorderMode(0);
50     gStyle->SetFrameBorderMode(0);
51    
52     gStyle->SetTitleFillColor(0);
53     gStyle->SetTitleBorderSize(0);
54     gStyle->SetTitleX(0.10);
55     gStyle->SetTitleY(0.98);
56     gStyle->SetTitleW(0.8);
57     gStyle->SetTitleH(0.06);
58    
59     gStyle->SetErrorX(0);
60     gStyle->SetStatColor(0);
61     gStyle->SetStatBorderSize(0);
62     gStyle->SetStatX(0);
63     gStyle->SetStatY(0);
64     gStyle->SetStatW(0);
65     gStyle->SetStatH(0);
66    
67     gStyle->SetTitleFont(22);
68     gStyle->SetLabelFont(22,"X");
69     gStyle->SetLabelFont(22,"Y");
70     gStyle->SetLabelFont(22,"Z");
71     gStyle->SetLabelSize(0.03,"X");
72     gStyle->SetLabelSize(0.03,"Y");
73     gStyle->SetLabelSize(0.03,"Z");
74 auterman 1.19 */
75     util::StyleSettings::paperNoTitle();
76     gStyle->SetPadBottomMargin(0.18);
77 auterman 1.1
78     //gROOT->SetStyle("MyStyle");
79     TCanvas * c1 = new TCanvas("c1","c1",600,600);
80 auterman 1.19 //c1->SetFillStyle ( 4000 );
81     //c1->SetLeftMargin ( 0.15 );
82     //c1->SetRightMargin ( 0.15 );
83     //c1->SetBottomMargin( 0.10 );
84 auterman 1.1 c1->cd();
85    
86 auterman 1.6 //Get limits and generator masses ---------------------------------------------------
87 auterman 1.4 TheLimits * genpoints = new TheLimits();
88 auterman 1.11 //genpoints->Fill(argc, argv);
89 auterman 1.20 //genpoints->Fill("limits_MHT_tb10_20110303/filelist.txt");
90     genpoints->Fill("limits_MHT_20110310/filelist_tb10.txt");
91 auterman 1.11
92     TheLimits * genpointsHT = new TheLimits();
93 auterman 1.20 //genpointsHT->Fill("limits_HT_tb10_20110303/filelist.txt");
94     genpointsHT->Fill("limits_HT_20110310/filelist_tb10.txt");
95 auterman 1.12
96 auterman 1.18 //TheLimits * genpointsHT = new TheLimits();
97     //genpointsHT->Fill("limits_HT_SigCont/filelist.txt");
98 auterman 1.13
99 auterman 1.18 //TheLimits * genpointsMHT = new TheLimits();
100     //genpointsMHT->Fill("limits_MHT_SigCont/filelist.txt");
101 auterman 1.11
102 auterman 1.10
103     //Replace read limits with specific numbers
104     //genpoints->OverwriteLimits("ABCD_MHT");
105    
106 auterman 1.17 genpoints->FillGeneratorMasses("GenScan_tb10.dat");
107 auterman 1.5 genpoints->match();
108 auterman 1.17 genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
109 auterman 1.11 genpointsHT->match();
110 auterman 1.18 //genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
111     //genpointsHT->match();
112     //genpointsMHT->FillGeneratorMasses("GenScan_tb10.dat");
113     //genpointsMHT->match();
114 auterman 1.2
115 auterman 1.6 //the plotting ----------------------------------------------------------------------
116 auterman 1.7 //plotting helper functions
117 auterman 1.4 PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
118 auterman 1.11 PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
119 auterman 1.18 //PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
120     //PlotTools<SusyScan> * plotToolsMHT = new PlotTools<SusyScan>(genpointsMHT->GetScan());
121 auterman 1.6 PlotTools<GeneratorMasses> * plotMasses = new PlotTools<GeneratorMasses>(genpoints->GetGeneratorMasses());
122 auterman 1.2
123 auterman 1.7 //iso mass lines
124 auterman 1.19 TGraph * gl500 = plotMasses->Line(Mzero, Mhalf, MGluino, 500, 0.8 );
125     TGraph * gl600 = plotMasses->Line(Mzero, Mhalf, MGluino, 600, 0.8 );
126     TGraph * gl700 = plotMasses->Line(Mzero, Mhalf, MGluino, 700, 0.8 );
127     TGraph * gl800 = plotMasses->Line(Mzero, Mhalf, MGluino, 800, 0.8 );
128     TGraph * gl900 = plotMasses->Line(Mzero, Mhalf, MGluino, 900, 0.8 );
129     TGraph * gl1000 = plotMasses->Line(Mzero, Mhalf, MGluino, 1000, 0.8);
130     TGraph * sq500 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 500, 1);
131     TGraph * sq600 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 600, 1);
132     TGraph * sq700 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 700, 1);
133     TGraph * sq800 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 800, 1);
134     TGraph * sq900 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 900, 1);
135     TGraph * sq1000 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 1000, 1);
136 auterman 1.7
137     TGraph * chi100 = plotMasses->Line(Mzero, Mhalf, MChi1, 50, 20);
138     TGraph * cha200 = plotMasses->Line(Mzero, Mhalf, MCha1, 200, 20);
139     cha200->SetLineColor(2);
140    
141     //the histograms
142     c1->SetLogz(1);
143     //h->SetMaximum(27);
144     //h->SetMinimum(0.01);
145 auterman 1.10
146     /**/
147 auterman 1.7 // cross-section in M0 - M1/2
148     TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
149 auterman 1.18 100,0,1009.9,35,100,450);
150 auterman 1.7 plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
151     hxsec->SetMinimum(0.01);
152     //sq500->Draw();
153     //gl500->Draw();
154     hxsec->Draw("colz");
155 auterman 1.10 c1->SaveAs("results/Xsection_m0_m12_tb10.pdf");
156 auterman 1.7 std::string wait;
157     //std::cin>>wait;
158 auterman 1.10
159 auterman 1.7 // Observed Limit in M0 - M1/2
160     TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
161 auterman 1.18 100,0,1009.9,35,100,450);
162 auterman 1.7 plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
163     hobslimit->SetMinimum(0.01);
164     hobslimit->Draw("colz");
165 auterman 1.14 c1->SaveAs("results/ObsLimit_m0_m12_tb10.pdf");
166    
167 auterman 1.15
168 auterman 1.7 // Expected Limit in M0 - M1/2
169     TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
170 auterman 1.18 100,0,1009.9,35,100,450);
171 auterman 1.7 plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
172     hexplimit->SetMinimum(0.01);
173     hexplimit->Draw("colz");
174 auterman 1.10 c1->SaveAs("results/ExpLimit_m0_m12_tb10.pdf");
175 auterman 1.7
176 auterman 1.13 // Signal Acceptance in M0 - M1/2
177 auterman 1.19 c1->SetRightMargin ( 0.15 );
178 auterman 1.7 TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
179 auterman 1.18 100,0,1009.9,35,100,450);
180 auterman 1.19 hsigacc->SetMinimum(0.01);
181     hsigacc->SetMaximum(1.0);
182 auterman 1.7 plotTools->Area(hsigacc, Mzero, Mhalf, SignalAcceptance);
183     hsigacc->SetMinimum(0.01);
184     hsigacc->SetMaximum(1.0);
185     hsigacc->Draw("colz");
186 auterman 1.14 //chi100->Draw();
187     //cha200->Draw();
188     //gl500 ->Draw();
189     // sq500 ->Draw();
190     c1->SaveAs("results/SigAccMHT_m0_m12_tb10.pdf");
191 auterman 1.18
192 auterman 1.14
193 auterman 1.19 // Signal Acceptance in M0 - M1/2
194     TH2F*hsigaccHT = new TH2F("sigaccHT",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
195     100,0,1009.9,35,100,450);
196     plotToolsHT->Area(hsigaccHT, Mzero, Mhalf, SignalAcceptance);
197     hsigaccHT->SetMinimum(0.01);
198     hsigaccHT->SetMaximum(1.0);
199     hsigaccHT->Draw("colz");
200     c1->SaveAs("results/SigAccHT_m0_m12_tb10.pdf");
201    
202    
203 auterman 1.7 // Exp. Limit on Number of Signal Events in M0 - M1/2
204 auterman 1.8 c1->SetLogz(0);
205 auterman 1.7 TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
206 auterman 1.18 100,0,1009.9,35,100,450);
207 auterman 1.7 plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
208 auterman 1.10 hexplimitnsig->SetMinimum(0.0);
209     hexplimitnsig->SetMaximum(20);
210 auterman 1.7 hexplimitnsig->Draw("colz");
211 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_m0_m12_tb10.pdf");
212 auterman 1.7
213     // Obs. Limit on Number of Signal Events in M0 - M1/2
214     TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
215 auterman 1.18 100,0,1009.9,35,100,450);
216 auterman 1.7 plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
217 auterman 1.10 hobslimitnsig->SetMinimum(0.0);
218     hobslimitnsig->SetMaximum(20);
219 auterman 1.7 hobslimitnsig->Draw("colz");
220 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_m0_m12_tb10.pdf");
221 auterman 1.7
222     c1->SetLogz(0);
223     // Expected Exclusion in M0 - M1/2
224     TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
225 auterman 1.8 100,0,1009.9,35,100,450);
226 auterman 1.11 plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
227 auterman 1.7 hexpexcl->Draw("colz");
228 auterman 1.10 c1->SaveAs("results/ExpExclusion_m0_m12_tb10.pdf");
229 auterman 1.7
230     // Observed Exclusion in M0 - M1/2
231     TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
232 auterman 1.8 100,0,1009.9,35,100,450);
233 auterman 1.11 plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
234 auterman 1.7 hobsexcl->Draw("colz");
235 auterman 1.10 c1->SaveAs("results/ObsExclusion_m0_m12_tb10.pdf");
236 auterman 1.7
237 auterman 1.8 // Observed Exclusion in M0 - M1/2
238     TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
239     100,0,1009.9,35,100,450);
240     plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
241     hPLobsexcl->Draw("colz");
242 auterman 1.10 c1->SaveAs("results/PL_ObsExclusion_m0_m12_tb10.pdf");
243 auterman 1.8
244 auterman 1.20 // TestContours in M0 - M1/2 /////////////////////////////////////////////////////////////////////
245 auterman 1.19 c1->SetRightMargin ( 0.1 );
246 auterman 1.14 TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
247     100,0,1009.9,35,100,450);
248 auterman 1.20 TH2F*testExclold=(TH2F*)texcl->Clone();
249     TH2F*testExclnew=(TH2F*)texcl->Clone();
250     plotTools->Area(testExclold, Mzero, Mhalf, NLOMCMCExpExclusion);
251     plotTools->Area(testExclnew, Mzero, Mhalf, NLOExpCL );
252    
253     std::vector<TGraph*> contours = plotTools->GetContours005(testExclnew,3);
254     //testExclold->Draw("colz");
255     testExclnew->Draw("colz");
256 auterman 1.14 int col=kBlue-10;
257 auterman 1.9 for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
258     if (! *cont) continue;
259     double x, y;
260     (*cont)->GetPoint(0, x, y);
261     (*cont)->SetLineColor(col);
262     (*cont)->Draw("l");
263     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col++);
264     char val[20];
265 auterman 1.13 sprintf(val,"%d",(int)(cont-contours.begin()));
266 auterman 1.9 l.DrawLatex(x,y,val);
267 auterman 1.11 //if (cont-contours.begin()>3) break;
268 auterman 1.9 }
269 auterman 1.10 c1->SaveAs("results/ExclusionTestContours_m0_m12_tb10.pdf");
270 auterman 1.9
271 auterman 1.20 /*
272 auterman 1.8 // Exclusion in M0 - M1/2
273     TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
274     100,0,1009.9,35,100,450);
275 auterman 1.19 TH2F*hs = new TH2F("hs","",100,0,1009.9,35,100,450);
276     TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
277     TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
278     TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
279 auterman 1.20 TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
280 auterman 1.19 TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
281     TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
282     TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
283     TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
284     TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
285     TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
286     TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
287     TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
288     TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
289     TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
290     TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
291     TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 1,1);
292     TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
293     TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
294     TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
295     TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
296     TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
297 auterman 1.20 TGraph * gCheck2 = (TGraph*)gCLsObsExclHT->Clone();
298     Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
299     Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
300     Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
301     Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
302     Smooth( gexpexcl, 27 );
303     Smooth( gobsexcl_LO, 27 );
304     Smooth( gobsexcl, 27 );
305     Smooth( gobsexcl, 27 );
306     Smooth( gexpexcl_LO, 27 );
307     Smooth( gFCexpexclHT, 27 );
308     Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
309     Smooth( gFCexpexclMHT, 27 );
310     Smooth( gFCexpexclMHT_LO, 27 );
311     Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
312     Smooth( gFCobsexclMHT, 27 );
313     Smooth( gCLsExpExclHTm1, 27 );
314     Smooth( gCLsExpExclHTp1, 27 );
315     Smooth( gCLsExpExclMHTm1, 27 );
316     Smooth( gCLsExpExclMHTp1, 27 );
317 auterman 1.19 TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
318 auterman 1.20 TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
319 auterman 1.19 TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
320     TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
321     TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
322 auterman 1.20 TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
323 auterman 1.19 TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
324     TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
325     TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
326     TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);
327     hexcl->GetYaxis()->SetTitleOffset(1.3);
328     hexcl->Draw("colz");
329 auterman 1.10 //set old exclusion Limits
330 auterman 1.16 TGraph* LEP_ch = set_lep_ch(10);
331     TGraph* LEP_sl = set_lep_sl(10);//slepton curve
332     TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
333     TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
334     //TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
335     //TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
336 auterman 1.20 TGraph* Atlas = Atlas_m0_m12_tb3_obs();
337     Atlas->SetLineColor( 28 );
338 auterman 1.16 TGraph* stau = set_tev_stau(10);//stau
339     TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
340     TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
341 auterman 1.15 TGraphErrors* First = getObserved_NLO_tanBeta10();
342     TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
343     TGraphErrors* Third = getLO_tanBeta10();
344 auterman 1.10 First->GetXaxis()->SetRangeUser(0,505);
345     First->GetYaxis()->SetRangeUser(80,500);
346     First->GetXaxis()->SetTitle("m_{0} (GeV)");
347     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
348     TSpline3 *sFirst = new TSpline3("sFirst",First);
349     sFirst->SetLineColor(kBlue);
350 auterman 1.18 sFirst->SetLineWidth(1);
351 auterman 1.10 TSpline3 *sSecond = new TSpline3("sSecond",Second);
352     sSecond->SetLineColor(kBlue);
353     sSecond->SetLineStyle(2);
354 auterman 1.18 sSecond->SetLineWidth(1);
355 auterman 1.10 TSpline3 *sThird = new TSpline3("sThird",Third);
356 auterman 1.18 sThird->SetLineColor(kBlue);
357 auterman 1.10 sThird->SetLineStyle(4);
358 auterman 1.18 sThird->SetLineWidth(1);
359 auterman 1.10 TEV_sn_d0_1->Draw("fsame");
360     TEV_sn_d0_2->Draw("fsame");
361 auterman 1.17 TEV_sg_d0->Draw("fsame");
362     TEV_sg_cdf->Draw("fsame");
363 auterman 1.10 LEP_ch->Draw("fsame");
364     LEP_sl->Draw("fsame");
365     stau->Draw("fsame");
366     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
367 auterman 1.19 b.DrawLatex( 10,380,"#tilde{#tau} LSP");
368 auterman 1.10 b.DrawLatex( 6,150,"LEP2");
369     b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
370     b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
371 auterman 1.17 b.DrawLatex(90,160,"CDF");
372     //b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
373 auterman 1.10 b.DrawLatex( 80,180,"D0 #tilde{#nu}");
374 auterman 1.18 if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
375     sFirst->Draw("same");
376 auterman 1.20 //sSecond->Draw("same");
377     //sThird->Draw("same");
378     Atlas->Draw("c,same");
379     TLatex ms; ms.SetTextSize(0.02); ms.SetTextColor(12);
380     ms.DrawLatex(150,433,"tan#beta=10, #mu>0, A_{0}=0");
381 auterman 1.19 //gCheck1->Draw("same");
382     //gCheck2->Draw("same");
383     //if (gexpexcl) gexpexcl->Draw("l");
384     //if (gobsexcl) gobsexcl->Draw("l");
385 auterman 1.18 //if (gexpexcl_LO) gexpexcl_LO->Draw("l");
386 auterman 1.19 //if (gMCMCexpexcl) gMCMCexpexcl->Draw("l");
387 auterman 1.18 //if (gCLsObsExclHT_LO) gCLsObsExclHT_LO->Draw("l");
388 auterman 1.19 //if (gCLsExpExclHT_LO) gCLsExpExclHT_LO->Draw("l");
389     //if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
390     //if (gCLsExpExclHT) gCLsExpExclHT->Draw("l");
391 auterman 1.18 //if (gobsexcl) gobsexcl->Draw("l");
392     //if (gPLobsexcl) gPLobsexcl->Draw("l");
393     //if (gFCobsexcl) gFCobsexcl->Draw("l");
394     //if (gMCMCobsexcl) gMCMCobsexcl->Draw("l");
395     //if (gFCexpexcl) gFCexpexcl->Draw("l");
396     //if (gPLexpexcl) gPLexpexcl->Draw("l");
397     //if (gFCexpexcl) gFCexpexcl->Draw("l");
398 auterman 1.19
399     CLsObsNLO->Draw("l");
400 auterman 1.20 CLsObsLO->Draw("l");
401 auterman 1.19 CLsExpNLO->Draw("l");
402 auterman 1.20 //CLsExpLO->Draw("l");
403     FCExpNLO->Draw("l");
404 auterman 1.19 //FCObsNLO->Draw("l");
405    
406     //gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c");
407     //sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c");
408 auterman 1.20 TLegend * leg = new TLegend(0.54,0.7,0.89,0.93);
409 auterman 1.19 leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);
410     TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
411     leg->SetHeader("CMS preliminary, 95%CL");
412 auterman 1.20 leg->AddEntry(expLeg, "Expected (CLs, NLO)","lf");
413 auterman 1.19 leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l");
414 auterman 1.20 leg->AddEntry(CLsObsLO, "Observed (CLs, LO)","l");
415     leg->AddEntry(FCExpLO, "Expected (FC, NLO)","l");
416     //leg->AddEntry(FCExpLO, "Expected (FC, LO)","l");
417 auterman 1.12 //if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
418 auterman 1.19 //if (gCLsExpExclHT) leg->AddEntry(expLeg ,"NLO Expected (HT, CLs)","l");
419     //if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
420 auterman 1.18 //if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
421 auterman 1.19 //if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
422 auterman 1.11 //if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
423     //if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
424 auterman 1.19 //if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l");
425     //if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l");
426 auterman 1.12 //if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
427 auterman 1.19 //if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
428     //if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
429     //if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l");
430 auterman 1.20 //if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)");
431     if (sFirst) leg->AddEntry(sFirst, "CMS #alpha_{T} Observed");
432     //if (sThird) leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)");
433     if (Atlas) leg->AddEntry(Atlas, "Atlas Observed (tan#beta=3)","l");
434 auterman 1.8 leg->Draw();
435 auterman 1.10 gPad->RedrawAxis();
436 auterman 1.13 c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
437 auterman 1.20 */
438     // Exclusion in M0 - M1/2
439     TH2F*hexcl = new TH2F("hexcl",";m_{0} (GeV); m_{1/2} (GeV); 95% CL Exclusion",
440     100,0,1009.9,35,110,450);
441     TH2F*hs = new TH2F("hs","",100,0,1009.9,35,100,450);
442     TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
443     TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
444     TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
445     TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
446     TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
447     TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
448     TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
449     TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
450     TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
451     TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
452     TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
453     TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
454     TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
455     TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
456     TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
457     TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
458     TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
459     TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
460     TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
461     TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
462     TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
463     TGraph * gCheck2 = (TGraph*)gexpexcl->Clone();
464     Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
465     Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
466     Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
467     Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
468     Smooth( gexpexcl, 27 );
469     Smooth( gobsexcl_LO, 27 );
470     //Smooth( gobsexcl, 27 );
471     Smooth( gobsexcl, 27 );
472     Smooth( gexpexcl_LO, 27 );
473     Smooth( gFCexpexclHT, 27 );
474     Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
475     Smooth( gFCexpexclMHT, 27 );
476     Smooth( gFCexpexclMHT_LO, 27 );
477     Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
478     Smooth( gFCobsexclMHT, 27 );
479     Smooth( gCLsExpExclHTm1, 27 );
480     Smooth( gCLsExpExclHTp1, 27 );
481     Smooth( gCLsExpExclMHTm1, 27 );
482     Smooth( gCLsExpExclMHTp1, 27 );
483     TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
484     TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
485     TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
486     TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
487     TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
488     TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
489     TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
490     TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
491     TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
492     TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);gCLsExp1Sigma->SetFillStyle(4010);
493     hexcl->GetYaxis()->SetTitleOffset(1.3);
494     hexcl->GetXaxis()->SetTitleOffset(0.92);
495     hexcl->Draw("colz");
496     //set old exclusion Limits
497     TGraph* LEP_ch = set_lep_ch(10);
498     TGraph* LEP_sl = set_lep_sl(10);//slepton curve
499     TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
500     TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
501     //TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
502     //TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
503     TGraph* Atlas = Atlas_m0_m12_tb3_obs();
504     Atlas->SetLineColor( 28 );
505     TGraph* stau = set_tev_stau(10);//stau
506     TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
507     TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
508     TGraphErrors* First = getObserved_NLO_tanBeta10();
509     TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
510     TGraphErrors* Third = getLO_tanBeta10();
511     First->GetXaxis()->SetRangeUser(0,505);
512     First->GetYaxis()->SetRangeUser(80,500);
513     First->GetXaxis()->SetTitle("m_{0} (GeV)");
514     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
515     TSpline3 *sFirst = new TSpline3("sFirst",First); sFirst->SetLineColor(kBlue);sFirst->SetLineWidth(1);
516     TSpline3 *sSecond = new TSpline3("sSecond",Second);
517     sSecond->SetLineColor(kBlue); sSecond->SetLineStyle(2); sSecond->SetLineWidth(1);
518     TSpline3 *sThird = new TSpline3("sThird",Third);
519     sThird->SetLineColor(kBlue); sThird->SetLineStyle(4); sThird->SetLineWidth(1);
520     //TEV_sn_d0_1->Draw("fsame");
521     //TEV_sn_d0_2->Draw("fsame"); //only for tb=3
522     LEP_ch->Draw("fsame");
523     LEP_sl->Draw("fsame");
524     TEV_sg_cdf->Draw("fsame");
525     TEV_sg_cdf->Draw("lsame");
526     TEV_sg_d0->Draw("fsame");
527     TEV_sg_d0->Draw("lsame");
528     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
529     //b.DrawLatex( 6,150,"LEP2");
530     //b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
531     //b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
532     //b.DrawLatex(90,160,"CDF");
533     //b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
534     //b.DrawLatex( 80,180,"D0 #tilde{#nu}");
535     if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
536     sFirst->Draw("same");
537     //sSecond->Draw("same");
538     //sThird->Draw("same");
539     Atlas->Draw("c,same");
540     TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
541     ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
542     //gCheck1->Draw("same");
543     //gCheck2->Draw("same");
544     //LM points
545     TMarker* LM0 = new TMarker(200.,160.,20);
546     TMarker* LM1 = new TMarker(60.,250.,20);
547     LM0->SetMarkerSize(0.7);
548     LM1->SetMarkerSize(0.7);
549     TLatex* tLM0 = new TLatex(205.,160.,"LM0");
550     tLM0->SetTextSize(0.02);
551     TLatex* tLM1 = new TLatex(65.,243.,"LM1");
552     tLM1->SetTextSize(0.02);
553     LM0->Draw("same");
554     tLM0->Draw("same");
555     LM1->Draw("same");
556     tLM1->Draw("same");
557     TLegend* legexp = new TLegend(0.61,0.62,0.93,0.93,NULL,"brNDC");
558     legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
559     //TEV_sg_cdf.SetLineColor(1);
560     legexp->SetHeader("CMS preliminary");
561     legexp->AddEntry(TEV_sg_cdf,"CDF #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
562     legexp->AddEntry(TEV_sg_d0,"D0 #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=3, #mu<0}","f");
563     //ch_gr.SetLineColor(1);
564     legexp->AddEntry(LEP_ch,"LEP2 #tilde{#chi}_{1}^{#pm}","f");
565     //sl_gr.SetLineColor(1);
566     legexp->AddEntry(LEP_sl,"LEP2 #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
567     //if(tanbeta == 3)
568     //legexp->AddEntry(TEV_sn_d0_1,"D0 #chi^{#pm}_{1}, #chi^{0}_{2}","f");
569     legexp->AddEntry(sFirst, "CMS #alpha_{T}");
570     legexp->AddEntry(Atlas, "Atlas, #scale[0.8]{tan#beta=3}","l");
571     legexp->Draw();
572    
573     CLsObsNLO->Draw("l");
574     CLsObsLO->Draw("l");
575     CLsExpNLO->Draw("l");
576     //FCExpNLO->Draw("l");
577     stau->Draw("fsame");
578     b.DrawLatex( 10,380,"#tilde{#tau} LSP");
579 auterman 1.15
580 auterman 1.20 //constant ssqquark and gluino lines
581     TF1* lnsq[4];
582     TF1* lngl[4];
583     TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
584     sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
585     sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
586     TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
587     glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
588     glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
589     glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
590     int tanBeta_=10;
591     for(int i = 0; i < 4; i++){
592     lnsq[i] = constant_squark(tanBeta_,i);
593     lngl[i] = constant_gluino(tanBeta_,i);
594     lngl[i]->Draw("same");
595     lnsq[i]->Draw("same");
596     }
597     TLegend * leg = new TLegend(0.25,0.80,0.6,0.93);
598     leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
599     TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
600     leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
601     leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
602     leg->AddEntry(CLsObsLO, "Observed, LO","l");
603     leg->AddEntry(expLeg, "Expected #pm 1#sigma, NLO","lf");
604     //leg->AddEntry(FCExpLO, "Expected, FC, NLO","l");
605     leg->Draw();
606     gPad->RedrawAxis();
607     c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
608    
609     hexcl->Draw("");
610     gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
611     CLsObsNLO->SetName("CLsNLO_ObservedTb10");
612     CLsExpNLO->SetName("CLsNLO_ExpectedTb10");
613     CLsExpLO->SetName("CLsLO_ExpectedTb10");
614     gCLsExp1Sigma->Draw("l");
615     CLsObsNLO->Draw("l");
616     CLsExpNLO->Draw("l");
617     CLsObsLO->Draw("l");
618     c1->SaveAs("results/LimitContours_tb10.C");
619     // ===============================================================================================================
620 auterman 1.8
621 auterman 1.13 // Signal Contamination in M0 - M1/2
622     c1->SetLogz(0);
623     TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
624 auterman 1.18 100,0,1009.9,35,100,450);
625     plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
626 auterman 1.13 hsigcont->SetMinimum(0.01);
627     hsigcont->SetMaximum(20);
628     hsigcont->Draw("colz");
629     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
630     c1->SaveAs("results/SignalContamination_HT_m0_m12_tb10.pdf");
631    
632     // Signal Contamination in M0 - M1/2
633     c1->SetLogz(0);
634     TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
635 auterman 1.18 100,0,1009.9,35,100,450);
636     plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
637 auterman 1.13 hsigcontMHT->SetMinimum(0.01);
638     hsigcontMHT->SetMaximum(20);
639     hsigcontMHT->Draw("colz");
640     if (gobsexcl) gobsexcl->Draw("l");
641     c1->SaveAs("results/SignalContamination_MHT_m0_m12_tb10.pdf");
642    
643 auterman 1.18 // Signal JEC Uncertainty MHT
644     c1->SetLogz(0);
645     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
646     100,0,1009.9,35,100,450);
647     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
648     hsig_jec_mht->SetMinimum(0.0);
649     hsig_jec_mht->SetMaximum(0.25);
650     hsig_jec_mht->SetContour(10);
651     hsig_jec_mht->Draw("colz");
652     if (gobsexcl) gobsexcl->Draw("l");
653     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
654    
655     // Signal JEC Uncertainty HT
656     c1->SetLogz(0);
657     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
658     100,0,1009.9,35,100,450);
659 auterman 1.19 plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
660 auterman 1.18 hsig_jec_ht->SetMinimum(0.0);
661     hsig_jec_ht->SetMaximum(0.25);
662     hsig_jec_ht->SetContour(10);
663     hsig_jec_ht->Draw("colz");
664     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
665     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
666    
667    
668     // Signal MuIso Uncertainty MHT
669     c1->SetLogz(0);
670     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
671     100,0,1009.9,35,100,450);
672     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
673 auterman 1.19 hsig_MuIso_mht->SetMinimum(0.0);hsig_MuIso_mht->GetZaxis()->SetTitleOffset(1.3);
674 auterman 1.18 hsig_MuIso_mht->SetMaximum(0.05);
675     hsig_MuIso_mht->SetContour(5);
676     hsig_MuIso_mht->Draw("colz");
677     if (gobsexcl) gobsexcl->Draw("l");
678     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
679    
680     // Signal MuIso Uncertainty HT
681     c1->SetLogz(0);
682     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
683     100,0,1009.9,35,100,450);
684 auterman 1.19 plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
685     hsig_MuIso_ht->SetMinimum(0.0);hsig_MuIso_ht->GetZaxis()->SetTitleOffset(1.3);
686 auterman 1.18 hsig_MuIso_ht->SetMaximum(0.05);
687     hsig_MuIso_ht->SetContour(5);
688     hsig_MuIso_ht->Draw("colz");
689     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
690     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
691    
692    
693     // Signal kFactor Uncertainty MHT
694     c1->SetLogz(0);
695 auterman 1.19 TH2F*hsig_kFactorUncert_mht = new TH2F("sigkFactorUncertmht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
696 auterman 1.18 100,0,1009.9,35,100,450);
697 auterman 1.19 plotTools->Area(hsig_kFactorUncert_mht, Mzero, Mhalf, SignalUncertKfactor);
698     hsig_kFactorUncert_mht->SetMinimum(0.1);hsig_kFactorUncert_mht->GetZaxis()->SetTitleOffset(1.3);
699     hsig_kFactorUncert_mht->SetMaximum(0.2);
700     hsig_kFactorUncert_mht->SetContour(10);
701     hsig_kFactorUncert_mht->Draw("colz");
702 auterman 1.18 if (gobsexcl) gobsexcl->Draw("l");
703 auterman 1.19 c1->SaveAs("results/SigkFactorUncert_MHT_m0_m12_tb10.pdf");
704 auterman 1.18
705     // Signal kFactor Uncertainty HT
706     c1->SetLogz(0);
707 auterman 1.19 TH2F*hsig_kFactorUncert_ht = new TH2F("sigkFactorUncertht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
708 auterman 1.18 100,0,1009.9,35,100,450);
709 auterman 1.19 plotToolsHT->Area(hsig_kFactorUncert_ht, Mzero, Mhalf, SignalUncertKfactor);
710     hsig_kFactorUncert_ht->SetMinimum(0.1);hsig_kFactorUncert_ht->GetZaxis()->SetTitleOffset(1.3);
711     hsig_kFactorUncert_ht->SetMaximum(0.2);
712     hsig_kFactorUncert_ht->SetContour(10);
713     hsig_kFactorUncert_ht->Draw("colz");
714     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
715     c1->SaveAs("results/SigkFactorUncert_HT_m0_m12_tb10.pdf");
716    
717     // Signal kFactor HT
718     c1->SetLogz(0);
719     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
720     100,0,1009.9,35,100,450);
721     plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalKfactor);
722     hsig_kFactor_ht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
723     hsig_kFactor_ht->SetMaximum(2.4);
724     hsig_kFactor_ht->SetContour(7);
725 auterman 1.18 hsig_kFactor_ht->Draw("colz");
726     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
727     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
728    
729 auterman 1.19 // Signal kFactor MHT
730     c1->SetLogz(0);
731     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
732     100,0,1009.9,35,100,450);
733     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalKfactor);
734     hsig_kFactor_mht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
735     hsig_kFactor_mht->SetMaximum(2.4);
736     hsig_kFactor_mht->SetContour(7);
737     hsig_kFactor_mht->Draw("colz");
738     if (gobsexcl) gobsexcl->Draw("l");
739     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
740    
741 auterman 1.1
742 auterman 1.2
743 auterman 1.7 //-----------------------------------------------------------------------------------
744     c1->SetLogz(1);
745    
746     TGraph * mz500 = plotMasses->Line( MSquarkL, MGluino, Mzero, 500, 1);
747     TGraph * mh250 = plotMasses->Line( MSquarkL, MGluino, Mhalf, 250, 1);
748    
749 auterman 1.8 // cross-section in squark - gluino mass
750 auterman 1.7 TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
751     60,200,1400,50,200,1200);
752     plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
753     hxsec_qg->SetMinimum(0.01);
754     hxsec_qg->Draw("colz");
755     mz500->Draw();
756     mh250->Draw();
757 auterman 1.10 c1->SaveAs("results/Xsection_mSql_mGl_tb10.pdf");
758 auterman 1.7
759 auterman 1.8 // Observed Limit in squark - gluino mass
760 auterman 1.7 TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
761     60,200,1400,50,200,1200);
762     plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
763     hobslimit_qg->SetMinimum(0.01);
764     hobslimit_qg->Draw("colz");
765 auterman 1.10 c1->SaveAs("results/ObsLimit_mSql_mGl_tb10.pdf");
766 auterman 1.7
767 auterman 1.8 // Expected Limit in squark - gluino mass
768 auterman 1.7 TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
769     60,200,1400,50,200,1200);
770     plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
771     hexplimit_qg->SetMinimum(0.01);
772     hexplimit_qg->Draw("colz");
773 auterman 1.10 c1->SaveAs("results/ExpLimit_mSql_mGl_tb10.pdf");
774 auterman 1.7
775 auterman 1.8 // Signal Acceptance in squark - gluino mass
776 auterman 1.7 TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
777     60,200,1400,50,200,1200);
778     plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
779     hsigacc_qg->SetMinimum(0.01);
780     hsigacc_qg->SetMaximum(1.0);
781     hsigacc_qg->Draw("colz");
782 auterman 1.14 c1->SaveAs("results/SigAccMHT_mSql_mGl_tb10.pdf");
783    
784     // Signal Acceptance in squark - gluino mass
785     TH2F*hsigaccHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
786     60,200,1400,50,200,1200);
787     plotToolsHT->Area(hsigaccHT_qg, MSquarkL, MGluino, SignalAcceptance);
788     hsigaccHT_qg->SetMinimum(0.01);
789     hsigaccHT_qg->SetMaximum(1.0);
790     hsigaccHT_qg->Draw("colz");
791     c1->SaveAs("results/SigAccHT_mSql_mGl_tb10.pdf");
792 auterman 1.19
793    
794     // Signal Contamination in squark - gluino mass
795     c1->SetLogz(0);
796     TH2F*hsigcont_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Contamination in bkg yield",
797     60,200,1400,50,200,1200);
798     plotTools->Area(hsigcont_qg, MSquarkL, MGluino, SignalContamination);
799     hsigcont_qg->SetMinimum(0.0);
800     hsigcont_qg->SetMaximum(20.0);
801     hsigcont_qg->Draw("colz");
802     c1->SaveAs("results/SignalContamination_MHT_mSql_mGl_tb10.pdf");
803    
804     // Signal Contamination in squark - gluino mass
805     TH2F*hsigcontHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance in bkg yield",
806     60,200,1400,50,200,1200);
807     plotToolsHT->Area(hsigcontHT_qg, MSquarkL, MGluino, SignalContamination);
808     hsigcontHT_qg->SetMinimum(0.0);
809     hsigcontHT_qg->SetMaximum(20.0);
810     hsigcontHT_qg->Draw("colz");
811     c1->SaveAs("results/SignalContamination_HT_mSql_mGl_tb10.pdf");
812 auterman 1.18 /*
813     // Signal JEC Uncertainty MHT
814     c1->SetLogz(0);
815     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
816     100,0,1009.9,35,100,450);
817     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
818     hsig_jec_mht->SetMinimum(0.0);
819     hsig_jec_mht->SetMaximum(0.25);
820     hsig_jec_mht->SetContour(10);
821     hsig_jec_mht->Draw("colz");
822     if (gobsexcl) gobsexcl->Draw("l");
823     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
824    
825     // Signal JEC Uncertainty HT
826     c1->SetLogz(0);
827     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
828     100,0,1009.9,35,100,450);
829     plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
830     hsig_jec_ht->SetMinimum(0.0);
831     hsig_jec_ht->SetMaximum(0.25);
832     hsig_jec_ht->SetContour(10);
833     hsig_jec_ht->Draw("colz");
834     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
835     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
836    
837 auterman 1.15
838 auterman 1.18 // Signal MuIso Uncertainty MHT
839     c1->SetLogz(0);
840     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
841     100,0,1009.9,35,100,450);
842     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
843     hsig_MuIso_mht->SetMinimum(0.0);
844     hsig_MuIso_mht->SetMaximum(0.05);
845     hsig_MuIso_mht->SetContour(5);
846     hsig_MuIso_mht->Draw("colz");
847     if (gobsexcl) gobsexcl->Draw("l");
848     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
849    
850     // Signal MuIso Uncertainty HT
851     c1->SetLogz(0);
852     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
853     100,0,1009.9,35,100,450);
854     plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
855     hsig_MuIso_ht->SetMinimum(0.0);
856     hsig_MuIso_ht->SetMaximum(0.05);
857     hsig_MuIso_ht->SetContour(5);
858     hsig_MuIso_ht->Draw("colz");
859     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
860     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
861    
862    
863     // Signal kFactor Uncertainty MHT
864     c1->SetLogz(0);
865     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
866     100,0,1009.9,35,100,450);
867     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
868     hsig_kFactor_mht->SetMinimum(0.1);
869     hsig_kFactor_mht->SetMaximum(0.2);
870     hsig_kFactor_mht->SetContour(10);
871     hsig_kFactor_mht->Draw("colz");
872     if (gobsexcl) gobsexcl->Draw("l");
873     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
874    
875     // Signal kFactor Uncertainty HT
876     c1->SetLogz(0);
877     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
878     100,0,1009.9,35,100,450);
879     plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
880     hsig_kFactor_ht->SetMinimum(0.1);
881     hsig_kFactor_ht->SetMaximum(0.2);
882     hsig_kFactor_ht->SetContour(10);
883     hsig_kFactor_ht->Draw("colz");
884     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
885     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
886 auterman 1.7
887 auterman 1.18 */
888 auterman 1.8 // Exp. Limit on Number of Signal Events in squark - gluino mass
889     c1->SetLogz(0);
890 auterman 1.7 TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
891     60,200,1400,50,200,1200);
892     plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
893 auterman 1.10 hexplimitnsig_qg->SetMinimum(0.0);
894     hexplimitnsig_qg->SetMaximum(20);
895 auterman 1.7 hexplimitnsig_qg->Draw("colz");
896 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_mSql_mGl_tb10.pdf");
897 auterman 1.6
898 auterman 1.8 // Obs. Limit on Number of Signal Events in squark - gluino mass
899 auterman 1.7 TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
900     60,200,1400,50,200,1200);
901     plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
902 auterman 1.10 hobslimitnsig_qg->SetMinimum(0.0);
903     hobslimitnsig_qg->SetMaximum(20);
904 auterman 1.7 hobslimitnsig_qg->Draw("colz");
905 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_mSql_mGl_tb10.pdf");
906 auterman 1.7
907     c1->SetLogz(0);
908 auterman 1.8 // Expected Exclusion in squark - gluino mass
909 auterman 1.7 TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
910 auterman 1.2 60,200,1400,50,200,1200);
911 auterman 1.11 plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
912 auterman 1.7 hexpexcl_qg->Draw("colz");
913 auterman 1.10 c1->SaveAs("results/ExpExclusion_mSql_mGl_tb10.pdf");
914 auterman 1.2
915 auterman 1.8 // Observed Exclusion in squark - gluino mass
916 auterman 1.7 TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
917     60,200,1400,50,200,1200);
918 auterman 1.11 plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
919 auterman 1.8 //std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
920 auterman 1.7 hobsexcl_qg->Draw("colz");
921 auterman 1.8 //for (std::vector<TGraph *>::iterator g=vobsexcl_qg.begin();g!=vobsexcl_qg.end();++g){
922     // if (! *g) continue;
923     // // (*g)->Draw("l");
924     //}
925 auterman 1.10 c1->SaveAs("results/ObsExclusion_mSql_mGl_tb10.pdf");
926 auterman 1.2
927 auterman 1.9 // TestContours in M0 - M1/2
928     TH2F*texcl_qg = new TH2F("texcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
929     60,200,1400,50,200,1200);
930     TH2F*tPLobsexcl_qg=(TH2F*)texcl_qg->Clone();
931     plotTools->Area(tPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
932 auterman 1.10 TH2F*tPLexpexcl_qg=(TH2F*)texcl_qg->Clone();
933     plotTools->Area(tPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
934 auterman 1.9 //std::vector<TGraph*> contours_qg = plotTools->GetContours(hexpexcl_qg,3);
935     //std::vector<TGraph*> contours_qg = plotTools->GetContours(hobsexcl_qg,3);
936     std::vector<TGraph*> contours_qg = plotTools->GetContours(tPLobsexcl_qg,3);
937 auterman 1.10 //std::vector<TGraph*> contours_qg = plotTools->GetContours(tPLexpexcl_qg,3);
938 auterman 1.9 //hPLexpexcl_qg
939     //hexcl_qg->Draw("colz");
940     //hexpexcl_qg->Draw("colz");
941     //hobsexcl_qg->Draw("colz");
942     tPLobsexcl_qg->Draw("colz");
943     int col_gl=kBlue-10;
944     for (std::vector<TGraph*>::iterator cont=contours_qg.begin(); cont!=contours_qg.end(); ++cont){
945     if (! *cont) continue;
946     double x, y;
947     (*cont)->GetPoint(0, x, y);
948     (*cont)->SetLineColor(col_gl);(*cont)->SetLineWidth(3);
949     (*cont)->Draw("l");
950     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col_gl++);
951     char val[20];
952 auterman 1.13 sprintf(val,"%d",(int)(cont-contours_qg.begin()));
953 auterman 1.9 l.DrawLatex(x,y,val);
954 auterman 1.10 if (cont-contours_qg.begin()>3) break;
955 auterman 1.9 }
956 auterman 1.10 c1->SaveAs("results/ExclusionTestContours_mSql_mGl_tb10.pdf");
957 auterman 1.9
958 auterman 1.8 // Observed Exclusion in squark - gluino mass
959     TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
960     60,200,1400,50,200,1200);
961     plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
962 auterman 1.10 //hPLobsexcl_qg->Draw("colz");
963     //c1->SaveAs("results/PL_ObsExclusion_mSql_mGl_tb10.pdf");
964 auterman 1.8
965 auterman 1.11
966    
967 auterman 1.10 // Exclusion in squark - gluino mass ----------------------------------------
968 auterman 1.19 TH2F*h_qg = new TH2F("hExcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
969     60,0,1200,50,0,1200);
970     TH2F*hExcl_qg = new TH2F("h_qg","",60,200,1400,50,200,1200);
971     TH2F*hexpFC_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
972     plotToolsHT->Area(hexpFC_HT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
973     TH2F*hexpFC_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
974     plotToolsHT->Area(hexpFC_HT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
975     TH2F*hexpFC_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
976     plotTools->Area(hexpFC_MHT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
977     TH2F*hexpFC_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
978     plotTools->Area(hexpFC_MHT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
979     //
980     TH2F*hexpCLs_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
981     plotToolsHT->Area(hexpCLs_HT_LO_qg, MSquarkL, MGluino, ExpExclCL);
982     TH2F*hexpCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
983     plotToolsHT->Area(hexpCLs_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
984     TH2F*hexpCLs_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
985     plotTools->Area(hexpCLs_MHT_LO_qg, MSquarkL, MGluino, ExpExclCL);
986     TH2F*hexpCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
987     plotTools->Area(hexpCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
988     //
989     TH2F*hexp1p_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
990     plotToolsHT->Area(hexp1p_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
991     TH2F*hexp1p_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
992     plotTools->Area(hexp1p_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
993     TH2F*hexp1m_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
994     plotToolsHT->Area(hexp1m_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
995     TH2F*hexp1m_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
996     plotTools->Area(hexp1m_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
997    
998     TH2F*hobsCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
999     plotToolsHT->Area(hobsCLs_HT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
1000     TH2F*hobsCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1001     plotTools->Area(hobsCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
1002     //
1003     TH2 * hexpFC_LO_qg = plotToolsHT->BinWiseOr(hexpFC_HT_LO_qg, hexpFC_MHT_LO_qg);
1004     TH2 * hexpCLs_LO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_qg, hexpCLs_MHT_LO_qg);
1005     TH2 * hexpCLs_NLO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_qg, hexpCLs_MHT_NLO_qg);
1006     TH2 * hobsCLs_NLO_qg = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_qg, hobsCLs_MHT_NLO_qg);
1007     TH2 * hexp1p_NLO_qg = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_qg, hexp1p_MHT_NLO_qg);
1008     TH2 * hexp1m_NLO_qg = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_qg, hexp1m_MHT_NLO_qg);
1009     //
1010     TGraph * gexpFC_LO_qg = plotToolsHT->GetContour(hexpFC_LO_qg,3,0);
1011     gexpFC_LO_qg->SetLineStyle(4); gexpFC_LO_qg->SetLineColor(3);gexpFC_LO_qg->SetLineWidth(2);
1012     TGraph * gexpCLs_LO_qg = plotToolsHT->GetContour(hexpCLs_LO_qg,3,0);
1013     gexpCLs_LO_qg->SetLineStyle(4); gexpCLs_LO_qg->SetLineColor(1);gexpCLs_LO_qg->SetLineWidth(2);
1014     TGraph * gexpCLs_NLO_qg = plotToolsHT->GetContour(hexpCLs_NLO_qg,3,0);
1015     gexpCLs_NLO_qg->SetLineStyle(2); gexpCLs_NLO_qg->SetLineColor(1);gexpCLs_NLO_qg->SetLineWidth(2);
1016     TGraph * gobsCLs_NLO_qg = plotToolsHT->GetContour(hobsCLs_NLO_qg,3,0);
1017     gobsCLs_NLO_qg->SetLineStyle(1); gobsCLs_NLO_qg->SetLineColor(1);gobsCLs_NLO_qg->SetLineWidth(2);
1018     TGraph * gexp1p_NLO_qg = plotToolsHT->GetContour(hexp1p_NLO_qg,3,0);
1019     TGraph * gexp1m_NLO_qg = plotToolsHT->GetContour(hexp1m_NLO_qg,3,0);
1020     Smooth(gexpFC_LO_qg, 11);
1021     Smooth(gexpCLs_LO_qg, 11);
1022     Smooth(gexpCLs_NLO_qg, 11);
1023     Smooth(gobsCLs_NLO_qg, 11);
1024     Smooth(gexp1p_NLO_qg, 11);
1025     Smooth(gexp1m_NLO_qg, 11);
1026    
1027     TGraph * g1Sigma_qg = MakeBand(gexp1p_NLO_qg, gexp1m_NLO_qg);
1028     g1Sigma_qg->SetFillColor(5);
1029     h_qg->GetYaxis()->SetTitleOffset(1.5);
1030     h_qg->Draw("");
1031     TGraph * sqLEP = sq_LEP(); sqLEP->Draw("f");
1032     TGraph * glTEV = gl_TEV(); glTEV->Draw("f");
1033     TGraph * glCDF = gl_CDF(); glCDF->Draw("f");
1034     TGraph * glDEZ = gl_DEZ(); glDEZ->Draw("f");
1035    
1036     g1Sigma_qg->Draw("f");
1037     gexpFC_LO_qg->Draw("l");
1038     gexpCLs_LO_qg->Draw("l");
1039     gexpCLs_NLO_qg->Draw("l");
1040     gobsCLs_NLO_qg->Draw("l");
1041     //TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
1042     b.DrawLatex( 800, 50 ,"Tevatron RunI");
1043     b.DrawLatex( 800, 200,"CDF RunII");
1044     b.DrawLatex( 800, 240,"Dzero RunII");
1045 auterman 1.20 b.DrawLatex( 100, 600,"LEP #tilde{q}#tilde{q}");
1046 auterman 1.19 //
1047     TLegend * l_qg = new TLegend(0.45,0.7,0.85,0.89);
1048     l_qg->SetBorderSize(0);l_qg->SetFillColor(0);l_qg->SetTextFont(42);
1049 auterman 1.20 l_qg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1050 auterman 1.19 TGraph * expL_qg = (TGraph*)gexpCLs_NLO_qg->Clone();expL_qg->SetFillStyle(g1Sigma_qg->GetFillStyle());expL_qg->SetFillColor(g1Sigma_qg->GetFillColor());
1051     if (expL_qg) l_qg->AddEntry(expL_qg, "Expected (CLs, NLO)","lf");
1052     if (gobsCLs_NLO_qg) l_qg->AddEntry(gobsCLs_NLO_qg,"Observed (CLs, NLO)","l");
1053     if (gexpCLs_LO_qg) l_qg->AddEntry(gexpCLs_LO_qg, "Expected (CLs, LO)","l");
1054     if (gexpFC_LO_qg) l_qg->AddEntry(gexpFC_LO_qg, "Expected (FC, LO)","l");
1055     l_qg->Draw();
1056     gPad->RedrawAxis();
1057 auterman 1.10 c1->SaveAs("results/Exclusion_mSql_mGl_tb10.pdf");
1058 auterman 1.19
1059 auterman 1.20 //====================================================================================================================
1060     // Exclusion in gluino - squark mass ----------------------------------------
1061     TH2F*h_gq = new TH2F("hExcl_gq",";m_{#tilde{g}} [GeV]; m_{#tilde{q}} [GeV]; 95% CL Observed Exclusion",
1062     50,0,900,60,0,1000);
1063     TH2F*hExcl_gq = new TH2F("h_gq","",50,200,1200,60,200,1400);
1064     TH2F*hexpFC_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
1065     plotToolsHT->Area(hexpFC_HT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
1066     TH2F*hexpFC_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1067     plotToolsHT->Area(hexpFC_HT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
1068     TH2F*hexpFC_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
1069     plotTools->Area(hexpFC_MHT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
1070     TH2F*hexpFC_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1071     plotTools->Area(hexpFC_MHT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
1072     //
1073     TH2F*hexpCLs_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
1074     plotToolsHT->Area(hexpCLs_HT_LO_gq, MGluino, MSquarkL, ExpExclCL);
1075     TH2F*hexpCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1076     plotToolsHT->Area(hexpCLs_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
1077     TH2F*hexpCLs_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
1078     plotTools->Area(hexpCLs_MHT_LO_gq, MGluino, MSquarkL, ExpExclCL);
1079     TH2F*hexpCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1080     plotTools->Area(hexpCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
1081     //
1082     TH2F*hexp1p_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1083     plotToolsHT->Area(hexp1p_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
1084     TH2F*hexp1p_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1085     plotTools->Area(hexp1p_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
1086     TH2F*hexp1m_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1087     plotToolsHT->Area(hexp1m_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
1088     TH2F*hexp1m_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1089     plotTools->Area(hexp1m_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
1090    
1091     TH2F*hobsCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1092     plotToolsHT->Area(hobsCLs_HT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
1093     TH2F*hobsCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1094     plotTools->Area(hobsCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
1095     //
1096     TH2 * hexpFC_LO_gq = plotToolsHT->BinWiseOr(hexpFC_HT_LO_gq, hexpFC_MHT_LO_gq);
1097     TH2 * hexpCLs_LO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_gq, hexpCLs_MHT_LO_gq);
1098     TH2 * hexpCLs_NLO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_gq, hexpCLs_MHT_NLO_gq);
1099     TH2 * hobsCLs_NLO_gq = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_gq, hobsCLs_MHT_NLO_gq);
1100     TH2 * hexp1p_NLO_gq = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_gq, hexp1p_MHT_NLO_gq);
1101     TH2 * hexp1m_NLO_gq = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_gq, hexp1m_MHT_NLO_gq);
1102     //
1103     // TGraph * gexpFC_LO_gq = plotToolsHT->GetContour(hexpFC_LO_gq,3,0);
1104     // gexpFC_LO_gq->SetLineStyle(4); gexpFC_LO_gq->SetLineColor(3);gexpFC_LO_gq->SetLineWidth(2);
1105     TGraph * gexpCLs_LO_gq = plotToolsHT->GetContour(hexpCLs_LO_gq,3,0);
1106     gexpCLs_LO_gq->SetLineStyle(4); gexpCLs_LO_gq->SetLineColor(1);gexpCLs_LO_gq->SetLineWidth(2);
1107     TGraph * gexpCLs_NLO_gq = plotToolsHT->GetContour(hexpCLs_NLO_gq,3,0);
1108     gexpCLs_NLO_gq->SetLineStyle(2); gexpCLs_NLO_gq->SetLineColor(1);gexpCLs_NLO_gq->SetLineWidth(2);
1109     TGraph * gobsCLs_NLO_gq = plotToolsHT->GetContour(hobsCLs_NLO_gq,3,0);
1110     gobsCLs_NLO_gq->SetLineStyle(1); gobsCLs_NLO_gq->SetLineColor(1);gobsCLs_NLO_gq->SetLineWidth(2);
1111     TGraph * gexp1p_NLO_gq = plotToolsHT->GetContour(hexp1p_NLO_gq,3,0);
1112     TGraph * gexp1m_NLO_gq = plotToolsHT->GetContour(hexp1m_NLO_gq,3,0);
1113     TGraph * ch1 = (TGraph*)gobsCLs_NLO_gq->Clone();
1114     // Smooth2D(gexpFC_LO_gq, 11);
1115     Smooth2D(gexpCLs_LO_gq, 11);
1116     Smooth2D(gexpCLs_NLO_gq, 11);
1117     Smooth2D(gobsCLs_NLO_gq, 11);
1118     Smooth2D(gexp1p_NLO_gq, 11);
1119     Smooth2D(gexp1m_NLO_gq, 11);
1120     //Close2D(gobsCLs_NLO_gq);
1121    
1122     TGraph * g1Sigma_gq = MakeBand(gexp1p_NLO_gq, gexp1m_NLO_gq);
1123     g1Sigma_gq->SetFillColor(5);
1124     h_gq->GetYaxis()->SetTitleOffset(1.5);
1125     h_gq->GetXaxis()->SetTitleOffset(0.98);
1126     h_gq->Draw("");
1127     TGraph * sqDEZ = sq_DEZ(); sqDEZ->Draw("f");
1128     TGraph * sqCDF = sq_CDF(); sqCDF->Draw("f");
1129     TGraph * sqTEV = sq_TEV(); sqTEV->Draw("f");
1130     TGraph * glLEP = gl_LEP(); glLEP->Draw("f");
1131     TGraph * noSol = glsq_NoSol(); noSol->Draw("f"); noSol->Draw("l");
1132    
1133     ms.DrawLatex(10,1013,"L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1134     ms.DrawLatex(660,1013,"tan#beta=10, #mu>0, A_{0}=0");
1135     g1Sigma_gq->Draw("f");
1136     //ch1->Draw("l");
1137     // gexpFC_LO_gq->Draw("l");
1138     gexpCLs_LO_gq->Draw("l");
1139     gexpCLs_NLO_gq->Draw("l");
1140     gobsCLs_NLO_gq->Draw("l");
1141     TGraph * expL_gq = (TGraph*)gexpCLs_NLO_gq->Clone();expL_gq->SetFillStyle(g1Sigma_gq->GetFillStyle());expL_gq->SetFillColor(g1Sigma_gq->GetFillColor());
1142     TLegend* legexp_gq = new TLegend(0.2,0.62,0.50,0.93,NULL,"brNDC");
1143     legexp_gq->SetFillColor(0);legexp_gq->SetShadowColor(0);//legexp_gq->SetFillStyle(4000);
1144     legexp_gq->SetTextFont(42);legexp_gq->SetTextSize(0.025);legexp_gq->SetBorderSize(1);
1145     //TEV_sg_cdf.SetLineColor(1);
1146     legexp_gq->SetHeader("CMS preliminary");
1147     if (expL_gq) legexp_gq->AddEntry(expL_gq, "Expected (CLs, NLO)","lf");
1148     if (gobsCLs_NLO_gq) legexp_gq->AddEntry(gobsCLs_NLO_gq,"Observed (CLs, NLO)","l");
1149     if (gexpCLs_LO_gq) legexp_gq->AddEntry(gexpCLs_LO_gq, "Expected (CLs, LO)","l");
1150     legexp_gq->AddEntry(sqTEV,"Tevatron RunI","f");
1151     legexp_gq->AddEntry(sqCDF,"CDF RunII","f");
1152     legexp_gq->AddEntry(sqDEZ,"Dzero RunII","f");
1153     legexp_gq->AddEntry(glLEP,"LEP2","f"); //NOT FOR tb=50!
1154     legexp_gq->Draw();
1155     //TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
1156     //b.DrawLatex( 800, 50,"LEP #tilde{q}#tilde{q}");
1157     //b.DrawLatex( 50, 1000 ,"Tevatron RunI");
1158     //b.DrawLatex( 200, 1100,"CDF RunII");
1159     //b.DrawLatex( 240, 1200,"Dzero RunII");
1160     TPaveText tpave(0.6,0.33,0.75,0.4,"NDC");
1161     tpave.AddText("no cMSSM");tpave.AddText("solution");
1162     tpave.SetBorderSize(0);tpave.SetFillColor(0);tpave.SetTextFont(42);
1163     tpave.Draw();
1164     //
1165     TLegend * l_gq = new TLegend(0.55,0.77,0.90,0.93);
1166     l_gq->SetBorderSize(0);l_gq->SetFillColor(0);l_gq->SetFillStyle(4000);l_gq->SetTextFont(42);
1167     l_gq->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1168     if (expL_gq) l_gq->AddEntry(expL_gq, "Expected (CLs, NLO)","lf");
1169     if (gobsCLs_NLO_gq) l_gq->AddEntry(gobsCLs_NLO_gq,"Observed (CLs, NLO)","l");
1170     if (gexpCLs_LO_gq) l_gq->AddEntry(gexpCLs_LO_gq, "Expected (CLs, LO)","l");
1171     // if (gexpFC_LO_gq) l_gq->AddEntry(gexpFC_LO_gq, "Expected (FC, LO)","l");
1172     //l_gq->Draw();
1173     gPad->RedrawAxis();
1174     c1->SaveAs("results/Exclusion_mGl_mSql_tb10.pdf");
1175     //====================================================================================================================
1176 auterman 1.16
1177     TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
1178 auterman 1.19 60,200,1400,50,200,1000);
1179 auterman 1.18 TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 2,2);
1180     TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,4);
1181     TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 2,1);
1182     TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
1183     TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
1184     TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2);
1185     TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4);
1186     TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1);
1187     TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
1188     TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
1189     Smooth(gexpexcl_qg, 11);
1190     Smooth(gexpexcl_qg_LO, 11);
1191     Smooth(gobsexcl_qg, 11);
1192     Smooth(gCLsExpExclHT_qg, 11);
1193     Smooth(gCLsExpExclHT_qg_LO, 11);
1194     Smooth(gCLsObsExclHT_qg, 11);
1195     Smooth(gCLsExpExclHTm1_qg, 11);
1196     Smooth(gCLsExpExclHTp1_qg, 11);
1197 auterman 1.16 TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
1198     hexcl_qg->Draw("");
1199     //
1200     if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
1201 auterman 1.19 if (gobsexcl_qg) gobsexcl_qg->Draw("l");
1202 auterman 1.16 if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
1203 auterman 1.19 if (gexpexcl_qg) gexpexcl_qg->Draw("l");
1204 auterman 1.16 if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
1205 auterman 1.18 if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l");
1206     if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
1207 auterman 1.19
1208 auterman 1.16 TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
1209     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
1210 auterman 1.18 TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
1211     if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
1212     if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
1213     if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
1214     if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
1215     if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg,"NLO Expected (HT, CLs)","lf");
1216     if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","lf");
1217 auterman 1.16 //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
1218     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
1219     leg_qg->Draw();
1220 auterman 1.19 c1->SaveAs("results/Exclusion_mSql_mGl_tb10_details.pdf");
1221 auterman 1.10 //c1->SaveAs("plot_tb10.pdf");
1222 auterman 1.2
1223 auterman 1.14 c1->SetLogy(1);
1224     c1->SetLogx(1);
1225     TH2F*hUncertScan = new TH2F("hUncertScan",";relative signal uncertainty; number of signal events",
1226     10,0.1,0.6,10,0.9,200);
1227     //hUncertScan->GetXaxis()->SetNdivisions(505, false);
1228     //hUncertScan->SetMinimum(1);
1229     TGraph * gCLsExpUncert = new TGraph(0);
1230     TGraph * gCLsObsUncert = new TGraph(0);
1231     TGraph * gPLExpUncert = new TGraph(0);
1232     TGraph * gPLObsUncert = new TGraph(0);
1233     TGraph * gFCExpUncert = new TGraph(0);
1234     TGraph * gFCObsUncert = new TGraph(0);
1235     TGraph * gMCMCExpUncert= new TGraph(0);
1236     TGraph * gMCMCObsUncert= new TGraph(0);
1237     plotTools->Graph(gCLsExpUncert, SignalRelUncertainty, ExpNSignLimit,1);
1238     plotTools->Graph(gCLsObsUncert, SignalRelUncertainty, ObsNSignLimit,1);
1239     plotTools->Graph(gPLExpUncert, SignalRelUncertainty, PLExpNSignLimit,1);
1240     plotTools->Graph(gPLObsUncert, SignalRelUncertainty, PLObsNSignLimit, 1);
1241     plotTools->Graph(gFCExpUncert, SignalRelUncertainty, FCExpNSignLimit,1);
1242     plotTools->Graph(gFCObsUncert, SignalRelUncertainty, FCObsNSignLimit, 1);
1243     plotTools->Graph(gMCMCExpUncert,SignalRelUncertainty, MCMCExpNSignLimit,1);
1244     plotTools->Graph(gMCMCObsUncert,SignalRelUncertainty, MCMCObsNSignLimit, 1);
1245     gCLsExpUncert->SetLineStyle(8);
1246     gPLObsUncert->SetLineColor(2);
1247     gPLExpUncert->SetLineColor(2); gPLExpUncert->SetLineStyle(8);
1248     gFCObsUncert->SetLineColor(3);
1249     gFCExpUncert->SetLineColor(3); gFCExpUncert->SetLineStyle(8);
1250     gMCMCObsUncert->SetLineColor(4);
1251     gMCMCExpUncert->SetLineColor(4); gMCMCExpUncert->SetLineStyle(8);
1252     hUncertScan->Draw("h");
1253     gCLsExpUncert->Draw("l");
1254     gCLsObsUncert->Draw("l");
1255     //gPLObsUncert->Draw("l");
1256     gPLExpUncert->Draw("l");
1257     //gFCObsUncert->Draw("l");
1258     //gFCExpUncert->Draw("l");
1259     //gMCMCObsUncert->Draw("l");
1260     //gMCMCExpUncert->Draw("l");
1261     TLatex l; l.SetTextSize(0.03); l.SetTextColor(1);
1262     l.DrawLatex(0.025,25,"excluded");
1263     TMarker m; m.SetMarkerStyle(29); m.SetMarkerColor(2);
1264     m.DrawMarker(0.128, 28.5);l.DrawLatex(0.14,28.5,"LM1");
1265     m.DrawMarker(0.147, 77.4);l.DrawLatex(0.16,77.4,"LM0");
1266    
1267     TLegend * leg_S = new TLegend(0.5,0.14,0.85,0.38);
1268     leg_S->SetBorderSize(0);leg_S->SetFillColor(0);
1269     leg_S->SetHeader("95% CL upper limits");
1270     leg_S->AddEntry(gCLsExpUncert,"Expected (CLs, TLimit)","l");
1271     leg_S->AddEntry(gCLsObsUncert,"Observed (CLs, TLimit)","l");
1272     //leg_S->AddEntry(gPLExpUncert,"Expected (PL, RooStat)","l");
1273     //leg_S->AddEntry(gPLObsUncert,"Observed (PL, RooStat)","l");
1274     //leg_S->AddEntry(gFCObsUncert,"Observed (FC, RooStat)","l");
1275     //leg_S->AddEntry(gFCExpUncert,"Expected (FC, RooStat)","l");
1276     //leg_S->AddEntry(gMCMCObsUncert,"Observed (MCMC, RooStat)","l");
1277     //leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
1278     leg_S->Draw();
1279     c1->SaveAs("results/UncertaintyScan.pdf");
1280 auterman 1.2
1281 auterman 1.7 //theApp.Run();
1282 auterman 1.1 }
1283    
1284    
1285     int main(int argc, char** argv)
1286     {
1287     return plot(argc, argv);
1288     }