ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/plot.cc
Revision: 1.21
Committed: Fri May 20 07:39:34 2011 UTC (13 years, 11 months ago) by auterman
Content type: text/plain
Branch: MAIN
Changes since 1.20: +314 -90 lines
Log Message: 
post CWR comments

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 auterman 1.21 genpoints->Fill("limits_moriond_MHT.old/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 auterman 1.21 genpointsHT->Fill("limits_moriond_HT.old/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.21 //genpoints->ExpandGrid(1);
116     //genpoints->ExpandGrid(1);
117     //genpoints->ExpandGrid(1);
118     //genpointsHT->ExpandGrid(1);
119     //genpointsHT->ExpandGrid(1);
120     //genpointsHT->ExpandGrid(1);
121    
122 auterman 1.6 //the plotting ----------------------------------------------------------------------
123 auterman 1.7 //plotting helper functions
124 auterman 1.4 PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
125 auterman 1.11 PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
126 auterman 1.18 //PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
127     //PlotTools<SusyScan> * plotToolsMHT = new PlotTools<SusyScan>(genpointsMHT->GetScan());
128 auterman 1.6 PlotTools<GeneratorMasses> * plotMasses = new PlotTools<GeneratorMasses>(genpoints->GetGeneratorMasses());
129 auterman 1.2
130 auterman 1.7 //iso mass lines
131 auterman 1.19 TGraph * gl500 = plotMasses->Line(Mzero, Mhalf, MGluino, 500, 0.8 );
132     TGraph * gl600 = plotMasses->Line(Mzero, Mhalf, MGluino, 600, 0.8 );
133     TGraph * gl700 = plotMasses->Line(Mzero, Mhalf, MGluino, 700, 0.8 );
134     TGraph * gl800 = plotMasses->Line(Mzero, Mhalf, MGluino, 800, 0.8 );
135     TGraph * gl900 = plotMasses->Line(Mzero, Mhalf, MGluino, 900, 0.8 );
136     TGraph * gl1000 = plotMasses->Line(Mzero, Mhalf, MGluino, 1000, 0.8);
137     TGraph * sq500 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 500, 1);
138     TGraph * sq600 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 600, 1);
139     TGraph * sq700 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 700, 1);
140     TGraph * sq800 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 800, 1);
141     TGraph * sq900 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 900, 1);
142     TGraph * sq1000 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 1000, 1);
143 auterman 1.7
144     TGraph * chi100 = plotMasses->Line(Mzero, Mhalf, MChi1, 50, 20);
145     TGraph * cha200 = plotMasses->Line(Mzero, Mhalf, MCha1, 200, 20);
146     cha200->SetLineColor(2);
147    
148     //the histograms
149     c1->SetLogz(1);
150     //h->SetMaximum(27);
151     //h->SetMinimum(0.01);
152 auterman 1.10
153     /**/
154 auterman 1.7 // cross-section in M0 - M1/2
155 auterman 1.21 c1->SetRightMargin ( 0.2 );
156 auterman 1.7 TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
157 auterman 1.18 100,0,1009.9,35,100,450);
158 auterman 1.7 plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
159     hxsec->SetMinimum(0.01);
160     //sq500->Draw();
161     //gl500->Draw();
162     hxsec->Draw("colz");
163 auterman 1.10 c1->SaveAs("results/Xsection_m0_m12_tb10.pdf");
164 auterman 1.7 std::string wait;
165     //std::cin>>wait;
166 auterman 1.10
167 auterman 1.7 // Observed Limit in M0 - M1/2
168     TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
169 auterman 1.18 100,0,1009.9,35,100,450);
170 auterman 1.7 plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
171     hobslimit->SetMinimum(0.01);
172     hobslimit->Draw("colz");
173 auterman 1.14 c1->SaveAs("results/ObsLimit_m0_m12_tb10.pdf");
174    
175 auterman 1.15
176 auterman 1.7 // Expected Limit in M0 - M1/2
177     TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
178 auterman 1.18 100,0,1009.9,35,100,450);
179 auterman 1.7 plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
180     hexplimit->SetMinimum(0.01);
181     hexplimit->Draw("colz");
182 auterman 1.10 c1->SaveAs("results/ExpLimit_m0_m12_tb10.pdf");
183 auterman 1.7
184 auterman 1.13 // Signal Acceptance in M0 - M1/2
185 auterman 1.21 c1->SetRightMargin ( 0.2 );
186     c1->SetTopMargin ( 0.11 );
187     c1->SetLogz(0);
188 auterman 1.7 TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
189 auterman 1.21 50,0,1009.9,35,100,450);
190     plotTools->Area(hsigacc, Mzero, Mhalf, NLOSignalAcceptance);
191     hsigacc->SetMinimum(0.0);
192     hsigacc->SetMaximum(0.35);
193     hsigacc->SetContour(14);
194     hsigacc->GetZaxis()->SetTitleOffset(1.5);
195 auterman 1.7 hsigacc->Draw("colz");
196 auterman 1.21 TLatex as; as.SetTextSize(0.025); as.SetTextFont(42);//ms.SetTextColor(12);
197     as.DrawLatex(520,461,"MHT selection, tan#beta=10, #mu>0, A_{0}=0");
198     as.SetTextSize(0.04);
199     as.DrawLatex( 50,461,"CMS");
200 auterman 1.14 c1->SaveAs("results/SigAccMHT_m0_m12_tb10.pdf");
201 auterman 1.18
202 auterman 1.14
203 auterman 1.19 // Signal Acceptance in M0 - M1/2
204     TH2F*hsigaccHT = new TH2F("sigaccHT",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
205 auterman 1.21 50,0,1009.9,35,100,450);
206     plotToolsHT->Area(hsigaccHT, Mzero, Mhalf, NLOSignalAcceptance);
207     hsigaccHT->SetMinimum(0.0);
208     hsigaccHT->SetMaximum(0.35);
209     hsigaccHT->SetContour(14);
210     hsigaccHT->GetZaxis()->SetTitleOffset(1.5);
211 auterman 1.19 hsigaccHT->Draw("colz");
212 auterman 1.21 as.SetTextSize(0.025); as.SetTextFont(42);//ms.SetTextColor(12);
213     as.DrawLatex(525,461,"HT selection, tan#beta=10, #mu>0, A_{0}=0");
214     as.SetTextSize(0.04);
215     as.DrawLatex( 50,461,"CMS");
216 auterman 1.19 c1->SaveAs("results/SigAccHT_m0_m12_tb10.pdf");
217 auterman 1.21 c1->SetTopMargin ( 0.1 );
218    
219 auterman 1.19
220 auterman 1.7 // Exp. Limit on Number of Signal Events in M0 - M1/2
221 auterman 1.8 c1->SetLogz(0);
222 auterman 1.7 TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
223 auterman 1.18 100,0,1009.9,35,100,450);
224 auterman 1.7 plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
225 auterman 1.10 hexplimitnsig->SetMinimum(0.0);
226     hexplimitnsig->SetMaximum(20);
227 auterman 1.7 hexplimitnsig->Draw("colz");
228 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_m0_m12_tb10.pdf");
229 auterman 1.7
230     // Obs. Limit on Number of Signal Events in M0 - M1/2
231     TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
232 auterman 1.18 100,0,1009.9,35,100,450);
233 auterman 1.7 plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
234 auterman 1.10 hobslimitnsig->SetMinimum(0.0);
235     hobslimitnsig->SetMaximum(20);
236 auterman 1.7 hobslimitnsig->Draw("colz");
237 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_m0_m12_tb10.pdf");
238 auterman 1.7
239     c1->SetLogz(0);
240     // Expected Exclusion in M0 - M1/2
241     TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
242 auterman 1.8 100,0,1009.9,35,100,450);
243 auterman 1.11 plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
244 auterman 1.7 hexpexcl->Draw("colz");
245 auterman 1.10 c1->SaveAs("results/ExpExclusion_m0_m12_tb10.pdf");
246 auterman 1.7
247     // Observed Exclusion in M0 - M1/2
248     TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
249 auterman 1.8 100,0,1009.9,35,100,450);
250 auterman 1.11 plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
251 auterman 1.7 hobsexcl->Draw("colz");
252 auterman 1.10 c1->SaveAs("results/ObsExclusion_m0_m12_tb10.pdf");
253 auterman 1.7
254 auterman 1.8 // Observed Exclusion in M0 - M1/2
255     TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
256     100,0,1009.9,35,100,450);
257     plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
258     hPLobsexcl->Draw("colz");
259 auterman 1.10 c1->SaveAs("results/PL_ObsExclusion_m0_m12_tb10.pdf");
260 auterman 1.8
261 auterman 1.20 // TestContours in M0 - M1/2 /////////////////////////////////////////////////////////////////////
262 auterman 1.19 c1->SetRightMargin ( 0.1 );
263 auterman 1.14 TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
264     100,0,1009.9,35,100,450);
265 auterman 1.20 TH2F*testExclold=(TH2F*)texcl->Clone();
266     TH2F*testExclnew=(TH2F*)texcl->Clone();
267 auterman 1.21 plotTools->Area(testExclold, Mzero, Mhalf, NLOHybridObsExclCL);
268     plotTools->Area(testExclnew, Mzero, Mhalf, NLOHybridObsCL );
269 auterman 1.20
270 auterman 1.21 //std::vector<TGraph*> contours = plotTools->GetContours005(testExclnew,3);
271     std::vector<TGraph*> contours = plotToolsHT->GetContours(testExclold,3);
272     testExclold->Draw("colz");
273     //testExclnew->Draw("colz");
274 auterman 1.14 int col=kBlue-10;
275 auterman 1.9 for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
276     if (! *cont) continue;
277     double x, y;
278     (*cont)->GetPoint(0, x, y);
279     (*cont)->SetLineColor(col);
280     (*cont)->Draw("l");
281     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col++);
282     char val[20];
283 auterman 1.13 sprintf(val,"%d",(int)(cont-contours.begin()));
284 auterman 1.9 l.DrawLatex(x,y,val);
285 auterman 1.11 //if (cont-contours.begin()>3) break;
286 auterman 1.9 }
287 auterman 1.10 c1->SaveAs("results/ExclusionTestContours_m0_m12_tb10.pdf");
288 auterman 1.9
289 auterman 1.20 /*
290 auterman 1.8 // Exclusion in M0 - M1/2
291     TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
292     100,0,1009.9,35,100,450);
293 auterman 1.19 TH2F*hs = new TH2F("hs","",100,0,1009.9,35,100,450);
294     TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
295     TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
296     TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
297 auterman 1.20 TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
298 auterman 1.19 TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
299     TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
300     TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
301     TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
302     TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
303     TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
304     TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
305     TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
306     TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
307     TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
308     TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
309     TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 1,1);
310     TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
311     TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
312     TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
313     TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
314     TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
315 auterman 1.20 TGraph * gCheck2 = (TGraph*)gCLsObsExclHT->Clone();
316     Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
317     Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
318     Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
319     Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
320     Smooth( gexpexcl, 27 );
321     Smooth( gobsexcl_LO, 27 );
322     Smooth( gobsexcl, 27 );
323     Smooth( gobsexcl, 27 );
324     Smooth( gexpexcl_LO, 27 );
325     Smooth( gFCexpexclHT, 27 );
326     Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
327     Smooth( gFCexpexclMHT, 27 );
328     Smooth( gFCexpexclMHT_LO, 27 );
329     Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
330     Smooth( gFCobsexclMHT, 27 );
331     Smooth( gCLsExpExclHTm1, 27 );
332     Smooth( gCLsExpExclHTp1, 27 );
333     Smooth( gCLsExpExclMHTm1, 27 );
334     Smooth( gCLsExpExclMHTp1, 27 );
335 auterman 1.19 TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
336 auterman 1.20 TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
337 auterman 1.19 TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
338     TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
339     TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
340 auterman 1.20 TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
341 auterman 1.19 TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
342     TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
343     TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
344     TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);
345     hexcl->GetYaxis()->SetTitleOffset(1.3);
346     hexcl->Draw("colz");
347 auterman 1.10 //set old exclusion Limits
348 auterman 1.16 TGraph* LEP_ch = set_lep_ch(10);
349     TGraph* LEP_sl = set_lep_sl(10);//slepton curve
350     TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
351     TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
352     //TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
353     //TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
354 auterman 1.20 TGraph* Atlas = Atlas_m0_m12_tb3_obs();
355     Atlas->SetLineColor( 28 );
356 auterman 1.16 TGraph* stau = set_tev_stau(10);//stau
357     TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
358     TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
359 auterman 1.15 TGraphErrors* First = getObserved_NLO_tanBeta10();
360     TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
361     TGraphErrors* Third = getLO_tanBeta10();
362 auterman 1.10 First->GetXaxis()->SetRangeUser(0,505);
363     First->GetYaxis()->SetRangeUser(80,500);
364     First->GetXaxis()->SetTitle("m_{0} (GeV)");
365     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
366     TSpline3 *sFirst = new TSpline3("sFirst",First);
367     sFirst->SetLineColor(kBlue);
368 auterman 1.18 sFirst->SetLineWidth(1);
369 auterman 1.10 TSpline3 *sSecond = new TSpline3("sSecond",Second);
370     sSecond->SetLineColor(kBlue);
371     sSecond->SetLineStyle(2);
372 auterman 1.18 sSecond->SetLineWidth(1);
373 auterman 1.10 TSpline3 *sThird = new TSpline3("sThird",Third);
374 auterman 1.18 sThird->SetLineColor(kBlue);
375 auterman 1.10 sThird->SetLineStyle(4);
376 auterman 1.18 sThird->SetLineWidth(1);
377 auterman 1.10 TEV_sn_d0_1->Draw("fsame");
378     TEV_sn_d0_2->Draw("fsame");
379 auterman 1.17 TEV_sg_d0->Draw("fsame");
380     TEV_sg_cdf->Draw("fsame");
381 auterman 1.10 LEP_ch->Draw("fsame");
382     LEP_sl->Draw("fsame");
383     stau->Draw("fsame");
384     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
385 auterman 1.19 b.DrawLatex( 10,380,"#tilde{#tau} LSP");
386 auterman 1.10 b.DrawLatex( 6,150,"LEP2");
387     b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
388     b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
389 auterman 1.17 b.DrawLatex(90,160,"CDF");
390     //b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
391 auterman 1.10 b.DrawLatex( 80,180,"D0 #tilde{#nu}");
392 auterman 1.18 if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
393     sFirst->Draw("same");
394 auterman 1.20 //sSecond->Draw("same");
395     //sThird->Draw("same");
396     Atlas->Draw("c,same");
397     TLatex ms; ms.SetTextSize(0.02); ms.SetTextColor(12);
398     ms.DrawLatex(150,433,"tan#beta=10, #mu>0, A_{0}=0");
399 auterman 1.19 //gCheck1->Draw("same");
400     //gCheck2->Draw("same");
401     //if (gexpexcl) gexpexcl->Draw("l");
402     //if (gobsexcl) gobsexcl->Draw("l");
403 auterman 1.18 //if (gexpexcl_LO) gexpexcl_LO->Draw("l");
404 auterman 1.19 //if (gMCMCexpexcl) gMCMCexpexcl->Draw("l");
405 auterman 1.18 //if (gCLsObsExclHT_LO) gCLsObsExclHT_LO->Draw("l");
406 auterman 1.19 //if (gCLsExpExclHT_LO) gCLsExpExclHT_LO->Draw("l");
407     //if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
408     //if (gCLsExpExclHT) gCLsExpExclHT->Draw("l");
409 auterman 1.18 //if (gobsexcl) gobsexcl->Draw("l");
410     //if (gPLobsexcl) gPLobsexcl->Draw("l");
411     //if (gFCobsexcl) gFCobsexcl->Draw("l");
412     //if (gMCMCobsexcl) gMCMCobsexcl->Draw("l");
413     //if (gFCexpexcl) gFCexpexcl->Draw("l");
414     //if (gPLexpexcl) gPLexpexcl->Draw("l");
415     //if (gFCexpexcl) gFCexpexcl->Draw("l");
416 auterman 1.19
417     CLsObsNLO->Draw("l");
418 auterman 1.20 CLsObsLO->Draw("l");
419 auterman 1.19 CLsExpNLO->Draw("l");
420 auterman 1.20 //CLsExpLO->Draw("l");
421     FCExpNLO->Draw("l");
422 auterman 1.19 //FCObsNLO->Draw("l");
423    
424     //gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c");
425     //sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c");
426 auterman 1.20 TLegend * leg = new TLegend(0.54,0.7,0.89,0.93);
427 auterman 1.19 leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);
428     TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
429     leg->SetHeader("CMS preliminary, 95%CL");
430 auterman 1.20 leg->AddEntry(expLeg, "Expected (CLs, NLO)","lf");
431 auterman 1.19 leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l");
432 auterman 1.20 leg->AddEntry(CLsObsLO, "Observed (CLs, LO)","l");
433     leg->AddEntry(FCExpLO, "Expected (FC, NLO)","l");
434     //leg->AddEntry(FCExpLO, "Expected (FC, LO)","l");
435 auterman 1.12 //if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
436 auterman 1.19 //if (gCLsExpExclHT) leg->AddEntry(expLeg ,"NLO Expected (HT, CLs)","l");
437     //if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
438 auterman 1.18 //if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
439 auterman 1.19 //if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
440 auterman 1.11 //if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
441     //if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
442 auterman 1.19 //if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l");
443     //if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l");
444 auterman 1.12 //if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
445 auterman 1.19 //if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
446     //if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
447     //if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l");
448 auterman 1.20 //if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)");
449     if (sFirst) leg->AddEntry(sFirst, "CMS #alpha_{T} Observed");
450     //if (sThird) leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)");
451     if (Atlas) leg->AddEntry(Atlas, "Atlas Observed (tan#beta=3)","l");
452 auterman 1.8 leg->Draw();
453 auterman 1.10 gPad->RedrawAxis();
454 auterman 1.13 c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
455 auterman 1.20 */
456     // Exclusion in M0 - M1/2
457     TH2F*hexcl = new TH2F("hexcl",";m_{0} (GeV); m_{1/2} (GeV); 95% CL Exclusion",
458     100,0,1009.9,35,110,450);
459     TH2F*hs = new TH2F("hs","",100,0,1009.9,35,100,450);
460     TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
461     TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
462     TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
463     TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
464     TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
465     TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
466     TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
467     TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
468     TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
469     TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
470     TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
471     TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
472     TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
473 auterman 1.21 //TGraph * gCLsExpNoSExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
474     //TGraph * gCLsExpNoSExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
475 auterman 1.20 TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
476     TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
477     TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
478     TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
479     TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
480     TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
481     TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
482 auterman 1.21 TGraph * gHybridCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
483     TGraph * gHybridCLsObsExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
484     TGraph * gCheck1 = (TGraph*)gCLsObsExclHT->Clone();
485     TGraph * gCheck2 = (TGraph*)gCLsExpExclHT->Clone();
486 auterman 1.20 Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
487 auterman 1.21 //Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
488     //Smooth( gCLsExpNoSExclHT, 27 ); gCLsExpNoSExclHT->SetLineWidth( 3 );
489 auterman 1.20 Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
490     Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
491     Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
492     Smooth( gexpexcl, 27 );
493     Smooth( gobsexcl_LO, 27 );
494     Smooth( gobsexcl, 27 );
495     Smooth( gexpexcl_LO, 27 );
496     Smooth( gFCexpexclHT, 27 );
497     Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
498     Smooth( gFCexpexclMHT, 27 );
499     Smooth( gFCexpexclMHT_LO, 27 );
500     Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
501     Smooth( gFCobsexclMHT, 27 );
502     Smooth( gCLsExpExclHTm1, 27 );
503     Smooth( gCLsExpExclHTp1, 27 );
504     Smooth( gCLsExpExclMHTm1, 27 );
505     Smooth( gCLsExpExclMHTp1, 27 );
506 auterman 1.21 // Smooth( gHybridCLsObsExclHT, 27); gHybridCLsObsExclHT->SetLineWidth( 3 );
507     // Smooth( gHybridCLsObsExclMHT, 27);gHybridCLsObsExclMHT->SetLineWidth( 3 );
508 auterman 1.20 TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
509     TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
510     TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
511 auterman 1.21 //TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
512 auterman 1.20 TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
513     TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
514     TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
515     TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
516 auterman 1.21 TGraph * HybridCLsObsNLO = plotToolsHT->ChooseBest(gHybridCLsObsExclHT,gHybridCLsObsExclMHT, gHybridCLsObsExclHT,gHybridCLsObsExclMHT);
517 auterman 1.20 TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
518     TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
519 auterman 1.21 TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);gCLsExp1Sigma->SetFillStyle(4010);
520 auterman 1.20 hexcl->GetYaxis()->SetTitleOffset(1.3);
521     hexcl->GetXaxis()->SetTitleOffset(0.92);
522     hexcl->Draw("colz");
523     //set old exclusion Limits
524     TGraph* LEP_ch = set_lep_ch(10);
525     TGraph* LEP_sl = set_lep_sl(10);//slepton curve
526     TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
527     TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
528     //TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
529     //TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
530     TGraph* Atlas = Atlas_m0_m12_tb3_obs();
531     Atlas->SetLineColor( 28 );
532     TGraph* stau = set_tev_stau(10);//stau
533     TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
534     TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
535     TGraphErrors* First = getObserved_NLO_tanBeta10();
536     TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
537     TGraphErrors* Third = getLO_tanBeta10();
538     First->GetXaxis()->SetRangeUser(0,505);
539     First->GetYaxis()->SetRangeUser(80,500);
540     First->GetXaxis()->SetTitle("m_{0} (GeV)");
541     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
542     TSpline3 *sFirst = new TSpline3("sFirst",First); sFirst->SetLineColor(kBlue);sFirst->SetLineWidth(1);
543     TSpline3 *sSecond = new TSpline3("sSecond",Second);
544     sSecond->SetLineColor(kBlue); sSecond->SetLineStyle(2); sSecond->SetLineWidth(1);
545     TSpline3 *sThird = new TSpline3("sThird",Third);
546     sThird->SetLineColor(kBlue); sThird->SetLineStyle(4); sThird->SetLineWidth(1);
547     //TEV_sn_d0_1->Draw("fsame");
548     //TEV_sn_d0_2->Draw("fsame"); //only for tb=3
549     LEP_ch->Draw("fsame");
550     LEP_sl->Draw("fsame");
551     TEV_sg_cdf->Draw("fsame");
552     TEV_sg_cdf->Draw("lsame");
553     TEV_sg_d0->Draw("fsame");
554     TEV_sg_d0->Draw("lsame");
555     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
556     //b.DrawLatex( 6,150,"LEP2");
557     //b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
558     //b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
559     //b.DrawLatex(90,160,"CDF");
560     //b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
561     //b.DrawLatex( 80,180,"D0 #tilde{#nu}");
562 auterman 1.21 // if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
563 auterman 1.20 sFirst->Draw("same");
564     //sSecond->Draw("same");
565     //sThird->Draw("same");
566     Atlas->Draw("c,same");
567     TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
568     ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
569 auterman 1.21
570     //gCheck1->Draw("same");
571     //gCheck2->Draw("same");
572    
573 auterman 1.20 //LM points
574     TMarker* LM0 = new TMarker(200.,160.,20);
575     TMarker* LM1 = new TMarker(60.,250.,20);
576     LM0->SetMarkerSize(0.7);
577     LM1->SetMarkerSize(0.7);
578     TLatex* tLM0 = new TLatex(205.,160.,"LM0");
579     tLM0->SetTextSize(0.02);
580     TLatex* tLM1 = new TLatex(65.,243.,"LM1");
581     tLM1->SetTextSize(0.02);
582     LM0->Draw("same");
583     tLM0->Draw("same");
584     LM1->Draw("same");
585     tLM1->Draw("same");
586 auterman 1.21 TLegend* legexp = new TLegend(0.61,0.61,0.93,0.88,NULL,"brNDC");
587 auterman 1.20 legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
588     //TEV_sg_cdf.SetLineColor(1);
589 auterman 1.21 legexp->SetHeader("CMS");
590 auterman 1.20 legexp->AddEntry(TEV_sg_cdf,"CDF #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
591     legexp->AddEntry(TEV_sg_d0,"D0 #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=3, #mu<0}","f");
592     //ch_gr.SetLineColor(1);
593     legexp->AddEntry(LEP_ch,"LEP2 #tilde{#chi}_{1}^{#pm}","f");
594     //sl_gr.SetLineColor(1);
595     legexp->AddEntry(LEP_sl,"LEP2 #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
596     //if(tanbeta == 3)
597     //legexp->AddEntry(TEV_sn_d0_1,"D0 #chi^{#pm}_{1}, #chi^{0}_{2}","f");
598     legexp->AddEntry(sFirst, "CMS #alpha_{T}");
599 auterman 1.21 //legexp->AddEntry(Atlas, "Atlas, #scale[0.8]{tan#beta=3}","l");
600 auterman 1.20 legexp->Draw();
601    
602     CLsObsNLO->Draw("l");
603     CLsObsLO->Draw("l");
604     CLsExpNLO->Draw("l");
605     //FCExpNLO->Draw("l");
606 auterman 1.21 //CLsExpNoSNLO->Draw("l");
607     // HybridCLsObsNLO->Draw("l");
608     //gobsexcl->Draw("l");
609     //gHybridCLsObsExclMHT->Draw("l");
610    
611    
612 auterman 1.20 stau->Draw("fsame");
613     b.DrawLatex( 10,380,"#tilde{#tau} LSP");
614     //constant ssqquark and gluino lines
615     TF1* lnsq[4];
616     TF1* lngl[4];
617     TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
618     sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
619     sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
620     TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
621     glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
622     glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
623     glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
624     int tanBeta_=10;
625     for(int i = 0; i < 4; i++){
626     lnsq[i] = constant_squark(tanBeta_,i);
627     lngl[i] = constant_gluino(tanBeta_,i);
628     lngl[i]->Draw("same");
629     lnsq[i]->Draw("same");
630     }
631 auterman 1.21 TLegend * leg = new TLegend(0.25,0.76,0.6,0.88);
632 auterman 1.20 leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
633     TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
634     leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
635     leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
636     leg->AddEntry(CLsObsLO, "Observed, LO","l");
637     leg->AddEntry(expLeg, "Expected #pm 1#sigma, NLO","lf");
638 auterman 1.21 //leg->AddEntry(CLsExpNoSNLO, "Expected, no-signal hyp., NLO","l");
639 auterman 1.20 leg->Draw();
640     gPad->RedrawAxis();
641     c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
642 auterman 1.21
643    
644     //plotToolsHT->Print(NLOSignal, Mzero,Mhalf, gCLsObsExclHT, 5);
645     plotToolsHT->Print(NLOXsection, Mzero,Mhalf, MGluino, MSquarkL, gCLsObsExclHT, 5);
646 auterman 1.20
647     hexcl->Draw("");
648     gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
649     CLsObsNLO->SetName("CLsNLO_ObservedTb10");
650     CLsExpNLO->SetName("CLsNLO_ExpectedTb10");
651     CLsExpLO->SetName("CLsLO_ExpectedTb10");
652     gCLsExp1Sigma->Draw("l");
653     CLsObsNLO->Draw("l");
654     CLsExpNLO->Draw("l");
655     CLsObsLO->Draw("l");
656     c1->SaveAs("results/LimitContours_tb10.C");
657     // ===============================================================================================================
658 auterman 1.8
659 auterman 1.21
660     //The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
661     TGraph * gCLsExpNoSExclMHT = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 1,2);
662     TGraph * gCLsExpNoSExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 1,2);
663     TGraph * gCLsExpNoSExclHTm1 = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
664     TGraph * gCLsExpNoSExclHTp1 = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
665     TGraph * gCLsExpNoSExclMHTm1 = plotTools ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
666     TGraph * gCLsExpNoSExclMHTp1 = plotTools ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
667     Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
668     Smooth( gCLsExpNoSExclHT, 27 ); gCLsExpNoSExclHT->SetLineWidth( 3 );
669     Smooth( gCLsExpNoSExclHTm1, 27 );
670     Smooth( gCLsExpNoSExclHTp1, 27 );
671     Smooth( gCLsExpNoSExclMHTm1, 27 );
672     Smooth( gCLsExpNoSExclMHTp1, 27 );
673     TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
674     TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
675     TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
676     TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
677     hexcl->Draw("colz");
678     LEP_ch->Draw("fsame");
679     LEP_sl->Draw("fsame");
680     TEV_sg_cdf->Draw("fsame");
681     TEV_sg_cdf->Draw("lsame");
682     TEV_sg_d0->Draw("fsame");
683     TEV_sg_d0->Draw("lsame");
684     ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
685     LM0->SetMarkerSize(0.7);
686     LM1->SetMarkerSize(0.7);
687     tLM0->SetTextSize(0.02);
688     tLM1->SetTextSize(0.02);
689     LM0->Draw("same");
690     tLM0->Draw("same");
691     LM1->Draw("same");
692     tLM1->Draw("same");
693     gCLsExpNoS1Sigma->Draw("lf,same");
694     sFirst->Draw("same");
695     //Atlas->Draw("c,same");
696     CLsObsNLO->Draw("l,same");
697     CLsObsLO->Draw("l,same");
698     CLsExpNoSNLO->Draw("l,same");
699     //FCExpNLO->Draw("l");
700     stau->Draw("fsame");
701     b.DrawLatex( 10,380,"#tilde{#tau} LSP");
702     //constant ssqquark and gluino lines
703     sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
704     sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
705     glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
706     glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
707     glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
708     for(int i = 0; i < 4; i++){
709     lngl[i]->Draw("same");
710     lnsq[i]->Draw("same");
711     }
712     legexp->Draw();
713     leg->Draw();
714     gPad->RedrawAxis();
715     c1->SaveAs("results/Exclusion_m0_m12_tb10_NoSigHypPseudoData.pdf");
716    
717    
718    
719     hexcl->Draw("colz");
720     LEP_ch->Draw("fsame");
721     LEP_sl->Draw("fsame");
722     TEV_sg_cdf->Draw("fsame");
723     TEV_sg_cdf->Draw("lsame");
724     TEV_sg_d0->Draw("fsame");
725     TEV_sg_d0->Draw("lsame");
726     gCLsExpNoS1Sigma->Draw("lf,same");
727     sFirst->Draw("same");
728     Atlas->Draw("c,same");
729     CLsObsNLO->Draw("l,same");
730     CLsExpNoSNLO->Draw("l,same");
731     ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
732     LM0->SetMarkerSize(0.7);
733     LM1->SetMarkerSize(0.7);
734     tLM0->SetTextSize(0.02);
735     tLM1->SetTextSize(0.02);
736     LM0->Draw("same");
737     tLM0->Draw("same");
738     LM1->Draw("same");
739     tLM1->Draw("same");
740     stau->Draw("fsame");
741     b.DrawLatex( 10,380,"#tilde{#tau} LSP");
742     //Jim's limits///////////////////////////////////////////////////////////////////////////////
743     TGraph* JimObsHT = Jim_ht_tb10(0);
744     TGraph* JimObsMHT = Jim_mht_tb10(0);
745     TGraph* JimExpHT = Jim_ht_tb10(1);
746     TGraph* JimExpMHT = Jim_mht_tb10(1);
747     TGraph* JimExpHTup = Jim_ht_tb10(2);
748     TGraph* JimExpMHTup = Jim_mht_tb10(2);
749     TGraph* JimExpHTdn = Jim_ht_tb10(3);
750     TGraph* JimExpMHTdn = Jim_mht_tb10(3);
751     Smooth(JimObsHT, 2);
752     Smooth(JimObsMHT,2);
753     Smooth(JimExpHT, 2);
754     Smooth(JimExpMHT,2);
755     Smooth(JimExpHTup, 2);
756     Smooth(JimExpMHTup,2);
757     Smooth(JimExpHTdn, 2);
758     Smooth(JimExpMHTdn,2);
759     TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
760     TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
761     TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
762     TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
763     TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
764     TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
765     JimExp1Sigma->Draw("f");
766     JimObs->Draw("c");
767     JimExp->Draw("c");
768     stau->Draw("fsame");
769     b.DrawLatex( 10,380,"#tilde{#tau} LSP");
770     //constant ssqquark and gluino lines
771     sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
772     sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
773     glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
774     glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
775     glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
776     for(int i = 0; i < 4; i++){
777     lngl[i]->Draw("same");
778     lnsq[i]->Draw("same");
779     }
780     legexp->Draw();
781     TLegend * legBayes = new TLegend(0.25,0.730,0.6,0.88);
782     legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
783     legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
784     legBayes->AddEntry(JimObs, "Observed, Bayes","l");
785     legBayes->AddEntry(JimLeg, "Expected #pm 1#sigma, Bayes","lf");
786     legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
787     legBayes->AddEntry(expLeg, "Expected #pm 1#sigma, CLs","lf");
788     legBayes->Draw();
789     gPad->RedrawAxis();
790     c1->SaveAs("results/Exclusion_m0_m12_tb10_Bayes.pdf");
791    
792    
793    
794 auterman 1.13 // Signal Contamination in M0 - M1/2
795     c1->SetLogz(0);
796 auterman 1.21 c1->SetRightMargin ( 0.2 );
797 auterman 1.13 TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
798 auterman 1.18 100,0,1009.9,35,100,450);
799     plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
800 auterman 1.13 hsigcont->SetMinimum(0.01);
801     hsigcont->SetMaximum(20);
802     hsigcont->Draw("colz");
803     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
804     c1->SaveAs("results/SignalContamination_HT_m0_m12_tb10.pdf");
805    
806     // Signal Contamination in M0 - M1/2
807     c1->SetLogz(0);
808     TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
809 auterman 1.18 100,0,1009.9,35,100,450);
810     plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
811 auterman 1.13 hsigcontMHT->SetMinimum(0.01);
812     hsigcontMHT->SetMaximum(20);
813     hsigcontMHT->Draw("colz");
814     if (gobsexcl) gobsexcl->Draw("l");
815     c1->SaveAs("results/SignalContamination_MHT_m0_m12_tb10.pdf");
816    
817 auterman 1.18 // Signal JEC Uncertainty MHT
818     c1->SetLogz(0);
819     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
820     100,0,1009.9,35,100,450);
821     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
822     hsig_jec_mht->SetMinimum(0.0);
823     hsig_jec_mht->SetMaximum(0.25);
824     hsig_jec_mht->SetContour(10);
825     hsig_jec_mht->Draw("colz");
826     if (gobsexcl) gobsexcl->Draw("l");
827     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
828    
829     // Signal JEC Uncertainty HT
830     c1->SetLogz(0);
831     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
832     100,0,1009.9,35,100,450);
833 auterman 1.19 plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
834 auterman 1.18 hsig_jec_ht->SetMinimum(0.0);
835     hsig_jec_ht->SetMaximum(0.25);
836     hsig_jec_ht->SetContour(10);
837     hsig_jec_ht->Draw("colz");
838     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
839     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
840    
841    
842     // Signal MuIso Uncertainty MHT
843     c1->SetLogz(0);
844     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
845     100,0,1009.9,35,100,450);
846     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
847 auterman 1.19 hsig_MuIso_mht->SetMinimum(0.0);hsig_MuIso_mht->GetZaxis()->SetTitleOffset(1.3);
848 auterman 1.18 hsig_MuIso_mht->SetMaximum(0.05);
849     hsig_MuIso_mht->SetContour(5);
850     hsig_MuIso_mht->Draw("colz");
851     if (gobsexcl) gobsexcl->Draw("l");
852     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
853    
854     // Signal MuIso Uncertainty HT
855     c1->SetLogz(0);
856     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
857     100,0,1009.9,35,100,450);
858 auterman 1.19 plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
859     hsig_MuIso_ht->SetMinimum(0.0);hsig_MuIso_ht->GetZaxis()->SetTitleOffset(1.3);
860 auterman 1.18 hsig_MuIso_ht->SetMaximum(0.05);
861     hsig_MuIso_ht->SetContour(5);
862     hsig_MuIso_ht->Draw("colz");
863     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
864     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
865    
866    
867     // Signal kFactor Uncertainty MHT
868     c1->SetLogz(0);
869 auterman 1.19 TH2F*hsig_kFactorUncert_mht = new TH2F("sigkFactorUncertmht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
870 auterman 1.18 100,0,1009.9,35,100,450);
871 auterman 1.19 plotTools->Area(hsig_kFactorUncert_mht, Mzero, Mhalf, SignalUncertKfactor);
872     hsig_kFactorUncert_mht->SetMinimum(0.1);hsig_kFactorUncert_mht->GetZaxis()->SetTitleOffset(1.3);
873     hsig_kFactorUncert_mht->SetMaximum(0.2);
874     hsig_kFactorUncert_mht->SetContour(10);
875     hsig_kFactorUncert_mht->Draw("colz");
876 auterman 1.18 if (gobsexcl) gobsexcl->Draw("l");
877 auterman 1.19 c1->SaveAs("results/SigkFactorUncert_MHT_m0_m12_tb10.pdf");
878 auterman 1.18
879     // Signal kFactor Uncertainty HT
880     c1->SetLogz(0);
881 auterman 1.19 TH2F*hsig_kFactorUncert_ht = new TH2F("sigkFactorUncertht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
882 auterman 1.18 100,0,1009.9,35,100,450);
883 auterman 1.19 plotToolsHT->Area(hsig_kFactorUncert_ht, Mzero, Mhalf, SignalUncertKfactor);
884     hsig_kFactorUncert_ht->SetMinimum(0.1);hsig_kFactorUncert_ht->GetZaxis()->SetTitleOffset(1.3);
885     hsig_kFactorUncert_ht->SetMaximum(0.2);
886     hsig_kFactorUncert_ht->SetContour(10);
887     hsig_kFactorUncert_ht->Draw("colz");
888     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
889     c1->SaveAs("results/SigkFactorUncert_HT_m0_m12_tb10.pdf");
890    
891     // Signal kFactor HT
892     c1->SetLogz(0);
893     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
894     100,0,1009.9,35,100,450);
895     plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalKfactor);
896     hsig_kFactor_ht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
897     hsig_kFactor_ht->SetMaximum(2.4);
898     hsig_kFactor_ht->SetContour(7);
899 auterman 1.18 hsig_kFactor_ht->Draw("colz");
900     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
901     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
902    
903 auterman 1.19 // Signal kFactor MHT
904     c1->SetLogz(0);
905     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
906     100,0,1009.9,35,100,450);
907     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalKfactor);
908     hsig_kFactor_mht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
909     hsig_kFactor_mht->SetMaximum(2.4);
910     hsig_kFactor_mht->SetContour(7);
911     hsig_kFactor_mht->Draw("colz");
912     if (gobsexcl) gobsexcl->Draw("l");
913     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
914    
915 auterman 1.1
916 auterman 1.2
917 auterman 1.7 //-----------------------------------------------------------------------------------
918     c1->SetLogz(1);
919    
920     TGraph * mz500 = plotMasses->Line( MSquarkL, MGluino, Mzero, 500, 1);
921     TGraph * mh250 = plotMasses->Line( MSquarkL, MGluino, Mhalf, 250, 1);
922    
923 auterman 1.8 // cross-section in squark - gluino mass
924 auterman 1.7 TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
925     60,200,1400,50,200,1200);
926     plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
927     hxsec_qg->SetMinimum(0.01);
928     hxsec_qg->Draw("colz");
929     mz500->Draw();
930     mh250->Draw();
931 auterman 1.10 c1->SaveAs("results/Xsection_mSql_mGl_tb10.pdf");
932 auterman 1.7
933 auterman 1.8 // Observed Limit in squark - gluino mass
934 auterman 1.7 TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
935     60,200,1400,50,200,1200);
936     plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
937     hobslimit_qg->SetMinimum(0.01);
938     hobslimit_qg->Draw("colz");
939 auterman 1.10 c1->SaveAs("results/ObsLimit_mSql_mGl_tb10.pdf");
940 auterman 1.7
941 auterman 1.8 // Expected Limit in squark - gluino mass
942 auterman 1.7 TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
943     60,200,1400,50,200,1200);
944     plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
945     hexplimit_qg->SetMinimum(0.01);
946     hexplimit_qg->Draw("colz");
947 auterman 1.10 c1->SaveAs("results/ExpLimit_mSql_mGl_tb10.pdf");
948 auterman 1.7
949 auterman 1.8 // Signal Acceptance in squark - gluino mass
950 auterman 1.7 TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
951     60,200,1400,50,200,1200);
952     plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
953     hsigacc_qg->SetMinimum(0.01);
954     hsigacc_qg->SetMaximum(1.0);
955     hsigacc_qg->Draw("colz");
956 auterman 1.14 c1->SaveAs("results/SigAccMHT_mSql_mGl_tb10.pdf");
957    
958     // Signal Acceptance in squark - gluino mass
959     TH2F*hsigaccHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
960     60,200,1400,50,200,1200);
961     plotToolsHT->Area(hsigaccHT_qg, MSquarkL, MGluino, SignalAcceptance);
962     hsigaccHT_qg->SetMinimum(0.01);
963     hsigaccHT_qg->SetMaximum(1.0);
964     hsigaccHT_qg->Draw("colz");
965     c1->SaveAs("results/SigAccHT_mSql_mGl_tb10.pdf");
966 auterman 1.19
967    
968     // Signal Contamination in squark - gluino mass
969     c1->SetLogz(0);
970     TH2F*hsigcont_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Contamination in bkg yield",
971     60,200,1400,50,200,1200);
972     plotTools->Area(hsigcont_qg, MSquarkL, MGluino, SignalContamination);
973     hsigcont_qg->SetMinimum(0.0);
974     hsigcont_qg->SetMaximum(20.0);
975     hsigcont_qg->Draw("colz");
976     c1->SaveAs("results/SignalContamination_MHT_mSql_mGl_tb10.pdf");
977    
978     // Signal Contamination in squark - gluino mass
979     TH2F*hsigcontHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance in bkg yield",
980     60,200,1400,50,200,1200);
981     plotToolsHT->Area(hsigcontHT_qg, MSquarkL, MGluino, SignalContamination);
982     hsigcontHT_qg->SetMinimum(0.0);
983     hsigcontHT_qg->SetMaximum(20.0);
984     hsigcontHT_qg->Draw("colz");
985     c1->SaveAs("results/SignalContamination_HT_mSql_mGl_tb10.pdf");
986 auterman 1.18 /*
987     // Signal JEC Uncertainty MHT
988     c1->SetLogz(0);
989     TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
990     100,0,1009.9,35,100,450);
991     plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
992     hsig_jec_mht->SetMinimum(0.0);
993     hsig_jec_mht->SetMaximum(0.25);
994     hsig_jec_mht->SetContour(10);
995     hsig_jec_mht->Draw("colz");
996     if (gobsexcl) gobsexcl->Draw("l");
997     c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
998    
999     // Signal JEC Uncertainty HT
1000     c1->SetLogz(0);
1001     TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
1002     100,0,1009.9,35,100,450);
1003     plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
1004     hsig_jec_ht->SetMinimum(0.0);
1005     hsig_jec_ht->SetMaximum(0.25);
1006     hsig_jec_ht->SetContour(10);
1007     hsig_jec_ht->Draw("colz");
1008     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
1009     c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
1010    
1011 auterman 1.15
1012 auterman 1.18 // Signal MuIso Uncertainty MHT
1013     c1->SetLogz(0);
1014     TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
1015     100,0,1009.9,35,100,450);
1016     plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
1017     hsig_MuIso_mht->SetMinimum(0.0);
1018     hsig_MuIso_mht->SetMaximum(0.05);
1019     hsig_MuIso_mht->SetContour(5);
1020     hsig_MuIso_mht->Draw("colz");
1021     if (gobsexcl) gobsexcl->Draw("l");
1022     c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
1023    
1024     // Signal MuIso Uncertainty HT
1025     c1->SetLogz(0);
1026     TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
1027     100,0,1009.9,35,100,450);
1028     plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
1029     hsig_MuIso_ht->SetMinimum(0.0);
1030     hsig_MuIso_ht->SetMaximum(0.05);
1031     hsig_MuIso_ht->SetContour(5);
1032     hsig_MuIso_ht->Draw("colz");
1033     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
1034     c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
1035    
1036    
1037     // Signal kFactor Uncertainty MHT
1038     c1->SetLogz(0);
1039     TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
1040     100,0,1009.9,35,100,450);
1041     plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
1042     hsig_kFactor_mht->SetMinimum(0.1);
1043     hsig_kFactor_mht->SetMaximum(0.2);
1044     hsig_kFactor_mht->SetContour(10);
1045     hsig_kFactor_mht->Draw("colz");
1046     if (gobsexcl) gobsexcl->Draw("l");
1047     c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
1048    
1049     // Signal kFactor Uncertainty HT
1050     c1->SetLogz(0);
1051     TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
1052     100,0,1009.9,35,100,450);
1053     plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
1054     hsig_kFactor_ht->SetMinimum(0.1);
1055     hsig_kFactor_ht->SetMaximum(0.2);
1056     hsig_kFactor_ht->SetContour(10);
1057     hsig_kFactor_ht->Draw("colz");
1058     if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
1059     c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
1060 auterman 1.7
1061 auterman 1.18 */
1062 auterman 1.8 // Exp. Limit on Number of Signal Events in squark - gluino mass
1063     c1->SetLogz(0);
1064 auterman 1.7 TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
1065     60,200,1400,50,200,1200);
1066     plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
1067 auterman 1.10 hexplimitnsig_qg->SetMinimum(0.0);
1068     hexplimitnsig_qg->SetMaximum(20);
1069 auterman 1.7 hexplimitnsig_qg->Draw("colz");
1070 auterman 1.10 c1->SaveAs("results/ExpLimitOnNSig_mSql_mGl_tb10.pdf");
1071 auterman 1.6
1072 auterman 1.8 // Obs. Limit on Number of Signal Events in squark - gluino mass
1073 auterman 1.7 TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
1074     60,200,1400,50,200,1200);
1075     plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
1076 auterman 1.10 hobslimitnsig_qg->SetMinimum(0.0);
1077     hobslimitnsig_qg->SetMaximum(20);
1078 auterman 1.7 hobslimitnsig_qg->Draw("colz");
1079 auterman 1.10 c1->SaveAs("results/ObsLimitOnNSig_mSql_mGl_tb10.pdf");
1080 auterman 1.7
1081     c1->SetLogz(0);
1082 auterman 1.8 // Expected Exclusion in squark - gluino mass
1083 auterman 1.7 TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
1084 auterman 1.2 60,200,1400,50,200,1200);
1085 auterman 1.11 plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
1086 auterman 1.7 hexpexcl_qg->Draw("colz");
1087 auterman 1.10 c1->SaveAs("results/ExpExclusion_mSql_mGl_tb10.pdf");
1088 auterman 1.2
1089 auterman 1.8 // Observed Exclusion in squark - gluino mass
1090 auterman 1.7 TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
1091     60,200,1400,50,200,1200);
1092 auterman 1.11 plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
1093 auterman 1.8 //std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
1094 auterman 1.7 hobsexcl_qg->Draw("colz");
1095 auterman 1.8 //for (std::vector<TGraph *>::iterator g=vobsexcl_qg.begin();g!=vobsexcl_qg.end();++g){
1096     // if (! *g) continue;
1097     // // (*g)->Draw("l");
1098     //}
1099 auterman 1.10 c1->SaveAs("results/ObsExclusion_mSql_mGl_tb10.pdf");
1100 auterman 1.2
1101 auterman 1.9
1102 auterman 1.8 // Observed Exclusion in squark - gluino mass
1103     TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
1104     60,200,1400,50,200,1200);
1105     plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
1106 auterman 1.10 //hPLobsexcl_qg->Draw("colz");
1107     //c1->SaveAs("results/PL_ObsExclusion_mSql_mGl_tb10.pdf");
1108 auterman 1.8
1109 auterman 1.11
1110    
1111 auterman 1.10 // Exclusion in squark - gluino mass ----------------------------------------
1112 auterman 1.21 c1->SetRightMargin ( 0.05 );
1113 auterman 1.19 TH2F*h_qg = new TH2F("hExcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
1114     60,0,1200,50,0,1200);
1115     TH2F*hExcl_qg = new TH2F("h_qg","",60,200,1400,50,200,1200);
1116     TH2F*hexpFC_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
1117     plotToolsHT->Area(hexpFC_HT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
1118     TH2F*hexpFC_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1119     plotToolsHT->Area(hexpFC_HT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
1120     TH2F*hexpFC_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
1121     plotTools->Area(hexpFC_MHT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
1122     TH2F*hexpFC_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1123     plotTools->Area(hexpFC_MHT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
1124     //
1125     TH2F*hexpCLs_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
1126     plotToolsHT->Area(hexpCLs_HT_LO_qg, MSquarkL, MGluino, ExpExclCL);
1127     TH2F*hexpCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1128     plotToolsHT->Area(hexpCLs_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
1129     TH2F*hexpCLs_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
1130     plotTools->Area(hexpCLs_MHT_LO_qg, MSquarkL, MGluino, ExpExclCL);
1131     TH2F*hexpCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1132     plotTools->Area(hexpCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
1133     //
1134     TH2F*hexp1p_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1135     plotToolsHT->Area(hexp1p_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
1136     TH2F*hexp1p_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1137     plotTools->Area(hexp1p_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
1138     TH2F*hexp1m_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1139     plotToolsHT->Area(hexp1m_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
1140     TH2F*hexp1m_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1141     plotTools->Area(hexp1m_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
1142    
1143     TH2F*hobsCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1144     plotToolsHT->Area(hobsCLs_HT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
1145     TH2F*hobsCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
1146     plotTools->Area(hobsCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
1147     //
1148     TH2 * hexpFC_LO_qg = plotToolsHT->BinWiseOr(hexpFC_HT_LO_qg, hexpFC_MHT_LO_qg);
1149     TH2 * hexpCLs_LO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_qg, hexpCLs_MHT_LO_qg);
1150     TH2 * hexpCLs_NLO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_qg, hexpCLs_MHT_NLO_qg);
1151     TH2 * hobsCLs_NLO_qg = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_qg, hobsCLs_MHT_NLO_qg);
1152     TH2 * hexp1p_NLO_qg = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_qg, hexp1p_MHT_NLO_qg);
1153     TH2 * hexp1m_NLO_qg = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_qg, hexp1m_MHT_NLO_qg);
1154     //
1155     TGraph * gexpFC_LO_qg = plotToolsHT->GetContour(hexpFC_LO_qg,3,0);
1156     gexpFC_LO_qg->SetLineStyle(4); gexpFC_LO_qg->SetLineColor(3);gexpFC_LO_qg->SetLineWidth(2);
1157     TGraph * gexpCLs_LO_qg = plotToolsHT->GetContour(hexpCLs_LO_qg,3,0);
1158     gexpCLs_LO_qg->SetLineStyle(4); gexpCLs_LO_qg->SetLineColor(1);gexpCLs_LO_qg->SetLineWidth(2);
1159     TGraph * gexpCLs_NLO_qg = plotToolsHT->GetContour(hexpCLs_NLO_qg,3,0);
1160     gexpCLs_NLO_qg->SetLineStyle(2); gexpCLs_NLO_qg->SetLineColor(1);gexpCLs_NLO_qg->SetLineWidth(2);
1161     TGraph * gobsCLs_NLO_qg = plotToolsHT->GetContour(hobsCLs_NLO_qg,3,0);
1162     gobsCLs_NLO_qg->SetLineStyle(1); gobsCLs_NLO_qg->SetLineColor(1);gobsCLs_NLO_qg->SetLineWidth(2);
1163     TGraph * gexp1p_NLO_qg = plotToolsHT->GetContour(hexp1p_NLO_qg,3,0);
1164     TGraph * gexp1m_NLO_qg = plotToolsHT->GetContour(hexp1m_NLO_qg,3,0);
1165     Smooth(gexpFC_LO_qg, 11);
1166     Smooth(gexpCLs_LO_qg, 11);
1167     Smooth(gexpCLs_NLO_qg, 11);
1168     Smooth(gobsCLs_NLO_qg, 11);
1169     Smooth(gexp1p_NLO_qg, 11);
1170     Smooth(gexp1m_NLO_qg, 11);
1171    
1172     TGraph * g1Sigma_qg = MakeBand(gexp1p_NLO_qg, gexp1m_NLO_qg);
1173     g1Sigma_qg->SetFillColor(5);
1174     h_qg->GetYaxis()->SetTitleOffset(1.5);
1175     h_qg->Draw("");
1176     TGraph * sqLEP = sq_LEP(); sqLEP->Draw("f");
1177     TGraph * glTEV = gl_TEV(); glTEV->Draw("f");
1178     TGraph * glCDF = gl_CDF(); glCDF->Draw("f");
1179     TGraph * glDEZ = gl_DEZ(); glDEZ->Draw("f");
1180    
1181     g1Sigma_qg->Draw("f");
1182     gexpFC_LO_qg->Draw("l");
1183     gexpCLs_LO_qg->Draw("l");
1184     gexpCLs_NLO_qg->Draw("l");
1185     gobsCLs_NLO_qg->Draw("l");
1186     //TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
1187     b.DrawLatex( 800, 50 ,"Tevatron RunI");
1188     b.DrawLatex( 800, 200,"CDF RunII");
1189     b.DrawLatex( 800, 240,"Dzero RunII");
1190 auterman 1.20 b.DrawLatex( 100, 600,"LEP #tilde{q}#tilde{q}");
1191 auterman 1.19 //
1192     TLegend * l_qg = new TLegend(0.45,0.7,0.85,0.89);
1193     l_qg->SetBorderSize(0);l_qg->SetFillColor(0);l_qg->SetTextFont(42);
1194 auterman 1.20 l_qg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1195 auterman 1.19 TGraph * expL_qg = (TGraph*)gexpCLs_NLO_qg->Clone();expL_qg->SetFillStyle(g1Sigma_qg->GetFillStyle());expL_qg->SetFillColor(g1Sigma_qg->GetFillColor());
1196     if (expL_qg) l_qg->AddEntry(expL_qg, "Expected (CLs, NLO)","lf");
1197     if (gobsCLs_NLO_qg) l_qg->AddEntry(gobsCLs_NLO_qg,"Observed (CLs, NLO)","l");
1198     if (gexpCLs_LO_qg) l_qg->AddEntry(gexpCLs_LO_qg, "Expected (CLs, LO)","l");
1199     if (gexpFC_LO_qg) l_qg->AddEntry(gexpFC_LO_qg, "Expected (FC, LO)","l");
1200     l_qg->Draw();
1201     gPad->RedrawAxis();
1202 auterman 1.10 c1->SaveAs("results/Exclusion_mSql_mGl_tb10.pdf");
1203 auterman 1.19
1204 auterman 1.20 //====================================================================================================================
1205     // Exclusion in gluino - squark mass ----------------------------------------
1206 auterman 1.21 c1->SetTopMargin ( 0.05 );
1207 auterman 1.20 TH2F*h_gq = new TH2F("hExcl_gq",";m_{#tilde{g}} [GeV]; m_{#tilde{q}} [GeV]; 95% CL Observed Exclusion",
1208     50,0,900,60,0,1000);
1209     TH2F*hExcl_gq = new TH2F("h_gq","",50,200,1200,60,200,1400);
1210 auterman 1.21 /* TH2F*hexpFC_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
1211 auterman 1.20 plotToolsHT->Area(hexpFC_HT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
1212     TH2F*hexpFC_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1213     plotToolsHT->Area(hexpFC_HT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
1214     TH2F*hexpFC_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
1215     plotTools->Area(hexpFC_MHT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
1216     TH2F*hexpFC_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1217     plotTools->Area(hexpFC_MHT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
1218 auterman 1.21 */ //
1219 auterman 1.20 TH2F*hexpCLs_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
1220     plotToolsHT->Area(hexpCLs_HT_LO_gq, MGluino, MSquarkL, ExpExclCL);
1221     TH2F*hexpCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1222     plotToolsHT->Area(hexpCLs_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
1223     TH2F*hexpCLs_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
1224     plotTools->Area(hexpCLs_MHT_LO_gq, MGluino, MSquarkL, ExpExclCL);
1225     TH2F*hexpCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1226     plotTools->Area(hexpCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
1227     //
1228     TH2F*hexp1p_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1229     plotToolsHT->Area(hexp1p_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
1230     TH2F*hexp1p_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1231     plotTools->Area(hexp1p_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
1232     TH2F*hexp1m_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1233     plotToolsHT->Area(hexp1m_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
1234     TH2F*hexp1m_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1235     plotTools->Area(hexp1m_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
1236    
1237     TH2F*hobsCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1238     plotToolsHT->Area(hobsCLs_HT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
1239     TH2F*hobsCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1240     plotTools->Area(hobsCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
1241 auterman 1.21 TH2F*hobsCLs_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
1242     plotToolsHT->Area(hobsCLs_HT_LO_gq, MGluino, MSquarkL, ObsExclCL);
1243     TH2F*hobsCLs_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
1244     plotTools->Area(hobsCLs_MHT_LO_gq, MGluino, MSquarkL, ObsExclCL);
1245 auterman 1.20 //
1246 auterman 1.21 //TH2 * hexpFC_LO_gq = plotToolsHT->BinWiseOr(hexpFC_HT_LO_gq, hexpFC_MHT_LO_gq);
1247 auterman 1.20 TH2 * hexpCLs_LO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_gq, hexpCLs_MHT_LO_gq);
1248     TH2 * hexpCLs_NLO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_gq, hexpCLs_MHT_NLO_gq);
1249     TH2 * hobsCLs_NLO_gq = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_gq, hobsCLs_MHT_NLO_gq);
1250 auterman 1.21 TH2 * hobsCLs_LO_gq = plotToolsHT->BinWiseOr(hobsCLs_HT_LO_gq, hobsCLs_MHT_LO_gq);
1251 auterman 1.20 TH2 * hexp1p_NLO_gq = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_gq, hexp1p_MHT_NLO_gq);
1252     TH2 * hexp1m_NLO_gq = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_gq, hexp1m_MHT_NLO_gq);
1253     //
1254 auterman 1.21 //TGraph * gexpFC_LO_gq = plotToolsHT->GetContour(hexpFC_LO_gq,3,0);
1255     //gexpFC_LO_gq->SetLineStyle(4); gexpFC_LO_gq->SetLineColor(3);gexpFC_LO_gq->SetLineWidth(2);
1256 auterman 1.20 TGraph * gexpCLs_LO_gq = plotToolsHT->GetContour(hexpCLs_LO_gq,3,0);
1257     gexpCLs_LO_gq->SetLineStyle(4); gexpCLs_LO_gq->SetLineColor(1);gexpCLs_LO_gq->SetLineWidth(2);
1258     TGraph * gexpCLs_NLO_gq = plotToolsHT->GetContour(hexpCLs_NLO_gq,3,0);
1259     gexpCLs_NLO_gq->SetLineStyle(2); gexpCLs_NLO_gq->SetLineColor(1);gexpCLs_NLO_gq->SetLineWidth(2);
1260     TGraph * gobsCLs_NLO_gq = plotToolsHT->GetContour(hobsCLs_NLO_gq,3,0);
1261     gobsCLs_NLO_gq->SetLineStyle(1); gobsCLs_NLO_gq->SetLineColor(1);gobsCLs_NLO_gq->SetLineWidth(2);
1262 auterman 1.21 TGraph * gobsCLs_LO_gq = plotToolsHT->GetContour(hobsCLs_LO_gq,3,0);
1263     gobsCLs_LO_gq->SetLineStyle(1); gobsCLs_LO_gq->SetLineColor(2);gobsCLs_LO_gq->SetLineWidth(2);
1264 auterman 1.20 TGraph * gexp1p_NLO_gq = plotToolsHT->GetContour(hexp1p_NLO_gq,3,0);
1265     TGraph * gexp1m_NLO_gq = plotToolsHT->GetContour(hexp1m_NLO_gq,3,0);
1266     TGraph * ch1 = (TGraph*)gobsCLs_NLO_gq->Clone();
1267 auterman 1.21 //Smooth2D(gexpFC_LO_gq, 11);
1268     //Smooth2D(gexpCLs_LO_gq, 11);
1269 auterman 1.20 Smooth2D(gexpCLs_NLO_gq, 11);
1270     Smooth2D(gobsCLs_NLO_gq, 11);
1271 auterman 1.21 Smooth2D(gobsCLs_LO_gq, 11);
1272 auterman 1.20 Smooth2D(gexp1p_NLO_gq, 11);
1273     Smooth2D(gexp1m_NLO_gq, 11);
1274     //Close2D(gobsCLs_NLO_gq);
1275     TGraph * g1Sigma_gq = MakeBand(gexp1p_NLO_gq, gexp1m_NLO_gq);
1276     g1Sigma_gq->SetFillColor(5);
1277     h_gq->GetYaxis()->SetTitleOffset(1.5);
1278     h_gq->GetXaxis()->SetTitleOffset(0.98);
1279     h_gq->Draw("");
1280     TGraph * sqDEZ = sq_DEZ(); sqDEZ->Draw("f");
1281     TGraph * sqCDF = sq_CDF(); sqCDF->Draw("f");
1282     TGraph * sqTEV = sq_TEV(); sqTEV->Draw("f");
1283     TGraph * glLEP = gl_LEP(); glLEP->Draw("f");
1284     TGraph * noSol = glsq_NoSol(); noSol->Draw("f"); noSol->Draw("l");
1285    
1286     ms.DrawLatex(10,1013,"L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1287     ms.DrawLatex(660,1013,"tan#beta=10, #mu>0, A_{0}=0");
1288     g1Sigma_gq->Draw("f");
1289     //ch1->Draw("l");
1290     // gexpFC_LO_gq->Draw("l");
1291 auterman 1.21 //gexpCLs_LO_gq->Draw("l");
1292 auterman 1.20 gexpCLs_NLO_gq->Draw("l");
1293     gobsCLs_NLO_gq->Draw("l");
1294 auterman 1.21 gobsCLs_LO_gq->Draw("l");
1295     TGraph * expL_gq = (TGraph*)gexpCLs_NLO_gq->Clone(); expL_gq->SetFillStyle(g1Sigma_gq->GetFillStyle());expL_gq->SetFillColor(g1Sigma_gq->GetFillColor());
1296 auterman 1.20 TLegend* legexp_gq = new TLegend(0.2,0.62,0.50,0.93,NULL,"brNDC");
1297     legexp_gq->SetFillColor(0);legexp_gq->SetShadowColor(0);//legexp_gq->SetFillStyle(4000);
1298     legexp_gq->SetTextFont(42);legexp_gq->SetTextSize(0.025);legexp_gq->SetBorderSize(1);
1299     //TEV_sg_cdf.SetLineColor(1);
1300 auterman 1.21 legexp_gq->SetHeader("CMS");
1301     if (gobsCLs_NLO_gq) legexp_gq->AddEntry(gobsCLs_NLO_gq,"Observed, NLO","l");
1302     if (gobsCLs_LO_gq) legexp_gq->AddEntry(gobsCLs_LO_gq,"Observed, LO","l");
1303     if (expL_gq) legexp_gq->AddEntry(expL_gq, "Expected #pm 1#sigma, NLO","lf");
1304     //if (gexpCLs_LO_gq) legexp_gq->AddEntry(gexpCLs_LO_gq, "Expected, LO","l");
1305 auterman 1.20 legexp_gq->AddEntry(sqTEV,"Tevatron RunI","f");
1306     legexp_gq->AddEntry(sqCDF,"CDF RunII","f");
1307 auterman 1.21 legexp_gq->AddEntry(sqDEZ,"D0 RunII","f");
1308 auterman 1.20 legexp_gq->AddEntry(glLEP,"LEP2","f"); //NOT FOR tb=50!
1309     legexp_gq->Draw();
1310     //TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
1311     //b.DrawLatex( 800, 50,"LEP #tilde{q}#tilde{q}");
1312     //b.DrawLatex( 50, 1000 ,"Tevatron RunI");
1313     //b.DrawLatex( 200, 1100,"CDF RunII");
1314     //b.DrawLatex( 240, 1200,"Dzero RunII");
1315 auterman 1.21
1316 auterman 1.20 TPaveText tpave(0.6,0.33,0.75,0.4,"NDC");
1317 auterman 1.21 tpave.AddText("no CMSSM");tpave.AddText("solution");
1318 auterman 1.20 tpave.SetBorderSize(0);tpave.SetFillColor(0);tpave.SetTextFont(42);
1319     tpave.Draw();
1320     //
1321     gPad->RedrawAxis();
1322     c1->SaveAs("results/Exclusion_mGl_mSql_tb10.pdf");
1323     //====================================================================================================================
1324 auterman 1.21
1325     // TestContours in gluino - squark mass
1326     TH2F*texcl_gq = new TH2F("texcl_gq",";m_{#tilde{g}} [GeV]; m_{#tilde{q}} [GeV]; 95% CL Observed Exclusion",
1327     50,200,1200,60,200,1400);
1328     TH2F*area_gq=(TH2F*)texcl_gq->Clone();
1329     plotTools->Area(area_gq, MGluino, MSquarkL, NLOObsExclCL);
1330     std::vector<TGraph*> contours_qg = plotTools->GetContours(hobsCLs_NLO_gq,3);
1331     area_gq->Draw("colz");
1332     //hobsCLs_MHT_NLO_gq->Draw("colz,same");
1333     int col_gl=kBlue-10;
1334     for (std::vector<TGraph*>::iterator cont=contours_qg.begin(); cont!=contours_qg.end(); ++cont){
1335     if (! *cont) continue;
1336     double x, y;
1337     (*cont)->GetPoint(0, x, y);
1338     (*cont)->SetLineColor(col_gl);(*cont)->SetLineWidth(3);
1339     (*cont)->Draw("l");
1340     TLatex l; l.SetTextSize(0.04); l.SetTextColor(col_gl++);
1341     char val[20];
1342     sprintf(val,"%d",(int)(cont-contours_qg.begin()));
1343     l.DrawLatex(x,y,val);
1344     if (cont-contours_qg.begin()>3) break;
1345     }
1346     c1->SaveAs("results/ExclusionTestContours_mGl_mSq_tb10.pdf");
1347 auterman 1.16
1348 auterman 1.21 //====================================================================================================================
1349     // Exclusion in gluino - squark mass ----------------------------------------
1350     //
1351     TH2F*hexpNoSCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1352     plotToolsHT->Area(hexpNoSCLs_HT_NLO_gq, MGluino, MSquarkL, NLOExpNoSigExclCL);
1353     TH2F*hexpNoSCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
1354     plotTools->Area(hexpNoSCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOExpNoSigExclCL);
1355     //
1356     TH2 * hexpNoSCLs_NLO_gq = plotToolsHT->BinWiseOr(hexpNoSCLs_HT_NLO_gq, hexpNoSCLs_MHT_NLO_gq);
1357     TGraph * gexpNoSCLs_NLO_gq = plotToolsHT->GetContour(hexpNoSCLs_NLO_gq,3,0);
1358     gexpNoSCLs_NLO_gq->SetLineStyle(2); gexpNoSCLs_NLO_gq->SetLineColor(1);gexpNoSCLs_NLO_gq->SetLineWidth(2);
1359     TGraph * gexpNoS1p_NLO_gq = plotToolsHT->ModifyExpSigmaY(gexp1p_NLO_gq, gexpCLs_NLO_gq, gexpNoSCLs_NLO_gq);
1360     TGraph * gexpNoS1m_NLO_gq = plotToolsHT->ModifyExpSigmaY(gexp1m_NLO_gq, gexpCLs_NLO_gq, gexpNoSCLs_NLO_gq);
1361     TGraph * ch2 = (TGraph*)gobsCLs_NLO_gq->Clone();
1362     h_gq->Draw("");
1363     Smooth(gexpNoS1p_NLO_gq, 11);
1364     Smooth(gexpNoS1m_NLO_gq, 11);
1365     Smooth2D(gexpNoSCLs_NLO_gq, 11);
1366     TGraph * g1SigmaNoS_gq = MakeBand(gexpNoS1p_NLO_gq, gexpNoS1m_NLO_gq);
1367     g1SigmaNoS_gq->SetFillColor(5);
1368     sqDEZ->Draw("f");
1369     sqCDF->Draw("f");
1370     sqTEV->Draw("f");
1371     glLEP->Draw("f");
1372     noSol->Draw("f"); noSol->Draw("l");
1373     ms.DrawLatex(10,1013,"L_{int} = 36/pb, #sqrt{s} = 7 TeV");
1374     ms.DrawLatex(660,1013,"tan#beta=10, #mu>0, A_{0}=0");
1375     g1SigmaNoS_gq->Draw("f");
1376     //gexpNoS1p_NLO_gq->Draw("l");
1377     //gexpNoS1m_NLO_gq->Draw("l");
1378     gexpNoSCLs_NLO_gq->Draw("l");
1379     //gexpCLs_NLO_gq->Draw("l");
1380     gobsCLs_NLO_gq->Draw("l");
1381     gobsCLs_LO_gq->Draw("l");
1382     legexp_gq->Draw();
1383     tpave.Draw();
1384     gPad->RedrawAxis();
1385     c1->SaveAs("results/Exclusion_mGl_mSql_tb10_NoSigHypPseudoData.pdf");
1386     //====================================================================================================================
1387    
1388    
1389    
1390    
1391    
1392    
1393 auterman 1.16 TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
1394 auterman 1.19 60,200,1400,50,200,1000);
1395 auterman 1.18 TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 2,2);
1396     TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,4);
1397     TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 2,1);
1398     TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
1399     TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
1400     TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2);
1401     TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4);
1402     TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1);
1403     TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
1404     TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
1405     Smooth(gexpexcl_qg, 11);
1406     Smooth(gexpexcl_qg_LO, 11);
1407     Smooth(gobsexcl_qg, 11);
1408     Smooth(gCLsExpExclHT_qg, 11);
1409     Smooth(gCLsExpExclHT_qg_LO, 11);
1410     Smooth(gCLsObsExclHT_qg, 11);
1411     Smooth(gCLsExpExclHTm1_qg, 11);
1412     Smooth(gCLsExpExclHTp1_qg, 11);
1413 auterman 1.16 TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
1414     hexcl_qg->Draw("");
1415     //
1416     if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
1417 auterman 1.19 if (gobsexcl_qg) gobsexcl_qg->Draw("l");
1418 auterman 1.16 if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
1419 auterman 1.19 if (gexpexcl_qg) gexpexcl_qg->Draw("l");
1420 auterman 1.16 if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
1421 auterman 1.18 if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l");
1422     if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
1423 auterman 1.19
1424 auterman 1.16 TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
1425     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
1426 auterman 1.18 TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
1427     if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
1428     if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
1429     if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
1430     if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
1431     if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg,"NLO Expected (HT, CLs)","lf");
1432     if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","lf");
1433 auterman 1.16 //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
1434     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
1435     leg_qg->Draw();
1436 auterman 1.19 c1->SaveAs("results/Exclusion_mSql_mGl_tb10_details.pdf");
1437 auterman 1.10 //c1->SaveAs("plot_tb10.pdf");
1438 auterman 1.2
1439 auterman 1.14 c1->SetLogy(1);
1440     c1->SetLogx(1);
1441     TH2F*hUncertScan = new TH2F("hUncertScan",";relative signal uncertainty; number of signal events",
1442     10,0.1,0.6,10,0.9,200);
1443     //hUncertScan->GetXaxis()->SetNdivisions(505, false);
1444     //hUncertScan->SetMinimum(1);
1445     TGraph * gCLsExpUncert = new TGraph(0);
1446     TGraph * gCLsObsUncert = new TGraph(0);
1447     TGraph * gPLExpUncert = new TGraph(0);
1448     TGraph * gPLObsUncert = new TGraph(0);
1449     TGraph * gFCExpUncert = new TGraph(0);
1450     TGraph * gFCObsUncert = new TGraph(0);
1451     TGraph * gMCMCExpUncert= new TGraph(0);
1452     TGraph * gMCMCObsUncert= new TGraph(0);
1453     plotTools->Graph(gCLsExpUncert, SignalRelUncertainty, ExpNSignLimit,1);
1454     plotTools->Graph(gCLsObsUncert, SignalRelUncertainty, ObsNSignLimit,1);
1455     plotTools->Graph(gPLExpUncert, SignalRelUncertainty, PLExpNSignLimit,1);
1456     plotTools->Graph(gPLObsUncert, SignalRelUncertainty, PLObsNSignLimit, 1);
1457     plotTools->Graph(gFCExpUncert, SignalRelUncertainty, FCExpNSignLimit,1);
1458     plotTools->Graph(gFCObsUncert, SignalRelUncertainty, FCObsNSignLimit, 1);
1459     plotTools->Graph(gMCMCExpUncert,SignalRelUncertainty, MCMCExpNSignLimit,1);
1460     plotTools->Graph(gMCMCObsUncert,SignalRelUncertainty, MCMCObsNSignLimit, 1);
1461     gCLsExpUncert->SetLineStyle(8);
1462     gPLObsUncert->SetLineColor(2);
1463     gPLExpUncert->SetLineColor(2); gPLExpUncert->SetLineStyle(8);
1464     gFCObsUncert->SetLineColor(3);
1465     gFCExpUncert->SetLineColor(3); gFCExpUncert->SetLineStyle(8);
1466     gMCMCObsUncert->SetLineColor(4);
1467     gMCMCExpUncert->SetLineColor(4); gMCMCExpUncert->SetLineStyle(8);
1468     hUncertScan->Draw("h");
1469     gCLsExpUncert->Draw("l");
1470     gCLsObsUncert->Draw("l");
1471     //gPLObsUncert->Draw("l");
1472     gPLExpUncert->Draw("l");
1473     //gFCObsUncert->Draw("l");
1474     //gFCExpUncert->Draw("l");
1475     //gMCMCObsUncert->Draw("l");
1476     //gMCMCExpUncert->Draw("l");
1477     TLatex l; l.SetTextSize(0.03); l.SetTextColor(1);
1478     l.DrawLatex(0.025,25,"excluded");
1479     TMarker m; m.SetMarkerStyle(29); m.SetMarkerColor(2);
1480     m.DrawMarker(0.128, 28.5);l.DrawLatex(0.14,28.5,"LM1");
1481     m.DrawMarker(0.147, 77.4);l.DrawLatex(0.16,77.4,"LM0");
1482    
1483     TLegend * leg_S = new TLegend(0.5,0.14,0.85,0.38);
1484     leg_S->SetBorderSize(0);leg_S->SetFillColor(0);
1485     leg_S->SetHeader("95% CL upper limits");
1486     leg_S->AddEntry(gCLsExpUncert,"Expected (CLs, TLimit)","l");
1487     leg_S->AddEntry(gCLsObsUncert,"Observed (CLs, TLimit)","l");
1488     //leg_S->AddEntry(gPLExpUncert,"Expected (PL, RooStat)","l");
1489     //leg_S->AddEntry(gPLObsUncert,"Observed (PL, RooStat)","l");
1490     //leg_S->AddEntry(gFCObsUncert,"Observed (FC, RooStat)","l");
1491     //leg_S->AddEntry(gFCExpUncert,"Expected (FC, RooStat)","l");
1492     //leg_S->AddEntry(gMCMCObsUncert,"Observed (MCMC, RooStat)","l");
1493     //leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
1494     leg_S->Draw();
1495     c1->SaveAs("results/UncertaintyScan.pdf");
1496 auterman 1.2
1497 auterman 1.21
1498    
1499    
1500    
1501    
1502    
1503    
1504    
1505 auterman 1.7 //theApp.Run();
1506 auterman 1.1 }
1507    
1508    
1509     int main(int argc, char** argv)
1510     {
1511     return plot(argc, argv);
1512     }