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.7 by auterman, Thu Mar 10 20:59:52 2011 UTC vs.
Revision 1.9 by auterman, Wed Jun 22 15:03:36 2011 UTC

#	Line 86 | Line 86 | int plot(int argc, char** argv)
86		TheLimits * genpoints = new TheLimits();
87		//genpoints->Fill(argc, argv);
88		//genpoints->Fill("limits_MHT_tb50_20110303/filelist.txt");
89	<	genpoints->Fill("limits_MHT_20110310/filelist_tb50.txt");
89	>	genpoints->Fill("limits_moriond_MHT/filelist_tb50.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92		//genpointsHT->Fill("limits_HT_tb50_20110303/filelist.txt");
93	<	genpointsHT->Fill("limits_HT_20110310/filelist_tb50.txt");
93	>	genpointsHT->Fill("limits_moriond_HT/filelist_tb50.txt");
94
95
96		//Replace read limits with specific numbers
#	Line 101 | 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 248 | 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);
#	Line 350 | Line 359 | int plot(int argc, char** argv)
359		sFirst->Draw("same");
360		//sSecond->Draw("same");
361		//sThird->Draw("same");
362	<	Atlas->Draw("c,same");
362	>	//Atlas->Draw("c,same");
363		TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
364		TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
365	<	ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
366	<	TLegend* legexp = new TLegend(0.69,0.77,0.93,0.92,NULL,"brNDC");
367	<	legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
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 preliminary");
371	<	//legexp->AddEntry(TEV_sg_cdf,"CDF  #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
372	<	//legexp->AddEntry(TEV_sg_d0,"D0   #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=3, #mu<0}","f");
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);
#	Line 367 | Line 378 | int plot(int argc, char** argv)
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");
381	>	//legexp->AddEntry(Atlas,  "Atlas, #scale[0.8]{tan #beta=3}","l");
382		legexp->Draw();
383
384		CLsObsNLO->Draw("l");
#	Line 382 | Line 393 | int plot(int argc, char** argv)
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,255,"#font[92]{#tilde{q}(650)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}");
#	Line 395 | Line 406 | int plot(int argc, char** argv)
406		lngl[i]->Draw("same");
407		lnsq[i]->Draw("same");
408		}
409	<	TLegend * leg = new TLegend(0.38,0.8,0.73,0.92);
410	<	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
409	>	TLegend * leg = new TLegend(0.38,0.72,0.81,0.88);
410	>	//TLegend * leg = new TLegend(0.25,0.72,0.68,0.88);
411	>	leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);leg->SetTextSize(0.025);
412		TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
413	<	leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
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");
#	Line 421 | Line 433 | int plot(int argc, char** argv)
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);
556		TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
#	Line 678 | 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 737 | 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