ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/auterman/SusyScan/PlotScript/plot50tb.cc (file contents):
Revision 1.3 by auterman, Wed Feb 9 21:00:29 2011 UTC vs.
Revision 1.10 by auterman, Wed Jun 22 16:02:08 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 37 | Line 38 | int plot(int argc, char** argv)
38		return 1;
39		}
40
41	+	/*
42		//Style stuff
43		gStyle->SetHistFillColor(0);
44		gStyle->SetPalette(1);
#	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::paper();
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_tb50/filelist.txt");
88	>	//genpoints->Fill("limits_MHT_tb50_20110303/filelist.txt");
89	>	genpoints->Fill("limits_moriond_MHT.old/filelist_tb50.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92	<	genpointsHT->Fill("limits_HT_tb50/filelist.txt");
92	>	//genpointsHT->Fill("limits_HT_tb50_20110303/filelist.txt");
93	>	genpointsHT->Fill("limits_moriond_HT.old/filelist_tb50.txt");
94
95
96		//Replace read limits with specific numbers
97		//genpoints->OverwriteLimits("ABCD_MHT");
98
99	<	genpoints->FillGeneratorMasses("tb50_mu1_a0_massscan.dat");
99	>	genpoints->FillGeneratorMasses("GenScan_tb50.dat");
100		genpoints->match();
101	<	genpointsHT->FillGeneratorMasses("tb50_mu1_a0_massscan.dat");
101	>	genpointsHT->FillGeneratorMasses("GenScan_tb50.dat");
102		genpointsHT->match();
103
104	+	//genpoints->ExpandGrid(1);
105	+	//genpoints->ExpandGrid(1);
106	+	//genpoints->ExpandGrid(1);
107	+	//genpointsHT->ExpandGrid(1);
108	+	//genpointsHT->ExpandGrid(1);
109	+	//genpointsHT->ExpandGrid(1);
110	+
111		//the plotting ----------------------------------------------------------------------
112		//plotting helper functions
113		PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
#	Line 166 | Line 180 | int plot(int argc, char** argv)
180		hexplimitnsig->Draw("colz");
181		c1->SaveAs("results_tb50/ExpLimitOnNSig_m0_m12_tb50.pdf");
182
183	+	// Charged LSP in M0 - M1/2
184	+	c1->SetRightMargin ( 0.15 );
185	+	c1->SetLogz(0);
186	+	TH2F*hchargedLSP = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; Charged LSP",
187	+	40,200,600,26,140,400);
188	+	plotTools->Area(hchargedLSP, Mzero, Mhalf, ChargedLSP);
189	+	//hchargedLSP->SetMinimum(-0.01);
190	+	//hchargedLSP->SetMaximum(1.1);
191	+	hchargedLSP->Draw("colz");
192	+	c1->SaveAs("results_tb50/ChargedLSP_m0_m12_tb50.pdf");
193	+
194	+	TH2F*hnLSP = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; Neutral LSP",
195	+	40,200,600,26,140,400);
196	+	plotTools->Area(hnLSP, Mzero, Mhalf, MChi1);
197	+	//hchargedLSP->SetMinimum(-0.01);
198	+	//hchargedLSP->SetMaximum(1.1);
199	+	hnLSP->Draw("colz");
200	+	c1->SaveAs("results_tb50/NeutralLSP_m0_m12_tb50.pdf");
201	+
202		// Obs. Limit on Number of Signal Events in M0 - M1/2
203		TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
204		40,200,600,26,140,400);
#	Line 197 | Line 230 | int plot(int argc, char** argv)
230		hPLobsexcl->Draw("colz");
231		c1->SaveAs("results_tb50/PL_ObsExclusion_m0_m12_tb50.pdf");
232
233	+	c1->SetRightMargin ( 0.1 );
234		// TestContours in M0 - M1/2
235		TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
236		40,200,600,26,140,400);
237		TH2F*texpexcl=(TH2F*)texcl->Clone();
238	<	plotTools->Area(texpexcl, Mzero, Mhalf, MCMCExpExclusion);
205	<	TH2F*tobsexcl=(TH2F*)texcl->Clone();
206	<	plotTools->Area(tobsexcl, Mzero, Mhalf, MCMCObsExclusion);
238	>	plotToolsHT->Area(texpexcl, Mzero, Mhalf, NLOExpExclCLp1sigma);
239		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");
240		texpexcl->Draw("colz");
241		int col=kBlue-10;
242		for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
#	Line 228 | Line 255 | int plot(int argc, char** argv)
255
256
257		// Exclusion in M0 - M1/2
258	<	TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
259	<	40,200,600,26,140,400);
260	<	TGraph * gexpexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 2,2);
261	<	TGraph * gobsexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 2,1);
262	<	TGraph * gMCMCexpexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2);
263	<	TGraph * gMCMCobsexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1);
264	<	TGraph * gCLsExpExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 1,2);
265	<	TGraph * gCLsObsExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 1,1);
266	<	TGraph * gCLsExpExclHTm1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLm1sigma,3,0, 5,2);
267	<	TGraph * gCLsExpExclHTp1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLp1sigma,3,0, 5,2);
268	<	TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1);
269	<	hexcl->Draw("colz");
270	<	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
271	<	if (gexpexcl)         gexpexcl->Draw("l");
272	<	if (gCLsExpExclHT)    gCLsExpExclHT->Draw("l");
273	<	if (gMCMCexpexcl)     gMCMCexpexcl->Draw("l");
274	<	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
275	<	//if (gobsexcl)       gobsexcl->Draw("l");
276	<	//if (gPLobsexcl)     gPLobsexcl->Draw("l");
277	<	//if (gFCobsexcl)     gFCobsexcl->Draw("l");
278	<	//if (gMCMCobsexcl)   gMCMCobsexcl->Draw("l");
279	<	//if (gFCexpexcl)     gFCexpexcl->Draw("l");
280	<	//if (gPLexpexcl)     gPLexpexcl->Draw("l");
281	<	//if (gFCexpexcl)     gFCexpexcl->Draw("l");
258	>	c1->SetRightMargin ( 0.1 );
259	>	c1->SetTopMargin ( 0.1 );
260	>	TH2F*hexcl = new TH2F("hexcl",";m_{0} (GeV); m_{1/2} (GeV); 95% CL Expected Exclusion",
261	>	40,200,600,26,140,450);
262	>	TH2F*hs = new TH2F("hs","",40,200,600,26,140,400);
263	>	TGraph * gexpexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,0, 2,2);
264	>	TGraph * gexpexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 2,4);
265	>	TGraph * gobsexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,0, 2,1);
266	>	TGraph * gobsexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
267	>	TGraph * gMCMCexpexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion,   3,0, 3,2);
268	>	TGraph * gMCMCobsexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion,   3,0, 3,1);
269	>	TGraph * gFCexpexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,  3,0, 3,2);
270	>	TGraph * gFCexpexclHT_LO  = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion,     3,0, 3,4);
271	>	TGraph * gFCexpexclMHT    = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,    3,0, 3,2);
272	>	TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion,       3,0, 3,4);
273	>	TGraph * gFCobsexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion,  3,0, 3,1);
274	>	TGraph * gFCobsexclMHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion,  3,0, 3,1);
275	>	TGraph * gCLsExpExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,0, 1,2);
276	>	TGraph * gCLsObsExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,0, 1,1);
277	>	TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 1,4);
278	>	TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
279	>	TGraph * gCLsExpExclHTm1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
280	>	TGraph * gCLsExpExclHTp1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
281	>	TGraph * gCLsExpExclMHTm1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,  3,0, 5,2);
282	>	TGraph * gCLsExpExclMHTp1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,  3,0, 5,2);
283	>	TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
284	>	TGraph * gCheck2 = (TGraph*)gexpexcl->Clone();
285	>	Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
286	>	Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
287	>	Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
288	>	Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
289	>	Smooth( gexpexcl, 27 );
290	>	Smooth( gobsexcl_LO, 27 );gobsexcl_LO->SetLineWidth( 2 );
291	>	//Smooth( gobsexcl, 27 );
292	>	Smooth( gobsexcl, 27 );
293	>	Smooth( gexpexcl_LO, 27 );
294	>	Smooth( gFCexpexclHT, 27 );
295	>	Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
296	>	Smooth( gFCexpexclMHT, 27 );
297	>	Smooth( gFCexpexclMHT_LO, 27 );
298	>	Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
299	>	Smooth( gFCobsexclMHT, 27 );
300	>	Smooth( gCLsExpExclHTm1, 27 );
301	>	Smooth( gCLsExpExclHTp1, 27 );
302	>	Smooth( gCLsExpExclMHTm1, 27 );
303	>	Smooth( gCLsExpExclMHTp1, 27 );
304	>	TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
305	>	TGraph * CLsObsLO  = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
306	>	TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
307	>	TGraph * CLsExpLO  = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
308	>	TGraph * FCExpLO   = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
309	>	TGraph * FCExpNLO  = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
310	>	TGraph * FCObsNLO  = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
311	>	TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
312	>	TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
313	>	TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclp1, gCLsExpExclm1);
314	>	hexcl->GetYaxis()->SetTitleOffset(1.3);
315	>	hexcl->GetXaxis()->SetTitleOffset(0.92);
316	>	hexcl->Draw("colz");
317	>	//gCLsExpExclHTp1->Draw("l");
318	>	//gCLsExpExclMHTp1->Draw("l");
319	>	//gCLsExpExclHTm1->Draw("l");
320	>	//gCLsExpExclMHTm1->Draw("l");
321	>	gCLsExp1Sigma->Draw("f");
322		//set old exclusion Limits
323	+	TGraph* Atlas = Atlas_m0_m12_tb3_obs();
324	+	Atlas->SetLineColor( 28 );
325		TGraph* LEP_ch = set_lep_ch(50);
326		TGraph* LEP_sl = set_lep_sl(50);//slepton curve
327		TGraph* TEV_sg_cdf = set_tev_sg_cdf(50);//squark gluino cdf
#	Line 271 | Line 340 | int plot(int argc, char** argv)
340		First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
341		TSpline3 *sFirst = new TSpline3("sFirst",First);
342		sFirst->SetLineColor(kBlue);
343	<	sFirst->SetLineWidth(3);
343	>	//sFirst->SetLineWidth(3);
344		TSpline3 *sSecond = new TSpline3("sSecond",Second);
345		sSecond->SetLineColor(kBlue);
346		sSecond->SetLineStyle(2);
347	<	sSecond->SetLineWidth(3);
347	>	//sSecond->SetLineWidth(3);
348		TSpline3 *sThird = new TSpline3("sThird",Third);
349	<	sThird->SetLineColor(kGreen+2);
349	>	sThird->SetLineColor(kBlue);
350		sThird->SetLineStyle(4);
351	<	sThird->SetLineWidth(3);
352	<	TEV_sn_d0_1->Draw("fsame");
353	<	TEV_sn_d0_2->Draw("fsame");
354	<	//TEV_sg_d0->Draw("fsame");
351	>	//sThird->SetLineWidth(3);
352	>	//TEV_sn_d0_1->Draw("fsame");
353	>	//TEV_sn_d0_2->Draw("fsame"); //only for tb=3
354	>	LEP_ch->Draw("fsame");
355	>	LEP_sl->Draw("fsame");
356		//TEV_sg_cdf->Draw("fsame");
357	+	//TEV_sg_cdf->Draw("lsame");
358	+	//TEV_sg_d0->Draw("fsame");
359	+	sFirst->Draw("same");
360	+	//sSecond->Draw("same");
361	+	//sThird->Draw("same");
362	+	//Atlas->Draw("c,same");
363	+	TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
364	+	TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
365	+	ms.DrawLatex(490,458,"tan #beta=50, #mu>0, A_{0}=0");
366	+
367	+	TLegend* legexp = new TLegend(0.7,0.764,0.93,0.88,NULL,"brNDC");
368	+	legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
369	+	//TEV_sg_cdf.SetLineColor(1);
370	+	//legexp->SetHeader("CMS   tan #beta=50, #mu>0, A_{0}=0");
371	+	legexp->SetHeader("CMS");
372	+	//legexp->AddEntry(TEV_sg_cdf,"CDF  #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan #beta=5, #mu<0}","f");
373	+	//legexp->AddEntry(TEV_sg_d0,"D0   #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan #beta=3, #mu<0}","f");
374	+	//ch_gr.SetLineColor(1);
375	+	legexp->AddEntry(LEP_ch,"LEP2   #tilde{#chi}_{1}^{#pm}","f");
376	+	//sl_gr.SetLineColor(1);
377	+	//legexp->AddEntry(LEP_sl,"LEP2   #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
378	+	//if(tanbeta == 3)
379	+	//legexp->AddEntry(TEV_sn_d0_1,"D0  #chi^{#pm}_{1}, #chi^{0}_{2}","f");
380	+	legexp->AddEntry(sFirst, "CMS #alpha_{T}");
381	+	//legexp->AddEntry(Atlas,  "Atlas, #scale[0.8]{tan #beta=3}","l");
382	+	legexp->Draw();
383	+
384	+	CLsObsNLO->Draw("l");
385	+	CLsObsLO->Draw("l");
386	+	CLsExpNLO->Draw("l");
387	+	//FCExpNLO->Draw("l");
388	+	stau->Draw("fsame");
389	+	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
390	+
391	+	//constant ssqquark and gluino lines
392	+	TF1* lnsq[4];
393	+	TF1* lngl[4];
394	+	TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
395	+	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
396	+	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
397	+	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
398	+	TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
399	+	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
400	+	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
401	+	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
402	+	int tanBeta_=50;
403	+	for(int i = 0; i < 4; i++){
404	+	lnsq[i] = constant_squark(tanBeta_,i);
405	+	lngl[i] = constant_gluino(tanBeta_,i);
406	+	lngl[i]->Draw("same");
407	+	lnsq[i]->Draw("same");
408	+	}
409	+	TLegend * leg = new TLegend(0.38,0.72,0.81,0.88);
410	+	//TLegend * leg = new TLegend(0.25,0.72,0.68,0.88);
411	+	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);leg->SetTextSize(0.025);
412	+	TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
413	+	leg->SetHeader("L_{int} = 36 pb^{-1}, #sqrt{s} = 7 TeV");
414	+	leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
415	+	leg->AddEntry(CLsObsLO, "Observed, LO","l");
416	+	leg->AddEntry(expLeg,   "Expected #pm 1#sigma, NLO","lf");
417	+	//leg->AddEntry(FCExpLO,  "Expected, FC, NLO","l");
418	+	leg->Draw();
419	+	gPad->RedrawAxis();
420	+	c1->SaveAs("results/Exclusion_m0_m12_tb50.pdf");
421	+
422	+
423	+	hexcl->Draw("");
424	+	gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb50");
425	+	CLsObsNLO->SetName("CLsNLO_ObservedTb50");
426	+	CLsExpNLO->SetName("CLsNLO_ExpectedTb50");
427	+	CLsExpLO->SetName("CLsLO_ExpectedTb50");
428	+	gCLsExp1Sigma->Draw("l");
429	+	CLsObsNLO->Draw("l");
430	+	CLsExpNLO->Draw("l");
431	+	CLsObsLO->Draw("l");
432	+	c1->SaveAs("results/LimitContours_tb50.C");
433	+
434	+
435	+
436	+	hexcl->Draw("colz");
437		LEP_ch->Draw("fsame");
438		LEP_sl->Draw("fsame");
439	+	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
440	+	sFirst->Draw("same");
441	+	Atlas->Draw("c,same");
442	+	ms.DrawLatex(490,458,"tan #beta=50, #mu>0, A_{0}=0");
443	+	CLsObsNLO->Draw("l");
444	+	//CLsObsLO->Draw("l");
445	+	CLsExpNLO->Draw("l");
446	+	//FCExpNLO->Draw("l");
447	+	//Jim's limits///////////////////////////////////////////////////////////////////////////////
448	+	TGraph* JimObsHT  = Jim_ht_tb50(0);
449	+	TGraph* JimObsMHT = Jim_mht_tb50(0);
450	+	TGraph* JimExpHT  = Jim_ht_tb50(1);
451	+	TGraph* JimExpMHT = Jim_mht_tb50(1);
452	+	TGraph* JimExpHTup  = Jim_ht_tb50(2);
453	+	TGraph* JimExpMHTup = Jim_mht_tb50(2);
454	+	TGraph* JimExpHTdn  = Jim_ht_tb50(3);
455	+	TGraph* JimExpMHTdn = Jim_mht_tb50(3);
456	+	Smooth(JimObsHT, 2);
457	+	Smooth(JimObsMHT,2);
458	+	Smooth(JimExpHT, 2);
459	+	Smooth(JimExpMHT,2);
460	+	Smooth(JimExpHTup, 2);
461	+	Smooth(JimExpMHTup,2);
462	+	Smooth(JimExpHTdn, 2);
463	+	Smooth(JimExpMHTdn,2);
464	+	TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
465	+	TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
466	+	TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
467	+	TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
468	+	TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
469	+	TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
470	+	JimExp1Sigma->Draw("f");
471	+	JimObs->Draw("c");
472	+	JimExp->Draw("c");
473		stau->Draw("fsame");
474	+	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
475	+	//constant ssqquark and gluino lines
476	+	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
477	+	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
478	+	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
479	+	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
480	+	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
481	+	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
482	+	for(int i = 0; i < 4; i++){
483	+	lngl[i]->Draw("same");
484	+	lnsq[i]->Draw("same");
485	+	}
486	+	legexp->Draw();
487	+	TLegend * legBayes = new TLegend(0.38,0.8,0.73,0.92);
488	+	legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
489	+	legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
490	+	legBayes->AddEntry(JimObs,   "Observed, Bayes","l");
491	+	legBayes->AddEntry(JimLeg,   "Expected #pm 1#sigma, Bayes","lf");
492	+	legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
493	+	legBayes->AddEntry(expLeg,   "Expected #pm 1#sigma, CLs","lf");
494	+	legBayes->Draw();
495	+	gPad->RedrawAxis();
496	+	c1->SaveAs("results/Exclusion_m0_m12_tb50_Bayes.pdf");
497	+
498	+
499	+	//The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
500	+	TGraph * gCLsExpNoSExclMHT   = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
501	+	TGraph * gCLsExpNoSExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
502	+	TGraph * gCLsExpNoSExclHTm1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
503	+	TGraph * gCLsExpNoSExclHTp1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
504	+	TGraph * gCLsExpNoSExclMHTm1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
505	+	TGraph * gCLsExpNoSExclMHTp1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
506	+	Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
507	+	Smooth( gCLsExpNoSExclHT, 27 );  gCLsExpNoSExclHT->SetLineWidth( 3 );
508	+	Smooth( gCLsExpNoSExclHTm1, 27 );
509	+	Smooth( gCLsExpNoSExclHTp1, 27 );
510	+	Smooth( gCLsExpNoSExclMHTm1, 27 );
511	+	Smooth( gCLsExpNoSExclMHTp1, 27 );
512	+	TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
513	+	TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
514	+	TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
515	+	TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
516	+	hexcl->Draw("colz");
517	+	LEP_ch->Draw("fsame");
518	+	LEP_sl->Draw("fsame");
519	+	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
520		sFirst->Draw("same");
521	<	sSecond->Draw("same");
522	<	sThird->Draw("same");
523	<	TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
524	<	//b.DrawLatex( 10,420,"#tilde{#tau} LSP");
525	<	//b.DrawLatex(  6,150,"LEP2");
526	<	//b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
527	<	//b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
528	<	//b.DrawLatex(100,140,"CDF");
529	<	b.DrawLatex(230,145,"D0 #tilde{g}, #tilde{q}");
530	<	//b.DrawLatex( 80,180,"D0 #tilde{#nu}");
531	<	TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
532	<	leg->SetBorderSize(0);leg->SetFillColor(0);
533	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
534	<	if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs)","l");
535	<	if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs)","l");
536	<	if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs)","l");
537	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
538	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
539	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
540	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
541	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
542	<	if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
543	<	if (sFirst)  leg->AddEntry(sFirst, "NLO RA1 Observed");
544	<	if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected");
545	<	if (sThird) leg->AddEntry(sThird,"LO RA1 Observed");
521	>	//Atlas->Draw("c,same");
522	>	ms.DrawLatex(490,380,"tan #beta=50, #mu>0, A_{0}=0");
523	>	gCLsExpNoS1Sigma->Draw("lf,same");
524	>	sFirst->Draw("same");
525	>	//Atlas->Draw("c,same");
526	>	CLsObsNLO->Draw("l,same");
527	>	CLsObsLO->Draw("l,same");
528	>	CLsExpNoSNLO->Draw("l,same");
529	>	//gCLsExpNoSExclMHT->Draw("l,same");
530	>	//gCLsExpNoSExclHT->Draw("l,same");
531	>	//gCLsObsExclHT->Draw("l,same");
532	>	//gobsexcl->Draw("l,same");
533	>	//FCExpNLO->Draw("l");
534	>	stau->Draw("fsame");
535	>	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
536	>	//constant ssqquark and gluino lines
537	>	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
538	>	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
539	>	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
540	>	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
541	>	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
542	>	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
543	>	for(int i = 0; i < 4; i++){
544	>	lngl[i]->Draw("same");
545	>	lnsq[i]->Draw("same");
546	>	}
547	>	legexp->Draw();
548		leg->Draw();
549		gPad->RedrawAxis();
550	<	c1->SaveAs("results_tb50/Exclusion_m0_m12_tb50.pdf");
550	>	c1->SaveAs("results/Exclusion_m0_m12_tb50_NoSigHypPseudoData.pdf");
551	>
552	>
553
554		// Signal Contamination in M0 - M1/2
555		c1->SetLogz(0);
#	Line 339 | Line 573 | int plot(int argc, char** argv)
573		if (gobsexcl)    gobsexcl->Draw("l");
574		c1->SaveAs("results_tb50/SignalContamination_MHT_m0_m12_tb50.pdf");
575
576	+	// Signal JEC Uncertainty  MHT
577	+	c1->SetLogz(0);
578	+	TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
579	+	40,200,600,26,140,400);
580	+	plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
581	+	hsig_jec_mht->SetMinimum(0.0);
582	+	hsig_jec_mht->SetMaximum(0.25);
583	+	hsig_jec_mht->SetContour(10);
584	+	hsig_jec_mht->Draw("colz");
585	+	if (gobsexcl)    gobsexcl->Draw("l");
586	+	c1->SaveAs("results_tb50/SigJEC_MHT_m0_m12_tb50.pdf");
587	+
588	+	// Signal JEC Uncertainty  HT
589	+	c1->SetLogz(0);
590	+	TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
591	+	40,200,600,26,140,400);
592	+	plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
593	+	hsig_jec_ht->SetMinimum(0.0);
594	+	hsig_jec_ht->SetMaximum(0.25);
595	+	hsig_jec_ht->SetContour(10);
596	+	hsig_jec_ht->Draw("colz");
597	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
598	+	c1->SaveAs("results_tb50/SigJEC_HT_m0_m12_tb50.pdf");
599	+
600	+
601	+	// Signal MuIso Uncertainty  MHT
602	+	c1->SetLogz(0);
603	+	TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
604	+	40,200,600,26,140,400);
605	+	plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
606	+	hsig_MuIso_mht->SetMinimum(0.0);
607	+	hsig_MuIso_mht->SetMaximum(0.05);
608	+	hsig_MuIso_mht->SetContour(5);
609	+	hsig_MuIso_mht->Draw("colz");
610	+	if (gobsexcl)    gobsexcl->Draw("l");
611	+	c1->SaveAs("results_tb50/SigMuIso_MHT_m0_m12_tb50.pdf");
612	+
613	+	// Signal MuIso Uncertainty  HT
614	+	c1->SetLogz(0);
615	+	TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
616	+	40,200,600,26,140,400);
617	+	plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
618	+	hsig_MuIso_ht->SetMinimum(0.0);
619	+	hsig_MuIso_ht->SetMaximum(0.05);
620	+	hsig_MuIso_ht->SetContour(5);
621	+	hsig_MuIso_ht->Draw("colz");
622	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
623	+	c1->SaveAs("results_tb50/SigMuIso_HT_m0_m12_tb50.pdf");
624	+
625	+
626	+	// Signal kFactor Uncertainty  MHT
627	+	c1->SetLogz(0);
628	+	TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
629	+	40,200,600,26,140,400);
630	+	plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
631	+	hsig_kFactor_mht->SetMinimum(0.1);
632	+	hsig_kFactor_mht->SetMaximum(0.2);
633	+	hsig_kFactor_mht->SetContour(10);
634	+	hsig_kFactor_mht->Draw("colz");
635	+	if (gobsexcl)    gobsexcl->Draw("l");
636	+	c1->SaveAs("results_tb50/SigkFactor_MHT_m0_m12_tb50.pdf");
637	+
638	+	// Signal kFactor Uncertainty  HT
639	+	c1->SetLogz(0);
640	+	TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
641	+	40,200,600,26,140,400);
642	+	plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
643	+	hsig_kFactor_ht->SetMinimum(0.1);
644	+	hsig_kFactor_ht->SetMaximum(0.2);
645	+	hsig_kFactor_ht->SetContour(10);
646	+	hsig_kFactor_ht->Draw("colz");
647	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
648	+	c1->SaveAs("results_tb50/SigkFactor_HT_m0_m12_tb50.pdf");
649	+
650
343	–	/*
651		//-----------------------------------------------------------------------------------
652		c1->SetLogz(1);
653
#	Line 349 | Line 656 | int plot(int argc, char** argv)
656
657		// cross-section in squark - gluino mass
658		TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
659	<	60,200,1400,50,200,1200);
659	>	40,200,900,30,200,900);
660		plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
661		hxsec_qg->SetMinimum(0.01);
662		hxsec_qg->Draw("colz");
#	Line 359 | Line 666 | int plot(int argc, char** argv)
666
667		// Observed Limit in squark - gluino mass
668		TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
669	<	60,200,1400,50,200,1200);
669	>	40,200,900,30,200,900);
670		plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
671		hobslimit_qg->SetMinimum(0.01);
672		hobslimit_qg->Draw("colz");
#	Line 367 | Line 674 | int plot(int argc, char** argv)
674
675		// Expected Limit in squark - gluino mass
676		TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
677	<	60,200,1400,50,200,1200);
677	>	40,200,900,30,200,900);
678		plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
679		hexplimit_qg->SetMinimum(0.01);
680		hexplimit_qg->Draw("colz");
#	Line 375 | Line 682 | int plot(int argc, char** argv)
682
683		// Signal Acceptance in squark - gluino mass
684		TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
685	<	60,200,1400,50,200,1200);
685	>	40,200,900,30,200,900);
686		plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
687		hsigacc_qg->SetMinimum(0.01);
688		hsigacc_qg->SetMaximum(1.0);
#	Line 385 | Line 692 | int plot(int argc, char** argv)
692		// Exp. Limit on Number of Signal Events in squark - gluino mass
693		c1->SetLogz(0);
694		TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
695	<	60,200,1400,50,200,1200);
695	>	40,200,900,30,200,900);
696		plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
697		hexplimitnsig_qg->SetMinimum(0.0);
698		hexplimitnsig_qg->SetMaximum(20);
#	Line 394 | Line 701 | int plot(int argc, char** argv)
701
702		// Obs. Limit on Number of Signal Events in squark - gluino mass
703		TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
704	<	60,200,1400,50,200,1200);
704	>	40,200,900,30,200,900);
705		plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
706		hobslimitnsig_qg->SetMinimum(0.0);
707		hobslimitnsig_qg->SetMaximum(20);
#	Line 404 | Line 711 | int plot(int argc, char** argv)
711		c1->SetLogz(0);
712		// Expected Exclusion in squark - gluino mass
713		TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
714	<	60,200,1400,50,200,1200);
714	>	40,200,900,30,200,900);
715		plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
716		hexpexcl_qg->Draw("colz");
717		c1->SaveAs("results_tb50/ExpExclusion_mSql_mGl_tb50.pdf");
718
719		// Observed Exclusion in squark - gluino mass
720		TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
721	<	60,200,1400,50,200,1200);
721	>	40,200,900,30,200,900);
722		plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
723		//std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
724		hobsexcl_qg->Draw("colz");
#	Line 423 | Line 730 | int plot(int argc, char** argv)
730
731		// TestContours in M0 - M1/2
732		TH2F*texcl_qg = new TH2F("texcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
733	<	60,200,1400,50,200,1200);
733	>	40,200,900,30,200,900);
734		TH2F*tPLobsexcl_qg=(TH2F*)texcl_qg->Clone();
735		plotTools->Area(tPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
736		TH2F*tPLexpexcl_qg=(TH2F*)texcl_qg->Clone();
#	Line 454 | Line 761 | int plot(int argc, char** argv)
761
762		// Observed Exclusion in squark - gluino mass
763		TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
764	<	60,200,1400,50,200,1200);
764	>	40,200,900,30,200,900);
765		plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
766		//hPLobsexcl_qg->Draw("colz");
767		//c1->SaveAs("results_tb50/PL_ObsExclusion_mSql_mGl_tb50.pdf");
768
769
770
464	–
465	–
771		// Exclusion in squark - gluino mass ----------------------------------------
772		TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
773	<	60,200,1400,50,200,1200);
774	<	TGraph * gexpexcl_qg = plotTools->GetContour(hexpexcl_qg,3,0);
775	<	if (gexpexcl_qg) gexpexcl_qg->SetLineStyle(2);
776	<	TGraph * gobsexcl_qg = plotTools->GetContour(hobsexcl_qg,3,0);
777	<	TH2F*hPLexpexcl_qg=(TH2F*)hexcl_qg->Clone();
778	<	plotTools->Area(hPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
779	<	TGraph * gPLexpexcl_qg = plotTools->GetContour(hPLexpexcl_qg,3,0);
780	<	if (gPLexpexcl_qg)   {gPLexpexcl_qg->SetLineStyle(2);gPLexpexcl_qg->SetLineColor(2);}
781	<	TGraph * gPLobsexcl_qg = plotTools->GetContour(hPLobsexcl_qg,3,0);
782	<	if (gPLobsexcl_qg) gPLobsexcl_qg->SetLineColor(2);
783	<	TH2F*hCLsExpHT_qg=(TH2F*)hexcl_qg->Clone();
784	<	TH2F*hCLsObsHT_qg=(TH2F*)hexcl_qg->Clone();
480	<	plotToolsHT->Area(hCLsExpHT_qg, MSquarkL, MGluino, ExpExclCL);
481	<	plotToolsHT->Area(hCLsObsHT_qg, MSquarkL, MGluino, ObsExclCL);
482	<	TGraph * gCLsExpExclHT_qg = plotTools->GetContour(hCLsExpHT_qg,3,0);
483	<	TGraph * gCLsObsExclHT_qg = plotTools->GetContour(hCLsObsHT_qg,3,0);
484	<	if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->SetLineColor(2);
485	<	if (gCLsExpExclHT_qg) {gCLsExpExclHT_qg->SetLineStyle(2);gCLsExpExclHT_qg->SetLineColor(2);}
773	>	40,200,900,30,200,900);
774	>	TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  NLOExpExclCL,       3,0, 2,2);
775	>	TGraph * gexpexcl_qg_LO      = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  ExpExclCL,       3,0, 2,4);
776	>	TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  NLOObsExclCL,       3,0, 2,1);
777	>	TGraph * gMCMCexpexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
778	>	TGraph * gMCMCobsexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
779	>	TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 1,2);
780	>	TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,4);
781	>	TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 1,1);
782	>	TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
783	>	TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
784	>	TGraph * gCLsExp1Sigma_qg    = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
785		hexcl_qg->Draw("");
786	<	if (gobsexcl_qg) gobsexcl_qg->Draw("l");
786	>	//
787	>	if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
788	>	if (gobsexcl_qg)      gobsexcl_qg->Draw("l");
789		if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
790	<	if (gexpexcl_qg) gexpexcl_qg->Draw("l");
791	<	//if (gPLobsexcl_qg) gPLobsexcl_qg->Draw("l");
491	<	//if (gPLexpexcl_qg) gPLexpexcl_qg->Draw("l");
790	>	if (gexpexcl_qg)      gexpexcl_qg->Draw("l");
791	>	if (gexpexcl_qg_LO)   gexpexcl_qg_LO->Draw("l");
792		if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
793	+	if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
794	+	//
795		TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
796	<	leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
797	<	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"Observed (MHT, CLs, TLimit)","l");
798	<	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"Expected (MHT, CLs, TLimit)","l");
799	<	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"Observed (HT, CLs, TLimit)","l");
800	<	if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"Expected (HT, CLs, TLimit)","l");
796	>	leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);leg->SetTextFont(42);
797	>	TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
798	>	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
799	>	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
800	>	if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
801	>	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
802	>	if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf");
803	>	if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l");
804		//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
805		//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
806		leg_qg->Draw();
807		c1->SaveAs("results_tb50/Exclusion_mSql_mGl_tb50.pdf");
808	<
808	>	//c1->SaveAs("plot_tb10.pdf");
809	>
810	>
811	>	/*
812		//c1->SaveAs("plot_tb50.pdf");
505	–	*/
813		c1->SetLogy(1);
814		c1->SetLogx(1);
815		TH2F*hUncertScan = new TH2F("hUncertScan",";relative signal uncertainty; number of signal events",
#	Line 560 | Line 867 | int plot(int argc, char** argv)
867		leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
868		leg_S->Draw();
869		c1->SaveAs("results_tb50/UncertaintyScan.pdf");
870	<
870	>	*/
871		//theApp.Run();
872		}
873

Diff Legend

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