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.10 by auterman, Wed Jun 22 16:02:08 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.old/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.old/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 226 | 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",
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);
#	Line 243 | Line 275 | int plot(int argc, char** argv)
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, 1,1);
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, 33); gCLsObsExclHT->SetLineWidth( 3 );
286	<	Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
287	<	Smooth( gCLsObsExclHT_LO, 33); gCLsObsExclHT_LO->SetLineWidth( 3 );
288	<	Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
289	<	Smooth( gexpexcl, 33 );
290	<	Smooth( gobsexcl, 33 );
291	<	Smooth( gobsexcl, 33 );
292	<	Smooth( gexpexcl_LO, 33 );
293	<	Smooth( gFCexpexclHT, 33 );
294	<	Smooth( gFCexpexclHT_LO, 33 );gFCexpexclHT_LO->SetLineWidth( 3 );
295	<	Smooth( gFCexpexclMHT, 33 );
296	<	Smooth( gFCexpexclMHT_LO, 33 );
297	<	Smooth( gFCobsexclHT, 33 );gFCobsexclHT->SetLineWidth( 3 );
298	<	Smooth( gFCobsexclMHT, 33 );
299	<	Smooth( gCLsExpExclHTm1, 33 );
300	<	Smooth( gCLsExpExclHTp1, 33 );
301	<	Smooth( gCLsExpExclMHTm1, 33 );
302	<	Smooth( gCLsExpExclMHTp1, 33 );
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");
#	Line 284 | Line 320 | int plot(int argc, char** argv)
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 311 | Line 349 | int plot(int argc, char** argv)
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");
316	<	TEV_sg_d0->Draw("fsame");
317	<	TEV_sg_cdf->Draw("fsame");
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	<	stau->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");
360	>	//sSecond->Draw("same");
361	>	//sThird->Draw("same");
362	>	//Atlas->Draw("c,same");
363		TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
364	<	//b.DrawLatex( 10,420,"#tilde{#tau} LSP");
365	<	b.DrawLatex(230,145,"LEP2");
366	<	//b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
367	<	//b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
368	<	//b.DrawLatex(100,140,"CDF");
369	<	//b.DrawLatex(230,145,"D0 #tilde{g}, #tilde{q}");
370	<	//b.DrawLatex( 80,180,"D0 #tilde{#nu}");
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	<	CLsExpLO->Draw("l");
388	<	FCExpLO->Draw("l");
389	<	//FCObsNLO->Draw("l");
390	<
391	<	//gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c");
392	<	//sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c");
393	<	//TLegend * leg = new TLegend(0.45,0.65,0.84,0.89);
394	<	//leg->SetBorderSize(0);leg->SetFillColor(0);
395	<	TLegend * leg = new TLegend(0.6,0.66,0.96,0.90);
396	<	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);
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("CMS preliminary, 95%CL");
414	<	leg->AddEntry(expLeg,   "Expected (CLs, NLO)","l");
415	<	leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l");
416	<	leg->AddEntry(CLsExpLO, "Expected (CLs, LO)","l");
417	<	leg->AddEntry(FCExpLO,  "Expected (FC, LO)","l");
350	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
351	<	//if (gCLsExpExclHT) leg->AddEntry(expLeg       ,"NLO Expected (HT, CLs)","l");
352	<	//if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
353	<	//if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
354	<	//if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
355	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
356	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
357	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l");
358	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l");
359	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
360	<	//if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
361	<	//if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
362	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l");
363	<	if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)");
364	<	if (sFirst)  leg->AddEntry(sFirst, "#alpha_{T} Observed (FC, NLO)");
365	<	if (sThird)  leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)");
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_tb50/Exclusion_m0_m12_tb50.pdf");
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	>	//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/Exclusion_m0_m12_tb50_NoSigHypPseudoData.pdf");
551	>
552	>
553
554		// Signal Contamination in M0 - M1/2
555		c1->SetLogz(0);
#	Line 624 | Line 808 | int plot(int argc, char** argv)
808		//c1->SaveAs("plot_tb10.pdf");
809
810
811	<
811	>	/*
812		//c1->SaveAs("plot_tb50.pdf");
813		c1->SetLogy(1);
814		c1->SetLogx(1);
#	Line 683 | 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