ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/auterman/SusyScan/PlotScript/plot3tb.cc

Comparing UserCode/auterman/SusyScan/PlotScript/plot3tb.cc (file contents):
Revision 1.2 by auterman, Wed Feb 9 09:32:31 2011 UTC vs.
Revision 1.8 by auterman, Fri May 20 07:39:34 2011 UTC

#	Line 4 | Line 4
4		#include "PlotTools.h"
5		#include "TheLimits.h"
6		#include "GlobalFunctions.h"
7	+	#include "StyleSettings.h"
8
9		#include "TRint.h"
10		#include "TROOT.h"
#	Line 38 | Line 39 | int plot(int argc, char** argv)
39		}
40
41		//Style stuff
42	+	/*
43		gStyle->SetHistFillColor(0);
44		gStyle->SetPalette(1);
45		gStyle->SetCanvasColor(0);
#	Line 68 | Line 70 | int plot(int argc, char** argv)
70		gStyle->SetLabelSize(0.03,"X");
71		gStyle->SetLabelSize(0.03,"Y");
72		gStyle->SetLabelSize(0.03,"Z");
73	+	*/
74	+	util::StyleSettings::paperNoTitle();
75	+	gStyle->SetPadBottomMargin(0.18);
76
77		//gROOT->SetStyle("MyStyle");
78		TCanvas * c1 = new TCanvas("c1","c1",600,600);
79	<	c1->SetFillStyle   ( 4000 );
80	<	c1->SetLeftMargin  ( 0.15 );
81	<	c1->SetRightMargin ( 0.15 );
82	<	c1->SetBottomMargin( 0.10 );
79	>	//c1->SetFillStyle   ( 4000 );
80	>	//c1->SetLeftMargin  ( 0.15 );
81	>	//c1->SetRightMargin ( 0.15 );
82	>	//c1->SetBottomMargin( 0.10 );
83		c1->cd();
84
85		//Get limits and generator masses ---------------------------------------------------
86		TheLimits * genpoints = new TheLimits();
87		//genpoints->Fill(argc, argv);
88	<	genpoints->Fill("limits_MHT_tb3/filelist.txt");
88	>	//genpoints->Fill("limits_MHT_tb3_20110303/filelist.txt");
89	>	genpoints->Fill("limits_MHT_tb3add/filelist_tb3.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92	<	genpointsHT->Fill("limits_HT_tb3/filelist.txt");
92	>	//genpointsHT->Fill("limits_HT_tb3_20110303/filelist.txt");
93	>	genpointsHT->Fill("limits_HT_tb3add/filelist_tb3.txt");
94
95
96		//Replace read limits with specific numbers
97		//genpoints->OverwriteLimits("ABCD_MHT");
98
99	<	genpoints->FillGeneratorMasses("tb3_mu1_a0_massscan.dat");
99	>	genpoints->FillGeneratorMasses("GenScan_tb3.dat");
100		genpoints->match();
101	<	genpointsHT->FillGeneratorMasses("tb3_mu1_a0_massscan.dat");
101	>	genpointsHT->FillGeneratorMasses("GenScan_tb3.dat");
102		genpointsHT->match();
103
104		//the plotting ----------------------------------------------------------------------
#	Line 117 | Line 124 | int plot(int argc, char** argv)
124		/**/
125		// cross-section in M0 - M1/2
126		TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
127	<	50,0,509.9,25,100,350);
127	>	150,0,1509.9,25,100,350);
128		plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
129		hxsec->SetMinimum(0.01);
130		//sq500->Draw();
#	Line 129 | Line 136 | int plot(int argc, char** argv)
136
137		// Observed Limit in M0 - M1/2
138		TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
139	<	50,0,509.9,25,100,350);
139	>	150,0,1509.9,25,100,350);
140		plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
141		hobslimit->SetMinimum(0.01);
142		hobslimit->Draw("colz");
#	Line 137 | Line 144 | int plot(int argc, char** argv)
144
145		// Expected Limit in M0 - M1/2
146		TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
147	<	50,0,509.9,25,100,350);
147	>	150,0,1509.9,25,100,350);
148		plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
149		hexplimit->SetMinimum(0.01);
150		hexplimit->Draw("colz");
#	Line 145 | Line 152 | int plot(int argc, char** argv)
152
153		// Signal Acceptance in M0 - M1/2
154		TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
155	<	50,0,509.9,25,100,350);
155	>	150,0,1509.9,25,100,350);
156		plotTools->Area(hsigacc, Mzero, Mhalf, SignalAcceptance);
157		hsigacc->SetMinimum(0.01);
158		hsigacc->SetMaximum(1.0);
#	Line 159 | Line 166 | int plot(int argc, char** argv)
166		// Exp. Limit on Number of Signal Events in M0 - M1/2
167		c1->SetLogz(0);
168		TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
169	<	50,0,509.9,25,100,350);
169	>	150,0,1509.9,25,100,350);
170		plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
171		hexplimitnsig->SetMinimum(0.0);
172		hexplimitnsig->SetMaximum(20);
#	Line 168 | Line 175 | int plot(int argc, char** argv)
175
176		// Obs. Limit on Number of Signal Events in M0 - M1/2
177		TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
178	<	50,0,509.9,25,100,350);
178	>	150,0,1509.9,25,100,350);
179		plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
180		hobslimitnsig->SetMinimum(0.0);
181		hobslimitnsig->SetMaximum(20);
#	Line 178 | Line 185 | int plot(int argc, char** argv)
185		c1->SetLogz(0);
186		// Expected Exclusion in M0 - M1/2
187		TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
188	<	50,0,509.9,35,100,450);
188	>	150,0,1509.9,35,100,450);
189		plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
190		hexpexcl->Draw("colz");
191		c1->SaveAs("results_tb3/ExpExclusion_m0_m12_tb3.pdf");
192
193		// Observed Exclusion in M0 - M1/2
194		TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
195	<	50,0,509.9,35,100,450);
195	>	150,0,1509.9,35,100,450);
196		plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
197		hobsexcl->Draw("colz");
198		c1->SaveAs("results_tb3/ObsExclusion_m0_m12_tb3.pdf");
199
200		// Observed Exclusion in M0 - M1/2
201		TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
202	<	50,0,509.9,35,100,450);
202	>	150,0,1509.9,35,100,450);
203		plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
204		hPLobsexcl->Draw("colz");
205		c1->SaveAs("results_tb3/PL_ObsExclusion_m0_m12_tb3.pdf");
206
207	<	// TestContours in M0 - M1/2
207	>	// TestContours in M0 - M1/2  ===================================================================
208		TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
209	<	50,0,509.9,35,100,450);
209	>	150,0,1509.9,35,100,450);
210		TH2F*texpexcl=(TH2F*)texcl->Clone();
211	<	plotTools->Area(texpexcl, Mzero, Mhalf, MCMCExpExclusion);
212	<	TH2F*tobsexcl=(TH2F*)texcl->Clone();
206	<	plotTools->Area(tobsexcl, Mzero, Mhalf, MCMCObsExclusion);
207	<	std::vector<TGraph*> contours = plotTools->GetContours(texpexcl,3);
208	<	//std::vector<TGraph*> contours = plotTools->GetContours(tobsexcl,3);
209	<	//hPLexpexcl
210	<	//hexcl->Draw("colz");
211	<	//hexpexcl->Draw("colz");
212	<	//hobsexcl->Draw("colz");
211	>	plotToolsHT->Area(texpexcl, Mzero, Mhalf, NLOExpNoSigExclCL);
212	>	std::vector<TGraph*> contours = plotToolsHT->GetContours(texpexcl,3);
213		texpexcl->Draw("colz");
214		int col=kBlue-10;
215		for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
#	Line 222 | Line 222 | int plot(int argc, char** argv)
222		char val[20];
223		sprintf(val,"%d",(int)(cont-contours.begin()));
224		l.DrawLatex(x,y,val);
225	–	//if (cont-contours.begin()>3) break;
225		}
226	<	c1->SaveAs("results_tb3/ExclusionTestContours_m0_m12_tb3.pdf");
226	>	c1->SaveAs("results/ExclusionTestContours_m0_m12_tb3.pdf");
227
228
229		// Exclusion in M0 - M1/2
230		TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
231	<	50,0,509.9,35,100,450);
232	<	TGraph * gexpexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 2,2);
233	<	TGraph * gobsexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 2,1);
234	<	TGraph * gMCMCexpexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2);
235	<	TGraph * gMCMCobsexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1);
236	<	TGraph * gCLsExpExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 1,2);
237	<	TGraph * gCLsObsExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 1,1);
238	<	TGraph * gCLsExpExclHTm1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLm1sigma,3,0, 5,2);
239	<	TGraph * gCLsExpExclHTp1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLp1sigma,3,0, 5,2);
240	<	TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1);
241	<	hexcl->Draw("colz");
242	<	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
243	<	if (gexpexcl)         gexpexcl->Draw("l");
244	<	if (gCLsExpExclHT)    gCLsExpExclHT->Draw("l");
245	<	if (gMCMCexpexcl)     gMCMCexpexcl->Draw("l");
246	<	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
247	<	//if (gobsexcl)       gobsexcl->Draw("l");
248	<	//if (gPLobsexcl)     gPLobsexcl->Draw("l");
249	<	//if (gFCobsexcl)     gFCobsexcl->Draw("l");
250	<	//if (gMCMCobsexcl)   gMCMCobsexcl->Draw("l");
251	<	//if (gFCexpexcl)     gFCexpexcl->Draw("l");
252	<	//if (gPLexpexcl)     gPLexpexcl->Draw("l");
253	<	//if (gFCexpexcl)     gFCexpexcl->Draw("l");
231	>	150,0,1509.9,35,100,500);
232	>	TH2F*hs = new TH2F("hs","",150,0,1509.9,35,100,450);
233	>	TGraph * gexpexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,3, 2,2);
234	>	TGraph * gexpexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 2,4);
235	>	TGraph * gobsexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,3, 2,1);
236	>	TGraph * gobsexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
237	>	TGraph * gMCMCexpexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion,   3,0, 3,2);
238	>	TGraph * gMCMCobsexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion,   3,0, 3,1);
239	>	TGraph * gFCexpexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,  3,0, 3,2);
240	>	TGraph * gFCexpexclHT_LO  = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion,     3,0, 3,4);
241	>	TGraph * gFCexpexclMHT    = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,    3,0, 3,2);
242	>	TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion,       3,0, 3,4);
243	>	TGraph * gFCobsexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion,  3,0, 3,1);
244	>	TGraph * gFCobsexclMHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion,  3,0, 3,1);
245	>	TGraph * gCLsExpExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,3, 1,2);
246	>	TGraph * gCLsObsExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,3, 1,1);
247	>	TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 1,4);
248	>	TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
249	>	TGraph * gCLsExpExclHTm1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,2, 5,2);
250	>	TGraph * gCLsExpExclHTp1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,2, 5,2);
251	>	TGraph * gCLsExpExclMHTm1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,  3,2, 5,2);
252	>	TGraph * gCLsExpExclMHTp1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,  3,2, 5,2);
253	>	TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
254	>	TGraph * gCheck2 = (TGraph*)gexpexcl->Clone();
255	>	Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
256	>	Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
257	>	Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
258	>	Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
259	>	Smooth( gexpexcl, 27 );
260	>	Smooth( gobsexcl_LO, 27 );
261	>	//Smooth( gobsexcl, 27 );
262	>	Smooth( gobsexcl, 27 );
263	>	Smooth( gexpexcl_LO, 27 );
264	>	Smooth( gFCexpexclHT, 27 );
265	>	Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
266	>	Smooth( gFCexpexclMHT, 27 );
267	>	Smooth( gFCexpexclMHT_LO, 27 );
268	>	Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 2 );
269	>	Smooth( gFCobsexclMHT, 27 );
270	>	Smooth( gCLsExpExclHTm1, 27 );
271	>	Smooth( gCLsExpExclHTp1, 27 );
272	>	Smooth( gCLsExpExclMHTm1, 27 );
273	>	Smooth( gCLsExpExclMHTp1, 27 );
274	>	TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
275	>	TGraph * CLsObsLO  = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
276	>	TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
277	>	TGraph * CLsExpLO  = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
278	>	TGraph * FCExpLO   = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
279	>	TGraph * FCExpNLO  = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
280	>	TGraph * FCObsNLO  = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
281	>	TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
282	>	TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
283	>	TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclp1, gCLsExpExclm1);
284	>	hexcl->GetYaxis()->SetTitleOffset(1.3);
285	>	hexcl->GetXaxis()->SetTitleOffset(0.92);
286	>	hexcl->Draw("colz");
287		//set old exclusion Limits
288	<	TGraph* LEP_ch = set_lep_ch();
289	<	TGraph* LEP_sl = set_lep_sl();//slepton curve
290	<	TGraph* TEV_sg_cdf = set_tev_sg_cdf();//squark gluino cdf
291	<	TGraph* TEV_sg_d0 = set_tev_sg_d0();//squark gluino d0
292	<	TGraph* TEV_tlp_cdf = set_tev_tlp_cdf();//trilepton cdf
293	<	TGraph* TEV_tlp_d0 = set_tev_tlp_d0();//trilepton d0
294	<	TGraph* stau = set_tev_stau();//stau
295	<	TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1();
296	<	TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2();
288	>	TGraph* Atlas = Atlas_m0_m12_tb3_obs();
289	>	Atlas->SetLineColor( 28 );
290	>	TGraph* LEP_ch = set_lep_ch(3);
291	>	TGraph* LEP_sl = set_lep_sl(3);//slepton curve
292	>	TGraph* TEV_sg_cdf = set_tev_sg_cdf(3);//squark gluino cdf
293	>	TGraph* TEV_sg_d0 = set_tev_sg_d0(3);//squark gluino d0
294	>	//TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
295	>	//TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
296	>	TGraph* stau = set_tev_stau(3);//stau
297	>	TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(3);
298	>	TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(3);
299		TGraphErrors* First  = getObserved_NLO_tanBeta3();
300		TGraphErrors* Second = getExpected_NLO_tanBeta3();//getLO_jetMultis();
301		TGraphErrors* Third  = getLO_tanBeta3();
#	Line 271 | Line 305 | int plot(int argc, char** argv)
305		First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
306		TSpline3 *sFirst = new TSpline3("sFirst",First);
307		sFirst->SetLineColor(kBlue);
308	<	sFirst->SetLineWidth(3);
308	>	//sFirst->SetLineWidth(3);
309		TSpline3 *sSecond = new TSpline3("sSecond",Second);
310		sSecond->SetLineColor(kBlue);
311		sSecond->SetLineStyle(2);
312	<	sSecond->SetLineWidth(3);
312	>	//sSecond->SetLineWidth(3);
313		TSpline3 *sThird = new TSpline3("sThird",Third);
314	<	sThird->SetLineColor(kGreen+2);
314	>	sThird->SetLineColor(kBlue);
315		sThird->SetLineStyle(4);
316	<	sThird->SetLineWidth(3);
316	>	//sThird->SetLineWidth(3);
317	>	//gCheck1->Draw("same");
318	>	//gCheck2->Draw("same");
319	>	//gexpexcl->Draw("l");
320	>	//gobsexcl->Draw("l");
321	>	//gexpexcl_LO->Draw("l");
322	>	//gCLsExpExclHT_LO->Draw("l");
323	>	//gCLsObsExclHT->Draw("l");
324	>	//gCLsExpExclHT->Draw("l");
325	>	//gFCexpexclHT_LO->Draw("l");
326	>	//gFCexpexclMHT_LO->Draw("l");
327	>	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
328		TEV_sn_d0_1->Draw("fsame");
329	<	TEV_sn_d0_2->Draw("fsame");
330	<	TEV_sg_d0->Draw("fsame");
329	>	TEV_sn_d0_2->Draw("fsame"); //only for tb=3
330	>	stau->Draw("fsame");
331	>	LEP_ch->Draw("fsame");
332	>	LEP_sl->Draw("fsame");
333		TEV_sg_cdf->Draw("fsame");
334	+	TEV_sg_cdf->Draw("lsame");
335	+	TEV_sg_d0->Draw("fsame");
336	+	TEV_sg_d0->Draw("lsame");
337	+	sFirst->Draw("same");
338	+	TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
339	+	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
340	+
341	+	//sSecond->Draw("same");
342	+	//sThird->Draw("same");
343	+	Atlas->Draw("c,same");
344	+	TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
345	+	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
346	+	//gCheck1->Draw("same");
347	+	//gCheck2->Draw("same");
348	+	TLegend* legexp = new TLegend(0.66,0.65,0.98,0.93,NULL,"brNDC");
349	+	legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
350	+	//TEV_sg_cdf.SetLineColor(1);
351	+	legexp->SetHeader("CMS preliminary");
352	+	legexp->AddEntry(TEV_sg_cdf,"CDF  #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
353	+	legexp->AddEntry(TEV_sg_d0,"D0   #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{#mu<0}","f");
354	+	//ch_gr.SetLineColor(1);
355	+	legexp->AddEntry(LEP_ch,"LEP2   #tilde{#chi}_{1}^{#pm}","f");
356	+	//sl_gr.SetLineColor(1);
357	+	legexp->AddEntry(LEP_sl,"LEP2   #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
358	+	//if(tanbeta == 3)
359	+	legexp->AddEntry(TEV_sn_d0_1,"D0  #chi^{#pm}_{1}, #chi^{0}_{2}","f");
360	+	legexp->AddEntry(sFirst, "CMS #alpha_{T}");
361	+	legexp->AddEntry(Atlas,  "Atlas","l");
362	+	legexp->Draw();
363	+
364	+	CLsObsNLO->Draw("l");
365	+	CLsObsLO->Draw("l");
366	+	CLsExpNLO->Draw("l");
367	+	//FCExpNLO->Draw("l");
368	+
369	+	//constant ssqquark and gluino lines
370	+	TF1* lnsq[4];
371	+	TF1* lngl[4];
372	+	TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-14);sqt.SetTextColor(kGray+2);
373	+	sqt.DrawLatex(148,218,"#font[92]{#tilde{q}(500)GeV}");
374	+	sqt.DrawLatex(220,385,"#font[92]{#tilde{q}(800)GeV}");
375	+	TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
376	+	glt.DrawLatex(430,184,"#font[92]{#tilde{g}(500)GeV}");
377	+	//glt.DrawLatex(450,235,"#font[92]{#tilde{g}(650)GeV}");
378	+	glt.DrawLatex(428,311,"#font[92]{#tilde{g}(800)GeV}");
379	+	int tanBeta_=3;
380	+	for(int i = 0; i < 4; i++){
381	+	lnsq[i] = constant_squark(tanBeta_,i);
382	+	lngl[i] = constant_gluino(tanBeta_,i);
383	+	lngl[i]->Draw("same");
384	+	lnsq[i]->Draw("same");
385	+	}
386	+	TLegend * leg = new TLegend(0.3,0.8,0.65,0.93);
387	+	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
388	+	TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
389	+	leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
390	+	leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
391	+	leg->AddEntry(CLsObsLO, "Observed, LO","l");
392	+	leg->AddEntry(expLeg,   "Expected #pm 1#sigma, NLO","lf");
393	+	//leg->AddEntry(FCExpLO,  "Expected, FC, NLO","l");
394	+	leg->Draw();
395	+	gPad->RedrawAxis();
396	+	c1->SaveAs("results/Exclusion_m0_m12_tb3.pdf");
397	+
398	+	hexcl->Draw("");
399	+	gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb3");
400	+	CLsObsNLO->SetName("CLsNLO_ObservedTb3");
401	+	CLsExpNLO->SetName("CLsNLO_ExpectedTb3");
402	+	CLsExpLO->SetName("CLsLO_ExpectedTb3");
403	+	gCLsExp1Sigma->Draw("l");
404	+	CLsObsNLO->Draw("l");
405	+	CLsExpNLO->Draw("l");
406	+	CLsObsLO->Draw("l");
407	+	c1->SaveAs("results/LimitContours_tb3.C");
408	+
409	+	//   plotToolsHT->Print(Xsection, Mzero,Mhalf,gCLsObsExclHT, 10);
410	+
411	+
412	+
413	+	hexcl->Draw("colz");
414		LEP_ch->Draw("fsame");
415		LEP_sl->Draw("fsame");
416	+	TEV_sg_cdf->Draw("fsame");
417	+	TEV_sg_cdf->Draw("lsame");
418	+	TEV_sg_d0->Draw("fsame");
419	+	TEV_sg_d0->Draw("lsame");
420	+	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
421	+	sFirst->Draw("same");
422	+	Atlas->Draw("c,same");
423	+	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
424	+	CLsObsNLO->Draw("l");
425	+	//CLsObsLO->Draw("l");
426	+	CLsExpNLO->Draw("l");
427	+	//FCExpNLO->Draw("l");
428	+	//Jim's limits///////////////////////////////////////////////////////////////////////////////
429	+	TGraph* JimObsHT  = Jim_ht_tb3(0);
430	+	TGraph* JimObsMHT = Jim_mht_tb3(0);
431	+	TGraph* JimExpHT  = Jim_ht_tb3(1);
432	+	TGraph* JimExpMHT = Jim_mht_tb3(1);
433	+	TGraph* JimExpHTup  = Jim_ht_tb3(2);
434	+	TGraph* JimExpMHTup = Jim_mht_tb3(2);
435	+	TGraph* JimExpHTdn  = Jim_ht_tb3(3);
436	+	TGraph* JimExpMHTdn = Jim_mht_tb3(3);
437	+	Smooth(JimObsHT, 2);
438	+	Smooth(JimObsMHT,2);
439	+	Smooth(JimExpHT, 2);
440	+	Smooth(JimExpMHT,2);
441	+	Smooth(JimExpHTup, 2);
442	+	Smooth(JimExpMHTup,2);
443	+	Smooth(JimExpHTdn, 2);
444	+	Smooth(JimExpMHTdn,2);
445	+	TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
446	+	TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
447	+	TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
448	+	TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
449	+	TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
450	+	TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
451	+	JimExp1Sigma->Draw("f");
452	+	JimObs->Draw("c");
453	+	JimExp->Draw("c");
454		stau->Draw("fsame");
455	+	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
456	+	//constant ssqquark and gluino lines
457	+	sqt.DrawLatex(148,218,"#font[92]{#tilde{q}(500)GeV}");
458	+	sqt.DrawLatex(220,385,"#font[92]{#tilde{q}(800)GeV}");
459	+	glt.DrawLatex(430,184,"#font[92]{#tilde{g}(500)GeV}");
460	+	glt.DrawLatex(428,311,"#font[92]{#tilde{g}(800)GeV}");
461	+	for(int i = 0; i < 4; i++){
462	+	lngl[i]->Draw("same");
463	+	lnsq[i]->Draw("same");
464	+	}
465	+	legexp->Draw();
466	+	TLegend * legBayes = new TLegend(0.3,0.8,0.65,0.93);
467	+	legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
468	+	legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
469	+	legBayes->AddEntry(JimObs,   "Observed, Bayes","l");
470	+	legBayes->AddEntry(JimLeg,   "Expected #pm 1#sigma, Bayes","lf");
471	+	legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
472	+	legBayes->AddEntry(expLeg,   "Expected #pm 1#sigma, CLs","lf");
473	+	legBayes->Draw();
474	+	gPad->RedrawAxis();
475	+	c1->SaveAs("results/Exclusion_m0_m12_tb3_Bayes.pdf");
476	+
477	+
478	+	//The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
479	+	TGraph * gCLsExpNoSExclMHT   = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,2, 1,2);
480	+	TGraph * gCLsExpNoSExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,2, 1,2);
481	+	TGraph * gCLsExpNoSExclHTm1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
482	+	TGraph * gCLsExpNoSExclHTp1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
483	+	TGraph * gCLsExpNoSExclMHTm1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
484	+	TGraph * gCLsExpNoSExclMHTp1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
485	+	Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
486	+	Smooth( gCLsExpNoSExclHT, 27 );  gCLsExpNoSExclHT->SetLineWidth( 3 );
487	+	Smooth( gCLsExpNoSExclHTm1, 27 );
488	+	Smooth( gCLsExpNoSExclHTp1, 27 );
489	+	Smooth( gCLsExpNoSExclMHTm1, 27 );
490	+	Smooth( gCLsExpNoSExclMHTp1, 27 );
491	+	TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
492	+	TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
493	+	TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
494	+	TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
495	+	hexcl->Draw("colz");
496	+	LEP_ch->Draw("fsame");
497	+	LEP_sl->Draw("fsame");
498	+	TEV_sg_cdf->Draw("fsame");
499	+	TEV_sg_cdf->Draw("lsame");
500	+	TEV_sg_d0->Draw("fsame");
501	+	TEV_sg_d0->Draw("lsame");
502	+	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
503	+	gCLsExpNoS1Sigma->Draw("lf,same");
504		sFirst->Draw("same");
505	<	sSecond->Draw("same");
506	<	sThird->Draw("same");
507	<	TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
508	<	b.DrawLatex( 10,420,"#tilde{#tau} LSP");
509	<	b.DrawLatex(  6,150,"LEP2");
510	<	b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
511	<	b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
512	<	b.DrawLatex(100,140,"CDF");
513	<	b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
514	<	b.DrawLatex( 80,180,"D0 #tilde{#nu}");
515	<	TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
516	<	leg->SetBorderSize(0);leg->SetFillColor(0);
517	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
518	<	if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs, TLimit)","l");
519	<	if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs, TLimit)","l");
520	<	if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs, TLimit)","l");
521	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
308	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
309	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
310	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
311	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
312	<	if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
313	<	if (sFirst)  leg->AddEntry(sFirst, "NLO RA1 Observed");
314	<	if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected");
315	<	if (sThird) leg->AddEntry(sThird,"LO RA1 Observed");
505	>	Atlas->Draw("c,same");
506	>	CLsObsNLO->Draw("l,same");
507	>	CLsObsLO->Draw("l,same");
508	>	CLsExpNoSNLO->Draw("l,same");
509	>	//FCExpNLO->Draw("l");
510	>	stau->Draw("fsame");
511	>	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
512	>	//constant ssqquark and gluino lines
513	>	sqt.DrawLatex(148,218,"#font[92]{#tilde{q}(500)GeV}");
514	>	sqt.DrawLatex(220,385,"#font[92]{#tilde{q}(800)GeV}");
515	>	glt.DrawLatex(430,184,"#font[92]{#tilde{g}(500)GeV}");
516	>	glt.DrawLatex(428,311,"#font[92]{#tilde{g}(800)GeV}");
517	>	for(int i = 0; i < 4; i++){
518	>	lngl[i]->Draw("same");
519	>	lnsq[i]->Draw("same");
520	>	}
521	>	legexp->Draw();
522		leg->Draw();
523		gPad->RedrawAxis();
524	<	c1->SaveAs("results_tb3/Exclusion_m0_m12_tb3.pdf");
524	>	c1->SaveAs("results/Exclusion_m0_m12_tb3_NoSigHypPseudoData.pdf");
525	>
526	>
527
528		// Signal Contamination in M0 - M1/2
529		c1->SetLogz(0);
530		TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
531	<	50,0,509.9,25,100,350);
531	>	150,0,1509.9,25,100,350);
532		plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
533		hsigcont->SetMinimum(0.01);
534		hsigcont->SetMaximum(20);
#	Line 331 | Line 539 | int plot(int argc, char** argv)
539		// Signal Contamination in M0 - M1/2
540		c1->SetLogz(0);
541		TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
542	<	50,0,509.9,25,100,350);
542	>	150,0,1509.9,25,100,350);
543		plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
544		hsigcontMHT->SetMinimum(0.01);
545		hsigcontMHT->SetMaximum(20);
#	Line 339 | Line 547 | int plot(int argc, char** argv)
547		if (gobsexcl)    gobsexcl->Draw("l");
548		c1->SaveAs("results_tb3/SignalContamination_MHT_m0_m12_tb3.pdf");
549
550	+
551	+	// Signal JEC Uncertainty  MHT
552	+	c1->SetLogz(0);
553	+	TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
554	+	150,0,1509.9,25,100,350);
555	+	plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
556	+	hsig_jec_mht->SetMinimum(0.0);
557	+	hsig_jec_mht->SetMaximum(0.25);
558	+	hsig_jec_mht->SetContour(10);
559	+	hsig_jec_mht->Draw("colz");
560	+	if (gobsexcl)    gobsexcl->Draw("l");
561	+	c1->SaveAs("results_tb3/SigJEC_MHT_m0_m12_tb3.pdf");
562	+
563	+	// Signal JEC Uncertainty  HT
564	+	c1->SetLogz(0);
565	+	TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
566	+	150,0,1509.9,25,100,350);
567	+	plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
568	+	hsig_jec_ht->SetMinimum(0.0);
569	+	hsig_jec_ht->SetMaximum(0.25);
570	+	hsig_jec_ht->SetContour(10);
571	+	hsig_jec_ht->Draw("colz");
572	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
573	+	c1->SaveAs("results_tb3/SigJEC_HT_m0_m12_tb3.pdf");
574	+
575	+
576	+	// Signal MuIso Uncertainty  MHT
577	+	c1->SetLogz(0);
578	+	TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
579	+	150,0,1509.9,25,100,350);
580	+	plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
581	+	hsig_MuIso_mht->SetMinimum(0.0);
582	+	hsig_MuIso_mht->SetMaximum(0.05);
583	+	hsig_MuIso_mht->SetContour(5);
584	+	hsig_MuIso_mht->Draw("colz");
585	+	if (gobsexcl)    gobsexcl->Draw("l");
586	+	c1->SaveAs("results_tb3/SigMuIso_MHT_m0_m12_tb3.pdf");
587	+
588	+	// Signal MuIso Uncertainty  HT
589	+	c1->SetLogz(0);
590	+	TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
591	+	150,0,1509.9,25,100,350);
592	+	plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
593	+	hsig_MuIso_ht->SetMinimum(0.0);
594	+	hsig_MuIso_ht->SetMaximum(0.05);
595	+	hsig_MuIso_ht->SetContour(5);
596	+	hsig_MuIso_ht->Draw("colz");
597	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
598	+	c1->SaveAs("results_tb3/SigMuIso_HT_m0_m12_tb3.pdf");
599	+
600
601	<	/*
601	>	// Signal kFactor Uncertainty  MHT
602	>	c1->SetLogz(0);
603	>	TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
604	>	150,0,1509.9,25,100,350);
605	>	plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
606	>	hsig_kFactor_mht->SetMinimum(0.1);
607	>	hsig_kFactor_mht->SetMaximum(0.2);
608	>	hsig_kFactor_mht->SetContour(10);
609	>	hsig_kFactor_mht->Draw("colz");
610	>	if (gobsexcl)    gobsexcl->Draw("l");
611	>	c1->SaveAs("results_tb3/SigkFactor_MHT_m0_m12_tb3.pdf");
612	>
613	>	// Signal kFactor Uncertainty  HT
614	>	c1->SetLogz(0);
615	>	TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
616	>	150,0,1509.9,25,100,350);
617	>	plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
618	>	hsig_kFactor_ht->SetMinimum(0.1);
619	>	hsig_kFactor_ht->SetMaximum(0.2);
620	>	hsig_kFactor_ht->SetContour(10);
621	>	hsig_kFactor_ht->Draw("colz");
622	>	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
623	>	c1->SaveAs("results_tb3/SigkFactor_HT_m0_m12_tb3.pdf");
624	>
625		//-----------------------------------------------------------------------------------
626		c1->SetLogz(1);
627
#	Line 349 | Line 630 | int plot(int argc, char** argv)
630
631		// cross-section in squark - gluino mass
632		TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
633	<	60,200,1400,50,200,1200);
633	>	60,200,800,40,200,900);
634		plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
635		hxsec_qg->SetMinimum(0.01);
636		hxsec_qg->Draw("colz");
#	Line 359 | Line 640 | int plot(int argc, char** argv)
640
641		// Observed Limit in squark - gluino mass
642		TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
643	<	60,200,1400,50,200,1200);
643	>	60,200,800,40,200,900);
644		plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
645		hobslimit_qg->SetMinimum(0.01);
646		hobslimit_qg->Draw("colz");
#	Line 367 | Line 648 | int plot(int argc, char** argv)
648
649		// Expected Limit in squark - gluino mass
650		TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
651	<	60,200,1400,50,200,1200);
651	>	60,200,800,40,200,900);
652		plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
653		hexplimit_qg->SetMinimum(0.01);
654		hexplimit_qg->Draw("colz");
#	Line 375 | Line 656 | int plot(int argc, char** argv)
656
657		// Signal Acceptance in squark - gluino mass
658		TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
659	<	60,200,1400,50,200,1200);
659	>	60,200,800,40,200,900);
660		plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
661		hsigacc_qg->SetMinimum(0.01);
662		hsigacc_qg->SetMaximum(1.0);
#	Line 385 | Line 666 | int plot(int argc, char** argv)
666		// Exp. Limit on Number of Signal Events in squark - gluino mass
667		c1->SetLogz(0);
668		TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
669	<	60,200,1400,50,200,1200);
669	>	60,200,800,40,200,900);
670		plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
671		hexplimitnsig_qg->SetMinimum(0.0);
672		hexplimitnsig_qg->SetMaximum(20);
#	Line 394 | Line 675 | int plot(int argc, char** argv)
675
676		// Obs. Limit on Number of Signal Events in squark - gluino mass
677		TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
678	<	60,200,1400,50,200,1200);
678	>	60,200,800,40,200,900);
679		plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
680		hobslimitnsig_qg->SetMinimum(0.0);
681		hobslimitnsig_qg->SetMaximum(20);
#	Line 404 | Line 685 | int plot(int argc, char** argv)
685		c1->SetLogz(0);
686		// Expected Exclusion in squark - gluino mass
687		TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
688	<	60,200,1400,50,200,1200);
688	>	60,200,800,40,200,900);
689		plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
690		hexpexcl_qg->Draw("colz");
691		c1->SaveAs("results_tb3/ExpExclusion_mSql_mGl_tb3.pdf");
692
693		// Observed Exclusion in squark - gluino mass
694		TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
695	<	60,200,1400,50,200,1200);
695	>	60,200,800,40,200,900);
696		plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
697		//std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
698		hobsexcl_qg->Draw("colz");
#	Line 423 | Line 704 | int plot(int argc, char** argv)
704
705		// TestContours in M0 - M1/2
706		TH2F*texcl_qg = new TH2F("texcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
707	<	60,200,1400,50,200,1200);
707	>	60,200,800,40,200,900);
708		TH2F*tPLobsexcl_qg=(TH2F*)texcl_qg->Clone();
709		plotTools->Area(tPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
710		TH2F*tPLexpexcl_qg=(TH2F*)texcl_qg->Clone();
#	Line 454 | Line 735 | int plot(int argc, char** argv)
735
736		// Observed Exclusion in squark - gluino mass
737		TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
738	<	60,200,1400,50,200,1200);
738	>	60,200,800,40,200,900);
739		plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
740		//hPLobsexcl_qg->Draw("colz");
741		//c1->SaveAs("results_tb3/PL_ObsExclusion_mSql_mGl_tb3.pdf");
742
743
744
464	–
465	–
466	–	// Exclusion in squark - gluino mass ----------------------------------------
745		TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
746	<	60,200,1400,50,200,1200);
747	<	TGraph * gexpexcl_qg = plotTools->GetContour(hexpexcl_qg,3,0);
748	<	if (gexpexcl_qg) gexpexcl_qg->SetLineStyle(2);
749	<	TGraph * gobsexcl_qg = plotTools->GetContour(hobsexcl_qg,3,0);
750	<	TH2F*hPLexpexcl_qg=(TH2F*)hexcl_qg->Clone();
751	<	plotTools->Area(hPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
752	<	TGraph * gPLexpexcl_qg = plotTools->GetContour(hPLexpexcl_qg,3,0);
753	<	if (gPLexpexcl_qg)   {gPLexpexcl_qg->SetLineStyle(2);gPLexpexcl_qg->SetLineColor(2);}
754	<	TGraph * gPLobsexcl_qg = plotTools->GetContour(hPLobsexcl_qg,3,0);
755	<	if (gPLobsexcl_qg) gPLobsexcl_qg->SetLineColor(2);
756	<	TH2F*hCLsExpHT_qg=(TH2F*)hexcl_qg->Clone();
757	<	TH2F*hCLsObsHT_qg=(TH2F*)hexcl_qg->Clone();
758	<	plotToolsHT->Area(hCLsExpHT_qg, MSquarkL, MGluino, ExpExclCL);
759	<	plotToolsHT->Area(hCLsObsHT_qg, MSquarkL, MGluino, ObsExclCL);
760	<	TGraph * gCLsExpExclHT_qg = plotTools->GetContour(hCLsExpHT_qg,3,0);
761	<	TGraph * gCLsObsExclHT_qg = plotTools->GetContour(hCLsObsHT_qg,3,0);
762	<	if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->SetLineColor(2);
763	<	if (gCLsExpExclHT_qg) {gCLsExpExclHT_qg->SetLineStyle(2);gCLsExpExclHT_qg->SetLineColor(2);}
746	>	50,200,800,30,200,900);
747	>	TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 2,2);
748	>	TGraph * gexpexcl_qg_LO      = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 2,4);
749	>	TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 2,1);
750	>	TGraph * gMCMCexpexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
751	>	TGraph * gMCMCobsexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
752	>	TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 1,2);
753	>	TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,4);
754	>	TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 1,1);
755	>	TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
756	>	TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
757	>	Smooth(gexpexcl_qg, 11);
758	>	Smooth(gobsexcl_qg, 11);
759	>	Smooth(gCLsExpExclHT_qg, 11);
760	>	Smooth(gCLsObsExclHT_qg, 11);
761	>	Smooth(gCLsExpExclHTm1_qg, 11);
762	>	Smooth(gCLsExpExclHTp1_qg, 11);
763	>	TGraph * gCLsExp1Sigma_qg    = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
764		hexcl_qg->Draw("");
765	<	if (gobsexcl_qg) gobsexcl_qg->Draw("l");
765	>	//
766	>	if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
767	>	if (gobsexcl_qg)      gobsexcl_qg->Draw("l");
768		if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
769	<	if (gexpexcl_qg) gexpexcl_qg->Draw("l");
490	<	//if (gPLobsexcl_qg) gPLobsexcl_qg->Draw("l");
491	<	//if (gPLexpexcl_qg) gPLexpexcl_qg->Draw("l");
769	>	if (gexpexcl_qg)      gexpexcl_qg->Draw("l");
770		if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
771	+	if (gexpexcl_qg_LO)      gexpexcl_qg_LO->Draw("l");
772	+	if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
773	+	//
774		TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
775		leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
776	<	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"Observed (MHT, CLs, TLimit)","l");
777	<	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"Expected (MHT, CLs, TLimit)","l");
778	<	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"Observed (HT, CLs, TLimit)","l");
779	<	if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"Expected (HT, CLs, TLimit)","l");
776	>	TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
777	>	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
778	>	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
779	>	if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
780	>	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
781	>	if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf");
782	>	if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l");
783		//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
784		//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
785		leg_qg->Draw();
786		c1->SaveAs("results_tb3/Exclusion_mSql_mGl_tb3.pdf");
787	<	*/
788	<	//c1->SaveAs("plot_tb3.pdf");
787	>	//c1->SaveAs("plot_tb10.pdf");
788	>
789	>
790
791		c1->SetLogy(1);
792		c1->SetLogx(1);

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines