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.7 by auterman, Thu Mar 10 20:59:52 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_20110310/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_20110310/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 197 | Line 204 | int plot(int argc, char** argv)
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);
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	>	plotTools->Area(texpexcl, Mzero, Mhalf, ObsExclCL);
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");
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	>	50,0,509.9,35,100,500);
232	>	TH2F*hs = new TH2F("hs","",50,0,509.9,35,100,450);
233	>	TGraph * gexpexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,2, 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,2, 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,2, 1,2);
246	>	TGraph * gCLsObsExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,0, 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");
286	<	TEV_sg_cdf->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	<	stau->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");
291	–	sSecond->Draw("same");
292	–	sThird->Draw("same");
338		TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
339	<	b.DrawLatex( 10,420,"#tilde{#tau} LSP");
340	<	b.DrawLatex(  6,150,"LEP2");
341	<	b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
342	<	b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
343	<	b.DrawLatex(100,140,"CDF");
344	<	b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
345	<	b.DrawLatex( 80,180,"D0 #tilde{#nu}");
346	<	TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
347	<	leg->SetBorderSize(0);leg->SetFillColor(0);
348	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
349	<	if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs, TLimit)","l");
350	<	if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs, TLimit)","l");
351	<	if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs, TLimit)","l");
352	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
353	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
354	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
355	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
356	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
357	<	if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
358	<	if (sFirst)  leg->AddEntry(sFirst, "NLO RA1 Observed");
359	<	if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected");
360	<	if (sThird) leg->AddEntry(sThird,"LO RA1 Observed");
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_tb3/Exclusion_m0_m12_tb3.pdf");
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		// Signal Contamination in M0 - M1/2
412		c1->SetLogz(0);
#	Line 339 | Line 430 | int plot(int argc, char** argv)
430		if (gobsexcl)    gobsexcl->Draw("l");
431		c1->SaveAs("results_tb3/SignalContamination_MHT_m0_m12_tb3.pdf");
432
433	+
434	+	// Signal JEC Uncertainty  MHT
435	+	c1->SetLogz(0);
436	+	TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
437	+	50,0,509.9,25,100,350);
438	+	plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
439	+	hsig_jec_mht->SetMinimum(0.0);
440	+	hsig_jec_mht->SetMaximum(0.25);
441	+	hsig_jec_mht->SetContour(10);
442	+	hsig_jec_mht->Draw("colz");
443	+	if (gobsexcl)    gobsexcl->Draw("l");
444	+	c1->SaveAs("results_tb3/SigJEC_MHT_m0_m12_tb3.pdf");
445	+
446	+	// Signal JEC Uncertainty  HT
447	+	c1->SetLogz(0);
448	+	TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
449	+	50,0,509.9,25,100,350);
450	+	plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
451	+	hsig_jec_ht->SetMinimum(0.0);
452	+	hsig_jec_ht->SetMaximum(0.25);
453	+	hsig_jec_ht->SetContour(10);
454	+	hsig_jec_ht->Draw("colz");
455	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
456	+	c1->SaveAs("results_tb3/SigJEC_HT_m0_m12_tb3.pdf");
457	+
458	+
459	+	// Signal MuIso Uncertainty  MHT
460	+	c1->SetLogz(0);
461	+	TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
462	+	50,0,509.9,25,100,350);
463	+	plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
464	+	hsig_MuIso_mht->SetMinimum(0.0);
465	+	hsig_MuIso_mht->SetMaximum(0.05);
466	+	hsig_MuIso_mht->SetContour(5);
467	+	hsig_MuIso_mht->Draw("colz");
468	+	if (gobsexcl)    gobsexcl->Draw("l");
469	+	c1->SaveAs("results_tb3/SigMuIso_MHT_m0_m12_tb3.pdf");
470	+
471	+	// Signal MuIso Uncertainty  HT
472	+	c1->SetLogz(0);
473	+	TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
474	+	50,0,509.9,25,100,350);
475	+	plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
476	+	hsig_MuIso_ht->SetMinimum(0.0);
477	+	hsig_MuIso_ht->SetMaximum(0.05);
478	+	hsig_MuIso_ht->SetContour(5);
479	+	hsig_MuIso_ht->Draw("colz");
480	+	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
481	+	c1->SaveAs("results_tb3/SigMuIso_HT_m0_m12_tb3.pdf");
482	+
483
484	<	/*
484	>	// Signal kFactor Uncertainty  MHT
485	>	c1->SetLogz(0);
486	>	TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
487	>	50,0,509.9,25,100,350);
488	>	plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
489	>	hsig_kFactor_mht->SetMinimum(0.1);
490	>	hsig_kFactor_mht->SetMaximum(0.2);
491	>	hsig_kFactor_mht->SetContour(10);
492	>	hsig_kFactor_mht->Draw("colz");
493	>	if (gobsexcl)    gobsexcl->Draw("l");
494	>	c1->SaveAs("results_tb3/SigkFactor_MHT_m0_m12_tb3.pdf");
495	>
496	>	// Signal kFactor Uncertainty  HT
497	>	c1->SetLogz(0);
498	>	TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
499	>	50,0,509.9,25,100,350);
500	>	plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
501	>	hsig_kFactor_ht->SetMinimum(0.1);
502	>	hsig_kFactor_ht->SetMaximum(0.2);
503	>	hsig_kFactor_ht->SetContour(10);
504	>	hsig_kFactor_ht->Draw("colz");
505	>	if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
506	>	c1->SaveAs("results_tb3/SigkFactor_HT_m0_m12_tb3.pdf");
507	>
508		//-----------------------------------------------------------------------------------
509		c1->SetLogz(1);
510
#	Line 349 | Line 513 | int plot(int argc, char** argv)
513
514		// cross-section in squark - gluino mass
515		TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
516	<	60,200,1400,50,200,1200);
516	>	60,200,800,40,200,900);
517		plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
518		hxsec_qg->SetMinimum(0.01);
519		hxsec_qg->Draw("colz");
#	Line 359 | Line 523 | int plot(int argc, char** argv)
523
524		// Observed Limit in squark - gluino mass
525		TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
526	<	60,200,1400,50,200,1200);
526	>	60,200,800,40,200,900);
527		plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
528		hobslimit_qg->SetMinimum(0.01);
529		hobslimit_qg->Draw("colz");
#	Line 367 | Line 531 | int plot(int argc, char** argv)
531
532		// Expected Limit in squark - gluino mass
533		TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
534	<	60,200,1400,50,200,1200);
534	>	60,200,800,40,200,900);
535		plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
536		hexplimit_qg->SetMinimum(0.01);
537		hexplimit_qg->Draw("colz");
#	Line 375 | Line 539 | int plot(int argc, char** argv)
539
540		// Signal Acceptance in squark - gluino mass
541		TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
542	<	60,200,1400,50,200,1200);
542	>	60,200,800,40,200,900);
543		plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
544		hsigacc_qg->SetMinimum(0.01);
545		hsigacc_qg->SetMaximum(1.0);
#	Line 385 | Line 549 | int plot(int argc, char** argv)
549		// Exp. Limit on Number of Signal Events in squark - gluino mass
550		c1->SetLogz(0);
551		TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
552	<	60,200,1400,50,200,1200);
552	>	60,200,800,40,200,900);
553		plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
554		hexplimitnsig_qg->SetMinimum(0.0);
555		hexplimitnsig_qg->SetMaximum(20);
#	Line 394 | Line 558 | int plot(int argc, char** argv)
558
559		// Obs. Limit on Number of Signal Events in squark - gluino mass
560		TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
561	<	60,200,1400,50,200,1200);
561	>	60,200,800,40,200,900);
562		plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
563		hobslimitnsig_qg->SetMinimum(0.0);
564		hobslimitnsig_qg->SetMaximum(20);
#	Line 404 | Line 568 | int plot(int argc, char** argv)
568		c1->SetLogz(0);
569		// Expected Exclusion in squark - gluino mass
570		TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
571	<	60,200,1400,50,200,1200);
571	>	60,200,800,40,200,900);
572		plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
573		hexpexcl_qg->Draw("colz");
574		c1->SaveAs("results_tb3/ExpExclusion_mSql_mGl_tb3.pdf");
575
576		// Observed Exclusion in squark - gluino mass
577		TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
578	<	60,200,1400,50,200,1200);
578	>	60,200,800,40,200,900);
579		plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
580		//std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
581		hobsexcl_qg->Draw("colz");
#	Line 423 | Line 587 | int plot(int argc, char** argv)
587
588		// TestContours in M0 - M1/2
589		TH2F*texcl_qg = new TH2F("texcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
590	<	60,200,1400,50,200,1200);
590	>	60,200,800,40,200,900);
591		TH2F*tPLobsexcl_qg=(TH2F*)texcl_qg->Clone();
592		plotTools->Area(tPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
593		TH2F*tPLexpexcl_qg=(TH2F*)texcl_qg->Clone();
#	Line 454 | Line 618 | int plot(int argc, char** argv)
618
619		// Observed Exclusion in squark - gluino mass
620		TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
621	<	60,200,1400,50,200,1200);
621	>	60,200,800,40,200,900);
622		plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
623		//hPLobsexcl_qg->Draw("colz");
624		//c1->SaveAs("results_tb3/PL_ObsExclusion_mSql_mGl_tb3.pdf");
625
626
627
464	–
465	–
466	–	// Exclusion in squark - gluino mass ----------------------------------------
628		TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
629	<	60,200,1400,50,200,1200);
630	<	TGraph * gexpexcl_qg = plotTools->GetContour(hexpexcl_qg,3,0);
631	<	if (gexpexcl_qg) gexpexcl_qg->SetLineStyle(2);
632	<	TGraph * gobsexcl_qg = plotTools->GetContour(hobsexcl_qg,3,0);
633	<	TH2F*hPLexpexcl_qg=(TH2F*)hexcl_qg->Clone();
634	<	plotTools->Area(hPLexpexcl_qg, MSquarkL, MGluino, PLExpExclusion);
635	<	TGraph * gPLexpexcl_qg = plotTools->GetContour(hPLexpexcl_qg,3,0);
636	<	if (gPLexpexcl_qg)   {gPLexpexcl_qg->SetLineStyle(2);gPLexpexcl_qg->SetLineColor(2);}
637	<	TGraph * gPLobsexcl_qg = plotTools->GetContour(hPLobsexcl_qg,3,0);
638	<	if (gPLobsexcl_qg) gPLobsexcl_qg->SetLineColor(2);
639	<	TH2F*hCLsExpHT_qg=(TH2F*)hexcl_qg->Clone();
640	<	TH2F*hCLsObsHT_qg=(TH2F*)hexcl_qg->Clone();
641	<	plotToolsHT->Area(hCLsExpHT_qg, MSquarkL, MGluino, ExpExclCL);
642	<	plotToolsHT->Area(hCLsObsHT_qg, MSquarkL, MGluino, ObsExclCL);
643	<	TGraph * gCLsExpExclHT_qg = plotTools->GetContour(hCLsExpHT_qg,3,0);
644	<	TGraph * gCLsObsExclHT_qg = plotTools->GetContour(hCLsObsHT_qg,3,0);
645	<	if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->SetLineColor(2);
646	<	if (gCLsExpExclHT_qg) {gCLsExpExclHT_qg->SetLineStyle(2);gCLsExpExclHT_qg->SetLineColor(2);}
629	>	50,200,800,30,200,900);
630	>	TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 2,2);
631	>	TGraph * gexpexcl_qg_LO      = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 2,4);
632	>	TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 2,1);
633	>	TGraph * gMCMCexpexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
634	>	TGraph * gMCMCobsexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
635	>	TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 1,2);
636	>	TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,4);
637	>	TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 1,1);
638	>	TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
639	>	TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
640	>	Smooth(gexpexcl_qg, 11);
641	>	Smooth(gobsexcl_qg, 11);
642	>	Smooth(gCLsExpExclHT_qg, 11);
643	>	Smooth(gCLsObsExclHT_qg, 11);
644	>	Smooth(gCLsExpExclHTm1_qg, 11);
645	>	Smooth(gCLsExpExclHTp1_qg, 11);
646	>	TGraph * gCLsExp1Sigma_qg    = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
647		hexcl_qg->Draw("");
648	<	if (gobsexcl_qg) gobsexcl_qg->Draw("l");
648	>	//
649	>	if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
650	>	if (gobsexcl_qg)      gobsexcl_qg->Draw("l");
651		if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
652	<	if (gexpexcl_qg) gexpexcl_qg->Draw("l");
490	<	//if (gPLobsexcl_qg) gPLobsexcl_qg->Draw("l");
491	<	//if (gPLexpexcl_qg) gPLexpexcl_qg->Draw("l");
652	>	if (gexpexcl_qg)      gexpexcl_qg->Draw("l");
653		if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
654	+	if (gexpexcl_qg_LO)      gexpexcl_qg_LO->Draw("l");
655	+	if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
656	+	//
657		TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
658		leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
659	<	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"Observed (MHT, CLs, TLimit)","l");
660	<	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"Expected (MHT, CLs, TLimit)","l");
661	<	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"Observed (HT, CLs, TLimit)","l");
662	<	if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"Expected (HT, CLs, TLimit)","l");
659	>	TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
660	>	if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
661	>	if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
662	>	if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
663	>	if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
664	>	if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf");
665	>	if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l");
666		//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
667		//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
668		leg_qg->Draw();
669		c1->SaveAs("results_tb3/Exclusion_mSql_mGl_tb3.pdf");
670	<	*/
671	<	//c1->SaveAs("plot_tb3.pdf");
670	>	//c1->SaveAs("plot_tb10.pdf");
671	>
672	>
673
674		c1->SetLogy(1);
675		c1->SetLogx(1);

Diff Legend

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