ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/plot.cc
Revision: 1.18
Committed: Tue Feb 22 13:59:19 2011 UTC (14 years, 2 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.17: +264 -82 lines
Log Message: 
NLO updates

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.1
8     #include "TRint.h"
9     #include "TROOT.h"
10     #include "TObjArray.h"
11     #include "TStyle.h"
12    
13     #include "TChain.h"
14     #include "TFile.h"
15 auterman 1.2 #include "TGraph.h"
16 auterman 1.1 #include "TH1.h"
17     #include "TH2.h"
18     #include "TH2F.h"
19     #include "TTree.h"
20     #include "TKey.h"
21     #include "Riostream.h"
22     #include "TCanvas.h"
23 auterman 1.8 #include "TLegend.h"
24 auterman 1.9 #include "TLatex.h"
25 auterman 1.10 #include "TMarker.h"
26 auterman 1.1
27     #include <string>
28 auterman 1.2 #include <cmath>
29 auterman 1.9 #include <stdio.h>
30 auterman 1.1
31     int plot(int argc, char** argv)
32     {
33 auterman 1.6 //interactive root session
34 auterman 1.7 //TApplication theApp("App", 0, 0);
35 auterman 1.2 if (gROOT->IsBatch()) {
36     fprintf(stderr, "%s: cannot run in batch mode\n", argv[0]);
37     return 1;
38     }
39    
40 auterman 1.6 //Style stuff
41 auterman 1.1 gStyle->SetHistFillColor(0);
42     gStyle->SetPalette(1);
43     gStyle->SetCanvasColor(0);
44     gStyle->SetCanvasBorderMode(0);
45     gStyle->SetPadColor(0);
46     gStyle->SetPadBorderMode(0);
47     gStyle->SetFrameBorderMode(0);
48    
49     gStyle->SetTitleFillColor(0);
50     gStyle->SetTitleBorderSize(0);
51     gStyle->SetTitleX(0.10);
52     gStyle->SetTitleY(0.98);
53     gStyle->SetTitleW(0.8);
54     gStyle->SetTitleH(0.06);
55    
56     gStyle->SetErrorX(0);
57     gStyle->SetStatColor(0);
58     gStyle->SetStatBorderSize(0);
59     gStyle->SetStatX(0);
60     gStyle->SetStatY(0);
61     gStyle->SetStatW(0);
62     gStyle->SetStatH(0);
63    
64     gStyle->SetTitleFont(22);
65     gStyle->SetLabelFont(22,"X");
66     gStyle->SetLabelFont(22,"Y");
67     gStyle->SetLabelFont(22,"Z");
68     gStyle->SetLabelSize(0.03,"X");
69     gStyle->SetLabelSize(0.03,"Y");
70     gStyle->SetLabelSize(0.03,"Z");
71    
72     //gROOT->SetStyle("MyStyle");
73     TCanvas * c1 = new TCanvas("c1","c1",600,600);
74     c1->SetFillStyle ( 4000 );
75     c1->SetLeftMargin ( 0.15 );
76     c1->SetRightMargin ( 0.15 );
77     c1->SetBottomMargin( 0.10 );
78     c1->cd();
79    
80 auterman 1.6 //Get limits and generator masses ---------------------------------------------------
81 auterman 1.4 TheLimits * genpoints = new TheLimits();
82 auterman 1.11 //genpoints->Fill(argc, argv);
83 auterman 1.18 genpoints->Fill("limits_MHT_tb10/filelist.txt");
84 auterman 1.11
85     TheLimits * genpointsHT = new TheLimits();
86 auterman 1.18 genpointsHT->Fill("limits_HT_tb10/filelist.txt");
87 auterman 1.12
88 auterman 1.18 //TheLimits * genpointsHT = new TheLimits();
89     //genpointsHT->Fill("limits_HT_SigCont/filelist.txt");
90 auterman 1.13
91 auterman 1.18 //TheLimits * genpointsMHT = new TheLimits();
92     //genpointsMHT->Fill("limits_MHT_SigCont/filelist.txt");
93 auterman 1.11
94 auterman 1.10
95     //Replace read limits with specific numbers
96     //genpoints->OverwriteLimits("ABCD_MHT");
97    
98 auterman 1.17 genpoints->FillGeneratorMasses("GenScan_tb10.dat");
99 auterman 1.5 genpoints->match();
100 auterman 1.17 genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
101 auterman 1.11 genpointsHT->match();
102 auterman 1.18 //genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
103     //genpointsHT->match();
104     //genpointsMHT->FillGeneratorMasses("GenScan_tb10.dat");
105     //genpointsMHT->match();
106 auterman 1.2
107 auterman 1.6 //the plotting ----------------------------------------------------------------------
108 auterman 1.7 //plotting helper functions
109 auterman 1.4 PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
110 auterman 1.11 PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
111 auterman 1.18 //PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
112     //PlotTools<SusyScan> * plotToolsMHT = new PlotTools<SusyScan>(genpointsMHT->GetScan());
113 auterman 1.6 PlotTools<GeneratorMasses> * plotMasses = new PlotTools<GeneratorMasses>(genpoints->GetGeneratorMasses());
114 auterman 1.2
115 auterman 1.7 //iso mass lines
116     TGraph * gl500 = plotMasses->Line(Mzero, Mhalf, MGluino, 500);
117     TGraph * sq500 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 500, 10);
118     sq500->SetLineWidth(2); sq500->SetLineColor(7);
119    
120     TGraph * chi100 = plotMasses->Line(Mzero, Mhalf, MChi1, 50, 20);
121     TGraph * cha200 = plotMasses->Line(Mzero, Mhalf, MCha1, 200, 20);
122     cha200->SetLineColor(2);
123    
124     //the histograms
125     c1->SetLogz(1);
126     //h->SetMaximum(27);
127     //h->SetMinimum(0.01);
128 auterman 1.10
129     /**/
130 auterman 1.7 // cross-section in M0 - M1/2
131     TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
132 auterman 1.18 100,0,1009.9,35,100,450);
133 auterman 1.7 plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
134     hxsec->SetMinimum(0.01);
135     //sq500->Draw();
136     //gl500->Draw();
137     hxsec->Draw("colz");
138 auterman 1.10 c1->SaveAs("results/Xsection_m0_m12_tb10.pdf");
139 auterman 1.7 std::string wait;
140     //std::cin>>wait;
141 auterman 1.10
142 auterman 1.7 // Observed Limit in M0 - M1/2
143     TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
144 auterman 1.18 100,0,1009.9,35,100,450);
145 auterman 1.7 plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
146     hobslimit->SetMinimum(0.01);
147     hobslimit->Draw("colz");
148 auterman 1.14 c1->SaveAs("results/ObsLimit_m0_m12_tb10.pdf");
149    
150 auterman 1.15
151 auterman 1.7 // Expected Limit in M0 - M1/2
152     TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
153 auterman 1.18 100,0,1009.9,35,100,450);
154 auterman 1.7 plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
155     hexplimit->SetMinimum(0.01);
156     hexplimit->Draw("colz");
157 auterman 1.10 c1->SaveAs("results/ExpLimit_m0_m12_tb10.pdf");
158 auterman 1.7
159 auterman 1.13 // Signal Acceptance in M0 - M1/2
160 auterman 1.7 TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
161 auterman 1.18 100,0,1009.9,35,100,450);
162 auterman 1.7 plotTools->Area(hsigacc, Mzero, Mhalf, SignalAcceptance);
163     hsigacc->SetMinimum(0.01);
164     hsigacc->SetMaximum(1.0);
165     hsigacc->Draw("colz");
166 auterman 1.14 //chi100->Draw();
167     //cha200->Draw();
168     //gl500 ->Draw();
169     // sq500 ->Draw();
170     c1->SaveAs("results/SigAccMHT_m0_m12_tb10.pdf");
171 auterman 1.18
172 auterman 1.14
173 auterman 1.7 // Exp. Limit on Number of Signal Events in M0 - M1/2
174 auterman 1.8 c1->SetLogz(0);
175 auterman 1.7 TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
176 auterman 1.18 100,0,1009.9,35,100,450);
177 auterman 1.7 plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
178 auterman 1.10 hexplimitnsig->SetMinimum(0.0);
179     hexplimitnsig->SetMaximum(20);
180 auterman 1.7 hexplimitnsig->Draw("colz");
181 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_m0_m12_tb10.pdf");
182 auterman 1.7
183     // Obs. Limit on Number of Signal Events in M0 - M1/2
184     TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
185 auterman 1.18 100,0,1009.9,35,100,450);
186 auterman 1.7 plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
187 auterman 1.10 hobslimitnsig->SetMinimum(0.0);
188     hobslimitnsig->SetMaximum(20);
189 auterman 1.7 hobslimitnsig->Draw("colz");
190 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_m0_m12_tb10.pdf");
191 auterman 1.7
192     c1->SetLogz(0);
193     // Expected Exclusion in M0 - M1/2
194     TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
195 auterman 1.8 100,0,1009.9,35,100,450);
196 auterman 1.11 plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
197 auterman 1.7 hexpexcl->Draw("colz");
198 auterman 1.10 c1->SaveAs("results/ExpExclusion_m0_m12_tb10.pdf");
199 auterman 1.7
200     // Observed Exclusion in M0 - M1/2
201     TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
202 auterman 1.8 100,0,1009.9,35,100,450);
203 auterman 1.11 plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
204 auterman 1.7 hobsexcl->Draw("colz");
205 auterman 1.10 c1->SaveAs("results/ObsExclusion_m0_m12_tb10.pdf");
206 auterman 1.7
207 auterman 1.8 // Observed Exclusion in M0 - M1/2
208     TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
209     100,0,1009.9,35,100,450);
210     plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
211     hPLobsexcl->Draw("colz");
212 auterman 1.10 c1->SaveAs("results/PL_ObsExclusion_m0_m12_tb10.pdf");
213 auterman 1.8
214 auterman 1.9 // TestContours in M0 - M1/2
215 auterman 1.14 TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
216     100,0,1009.9,35,100,450);
217     TH2F*tPLexpexcl=(TH2F*)texcl->Clone();
218     plotTools->Area(tPLexpexcl, Mzero, Mhalf, PLExpExclusion);
219     TH2F*tFCexpexcl=(TH2F*)texcl->Clone();
220     plotTools->Area(tFCexpexcl, Mzero, Mhalf, FCExpExclusion);
221 auterman 1.18 TH2F*testExcl=(TH2F*)texcl->Clone();
222     plotTools->Area(testExcl, Mzero, Mhalf, NLOMCMCExpExclusion);
223     //testExcl->SetMinimum(0.9);
224     //testExcl->SetMaximum(1.1);
225     std::vector<TGraph*> contours = plotTools->GetContours(testExcl,3);
226 auterman 1.14 //std::vector<TGraph*> contours = plotTools->GetContours(hobsexcl,3);
227     //std::vector<TGraph*> contours = plotTools->GetContours(tPLexpexcl,3);
228     //std::vector<TGraph*> contours = plotTools->GetContours(hPLobsexcl,3);
229 auterman 1.9 //hPLexpexcl
230     //hexcl->Draw("colz");
231 auterman 1.18 testExcl->Draw("colz");
232 auterman 1.9 //hobsexcl->Draw("colz");
233 auterman 1.14 int col=kBlue-10;
234 auterman 1.9 for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
235     if (! *cont) continue;
236     double x, y;
237     (*cont)->GetPoint(0, x, y);
238     (*cont)->SetLineColor(col);
239     (*cont)->Draw("l");
240     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col++);
241     char val[20];
242 auterman 1.13 sprintf(val,"%d",(int)(cont-contours.begin()));
243 auterman 1.9 l.DrawLatex(x,y,val);
244 auterman 1.11 //if (cont-contours.begin()>3) break;
245 auterman 1.9 }
246 auterman 1.10 c1->SaveAs("results/ExclusionTestContours_m0_m12_tb10.pdf");
247 auterman 1.9
248 auterman 1.13
249 auterman 1.8 // Exclusion in M0 - M1/2
250     TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
251     100,0,1009.9,35,100,450);
252 auterman 1.14 hexcl->Draw("colz");
253 auterman 1.18 TGraph * gexpexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
254     TGraph * gexpexcl_LO = plotTools ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
255     TGraph * gobsexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
256     TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2);
257     TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1);
258     TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
259     TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
260     TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
261     TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,ObsExclCL, 3,0, 1,1);
262     TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
263     TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
264     Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 );
265     Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
266     Smooth( gCLsObsExclHT_LO, 33); gCLsObsExclHT_LO->SetLineWidth( 3 );
267     Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
268     Smooth( gexpexcl, 33 );
269     Smooth( gexpexcl_LO, 33 );
270     Smooth( gMCMCexpexcl, 33 );
271     Smooth( gCLsExpExclHTm1, 33 );
272     Smooth( gCLsExpExclHTp1, 33 );
273 auterman 1.13 TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1);
274     hexcl->Draw("colz");
275 auterman 1.10 //set old exclusion Limits
276 auterman 1.16 TGraph* LEP_ch = set_lep_ch(10);
277     TGraph* LEP_sl = set_lep_sl(10);//slepton curve
278     TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
279     TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
280     //TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
281     //TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
282     TGraph* stau = set_tev_stau(10);//stau
283     TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
284     TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
285 auterman 1.15 TGraphErrors* First = getObserved_NLO_tanBeta10();
286     TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
287     TGraphErrors* Third = getLO_tanBeta10();
288 auterman 1.10 First->GetXaxis()->SetRangeUser(0,505);
289     First->GetYaxis()->SetRangeUser(80,500);
290     First->GetXaxis()->SetTitle("m_{0} (GeV)");
291     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
292     TSpline3 *sFirst = new TSpline3("sFirst",First);
293     sFirst->SetLineColor(kBlue);
294 auterman 1.18 sFirst->SetLineWidth(1);
295 auterman 1.10 TSpline3 *sSecond = new TSpline3("sSecond",Second);
296     sSecond->SetLineColor(kBlue);
297     sSecond->SetLineStyle(2);
298 auterman 1.18 sSecond->SetLineWidth(1);
299 auterman 1.10 TSpline3 *sThird = new TSpline3("sThird",Third);
300 auterman 1.18 sThird->SetLineColor(kBlue);
301 auterman 1.10 sThird->SetLineStyle(4);
302 auterman 1.18 sThird->SetLineWidth(1);
303 auterman 1.10 TEV_sn_d0_1->Draw("fsame");
304     TEV_sn_d0_2->Draw("fsame");
305 auterman 1.17 TEV_sg_d0->Draw("fsame");
306     TEV_sg_cdf->Draw("fsame");
307 auterman 1.10 LEP_ch->Draw("fsame");
308     LEP_sl->Draw("fsame");
309     stau->Draw("fsame");
310     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
311     b.DrawLatex( 10,420,"#tilde{#tau} LSP");
312     b.DrawLatex( 6,150,"LEP2");
313     b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
314     b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
315 auterman 1.17 b.DrawLatex(90,160,"CDF");
316     //b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
317 auterman 1.10 b.DrawLatex( 80,180,"D0 #tilde{#nu}");
318 auterman 1.18 if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
319     sFirst->Draw("same");
320     sSecond->Draw("same");
321     sThird->Draw("same");
322     if (gexpexcl) gexpexcl->Draw("l");
323     //if (gexpexcl_LO) gexpexcl_LO->Draw("l");
324     if (gMCMCexpexcl) gMCMCexpexcl->Draw("l");
325     //if (gCLsObsExclHT_LO) gCLsObsExclHT_LO->Draw("l");
326     if (gCLsExpExclHT_LO) gCLsExpExclHT_LO->Draw("l");
327     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
328     if (gCLsExpExclHT) gCLsExpExclHT->Draw("l");
329     //if (gobsexcl) gobsexcl->Draw("l");
330     //if (gPLobsexcl) gPLobsexcl->Draw("l");
331     //if (gFCobsexcl) gFCobsexcl->Draw("l");
332     //if (gMCMCobsexcl) gMCMCobsexcl->Draw("l");
333     //if (gFCexpexcl) gFCexpexcl->Draw("l");
334     //if (gPLexpexcl) gPLexpexcl->Draw("l");
335     //if (gFCexpexcl) gFCexpexcl->Draw("l");
336 auterman 1.14 TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
337 auterman 1.8 leg->SetBorderSize(0);leg->SetFillColor(0);
338 auterman 1.18 TGraph * expLeg = (TGraph*)gCLsExpExclHT->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
339 auterman 1.12 //if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
340 auterman 1.18 if (gCLsExpExclHT) leg->AddEntry(expLeg ,"NLO Expected (HT, CLs)","l");
341     if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
342     //if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
343     if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
344 auterman 1.11 //if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
345     //if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
346 auterman 1.10 //if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
347     //if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
348 auterman 1.12 //if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
349     if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
350 auterman 1.18 if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
351     if (sSecond) leg->AddEntry(sSecond,"NLO #alpha_{T} Expected");
352     if (sFirst) leg->AddEntry(sFirst, "NLO #alpha_{T} Observed");
353     if (sThird) leg->AddEntry(sThird,"LO #alpha_{T} Expected");
354 auterman 1.8 leg->Draw();
355 auterman 1.10 gPad->RedrawAxis();
356 auterman 1.13 c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
357 auterman 1.15
358 auterman 1.8
359 auterman 1.13 // Signal Contamination in M0 - M1/2
360     c1->SetLogz(0);
361     TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
362 auterman 1.18 100,0,1009.9,35,100,450);
363     plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
364 auterman 1.13 hsigcont->SetMinimum(0.01);
365     hsigcont->SetMaximum(20);
366     hsigcont->Draw("colz");
367     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
368     c1->SaveAs("results/SignalContamination_HT_m0_m12_tb10.pdf");
369    
370     // Signal Contamination in M0 - M1/2
371     c1->SetLogz(0);
372     TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
373 auterman 1.18 100,0,1009.9,35,100,450);
374     plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
375 auterman 1.13 hsigcontMHT->SetMinimum(0.01);
376     hsigcontMHT->SetMaximum(20);
377     hsigcontMHT->Draw("colz");
378     if (gobsexcl) gobsexcl->Draw("l");
379     c1->SaveAs("results/SignalContamination_MHT_m0_m12_tb10.pdf");
380    
381 auterman 1.18 // Signal JEC Uncertainty MHT
382     c1->SetLogz(0);
383     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
384     100,0,1009.9,35,100,450);
385     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
386     hsig_jec_mht->SetMinimum(0.0);
387     hsig_jec_mht->SetMaximum(0.25);
388     hsig_jec_mht->SetContour(10);
389     hsig_jec_mht->Draw("colz");
390     if (gobsexcl) gobsexcl->Draw("l");
391     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
392    
393     // Signal JEC Uncertainty HT
394     c1->SetLogz(0);
395     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
396     100,0,1009.9,35,100,450);
397     plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
398     hsig_jec_ht->SetMinimum(0.0);
399     hsig_jec_ht->SetMaximum(0.25);
400     hsig_jec_ht->SetContour(10);
401     hsig_jec_ht->Draw("colz");
402     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
403     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
404    
405    
406     // Signal MuIso Uncertainty MHT
407     c1->SetLogz(0);
408     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
409     100,0,1009.9,35,100,450);
410     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
411     hsig_MuIso_mht->SetMinimum(0.0);
412     hsig_MuIso_mht->SetMaximum(0.05);
413     hsig_MuIso_mht->SetContour(5);
414     hsig_MuIso_mht->Draw("colz");
415     if (gobsexcl) gobsexcl->Draw("l");
416     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
417    
418     // Signal MuIso Uncertainty HT
419     c1->SetLogz(0);
420     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
421     100,0,1009.9,35,100,450);
422     plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
423     hsig_MuIso_ht->SetMinimum(0.0);
424     hsig_MuIso_ht->SetMaximum(0.05);
425     hsig_MuIso_ht->SetContour(5);
426     hsig_MuIso_ht->Draw("colz");
427     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
428     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
429    
430    
431     // Signal kFactor Uncertainty MHT
432     c1->SetLogz(0);
433     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
434     100,0,1009.9,35,100,450);
435     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
436     hsig_kFactor_mht->SetMinimum(0.1);hsig_kFactor_mht->GetZaxis()->SetTitleOffset(1.3);
437     hsig_kFactor_mht->SetMaximum(0.2);
438     hsig_kFactor_mht->SetContour(10);
439     hsig_kFactor_mht->Draw("colz");
440     if (gobsexcl) gobsexcl->Draw("l");
441     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
442    
443     // Signal kFactor Uncertainty HT
444     c1->SetLogz(0);
445     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
446     100,0,1009.9,35,100,450);
447     plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
448     hsig_kFactor_ht->SetMinimum(0.1);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
449     hsig_kFactor_ht->SetMaximum(0.2);
450     hsig_kFactor_ht->SetContour(10);
451     hsig_kFactor_ht->Draw("colz");
452     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
453     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
454    
455 auterman 1.1
456 auterman 1.2
457 auterman 1.7 //-----------------------------------------------------------------------------------
458     c1->SetLogz(1);
459    
460     TGraph * mz500 = plotMasses->Line( MSquarkL, MGluino, Mzero, 500, 1);
461     TGraph * mh250 = plotMasses->Line( MSquarkL, MGluino, Mhalf, 250, 1);
462    
463 auterman 1.8 // cross-section in squark - gluino mass
464 auterman 1.7 TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
465     60,200,1400,50,200,1200);
466     plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
467     hxsec_qg->SetMinimum(0.01);
468     hxsec_qg->Draw("colz");
469     mz500->Draw();
470     mh250->Draw();
471 auterman 1.10 c1->SaveAs("results/Xsection_mSql_mGl_tb10.pdf");
472 auterman 1.7
473 auterman 1.8 // Observed Limit in squark - gluino mass
474 auterman 1.7 TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
475     60,200,1400,50,200,1200);
476     plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
477     hobslimit_qg->SetMinimum(0.01);
478     hobslimit_qg->Draw("colz");
479 auterman 1.10 c1->SaveAs("results/ObsLimit_mSql_mGl_tb10.pdf");
480 auterman 1.7
481 auterman 1.8 // Expected Limit in squark - gluino mass
482 auterman 1.7 TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
483     60,200,1400,50,200,1200);
484     plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
485     hexplimit_qg->SetMinimum(0.01);
486     hexplimit_qg->Draw("colz");
487 auterman 1.10 c1->SaveAs("results/ExpLimit_mSql_mGl_tb10.pdf");
488 auterman 1.7
489 auterman 1.8 // Signal Acceptance in squark - gluino mass
490 auterman 1.7 TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
491     60,200,1400,50,200,1200);
492     plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
493     hsigacc_qg->SetMinimum(0.01);
494     hsigacc_qg->SetMaximum(1.0);
495     hsigacc_qg->Draw("colz");
496 auterman 1.14 c1->SaveAs("results/SigAccMHT_mSql_mGl_tb10.pdf");
497    
498     // Signal Acceptance in squark - gluino mass
499     TH2F*hsigaccHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
500     60,200,1400,50,200,1200);
501     plotToolsHT->Area(hsigaccHT_qg, MSquarkL, MGluino, SignalAcceptance);
502     hsigaccHT_qg->SetMinimum(0.01);
503     hsigaccHT_qg->SetMaximum(1.0);
504     hsigaccHT_qg->Draw("colz");
505     c1->SaveAs("results/SigAccHT_mSql_mGl_tb10.pdf");
506 auterman 1.18 /*
507     // Signal JEC Uncertainty MHT
508     c1->SetLogz(0);
509     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
510     100,0,1009.9,35,100,450);
511     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
512     hsig_jec_mht->SetMinimum(0.0);
513     hsig_jec_mht->SetMaximum(0.25);
514     hsig_jec_mht->SetContour(10);
515     hsig_jec_mht->Draw("colz");
516     if (gobsexcl) gobsexcl->Draw("l");
517     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
518    
519     // Signal JEC Uncertainty HT
520     c1->SetLogz(0);
521     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
522     100,0,1009.9,35,100,450);
523     plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
524     hsig_jec_ht->SetMinimum(0.0);
525     hsig_jec_ht->SetMaximum(0.25);
526     hsig_jec_ht->SetContour(10);
527     hsig_jec_ht->Draw("colz");
528     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
529     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
530    
531 auterman 1.15
532 auterman 1.18 // Signal MuIso Uncertainty MHT
533     c1->SetLogz(0);
534     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
535     100,0,1009.9,35,100,450);
536     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
537     hsig_MuIso_mht->SetMinimum(0.0);
538     hsig_MuIso_mht->SetMaximum(0.05);
539     hsig_MuIso_mht->SetContour(5);
540     hsig_MuIso_mht->Draw("colz");
541     if (gobsexcl) gobsexcl->Draw("l");
542     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
543    
544     // Signal MuIso Uncertainty HT
545     c1->SetLogz(0);
546     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
547     100,0,1009.9,35,100,450);
548     plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
549     hsig_MuIso_ht->SetMinimum(0.0);
550     hsig_MuIso_ht->SetMaximum(0.05);
551     hsig_MuIso_ht->SetContour(5);
552     hsig_MuIso_ht->Draw("colz");
553     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
554     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
555    
556    
557     // Signal kFactor Uncertainty MHT
558     c1->SetLogz(0);
559     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
560     100,0,1009.9,35,100,450);
561     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
562     hsig_kFactor_mht->SetMinimum(0.1);
563     hsig_kFactor_mht->SetMaximum(0.2);
564     hsig_kFactor_mht->SetContour(10);
565     hsig_kFactor_mht->Draw("colz");
566     if (gobsexcl) gobsexcl->Draw("l");
567     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
568    
569     // Signal kFactor Uncertainty HT
570     c1->SetLogz(0);
571     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
572     100,0,1009.9,35,100,450);
573     plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
574     hsig_kFactor_ht->SetMinimum(0.1);
575     hsig_kFactor_ht->SetMaximum(0.2);
576     hsig_kFactor_ht->SetContour(10);
577     hsig_kFactor_ht->Draw("colz");
578     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
579     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
580 auterman 1.7
581 auterman 1.18 */
582 auterman 1.8 // Exp. Limit on Number of Signal Events in squark - gluino mass
583     c1->SetLogz(0);
584 auterman 1.7 TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
585     60,200,1400,50,200,1200);
586     plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
587 auterman 1.10 hexplimitnsig_qg->SetMinimum(0.0);
588     hexplimitnsig_qg->SetMaximum(20);
589 auterman 1.7 hexplimitnsig_qg->Draw("colz");
590 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_mSql_mGl_tb10.pdf");
591 auterman 1.6
592 auterman 1.8 // Obs. Limit on Number of Signal Events in squark - gluino mass
593 auterman 1.7 TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
594     60,200,1400,50,200,1200);
595     plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
596 auterman 1.10 hobslimitnsig_qg->SetMinimum(0.0);
597     hobslimitnsig_qg->SetMaximum(20);
598 auterman 1.7 hobslimitnsig_qg->Draw("colz");
599 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_mSql_mGl_tb10.pdf");
600 auterman 1.7
601     c1->SetLogz(0);
602 auterman 1.8 // Expected Exclusion in squark - gluino mass
603 auterman 1.7 TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
604 auterman 1.2 60,200,1400,50,200,1200);
605 auterman 1.11 plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
606 auterman 1.7 hexpexcl_qg->Draw("colz");
607 auterman 1.10 c1->SaveAs("results/ExpExclusion_mSql_mGl_tb10.pdf");
608 auterman 1.2
609 auterman 1.8 // Observed Exclusion in squark - gluino mass
610 auterman 1.7 TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
611     60,200,1400,50,200,1200);
612 auterman 1.11 plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
613 auterman 1.8 //std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
614 auterman 1.7 hobsexcl_qg->Draw("colz");
615 auterman 1.8 //for (std::vector<TGraph *>::iterator g=vobsexcl_qg.begin();g!=vobsexcl_qg.end();++g){
616     // if (! *g) continue;
617     // // (*g)->Draw("l");
618     //}
619 auterman 1.10 c1->SaveAs("results/ObsExclusion_mSql_mGl_tb10.pdf");
620 auterman 1.2
621 auterman 1.9 // TestContours in M0 - M1/2
622     TH2F*texcl_qg = new TH2F("texcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
623     60,200,1400,50,200,1200);
624     TH2F*tPLobsexcl_qg=(TH2F*)texcl_qg->Clone();
625     plotTools->Area(tPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
626 auterman 1.10 TH2F*tPLexpexcl_qg=(TH2F*)texcl_qg->Clone();
627     plotTools->Area(tPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
628 auterman 1.9 //std::vector<TGraph*> contours_qg = plotTools->GetContours(hexpexcl_qg,3);
629     //std::vector<TGraph*> contours_qg = plotTools->GetContours(hobsexcl_qg,3);
630     std::vector<TGraph*> contours_qg = plotTools->GetContours(tPLobsexcl_qg,3);
631 auterman 1.10 //std::vector<TGraph*> contours_qg = plotTools->GetContours(tPLexpexcl_qg,3);
632 auterman 1.9 //hPLexpexcl_qg
633     //hexcl_qg->Draw("colz");
634     //hexpexcl_qg->Draw("colz");
635     //hobsexcl_qg->Draw("colz");
636     tPLobsexcl_qg->Draw("colz");
637     int col_gl=kBlue-10;
638     for (std::vector<TGraph*>::iterator cont=contours_qg.begin(); cont!=contours_qg.end(); ++cont){
639     if (! *cont) continue;
640     double x, y;
641     (*cont)->GetPoint(0, x, y);
642     (*cont)->SetLineColor(col_gl);(*cont)->SetLineWidth(3);
643     (*cont)->Draw("l");
644     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col_gl++);
645     char val[20];
646 auterman 1.13 sprintf(val,"%d",(int)(cont-contours_qg.begin()));
647 auterman 1.9 l.DrawLatex(x,y,val);
648 auterman 1.10 if (cont-contours_qg.begin()>3) break;
649 auterman 1.9 }
650 auterman 1.10 c1->SaveAs("results/ExclusionTestContours_mSql_mGl_tb10.pdf");
651 auterman 1.9
652 auterman 1.8 // Observed Exclusion in squark - gluino mass
653     TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
654     60,200,1400,50,200,1200);
655     plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
656 auterman 1.10 //hPLobsexcl_qg->Draw("colz");
657     //c1->SaveAs("results/PL_ObsExclusion_mSql_mGl_tb10.pdf");
658 auterman 1.8
659 auterman 1.11
660    
661 auterman 1.16 /*
662 auterman 1.10 // Exclusion in squark - gluino mass ----------------------------------------
663 auterman 1.8 TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
664     60,200,1400,50,200,1200);
665 auterman 1.18 TGraph * gexpexcl_qg = plotToolsMHT->GetContour(hexpexcl_qg,3,0);
666 auterman 1.8 if (gexpexcl_qg) gexpexcl_qg->SetLineStyle(2);
667 auterman 1.18 TGraph * gobsexcl_qg = plotToolsMHT->GetContour(hobsexcl_qg,3,0);
668 auterman 1.8 TH2F*hPLexpexcl_qg=(TH2F*)hexcl_qg->Clone();
669 auterman 1.18 plotToolsMHT->Area(hPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
670     TGraph * gPLexpexcl_qg = plotToolsMHT->GetContour(hPLexpexcl_qg,3,0);
671 auterman 1.10 if (gPLexpexcl_qg) {gPLexpexcl_qg->SetLineStyle(2);gPLexpexcl_qg->SetLineColor(2);}
672 auterman 1.18 TGraph * gPLobsexcl_qg = plotToolsMHT->GetContour(hPLobsexcl_qg,3,0);
673 auterman 1.10 if (gPLobsexcl_qg) gPLobsexcl_qg->SetLineColor(2);
674 auterman 1.11 TH2F*hCLsExpHT_qg=(TH2F*)hexcl_qg->Clone();
675     TH2F*hCLsObsHT_qg=(TH2F*)hexcl_qg->Clone();
676 auterman 1.18 plotToolsHT->Area(hCLsExpHT_qg, MSquarkL, MGluino, ExpExclCL);
677     plotToolsHT->Area(hCLsObsHT_qg, MSquarkL, MGluino, ObsExclCL);
678     TGraph * gCLsExpExclHT_qg = plotToolsMHT->GetContour(hCLsExpHT_qg,3,0);
679     TGraph * gCLsObsExclHT_qg = plotToolsMHT->GetContour(hCLsObsHT_qg,3,0);
680 auterman 1.11 if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->SetLineColor(2);
681     if (gCLsExpExclHT_qg) {gCLsExpExclHT_qg->SetLineStyle(2);gCLsExpExclHT_qg->SetLineColor(2);}
682 auterman 1.10 hexcl_qg->Draw("");
683 auterman 1.9 if (gobsexcl_qg) gobsexcl_qg->Draw("l");
684 auterman 1.11 if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
685 auterman 1.9 if (gexpexcl_qg) gexpexcl_qg->Draw("l");
686 auterman 1.11 //if (gPLobsexcl_qg) gPLobsexcl_qg->Draw("l");
687     //if (gPLexpexcl_qg) gPLexpexcl_qg->Draw("l");
688     if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
689 auterman 1.8 TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
690     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
691 auterman 1.11 if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"Observed (MHT, CLs, TLimit)","l");
692     if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"Expected (MHT, CLs, TLimit)","l");
693     if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"Observed (HT, CLs, TLimit)","l");
694     if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"Expected (HT, CLs, TLimit)","l");
695     //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
696     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
697 auterman 1.8 leg_qg->Draw();
698 auterman 1.10 c1->SaveAs("results/Exclusion_mSql_mGl_tb10.pdf");
699 auterman 1.16 */
700    
701     TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
702     60,200,1400,50,200,1200);
703 auterman 1.18 TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 2,2);
704     TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,4);
705     TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 2,1);
706     TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
707     TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
708     TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2);
709     TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4);
710     TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1);
711     TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
712     TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
713     Smooth(gexpexcl_qg, 11);
714     Smooth(gexpexcl_qg_LO, 11);
715     Smooth(gobsexcl_qg, 11);
716     Smooth(gCLsExpExclHT_qg, 11);
717     Smooth(gCLsExpExclHT_qg_LO, 11);
718     Smooth(gCLsObsExclHT_qg, 11);
719     Smooth(gCLsExpExclHTm1_qg, 11);
720     Smooth(gCLsExpExclHTp1_qg, 11);
721 auterman 1.16 TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
722     hexcl_qg->Draw("");
723     //
724     if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
725     if (gobsexcl_qg) gobsexcl_qg->Draw("l");
726     if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
727     if (gexpexcl_qg) gexpexcl_qg->Draw("l");
728     if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
729 auterman 1.18 if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l");
730     if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
731 auterman 1.16 //
732     TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
733     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
734 auterman 1.18 TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
735     if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
736     if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
737     if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
738     if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
739     if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg,"NLO Expected (HT, CLs)","lf");
740     if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","lf");
741 auterman 1.16 //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
742     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
743     leg_qg->Draw();
744     c1->SaveAs("results/Exclusion_mSql_mGl_tb10.pdf");
745 auterman 1.10 //c1->SaveAs("plot_tb10.pdf");
746 auterman 1.2
747 auterman 1.14 c1->SetLogy(1);
748     c1->SetLogx(1);
749     TH2F*hUncertScan = new TH2F("hUncertScan",";relative signal uncertainty; number of signal events",
750     10,0.1,0.6,10,0.9,200);
751     //hUncertScan->GetXaxis()->SetNdivisions(505, false);
752     //hUncertScan->SetMinimum(1);
753     TGraph * gCLsExpUncert = new TGraph(0);
754     TGraph * gCLsObsUncert = new TGraph(0);
755     TGraph * gPLExpUncert = new TGraph(0);
756     TGraph * gPLObsUncert = new TGraph(0);
757     TGraph * gFCExpUncert = new TGraph(0);
758     TGraph * gFCObsUncert = new TGraph(0);
759     TGraph * gMCMCExpUncert= new TGraph(0);
760     TGraph * gMCMCObsUncert= new TGraph(0);
761     plotTools->Graph(gCLsExpUncert, SignalRelUncertainty, ExpNSignLimit,1);
762     plotTools->Graph(gCLsObsUncert, SignalRelUncertainty, ObsNSignLimit,1);
763     plotTools->Graph(gPLExpUncert, SignalRelUncertainty, PLExpNSignLimit,1);
764     plotTools->Graph(gPLObsUncert, SignalRelUncertainty, PLObsNSignLimit, 1);
765     plotTools->Graph(gFCExpUncert, SignalRelUncertainty, FCExpNSignLimit,1);
766     plotTools->Graph(gFCObsUncert, SignalRelUncertainty, FCObsNSignLimit, 1);
767     plotTools->Graph(gMCMCExpUncert,SignalRelUncertainty, MCMCExpNSignLimit,1);
768     plotTools->Graph(gMCMCObsUncert,SignalRelUncertainty, MCMCObsNSignLimit, 1);
769     gCLsExpUncert->SetLineStyle(8);
770     gPLObsUncert->SetLineColor(2);
771     gPLExpUncert->SetLineColor(2); gPLExpUncert->SetLineStyle(8);
772     gFCObsUncert->SetLineColor(3);
773     gFCExpUncert->SetLineColor(3); gFCExpUncert->SetLineStyle(8);
774     gMCMCObsUncert->SetLineColor(4);
775     gMCMCExpUncert->SetLineColor(4); gMCMCExpUncert->SetLineStyle(8);
776     hUncertScan->Draw("h");
777     gCLsExpUncert->Draw("l");
778     gCLsObsUncert->Draw("l");
779     //gPLObsUncert->Draw("l");
780     gPLExpUncert->Draw("l");
781     //gFCObsUncert->Draw("l");
782     //gFCExpUncert->Draw("l");
783     //gMCMCObsUncert->Draw("l");
784     //gMCMCExpUncert->Draw("l");
785     TLatex l; l.SetTextSize(0.03); l.SetTextColor(1);
786     l.DrawLatex(0.025,25,"excluded");
787     TMarker m; m.SetMarkerStyle(29); m.SetMarkerColor(2);
788     m.DrawMarker(0.128, 28.5);l.DrawLatex(0.14,28.5,"LM1");
789     m.DrawMarker(0.147, 77.4);l.DrawLatex(0.16,77.4,"LM0");
790    
791     TLegend * leg_S = new TLegend(0.5,0.14,0.85,0.38);
792     leg_S->SetBorderSize(0);leg_S->SetFillColor(0);
793     leg_S->SetHeader("95% CL upper limits");
794     leg_S->AddEntry(gCLsExpUncert,"Expected (CLs, TLimit)","l");
795     leg_S->AddEntry(gCLsObsUncert,"Observed (CLs, TLimit)","l");
796     //leg_S->AddEntry(gPLExpUncert,"Expected (PL, RooStat)","l");
797     //leg_S->AddEntry(gPLObsUncert,"Observed (PL, RooStat)","l");
798     //leg_S->AddEntry(gFCObsUncert,"Observed (FC, RooStat)","l");
799     //leg_S->AddEntry(gFCExpUncert,"Expected (FC, RooStat)","l");
800     //leg_S->AddEntry(gMCMCObsUncert,"Observed (MCMC, RooStat)","l");
801     //leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
802     leg_S->Draw();
803     c1->SaveAs("results/UncertaintyScan.pdf");
804 auterman 1.2
805 auterman 1.7 //theApp.Run();
806 auterman 1.1 }
807    
808    
809     int main(int argc, char** argv)
810     {
811     return plot(argc, argv);
812     }