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.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_tb3-20110301/filelist.txt");
88	>	//genpoints->Fill("limits_MHT_tb3_20110303/filelist.txt");
89	>	genpoints->Fill("limits_MHT_tb3add/filelist_tb3.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92	<	genpointsHT->Fill("limits_HT_tb3-20110301/filelist.txt");
92	>	//genpointsHT->Fill("limits_HT_tb3_20110303/filelist.txt");
93	>	genpointsHT->Fill("limits_HT_tb3add/filelist_tb3.txt");
94
95
96		//Replace read limits with specific numbers
#	Line 122 | Line 124 | int plot(int argc, char** argv)
124		/**/
125		// cross-section in M0 - M1/2
126		TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
127	<	50,0,509.9,25,100,350);
127	>	150,0,1509.9,25,100,350);
128		plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
129		hxsec->SetMinimum(0.01);
130		//sq500->Draw();
#	Line 134 | Line 136 | int plot(int argc, char** argv)
136
137		// Observed Limit in M0 - M1/2
138		TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
139	<	50,0,509.9,25,100,350);
139	>	150,0,1509.9,25,100,350);
140		plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
141		hobslimit->SetMinimum(0.01);
142		hobslimit->Draw("colz");
#	Line 142 | Line 144 | int plot(int argc, char** argv)
144
145		// Expected Limit in M0 - M1/2
146		TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
147	<	50,0,509.9,25,100,350);
147	>	150,0,1509.9,25,100,350);
148		plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
149		hexplimit->SetMinimum(0.01);
150		hexplimit->Draw("colz");
#	Line 150 | Line 152 | int plot(int argc, char** argv)
152
153		// Signal Acceptance in M0 - M1/2
154		TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
155	<	50,0,509.9,25,100,350);
155	>	150,0,1509.9,25,100,350);
156		plotTools->Area(hsigacc, Mzero, Mhalf, SignalAcceptance);
157		hsigacc->SetMinimum(0.01);
158		hsigacc->SetMaximum(1.0);
#	Line 164 | Line 166 | int plot(int argc, char** argv)
166		// Exp. Limit on Number of Signal Events in M0 - M1/2
167		c1->SetLogz(0);
168		TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
169	<	50,0,509.9,25,100,350);
169	>	150,0,1509.9,25,100,350);
170		plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
171		hexplimitnsig->SetMinimum(0.0);
172		hexplimitnsig->SetMaximum(20);
#	Line 173 | Line 175 | int plot(int argc, char** argv)
175
176		// Obs. Limit on Number of Signal Events in M0 - M1/2
177		TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
178	<	50,0,509.9,25,100,350);
178	>	150,0,1509.9,25,100,350);
179		plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
180		hobslimitnsig->SetMinimum(0.0);
181		hobslimitnsig->SetMaximum(20);
#	Line 183 | Line 185 | int plot(int argc, char** argv)
185		c1->SetLogz(0);
186		// Expected Exclusion in M0 - M1/2
187		TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
188	<	50,0,509.9,35,100,450);
188	>	150,0,1509.9,35,100,450);
189		plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
190		hexpexcl->Draw("colz");
191		c1->SaveAs("results_tb3/ExpExclusion_m0_m12_tb3.pdf");
192
193		// Observed Exclusion in M0 - M1/2
194		TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
195	<	50,0,509.9,35,100,450);
195	>	150,0,1509.9,35,100,450);
196		plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
197		hobsexcl->Draw("colz");
198		c1->SaveAs("results_tb3/ObsExclusion_m0_m12_tb3.pdf");
199
200		// Observed Exclusion in M0 - M1/2
201		TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
202	<	50,0,509.9,35,100,450);
202	>	150,0,1509.9,35,100,450);
203		plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
204		hPLobsexcl->Draw("colz");
205		c1->SaveAs("results_tb3/PL_ObsExclusion_m0_m12_tb3.pdf");
206
207		// TestContours in M0 - M1/2  ===================================================================
208		TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
209	<	50,0,509.9,35,100,450);
209	>	150,0,1509.9,35,100,450);
210		TH2F*texpexcl=(TH2F*)texcl->Clone();
211	<	plotToolsHT->Area(texpexcl, Mzero, Mhalf, NLOExpExclCLm1sigma);
211	>	plotToolsHT->Area(texpexcl, Mzero, Mhalf, NLOExpNoSigExclCL);
212		std::vector<TGraph*> contours = plotToolsHT->GetContours(texpexcl,3);
213		texpexcl->Draw("colz");
214		int col=kBlue-10;
#	Line 221 | Line 223 | int plot(int argc, char** argv)
223		sprintf(val,"%d",(int)(cont-contours.begin()));
224		l.DrawLatex(x,y,val);
225		}
226	<	c1->SaveAs("results_tb3/ExclusionTestContours_m0_m12_tb3.pdf");
226	>	c1->SaveAs("results/ExclusionTestContours_m0_m12_tb3.pdf");
227
228
229		// Exclusion in M0 - M1/2
230		TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
231	<	50,0,509.9,35,100,450);
232	<	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);
231	>	150,0,1509.9,35,100,500);
232	>	TH2F*hs = new TH2F("hs","",150,0,1509.9,35,100,450);
233	>	TGraph * gexpexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,3, 2,2);
234		TGraph * gexpexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 2,4);
235	<	TGraph * gobsexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,2, 2,1);
235	>	TGraph * gobsexcl         = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,3, 2,1);
236	>	TGraph * gobsexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 2,1);
237		TGraph * gMCMCexpexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion,   3,0, 3,2);
238		TGraph * gMCMCobsexcl     = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion,   3,0, 3,1);
239		TGraph * gFCexpexclHT     = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion,  3,0, 3,2);
#	Line 239 | Line 242 | int plot(int argc, char** argv)
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);
245	>	TGraph * gCLsExpExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL,       3,3, 1,2);
246	>	TGraph * gCLsObsExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL,       3,3, 1,1);
247		TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL,          3,0, 1,4);
248	<	TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL,          3,0, 1,1);
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, 33); gCLsObsExclHT->SetLineWidth( 3 );
256	<	Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
257	<	Smooth( gCLsObsExclHT_LO, 33); gCLsObsExclHT_LO->SetLineWidth( 3 );
258	<	Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
259	<	Smooth( gexpexcl, 33 );
260	<	Smooth( gobsexcl, 33 );
261	<	Smooth( gobsexcl, 33 );
262	<	Smooth( gexpexcl_LO, 33 );
263	<	Smooth( gFCexpexclHT, 33 );
264	<	Smooth( gFCexpexclHT_LO, 33 );gFCexpexclHT_LO->SetLineWidth( 3 );
265	<	Smooth( gFCexpexclMHT, 33 );
266	<	Smooth( gFCexpexclMHT_LO, 33 );
267	<	Smooth( gFCobsexclHT, 33 );gFCobsexclHT->SetLineWidth( 3 );
268	<	Smooth( gFCobsexclMHT, 33 );
269	<	Smooth( gCLsExpExclHTm1, 33 );
270	<	Smooth( gCLsExpExclHTp1, 33 );
271	<	Smooth( gCLsExpExclMHTm1, 33 );
272	<	Smooth( gCLsExpExclMHTp1, 33 );
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* 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
#	Line 317 | Line 326 | int plot(int argc, char** argv)
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");
322	<	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");
338	<	sSecond->Draw("same");
339	<	sThird->Draw("same");
338	>	TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
339	>	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
340	>
341	>	//sSecond->Draw("same");
342	>	//sThird->Draw("same");
343	>	Atlas->Draw("c,same");
344	>	TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
345	>	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
346	>	//gCheck1->Draw("same");
347	>	//gCheck2->Draw("same");
348	>	TLegend* legexp = new TLegend(0.66,0.65,0.98,0.93,NULL,"brNDC");
349	>	legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
350	>	//TEV_sg_cdf.SetLineColor(1);
351	>	legexp->SetHeader("CMS preliminary");
352	>	legexp->AddEntry(TEV_sg_cdf,"CDF  #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
353	>	legexp->AddEntry(TEV_sg_d0,"D0   #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{#mu<0}","f");
354	>	//ch_gr.SetLineColor(1);
355	>	legexp->AddEntry(LEP_ch,"LEP2   #tilde{#chi}_{1}^{#pm}","f");
356	>	//sl_gr.SetLineColor(1);
357	>	legexp->AddEntry(LEP_sl,"LEP2   #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
358	>	//if(tanbeta == 3)
359	>	legexp->AddEntry(TEV_sn_d0_1,"D0  #chi^{#pm}_{1}, #chi^{0}_{2}","f");
360	>	legexp->AddEntry(sFirst, "CMS #alpha_{T}");
361	>	legexp->AddEntry(Atlas,  "Atlas","l");
362	>	legexp->Draw();
363	>
364		CLsObsNLO->Draw("l");
365	+	CLsObsLO->Draw("l");
366		CLsExpNLO->Draw("l");
367	<	CLsExpLO->Draw("l");
332	<	FCExpLO->Draw("l");
333	<	//FCObsNLO->Draw("l");
334	<
335	<	//gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c");
336	<	//sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c");
337	<	//TLegend * leg = new TLegend(0.55,0.7,0.84,0.89);
338	<	//leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetTextFont(42);
339	<	TLegend * leg = new TLegend(0.6,0.66,0.96,0.90);
340	<	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);
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("CMS preliminary, 95%CL");
390	<	leg->AddEntry(expLeg,   "Expected (CLs, NLO)","l");
391	<	leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l");
392	<	leg->AddEntry(CLsExpLO, "Expected (CLs, LO)","l");
393	<	leg->AddEntry(FCExpLO,  "Expected (FC, LO)","l");
348	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
349	<	//if (gCLsExpExclHT) leg->AddEntry(expLeg       ,"NLO Expected (HT, CLs)","l");
350	<	//if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
351	<	//if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l");
352	<	//if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
353	<	//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
354	<	//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
355	<	//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l");
356	<	//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l");
357	<	//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
358	<	//if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
359	<	//if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l");
360	<	//if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l");
361	<	if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)");
362	<	if (sFirst)  leg->AddEntry(sFirst, "#alpha_{T} Observed (FC, NLO)");
363	<	if (sThird)  leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)");
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	>
412	>
413	>	hexcl->Draw("colz");
414	>	LEP_ch->Draw("fsame");
415	>	LEP_sl->Draw("fsame");
416	>	TEV_sg_cdf->Draw("fsame");
417	>	TEV_sg_cdf->Draw("lsame");
418	>	TEV_sg_d0->Draw("fsame");
419	>	TEV_sg_d0->Draw("lsame");
420	>	if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
421	>	sFirst->Draw("same");
422	>	Atlas->Draw("c,same");
423	>	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
424	>	CLsObsNLO->Draw("l");
425	>	//CLsObsLO->Draw("l");
426	>	CLsExpNLO->Draw("l");
427	>	//FCExpNLO->Draw("l");
428	>	//Jim's limits///////////////////////////////////////////////////////////////////////////////
429	>	TGraph* JimObsHT  = Jim_ht_tb3(0);
430	>	TGraph* JimObsMHT = Jim_mht_tb3(0);
431	>	TGraph* JimExpHT  = Jim_ht_tb3(1);
432	>	TGraph* JimExpMHT = Jim_mht_tb3(1);
433	>	TGraph* JimExpHTup  = Jim_ht_tb3(2);
434	>	TGraph* JimExpMHTup = Jim_mht_tb3(2);
435	>	TGraph* JimExpHTdn  = Jim_ht_tb3(3);
436	>	TGraph* JimExpMHTdn = Jim_mht_tb3(3);
437	>	Smooth(JimObsHT, 2);
438	>	Smooth(JimObsMHT,2);
439	>	Smooth(JimExpHT, 2);
440	>	Smooth(JimExpMHT,2);
441	>	Smooth(JimExpHTup, 2);
442	>	Smooth(JimExpMHTup,2);
443	>	Smooth(JimExpHTdn, 2);
444	>	Smooth(JimExpMHTdn,2);
445	>	TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
446	>	TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
447	>	TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
448	>	TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
449	>	TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
450	>	TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
451	>	JimExp1Sigma->Draw("f");
452	>	JimObs->Draw("c");
453	>	JimExp->Draw("c");
454	>	stau->Draw("fsame");
455	>	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
456	>	//constant ssqquark and gluino lines
457	>	sqt.DrawLatex(148,218,"#font[92]{#tilde{q}(500)GeV}");
458	>	sqt.DrawLatex(220,385,"#font[92]{#tilde{q}(800)GeV}");
459	>	glt.DrawLatex(430,184,"#font[92]{#tilde{g}(500)GeV}");
460	>	glt.DrawLatex(428,311,"#font[92]{#tilde{g}(800)GeV}");
461	>	for(int i = 0; i < 4; i++){
462	>	lngl[i]->Draw("same");
463	>	lnsq[i]->Draw("same");
464	>	}
465	>	legexp->Draw();
466	>	TLegend * legBayes = new TLegend(0.3,0.8,0.65,0.93);
467	>	legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
468	>	legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
469	>	legBayes->AddEntry(JimObs,   "Observed, Bayes","l");
470	>	legBayes->AddEntry(JimLeg,   "Expected #pm 1#sigma, Bayes","lf");
471	>	legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
472	>	legBayes->AddEntry(expLeg,   "Expected #pm 1#sigma, CLs","lf");
473	>	legBayes->Draw();
474	>	gPad->RedrawAxis();
475	>	c1->SaveAs("results/Exclusion_m0_m12_tb3_Bayes.pdf");
476	>
477	>
478	>	//The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
479	>	TGraph * gCLsExpNoSExclMHT   = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,2, 1,2);
480	>	TGraph * gCLsExpNoSExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,2, 1,2);
481	>	TGraph * gCLsExpNoSExclHTm1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
482	>	TGraph * gCLsExpNoSExclHTp1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
483	>	TGraph * gCLsExpNoSExclMHTm1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
484	>	TGraph * gCLsExpNoSExclMHTp1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
485	>	Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
486	>	Smooth( gCLsExpNoSExclHT, 27 );  gCLsExpNoSExclHT->SetLineWidth( 3 );
487	>	Smooth( gCLsExpNoSExclHTm1, 27 );
488	>	Smooth( gCLsExpNoSExclHTp1, 27 );
489	>	Smooth( gCLsExpNoSExclMHTm1, 27 );
490	>	Smooth( gCLsExpNoSExclMHTp1, 27 );
491	>	TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
492	>	TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
493	>	TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
494	>	TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
495	>	hexcl->Draw("colz");
496	>	LEP_ch->Draw("fsame");
497	>	LEP_sl->Draw("fsame");
498	>	TEV_sg_cdf->Draw("fsame");
499	>	TEV_sg_cdf->Draw("lsame");
500	>	TEV_sg_d0->Draw("fsame");
501	>	TEV_sg_d0->Draw("lsame");
502	>	ms.DrawLatex(395,508,"tan#beta=3, #mu>0, A_{0}=0");
503	>	gCLsExpNoS1Sigma->Draw("lf,same");
504	>	sFirst->Draw("same");
505	>	Atlas->Draw("c,same");
506	>	CLsObsNLO->Draw("l,same");
507	>	CLsObsLO->Draw("l,same");
508	>	CLsExpNoSNLO->Draw("l,same");
509	>	//FCExpNLO->Draw("l");
510	>	stau->Draw("fsame");
511	>	b.DrawLatex( 10,435,"#tilde{#tau} LSP");
512	>	//constant ssqquark and gluino lines
513	>	sqt.DrawLatex(148,218,"#font[92]{#tilde{q}(500)GeV}");
514	>	sqt.DrawLatex(220,385,"#font[92]{#tilde{q}(800)GeV}");
515	>	glt.DrawLatex(430,184,"#font[92]{#tilde{g}(500)GeV}");
516	>	glt.DrawLatex(428,311,"#font[92]{#tilde{g}(800)GeV}");
517	>	for(int i = 0; i < 4; i++){
518	>	lngl[i]->Draw("same");
519	>	lnsq[i]->Draw("same");
520	>	}
521	>	legexp->Draw();
522	>	leg->Draw();
523	>	gPad->RedrawAxis();
524	>	c1->SaveAs("results/Exclusion_m0_m12_tb3_NoSigHypPseudoData.pdf");
525	>
526
368	–	plotToolsHT->Print(Xsection, Mzero,Mhalf,gCLsObsExclHT, 10);
527
528		// Signal Contamination in M0 - M1/2
529		c1->SetLogz(0);
530		TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
531	<	50,0,509.9,25,100,350);
531	>	150,0,1509.9,25,100,350);
532		plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
533		hsigcont->SetMinimum(0.01);
534		hsigcont->SetMaximum(20);
#	Line 381 | Line 539 | int plot(int argc, char** argv)
539		// Signal Contamination in M0 - M1/2
540		c1->SetLogz(0);
541		TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
542	<	50,0,509.9,25,100,350);
542	>	150,0,1509.9,25,100,350);
543		plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
544		hsigcontMHT->SetMinimum(0.01);
545		hsigcontMHT->SetMaximum(20);
#	Line 393 | Line 551 | int plot(int argc, char** argv)
551		// Signal JEC Uncertainty  MHT
552		c1->SetLogz(0);
553		TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
554	<	50,0,509.9,25,100,350);
554	>	150,0,1509.9,25,100,350);
555		plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
556		hsig_jec_mht->SetMinimum(0.0);
557		hsig_jec_mht->SetMaximum(0.25);
#	Line 405 | Line 563 | int plot(int argc, char** argv)
563		// Signal JEC Uncertainty  HT
564		c1->SetLogz(0);
565		TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
566	<	50,0,509.9,25,100,350);
566	>	150,0,1509.9,25,100,350);
567		plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
568		hsig_jec_ht->SetMinimum(0.0);
569		hsig_jec_ht->SetMaximum(0.25);
#	Line 418 | Line 576 | int plot(int argc, char** argv)
576		// Signal MuIso Uncertainty  MHT
577		c1->SetLogz(0);
578		TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
579	<	50,0,509.9,25,100,350);
579	>	150,0,1509.9,25,100,350);
580		plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
581		hsig_MuIso_mht->SetMinimum(0.0);
582		hsig_MuIso_mht->SetMaximum(0.05);
#	Line 430 | Line 588 | int plot(int argc, char** argv)
588		// Signal MuIso Uncertainty  HT
589		c1->SetLogz(0);
590		TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
591	<	50,0,509.9,25,100,350);
591	>	150,0,1509.9,25,100,350);
592		plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
593		hsig_MuIso_ht->SetMinimum(0.0);
594		hsig_MuIso_ht->SetMaximum(0.05);
#	Line 443 | Line 601 | int plot(int argc, char** argv)
601		// Signal kFactor Uncertainty  MHT
602		c1->SetLogz(0);
603		TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
604	<	50,0,509.9,25,100,350);
604	>	150,0,1509.9,25,100,350);
605		plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
606		hsig_kFactor_mht->SetMinimum(0.1);
607		hsig_kFactor_mht->SetMaximum(0.2);
#	Line 455 | Line 613 | int plot(int argc, char** argv)
613		// Signal kFactor Uncertainty  HT
614		c1->SetLogz(0);
615		TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
616	<	50,0,509.9,25,100,350);
616	>	150,0,1509.9,25,100,350);
617		plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
618		hsig_kFactor_ht->SetMinimum(0.1);
619		hsig_kFactor_ht->SetMaximum(0.2);

Diff Legend

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