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.7 by auterman, Thu Mar 10 20:59:52 2011 UTC vs.
Revision 1.8 by auterman, Fri May 20 07:39:34 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_tb3_20110303/filelist.txt");
89	<	genpoints->Fill("limits_MHT_20110310/filelist_tb3.txt");
89	>	genpoints->Fill("limits_MHT_tb3add/filelist_tb3.txt");
90
91		TheLimits * genpointsHT = new TheLimits();
92		//genpointsHT->Fill("limits_HT_tb3_20110303/filelist.txt");
93	<	genpointsHT->Fill("limits_HT_20110310/filelist_tb3.txt");
93	>	genpointsHT->Fill("limits_HT_tb3add/filelist_tb3.txt");
94
95
96		//Replace read limits with specific numbers
#	Line 124 | 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 136 | 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 144 | 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 152 | 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 166 | 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 175 | 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 185 | 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	<	plotTools->Area(texpexcl, Mzero, Mhalf, ObsExclCL);
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 223 | 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,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);
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);
236	<	TGraph * gobsexcl_LO      = plotTools  ->GetContour(hs,Mzero,Mhalf,ObsExclCL,       3,0, 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 242 | 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, 2,1);
249		TGraph * gCLsExpExclHTm1  = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,2, 5,2);
#	Line 408 | Line 408 | int plot(int argc, char** argv)
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	+
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 422 | 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 434 | 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 446 | 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 459 | 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 471 | 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 484 | 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 496 | 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