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.6 by auterman, Wed Mar 2 13:31:33 2011 UTC vs.
Revision 1.8 by auterman, Fri May 20 07:39:34 2011 UTC

#	Line 85 | Line 85 | int plot(int argc, char** argv)
85		//Get limits and generator masses ---------------------------------------------------
86		TheLimits * genpoints = new TheLimits();
87		//genpoints->Fill(argc, argv);
88	<	genpoints->Fill("limits_MHT_tb50-20110301/filelist.txt");
88	>	//genpoints->Fill("limits_MHT_tb50_20110303/filelist.txt");
89	>	genpoints->Fill("limits_moriond_MHT/filelist_tb50.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92	<	genpointsHT->Fill("limits_HT_tb50-20110301/filelist.txt");
92	>	//genpointsHT->Fill("limits_HT_tb50_20110303/filelist.txt");
93	>	genpointsHT->Fill("limits_moriond_HT/filelist_tb50.txt");
94
95
96		//Replace read limits with specific numbers
#	Line 99 | Line 101 | int plot(int argc, char** argv)
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 171 | 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 202 | 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);
#	Line 232 | Line 261 | int plot(int argc, char** argv)
261		TGraph * gexpexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,0, 2,2);
262		TGraph * gexpexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 2,4);
263		TGraph * gobsexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,0, 2,1);
264	+	TGraph * gobsexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
265		TGraph * gMCMCexpexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion,   3,0, 3,2);
266		TGraph * gMCMCobsexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion,   3,0, 3,1);
267		TGraph * gFCexpexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,  3,0, 3,2);
#	Line 243 | Line 273 | int plot(int argc, char** argv)
273		TGraph * gCLsExpExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,0, 1,2);
274		TGraph * gCLsObsExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,0, 1,1);
275		TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 1,4);
276	<	TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 1,1);
276	>	TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
277		TGraph * gCLsExpExclHTm1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
278		TGraph * gCLsExpExclHTp1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
279		TGraph * gCLsExpExclMHTm1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,  3,0, 5,2);
280		TGraph * gCLsExpExclMHTp1 = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,  3,0, 5,2);
281		TGraph * gCheck1 = (TGraph*)gobsexcl->Clone();
282		TGraph * gCheck2 = (TGraph*)gexpexcl->Clone();
283	<	Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 );
284	<	Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
285	<	Smooth( gCLsObsExclHT_LO, 33); gCLsObsExclHT_LO->SetLineWidth( 3 );
286	<	Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
287	<	Smooth( gexpexcl, 33 );
288	<	Smooth( gobsexcl, 33 );
289	<	Smooth( gobsexcl, 33 );
290	<	Smooth( gexpexcl_LO, 33 );
291	<	Smooth( gFCexpexclHT, 33 );
292	<	Smooth( gFCexpexclHT_LO, 33 );gFCexpexclHT_LO->SetLineWidth( 3 );
293	<	Smooth( gFCexpexclMHT, 33 );
294	<	Smooth( gFCexpexclMHT_LO, 33 );
295	<	Smooth( gFCobsexclHT, 33 );gFCobsexclHT->SetLineWidth( 3 );
296	<	Smooth( gFCobsexclMHT, 33 );
297	<	Smooth( gCLsExpExclHTm1, 33 );
298	<	Smooth( gCLsExpExclHTp1, 33 );
299	<	Smooth( gCLsExpExclMHTm1, 33 );
300	<	Smooth( gCLsExpExclMHTp1, 33 );
283	>	Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
284	>	Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
285	>	Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
286	>	Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
287	>	Smooth( gexpexcl, 27 );
288	>	Smooth( gobsexcl_LO, 27 );gobsexcl_LO->SetLineWidth( 2 );
289	>	//Smooth( gobsexcl, 27 );
290	>	Smooth( gobsexcl, 27 );
291	>	Smooth( gexpexcl_LO, 27 );
292	>	Smooth( gFCexpexclHT, 27 );
293	>	Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
294	>	Smooth( gFCexpexclMHT, 27 );
295	>	Smooth( gFCexpexclMHT_LO, 27 );
296	>	Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
297	>	Smooth( gFCobsexclMHT, 27 );
298	>	Smooth( gCLsExpExclHTm1, 27 );
299	>	Smooth( gCLsExpExclHTp1, 27 );
300	>	Smooth( gCLsExpExclMHTm1, 27 );
301	>	Smooth( gCLsExpExclMHTp1, 27 );
302		TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
303	+	TGraph * CLsObsLO  = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
304		TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
305		TGraph * CLsExpLO  = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
306		TGraph * FCExpLO   = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
307	+	TGraph * FCExpNLO  = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
308		TGraph * FCObsNLO  = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
309		TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
310		TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
311		TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclp1, gCLsExpExclm1);
312		hexcl->GetYaxis()->SetTitleOffset(1.3);
313	+	hexcl->GetXaxis()->SetTitleOffset(0.92);
314		hexcl->Draw("colz");
315		//gCLsExpExclHTp1->Draw("l");
316		//gCLsExpExclMHTp1->Draw("l");
#	Line 284 | Line 318 | int plot(int argc, char** argv)
318		//gCLsExpExclMHTm1->Draw("l");
319		gCLsExp1Sigma->Draw("f");
320		//set old exclusion Limits
321	+	TGraph* Atlas = Atlas_m0_m12_tb3_obs();
322	+	Atlas->SetLineColor( 28 );
323		TGraph* LEP_ch = set_lep_ch(50);
324		TGraph* LEP_sl = set_lep_sl(50);//slepton curve
325		TGraph* TEV_sg_cdf = set_tev_sg_cdf(50);//squark gluino cdf
#	Line 311 | Line 347 | int plot(int argc, char** argv)
347		sThird->SetLineColor(kBlue);
348		sThird->SetLineStyle(4);
349		//sThird->SetLineWidth(3);
350	<	TEV_sn_d0_1->Draw("fsame");
351	<	TEV_sn_d0_2->Draw("fsame");
316	<	TEV_sg_d0->Draw("fsame");
317	<	TEV_sg_cdf->Draw("fsame");
350	>	//TEV_sn_d0_1->Draw("fsame");
351	>	//TEV_sn_d0_2->Draw("fsame"); //only for tb=3
352		LEP_ch->Draw("fsame");
353		LEP_sl->Draw("fsame");
354	<	stau->Draw("fsame");
354	>	//TEV_sg_cdf->Draw("fsame");
355	>	//TEV_sg_cdf->Draw("lsame");
356	>	//TEV_sg_d0->Draw("fsame");
357		sFirst->Draw("same");
358	<	sSecond->Draw("same");
359	<	sThird->Draw("same");
358	>	//sSecond->Draw("same");
359	>	//sThird->Draw("same");
360	>	Atlas->Draw("c,same");
361		TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
362	<	//b.DrawLatex( 10,420,"#tilde{#tau} LSP");
363	<	b.DrawLatex(230,145,"LEP2");
364	<	//b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
365	<	//b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
366	<	//b.DrawLatex(100,140,"CDF");
367	<	//b.DrawLatex(230,145,"D0 #tilde{g}, #tilde{q}");
368	<	//b.DrawLatex( 80,180,"D0 #tilde{#nu}");
362	>	TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
363	>	ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
364	>	TLegend* legexp = new TLegend(0.69,0.77,0.93,0.92,NULL,"brNDC");
365	>	legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
366	>	//TEV_sg_cdf.SetLineColor(1);
367	>	legexp->SetHeader("CMS preliminary");
368	>	//legexp->AddEntry(TEV_sg_cdf,"CDF  #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
369	>	//legexp->AddEntry(TEV_sg_d0,"D0   #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=3, #mu<0}","f");
370	>	//ch_gr.SetLineColor(1);
371	>	legexp->AddEntry(LEP_ch,"LEP2   #tilde{#chi}_{1}^{#pm}","f");
372	>	//sl_gr.SetLineColor(1);
373	>	//legexp->AddEntry(LEP_sl,"LEP2   #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
374	>	//if(tanbeta == 3)
375	>	//legexp->AddEntry(TEV_sn_d0_1,"D0  #chi^{#pm}_{1}, #chi^{0}_{2}","f");
376	>	legexp->AddEntry(sFirst, "CMS #alpha_{T}");
377	>	legexp->AddEntry(Atlas,  "Atlas, #scale[0.8]{tan#beta=3}","l");
378	>	legexp->Draw();
379	>
380		CLsObsNLO->Draw("l");
381	+	CLsObsLO->Draw("l");
382		CLsExpNLO->Draw("l");
383	<	CLsExpLO->Draw("l");
384	<	FCExpLO->Draw("l");
385	<	//FCObsNLO->Draw("l");
386	<
387	<	//gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c");
388	<	//sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c");
389	<	//TLegend * leg = new TLegend(0.45,0.65,0.84,0.89);
390	<	//leg->SetBorderSize(0);leg->SetFillColor(0);
391	<	TLegend * leg = new TLegend(0.6,0.66,0.96,0.90);
392	<	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);
383	>	//FCExpNLO->Draw("l");
384	>	stau->Draw("fsame");
385	>	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
386	>
387	>	//constant ssqquark and gluino lines
388	>	TF1* lnsq[4];
389	>	TF1* lngl[4];
390	>	TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
391	>	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
392	>	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
393	>	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
394	>	TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
395	>	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
396	>	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
397	>	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
398	>	int tanBeta_=50;
399	>	for(int i = 0; i < 4; i++){
400	>	lnsq[i] = constant_squark(tanBeta_,i);
401	>	lngl[i] = constant_gluino(tanBeta_,i);
402	>	lngl[i]->Draw("same");
403	>	lnsq[i]->Draw("same");
404	>	}
405	>	TLegend * leg = new TLegend(0.38,0.8,0.73,0.92);
406	>	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
407		TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
408	<	leg->SetHeader("CMS preliminary, 95%CL");
409	<	leg->AddEntry(expLeg,   "Expected (CLs, NLO)","l");
410	<	leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l");
411	<	leg->AddEntry(CLsExpLO, "Expected (CLs, LO)","l");
412	<	leg->AddEntry(FCExpLO,  "Expected (FC, LO)","l");
413	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
414	<	//if (gCLsExpExclHT) leg->AddEntry(expLeg       ,"NLO Expected (HT, CLs)","l");
415	<	//if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
416	<	//if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
417	<	//if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
418	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
419	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
420	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l");
421	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l");
422	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
423	<	//if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
424	<	//if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
425	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l");
426	<	if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)");
427	<	if (sFirst)  leg->AddEntry(sFirst, "#alpha_{T} Observed (FC, NLO)");
428	<	if (sThird)  leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)");
408	>	leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
409	>	leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
410	>	leg->AddEntry(CLsObsLO, "Observed, LO","l");
411	>	leg->AddEntry(expLeg,   "Expected #pm 1#sigma, NLO","lf");
412	>	//leg->AddEntry(FCExpLO,  "Expected, FC, NLO","l");
413	>	leg->Draw();
414	>	gPad->RedrawAxis();
415	>	c1->SaveAs("results/Exclusion_m0_m12_tb50.pdf");
416	>
417	>
418	>	hexcl->Draw("");
419	>	gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb50");
420	>	CLsObsNLO->SetName("CLsNLO_ObservedTb50");
421	>	CLsExpNLO->SetName("CLsNLO_ExpectedTb50");
422	>	CLsExpLO->SetName("CLsLO_ExpectedTb50");
423	>	gCLsExp1Sigma->Draw("l");
424	>	CLsObsNLO->Draw("l");
425	>	CLsExpNLO->Draw("l");
426	>	CLsObsLO->Draw("l");
427	>	c1->SaveAs("results/LimitContours_tb50.C");
428	>
429	>
430	>
431	>	hexcl->Draw("colz");
432	>	LEP_ch->Draw("fsame");
433	>	LEP_sl->Draw("fsame");
434	>	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
435	>	sFirst->Draw("same");
436	>	Atlas->Draw("c,same");
437	>	ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
438	>	CLsObsNLO->Draw("l");
439	>	//CLsObsLO->Draw("l");
440	>	CLsExpNLO->Draw("l");
441	>	//FCExpNLO->Draw("l");
442	>	//Jim's limits///////////////////////////////////////////////////////////////////////////////
443	>	TGraph* JimObsHT  = Jim_ht_tb50(0);
444	>	TGraph* JimObsMHT = Jim_mht_tb50(0);
445	>	TGraph* JimExpHT  = Jim_ht_tb50(1);
446	>	TGraph* JimExpMHT = Jim_mht_tb50(1);
447	>	TGraph* JimExpHTup  = Jim_ht_tb50(2);
448	>	TGraph* JimExpMHTup = Jim_mht_tb50(2);
449	>	TGraph* JimExpHTdn  = Jim_ht_tb50(3);
450	>	TGraph* JimExpMHTdn = Jim_mht_tb50(3);
451	>	Smooth(JimObsHT, 2);
452	>	Smooth(JimObsMHT,2);
453	>	Smooth(JimExpHT, 2);
454	>	Smooth(JimExpMHT,2);
455	>	Smooth(JimExpHTup, 2);
456	>	Smooth(JimExpMHTup,2);
457	>	Smooth(JimExpHTdn, 2);
458	>	Smooth(JimExpMHTdn,2);
459	>	TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
460	>	TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
461	>	TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
462	>	TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
463	>	TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
464	>	TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
465	>	JimExp1Sigma->Draw("f");
466	>	JimObs->Draw("c");
467	>	JimExp->Draw("c");
468	>	stau->Draw("fsame");
469	>	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
470	>	//constant ssqquark and gluino lines
471	>	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
472	>	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
473	>	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
474	>	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
475	>	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
476	>	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
477	>	for(int i = 0; i < 4; i++){
478	>	lngl[i]->Draw("same");
479	>	lnsq[i]->Draw("same");
480	>	}
481	>	legexp->Draw();
482	>	TLegend * legBayes = new TLegend(0.38,0.8,0.73,0.92);
483	>	legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
484	>	legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
485	>	legBayes->AddEntry(JimObs,   "Observed, Bayes","l");
486	>	legBayes->AddEntry(JimLeg,   "Expected #pm 1#sigma, Bayes","lf");
487	>	legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
488	>	legBayes->AddEntry(expLeg,   "Expected #pm 1#sigma, CLs","lf");
489	>	legBayes->Draw();
490	>	gPad->RedrawAxis();
491	>	c1->SaveAs("results/Exclusion_m0_m12_tb50_Bayes.pdf");
492	>
493	>
494	>	//The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
495	>	TGraph * gCLsExpNoSExclMHT   = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
496	>	TGraph * gCLsExpNoSExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
497	>	TGraph * gCLsExpNoSExclHTm1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
498	>	TGraph * gCLsExpNoSExclHTp1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
499	>	TGraph * gCLsExpNoSExclMHTm1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
500	>	TGraph * gCLsExpNoSExclMHTp1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
501	>	Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
502	>	Smooth( gCLsExpNoSExclHT, 27 );  gCLsExpNoSExclHT->SetLineWidth( 3 );
503	>	Smooth( gCLsExpNoSExclHTm1, 27 );
504	>	Smooth( gCLsExpNoSExclHTp1, 27 );
505	>	Smooth( gCLsExpNoSExclMHTm1, 27 );
506	>	Smooth( gCLsExpNoSExclMHTp1, 27 );
507	>	TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
508	>	TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
509	>	TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
510	>	TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
511	>	hexcl->Draw("colz");
512	>	LEP_ch->Draw("fsame");
513	>	LEP_sl->Draw("fsame");
514	>	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
515	>	sFirst->Draw("same");
516	>	Atlas->Draw("c,same");
517	>	ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
518	>	gCLsExpNoS1Sigma->Draw("lf,same");
519	>	sFirst->Draw("same");
520	>	Atlas->Draw("c,same");
521	>	CLsObsNLO->Draw("l,same");
522	>	CLsObsLO->Draw("l,same");
523	>	CLsExpNoSNLO->Draw("l,same");
524	>	//gCLsExpNoSExclMHT->Draw("l,same");
525	>	//gCLsExpNoSExclHT->Draw("l,same");
526	>	//gCLsObsExclHT->Draw("l,same");
527	>	//gobsexcl->Draw("l,same");
528	>	//FCExpNLO->Draw("l");
529	>	stau->Draw("fsame");
530	>	b.DrawLatex( 220,420,"#tilde{#tau} LSP");
531	>	//constant ssqquark and gluino lines
532	>	sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
533	>	sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
534	>	sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
535	>	glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
536	>	glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
537	>	glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
538	>	for(int i = 0; i < 4; i++){
539	>	lngl[i]->Draw("same");
540	>	lnsq[i]->Draw("same");
541	>	}
542	>	legexp->Draw();
543		leg->Draw();
544		gPad->RedrawAxis();
545	<	c1->SaveAs("results_tb50/Exclusion_m0_m12_tb50.pdf");
545	>	c1->SaveAs("results/Exclusion_m0_m12_tb50_NoSigHypPseudoData.pdf");
546	>
547	>
548
549		// Signal Contamination in M0 - M1/2
550		c1->SetLogz(0);
#	Line 624 | Line 803 | int plot(int argc, char** argv)
803		//c1->SaveAs("plot_tb10.pdf");
804
805
806	<
806	>	/*
807		//c1->SaveAs("plot_tb50.pdf");
808		c1->SetLogy(1);
809		c1->SetLogx(1);
#	Line 683 | Line 862 | int plot(int argc, char** argv)
862		leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
863		leg_S->Draw();
864		c1->SaveAs("results_tb50/UncertaintyScan.pdf");
865	<
865	>	*/
866		//theApp.Run();
867		}
868

Diff Legend

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