ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents):
Revision 1.56 by buchmann, Wed Sep 5 20:16:25 2012 UTC vs.
Revision 1.92 by buchmann, Tue Apr 30 09:49:05 2013 UTC

#	Line 33 | Line 33 | void todo() {
33
34
35
36	<	void find_one_peak_combination(TCut specialcut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="")
36	>	void find_one_peak_combination(TCut specialcut, bool SwitchOffNJetsCut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="")
37		{
38		// Temporarily switch off PU reweighting, if asked
39		TCut weightbackup=cutWeight;
#	Line 42 | Line 42 | void find_one_peak_combination(TCut spec
42		int nbins=100;
43		if(PlottingSetup::DoBTag) nbins=25;
44
45	+	TCut nJetsCut(cutnJets);
46	+	if(SwitchOffNJetsCut) nJetsCut=specialcut;
47	+
48		TCanvas *tempcan = new TCanvas("tempcan","Temporary canvas for peak finding preparations");
49	<	TH1F *rawJZBeemmMC      = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,mc, luminosity);
50	<	TH1F *rawJZBeemmData    = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,data, luminosity);
51	<	TH1F *rawJZBemMC      = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,mc, luminosity);
52	<	TH1F *rawJZBemData    = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,data, luminosity);
49	>	TH1F *rawJZBeemmMC      = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,mc, luminosity);
50	>	TH1F *rawJZBeemmData    = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,data, luminosity);
51	>	TH1F *rawJZBemMC      = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,mc, luminosity);
52	>	TH1F *rawJZBemData    = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,data, luminosity);
53		TH1F *rawttbarjzbeemmMC;
54
55		if(method==doKM) {
56		//we only need this histo for the KM fitting...
57	<	rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&specialcut,mc,luminosity,allsamples.FindSample("TTJet"));
57	>	rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&nJetsCut&&specialcut,mc,luminosity,allsamples.FindSample("TTJet"));
58		MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas);
59		DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas);
60		delete rawttbarjzbeemmMC;
#	Line 84 | Line 87 | void find_one_peak_combination(TCut spec
87		delete tempcan;
88		}
89
90	<	void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb)
90	>	void DoPeakVsPU() {
91	>	cout << "Trying to find peak position as a function of PU" << endl;
92	>	float MCPeak, MCPeakError, DataPeak, DataPeakError;
93	>	stringstream result, datajzb, mcjzb;
94	>	bool doPUreweighting=true;
95	>	bool SwitchOffNJetsCut=false;
96	>
97	>
98	>	TCanvas *can = new TCanvas("can","can");
99	>	float mcSigma,mcSigmaError,dataSigma,dataSigmaError;
100	>
101	>	float NumVtxBin[6] = {0,5,10,15,20,30};
102	>
103	>	stringstream NowCut;
104	>
105	>	TGraphErrors *gMCPeak = new TGraphErrors();
106	>	gMCPeak->SetTitle("gMCPeak");
107	>	gMCPeak->SetName("gMCPeak");
108	>	TGraphErrors *gDataPeak = new TGraphErrors();
109	>	gDataPeak->SetTitle("gDataPeak");
110	>	gDataPeak->SetName("gDataPeak");
111	>	TGraphErrors *gMCWidth = new TGraphErrors();
112	>	gMCWidth->SetTitle("gMCWidth");
113	>	gMCWidth->SetName("gMCWidth");
114	>	TGraphErrors *gDataWidth = new TGraphErrors();
115	>	gDataWidth->SetTitle("gDataWidth");
116	>	gDataWidth->SetName("gDataWidth");
117	>
118	>	for(int i=0;i<5;i++) {
119	>	NowCut.str("");
120	>	NowCut << "numVtx>=" << NumVtxBin[i] << "&&numVtx<" << NumVtxBin[i+1];
121	>	find_one_peak_combination(TCut(NowCut.str().c_str()),SwitchOffNJetsCut,MCPeak,MCPeakError, DataPeak,DataPeakError,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
122	>	cout << "    " << MCPeak << " +/- " << MCPeakError << endl;
123	>	cout << "    " << DataPeak << " +/- " << DataPeakError << endl;
124	>
125	>	gMCPeak->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),MCPeak);
126	>	gMCPeak->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),MCPeakError);
127	>
128	>	gDataPeak->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),DataPeak);
129	>	gDataPeak->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),DataPeakError);
130	>
131	>	gMCWidth->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),mcSigma);
132	>	gMCWidth->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),mcSigmaError);
133	>
134	>	gDataWidth->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),dataSigma);
135	>	gDataWidth->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),dataSigmaError);
136	>	}
137	>
138	>	can->cd();
139	>	gMCPeak->GetXaxis()->SetTitle("N(Vertices)");
140	>	gMCPeak->GetYaxis()->SetTitle("Peak position");
141	>	gMCPeak->GetXaxis()->CenterTitle();
142	>	gMCPeak->GetYaxis()->CenterTitle();
143	>	gDataPeak->GetXaxis()->SetTitle("N(Vertices)");
144	>	gDataPeak->GetYaxis()->SetTitle("Peak position");
145	>	gDataPeak->GetXaxis()->CenterTitle();
146	>	gDataPeak->GetYaxis()->CenterTitle();
147	>
148	>	gDataWidth->GetXaxis()->SetTitle("N(Vertices)");
149	>	gDataWidth->GetYaxis()->SetTitle("#sigma_{JZB}");
150	>	gDataWidth->GetXaxis()->CenterTitle();
151	>	gDataWidth->GetYaxis()->CenterTitle();
152	>	gMCWidth->GetXaxis()->SetTitle("N(Vertices)");
153	>	gMCWidth->GetYaxis()->SetTitle("#sigma_{JZB}");
154	>	gMCWidth->GetXaxis()->CenterTitle();
155	>	gMCWidth->GetYaxis()->CenterTitle();
156	>
157	>	gMCPeak->SetFillColor(TColor::GetColor("#2E9AFE"));
158	>	gDataPeak->SetFillColor(TColor::GetColor("#2E9AFE"));
159	>	gDataWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
160	>	gMCWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
161	>
162	>	gMCPeak->Draw("A2");
163	>	DrawMCPrelim();
164	>	CompleteSave(can,"PUStudy/MCPeak");
165	>	can->cd();
166	>	gMCWidth->Draw("A2");
167	>	DrawMCPrelim();
168	>	CompleteSave(can,"PUStudy/MCWidth");
169	>	can->cd();
170	>
171	>	gDataPeak->Draw("A2");
172	>	DrawPrelim();
173	>	CompleteSave(can,"PUStudy/DataPeak");
174	>	can->cd();
175	>	gDataWidth->Draw("A2");
176	>	DrawPrelim();
177	>	CompleteSave(can,"PUStudy/DataWidth");
178	>
179	>
180	>
181	>	delete can;
182	>
183	>	}
184	>
185	>	void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb, string sSpecialCut="", bool SwitchOffNJetsCut=false)
186		{
187	+	bool overunderflowstatus=addoverunderflowbins;
188		switch_overunderflow(false);
189	<
189	>
190	>	TCut SpecialCut("mll>=0");
191	>	if(sSpecialCut!="") SpecialCut=TCut(sSpecialCut.c_str());
192	>
193		bool DoInvidualeemmPeaks=false;
194
195		float mcpeak, datapeak;
#	Line 102 | Line 204 | void find_peaks(float &MCPeak,float &MCP
204		float mcSigma,mcSigmaError, dataSigma, dataSigmaError;
205
206		dout << "Finding global peak : " << endl;
207	<	find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
207	>	find_one_peak_combination(SpecialCut,SwitchOffNJetsCut,mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
208
209		if(DoInvidualeemmPeaks) {
210		dout << "Finding peak for electrons : " << endl;
211	<	find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele");
211	>	find_one_peak_combination(SpecialCut&&TCut("id1==0"),SwitchOffNJetsCut,mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele");
212		dout << "Finding peak for muons : " << endl;
213	<	find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu");
213	>	find_one_peak_combination(SpecialCut&&TCut("id1==1"),SwitchOffNJetsCut,mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu");
214
215		datajzb << "(" << jzbvariabledata;
216		mcjzb << "(" << jzbvariablemc;
#	Line 154 | Line 256 | void find_peaks(float &MCPeak,float &MCP
256		MCPeakError=mcpeakerr;
257		DataPeak=datapeak;
258		DataPeakError=datapeakerr;
259	<
158	<	switch_overunderflow(true);
159	<
259	>	switch_overunderflow(overunderflowstatus);
260		}
261
262		void make_special_obs_pred_mll_plot(string datajzb, string mcjzb, float jzbthreshold, float binWidth = 5.0) {
#	Line 223 | Line 323 | void make_special_obs_pred_mll_plot(stri
323
324		RcorrJZBeemm->SetTitleOffset(1.3,"y");
325		RcorrJZBeemm->Draw();
326	<	mcRcorrJZBeemm.Draw("same");
326	>	mcRcorrJZBeemm.Draw("histo,same");
327		Bpred->Draw("histo,same");
328		RcorrJZBeemm->Draw("same");
329
330		TLegend *leg = allsamples.allbglegend();
331		leg->SetX1(0.58);
332	<	leg->AddEntry(RcorrJZBeemm,"observed (data)","l");
332	>	//  leg->AddEntry(RcorrJZBeemm,"observed (data)","l");
333		leg->AddEntry(Bpred,"predicted (data)","l");
334		leg->Draw("same");
335
#	Line 257 | Line 357 | void make_special_obs_pred_mll_plot(stri
357		RcorrJZBeemm->Draw("same");
358
359		TLegend *leg2 = make_legend("",0.52,0.7);
360	<	leg2->AddEntry(RcorrJZBeemm,"observed (data)","lp");
360	>	//  leg2->AddEntry(RcorrJZBeemm,"observed (data)","lp");
361		leg2->AddEntry(Bpred,"predicted (data)","l");
362		leg2->AddEntry(SFN,  "  SF JZB<0","l");
363		leg2->AddEntry(OFN,  "  OF JZB<0","l");
#	Line 280 | Line 380 | void make_special_obs_pred_mll_plot(stri
380		}
381		delete Bpred;
382		delete ckin;
383	+	CleanLegends();
384		}
385
386		void make_special_mll_plot(int nbins, float min, float max, bool logscale,string xlabel) {
#	Line 300 | Line 401 | void make_special_mll_plot(int nbins, fl
401		datahistoOSSF->SetMarkerSize(DataMarkerSize);
402		datahistoOSSF->Draw();
403
404	<	mcstackOSSF.Draw("same");
404	>	mcstackOSSF.Draw("histo,same");
405		datahistoOSSF->Draw("same");
406
407		datahistoOSOF->SetMarkerColor(TColor::GetColor("#FE642E"));
#	Line 326 | Line 427 | void make_special_mll_plot(int nbins, fl
427		delete datahistoOSOF;
428		delete datahistoOSSF;
429		delete ckin;
430	+	CleanLegends();
431		}
432
433
#	Line 367 | Line 469 | void draw_ratio_plot(TH1* hdata, THStack
469		}
470
471		float make_one_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, float ymax, bool logscale,
472	<	string xlabel, string filename, bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) {
472	>	string xlabel, string filename, float legendPosition=0.55) {
473
474		TCut ibasiccut=basiccut;
475		bool draw_separation_lines=false;
#	Line 377 | Line 479 | float make_one_OFSF_plot(string variable
479		TCut cutSF;
480		TCut cutOF;
481
482	+	cout << "Going to use cutOSSF " << (const char*) cutOSSF << " and cutnJets " << (const char*) cutnJets << " and ibasiccut " << (const char*) ibasiccut << endl;
483		cutOF = cutOSOF&&cutnJets&&ibasiccut;
484		cutSF = cutOSSF&&cutnJets&&ibasiccut;
485	+
486	+	TCanvas *ofsf_can = new TCanvas("ofsf_can","ofsf_can");
487
383	–	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
488		TPad* rcan = new TPad("rcan","rcan",0,0,1,1);
489		rcan->SetLogy(logscale);
490		rcan->cd();
491
492	+	std::cout << "OF/SF comparison: variable = " << variable << ", cut = " << cutSF.GetTitle() << std::endl;
493		TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cutSF,data,luminosity);
494		TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cutOF,data,luminosity);
495	<	string signal("LM3");
496	<	TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
497	<	int idx = signalsamples.FindSample(signal)[0];
498	<	(signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF);
499	<	signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
500	<	signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor);
501	<	signalhisto->SetLineStyle(2);
495	>	//   string signal("LM3");
496	>	//   TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
497	>	//   int idx = signalsamples.FindSample(signal)[0];
498	>	//   (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF);
499	>	//   signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
500	>	//   signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor);
501	>	//   signalhisto->SetLineStyle(2);
502		datahistoSF->SetMarkerSize(DataMarkerSize);
503		datahistoOF->SetLineColor(kRed);
504
#	Line 404 | Line 509 | float make_one_OFSF_plot(string variable
509		}
510		if (ymax<0) {
511		if ( logscale ) datahistoSF->SetMaximum(5.3*datahistoSF->GetMaximum());
512	<	else datahistoSF->SetMaximum(1.5*datahistoSF->GetMaximum());
512	>	else datahistoSF->SetMaximum(0.8*datahistoSF->GetMaximum());
513		} else {
514		datahistoSF->SetMaximum(ymax);
515		}
#	Line 423 | Line 528 | float make_one_OFSF_plot(string variable
528		} else {
529		mleg->AddEntry((TObject*)0, "", "");
530		}
531	<	mleg->AddEntry(signalhisto, "LM3", "L");
531	>	//mleg->AddEntry(signalhisto, "LM3", "L");
532
533		datahistoSF->Draw("E1");
534		if (datahistoOF->Integral()>0) datahistoOF->Draw("HIST,SAMES");
535	<	signalhisto->Draw("HIST,SAMES");
535	>	//signalhisto->Draw("HIST,SAMES");
536		mleg->Draw();
537		DrawPrelim();
538	<	if (datahistoOF->Integral()>0) save_with_ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" );
539	<	else CompleteSave(ckin, "SFOF/" + filename);
538	>	if (datahistoOF->Integral()>0) {
539	>	Save_With_Ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" );
540	>	} else {
541	>	CompleteSave(rcan, "SFOF/" + filename);
542	>	}
543
544		datahistoSF->Delete();
545		datahistoOF->Delete();
546	<	signalhisto->Delete();
546	>	//signalhisto->Delete();
547		delete mleg;
548	<	delete ckin;
548	>	delete rcan;
549	>	delete ofsf_can;
550
551		return ymaxSet;
552
#	Line 458 | Line 567 | float make_data_comparison_plot(string v
567
568		TH1F *data1 = allsamples.Draw("data1",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity);
569		TH1F *data2 = comparesamples.Draw("data2",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity);
570	+	data1->Scale(5.0/9.2);
571
572		data1->SetMarkerSize(DataMarkerSize);
573		data2->SetLineColor(kRed);
#	Line 482 | Line 592 | float make_data_comparison_plot(string v
592		data1->GetYaxis()->CenterTitle();
593
594		TLegend *mleg = make_legend("",legendPosition,0.7,false,legendPosition+0.2);
595	<	mleg->AddEntry(data1, "Old 5/fb", "PL");
596	<	mleg->AddEntry(data2, "New x/fb", "L");
595	>	mleg->AddEntry(data1, "2012 data (scaled)", "PL");
596	>	mleg->AddEntry(data2, "2011 data", "L");
597
598		data1->Draw("E1");
599		data2->Draw("HIST,SAMES");
600		mleg->Draw();
601		DrawPrelim();
602	<	save_with_ratio( data1, data2, rcan, "compareData/" + filename, false, false, "old/new" );
602	>	Save_With_Ratio( data1, data2, rcan, "compareData/" + filename, false, false, "new/old" );
603
604		data1->Delete();
605		data2->Delete();
606		delete mleg;
607	+	delete rcan;
608		delete ckin;
609
610		}
#	Line 501 | Line 612 | float make_data_comparison_plot(string v
612		void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale,
613		string xlabel, string filename, float legendPosition=0.55) {
614
615	<	string mllcuts[] = { "mll>20","mll>20&&mll<70", "mll>75&&mll<105", "mll>120" };
616	<	string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "75 < m_{ll} < 105 GeV", "m_{ll} > 120 GeV" };
615	>	string mllcuts[] = { "mll>20","mll>15&&mll<70", "mll>75&&mll<105", "mll>120" };
616	>	string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 120 GeV" };
617		string plotname[] = {"_all","_low","_peak","_high"};
618		float ymax;
619
#	Line 512 | Line 623 | void make_OFSF_plots(string variable, st
623		for ( int i=start; i<4; ++i ) {
624		if ( addcut != "" ) mllcuts[i] += "&&"+addcut;
625		if ( i==start ) {
626	<	ymax = make_one_OFSF_plot(variable, mllcuts[i],            mllcutname[i], nbins, min, max, -1, logscale, xlabel,
626	>	ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel,
627		filename+plotname[i], legendPosition );
628		} else {
629	<	make_one_OFSF_plot(variable, mllcuts[i],            mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
629	>	make_one_OFSF_plot(variable, mllcuts[i],        mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
630		filename+plotname[i], legendPosition );
631		}
632		make_one_OFSF_plot(variable, "id1==1&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel,
#	Line 526 | Line 637 | void make_OFSF_plots(string variable, st
637
638		}
639
640	+	TCut ReplaceInCut(TCut cut, string replacethis, string withthis) {
641	+	string scut=(string)(const char*)cut;
642	+	string acut=ReplaceAll(scut,replacethis,withthis);
643	+	return TCut(acut.c_str());
644	+	}
645	+
646	+	void DoMCSystPlot(string datavariable, string mcvariable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) {
647	+	TCut cut=addcut&&basiccut;
648	+
649	+	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
650	+	ckin->SetLogy(dolog);
651	+	cout << "\r Creating " << filename << " : data   (1/6)" << std::flush;
652	+	TH1F *datahisto = allsamples.Draw("datahisto",datavariable,nbins,min,max, xlabel, "events",cut,data,luminosity);
653	+	datahisto->SetMarkerSize(DataMarkerSize);
654	+	if ( !dolog ) datahisto->SetMinimum(0); // Defaults
655	+	else datahisto->SetMinimum(0.5);
656	+	if (dolog) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
657	+	else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
658	+	cout << "\r Creating " << filename << " : MC central   (2/6)" << std::flush;
659	+	THStack mcstack   = allsamples.DrawStack("mcstack",  mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
660	+	TH1F *MCcentral = CollapseStack(mcstack);
661	+
662	+	TCut bkpcut = cut;
663	+	cut = ReplaceInCut(cut,"pfJetGoodNum40","pfJetGoodNum40p1sigma");
664	+	cout << "\r Creating " << filename << " : MC JES up   (3/6)" << std::flush;
665	+	TH1F *MCJESup   = allsamples.Draw("MCJESup",  mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
666	+
667	+	cut = ReplaceInCut(cut,"pfJetGoodNum40p1sigma","pfJetGoodNum40n1sigma");
668	+	cout << "\r Creating " << filename << " : MC JES down   (4/6)" << std::flush;
669	+	TH1F *MCJESdn   = allsamples.Draw("MCJESdn",  mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
670	+
671	+	cut=bkpcut;
672	+	TCut bkpweight = cutWeight;
673	+	cutWeight= ReplaceInCut(cutWeight,"weight","PUweightUP");
674	+	cout << "\r Creating " << filename << " : MC PU up   (5/6)" << std::flush;
675	+	TH1F *MCPUup    = allsamples.Draw("MCPUup",   mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
676	+
677	+	cutWeight= ReplaceInCut(cutWeight,"PUweightUP","PUweightDOWN");
678	+	cout << "\r Creating " << filename << " : MC PU down   (6/6)" << std::flush;
679	+	TH1F *MCPUdn    = allsamples.Draw("MCPUdn",   mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
680	+	cutWeight=bkpweight;
681	+	cout << "\r Creating " << filename << " : Processing systematics ..." << std::flush;
682	+	TH1F *Systematic = (TH1F*)MCPUdn->Clone("Systematic");
683	+	for(int i=1;i<=Systematic->GetNbinsX();i++) {
684	+	float jesdn = abs(MCcentral->GetBinContent(i)-MCJESdn->GetBinContent(i));
685	+	float jesup = abs(MCcentral->GetBinContent(i)-MCJESup->GetBinContent(i));
686	+	float jes   = jesdn>jesup?jesdn:jesup;
687	+
688	+	float PUup = abs(MCcentral->GetBinContent(i)-MCPUup->GetBinContent(i));
689	+	float PUdn = abs(MCcentral->GetBinContent(i)-MCPUdn->GetBinContent(i));
690	+	float pu   = PUdn>PUup?PUdn:PUup;
691	+
692	+	float sys=sqrt(jes*jes+pu*pu);
693	+	sys/=MCcentral->GetBinContent(i);
694	+	if(MCcentral->GetBinContent(i)==0) sys=0;
695	+	Systematic->SetBinContent(i,sys);
696	+	}
697	+
698	+	datahisto->Draw("e1");
699	+	ckin->Update();
700	+	mcstack.Draw("histo,same");
701	+
702	+	datahisto->Draw("same,e1");
703	+	TLegend *kinleg = allsamples.allbglegend();
704	+	kinleg->Draw();
705
706	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
707	+	kinpad->SetLogy(dolog);
708	+	kinpad->cd();
709	+	datahisto->Draw("e1");
710	+	mcstack.Draw("histo,same");
711	+	datahisto->Draw("same,e1");
712	+	datahisto->Draw("same,axis");
713	+
714	+	kinleg->Draw("same");
715	+	DrawPrelim();
716	+	string saveas="kin/"+filename;
717	+
718	+	cout << "Passing filename to be saved : " << filename << " as " << saveas << endl;
719	+	save_with_ratio_and_sys_band(0.0, datahisto, CollapseStack(mcstack), kinpad->cd(), saveas, false, false, "data/mc",Systematic );
720	+
721	+	ckin->cd();
722	+	MCcentral->SetFillColor(kWhite);
723	+	MCcentral->SetLineColor(kBlack);
724	+	MCcentral->SetLineWidth(2);
725	+	MCPUdn->SetLineColor(kRed);
726	+	MCPUup->SetLineColor(kYellow);
727	+	MCJESdn->SetLineColor(kBlue);
728	+	MCJESup->SetLineColor(kGreen);
729	+
730	+	TLegend *legsmall = new TLegend(0.0,0.85,1.0,1.00);
731	+	legsmall->SetNColumns(5);
732	+	legsmall->AddEntry(MCcentral,"central","L");
733	+	legsmall->AddEntry(MCPUdn,"PU down","L");
734	+	legsmall->AddEntry(MCPUup,"PU up","L");
735	+	legsmall->AddEntry(MCJESdn,"JES down","L");
736	+	legsmall->AddEntry(MCJESup,"JES up","L");
737	+	legsmall->SetFillColor(kWhite);
738	+	legsmall->SetBorderSize(0);
739	+
740	+	datahisto->Draw("e1");
741	+	MCcentral->Draw("histo,same");
742	+	MCPUdn->Draw("histo,same");
743	+	MCPUup->Draw("histo,same");
744	+	MCJESdn->Draw("histo,same");
745	+	MCJESup->Draw("histo,same");
746	+	datahisto->Draw("e1,same");
747	+	legsmall->Draw();
748	+	CompleteSave(ckin,saveas+"___AllLines");
749	+
750	+	delete datahisto;
751	+	delete MCcentral;
752	+	delete MCPUdn;
753	+	delete MCPUup;
754	+	delete MCJESdn;
755	+	delete MCJESup;
756	+	delete ckin;
757	+	}
758	+
759	+	void DoMCSystPlot(string variable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) {
760	+	DoMCSystPlot(variable, variable, nbins, min, max, addcut, dolog, xlabel, filename);
761	+	}
762	+
763		float lastrange_min=0;
764		float lastrange_max=0;
765
#	Line 535 | Line 768 | void make_kin_plot(string variable, stri
768		//  TCut basiccut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])&&(passed_triggers||!is_data)");
769		TCut ibasiccut=basiccut;
770		bool draw_separation_lines=false;
771	+	bool drawsignal = false;
772
773		if(isPF) ibasiccut=basiccut&&"pfjzb[0]>-998";
774
#	Line 556 | Line 790 | void make_kin_plot(string variable, stri
790		//Step 2: Refine the cut
791		TCut cut;
792		cut=cutmass&&cutOSSF&&cutnJets&&ibasiccut;
793	<	if(filename=="nJets" || filename=="nJets_inclusive") cut=cutmass&&cutOSSF&&ibasiccut;
794	<	if(filename=="nJets_osof" || filename=="nJets_osof_inclusive") cut=cutmass&&cutOSOF&&ibasiccut;
793	>	if(filename=="nJets" || filename=="nJets_inclusive" || filename=="nJets_met100" || filename=="nJets_inclusive_met100") cut=cutmass&&cutOSSF&&ibasiccut;
794	>	if(filename=="nJets_osof" || filename=="nJets_osof_inclusive" || filename=="nJets_osof_met100" || filename=="nJets_osof_inclusive_met100") cut=cutmass&&cutOSOF&&ibasiccut;
795		if(filename=="nJets_nocuts_except_mll_ossf") cut=cutmass&&cutOSSF;
796	<	if(filename=="mll") {
796	>	if(filename=="mll"||filename=="mll_met100") {
797		cut=cutOSSF&&cutnJets&&ibasiccut;
798		draw_separation_lines=true;
799		}
800	<	if(filename=="mll_ee") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0";
801	<	if(filename=="mll_mm") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1";
802	<	if(filename=="mll_osof") {
800	>	if(filename=="mll_ee"||filename=="mll_ee_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0";
801	>	if(filename=="mll_mm"||filename=="mll_mm_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1";
802	>	if(filename=="mll_osof"||filename=="mll_osof_met100") {
803		cut=cutOSOF&&cutnJets&&ibasiccut;
804		draw_separation_lines=true;
805		}
#	Line 573 | Line 807 | void make_kin_plot(string variable, stri
807		if(Contains(filename,"mll_ee_above")) cut=cut&&"id1==0";
808		if(Contains(filename,"mll_mm_above")) cut=cut&&"id1==1";
809		if(Contains(filename,"mll_osof_aboveJZB")) cut=cutOSOF&&cutnJets&&ibasiccut;
810	<	if(filename=="mll_inclusive"||filename=="mll_inclusive_highrange") cut=cutmass&&cutOSSF;
811	<	if(filename=="mll_inclusive_osof") cut=cutmass&&cutOSOF;
812	<	if(filename=="mll_inclusive_ee") cut=cutmass&&cutOSSF&&"id1==0";
813	<	if(filename=="mll_inclusive_mm") cut=cutmass&&cutOSSF&&"id1==1";
810	>	if(filename=="mll_inclusive"||filename=="mll_inclusive_highrange") cut=cutOSSF;
811	>	if(filename=="mll_inclusive_osof") cut=cutOSOF;
812	>	if(filename=="mll_inclusive_ee") cut=cutOSSF&&"id1==0";
813	>	if(filename=="mll_inclusive_mm") cut=cutOSSF&&"id1==1";
814		if(filename=="pfJetGoodEta_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets;
815		if(filename=="pfJetGoodPt_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets;
816		if(filename=="numVtx") cut=cutmass&&ibasiccut;
#	Line 597 | Line 831 | void make_kin_plot(string variable, stri
831		if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
832		else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
833
834	<	cout << "******** Cut used : " << (const char*) cut << endl;
834	>	cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl;
835		if(loadlastminmax) {
836		datahisto->SetMinimum(lastrange_min);
837		datahisto->SetMaximum(lastrange_max);
#	Line 611 | Line 845 | void make_kin_plot(string variable, stri
845		string signal("LM3");
846		TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
847		int idx = signalsamples.FindSample(signal)[0];
848	<	(signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut);
849	<	signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
850	<	signalhisto->SetLineColor(kOrange);
848	>	if(drawsignal) (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut);
849	>	if(drawsignal) signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
850	>	if(drawsignal) signalhisto->SetLineColor(kOrange);
851
852		THStack mcstack   = allsamples.DrawStack("mcstack",  variable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
853		datahisto->Draw("e1");
854		ckin->Update();
855	<	mcstack.Draw("same");
855	>	mcstack.Draw("histo,same");
856
857		datahisto->Draw("same,e1");
858		TLegend *kinleg = allsamples.allbglegend();
#	Line 661 | Line 895 | void make_kin_plot(string variable, stri
895		kinpad->cd();
896		kinpad->SetLogy(logscale);
897		datahisto->Draw("e1");
898	<	mcstack.Draw("same");
899	<	signalhisto->Draw("same");
898	>	mcstack.Draw("histo,same");
899	>	if(drawsignal) signalhisto->Draw("same");
900		datahisto->Draw("same,e1");
901		datahisto->Draw("same,axis");
902		if(RestrictToMassPeak&&draw_separation_lines) {
#	Line 670 | Line 904 | void make_kin_plot(string variable, stri
904		upperboundary->Draw("same");
905		}
906
907	<	kinleg->AddEntry("signalhisto",signal.c_str(),"l");
907	>	if(drawsignal) kinleg->AddEntry("signalhisto",signal.c_str(),"l");
908		kinleg->Draw();
909		write_cut->Draw();
910		DrawPrelim();
911		string saveas="kin/"+filename;
912		if(isPF) saveas="kin/"+filename+"__PF";
913	<	save_with_ratio(datahisto,mcstack,kinpad->cd(),saveas);
913	>	Save_With_Ratio(datahisto,mcstack,kinpad->cd(),saveas);
914		//      if(isPF) CompleteSave(with_ratio,"kin/"+filename+"__PF_withratio");
915		//      else CompleteSave(with_ratio,"kin/"+filename+"_withratio");
916		//      delete with_ratio;
#	Line 685 | Line 919 | void make_kin_plot(string variable, stri
919		else CompleteSave(ckin,"kin/"+filename);
920		}
921		datahisto->Delete();
922	+	delete signalhisto;
923	+	delete ckin;
924	+	}
925	+
926	+
927	+	void make_plain_kin_plot(string variable, string addcut, int nbins, float min, float max, bool logscale,
928	+	string xlabel, string filename, bool isPF=true, bool plotratio=true, bool loadlastminmax=false ) {
929	+	TCut cut=TCut(addcut.c_str())&&basiccut;
930	+
931	+	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
932	+	ckin->SetLogy(logscale);
933	+	TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity);
934	+	datahisto->SetMarkerSize(DataMarkerSize);
935	+	if ( !logscale ) datahisto->SetMinimum(0); // Defaults
936	+	else datahisto->SetMinimum(0.5);
937	+	// Custom max.
938	+	if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
939	+	else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
940	+
941	+	cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl;
942	+	if(loadlastminmax) {
943	+	datahisto->SetMinimum(lastrange_min);
944	+	datahisto->SetMaximum(lastrange_max);
945	+	if(logscale) {
946	+	datahisto->SetMinimum(pow(10,lastrange_min));
947	+	datahisto->SetMaximum(pow(10,lastrange_max));
948	+	}
949	+	}
950	+
951	+	THStack mcstack   = allsamples.DrawStack("mcstack",  variable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
952	+	datahisto->Draw("e1");
953	+	ckin->Update();
954	+	mcstack.Draw("histo,same");
955	+
956	+	datahisto->Draw("same,e1");
957	+	TLegend *kinleg = allsamples.allbglegend();
958	+	kinleg->Draw();
959	+
960	+	lastrange_min=ckin->GetUymin();
961	+	lastrange_max=ckin->GetUymax();
962	+
963	+	if ( plotratio ) {
964	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
965	+	kinpad->cd();
966	+	kinpad->SetLogy(logscale);
967	+	datahisto->Draw("e1");
968	+	mcstack.Draw("histo,same");
969	+	datahisto->Draw("same,e1");
970	+	datahisto->Draw("same,axis");
971	+
972	+	kinleg->Draw();
973	+	DrawPrelim();
974	+	string saveas="kin/"+filename;
975	+	if(isPF) saveas="kin/"+filename+"__PF";
976	+	Save_With_Ratio(datahisto,mcstack,kinpad->cd(),saveas);
977	+	} else {
978	+	if(isPF) CompleteSave(ckin,"kin/"+filename+"__PF");
979	+	else CompleteSave(ckin,"kin/"+filename);
980	+	}
981	+	datahisto->Delete();
982		delete ckin;
983		}
984
985	+
986		void make_JES_plot(TCut cut, string name) {
987
988		int nbins=10;
#	Line 721 | Line 1016 | void make_JES_plot(TCut cut, string name
1016		JESunc->SetFillColor(TColor::GetColor("#00ADE1"));
1017		JESunc->SetFillStyle(3002);
1018		datahisto->Draw("e1");
1019	<	mcstack.Draw("same");
1019	>	mcstack.Draw("histo,same");
1020		JESunc->Draw("2");
1021		datahisto->Draw("same,e1");
1022		TLegend *kinleg = allsamples.allbglegend();
#	Line 730 | Line 1025 | void make_JES_plot(TCut cut, string name
1025		CompleteSave(ckin,"Systematics/JES"+name);
1026		datahisto->Delete();
1027		delete ckin;
1028	+	delete JESunc;
1029	+	delete JESup;
1030	+	delete JESdn;
1031	+	CleanLegends();
1032	+	}
1033	+
1034	+	string Cut2Str(TCut cut) {
1035	+	return ((string)(const char*)cut);
1036	+	}
1037	+
1038	+	void MakeElegantTwoThreeComparisons() {
1039	+
1040	+	TCut signal("met[4]>100&&met[4]<150&&mll>20"&&PlottingSetup::basicqualitycut);
1041	+	TCut TwoJets("pfJetGoodNum40==2");
1042	+	TCut ThreeJets("pfJetGoodNum40>=3");
1043	+
1044	+	float min=20;
1045	+	float max=300;
1046	+	int nbins=int((max-min)/5.0);
1047	+	string xlabel="m_{ll} [GeV]";
1048	+
1049	+	TCanvas *can2 = new TCanvas("can2","can2");
1050	+	TH1F *TwoOF = allsamples.Draw("TwoOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_"));
1051	+	TH1F *TwoSF = allsamples.Draw("TwoSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_"));
1052	+	TH1F *ThreeOF = allsamples.Draw("ThreeOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_"));
1053	+	TH1F *ThreeSF = allsamples.Draw("ThreeSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_"));
1054	+
1055	+	can2->cd(1);
1056	+
1057	+	TwoSF->SetFillColor(kWhite);
1058	+	TwoOF->SetFillColor(kWhite);
1059	+	ThreeSF->SetFillColor(kWhite);
1060	+	ThreeOF->SetFillColor(kWhite);
1061	+
1062	+	TPad *kin3 = new TPad("kin3","kin3",0,0,1,1);
1063	+	kin3->cd();
1064	+	TwoOF->Draw("histo");
1065	+	TwoSF->Draw("e1,same");
1066	+	TLegend *leg_SFOF_Two = make_legend();
1067	+	leg_SFOF_Two->SetHeader("t#bar{t} MC:");
1068	+	leg_SFOF_Two->AddEntry(TwoOF,"OF, nJets==2","l");
1069	+	leg_SFOF_Two->AddEntry(TwoSF,"SF, nJets==2","lp");
1070	+	leg_SFOF_Two->Draw();
1071	+	Save_With_Ratio(TwoSF,TwoOF,kin3->cd(),"DoubleRatios/SFOF_Two");
1072	+
1073	+	TPad *kin4 = new TPad("kin4","kin4",0,0,1,1);
1074	+	kin4->cd();
1075	+	ThreeOF->Draw("histo");
1076	+	ThreeSF->Draw("e1,same");
1077	+	TLegend *leg_SFOF_Three = make_legend();
1078	+	leg_SFOF_Three->SetHeader("t#bar{t} MC:");
1079	+	leg_SFOF_Three->AddEntry(TwoOF,"OF, nJet#geq3","l");
1080	+	leg_SFOF_Three->AddEntry(TwoSF,"SF, nJets#geq3","lp");
1081	+	leg_SFOF_Three->Draw();
1082	+	Save_With_Ratio(ThreeSF,ThreeOF,kin4->cd(),"DoubleRatios/SFOF_Three");
1083	+
1084	+
1085	+	TwoSF->Scale(ThreeSF->Integral()/TwoSF->Integral());
1086	+	TPad *kin = new TPad("kin","kin",0,0,1,1);
1087	+	kin->cd();
1088	+	TwoSF->Draw("histo");
1089	+	ThreeSF->Draw("e1,same");
1090	+	TLegend *leg_SF_Two_Three = make_legend();
1091	+	leg_SF_Two_Three->SetHeader("t#bar{t} MC:");
1092	+	leg_SF_Two_Three->AddEntry(TwoSF,"#splitline{SF, nJets==2}{ (scaled)}","l");
1093	+	leg_SF_Two_Three->AddEntry(ThreeSF,"SF, nJets#geq3","lp");
1094	+	leg_SF_Two_Three->Draw();
1095	+	Save_With_Ratio(TwoSF,ThreeSF,kin->cd(),"DoubleRatios/SF_Two_Three");
1096	+
1097	+	TwoOF->Scale(ThreeOF->Integral()/TwoOF->Integral());
1098	+	TPad *kin2 = new TPad("kin2","kin2",0,0,1,1);
1099	+	kin2->cd();
1100	+	TwoOF->Draw("histo");
1101	+	ThreeOF->Draw("e1,same");
1102	+	TLegend *leg_OF_Two_Three = make_legend();
1103	+	leg_OF_Two_Three->SetHeader("t#bar{t} MC:");
1104	+	leg_OF_Two_Three->AddEntry(TwoSF,"#splitline{OF, nJets==2}{ (scaled)}","l");
1105	+	leg_OF_Two_Three->AddEntry(ThreeSF,"OF, nJets#geq3","lp");
1106	+	leg_OF_Two_Three->Draw();
1107	+	Save_With_Ratio(TwoOF,ThreeOF,kin2->cd(),"DoubleRatios/OF_Two_Three");
1108	+
1109	+	}
1110	+
1111	+	void PresentRange(TH1F *Down, TH1F *Up, TH1F *central, TVirtualPad *pad, string title) {
1112	+	pad->cd();
1113	+	TGraphErrors *gr = new TGraphErrors(Down->GetNbinsX());
1114	+	for(int i=1;i<=Down->GetNbinsX();i++) {
1115	+	float average=0.5*(Down->GetBinContent(i)+Up->GetBinContent(i));
1116	+	float error=abs(Down->GetBinContent(i)-average);
1117	+	gr->SetPoint(i-1,Down->GetBinCenter(i),average);
1118	+	gr->SetPointError(i-1,Down->GetBinWidth(i)/2.0,error);
1119	+	}
1120	+	gr->SetFillColor(TColor::GetColor("#F79F81"));
1121	+	gr->GetYaxis()->SetRangeUser(0,2);
1122	+	gr->GetXaxis()->SetTitle(central->GetXaxis()->GetTitle());
1123	+	gr->GetYaxis()->SetTitle(central->GetYaxis()->GetTitle());
1124	+	gr->GetXaxis()->CenterTitle();
1125	+	gr->GetYaxis()->CenterTitle();
1126	+
1127	+	gr->Draw("A2");
1128	+	central->Draw("e1,same");
1129	+	DrawMCPrelim();
1130	+	TText *rtitle = write_text(0.8,0.8,title);
1131	+	TLine *line = new TLine(Down->GetBinLowEdge(1),1.0,Down->GetBinLowEdge(Down->GetNbinsX())+Down->GetBinWidth(Down->GetNbinsX()),1.0);
1132	+	line->SetLineColor(kBlue);
1133	+	line->SetLineStyle(2);
1134	+	line->Draw("same");
1135	+	rtitle->Draw("same");
1136	+	}
1137	+
1138	+	TH1F* ProduceTTbarRatio(TH1F *datahisto, TH1F *fullmc, TH1F *originalttbar, TH1F *httbar) {
1139	+	TH1F *basehisto = (TH1F*)fullmc->Clone("basehisto");
1140	+	basehisto->Add(originalttbar,-1);
1141	+	basehisto->Add(httbar);
1142	+
1143	+	string sname=httbar->GetName();
1144	+	sname="R"+sname.substr(1,sname.length());
1145	+	TH1F *ratio = (TH1F*)datahisto->Clone(sname.c_str());
1146	+	ratio->Divide(basehisto);
1147	+	delete basehisto;
1148	+	return ratio;
1149	+	}
1150	+
1151	+	void make_ttbar_comparison(string datavariable, string mcvariable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) {
1152	+
1153	+	TCut cut=addcut&&basiccut;
1154	+
1155	+	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
1156	+	ckin->SetLogy(dolog);
1157	+	cout << "\r Creating " << filename << " : data   (1/7)" << std::flush;
1158	+	TH1F *datahisto = allsamples.Draw("datahisto",datavariable,nbins,min,max, xlabel, "events",cut,data,luminosity);
1159	+	datahisto->SetMarkerSize(DataMarkerSize);
1160	+	if ( !dolog ) datahisto->SetMinimum(0); // Defaults
1161	+	else datahisto->SetMinimum(0.5);
1162	+	if (dolog) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
1163	+	else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
1164	+	cout << "\r Creating " << filename << " : MC central  (stack) (2/7)" << std::flush;
1165	+	THStack mcstack   = allsamples.DrawStack("mcstack",  mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
1166	+	TH1F *MCcentral = CollapseStack(mcstack);
1167	+
1168	+	cout << "\r Creating " << filename << " : MC central (histo)   (3/7)" << std::flush;
1169	+	TH1F *TCentral        = allsamples.Draw ("TCentral",     mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,allsamples.FindSample("TTJets_"));
1170	+	TH1F *RCentral        = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TCentral);
1171	+
1172	+	cout << "\r Creating " << filename << " : MC Matching up   (4/7)" << std::flush;
1173	+	TH1F *TMatchingUp     = systsamples.Draw("TMatchingUp",  mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_matchingup"));
1174	+	TH1F *RMatchingUp     = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TMatchingUp);
1175	+
1176	+	cout << "\r Creating " << filename << " : MC Matching down   (5/7)" << std::flush;
1177	+	TH1F *TMatchingDown   = systsamples.Draw("TMatchingDown",mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_matchingdown"));
1178	+	TH1F *RMatchingDown   = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TMatchingDown);
1179	+
1180	+	cout << "\r Creating " << filename << " : MC Scale up   (6/7)" << std::flush;
1181	+	TH1F *TScaleUp        = systsamples.Draw("TScaleUp",     mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_scaleup"));
1182	+	TH1F *RScaleUp        = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TScaleUp);
1183	+
1184	+	cout << "\r Creating " << filename << " : MC Scale down   (7/7)" << std::flush;
1185	+	TH1F *TScaleDown      = systsamples.Draw("TScaleDown",   mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_scaledown"));
1186	+	TH1F *RScaleDown      = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TScaleDown);
1187	+
1188	+	datahisto->Draw("e1");
1189	+	ckin->Update();
1190	+	mcstack.Draw("histo,same");
1191	+
1192	+	datahisto->Draw("same,e1");
1193	+	TLegend *kinleg = allsamples.allbglegend();
1194	+	kinleg->Draw();
1195	+
1196	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
1197	+	kinpad->SetLogy(dolog);
1198	+	kinpad->cd();
1199	+	datahisto->Draw("e1");
1200	+	mcstack.Draw("histo,same");
1201	+	datahisto->Draw("same,e1");
1202	+	datahisto->Draw("same,axis");
1203	+
1204	+	kinleg->Draw("same");
1205	+	DrawPrelim();
1206	+	string saveas="kin/"+filename;
1207	+
1208	+	cout << "Passing filename to be saved : " << filename << " as " << saveas << endl;
1209	+	Save_With_Ratio( datahisto, CollapseStack(mcstack), kinpad->cd(), saveas, false, false, "data/mc");
1210	+	delete kinpad;
1211	+
1212	+	ckin->cd();
1213	+	TCentral->SetLineColor(kBlack);
1214	+
1215	+	TMatchingUp->SetLineColor(TColor::GetColor("#2c17b1"));//blue
1216	+	TMatchingUp->SetLineStyle(2);//dashed
1217	+	TMatchingDown->SetLineColor(TColor::GetColor("#2c17b1"));
1218	+	TMatchingDown->SetLineStyle(3);//dotted
1219	+
1220	+	TScaleUp->SetLineColor(TColor::GetColor("#FF8500"));//orange
1221	+	TScaleUp->SetLineStyle(2);//dashed
1222	+	TScaleDown->SetLineColor(TColor::GetColor("#FF8500"));
1223	+	TScaleDown->SetLineStyle(3);//dotted
1224	+
1225	+	MCcentral->SetLineStyle(3);
1226	+	MCcentral->SetLineColor(kRed);
1227	+
1228	+	TCentral->Draw("e1,same");
1229	+	TCentral->Draw("histo,same");
1230	+	TMatchingUp->Draw("histo,same");
1231	+	TMatchingDown->Draw("histo,same");
1232	+	TScaleUp->Draw("histo,same");
1233	+	TScaleDown->Draw("histo,same");
1234	+
1235	+	TLegend *sleg = make_legend();
1236	+	sleg->AddEntry(TCentral,"Central","L");
1237	+	sleg->AddEntry(TMatchingUp,"matching up","L");
1238	+	sleg->AddEntry(TMatchingDown,"matching down","L");
1239	+	sleg->AddEntry(TScaleUp,"scale up","L");
1240	+	sleg->AddEntry(TScaleDown,"scale down","L");
1241	+	sleg->Draw();
1242	+
1243	+	CompleteSave(ckin,saveas+"___AllLines");
1244	+	ckin->cd();
1245	+	ckin->SetLogy(0);
1246	+
1247	+	RCentral->SetLineColor(kBlack);
1248	+
1249	+	RMatchingUp->SetLineColor(TColor::GetColor("#2c17b1"));//blue
1250	+	RMatchingUp->SetLineStyle(2);//dashed
1251	+	RMatchingDown->SetLineColor(TColor::GetColor("#2c17b1"));
1252	+	RMatchingDown->SetLineStyle(3);//dotted
1253	+
1254	+	RScaleUp->SetLineColor(TColor::GetColor("#FF8500"));//orange
1255	+	RScaleUp->SetLineStyle(2);//dashed
1256	+	RScaleDown->SetLineColor(TColor::GetColor("#FF8500"));
1257	+	RScaleDown->SetLineStyle(3);//dotted
1258	+
1259	+	RCentral->GetYaxis()->SetRangeUser(0,2);
1260	+	RCentral->Draw("e1");
1261	+	RMatchingUp->Draw("histo,same");
1262	+	RMatchingDown->Draw("histo,same");
1263	+	RScaleUp->Draw("histo,same");
1264	+	RScaleDown->Draw("histo,same");
1265	+
1266	+	CompleteSave(ckin,saveas+"___AllRatios");
1267	+
1268	+	TCanvas *multicanvas = new TCanvas("multicanvas","multicanvas",1400,700);
1269	+	multicanvas->Divide(2,1);
1270	+	PresentRange(RMatchingUp,RMatchingDown,RCentral,multicanvas->cd(1),"Matching");
1271	+	PresentRange(RScaleUp,RScaleDown,RCentral,multicanvas->cd(2),"Scale");
1272	+	CompleteSave(multicanvas,saveas+"___RangeIllustration");
1273	+
1274	+
1275	+	delete datahisto;
1276	+	delete   TCentral;
1277	+	delete   TMatchingUp;
1278	+	delete   TMatchingDown;
1279	+	delete   TScaleUp;
1280	+	delete   TScaleDown;
1281	+	//   delete   Ttmass166;
1282	+	//   delete   Ttmass169;
1283	+	//   delete   Ttmass175;
1284	+	//   delete   Ttmass178;
1285	+
1286	+	delete ckin;
1287	+
1288	+	}
1289	+
1290	+	void make_ttbar_comparison(string variable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) {
1291	+	make_ttbar_comparison(variable, variable, nbins, min, max, addcut, dolog, xlabel, filename);
1292	+	}
1293	+
1294	+	void ProduceJanPlots() {
1295	+	TCanvas *c1 = new TCanvas("c1","c1");
1296	+	TH1F *metSF = allsamples.Draw("metSF","met[4]",30,100,400, "E_{T}^{miss} [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data, luminosity);
1297	+	TH1F *metOF = allsamples.Draw("metOF","met[4]",30,100,400, "E_{T}^{miss} [GeV]", "events", cutmass&&cutOSOF&&cutnJets,data, luminosity);
1298	+
1299	+	TPad* rcan = new TPad("rcan","rcan",0,0,1,1);
1300	+	rcan->cd();
1301	+
1302	+	metOF->SetLineColor(kRed);
1303	+	metSF->Draw("e1");
1304	+	metOF->Draw("histo,same");
1305	+	metSF->Draw("e1,same");
1306	+
1307	+	TLegend *leg = make_legend();
1308	+
1309	+	leg->AddEntry(metSF,"Data SF","p");
1310	+	leg->AddEntry(metOF,"Data OF","l");
1311	+	leg->Draw();
1312	+
1313	+	DrawPrelim();
1314	+
1315	+	Save_With_Ratio(metSF,metOF,rcan,"JanPlots/ETHConfig",false,false,"SF/OF");
1316	+
1317	+	delete rcan;
1318	+	delete metSF;
1319	+	delete metOF;
1320	+	delete c1;
1321	+	}
1322	+
1323	+	THStack MakeOneSystematicsPlot(TCut cut, string saveas, string variation, TH1F *hdata, string variable, int nbins, float bmin, float bmax, string label, TH1F* &thisto) {
1324	+	THStack SystPlot = allsamples.DrawStack(variation,variable,nbins,bmin,bmax,label,"events",cut,mc,PlottingSetup::luminosity);
1325	+
1326	+	//now need to process the plot (calculate the variation and set the member of thstack accordingly!)
1327	+	if(variation!="Central") {
1328	+	TH1F *varttbar;
1329	+	if(variation=="MatchingUp") varttbar = systsamples.Draw("varttbar",variable,nbins,bmin,bmax,label,"events",cut,mc,PlottingSetup::luminosity,systsamples.FindSample("TTJets_matchingup_TuneZ2s"));
1330	+	if(variation=="MatchingDown") varttbar = systsamples.Draw("varttbar",variable,nbins,bmin,bmax,label,"events",cut,mc,PlottingSetup::luminosity,systsamples.FindSample("TTJets_matchingdown"));
1331	+	if(variation=="ScaleUp") varttbar = systsamples.Draw("varttbar",variable,nbins,bmin,bmax,label,"events",cut,mc,PlottingSetup::luminosity,systsamples.FindSample("TTJets_scaleup"));
1332	+	if(variation=="ScaleDown") varttbar = systsamples.Draw("varttbar",variable,nbins,bmin,bmax,label,"events",cut,mc,PlottingSetup::luminosity,systsamples.FindSample("TTJets_scaledown"));
1333	+	assert(varttbar);
1334	+
1335	+	TIter nextHisto(SystPlot.GetHists());
1336	+	TH1F* h;
1337	+	while ( h = (TH1F*)nextHisto() ) {
1338	+	if(Contains(h->GetName(),"t_bar_t")) {
1339	+	varttbar->Scale(h->Integral()/varttbar->Integral());
1340	+	for(int i=0;i<=h->GetNbinsX()+1;i++) {//note that we deliberatly consider the under/overflow bin as well!
1341	+	h->SetBinContent(i,varttbar->GetBinContent(i));
1342	+	h->SetBinError(i,varttbar->GetBinError(i));
1343	+	}
1344	+	thisto=(TH1F*)varttbar->Clone(variation.c_str());
1345	+	SystPlot.Modified();
1346	+	}
1347	+	}
1348	+	delete varttbar;
1349	+	} else {
1350	+	TIter nextHisto(SystPlot.GetHists());
1351	+	TH1F* h;
1352	+	while ( h = (TH1F*)nextHisto() ) {
1353	+	if(Contains(h->GetName(),"t_bar_t")) thisto=(TH1F*)h->Clone(variation.c_str());
1354	+	}
1355	+	}
1356	+
1357	+	TLegend *fullleg = allsamples.allbglegend();
1358	+	fullleg->SetHeader(variation.c_str());
1359	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
1360	+	kinpad->SetLogy(1);
1361	+	kinpad->cd();
1362	+
1363	+	hdata->Draw("e1");
1364	+	SystPlot.Draw("histo,same");
1365	+	hdata->Draw("e1,same");
1366	+	fullleg->Draw();
1367	+	DrawPrelim();
1368	+
1369	+	Save_With_Ratio(hdata,SystPlot,kinpad,saveas+"_"+variation);
1370	+	CleanLegends();
1371	+
1372	+	delete kinpad;
1373	+	return SystPlot;
1374	+	}
1375	+
1376	+
1377	+
1378	+	void ProduceMCSystematicPlot(string variable, int nbins, float bmin, float bmax, string label, TCut thiscut, string saveas) {
1379	+	TCanvas *ca = new TCanvas("ca","ca");
1380	+	TH1F *hdata = allsamples.Draw("hdata",variable,nbins,bmin,bmax,label,"events",thiscut,data,luminosity);
1381	+
1382	+	TH1F *tScaleUp,*tScaleDown,*tMatchingUp,*tMatchingDown,*tCentral;
1383	+
1384	+	THStack ScaleUp      = MakeOneSystematicsPlot(thiscut,saveas,"ScaleUp",hdata,variable, nbins, bmin, bmax, label,tScaleUp);
1385	+	THStack ScaleDown    = MakeOneSystematicsPlot(thiscut,saveas,"ScaleDown",hdata,variable, nbins, bmin, bmax, label,tScaleDown);
1386	+	THStack MatchingUp   = MakeOneSystematicsPlot(thiscut,saveas,"MatchingUp",hdata,variable, nbins, bmin, bmax, label,tMatchingUp);
1387	+	THStack MatchingDown = MakeOneSystematicsPlot(thiscut,saveas,"MatchingDown",hdata,variable, nbins, bmin, bmax, label,tMatchingDown);
1388	+
1389	+	TH1F *HScaleUp = CollapseStack(ScaleUp);
1390	+	TH1F *HScaleDown = CollapseStack(ScaleDown);
1391	+	TH1F *HMatchingUp = CollapseStack(MatchingUp);
1392	+	TH1F *HMatchingDown = CollapseStack(MatchingDown);
1393	+
1394	+	THStack Central = MakeOneSystematicsPlot(thiscut,saveas,"Central",hdata,variable, nbins, bmin, bmax, label,tCentral);
1395	+	TH1F *HCentral = CollapseStack(Central);
1396	+
1397	+	TH1F *Systematic = (TH1F*)hdata->Clone("Systematic");
1398	+	for(int i=1;i<=Systematic->GetNbinsX();i++) {
1399	+	float matchingup   = abs(HMatchingUp->GetBinContent(i)-HCentral->GetBinContent(i));
1400	+	float matchingdown = abs(HMatchingDown->GetBinContent(i)-HCentral->GetBinContent(i));
1401	+	float scaleup   = abs(HScaleUp->GetBinContent(i)-HCentral->GetBinContent(i));
1402	+	float scaledown = abs(HScaleDown->GetBinContent(i)-HCentral->GetBinContent(i));
1403	+
1404	+	float matching = matchingup>matchingdown?matchingup:matchingdown;
1405	+	float scale    = scaleup>scaledown?scaleup:scaledown;
1406	+
1407	+	float sys=sqrt(matching*matching+scale*scale);
1408	+	sys/=HCentral->GetBinContent(i);
1409	+	if(HCentral->GetBinContent(i)==0) sys=0;
1410	+	Systematic->SetBinContent(i,sys);
1411	+	}
1412	+
1413	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
1414	+	kinpad->SetLogy(1);
1415	+	kinpad->cd();
1416	+	if(variable=="genMET") {
1417	+	Central.Draw("histo");
1418	+	} else {
1419	+	hdata->Draw("e1");
1420	+	Central.Draw("histo,same");
1421	+	hdata->Draw("e1,same");
1422	+	}
1423	+
1424	+	DrawPrelim();
1425	+
1426	+	save_with_ratio_and_sys_band( 0.0, hdata, HCentral, kinpad->cd(), saveas, false, false, "data/mc",Systematic );
1427	+
1428	+	kinpad->cd();
1429	+	tCentral->SetFillColor(kWhite);
1430	+	tCentral->SetLineColor(kBlack);
1431	+	tScaleUp->SetLineColor(kRed);
1432	+	tScaleDown->SetLineColor(kRed);
1433	+	tMatchingUp->SetLineColor(kBlue);
1434	+	tMatchingDown->SetLineColor(kBlue);
1435	+	tScaleUp->SetLineStyle(2);
1436	+	tScaleDown->SetLineStyle(3);
1437	+	tMatchingUp->SetLineStyle(2);
1438	+	tMatchingDown->SetLineStyle(3);
1439	+
1440	+	TLegend *leg2 = make_legend();
1441	+	leg2->AddEntry(hdata,"Data","p");
1442	+	leg2->AddEntry(tCentral,"Central (ttbar)","l");
1443	+	leg2->AddEntry(tScaleUp,"ScaleUp (ttbar)","l");
1444	+	leg2->AddEntry(tScaleDown,"ScaleDown (ttbar)","l");
1445	+	leg2->AddEntry(tMatchingUp,"MatchingUp (ttbar)","l");
1446	+	leg2->AddEntry(tMatchingDown,"MatchingDown (ttbar)","l");
1447	+	if(variable=="genMET") {
1448	+	Central.Draw("histo");
1449	+	}else {
1450	+	hdata->Draw("e1");
1451	+	Central.Draw("histo,same");
1452	+	hdata->Draw("e1,same");
1453	+	}
1454	+	tCentral->Draw("histo,same");
1455	+	tScaleUp->Draw("histo,same");
1456	+	tScaleDown->Draw("histo,same");
1457	+	tMatchingUp->Draw("histo,same");
1458	+	tMatchingDown->Draw("histo,same");
1459	+	leg2->Draw();
1460	+
1461	+	CompleteSave(kinpad,saveas+"__TTbarComparison");
1462	+
1463	+	gStyle->SetOptFit(0);
1464	+
1465	+	kinpad->cd();
1466	+	kinpad->SetLogy(0);
1467	+	TH1F *MatchingRatio = (TH1F*)tMatchingUp->Clone("MatchingRatio");
1468	+	MatchingRatio->Divide(tMatchingDown);
1469	+	TLine *lone = new TLine(tScaleUp->GetBinLowEdge(1),1,tScaleUp->GetBinLowEdge(tScaleUp->GetNbinsX())+tScaleUp->GetBinWidth(tScaleUp->GetNbinsX()),1);
1470	+	lone->SetLineColor(TColor::GetColor("#01DF01"));
1471	+	lone->SetLineStyle(2);
1472	+	TH1F *ScaleRatio = (TH1F*)tScaleUp->Clone("ScaleRatio");
1473	+	ScaleRatio->Divide(tScaleDown);
1474	+	MatchingRatio->GetYaxis()->SetRangeUser(0,3);
1475	+	MatchingRatio->Draw("e1");
1476	+	TF1 *QP1 = new TF1("QP1","[0]+[1]*x",50,200);//simple linear function ranging from 50 to 200
1477	+	MatchingRatio->Fit("QP1","R");
1478	+	lone->Draw();
1479	+	stringstream summary;
1480	+	summary << " #splitline{Fit result for f(x) = a+bx :}{#splitline{a=" << DigitsAfterComma(QP1->GetParameter(0),4) << " +/- " << DigitsAfterComma(QP1->GetParError(0),4) << "}{b=" << DigitsAfterComma(QP1->GetParameter(1),4);
1481	+	summary << " +/- " << DigitsAfterComma(QP1->GetParError(1),4) << "}}";
1482	+	TText *infobox = write_title(summary.str());
1483	+	infobox->SetX(0.75);
1484	+	infobox->SetTextSize(0.03);
1485	+	infobox->SetY(0.75);
1486	+	infobox->Draw();
1487	+	CompleteSave(kinpad,saveas+"__TTbarComparison_MatchingUpDividedMatchingDown");
1488	+	kinpad->cd();
1489	+	ScaleRatio->GetYaxis()->SetRangeUser(0,3);
1490	+	ScaleRatio->Draw("e1");
1491	+	ScaleRatio->Fit("QP1","R");
1492	+	summary.str("");
1493	+	summary << " #splitline{Fit result for f(x) = a+bx :}{#splitline{a=" << DigitsAfterComma(QP1->GetParameter(0),4) << " +/- " << DigitsAfterComma(QP1->GetParError(0),4) << "}{b=" << DigitsAfterComma(QP1->GetParameter(1),4);
1494	+	summary << " +/- " << DigitsAfterComma(QP1->GetParError(1),4) << "}}";
1495	+	TText *infobox2 = write_title(summary.str());
1496	+	infobox2->SetX(0.75);
1497	+	infobox2->SetTextSize(0.03);
1498	+	infobox2->SetY(0.75);
1499	+	infobox2->Draw();
1500	+
1501	+	lone->Draw();
1502	+	CompleteSave(kinpad,saveas+"__TTbarComparison_ScaleUpDividedScaleDown");
1503	+
1504	+
1505	+	delete QP1;
1506	+	delete infobox;
1507	+	delete infobox2;
1508	+	delete MatchingRatio;
1509	+	delete ScaleRatio;
1510	+	DeleteStack(ScaleUp);
1511	+	DeleteStack(ScaleDown);
1512	+	DeleteStack(MatchingUp);
1513	+	DeleteStack(MatchingDown);
1514	+	delete leg2;
1515	+	CleanLegends();
1516	+	DeleteStack(Central);
1517	+	delete Systematic;
1518	+	delete HScaleUp;
1519	+	delete HScaleDown;
1520	+	delete HMatchingUp;
1521	+	delete HMatchingDown;
1522	+	delete hdata;
1523	+	delete HCentral;
1524	+	delete ca;
1525	+	}
1526	+
1527	+	TH1F* ImposeBinning(TH1F *binninghisto, TH1F* histo) {
1528	+	float val=0,err=0;
1529	+	vector<float> bins;
1530	+	vector<float> vals;
1531	+	vector<float> errs;
1532	+
1533	+	bins.push_back(binninghisto->GetBinLowEdge(1));
1534	+	int iBin=1;
1535	+
1536	+	for(unsigned int i=1;i<histo->GetNbinsX();i++) {//going to second last one on purpose!
1537	+	if((histo->GetBinLowEdge(i)+0.00001>=binninghisto->GetBinLowEdge(iBin)+binninghisto->GetBinWidth(iBin))) {
1538	+	bins.push_back(histo->GetBinLowEdge(i));
1539	+	vals.push_back(val);val=0;
1540	+	errs.push_back(err);err=0;
1541	+	iBin++;
1542	+	}
1543	+	val+=histo->GetBinContent(i);
1544	+	err=sqrt(err*err+histo->GetBinError(i)*histo->GetBinError(i));
1545	+	}
1546	+	bins.push_back(histo->GetBinLowEdge(histo->GetNbinsX())+histo->GetBinWidth(histo->GetNbinsX()));
1547	+	vals.push_back(val);val=0;
1548	+	errs.push_back(err);err=0;
1549	+
1550	+	TH1F *h = new TH1F(("r"+(string)histo->GetName()).c_str(),("r"+(string)histo->GetName()).c_str(),bins.size()-1,&bins[0]);
1551	+	for(unsigned int i=0;i<vals.size();i++) {
1552	+	h->SetBinContent(i+1,vals[i]);
1553	+	h->SetBinError(i+1,errs[i]);
1554	+	}
1555	+
1556	+	h->GetXaxis()->SetTitle(histo->GetXaxis()->GetTitle());
1557	+	h->GetYaxis()->SetTitle(histo->GetYaxis()->GetTitle());
1558	+
1559	+	h->GetXaxis()->CenterTitle();
1560	+	h->GetYaxis()->CenterTitle();
1561	+
1562	+	return h;
1563	+	}
1564	+
1565	+
1566	+	TH1F* ReBinOptimizingStats(int nbins, TH1F *histo) {
1567	+	float statsperbin = (1/(float)nbins) * histo->Integral();
1568	+	float val=0,err=0;
1569	+	vector<float> bins;
1570	+	vector<float> vals;
1571	+	vector<float> errs;
1572	+
1573	+	bins.push_back(histo->GetBinLowEdge(1));
1574	+	for(unsigned int i=1;i<=histo->GetNbinsX();i++) {
1575	+	val+=histo->GetBinContent(i);
1576	+	err=sqrt(err*err+histo->GetBinError(i)*histo->GetBinError(i));
1577	+	if( val/statsperbin > 0.85) {
1578	+	if(bins.size()<nbins) {
1579	+	bins.push_back(histo->GetBinLowEdge(i)+histo->GetBinWidth(i));
1580	+	vals.push_back(val);val=0;
1581	+	errs.push_back(err);err=0;
1582	+	}
1583	+	}
1584	+	if(i==histo->GetNbinsX()) {
1585	+	bins.push_back(histo->GetBinLowEdge(i)+histo->GetBinWidth(i));
1586	+	vals.push_back(val);val=0;
1587	+	errs.push_back(err);err=0;
1588	+	}
1589	+	}
1590	+	TH1F *h = new TH1F(("r"+(string)histo->GetName()).c_str(),("r"+(string)histo->GetName()).c_str(),bins.size()-1,&bins[0]);
1591	+
1592	+	h->GetXaxis()->SetTitle(histo->GetXaxis()->GetTitle());
1593	+	h->GetYaxis()->SetTitle(histo->GetYaxis()->GetTitle());
1594	+
1595	+	h->GetXaxis()->CenterTitle();
1596	+	h->GetYaxis()->CenterTitle();
1597	+
1598	+	for(unsigned int i=0;i<vals.size();i++) {
1599	+	h->SetBinContent(i+1,vals[i]);
1600	+	h->SetBinError(i+1,errs[i]);
1601	+	}
1602	+
1603	+	return h;
1604	+	}
1605
1606	+
1607	+
1608	+	void ShowBinning(TH1F *histo) {
1609	+	cout << "Showing binning for " << histo->GetName() << " (\" " << histo->GetTitle() << "\")" << endl;
1610	+	for(unsigned int i=1;i<=histo->GetNbinsX();i++) cout << "   Bin " << i << " : " << histo->GetBinLowEdge(i) << " , " << histo->GetBinLowEdge(i) + histo->GetBinWidth(i) << "  : " << histo->GetBinContent(i) << endl;
1611	+	}
1612	+
1613	+	void QuickProduceMCSystematicPlot(string variable, int nbins, float bmin, float bmax, string label, TCut thiscut, string saveas) {
1614	+	TCanvas *ca = new TCanvas("ca","ca");
1615	+	TH1F *tScaleUp = systsamples.Draw("tScaleUp",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity,systsamples.FindSample("scaleup"));
1616	+	TH1F *tScaleDown = systsamples.Draw("tScaleDown",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity,systsamples.FindSample("scaledown"));
1617	+	TH1F *tMatchingUp = systsamples.Draw("tMatchingUp",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity,systsamples.FindSample("matchingup"));
1618	+	TH1F *tMatchingDown = systsamples.Draw("tMatchingDown",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity,systsamples.FindSample("matchingdown"));
1619	+	TH1F *tCentral = systsamples.Draw("tCentral",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity,systsamples.FindSample("TTJets_MassiveBinDECAY_TuneZ2sta"));
1620	+
1621	+	tScaleUp->Scale(1.0/tScaleUp->Integral());
1622	+	tScaleDown->Scale(1.0/tScaleDown->Integral());
1623	+	tMatchingDown->Scale(1.0/tMatchingDown->Integral());
1624	+	tMatchingUp->Scale(1.0/tMatchingUp->Integral());
1625	+	tCentral->Scale(1.0/tCentral->Integral());
1626	+
1627	+	cout << "At this point we want to rebin the scale/matching histograms so we get decent stats!" << endl;
1628	+
1629	+
1630	+	TH1F *rtScaleUp = ReBinOptimizingStats(7,tScaleUp); // using the scale up sample to get the binning (shouldn't use the central one for obvious reasons)
1631	+	TH1F *rtScaleDown = ImposeBinning(rtScaleUp,tScaleDown);
1632	+	TH1F *rtMatchingDown = ImposeBinning(rtScaleUp,tMatchingDown);
1633	+	TH1F *rtMatchingUp = ImposeBinning(rtScaleUp,tMatchingUp);
1634	+	TH1F *rtCentral = ImposeBinning(rtScaleUp,tCentral);
1635	+
1636	+	float min=rtScaleUp->GetMinimum();
1637	+	float max=rtScaleUp->GetMaximum();
1638	+
1639	+	if(rtScaleDown->GetMinimum()<min) min=rtScaleDown->GetMinimum();
1640	+	if(rtMatchingDown->GetMinimum()<min) min=rtMatchingDown->GetMinimum();
1641	+	if(rtMatchingDown->GetMinimum()<min) min=rtMatchingDown->GetMinimum();
1642	+	if(rtMatchingUp->GetMinimum()<min) min=rtMatchingUp->GetMinimum();
1643	+	if(rtCentral->GetMinimum()<min) min=rtCentral->GetMinimum();
1644	+
1645	+	if(rtScaleDown->GetMaximum()<min) max=rtScaleDown->GetMaximum();
1646	+	if(rtMatchingDown->GetMaximum()<min) max=rtMatchingDown->GetMaximum();
1647	+	if(rtMatchingDown->GetMaximum()<min) max=rtMatchingDown->GetMaximum();
1648	+	if(rtMatchingUp->GetMaximum()<min) max=rtMatchingUp->GetMaximum();
1649	+	if(rtCentral->GetMaximum()<min) max=rtCentral->GetMaximum();
1650	+
1651	+	rtCentral->SetMaximum(1.2*max);
1652	+	rtCentral->SetMinimum(0.8*min);
1653	+
1654	+	/*  tScaleUp->Rebin(40);
1655	+	tScaleDown->Rebin(40);
1656	+	tMatchingDown->Rebin(40);
1657	+	tMatchingUp->Rebin(40);
1658	+	tCentral->Rebin(40);*/
1659	+
1660	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
1661	+	kinpad->SetLogy(1);
1662	+	kinpad->cd();
1663	+	tCentral->SetFillColor(kWhite);
1664	+	tCentral->SetLineColor(kBlack);
1665	+	tScaleUp->SetLineColor(kRed);
1666	+	tScaleDown->SetLineColor(kRed);
1667	+	tMatchingUp->SetLineColor(kBlue);
1668	+	tMatchingDown->SetLineColor(kBlue);
1669	+	tScaleUp->SetLineStyle(2);
1670	+	tScaleDown->SetLineStyle(3);
1671	+	tMatchingUp->SetLineStyle(2);
1672	+	tMatchingDown->SetLineStyle(3);
1673	+
1674	+	rtScaleDown->SetLineColor(kRed);
1675	+	rtMatchingUp->SetLineColor(kBlue);
1676	+	rtMatchingDown->SetLineColor(kBlue);
1677	+	rtScaleUp->SetLineStyle(2);
1678	+	rtScaleDown->SetLineStyle(3);
1679	+	rtMatchingUp->SetLineStyle(2);
1680	+	rtMatchingDown->SetLineStyle(3);
1681	+
1682	+	TLegend *leg2 = make_legend();
1683	+	leg2->AddEntry(tCentral,"Central (ttbar)","l");
1684	+	leg2->AddEntry(tScaleUp,"ScaleUp (ttbar)","l");
1685	+	leg2->AddEntry(tScaleDown,"ScaleDown (ttbar)","l");
1686	+	leg2->AddEntry(tMatchingUp,"MatchingUp (ttbar)","l");
1687	+	leg2->AddEntry(tMatchingDown,"MatchingDown (ttbar)","l");
1688	+
1689	+	tCentral->Draw("histo");
1690	+	tScaleUp->Draw("histo,same");
1691	+	tScaleDown->Draw("histo,same");
1692	+	tMatchingUp->Draw("histo,same");
1693	+	tMatchingDown->Draw("histo,same");
1694	+	leg2->Draw();
1695	+	DrawMCPrelim();
1696	+	CompleteSave(kinpad,saveas+"__TTbarComparison");
1697	+	kinpad->cd();
1698	+	kinpad->SetLogy(0);
1699	+	rtCentral->Draw("histo");
1700	+	rtScaleUp->Draw("histo,same");
1701	+	rtScaleDown->Draw("histo,same");
1702	+	rtMatchingUp->Draw("histo,same");
1703	+	rtMatchingDown->Draw("histo,same");
1704	+	leg2->Draw();
1705	+	DrawMCPrelim();
1706	+	CompleteSave(kinpad,saveas+"__TTbarComparison__REBINNED");
1707	+
1708	+
1709	+	gStyle->SetOptFit(0);
1710	+
1711	+	kinpad->cd();
1712	+	kinpad->SetLogy(0);
1713	+	TH1F *MatchingRatio = (TH1F*)rtMatchingUp->Clone("MatchingRatio");
1714	+	MatchingRatio->Divide(rtMatchingDown);
1715	+	TLine *lone = new TLine(tScaleUp->GetBinLowEdge(1),1,tScaleUp->GetBinLowEdge(tScaleUp->GetNbinsX())+tScaleUp->GetBinWidth(tScaleUp->GetNbinsX()),1);
1716	+	lone->SetLineColor(TColor::GetColor("#01DF01"));
1717	+	lone->SetLineStyle(2);
1718	+	TH1F *ScaleRatio = (TH1F*)rtScaleUp->Clone("ScaleRatio");
1719	+	ScaleRatio->Divide(rtScaleDown);
1720	+	MatchingRatio->GetYaxis()->SetRangeUser(0,3);
1721	+	MatchingRatio->Draw("e1");
1722	+	TF1 *QP1 = new TF1("QP1","[0]+[1]*x",50,200);//simple linear function ranging from 50 to 200
1723	+	if(variable=="genMET" || variable=="met[4]" || variable=="met[1]") MatchingRatio->Fit("QP1","RQ");
1724	+	lone->Draw();
1725	+	stringstream summary;
1726	+	summary << " #splitline{Fit result for f(x) = a+bx :}{#splitline{a=" << DigitsAfterComma(QP1->GetParameter(0),4) << " +/- " << DigitsAfterComma(QP1->GetParError(0),4) << "}{b=" << DigitsAfterComma(QP1->GetParameter(1),4);
1727	+	summary << " +/- " << DigitsAfterComma(QP1->GetParError(1),4) << "}}";
1728	+	TText *infobox = write_title(summary.str());
1729	+	infobox->SetX(0.75);
1730	+	infobox->SetTextSize(0.03);
1731	+	infobox->SetY(0.75);
1732	+	if(variable=="genMET" || variable=="met[4]" || variable=="met[1]") infobox->Draw();
1733	+	DrawMCPrelim();
1734	+	CompleteSave(kinpad,saveas+"__TTbarComparison_MatchingUpDividedMatchingDown");
1735	+	kinpad->cd();
1736	+	ScaleRatio->GetYaxis()->SetRangeUser(0,3);
1737	+	ScaleRatio->Draw("e1");
1738	+	if(variable=="genMET" || variable=="met[4]" || variable=="met[1]") ScaleRatio->Fit("QP1","RQ");
1739	+	summary.str("");
1740	+	summary << " #splitline{Fit result for f(x) = a+bx :}{#splitline{a=" << DigitsAfterComma(QP1->GetParameter(0),4) << " +/- " << DigitsAfterComma(QP1->GetParError(0),4) << "}{b=" << DigitsAfterComma(QP1->GetParameter(1),4);
1741	+	summary << " +/- " << DigitsAfterComma(QP1->GetParError(1),4) << "}}";
1742	+	TText *infobox2 = write_title(summary.str());
1743	+	infobox2->SetX(0.75);
1744	+	infobox2->SetTextSize(0.03);
1745	+	infobox2->SetY(0.75);
1746	+	if(variable=="genMET" || variable=="met[4]" || variable=="met[1]") infobox2->Draw();
1747	+
1748	+	lone->Draw();
1749	+	DrawMCPrelim();
1750	+	CompleteSave(kinpad,saveas+"__TTbarComparison_ScaleUpDividedScaleDown");
1751	+
1752	+	kinpad->cd();
1753	+	TH1F *SysMatching = (TH1F*)rtMatchingDown->Clone("SysMatching");
1754	+	TH1F *SysScale = (TH1F*)rtScaleDown->Clone("SysMatching");
1755	+	TH1F *SysHisto = (TH1F*)rtScaleDown->Clone("SysHisto");
1756	+
1757	+	for(int i=1;i<=SysScale->GetNbinsX();i++) {
1758	+	float matching, scale;
1759	+	bool AssumeFactor=false;
1760	+	if(AssumeFactor) {
1761	+	// assume that an upward/downward variation means a change by a constant factor, i.e. f(up) = alpha * f(central)
1762	+	matching=0.;
1763	+	scale=0.;
1764	+	} else {
1765	+	// assume that the central value is exactly between up & down, i.e. central = (up-down)/2
1766	+	matching = fabs(    ( rtMatchingUp->GetBinContent(i) - rtMatchingDown->GetBinContent(i)) / (rtMatchingUp->GetBinContent(i) + rtMatchingDown->GetBinContent(i)) );
1767	+	scale    = fabs(    ( rtScaleUp->GetBinContent(i) - rtScaleDown->GetBinContent(i)) / (rtScaleUp->GetBinContent(i) + rtScaleDown->GetBinContent(i)) );
1768	+	}
1769	+
1770	+	SysMatching->SetBinContent(i,1+matching);
1771	+	SysScale->SetBinContent(i,1+scale);
1772	+	SysHisto->SetBinContent(i,sqrt(matching*matching+scale*scale)+1);
1773	+	}
1774	+
1775	+	SysHisto->SetLineColor(kGreen);
1776	+
1777	+	SysMatching->SetMaximum(3.0);
1778	+	SysMatching->SetMinimum(0.0);
1779	+	SysMatching->Draw("histo");
1780	+	SysScale->Draw("histo,same");
1781	+	SysHisto->Draw("histo,same");
1782	+	CompleteSave(kinpad,saveas+"__TTbarComparison_DerivedSystematic");
1783	+
1784	+	delete SysMatching;
1785	+	delete SysScale;
1786	+
1787	+	if(!Contains(variable,"gen")) {
1788	+	TH1F *hdata  = allsamples.Draw("hdata",   variable,nbins,bmin,bmax,label,"events",thiscut,data,luminosity);
1789	+	THStack smc = allsamples.DrawStack("smc",variable,nbins,bmin,bmax,label,"events",thiscut,mc,luminosity);
1790	+
1791	+	TH1F *Systematic = (TH1F*)hdata->Clone("Systematic");
1792	+
1793	+	for(int i=1;i<=Systematic->GetNbinsX();i++) {
1794	+	int iBin=SysHisto->FindBin(Systematic->GetBinCenter(i));
1795	+	if(iBin>SysHisto->GetNbinsX()) --iBin; // if we're over the end of the histo, use the last bin
1796	+	float sys = SysHisto->GetBinContent(iBin)-1;
1797	+	Systematic->SetBinContent(i,sys);
1798	+	}
1799	+
1800	+	kinpad->cd();
1801	+	TLegend *leg = allsamples.allbglegend();
1802	+	hdata->Draw();
1803	+	smc.Draw("histo,same");
1804	+	hdata->Draw("same");
1805	+	leg->Draw();
1806	+	DrawPrelim();
1807	+	save_with_ratio_and_sys_band( 0.0, hdata, CollapseStack(smc), kinpad->cd(), saveas+"__DataVsMC", false, false, "data/mc",Systematic );
1808	+
1809	+	delete Systematic;
1810	+	delete hdata;
1811	+	}
1812	+
1813	+
1814	+	delete SysHisto;
1815	+	delete tScaleUp;
1816	+	delete tScaleDown;
1817	+	delete tCentral;
1818	+	delete tMatchingUp;
1819	+	delete tMatchingDown;
1820	+	delete rtScaleUp;
1821	+	delete rtScaleDown;
1822	+	delete rtCentral;
1823	+	delete rtMatchingUp;
1824	+	delete rtMatchingDown;
1825	+	delete QP1;
1826	+	delete infobox;
1827	+	delete infobox2;
1828	+	delete MatchingRatio;
1829	+	delete ScaleRatio;
1830	+	delete leg2;
1831	+	CleanLegends();
1832	+	delete ca;
1833	+	}
1834	+
1835	+	void ProduceMCSystematicPlots() {
1836	+	cout << "Getting ready to produce systematic plots " << endl;
1837	+	TCut cutweightBKP = cutWeight;
1838	+
1839	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_SameFlavor");
1840	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_OppositeFlavor");
1841	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_SameFlavor");
1842	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_OppositeFlavor");
1843	+
1844	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE");
1845	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE");
1846	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE_HighMass");
1847	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE_HighMass");
1848	+
1849	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_SameFlavor_INCLUSIVE");
1850	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE");
1851	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/MET/ETH_SameFlavor_INCLUSIVE_HighMass");
1852	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE_HighMass");
1853	+	QuickProduceMCSystematicPlot("pfJetGoodNum40",8,-0.5,7.5,"NJets",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70&&met[4]>100"),"MCSystPlots/NJets/ETH_SameFlavor");
1854	+	QuickProduceMCSystematicPlot("pfJetGoodNum40",8,-0.5,7.5,"NJets",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70&&met[4]>100"),"MCSystPlots/NJets/ETH_OppositeFlavor");
1855	+
1856	+
1857	+
1858	+	cout << "Kicking cutWeight  " << (const char*) cutWeight << endl;
1859	+	cout << "Keeping OSSF cut " << (const char*)cutOSSF << endl;
1860	+	cutWeight="1.0";
1861	+
1862	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE__NOPURW");
1863	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE__NOPURW");
1864	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1865	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/genMET/ETH_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");
1866	+
1867	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_SameFlavor_INCLUSIVE__NOPURW");
1868	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE__NOPURW");
1869	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/MET/ETH_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1870	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/MET/ETH_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");
1871	+
1872	+
1873	+	// --------------------------------------    ***** AACHEN *****      --------------------------------------
1874	+
1875	+	/*
1876	+	cutWeight=cutweightBKP;
1877	+	TCut essentialcutBKP = essentialcut;
1878	+	cout << (const char*) essentialcut << endl;
1879	+
1880	+	essentialcut = TCut((ReplaceAll((const char*)essentialcut,"pt2>20","pt2>10")).c_str());
1881	+	essentialcut = TCut((ReplaceAll((const char*)essentialcut,"abs(eta1)<1.4","abs(eta1)<2.4")).c_str());
1882	+	essentialcut = TCut((ReplaceAll((const char*)essentialcut,"abs(eta2)<1.4","abs(eta2)<2.4")).c_str());
1883	+
1884	+	TCut cutnJetsBKP = cutnJets;
1885	+	cutnJets = TCut((ReplaceAll((const char*)cutnJets,"pt2>20","pt2>10")).c_str());
1886	+	cutnJets = TCut((ReplaceAll((const char*)cutnJets,"pfJetGoodNum40>=3","pfJetGoodNum40>=2")).c_str());
1887	+	cutnJets = TCut((ReplaceAll((const char*)cutnJets,"abs(eta2)<1.4","abs(eta2)<2.4")).c_str());
1888	+	cutnJets = TCut((ReplaceAll((const char*)cutnJets,"abs(eta1)<1.4","abs(eta1)<2.4")).c_str());
1889	+
1890	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_SameFlavor");
1891	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_OppositeFlavor");
1892	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_SameFlavor");
1893	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_OppositeFlavor");
1894	+
1895	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_SameFlavor_INCLUSIVE");
1896	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_OppositeFlavor_INCLUSIVE");
1897	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/genMET/Aachen_SameFlavor_INCLUSIVE_HighMass");
1898	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/genMET/Aachen_OppositeFlavor_INCLUSIVE_HighMass");
1899	+
1900	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_SameFlavor_INCLUSIVE");
1901	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_OppositeFlavor_INCLUSIVE");
1902	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/MET/Aachen_SameFlavor_INCLUSIVE_HighMass");
1903	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/MET/Aachen_OppositeFlavor_INCLUSIVE_HighMass");
1904	+	QuickProduceMCSystematicPlot("pfJetGoodNum40",8,-0.5,7.5,"NJets",cutOSSF&&cutnJets&&TCut("mll>20&&mll<70&&met[4]>100"),"MCSystPlots/NJets/Aachen_SameFlavor");
1905	+	QuickProduceMCSystematicPlot("pfJetGoodNum40",8,-0.5,7.5,"NJets",cutOSOF&&cutnJets&&TCut("mll>20&&mll<70&&met[4]>100"),"MCSystPlots/NJets/Aachen_OppositeFlavor");
1906	+
1907	+
1908	+
1909	+	cout << "Kicking cutWeight  " << (const char*) cutWeight << endl;
1910	+	cout << "Keeping OSSF cut " << (const char*)cutOSSF << endl;
1911	+	cutWeight="1.0";
1912	+
1913	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_SameFlavor_INCLUSIVE__NOPURW");
1914	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/genMET/Aachen_OppositeFlavor_INCLUSIVE__NOPURW");
1915	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/genMET/Aachen_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1916	+	QuickProduceMCSystematicPlot("genMET",20,0,200,"gen#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/genMET/Aachen_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");
1917	+
1918	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_SameFlavor_INCLUSIVE__NOPURW");
1919	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>20&&mll<70"),"MCSystPlots/MET/Aachen_OppositeFlavor_INCLUSIVE__NOPURW");
1920	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSSF&&TCut("mll>120"),"MCSystPlots/MET/Aachen_SameFlavor_INCLUSIVE_HighMass__NOPURW");
1921	+	QuickProduceMCSystematicPlot("met[4]",20,0,200,"PF#slash{E}_{T} [GeV]",cutOSOF&&TCut("mll>120"),"MCSystPlots/MET/Aachen_OppositeFlavor_INCLUSIVE_HighMass__NOPURW");
1922	+
1923	+	cutnJets  = cutnJetsBKP;
1924	+	cutWeight = cutweightBKP;
1925	+	essentialcut = essentialcutBKP;
1926	+	*/
1927	+	write_error(__FUNCTION__,"Still need to add systematic shape");assert(0);
1928		}
1929
1930	+	void MakeTauPlot(string label, TCut pcut, string filename, bool TauOnly) {
1931	+	TCut cut=pcut&&basiccut;
1932	+	string variable="mll";
1933	+	string xlabel="m_{ll} [GeV]";
1934	+	int nbins=40;
1935	+	float min=0;
1936	+	float max=200;
1937	+
1938	+	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
1939	+	TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cut&&cutOSSF,data,luminosity);
1940	+	TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cut&&cutOSOF,data,luminosity);
1941	+	TH1F *tauSF = allsamples.Draw("tauSF",variable,nbins,min,max, xlabel, "events",cut&&cutOSSF,data,luminosity,systsamples.FindSample("DYToTauTau"));
1942	+	TH1F *tauOF = allsamples.Draw("tauOF",variable,nbins,min,max, xlabel, "events",cut&&cutOSOF,data,luminosity,systsamples.FindSample("DYToTauTau"));
1943	+	datahistoSF->SetMarkerSize(DataMarkerSize);
1944	+	datahistoSF->SetMaximum(1.3*datahistoSF->GetMaximum());
1945	+	datahistoOF->SetMarkerSize(DataMarkerSize);
1946	+	datahistoOF->SetMaximum(1.3*datahistoOF->GetMaximum());
1947	+
1948	+	THStack mcstackSF   = allsamples.DrawStack("mcstackSF",  variable,nbins,min,max,xlabel,"events",cut&&cutOSSF,mc,luminosity);
1949	+	THStack mcstackOF   = allsamples.DrawStack("mcstackOF",  variable,nbins,min,max,xlabel,"events",cut&&cutOSOF,mc,luminosity);
1950	+	datahistoSF->Draw("e1");
1951	+	ckin->Update();
1952	+	mcstackSF.Draw("histo,same");
1953	+
1954	+	datahistoSF->Draw("same,e1");
1955	+	TLegend *kinleg = allsamples.allbglegend();
1956	+	kinleg->Draw();
1957	+
1958	+	TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
1959	+	kinpad->cd();
1960	+	datahistoSF->Draw("e1");
1961	+	mcstackSF.Draw("histo,same");
1962	+	datahistoSF->Draw("same,e1");
1963	+	datahistoSF->Draw("same,axis");
1964	+	kinleg->Draw();
1965	+	DrawPrelim();
1966	+	Save_With_Ratio(datahistoSF,mcstackSF,kinpad->cd(),filename+"__SF_mc");
1967	+
1968	+	kinpad->cd();
1969	+	datahistoOF->Draw("e1");
1970	+	mcstackOF.Draw("histo,same");
1971	+	datahistoOF->Draw("same,e1");
1972	+	datahistoOF->Draw("same,axis");
1973	+	kinleg->Draw();
1974	+	DrawPrelim();
1975	+	Save_With_Ratio(datahistoOF,mcstackOF,kinpad->cd(),filename+"__OF_mc");
1976	+
1977	+	kinpad->cd();
1978	+	tauSF->Draw("e1");
1979	+	tauOF->Draw("histo,same");
1980	+	tauSF->Draw("e1,same");
1981	+	TLegend *legtau = make_legend();
1982	+	legtau->AddEntry(tauSF,"DY->#tau#tau, SF","p");
1983	+	legtau->AddEntry(tauSF,"DY->#tau#tau, OF","l");
1984	+	legtau->Draw();
1985	+	DrawPrelim();
1986	+	Save_With_Ratio(tauSF,tauOF,kinpad->cd(),filename+"__Tau_SF_vs_OF");
1987	+
1988	+	delete datahistoSF;
1989	+	delete datahistoOF;
1990	+	delete tauSF;
1991	+	delete tauOF;
1992	+	delete kinpad;
1993	+	delete ckin;
1994	+	CleanLegends();
1995	+	}
1996	+
1997	+
1998	+	void TauQuestion() {
1999	+	//   MakeTauPlot("MET>100 GeV, #geq 3 jets",cutnJets&&TCut("met[4]>100"),"TauQ/MET100_3Jets",true);
2000	+	//   MakeTauPlot("MET>100 GeV",TCut("met[4]>100"),"TauQ/MET100",true);
2001	+	//   MakeTauPlot("MET>0 GeV",TCut("met[4]>0"),"TauQ/MET0",true);
2002	+	MakeTauPlot("b-tag veto, 50<MET<100",TCut("pfJetGoodNumBtag30==0&&met[4]>50&&met[4]<100"),"TauQ/ControlPlots",false);
2003	+	}
2004	+
2005	+	void ProduceOFSFPlots(string mcjzb, string datajzb, bool isAachen) {
2006	+	string BaseMetCut="met[4]>100";
2007	+	string Prefix="ETH/";
2008	+	if(isAachen) {
2009	+	BaseMetCut="met[4]>150";
2010	+	Prefix="Aachen/";
2011	+	}
2012	+	//        make_OFSF_plots("st", BaseMetCut, 20, 0, 2000., false, "S_{T}", Prefix+"FrankT");
2013	+	make_OFSF_plots("mt2j", BaseMetCut, 35, 0, 700., false, "M_{T2J}", Prefix+"MT2J");
2014	+
2015	+	make_OFSF_plots("ml1b", BaseMetCut, 20, 0, 800., false, "m_{l_{1}b}", Prefix+"ML1B");
2016	+	make_OFSF_plots("ml2b", BaseMetCut, 20, 0, 400., false, "m_{l_{2}b}", Prefix+"ML2B");
2017	+	make_OFSF_plots("min(ml1b,ml2b)", BaseMetCut, 20, 0, 400., false, "m_{lb}", Prefix+"MLB");
2018	+
2019	+	//      make_OFSF_plots(datajzb, BaseMetCut, 55, -100, 450., false, "JZB", Prefix+"JZB_In_Signal_Region");
2020	+	make_OFSF_plots("mt2", BaseMetCut, 15, 0, 150., false, "M_{T2}", Prefix+"MT2");
2021	+	make_OFSF_plots("d2", BaseMetCut, 20, 0, 200., true, "D2", Prefix+"D2");
2022	+	make_OFSF_plots("pfJetGoodNum40", BaseMetCut, 10, -0.5, 9.5, true, "N_{jets}", Prefix+"nJets");
2023	+
2024	+	//      make_OFSF_plots("mll", BaseMetCut, 60, 20., 320., false, "m_{ll}", Prefix+"mll");
2025	+	//      make_OFSF_plots("leptonNum", BaseMetCut, 3, 2, 5., false, "N(leptons)", Prefix+"NLeptons");
2026	+
2027	+	//      make_OFSF_plots("pfJetGoodNum40", BaseMetCut,                        7, 3, 10, true, "#(jets)", Prefix+"njets");
2028	+	//      make_OFSF_plots("pfJetGoodNum40", BaseMetCut+"&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", Prefix+"njets_btagVeto");
2029	+	//      make_OFSF_plots("pfJetGoodNum40", BaseMetCut"&&pfJetGoodNumBtag30>0",  7, 3, 10, true, "#(jets)", Prefix+"njets_AtLeastOneBJet30");
2030	+
2031	+	//      make_OFSF_plots("pfJetGoodNumBtag30", BaseMetCut, 5, -0.5, 4.5, true, "N_{b-jets}", Prefix+"nbjets");
2032	+	//      make_OFSF_plots("pfJetGoodPtBtag[0]", BaseMetCut+"&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", Prefix+"ptb1");
2033	+	//      make_OFSF_plots("pfJetGoodPt",  BaseMetCut, 30, 0., 400., true, "p_{T}^{J}", Prefix+"ptJet");
2034	+
2035	+	//      make_OFSF_plots("iso1", BaseMetCut, 20, 0, 0.3, true, "lepton 1 isolation", Prefix+"iso1");
2036	+	//      make_OFSF_plots("iso2", BaseMetCut, 20, 0, 0.3, true, "lepton 2 isolation", Prefix+"iso2");
2037	+	make_OFSF_plots("pt",  BaseMetCut, 30, 0., 300., true, "p_{ll}", Prefix+"pt");
2038	+	make_OFSF_plots("pt1",  BaseMetCut, 30, 0., 300., true, "p_{T,1}", Prefix+"pt1");
2039	+	make_OFSF_plots("pt2",  BaseMetCut, 22, 0., 220., true, "p_{T,2}", Prefix+"pt2");
2040	+	//      make_OFSF_plots("eta1", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{1}", Prefix+"eta1", 0.15);
2041	+	//      make_OFSF_plots("eta2", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{2}", Prefix+"eta2", 0.15);
2042	+	//      make_OFSF_plots("phi1", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", Prefix+"phi1", 0.2);
2043	+	//      make_OFSF_plots("phi2", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", Prefix+"phi2", 0.2);
2044	+	//      make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", BaseMetCut, 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", Prefix+"jpt1pt2", 0.2);
2045	+	//      make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", BaseMetCut, 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", Prefix+"dphij1met", 0.2);
2046	+
2047	+	//      make_OFSF_plots("met[4]",  BaseMetCut, 30, 100., 400., true, "E_{T}^{miss}", Prefix+"MET");
2048	+	//      make_OFSF_plots("TMath::Abs(dphi)", BaseMetCut, 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", Prefix+"dphi", 0.2);
2049	+	//      make_OFSF_plots("TMath::Abs(dphiMet1)", BaseMetCut, 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", Prefix+"dphiMet1", 0.2);
2050	+	//      make_OFSF_plots("TMath::Abs(dphiMet2)", BaseMetCut, 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", Prefix+"dphiMet2", 0.2);
2051	+	//      make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", BaseMetCut, 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", Prefix+"dphilc");
2052	+	//      make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", BaseMetCut, 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", Prefix+"dphijc");
2053	+	//      make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", BaseMetCut, 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", Prefix+"dphilco");
2054	+	//      make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", BaseMetCut, 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(jet,MET)|)", Prefix+"dphijco");
2055	+
2056	+	}
2057	+
2058	+	void ProduceOFSFPlots(string mcjzb, string datajzb) {
2059	+
2060	+
2061	+	//FIRST : ETH SELECTION
2062	+	ProduceOFSFPlots(mcjzb,datajzb,0); // this is the ETH selection
2063	+
2064	+
2065	+
2066	+	string backup_cutnJets = (const char*)cutnJets;
2067	+	cutnJets = ReplaceInCut(cutnJets,"pfJetGoodNum40>=3","pfJetGoodNum40>=2");
2068	+	cutnJets = ReplaceInCut(cutnJets,"1.4","2.4");
2069	+
2070	+	string backup_basiccut = (const char*)basiccut;
2071	+	basiccut = ReplaceInCut(basiccut,"1.4","2.4");
2072	+
2073	+	string backup_essential = (const char*) essentialcut;
2074	+	essentialcut = ReplaceInCut(essentialcut,"1.4","2.4");
2075	+
2076	+	ProduceOFSFPlots(mcjzb,datajzb,1); // this is the Aachen selection
2077	+
2078	+	cutnJets     = TCut(backup_cutnJets.c_str());
2079	+	basiccut     = TCut(backup_basiccut.c_str());
2080	+	essentialcut = TCut(backup_essential.c_str());
2081	+	}
2082	+
2083	+
2084		void do_kinematic_plots(string mcjzb, string datajzb, bool doPF=false)
2085		{
2086	+	//  switch_overunderflow(true);
2087		bool dolog=true;
2088		bool nolog=false;
2089
2090	<	bool doOFSF = false;
2090	>	bool doOFSF = true;
2091		bool doKin  = false;
2092	<	bool doDataComp = true;
2092	>	bool doDataComp = false;
2093	>	bool MakeTwoThreeJetComparison = false;
2094
2095
2096		if(doPF) write_warning(__FUNCTION__,"Please use caution when trying to produce PF plots; not all versions of the JZB trees have these variables!");
#	Line 750 | Line 2100 | void do_kinematic_plots(string mcjzb, st
2100		mll_low=20;
2101		mll_hi=320;
2102		}
2103	+
2104	+	if(MakeTwoThreeJetComparison) MakeElegantTwoThreeComparisons();
2105	+
2106	+	//TauQuestion();
2107
2108	<	if ( doOFSF ) {
2109	<	make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll");
2110	<
2111	<	make_OFSF_plots("pfJetGoodNum40",     "met[4]>100", 7, 3, 10, true, "#(jets)", "njets");
2112	<	make_OFSF_plots("pfJetGoodNum40",     "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto");
2113	<	make_OFSF_plots("pfJetGoodNum40",     "met[4]>100&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30");
2114	<
761	<	make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets");
762	<	make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1");
763	<	make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2");
764	<	//  make_OFSF_plots("pt1",  "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1");
765	<	//  make_OFSF_plots("pt2",  "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2");
766	<	make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", 0.15);
767	<	make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", 0.15);
768	<	//  make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", 0.2);
769	<	//  make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", 0.2);
770	<	//  make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", 0.2);
771	<	make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", 0.2);
772	<	make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", 0.2);
773	<	//  make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", 0.2);
774	<	//  make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", 0.2);
775	<	make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilc");
776	<	make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijc");
777	<	make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", "dphilco");
778	<	make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(jet,MET)|)", "dphijco");
779	<	}
2108	>	//ProduceJanPlots();
2109	>	//ProduceMCSystematicPlots();
2110	>	//assert(0);
2111	>	//  make_plain_kin_plot("pt",Cut2Str(cutOSSF&&TCut("mll>20&&pfJetGoodNumBtag30>=2")),40,0,200,nolog,"Z p_{T}","PlotForKostas",doPF,true);
2112	>
2113	>	if ( doOFSF ) ProduceOFSFPlots(mcjzb,datajzb);
2114	>	assert(0);
2115
2116	<	if ( doDataComp && !PlottingSetup::openBox ) {
2116	>	if ( doDataComp) {
2117		TCut mllCut("");
2118	<	if ( !PlottingSetup::openBox ) mllCut = "mll>120";
2118	>	float massmin = 15.;
2119	>	float massmax = 315;
2120	>	int massnbins = 60;
2121	>	/*    if ( !PlottingSetup::openBox ) {
2122	>	mllCut = "mll>70";
2123	>	massmin = 120;
2124	>	massmax = 360;
2125	>	massnbins = 14;
2126	>	}*/
2127
2128		TCut cutSignal = cutmass&&cutnJets&&"met[4]>100";
2129	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut,           60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig");
2130	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee");
2131	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm");
2132	<
2133	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0",           60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto");
2134	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_btagVeto");
2135	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_btagVeto");
2136	<
2137	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0",           60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet");
2138	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_AtLeastOneBJet");
2139	<	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_AtLeastOneBJet");
2129	>
2130	>	make_data_comparison_plot("numVtx", "",40, -0.5, 39.5,-1., true, "N(Vertices)", "numVtx");
2131	>	make_data_comparison_plot("pfJetGoodEta[0]", "",50, -3.0, 3.0,-1., true, "N(Jets)", "pfJetGoodEta0");
2132	>
2133	>
2134	>	make_data_comparison_plot("pfJetGoodNum40", cutOSOF,10, -0.5, 9.5,-1., true, "N(Jets)", "njets_OF");
2135	>	make_data_comparison_plot("met[4]", cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_inclusive");
2136	>	make_data_comparison_plot("met[1]", cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_inclusive");
2137	>	make_data_comparison_plot("met[4]", cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_inclusive");
2138	>	make_data_comparison_plot("met[1]", cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_inclusive");
2139	>
2140	>	make_data_comparison_plot("met[4]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_exclusive");
2141	>	make_data_comparison_plot("met[1]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_exclusive");
2142	>	make_data_comparison_plot("met[4]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_exclusive");
2143	>	make_data_comparison_plot("met[1]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_exclusive");
2144	>
2145	>
2146	>	string basiccutsaved = (const char*)basiccut;
2147	>
2148	>	string newbasiccut=basiccutsaved;
2149	>
2150	>	size_t found;
2151	>	found = newbasiccut.find( "pt2>20" );
2152	>	while (found!=string::npos){
2153	>	newbasiccut.replace( found, string( "pt2>20" ).length(), "pt2>10" );
2154	>	found = newbasiccut.find( "pt2>20" );
2155	>	}
2156	>	basiccut=TCut(newbasiccut.c_str());
2157	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSSF,60, 0., 300.,-1., false, "mll", "mll_SF_Aachen_pt2010_met");
2158	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSOF,60, 0., 300.,-1., false, "mll", "mll_OF_Aachen_pt2010_met");
2159	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_Aachen_pt2010_t1cpfmet");
2160	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_Aachen_pt2010_t1cpfmet");
2161	>	basiccut=TCut(basiccutsaved.c_str());
2162	>
2163	>	make_data_comparison_plot("pfJetGoodNum40", cutOSSF, 8, 0., 8.,-1., true, "#(jets)", "njets_SF");
2164	>	make_data_comparison_plot("pfJetGoodNum40", cutOSOF, 8, 0., 8.,-1., true, "#(jets)", "njets_OF");
2165	>
2166	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_met");
2167	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_met");
2168	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_t1cpfmet");
2169	>	make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_t1cpfmet");
2170	>
2171	>
2172	>	make_data_comparison_plot("mll", cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive");
2173	>	make_data_comparison_plot("mll", cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_ee");
2174	>	make_data_comparison_plot("mll", cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_mm");
2175	>	make_data_comparison_plot("mll", cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_inclusive");
2176	>	make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_ee");
2177	>	make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_mm");
2178	>	make_data_comparison_plot("mll", cutnJets&&cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive");
2179	>	make_data_comparison_plot("mll", cutnJets&&cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_exclusive");
2180	>
2181	>	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig");
2182	>	make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig");
2183	>
2184	>	make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSSF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_central");
2185	>	make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSOF&&cutSignal&&mllCut,           massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_central");
2186	>
2187	>	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto");
2188	>	make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_btagVeto");
2189	>
2190	>	make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0",  massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet");
2191	>	make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0",  massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_AtLeastOneBJet");
2192	>
2193	>
2194	>	make_data_comparison_plot("mll", mllCut&&cutOSOF&&cutSignal, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig");
2195	>	make_data_comparison_plot("mll", cutmass&&cutOSSF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_CR");
2196	>	make_data_comparison_plot("mll", cutmass&&cutOSOF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_CR");
2197	>
2198	>	make_data_comparison_plot("pfJetGoodNum40", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_SF_sig");
2199	>	make_data_comparison_plot("pfJetGoodNum40", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_OF_sig");
2200	>	make_data_comparison_plot("pfJetGoodNumBtag30", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_SF_sig");
2201	>	make_data_comparison_plot("pfJetGoodNumBtag30", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_OF_sig");
2202
798	–	make_data_comparison_plot("mll", cutOSOF&&cutSignal, 60, 20., 320.,-1., false, "m_{ll}", "mll_OF_sig");
799	–	make_data_comparison_plot("mll", cutmass&&"met[4]>100&&met[4]<150&&pfJetGoodNum40<3", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_CR");
2203		}
2204
2205
2206		if ( doKin ) {
2207		string mllCut("");
2208	<	if ( !PlottingSetup::openBox ) mllCut = "&&mll>120";
2208	>
2209	>	make_plain_kin_plot("mll",Cut2Str(cutOSSF&&TCut("mll>20")),56,20,300,dolog,"m_{ll}","NoPreselection/mll_SF",doPF,true);
2210	>	make_plain_kin_plot("mll",Cut2Str(cutOSOF&&TCut("mll>20")),56,20,300,dolog,"m_{ll}","NoPreselection/mll_OF",doPF,true);
2211
2212	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("mll>20")),80,0,400,dolog,"E_{T}^{miss}","NoPreselection/met_SF",doPF,true);
2213	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>20")),80,0,400,dolog,"E_{T}^{miss}","NoPreselection/met_OF",doPF,true);
2214	+
2215	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&TCut("mll>20")),20,0,800,dolog,"H_{T}","NoPreselection/HT_SF",doPF,true);
2216	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&TCut("mll>20")),20,0,800,dolog,"H_{T}","NoPreselection/HT_OF",doPF,true);
2217	+
2218	+	make_plain_kin_plot("pfJetGoodNum40",Cut2Str(cutOSSF&&TCut("mll>20")),10,-0.5,9.5,dolog,"n_{jets}","NoPreselection/NJets_SF",doPF,true);
2219	+	make_plain_kin_plot("pfJetGoodNum40",Cut2Str(cutOSOF&&TCut("mll>20")),10,-0.5,9.5,dolog,"n_{jets}","NoPreselection/NJets_OF",doPF,true);
2220	+
2221	+	make_plain_kin_plot("pfJetGoodNumBtag30",Cut2Str(cutOSSF&&TCut("mll>20")),10,-0.5,9.5,dolog,"n_{btags}","NoPreselection/NBTags_SF",doPF,true);
2222	+	make_plain_kin_plot("pfJetGoodNumBtag30",Cut2Str(cutOSOF&&TCut("mll>20")),10,-0.5,9.5,dolog,"n_{btags}","NoPreselection/NBTags_OF",doPF,true);
2223	+
2224	+	make_plain_kin_plot("pt1",Cut2Str(cutOSSF&&TCut("mll>20")),80,0,400,dolog,"p_{T}^{1}","NoPreselection/Pt1_SF",doPF,true);
2225	+	make_plain_kin_plot("pt1",Cut2Str(cutOSOF&&TCut("mll>20")),80,0,400,dolog,"p_{T}^{1}","NoPreselection/Pt1_OF",doPF,true);
2226	+
2227	+	make_plain_kin_plot("pt2",Cut2Str(cutOSSF&&TCut("mll>20")),80,0,400,dolog,"p_{T}^{2}","NoPreselection/Pt2_SF",doPF,true);
2228	+	make_plain_kin_plot("pt2",Cut2Str(cutOSOF&&TCut("mll>20")),80,0,400,dolog,"p_{T}^{2}","NoPreselection/Pt2_OF",doPF,true);
2229	+
2230	+
2231	+	assert(0);
2232	+	make_plain_kin_plot("mll",Cut2Str(cutOSSF&&TCut("mll>60&&mll<120")),60,60,120,dolog,"m_{ll}","TestingLuminosity_SF",doPF,true);
2233	+	make_plain_kin_plot("mll",Cut2Str(cutOSOF&&TCut("mll>60&&mll<120")),60,60,120,dolog,"m_{ll}","TestingLuminosity_OF",doPF,true);
2234	+
2235	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_0j",doPF,true);
2236	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_1j",doPF,true);
2237	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_2j",doPF,true);
2238	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_3j",doPF,true);
2239	+
2240	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_0j",doPF,true);
2241	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_1j",doPF,true);
2242	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_2j",doPF,true);
2243	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_3j",doPF,true);
2244	+
2245	+
2246	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets),120,0,1200,dolog,"H_{T}","HT_3jets_OF",doPF,true);
2247	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF",doPF,true);
2248	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets&&TCut("met[4]>100")),120,0,1200,dolog,"H_{T}","HT_3jets_OF_MET100",doPF,true);
2249	+	make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets&&TCut("met[4]>100")) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF_MET100",doPF,true);
2250	+
2251	+
2252	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_0j",doPF,true);
2253	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_1j",doPF,true);
2254	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_2j",doPF,true);
2255	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_3j",doPF,true);
2256	+
2257	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_0j",doPF,true);
2258	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_1j",doPF,true);
2259	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_2j",doPF,true);
2260	+	make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_3j",doPF,true);
2261	+
2262	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70",doPF,true);
2263	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70",doPF,true);
2264	+
2265	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120",doPF,true);
2266	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120",doPF,true);
2267	+
2268	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70___2jetsonly",doPF,true);
2269	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70___2jetsonly",doPF,true);
2270	+
2271	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120___2jetsonly",doPF,true);
2272	+	make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120___2jetsonly",doPF,true);
2273	+
2274	+
2275	+	make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true);
2276	+	make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true);
2277		// Plots in signal region
2278	<	make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true);
2279	<	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true);
2280	<	make_kin_plot("mll","mll>20&&met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true);
2281	<	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true);
2282	<	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true);
2278	>	make_kin_plot("met[4]","",70,0,350,dolog,"#slash{E}_{T} [GeV]","met",doPF,true);
2279	>
2280	>
2281	>	make_kin_plot("MetFactor","",20,0,2,nolog,"MetFactor","MetFactor",doPF,true);
2282	>
2283	>	make_ttbar_comparison("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j__ALLMCSAMPLES");
2284	>	/*
2285	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_OF_30to70_0j");
2286	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_OF_30to70_1j");
2287	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j");
2288	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),dolog,"PFMET","Debunking/MET_OF_30to70_3j");
2289	>
2290	>	DoMCSystPlot("met[4]",40,0,400,(cutOSSF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_SF_30to70_0j");
2291	>	DoMCSystPlot("met[4]",40,0,400,(cutOSSF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_SF_30to70_1j");
2292	>	DoMCSystPlot("met[4]",40,0,400,(cutOSSF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_SF_30to70_2j");
2293	>	DoMCSystPlot("met[4]",40,0,400,(cutOSSF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),dolog,"PFMET","Debunking/MET_SF_30to70_3j");
2294	>
2295	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_OF_Z_0j");
2296	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_OF_Z_1j");
2297	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_Z_2j");
2298	>	DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),dolog,"PFMET","Debunking/MET_OF_Z_3j");
2299	>
2300	>
2301	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF,dolog,"m_{ll} [GeV]","kinSys/mll");
2302	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF&&TCut("id1==0"),dolog,"m_{ll} [GeV]","kinSys/mll_ee");
2303	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF&&TCut("id1==1"),dolog,"m_{ll} [GeV]","kinSys/mll_mm");
2304	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSOF,dolog,"m_{ll} [GeV]","kinSys/mll_osof");
2305	>
2306	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutOSSF,dolog,"m_{ll} [GeV]","kinSys/mll_inclusive");
2307	>	DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutOSOF,dolog,"m_{ll} [GeV]","kinSys/mll_inclusive_osof");
2308	>
2309	>	DoMCSystPlot("pfJetGoodNum40", 9,-0.5,8.5,cutOSSF,dolog,"nJets","kinSys/nJets_ossf");
2310	>	DoMCSystPlot("pfJetGoodNum40", 9,-0.5,8.5,cutOSOF,dolog,"nJets","kinSys/nJets_osof");
2311	>
2312	>	DoMCSystPlot("pfJetGoodPt[0]/pfJetGoodPt[1]",45,1,10,cutnJets&&cutOSSF,dolog,"p_{T}^{J1} / p_{T}^{J2}","kinSys/j1j2ratio");
2313	>	DoMCSystPlot("TMath::Abs(pfJetDphiMet[0])",32,0,3.2,cutnJets&&cutOSSF,dolog,"|#Delta#phi(jet1,MET)|","kinSys/dphiJ1MET");
2314	>	DoMCSystPlot(datajzb,mcjzb, (int)((PlottingSetup::jzbHigh+110)/10),-110,PlottingSetup::jzbHigh,cutnJets&&cutOSSF&&Restrmasscut,dolog,"m_{ll} [GeV]","kinSys/jzb_OS_SFZP");
2315	>	DoMCSystPlot(datajzb,mcjzb, (int)((PlottingSetup::jzbHigh+110)/10),-110,PlottingSetup::jzbHigh,cutnJets&&cutOSOF&&Restrmasscut,dolog,"m_{ll} [GeV]","kinSys/jzb_OS_OFZP");
2316	>
2317	>	make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true);
2318	>	make_kin_plot("mll","mll>20",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true);
2319	>	make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true);
2320	>	make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true);
2321	>
2322	>	make_kin_plot("pfJetGoodNum40",mllCut,9,-0.5,8.5,dolog,"nJets","nJets",doPF);
2323	>	make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF);
2324	>
2325	>	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_met100",doPF,true);
2326	>	make_kin_plot("mll","mll>20&&met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_met100",doPF,true,true);
2327	>	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee_met100",doPF,true);
2328	>	make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm_met100",doPF,true);
2329	>
2330
2331	<	make_kin_plot("pfJetGoodNum40","met[4]>100"+mllCut,9,-0.5,8.5,dolog,"nJets","nJets",doPF);
2332	<	make_kin_plot("pfJetGoodNum40","met[4]>100",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF);
2331	>	make_kin_plot("pfJetGoodNum40","met[4]>100"+mllCut,9,-0.5,8.5,dolog,"nJets","nJets_met100",doPF);
2332	>	make_kin_plot("pfJetGoodNum40","met[4]>100",9,-0.5,8.5,dolog,"nJets","nJets_osof_met100",doPF);
2333
2334		// Further inclusive invariant mass plots
2335	<	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
2335	>	make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
2336		make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true);
2337		make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true);
2338	<	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
2338	>	make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
2339
2340	<	//make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
2340	>	make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
2341		//if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]");
2342
826	–	//make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true);
827	–	//make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true);
2343
2344		// Number of jets
2345		make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_inclusive",doPF);
2346		make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof_inclusive",doPF);
2347	<	//make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF);
2347	>	make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF);
2348
2349		// Others
2350		make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF);
#	Line 861 | Line 2376 | void do_kinematic_plots(string mcjzb, st
2376
2377		//   make_kin_plot("mll","met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV] (MET>100GeV)","mll_met100_ll",doPF,true);
2378		//make_kin_plot("mll","met[4]>150&&id1==0",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ee} [GeV] (MET>150GeV)","mll_met150_ee",doPF,true);
2379	<	//make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);
2379	>	//make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true);*/
2380		}
2381
2382	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
2383	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
2384	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
2385	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,100);
2386	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,150);
2387	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,200);
2388	<	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,250);
2389	<
2390	<	//   make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
2391	<	//   make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
2392	<
2382	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
2383	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
2384	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
2385	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,100);
2386	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,150);
2387	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,200);
2388	>	make_special_obs_pred_mll_plot(datajzb,mcjzb,250);
2389	>	/*
2390	>	make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
2391	>	make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
2392	>	*/
2393	>	switch_overunderflow(false);
2394		}
2395
2396		void make_comp_plot( string var, string xlabel, string filename, float jzbcut, string mcjzb, string datajzb,
#	Line 914 | Line 2430 | void make_comp_plot( string var, string
2430		if (ymax>ymin) datahisto->SetMaximum(ymax);
2431		lm3histo->SetLineStyle(2);
2432		datahisto->Draw("e1");
2433	<	mcstack.Draw("same");
2433	>	mcstack.Draw("histo,same");
2434		datahisto->Draw("same,e1");
2435		lm3histo->Draw("hist,same");
2436		TLegend *kinleg = allsamples.allbglegend((sRegions[iregion]+(ijzb?"neg":"pos")).c_str());
#	Line 944 | Line 2460 | void region_comparison_plots(string mcjz
2460
2461		//  switch_overunderflow(true);  // switching overflow/underflow bins on
2462
2463	<
2463	>	switch_overunderflow(true);
2464		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- the arguments changed
2465		for(int ijzb=0;ijzb<(int)jzb_cuts.size();ijzb++) {
2466		float jzbcut=jzb_cuts[ijzb]; // Comparison plots are done for this JZB cut
#	Line 957 | Line 2473 | void region_comparison_plots(string mcjz
2473		make_comp_plot("TMath::Abs(pfJetDphiMet[0])","|#Delta#phi(jet1,MET)|","dphiJ1MET",jzbcut,mcjzb,datajzb,32,0,3.2,false,0,0,true);
2474
2475		make_comp_plot("mll","m_{ll} [GeV]","mll",jzbcut,mcjzb,datajzb,56,mll_low,mll_high,false,0,16.);
2476	<	make_comp_plot("met[4]","pfMET [GeV]","pfmet",jzbcut,mcjzb,datajzb,18,0,360,false,0,16.);
2476	>	make_comp_plot("met[4]","E_{T}^{miss} [GeV]","pfmet",jzbcut,mcjzb,datajzb,18,0,360,false,0,16.);
2477		make_comp_plot("pfJetGoodNum40","#(jets)","njets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.);
2478		make_comp_plot("pfJetGoodNumBtag","#(b-jets)","nBjets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.);
2479		make_comp_plot("pt","Z p_{T} [GeV]","Zpt",jzbcut,mcjzb,datajzb,26,0,525,false,0.,21.);
#	Line 989 | Line 2505 | void signal_bg_comparison()
2505
2506		float simulatedlumi=luminosity;//in pb please - adjust to your likings
2507
2508	<	TH1F *JZBplotZJETs = allsamples.Draw("JZBplotZJETs",jzbvariablemc,sbg_nbins,sbg_min,sbg_max, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/DY"));
2508	>	TH1F *JZBplotZJETs = allsamples.Draw("JZBplotZJETs",jzbvariablemc,sbg_nbins,sbg_min,sbg_max, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
2509		TH1F *JZBplotLM4;
2510		if(PlottingSetup::RestrictToMassPeak) JZBplotLM4 = allsamples.Draw("JZBplotLM4",jzbvariablemc,sbg_nbins,sbg_min,sbg_max, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("LM4"));
2511		else JZBplotLM4 = allsamples.Draw("JZBplotLM4",jzbvariablemc,sbg_nbins,sbg_min,sbg_max, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("LM3"));
#	Line 1261 | Line 2777 | if you want to start from scratch (witho
2777
2778		void do_prediction_plot(string jzb, TCanvas *globalcanvas, float high, int use_data, bool overlay_signal = false,string subdir="" )
2779		{
2780	<	//  switch_overunderflow(true);
2780	>
2781		bool is_data=false;
2782		bool use_signal=false;
2783		if(use_data==1) is_data=true;
2784		if(use_data==2) use_signal=true;
2785	<	int nbins=50;//100;
2786	<	if(is_data) nbins=50;
2785	>	int nbins=int(high/10);//100;
2786	>	if(is_data) nbins=int(high/10);
2787		float low=0;
2788	<	float hi=500;
2788	>	float hi=high;
2789	>
2790	>	stringstream cutpositiveS;
2791	>	cutpositiveS << "(" << jzb << ">0)";
2792	>	TCut cutpositive(cutpositiveS.str().c_str());
2793	>	stringstream cutnegativeS;
2794	>	cutnegativeS << "(" << jzb << "<0)";
2795	>	TCut cutnegative(cutnegativeS.str().c_str());
2796	>
2797
2798		TH1F *blankback = new TH1F("blankback","blankback",int(high/10),0,high);
2799	<	TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2800	<	TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2801	<	TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2802	<	TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2799	>	TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2800	>	TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2801	>	TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2802	>	TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2803
2804		blankback->GetXaxis()->SetTitle(RcorrJZBeemm->GetXaxis()->GetTitle());
2805		blankback->GetYaxis()->SetTitle(RcorrJZBeemm->GetYaxis()->GetTitle());
#	Line 1290 | Line 2814 | void do_prediction_plot(string jzb, TCan
2814
2815		TH1F *RcorrJZBeemmNoS;
2816
2817	<	//these are for the ratio
2818	<
2819	<	TH1F *JRcorrJZBeemm   = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2820	<	TH1F *JLcorrJZBeemm   = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2821	<	TH1F *JRcorrJZBem     = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
1298	<	TH1F *JLcorrJZBem     = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2817	>	//these are for the ratio
2818	>	TH1F *JRcorrJZBeemm   = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2819	>	TH1F *JLcorrJZBeemm   = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2820	>	TH1F *JRcorrJZBem     = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2821	>	TH1F *JLcorrJZBem     = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2822
2823		TH1F *JRcorrJZBSBem;
2824		TH1F *JLcorrJZBSBem;
2825		TH1F *JRcorrJZBSBeemm;
2826		TH1F *JLcorrJZBSBeemm;
2827
2828	<	if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false);
2828	>	if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false);
2829
2830
2831		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
2832	<	RcorrJZBSBem   = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2833	<	LcorrJZBSBem   = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2834	<	RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2835	<	LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2832	>	RcorrJZBSBem   = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2833	>	LcorrJZBSBem   = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2834	>	RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2835	>	LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2836
2837		//these are for the ratio
2838	<	JRcorrJZBSBem   = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2839	<	JLcorrJZBSBem   = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2840	<	JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2841	<	JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2838	>	JRcorrJZBSBem   = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2839	>	JLcorrJZBSBem   = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal);
2840	>	JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2841	>	JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal);
2842		}
2843
2844		TH1F *lm4RcorrJZBeemm;
2845	<	if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM"));
2845	>	if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM"));
2846
2847		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed.
2848
#	Line 1383 | Line 2906 | void do_prediction_plot(string jzb, TCan
2906		SQRT(BpredSys);
2907		BpredSys->Divide(JBpred);
2908
1386	–
2909		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed
2910		TH1F *Tpred = (TH1F*)RcorrJZBem->Clone("Bpred");
2911		Tpred->Add(LcorrJZBem,-1.0);
#	Line 1436 | Line 2958 | void do_prediction_plot(string jzb, TCan
2958		Bpred->Draw("hist,same");
2959		//analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same");
2960		RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP!
2961	<	lm4RcorrJZBeemm->Draw("hist,same");
2961	>	//lm4RcorrJZBeemm->Draw("hist,same");
2962		legBpred->AddEntry(RcorrJZBeemm,"observed","p");
2963		legBpred->AddEntry(Bpred,"predicted","l");
2964		//    legBpred->AddEntry(analytical_function[1],"predicted fit","l");
2965		//    legBpred->AddEntry(analytical_function[2],"stat. uncert.","l");
2966	<	legBpred->AddEntry(lm4RcorrJZBeemm,(allsamples.collection[allsamples.FindSample("LM")[0]].samplename).c_str(),"l");
2966	>	//     legBpred->AddEntry(lm4RcorrJZBeemm,(allsamples.collection[allsamples.FindSample("LM")[0]].samplename).c_str(),"l");
2967		legBpred->Draw();
2968		DrawPrelim();
2969
#	Line 1451 | Line 2973 | void do_prediction_plot(string jzb, TCan
2973		Bpred->SetLineWidth(2);
2974		predcomppad->cd();
2975		predcomppad->SetLogy(1);
1454	–
2976		TH1F *jzbnegative = (TH1F*)LcorrJZBeemm->Clone("jzbnegative");
2977		TH1F *sidebandsemu = (TH1F*)Bpred->Clone("sidebandsemu");
2978		sidebandsemu->Add(jzbnegative,-1);
#	Line 1473 | Line 2994 | void do_prediction_plot(string jzb, TCan
2994		//    lm4RcorrJZBeemm->SetLineColor(kOrange+1);
2995		lm4RcorrJZBeemm->SetLineWidth(2);
2996		//lm4RcorrJZBeemm->SetLineWidth(2); // paper style. overruled.
2997	<	lm4RcorrJZBeemm->Draw("histo,same");
2997	>	//     lm4RcorrJZBeemm->Draw("histo,same");
2998		DrawPrelim();
2999		TLegend *speciallegBpred = make_legend("",0.45,0.55);
3000		//TLegend *speciallegBpred = make_legend("",0.35,0.55); // paper style. overruled.
#	Line 1483 | Line 3004 | void do_prediction_plot(string jzb, TCan
3004		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) speciallegBpred->AddEntry(sidebandsemu,"Sidebands/e#mu (data)","f");
3005		else speciallegBpred->AddEntry(sidebandsemu,"e#mu (data)","f");
3006		//    speciallegBpred->AddEntry(lm4RcorrJZBeemmC,"LM4","l");
3007	<	speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3007	>	//     speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3008		speciallegBpred->Draw();
3009	<	save_with_ratio(JRcorrJZBeemm,JBpred,predcomppad,subdir+"Bpred_Data_____PredictionComposition",true,true,"data/pred",BpredSys);
3009	>	Save_With_Ratio(JRcorrJZBeemm,JBpred,predcomppad,subdir+"Bpred_Data_____PredictionComposition",true,true,"data/pred",BpredSys);
3010	>	delete predcomppad;
3011
3012		TCanvas *specialcanv = new TCanvas("specialcanv","specialcanv");
3013		specialcanv->SetLogy(1);
3014		//    THStack kostack = allsamples.DrawStack("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,!is_data, luminosity,use_signal);
3015		blankback->Draw();
3016		//    kostack.Draw("same");
3017	<	predcomposition.Draw();
3017	>	predcomposition.Draw("hist");
3018		Bpred->Draw("hist,same");
3019		//analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same");
3020		RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP!
#	Line 1519 | Line 3041 | void do_prediction_plot(string jzb, TCan
3041		legBpred->AddEntry(Bpred,"MC predicted","l");
3042		if(versok) legBpred->AddEntry((TObject*)0,"",""); // Just for alignment // causes seg fault on root v5.18
3043		if(versok) legBpred->AddEntry((TObject*)0,"",""); // causes seg fault on root v5.18
3044	<	if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3044	>	//     if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3045		legBpred->Draw();
3046		DrawMCPrelim();
3047		Bpredsaveas="Bpred_MC";
#	Line 1551 | Line 3073 | void do_prediction_plot(string jzb, TCan
3073
3074		RcorrJZBeemmNoS->SetLineStyle(2);
3075		legBpred2->AddEntry(RcorrJZBeemmNoS,"MC B","l");
3076	<	legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l");
3076	>	//       legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l");
3077		legBpred2->Draw();
3078		RcorrJZBeemmNoS->SetLineColor(TColor::GetColor("#61210B"));
3079		RcorrJZBeemmNoS->Draw("histo,same");
#	Line 1568 | Line 3090 | void do_prediction_plot(string jzb, TCan
3090		string ytitle("ratio");
3091		if ( use_data==1 ) ytitle = "data/pred";
3092		//save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,Bpredsaveas,true,use_data!=1,ytitle);
3093	<	save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,false,ytitle,BpredSys);//not extending the y range anymore up to 4
3093	>	Save_With_Ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,false,ytitle,BpredSys);//not extending the y range anymore up to 4
3094	>
3095	>	delete kinpad;
3096	>
3097
3098
3099		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
#	Line 1628 | Line 3153 | void do_prediction_plot(string jzb, TCan
3153		legBpredc->AddEntry(Bpredem,"MC OFZP","l");
3154		legBpredc->AddEntry(BpredSBem,"MC OFSB","l");
3155		legBpredc->AddEntry(BpredSBeemm,"MC SFSB","l");
3156	<	if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3156	>	//       if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l");
3157		legBpredc->Draw();
3158		CompleteSave(globalcanvas,subdir+"Bpred_MCwithS_comparison");
3159		}
#	Line 1664 | Line 3189 | void do_prediction_plot(string jzb, TCan
3189		delete JLcorrJZBSBeemm;
3190		}
3191		if(overlay_signal || use_data==2) delete lm4RcorrJZBeemm;
1667	–	switch_overunderflow(false);
3192		}
3193
3194		void do_prediction_plots(string mcjzb, string datajzb, float DataSigma, float MCSigma, bool overlay_signal ) {
3195	+	switch_overunderflow(true);
3196		TCanvas *globalcanvas = new TCanvas("globalcanvas","Prediction Canvas");
3197		do_prediction_plot(datajzb,globalcanvas,jzbHigh ,data,overlay_signal);
3198		if ( !PlottingSetup::Approved ) {
#	Line 1676 | Line 3201 | void do_prediction_plots(string mcjzb, s
3201		} else {
3202		write_info(__FUNCTION__,"You set approved to true, therefore not producing prediction/observation plots for MC with and without signal.");
3203		}
3204	<	}
3205	<
1681	<	void do_ratio_plot(int is_data,vector<float> binning, string jzb, TCanvas *can, float high=-9999) {
1682	<	bool do_data=0;
1683	<	bool dosignal=0;
1684	<	if(is_data==1) do_data=1;
1685	<	if(is_data==2) dosignal=1;
1686	<	TH1F *RcorrJZBeemm   = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1687	<	TH1F *LcorrJZBeemm   = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1688	<	TH1F *RcorrJZBem     = allsamples.Draw("RcorrJZBem",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1689	<	TH1F *LcorrJZBem     = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1690	<
1691	<	TH1F *RcorrJZBSBem;
1692	<	TH1F *LcorrJZBSBem;
1693	<	TH1F *RcorrJZBSBeemm;
1694	<	TH1F *LcorrJZBSbeemm;
1695	<
1696	<	if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
1697	<	RcorrJZBSBem   = allsamples.Draw("RcorrJZBSbem",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1698	<	LcorrJZBSBem   = allsamples.Draw("LcorrJZBSbem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal);
1699	<	RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSbeemm",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1700	<	LcorrJZBSbeemm = allsamples.Draw("LcorrJZBSbeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal);
1701	<	}
1702	<
1703	<
1704	<
1705	<
1706	<	flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
1707	<	TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
1708	<	if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
1709	<	Bpred->Add(RcorrJZBem,1.0/3);
1710	<	Bpred->Add(LcorrJZBem,-1.0/3);
1711	<	Bpred->Add(RcorrJZBSBem,1.0/3);
1712	<	Bpred->Add(LcorrJZBSBem,-1.0/3);
1713	<	Bpred->Add(RcorrJZBSBeemm,1.0/3);
1714	<	Bpred->Add(LcorrJZBSbeemm,-1.0/3);
1715	<	} else {
1716	<	Bpred->Add(RcorrJZBem,1.0);
1717	<	Bpred->Add(LcorrJZBem,-1.0);
1718	<	}
1719	<
1720	<	can->cd();
1721	<	can->SetLogy(0);
1722	<	Bpred->SetLineColor(kRed);
1723	<	Bpred->SetStats(0);
1724	<	if(high>0) Bpred->GetXaxis()->SetRangeUser(0,high);
1725	<	TH1F *JZBratioforfitting=(TH1F*)RcorrJZBeemm->Clone("JZBratioforfitting");
1726	<	JZBratioforfitting->Divide(Bpred);
1727	<	TGraphAsymmErrors *JZBratio = histRatio(RcorrJZBeemm,Bpred,is_data,binning,false);
1728	<
1729	<
1730	<	JZBratio->SetTitle("");
1731	<	JZBratio->GetYaxis()->SetRangeUser(0.0,9.0);
1732	<	//  if(is_data==1) JZBratio->GetXaxis()->SetRangeUser(0,jzbHigh );
1733	<
1734	<	TF1 *pol0 = new TF1("pol0","[0]",0,1000);
1735	<	TF1 *pol0d = new TF1("pol0","[0]",0,1000);
1736	<	//  straightline_fit->SetParameter(0,1);
1737	<	JZBratioforfitting->Fit(pol0,"Q0R","",0,30);
1738	<	pol0d->SetParameter(0,pol0->GetParameter(0));
1739	<
1740	<	JZBratio->GetYaxis()->SetTitle("Observed/Predicted");
1741	<	JZBratio->GetXaxis()->SetTitle("JZB [GeV]");
1742	<	if ( high>0 ) JZBratio->GetXaxis()->SetRangeUser(0.0,high);
1743	<	JZBratio->GetYaxis()->SetNdivisions(519);
1744	<	JZBratio->GetYaxis()->SetRangeUser(0.0,4.0);
1745	<	JZBratio->GetYaxis()->CenterTitle();
1746	<	JZBratio->GetXaxis()->CenterTitle();
1747	<	JZBratio->SetMarkerSize(DataMarkerSize);
1748	<	JZBratio->Draw("AP");
1749	<	/////----------------------------
1750	<	TPaveText *writeresult = new TPaveText(0.15,0.78,0.49,0.91,"blNDC");
1751	<	writeresult->SetFillStyle(4000);
1752	<	writeresult->SetFillColor(kWhite);
1753	<	writeresult->SetTextFont(42);
1754	<	writeresult->SetTextSize(0.03);
1755	<	writeresult->SetTextAlign(12);
1756	<	ostringstream jzb_agreement_data_text;
1757	<	jzb_agreement_data_text<< setprecision(2) << "mean =" << pol0->GetParameter(0) << " #pm " << setprecision(1) <<  pol0->GetParError(0);
1758	<	if(is_data==1) fitresultconstdata=pol0->GetParameter(0);// data
1759	<	if(is_data==0) fitresultconstmc=pol0->GetParameter(0); // monte carlo, no signal
1760	<	/*  if(is_data) writeresult->AddText("Data closure test");
1761	<	else writeresult->AddText("MC closure test");
1762	<	*/
1763	<	writeresult->AddText(jzb_agreement_data_text.str().c_str());
1764	<	//  writeresult->Draw("same");
1765	<	//  pol0d->Draw("same");
1766	<	TF1 *topline =  new TF1("","1.5",0,1000);
1767	<	TF1 *bottomline =  new TF1("","0.5",0,1000);
1768	<	topline->SetLineColor(kBlue);
1769	<	topline->SetLineStyle(2);
1770	<	bottomline->SetLineColor(kBlue);
1771	<	bottomline->SetLineStyle(2);
1772	<	//  topline->Draw("same");
1773	<	//  bottomline->Draw("same");
1774	<	TF1 *oneline = new TF1("","1.0",0,1000);
1775	<	oneline->SetLineColor(kBlue);
1776	<	oneline->SetLineStyle(1);
1777	<	oneline->Draw("same");
1778	<	TLegend *phony_leg = make_legend("ratio",0.6,0.55,false);//this line is just to have the default CMS Preliminary (...) on the canvas as well.
1779	<	if(is_data==1) DrawPrelim();
1780	<	else DrawMCPrelim();
1781	<	TLegend *leg = new TLegend(0.55,0.75,0.89,0.89);
1782	<	leg->SetTextFont(42);
1783	<	leg->SetTextSize(0.04);
1784	<	//   if(is_data==1) leg->SetHeader("Ratio (data)");
1785	<	//   else leg->SetHeader("Ratio (MC)");
1786	<
1787	<	TString MCtitle("MC ");
1788	<	if (is_data==1) MCtitle = "";
1789	<
1790	<	leg->SetFillStyle(4000);
1791	<	leg->SetFillColor(kWhite);
1792	<	leg->SetTextFont(42);
1793	<	//  leg->AddEntry(topline,"+20\% sys envelope","l");
1794	<	leg->AddEntry(JZBratio,MCtitle+"obs / "+MCtitle+"pred","p");
1795	<	leg->AddEntry(oneline,"ratio = 1","l");
1796	<	//  leg->AddEntry(pol0d,"fit in [0,30] GeV","l");
1797	<	//  leg->AddEntry(bottomline,"#pm50% envelope","l");
1798	<
1799	<
1800	<	//leg->Draw("same"); // no longer drawing legend
1801	<
1802	<	if(is_data==1) CompleteSave(can, "jzb_ratio_data");
1803	<	if(is_data==0) CompleteSave(can, "jzb_ratio_mc");
1804	<	if(is_data==2) CompleteSave(can, "jzb_ratio_mc_BandS");//special case, MC with signal!
1805	<
1806	<	delete RcorrJZBeemm;
1807	<	delete LcorrJZBeemm;
1808	<	delete RcorrJZBem;
1809	<	delete LcorrJZBem;
1810	<
1811	<	delete RcorrJZBSBem;
1812	<	delete LcorrJZBSBem;
1813	<	delete RcorrJZBSBeemm;
1814	<	delete LcorrJZBSbeemm;
1815	<	}
1816	<
1817	<	void do_ratio_plots(string mcjzb,string datajzb,vector<float> ratio_binning) {
1818	<	TCanvas *globalc = new TCanvas("globalc","Ratio Plot Canvas");
1819	<	globalc->SetLogy(0);
1820	<
1821	<	do_ratio_plot(mc,ratio_binning,mcjzb,globalc, jzbHigh  );
1822	<	do_ratio_plot(data,ratio_binning,datajzb,globalc, jzbHigh  );
1823	<	do_ratio_plot(mcwithsignal,ratio_binning,mcjzb,globalc, jzbHigh  );
3204	>	delete globalcanvas;
3205	>	switch_overunderflow(false);
3206		}
3207
3208		string give_jzb_expression(float peak, int type) {
#	Line 1843 | Line 3225 | void lepton_comparison_plots() {
3225		Float_t ymin = 1.e-5, ymax = 0.25;
3226		TCanvas *can = new TCanvas("can","Lepton Comparison Canvas");
3227		can->SetLogy(1);
3228	<	TH1F *eemc = allsamples.Draw("eemc","mll",50,50,150, "mll [GeV]", "events", cutOSSF&&cutnJets&&"(id1==0)",mc, luminosity,allsamples.FindSample("/DY"));
3229	<	TH1F *mmmc = allsamples.Draw("mmmc","mll",50,50,150, "mll [GeV]", "events", cutOSSF&&cutnJets&&"(id1==1)",mc, luminosity,allsamples.FindSample("/DY"));
3228	>	TH1F *eemc = allsamples.Draw("eemc","mll",50,50,150, "mll [GeV]", "events", cutOSSF&&cutnJets&&"(id1==0)",mc, luminosity,allsamples.FindSample("DYJetsToLL"));
3229	>	TH1F *mmmc = allsamples.Draw("mmmc","mll",50,50,150, "mll [GeV]", "events", cutOSSF&&cutnJets&&"(id1==1)",mc, luminosity,allsamples.FindSample("DYJetsToLL"));
3230		eemc->SetLineColor(kBlue);
3231		mmmc->SetLineColor(kRed);
3232		eemc->SetMinimum(0.1);
#	Line 1907 | Line 3289 | void lepton_comparison_plots() {
3289		eemmlegend->AddEntry(jmmd,"#mu#mu","l");
3290		eemmlegend->Draw();
3291		DrawPrelim();
3292	<	save_with_ratio(jeed,jmmd,eemmpad->cd(),"lepton_comparison/jzb_Comparing_ee_mm_data");
3292	>	Save_With_Ratio(jeed,jmmd,eemmpad->cd(),"lepton_comparison/jzb_Comparing_ee_mm_data");
3293
3294	<	TH1F *zjeed = allsamples.Draw("zjeed",jzbvariablemc,    int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&"(id1==0)",mc,  luminosity,allsamples.FindSample("/DY"));
3295	<	TH1F *zjmmd = allsamples.Draw("zjmmd",jzbvariablemc,    int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&"(id1==1)",mc,  luminosity,allsamples.FindSample("/DY"));
3296	<	TH1F *zjeemmd = allsamples.Draw("zjeemmd",jzbvariablemc,int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,            mc,  luminosity,allsamples.FindSample("/DY"));
3294	>	TH1F *zjeed = allsamples.Draw("zjeed",jzbvariablemc,    int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&"(id1==0)",mc,  luminosity,allsamples.FindSample("DYJetsToLL"));
3295	>	TH1F *zjmmd = allsamples.Draw("zjmmd",jzbvariablemc,    int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&"(id1==1)",mc,  luminosity,allsamples.FindSample("DYJetsToLL"));
3296	>	TH1F *zjeemmd = allsamples.Draw("zjeemmd",jzbvariablemc,int((jzbHigh+150)/5),-150,jzbHigh , "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,            mc,  luminosity,allsamples.FindSample("DYJetsToLL"));
3297		dout << "Z+Jets ee : " << zjeed->GetMean() << "+/-" << zjeed->GetMeanError() << endl;
3298		dout << "Z+Jets ee : " << zjmmd->GetMean() << "+/-" << zjmmd->GetMeanError() << endl;
3299		dout << "Z+Jets eemd : " << zjeemmd->GetMean() << "+/-" << zjeemmd->GetMeanError() << endl;
#	Line 2008 | Line 3390 | void draw_pure_jzb_histo(TCut cut,string
3390		datahisto->SetMinimum(0.1);
3391		datahisto->SetMarkerSize(DataMarkerSize);
3392		datahisto->Draw("e1");
3393	<	mcstack.Draw("same");
3393	>	mcstack.Draw("histo,same");
3394		datahisto->Draw("same,e1");
3395
3396		TLegend *leg;
#	Line 2030 | Line 3412 | void draw_pure_jzb_histo(TCut cut,string
3412		TText *writeline1 = write_cut_on_canvas(write_cut.c_str());
3413		writeline1->SetTextSize(0.035);
3414		writeline1->Draw();
3415	<	if(!Contains(savename,"Dibosons")) save_with_ratio(datahisto,mcstack,jzbpad->cd(),("jzb/"+savename));
3416	<	else save_with_ratio(datahisto,mcstack,jzbpad->cd(),savename);
3415	>	if(!Contains(savename,"Dibosons")) Save_With_Ratio(datahisto,mcstack,jzbpad->cd(),("jzb/"+savename));
3416	>	else Save_With_Ratio(datahisto,mcstack,jzbpad->cd(),savename);
3417		TPad *jzbpad2 = new TPad("jzbpad2","jzbpad2",0,0,1,1);
3418		jzbpad2->cd();
3419		jzbpad2->SetLogy(1);
#	Line 2039 | Line 3421 | void draw_pure_jzb_histo(TCut cut,string
3421		datahisto->SetMinimum(0.1);
3422		datahisto->SetMarkerSize(DataMarkerSize);
3423		datahisto->Draw("e1");
3424	<	mcstack.Draw("same");
3424	>	mcstack.Draw("histo,same");
3425		datahisto->Draw("same,e1");
3426		leg->SetHeader("");
3427		leg->Draw();
3428		writeline1->SetTextSize(0.035);
3429		writeline1->Draw();
3430		DrawPrelim();
3431	<	if(!Contains(savename,"Dibosons")) save_with_ratio(datahisto,mcstack,jzbpad2->cd(),("jzb/PositiveSideOnly/"+savename+""));
3432	<	else save_with_ratio(datahisto,mcstack,jzbpad2->cd(),(savename+"__PosOnly"));
3431	>	if(!Contains(savename,"Dibosons")) Save_With_Ratio(datahisto,mcstack,jzbpad2->cd(),("jzb/PositiveSideOnly/"+savename+""));
3432	>	else Save_With_Ratio(datahisto,mcstack,jzbpad2->cd(),(savename+"__PosOnly"));
3433		datahisto->Delete();
3434		mcstack.Delete();
3435		}
#	Line 2088 | Line 3470 | void produce_stretched_jzb_plots(string
3470
3471
3472		void diboson_plots(string mcjzb, string datajzb,vector<float> ratio_binning) {
3473	+	switch_overunderflow(true);
3474		vector<int> SamplesToBeModified = allsamples.FindSampleBySampleName("WW/WZ/ZZ");
3475
3476		if(SamplesToBeModified.size()==0 || SamplesToBeModified[0]==-1) {
#	Line 2127 | Line 3510 | void diboson_plots(string mcjzb, string
3510		dout << "     Reset xs for sample " << (allsamples.collection)[SamplesToBeModified[i]].samplename << " from " << Upxs << " to " << (allsamples.collection)[SamplesToBeModified[i]].xs << " (by a factor of " << stretchfactor << ") and reset the correct name (from " << Upname << ")" << endl;
3511
3512		}
3513	<
3513	>	//   switch_overunderflow(false);
3514		}
3515
3516
#	Line 2154 | Line 3537 | void draw_normalized_data_vs_data_histo(
3537		leg->AddEntry(datahisto2,legentry2.c_str());
3538		leg->Draw();
3539
3540	<	save_with_ratio(datahisto1,datahisto2,jzbpad->cd(),("jzb/"+savename));
3540	>	Save_With_Ratio(datahisto1,datahisto2,jzbpad->cd(),("jzb/"+savename));
3541
3542		datahisto1->Delete();
3543		datahisto2->Delete();
#	Line 2162 | Line 3545 | void draw_normalized_data_vs_data_histo(
3545
3546
3547		void jzb_plots(string mcjzb, string datajzb,vector<float> ratio_binning) {
3548	+	switch_overunderflow(true);
3549		TCanvas *can = new TCanvas("can","JZB Plots Canvas");
3550		float max=jzbHigh ;
3551		float min=-120;
#	Line 2198 | Line 3582 | void jzb_plots(string mcjzb, string data
3582		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB_coarse",can,coarse_binning);
3583		}
3584		delete can;
3585	+	switch_overunderflow(false);
3586		}
3587
3588
#	Line 2646 | Line 4031 | void find_correction_factors(string &jzb
4031
4032		}
4033
4034	<	void pick_up_events(string cut) {
4035	<	dout << "Picking up events with cut " << cut << endl;
2651	<	allsamples.PickUpEvents(cut);
4034	>	void pick_up_events(string cut, string filename, bool QuietMode=false) {
4035	>	allsamples.PickUpEvents(cut,filename,QuietMode);
4036		}
4037
4038		void save_template(string mcjzb, string datajzb,vector<float> jzb_cuts,float MCPeakError,float DataPeakError, vector<float> jzb_shape_limit_bins) {
#	Line 2738 | Line 4122 | void ttbar_sidebands_comparison(string m
4122		//in the case of the on peak analysis, we compare the 3 control regions to the real value
4123		//in the case of the OFF peak analysis, we compare our control region to the real value
4124		TCut weightbackup=cutWeight;
4125	+	switch_overunderflow(true);
4126
4127		bool doPURW=false;
4128
#	Line 2759 | Line 4144 | void ttbar_sidebands_comparison(string m
4144
4145		float simulatedlumi = luminosity; //in pb please - adjust to your likings
4146
4147	<	TH1F *TZem  = systsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4148	<	TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4147	>	TH1F *TZem  = allsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4148	>	TH1F *nTZem = allsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4149		TH1F *TSem;
4150		TH1F *nTSem;
4151	<	TH1F *TZeemm  = systsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4152	<	TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4151	>	TH1F *TZeemm  = allsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4152	>	TH1F *nTZeemm = allsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4153		TH1F *TSeemm;
4154		TH1F *nTSeemm;
4155
4156		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4157	<	TSem    = systsamples.Draw("TSem",       mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4158	<	nTSem   = systsamples.Draw("nTSem",  "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4159	<	TSeemm  = systsamples.Draw("TSeemm",     mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4160	<	nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
4157	>	TSem    = allsamples.Draw("TSem",       mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4158	>	nTSem   = allsamples.Draw("nTSem",  "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4159	>	TSeemm  = allsamples.Draw("TSeemm",     mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4160	>	nTSeemm = allsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,allsamples.FindSample("/TT"));
4161		}
4162
4163		TCanvas *tcan = new TCanvas("tcan","tcan");
#	Line 2804 | Line 4189 | void ttbar_sidebands_comparison(string m
4189		histos.push_back(TSem);
4190		histos.push_back(TSeemm);
4191		}
4192	<	draw_all_ttbar_histos(tcan,histos,"histo",1);
4192	>	draw_all_ttbar_histos(tcan,histos,"histo",8);
4193
4194		TLegend *leg = make_legend("MC t#bar{t}",0.6,0.65,false);
4195		leg->AddEntry(TZeemm,"SFZP","l");
#	Line 2816 | Line 4201 | void ttbar_sidebands_comparison(string m
4201		leg->Draw("same");
4202		DrawMCPrelim(simulatedlumi);
4203		CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison");
4204	+
4205		TH1F *TZemcopy = (TH1F*)TZem->Clone("TZemcopy");
4206		TH1F *TZeemmcopy = (TH1F*)TZeemm->Clone("TZeemmcopy");
4207		TH1F *TSeemmcopy;
#	Line 2881 | Line 4267 | void ttbar_sidebands_comparison(string m
4267		DrawMCPrelim(simulatedlumi);
4268		CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison_ratio");
4269
4270	+
4271		if (0) { // Turn this off: we don't need this
4272
4273		///-------------- second part: only look at the quantity we actually care about!
#	Line 2986 | Line 4373 | void ttbar_sidebands_comparison(string m
4373
4374		delete tcan;
4375		cutWeight=weightbackup;
4376	+	switch_overunderflow(false);
4377		}
4378
4379		void ttbar_sidebands_comparison(string mcjzb, vector<float> jzb_binning) {
#	Line 3020 | Line 4408 | void ttbar_sidebands_comparison(string m
4408		nicer_binning.push_back(150);
4409		//nicer_binning.push_back(175);
4410		nicer_binning.push_back(200);
4411	<	nicer_binning.push_back(250);
4412	<	nicer_binning.push_back(300);
4413	<	nicer_binning.push_back(400);
4411	>	//  nicer_binning.push_back(250);
4412	>	//  nicer_binning.push_back(300);
4413	>	//  nicer_binning.push_back(400);
4414
4415		ttbar_sidebands_comparison(mcjzb,nicer_binning, "ttbar/");
4416		}
4417
4418
4419	<	void zjets_prediction_comparison(string mcjzbWithPUa) {
4419	>	void zjets_prediction_comparison(string mcjzbWithPUa, TCut massregioncut, string massregionname) {
4420	>	cout << "****************************************************************************************************************************************************************" << endl;
4421		TCanvas *zcan = new TCanvas("zcan","zcan");
4422		//  zcan->SetLogy(1);
4423		TCut weightbackup=cutWeight;
4424
4425	<	bool UsePURW=false;
4425	>	bool UsePURW=true;
4426
4427
4428		string mcjzb;
4429		if(UsePURW) {
4430		mcjzb=mcjzbWithPUa;
4431	+	cout << "Using PURW peak positions" << endl;
4432		} else {
4433		// Do it without PU re-weighting
4434	+	cutWeight="1.0";
4435		float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0;
4436		stringstream resultsNoPU;
4437		stringstream noPUdatajzb;
#	Line 3051 | Line 4442 | void zjets_prediction_comparison(string
4442
4443		mcjzb = noPUmcjzb.str();
4444
3054	–	cutWeight="1.0";
4445		}
4446
4447
4448	+
4449		vector<float> binning;
4450		binning.push_back(0);
4451		binning.push_back(10);
#	Line 3072 | Line 4463 | void zjets_prediction_comparison(string
4463
4464		TCut kPos((mcjzb+">0").c_str());
4465		TCut kNeg((mcjzb+"<0").c_str());
4466	+	TCut reducedNJets(cutnJets);
4467		string var( "abs("+mcjzb+")" );
4468
4469	+
4470		TCut notTau("abs(genMID1)!=15");
4471	<	TCut ONLYTau("mll>20");
4471	>	TCut ee_mm_tautau("mll>0");
4472
4473
4474	<	TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
4475	<	TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
4474	>	TH1F *hJZBpos = allsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kPos&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4475	>	TH1F *hJZBneg = allsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kNeg&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4476
4477		hJZBpos->SetLineColor(kBlack);
4478		hJZBneg->SetLineColor(kRed);
#	Line 3094 | Line 4487 | void zjets_prediction_comparison(string
4487		leg->AddEntry(hJZBneg,"Predicted","l");
4488		leg->Draw("same");
4489		DrawMCPrelim(simulatedlumi);
4490	<	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction");
4490	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction"+any2string(massregionname));
4491
4492		TH1F* hratio = (TH1F*)hJZBpos->Clone("hratio");
4493		hratio->Divide(hJZBneg);
4494
4495		for(int i=1;i<=hJZBpos->GetNbinsX();i++) {
4496	<	cout << "Positive: " << hJZBpos->GetBinContent(i) << "   vs    Negative : " << hJZBneg->GetBinContent(i) << "   (ratio : " << hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i) << endl;
4496	>	cout << "Positive: " << hJZBpos->GetBinContent(i) << "   vs    Negative : " << hJZBneg->GetBinContent(i) << "   (ratio : " << hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i) << " +/- " << (hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i)) * sqrt(1.0/hJZBpos->GetBinContent(i) +  1.0/hJZBneg->GetBinContent(i)) << ")" << endl;
4497		}
4498
4499	<	zcan->SetLogy(0);
4500	<	hratio->GetYaxis()->SetRangeUser(0,2.0);
4499	>	//  zcan->SetLogy(0);
4500	>	hratio->GetYaxis()->SetRangeUser(0,2.5);
4501		hratio->GetYaxis()->SetTitle("Observed/Predicted");
4502		hratio->Draw("e1");
4503
4504	<	TLine *top = new TLine(binning[0],1.25,binning[binning.size()-1],1.25);
4504	>	TLine *top = new TLine(binning[0],1.2,binning[binning.size()-1],1.2);
4505		TLine *center = new TLine(binning[0],1.0,binning[binning.size()-1],1.0);
4506	<	TLine *bottom = new TLine(binning[0],0.75,binning[binning.size()-1],0.75);
4506	>	TLine *bottom = new TLine(binning[0],0.8,binning[binning.size()-1],0.8);
4507
4508
4509		top->SetLineColor(kBlue);top->SetLineStyle(2);
#	Line 3126 | Line 4519 | void zjets_prediction_comparison(string
4519		leg2->AddEntry(bottom,"syst. envelope","l");
4520		leg2->Draw("same");
4521		DrawMCPrelim(simulatedlumi);
4522	<	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio");
4522	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio"+any2string(massregionname));
4523
3131	–	TCut reducedNJets(cutnJets);
4524
4525	<	TH1F *TAUhJZBpos = systsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
4526	<	TH1F *LcorrJZBeemm     = systsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
4527	<	TH1F *RcorrJZBem       = systsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
4528	<	TH1F *LcorrJZBem       = systsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
4525	>	TH1F *TAUhJZBpos       = allsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4526	>	TH1F *LcorrJZBeemm     = allsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4527	>	TH1F *RcorrJZBem       = allsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4528	>	TH1F *LcorrJZBem       = allsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4529
4530		TH1F *RcorrJZBSBem;
4531		TH1F *LcorrJZBSBem;
#	Line 3141 | Line 4533 | void zjets_prediction_comparison(string
4533		TH1F *LcorrJZBSBeemm;
4534
4535		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4536	<	RcorrJZBSBem   = systsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
4537	<	LcorrJZBSBem   = systsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
4538	<	RcorrJZBSBeemm = systsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
4539	<	LcorrJZBSBeemm = systsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
4536	>	RcorrJZBSBem   = allsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4537	>	LcorrJZBSBem   = allsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4538	>	RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4539	>	LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4540		}
4541
4542		TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
#	Line 3175 | Line 4567 | void zjets_prediction_comparison(string
4567		TAUleg->AddEntry(Bpred,"Predicted","l");
4568		TAUleg->Draw("same");
4569		DrawMCPrelim(simulatedlumi);
4570	<	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction");
4570	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction__"+any2string(massregionname));
4571
4572		TH1F* TAUhratio = (TH1F*)TAUhJZBpos->Clone("TAUhratio");
4573		TAUhratio->Divide(Bpred);
#	Line 3193 | Line 4585 | void zjets_prediction_comparison(string
4585		center->Draw("same");
4586		bottom->Draw("same");
4587
4588	<	TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow #tau#tau",0.25,0.75,false);
4588	>	TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false);
4589		TAUleg2->AddEntry(TAUhratio,"obs / pred","pe");
4590		TAUleg2->AddEntry(bottom,"syst. envelope","l");
4591		TAUleg2->Draw("same");
4592		DrawMCPrelim(simulatedlumi);
4593	<	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio");
4593	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio"+any2string(massregionname));
4594
4595		delete Bpred;
4596		delete TAUhJZBpos;
4597		delete LcorrJZBeemm;
4598		delete RcorrJZBem;
4599		delete LcorrJZBem;
4600	+	delete hJZBpos;
4601	+	delete hJZBneg;
4602
4603		if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4604		delete RcorrJZBSBem;
#	Line 3213 | Line 4607 | void zjets_prediction_comparison(string
4607		delete LcorrJZBSBeemm;
4608		}
4609
4610	+	//***************************************************************************************************
4611	+
4612	+	string jzbvar=var;
4613	+	for(int i=0;i<60;i+=10) {
4614	+	stringstream sSpecialJZBCut;
4615	+	sSpecialJZBCut << jzbvar << ">" << i;
4616	+	TCut SpecialJZBCut(sSpecialJZBCut.str().c_str());
4617	+
4618	+	var="mll";
4619	+
4620	+	TH1F *FullJZBhJZBpos       = allsamples.Draw("FullJZBhJZBpos",var,10,0,200, "m_{ll} [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&notTau&&kPos&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4621	+	TH1F *LcorrJZBPluseemm     = allsamples.Draw("LcorrJZBPluseemm",var,10,0,200, "m_{ll} [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&notTau&&kNeg&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4622	+	TH1F *RcorrJZBPlusem       = allsamples.Draw("RcorrJZBPlusem",var,10,0,200, "m_{ll} [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&notTau&&kPos&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4623	+	TH1F *LcorrJZBPlusem       = allsamples.Draw("LcorrJZBPlusem",var,10,0,200, "m_{ll} [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&notTau&&kNeg&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4624	+
4625	+	TH1F *RcorrJZBPlusSBem;
4626	+	TH1F *LcorrJZBPlusSBem;
4627	+	TH1F *RcorrJZBPlusSBeemm;
4628	+	TH1F *LcorrJZBPlusSBeemm;
4629	+
4630	+
4631	+	if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4632	+	RcorrJZBPlusSBem   = allsamples.Draw("RcorrJZBPlusSBem",var,10,0,200, "m_{ll} [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&notTau&&kPos&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4633	+	LcorrJZBPlusSBem   = allsamples.Draw("LcorrJZBPlusSBem",var,10,0,200, "m_{ll} [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&notTau&&kNeg&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4634	+	RcorrJZBPlusSBeemm = allsamples.Draw("RcorrJZBPlusSBeemm",var,10,0,200, "m_{ll} [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&notTau&&kPos&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4635	+	LcorrJZBPlusSBeemm = allsamples.Draw("LcorrJZBPlusSBeemm",var,10,0,200, "m_{ll} [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&notTau&&kNeg&&massregioncut&&SpecialJZBCut,mc,simulatedlumi,allsamples.FindSample("DYJetsToLL"));
4636	+	}
4637	+
4638	+	TH1F *BpredPrime = (TH1F*)LcorrJZBPluseemm->Clone("BpredPrime");
4639	+	if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4640	+	BpredPrime->Add(RcorrJZBPlusem,1.0/3.);
4641	+	BpredPrime->Add(LcorrJZBPlusem,-1.0/3.);
4642	+	BpredPrime->Add(RcorrJZBPlusSBem,1.0/3.);
4643	+	BpredPrime->Add(LcorrJZBPlusSBem,-1.0/3.);
4644	+	BpredPrime->Add(RcorrJZBPlusSBeemm,1.0/3.);
4645	+	BpredPrime->Add(LcorrJZBPlusSBeemm,-1.0/3.);
4646	+	} else {
4647	+	BpredPrime->Add(RcorrJZBPlusem,1.0);
4648	+	BpredPrime->Add(LcorrJZBPlusem,-1.0);
4649	+	}
4650	+
4651	+	BpredPrime->SetLineColor(kRed);
4652	+
4653	+	FullJZBhJZBpos->SetLineColor(kBlack);
4654	+	BpredPrime->SetLineColor(kRed);
4655	+
4656	+	FullJZBhJZBpos->SetMinimum(1.0);
4657	+	FullJZBhJZBpos->Draw("e1");
4658	+	BpredPrime->Draw("same,hist");
4659	+	FullJZBhJZBpos->Draw("same,e1");
4660	+
4661	+	TLegend *FullJZBleg = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.55,0.75,false);
4662	+	FullJZBleg->AddEntry(FullJZBhJZBpos,"Observed","pe");
4663	+	FullJZBleg->AddEntry(BpredPrime,"Predicted","l");
4664	+	FullJZBleg->Draw("same");
4665	+	DrawMCPrelim(simulatedlumi);
4666	+	CompleteSave(zcan,"Systematics/ZJets/Mll_Distribution_JZB"+any2string(i)+"_zjets_eemumu_prediction__"+any2string(massregionname));
4667	+
4668	+	TH1F* FullJZBhratio = (TH1F*)FullJZBhJZBpos->Clone("FullJZBhratio");
4669	+	FullJZBhratio->Divide(BpredPrime);
4670	+
4671	+	zcan->SetLogy(0);
4672	+	FullJZBhratio->GetYaxis()->SetRangeUser(0,2.5);
4673	+	FullJZBhratio->GetYaxis()->SetTitle("Observed/Predicted");
4674	+	FullJZBhratio->Draw("e1");
4675	+
4676	+	TLine *atop = new TLine(0,1.2,200,1.2);
4677	+	TLine *acenter = new TLine(0,1.0,200,1.0);
4678	+	TLine *abottom = new TLine(0,0.8,200,0.8);
4679	+
4680	+	atop->SetLineColor(kBlue);
4681	+	atop->SetLineStyle(2);
4682	+	acenter->SetLineColor(kBlue);
4683	+	abottom->SetLineColor(kBlue);
4684	+	abottom->SetLineStyle(2);
4685	+
4686	+	atop->Draw("same");
4687	+	acenter->Draw("same");
4688	+	abottom->Draw("same");
4689	+
4690	+	TLegend *FullJZBleg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false);
4691	+	FullJZBleg2->AddEntry(FullJZBhratio,"obs / pred","pe");
4692	+	FullJZBleg2->AddEntry(abottom,"syst. envelope","l");
4693	+	FullJZBleg2->Draw("same");
4694	+	DrawMCPrelim(simulatedlumi);
4695	+	CompleteSave(zcan,"Systematics/ZJets/Mll_Distribution_JZB"+any2string(i)+"zjets_eemumu_prediction_ratio"+any2string(massregionname));
4696	+
4697	+
4698	+	delete BpredPrime;
4699	+	delete FullJZBhJZBpos;
4700	+	delete LcorrJZBPluseemm;
4701	+	delete RcorrJZBPlusem;
4702	+	delete LcorrJZBPlusem;
4703	+	//  delete FullJZBhJZBpos;
4704	+
4705	+	if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) {
4706	+	delete RcorrJZBPlusSBem;
4707	+	delete LcorrJZBPlusSBem;
4708	+	delete RcorrJZBPlusSBeemm;
4709	+	delete LcorrJZBPlusSBeemm;
4710	+	}
4711	+
4712	+	}
4713	+	//***************************************************************************************************
4714
4715		delete zcan;
4716		cutWeight=weightbackup;
3219	–
4717		}
4718
4719
4720	<
4720	>	void zjets_prediction_comparison(string mcjzbWithPUa) {
4721	>	zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>20"), "nomasscut");
4722	>	zjets_prediction_comparison(mcjzbWithPUa, TCut("abs(mll-91)<10"), "Zwindow_10");
4723	>	zjets_prediction_comparison(mcjzbWithPUa, TCut("abs(mll-91)>10&&mll>20"), "Outside_Zwindow_10");
4724	>	zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>20&&mll<70"), "LowMassRegion2070");
4725	>	zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>110"), "HighMassRegion110");
4726	>	}
4727	>
4728		void sideband_assessment(string datajzb, float min=30.0, float max=50.0) {
4729		tout << endl << endl;
4730		stringstream bordercut;
#	Line 3328 | Line 4832 | void make_table(samplecollection &coll,
4832		}
4833
4834		void met_jzb_cut(string datajzb, string mcjzb, vector<float> jzb_cut) {
4835	+	cout << "You probably don't want --met, you want --metplots ... " << endl;
4836	+	assert(0);
4837		/*we want a table like this:
4838		__________________     50   |   100  | ...
4839		| Data prediction |  a vs b | a vs b | ...
#	Line 3561 | Line 5067 | void qcd_plots(string datajzb, string mc
5067		//3rd last argument: do special bpred ratio, 2nd last argument: extended range!, last: y-axis title
5068		string ytitle("ratio");
5069		if ( use_data==1 ) ytitle = "data/pred";
5070	<	save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,"QCD/Bpred",true,false,ytitle);
5070	>	Save_With_Ratio(JRcorrJZBeemm,JBpred,kinpad,"QCD/Bpred",true,false,ytitle);
5071	>	delete kinpad;
5072
5073		TH1F *allevents = qcdsamples.Draw("allevents","pfJetGoodNum",1,0,100, "internal code", "events", "" ,mc, luminosity);
5074		TH1F *ossf = qcdsamples.Draw("ossf","pfJetGoodNum",1,0,100, "internal code", "events", cutOSSF ,mc, luminosity);
#	Line 3769 | Line 5276 | void met_vs_jzb_plots(string datajzb, st
5276		sRIGHT << "((" << mcjzb << ")>0)";
5277		TCut RIGHT(sRIGHT.str().c_str());
5278
5279	<	TH1F *metleft   = allsamples.Draw("metleft","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5280	<	TH1F *metleftO  = allsamples.Draw("metleftO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5281	<	TH1F *metright  = allsamples.Draw("metright","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5282	<	TH1F *metrightO = allsamples.Draw("metrightO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5279	>	TH1F *metleft   = allsamples.Draw("metleft","met[4]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5280	>	TH1F *metleftO  = allsamples.Draw("metleftO","met[4]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5281	>	TH1F *metright  = allsamples.Draw("metright","met[4]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5282	>	TH1F *metrightO = allsamples.Draw("metrightO","met[4]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5283
5284
5285		TH1F *Bpred =  (TH1F*)metleft->Clone("Bpred");
#	Line 3798 | Line 5305 | void met_vs_jzb_plots(string datajzb, st
5305		lg->SetHeader("DY");
5306
5307		lg->Draw();
5308	<	save_with_ratio(metright,metleft,metpad->cd(),"METvsJZBplots/ComparingLeftToRightinMETspectrum");
5308	>	Save_With_Ratio(metright,metleft,metpad->cd(),"METvsJZBplots/ComparingLeftToRightinMETspectrum");
5309
5310		TPad *metpad3 = new TPad("metpad3","metpad3",0,0,1,1);
5311		metpad3->cd();
#	Line 3816 | Line 5323 | void met_vs_jzb_plots(string datajzb, st
5323
5324		lg2->Draw();
5325
5326	<	save_with_ratio(obs,Bpred,metpad3->cd(),"METvsJZBplots/ComparingPredObsinMET");
5326	>	Save_With_Ratio(obs,Bpred,metpad3->cd(),"METvsJZBplots/ComparingPredObsinMET");
5327
5328		TPad *metpad2 = new TPad("metpad2","metpad2",0,0,1,1);
5329		metpad2->cd();
5330		metpad2->SetLogy(1);
5331
5332	<	TH1F *metlefta   = allsamples.Draw("metlefta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5333	<	TH1F *metleftOa  = allsamples.Draw("metleftOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5334	<	TH1F *metrighta  = allsamples.Draw("metrighta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5335	<	TH1F *metrightOa = allsamples.Draw("metrightOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5332	>	TH1F *metlefta   = allsamples.Draw("metlefta","met[2]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5333	>	TH1F *metleftOa  = allsamples.Draw("metleftOa","met[2]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
5334	>	TH1F *metrighta  = allsamples.Draw("metrighta","met[2]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5335	>	TH1F *metrightOa = allsamples.Draw("metrightOa","met[2]",nbins,low,high, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
5336
5337		metlefta->Add(metleftOa,-1);
5338		metrighta->Add(metrightOa,-1);
#	Line 3835 | Line 5342 | void met_vs_jzb_plots(string datajzb, st
5342		metlefta->Draw("histo");
5343		metrighta->Draw("same");
5344		lg->Draw();
5345	<	save_with_ratio(metrighta,metlefta,metpad2->cd(),"METvsJZBplots/ComparingLeftToRightinMET_type1_spectrum");
5345	>	Save_With_Ratio(metrighta,metlefta,metpad2->cd(),"METvsJZBplots/ComparingLeftToRightinMET_type1_spectrum");
5346
5347		delete Bpred;
5348		delete obs;
#	Line 3844 | Line 5351 | void met_vs_jzb_plots(string datajzb, st
5351		int newNBins=30;
5352
5353		TPad *metpad4 = new TPad("metpad4","metpad4",0,0,1,1);
5354	<	TH1F *Ametleft   = allsamples.Draw("Ametleft","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity);
5355	<	TH1F *AmetleftO  = allsamples.Draw("AmetleftO","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity);
5356	<	TH1F *Ametright  = allsamples.Draw("Ametright","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity);
5357	<	TH1F *AmetrightO = allsamples.Draw("AmetrightO","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity);
5354	>	TH1F *Ametleft   = allsamples.Draw("Ametleft","met[4]",newNBins,low,newhigh, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity);
5355	>	TH1F *AmetleftO  = allsamples.Draw("AmetleftO","met[4]",newNBins,low,newhigh, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity);
5356	>	TH1F *Ametright  = allsamples.Draw("Ametright","met[4]",newNBins,low,newhigh, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity);
5357	>	TH1F *AmetrightO = allsamples.Draw("AmetrightO","met[4]",newNBins,low,newhigh, "#slash{E}_{T} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity);
5358
5359		TH1F *aBpred = (TH1F*)Ametleft->Clone("aBpred");
5360		aBpred->Add(AmetleftO,-1);
#	Line 3862 | Line 5369 | void met_vs_jzb_plots(string datajzb, st
5369		aobs->Draw("same");
5370		lg->SetHeader("All MC");
5371		lg->Draw();
5372	<	save_with_ratio(aobs,aBpred,metpad4->cd(),"METvsJZBplots/ComparingPredObsinMET_ALLSAMPLES");
5372	>	Save_With_Ratio(aobs,aBpred,metpad4->cd(),"METvsJZBplots/ComparingPredObsinMET_ALLSAMPLES");
5373
5374
5375		delete lg;
#	Line 3872 | Line 5379 | void met_vs_jzb_plots(string datajzb, st
5379		delete metright;
5380		delete metrightO;
5381		}
5382	+
5383	+
5384	+	void do_one_ttbar_test(TH1F* &observed, TH1F* &predicted, TH1F* &ratio, TH1F* &ratio2, string mcjzb, string prestring) {
5385	+
5386	+	if(PlottingSetup::RestrictToMassPeak) {
5387	+	write_error(__FUNCTION__,"This function (ttbar_tests) was not written for on-peak studies - sorry.");
5388	+	assert(!PlottingSetup::RestrictToMassPeak);
5389	+	}
5390	+
5391	+	vector<float> binning;
5392	+
5393	+	binning.push_back(-200);
5394	+	binning.push_back(-150);
5395	+	binning.push_back(-125);
5396	+	binning.push_back(-100);
5397	+	binning.push_back(-75);
5398	+	binning.push_back(-50);
5399	+	binning.push_back(-25);
5400	+	binning.push_back(0);
5401	+	binning.push_back(25);
5402	+	binning.push_back(50);
5403	+	binning.push_back(75);
5404	+	binning.push_back(100);
5405	+	binning.push_back(125);
5406	+	binning.push_back(150);
5407	+	binning.push_back(200);
5408	+
5409	+	switch_overunderflow(true);
5410	+
5411	+	TH1F *TZem  = allsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,luminosity,allsamples.FindSample("/TT"));
5412	+	TH1F *nTZem = allsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,luminosity,allsamples.FindSample("/TT"));
5413	+	TH1F *TZeemm  = allsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,luminosity,allsamples.FindSample("/TT"));
5414	+	TH1F *nTZeemm = allsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,luminosity,allsamples.FindSample("/TT"));
5415	+
5416	+	TCanvas *tcan = new TCanvas("tcan","tcan");
5417	+	tcan->SetLogy(1);
5418	+
5419	+	TZeemm->SetLineColor(kBlack);
5420	+	TZem->SetLineColor(kRed);
5421	+	TZeemm->SetMarkerColor(kBlack);
5422	+	TZem->SetMarkerColor(kRed);
5423	+
5424	+	vector<TH1F*> histos;
5425	+	//  TZem->GetXaxis()->SetRangeUser(-100,binning[binning.size()-1]);
5426	+	//  TZeemm->GetXaxis()->SetRangeUser(-100,binning[binning.size()-1]);
5427	+	histos.push_back(TZem);
5428	+	histos.push_back(TZeemm);
5429	+	draw_all_ttbar_histos(tcan,histos,"histo",8);
5430	+
5431	+	TLegend *leg = make_legend("MC t#bar{t}",0.6,0.65,false);
5432	+	leg->AddEntry(TZeemm,"SFZP","l");
5433	+	leg->AddEntry(TZem,"OFZP","l");
5434	+	leg->Draw("same");
5435	+	DrawMCPrelim();
5436	+	CompleteSave(tcan,"Systematics/TtbarTests/"+prestring+"/ttbar_shape_comparison");
5437	+
5438	+	TH1F *TZemcopy = (TH1F*)TZem->Clone("TZemcopy");
5439	+	TH1F *TZeemmcopy = (TH1F*)TZeemm->Clone("TZeemmcopy");
5440	+	TH1F *TSeemmcopy;
5441	+	TH1F *TSemcopy;
5442	+
5443	+	predicted=(TH1F*)TZem->Clone(((string)"predicted_"+prestring).c_str());
5444	+	observed=(TH1F*)TZeemm->Clone(((string)"observed_"+prestring).c_str());
5445	+
5446	+	vector<float> posbinning;
5447	+	for(unsigned int i=0;i<binning.size();i++) {
5448	+	if(binning[i]<0) continue;
5449	+	posbinning.push_back(binning[i]);
5450	+	}
5451	+
5452	+	TH1F *pred2 = new TH1F("pred2","pred2",posbinning.size()-1,&posbinning[0]);pred2->Sumw2();
5453	+	TH1F *obs2 = new TH1F("obs2","obs2",posbinning.size()-1,&posbinning[0]);obs2->Sumw2();
5454	+
5455	+	for(unsigned int i=1;i<=predicted->GetNbinsX();i++) {
5456	+	int index=pred2->FindBin(abs(TZem->GetBinCenter(i)));
5457	+	if(TZem->GetBinCenter(i)<0) {
5458	+	//we're on the left: subtract em!
5459	+	pred2->SetBinContent(index,pred2->GetBinContent(index)-TZem->GetBinContent(i));
5460	+	pred2->SetBinError(index,sqrt(pow(pred2->GetBinError(index),2)+pow(TZem->GetBinError(i),2)));
5461	+	//we're on the left: add eemm!
5462	+	pred2->SetBinContent(index,pred2->GetBinContent(index)+TZeemm->GetBinContent(i));
5463	+	pred2->SetBinError(index,sqrt(pow(pred2->GetBinError(index),2)+pow(TZeemm->GetBinError(i),2)));
5464	+	} else {
5465	+	//we're on the right: add em!
5466	+	pred2->SetBinContent(index,pred2->GetBinContent(index)+TZem->GetBinContent(i));
5467	+	pred2->SetBinError(index,sqrt(pow(pred2->GetBinError(index),2)+pow(TZem->GetBinError(i),2)));
5468	+	//we're on the left: add eemm!
5469	+	obs2->SetBinContent(index,obs2->GetBinContent(index)+TZeemm->GetBinContent(i));
5470	+	obs2->SetBinError(index,sqrt(pow(pred2->GetBinError(index),2)+pow(TZeemm->GetBinError(i),2)));
5471	+	}
5472	+	}
5473	+
5474	+	ratio2 = (TH1F*)pred2->Clone(((string)"ratio2_"+prestring).c_str());
5475	+	ratio2->Divide(obs2);
5476	+
5477	+	TZem->Divide(TZeemm);
5478	+	TZem->SetMarkerStyle(21);
5479	+
5480	+	TZem->SetTitle("OF / SF");
5481	+	ratio = (TH1F*)TZem->Clone(((string)"ratio_"+prestring).c_str());
5482	+
5483	+	tcan->SetLogy(0);
5484	+	TZem->GetYaxis()->SetRangeUser(0,2.5);
5485	+	TZem->GetYaxis()->SetTitle("ratio");
5486	+	TZem->Draw();
5487	+
5488	+	float linepos=emuncertOFFPEAK;
5489	+
5490	+	TLine *top = new TLine(binning[0],1.0+linepos,binning[binning.size()-1],1.0+linepos);
5491	+	TLine *center = new TLine(binning[0],1.0,binning[binning.size()-1],1.0);
5492	+	TLine *bottom = new TLine(binning[0],1.0-linepos,binning[binning.size()-1],1.0-linepos);
5493	+
5494	+	top->SetLineColor(kBlue);top->SetLineStyle(2);
5495	+	bottom->SetLineColor(kBlue);bottom->SetLineStyle(2);
5496	+	center->SetLineColor(kBlue);
5497	+
5498	+	top->Draw("same");
5499	+	center->Draw("same");
5500	+	bottom->Draw("same");
5501	+
5502	+	TLegend *leg2 = make_legend("MC t#bar{t}",0.55,0.75,false);
5503	+	leg2->AddEntry(TZem,"OFZP / SFZP","ple");
5504	+	leg2->AddEntry(bottom,"syst. envelope","l");
5505	+	leg2->SetX1(0.25);leg2->SetX2(0.6);
5506	+	leg2->SetY1(0.65);
5507	+
5508	+	leg2->Draw("same");
5509	+
5510	+	DrawMCPrelim();
5511	+	CompleteSave(tcan,"Systematics/TtbarTests/"+prestring+"/ttbar_shape_comparison_ratio");
5512	+
5513	+	ratio2->SetLineColor(TZem->GetLineColor());
5514	+	ratio2->SetMarkerColor(TZem->GetMarkerColor());
5515	+	ratio2->GetYaxis()->SetRangeUser(0,2.5);
5516	+	ratio2->GetYaxis()->SetTitle("ratio");
5517	+	ratio2->Draw();
5518	+
5519	+	TLine *top2 = new TLine(0,1.0+linepos,binning[binning.size()-1],1.0+linepos);
5520	+	TLine *center2 = new TLine(0,1.0,binning[binning.size()-1],1.0);
5521	+	TLine *bottom2 = new TLine(0,1.0-linepos,binning[binning.size()-1],1.0-linepos);
5522	+
5523	+	top2->SetLineColor(kBlue);top->SetLineStyle(2);
5524	+	bottom2->SetLineColor(kBlue);bottom->SetLineStyle(2);
5525	+	center2->SetLineColor(kBlue);
5526	+
5527	+	top2->Draw("same");
5528	+	center2->Draw("same");
5529	+	bottom2->Draw("same");
5530	+
5531	+	leg2->Draw("same");
5532	+
5533	+	DrawMCPrelim();
5534	+	CompleteSave(tcan,"Systematics/TtbarTests/"+prestring+"/ttbar_shape_comparison_FINALratio");
5535	+
5536	+	delete TZem;
5537	+	delete nTZem;
5538	+	delete TZeemm;
5539	+	delete nTZeemm;
5540	+	delete obs2;
5541	+	delete pred2;
5542	+	delete tcan;
5543	+	switch_overunderflow(false);
5544	+	}
5545	+
5546	+	void IllustrateVariation(TH1F *vanilla_ratio, TH1F *response_only_up_ratio, TH1F *response_only_down_ratio, TH1F *reponse_plus_peak_up_ratio, TH1F *reponse_plus_peak_down_ratio, string PostString) {
5547
5548	+
5549	+	TCanvas *ttcan2 = new TCanvas("ttcan2","ttcan2");
5550	+	vanilla_ratio->GetYaxis()->SetRangeUser(0.6,1.4);
5551	+	vanilla_ratio->SetFillColor(kGreen);
5552	+	vanilla_ratio->SetLineColor(kGreen);
5553	+	vanilla_ratio->SetMarkerColor(kGreen);
5554	+	vanilla_ratio->SetFillStyle(3001);
5555	+	vanilla_ratio->Draw("e3");
5556	+
5557	+	response_only_up_ratio->SetFillColor(kRed);
5558	+	response_only_up_ratio->SetMarkerColor(kRed);
5559	+	response_only_up_ratio->SetLineColor(kRed);
5560	+	response_only_up_ratio->SetFillStyle(3003);
5561	+	response_only_up_ratio->Draw("e3,same");
5562	+
5563	+	response_only_down_ratio->SetFillColor(kOrange);
5564	+	response_only_down_ratio->SetMarkerColor(kOrange);
5565	+	response_only_down_ratio->SetLineColor(kOrange);
5566	+	response_only_down_ratio->SetFillStyle(3004);
5567	+	response_only_down_ratio->Draw("e3,same");
5568	+
5569	+	reponse_plus_peak_up_ratio->SetFillColor(TColor::GetColor("#BDBDBD"));
5570	+	reponse_plus_peak_up_ratio->SetMarkerColor(TColor::GetColor("#BDBDBD"));
5571	+	reponse_plus_peak_up_ratio->SetLineColor(TColor::GetColor("#BDBDBD"));
5572	+	reponse_plus_peak_up_ratio->SetFillStyle(3005);
5573	+	reponse_plus_peak_up_ratio->Draw("e3,same");
5574	+
5575	+	reponse_plus_peak_down_ratio->SetFillColor(kBlue);
5576	+	reponse_plus_peak_down_ratio->SetMarkerColor(kBlue);
5577	+	reponse_plus_peak_down_ratio->SetLineColor(kBlue);
5578	+	reponse_plus_peak_down_ratio->SetFillStyle(3006);
5579	+	reponse_plus_peak_down_ratio->Draw("e3,same");
5580	+
5581	+	TLine *llow = new TLine(vanilla_ratio->GetBinLowEdge(1),1.0-emuncertOFFPEAK,vanilla_ratio->GetBinLowEdge(vanilla_ratio->GetNbinsX())+vanilla_ratio->GetBinWidth(vanilla_ratio->GetNbinsX()),1.0-emuncertOFFPEAK);
5582	+	TLine *lhi = new TLine(vanilla_ratio->GetBinLowEdge(1),1.0+emuncertOFFPEAK,vanilla_ratio->GetBinLowEdge(vanilla_ratio->GetNbinsX())+vanilla_ratio->GetBinWidth(vanilla_ratio->GetNbinsX()),1.0+emuncertOFFPEAK);
5583	+	TLine *lOne = new TLine(vanilla_ratio->GetBinLowEdge(1),1.0,vanilla_ratio->GetBinLowEdge(vanilla_ratio->GetNbinsX())+vanilla_ratio->GetBinWidth(vanilla_ratio->GetNbinsX()),1.0);
5584	+
5585	+	TLine *llowA = new TLine(vanilla_ratio->GetBinLowEdge(1),1.0-0.07,vanilla_ratio->GetBinLowEdge(vanilla_ratio->GetNbinsX())+vanilla_ratio->GetBinWidth(vanilla_ratio->GetNbinsX()),1.0-0.07);
5586	+	TLine *lhiA = new TLine(vanilla_ratio->GetBinLowEdge(1),1.0+0.07,vanilla_ratio->GetBinLowEdge(vanilla_ratio->GetNbinsX())+vanilla_ratio->GetBinWidth(vanilla_ratio->GetNbinsX()),1.0+0.07);
5587	+
5588	+	lhiA->SetLineColor(kRed);
5589	+	lhiA->SetLineStyle(2);
5590	+	llowA->SetLineColor(kRed);
5591	+	llowA->SetLineStyle(2);
5592	+
5593	+	llow->SetLineColor(kBlue);
5594	+	llow->SetLineStyle(2);
5595	+	lhi->SetLineColor(kBlue);
5596	+	lhi->SetLineStyle(2);
5597	+	lOne->SetLineColor(kBlue);
5598	+
5599	+	llow->Draw();
5600	+	lhi->Draw();
5601	+	lOne->Draw();
5602	+
5603	+	llowA->Draw();
5604	+	lhiA->Draw();
5605	+
5606	+	TLegend *leg = make_legend();
5607	+	leg->SetX1(0.15);
5608	+	leg->SetY1(0.7);
5609	+	leg->AddEntry(vanilla_ratio,"Raw JZB","pf");
5610	+	leg->AddEntry(response_only_up_ratio,"With response correction","pf");
5611	+	leg->AddEntry(response_only_down_ratio,"With inverted response correction","pf");
5612	+	leg->AddEntry(reponse_plus_peak_up_ratio,"With response and peak correction","pf");
5613	+	leg->AddEntry(reponse_plus_peak_down_ratio,"With response and inverted peak corr.","pf");
5614	+	leg->AddEntry(lhi,"Systematic envelope","l");
5615	+	//  leg->SetNColumns(2);
5616	+	leg->Draw();
5617	+
5618	+	CompleteSave(ttcan2,"Systematics/TtbarTests/Comparison"+PostString);
5619	+
5620	+	vanilla_ratio->GetYaxis()->SetRangeUser(0.6,1.4);
5621	+	vanilla_ratio->Draw("e3");
5622	+	response_only_up_ratio->Draw("e3 same");
5623	+	response_only_down_ratio->Draw("e3 same");
5624	+	TLegend *leg2 = make_legend();
5625	+	leg2->SetX1(0.15);
5626	+	leg2->SetY1(0.7);
5627	+	leg2->AddEntry(vanilla_ratio,"Raw JZB","pf");
5628	+	leg2->AddEntry(response_only_up_ratio,"With response correction","pf");
5629	+	leg2->AddEntry(response_only_down_ratio,"With inverted response correction","pf");
5630	+	leg2->AddEntry(lhi,"Systematic envelope","l");
5631	+	leg2->Draw();
5632	+	llow->Draw();
5633	+	lOne->Draw();
5634	+	llowA->Draw();
5635	+	lhiA->Draw();
5636	+
5637	+	CompleteSave(ttcan2,"Systematics/TtbarTests/Comparison_Response"+PostString);
5638	+
5639	+	response_only_up_ratio->GetYaxis()->SetRangeUser(0.6,1.4);
5640	+	response_only_up_ratio->Draw("e3");
5641	+	reponse_plus_peak_up_ratio->Draw("e3 same");
5642	+	reponse_plus_peak_down_ratio->Draw("e3 same");
5643	+	TLegend *leg3 = make_legend();
5644	+	leg3->SetX1(0.15);
5645	+	leg3->SetY1(0.7);
5646	+	leg3->AddEntry(response_only_up_ratio,"With response correction","pf");
5647	+	leg3->AddEntry(reponse_plus_peak_up_ratio,"With response and peak correction","pf");
5648	+	leg3->AddEntry(reponse_plus_peak_down_ratio,"With response and inverted peak corr.","pf");
5649	+	leg3->AddEntry(lhi,"Systematic envelope","l");
5650	+	leg3->Draw();
5651	+	llow->Draw();
5652	+	lOne->Draw();
5653	+	llowA->Draw();
5654	+	lhiA->Draw();
5655	+
5656	+	CompleteSave(ttcan2,"Systematics/TtbarTests/Comparison_Peak"+PostString);
5657	+
5658	+	delete ttcan2;
5659	+	}
5660	+
5661	+	void ttbar_correction_tests() {
5662	+	TCanvas *ttcan = new TCanvas("ttcan","ttcan");
5663	+
5664	+	TH1F *vanilla_observed, *response_only_up_observed, *response_only_down_observed, *reponse_plus_peak_up_observed, *reponse_plus_peak_down_observed;
5665	+
5666	+	TH1F *vanilla_predicted, *response_only_up_predicted, *response_only_down_predicted, *reponse_plus_peak_up_predicted, *reponse_plus_peak_down_predicted;
5667	+
5668	+	TH1F *vanilla_ratio, *response_only_up_ratio, *response_only_down_ratio, *reponse_plus_peak_up_ratio, *reponse_plus_peak_down_ratio;
5669	+
5670	+	TH1F *vanilla_ratio2, *response_only_up_ratio2, *response_only_down_ratio2, *reponse_plus_peak_up_ratio2, *reponse_plus_peak_down_ratio2;
5671	+
5672	+	string vanilla="jzb[1]";
5673	+	string response_only_up="((jzb[1]+0.034665*pt))";
5674	+	string response_only_down="((jzb[1]-0.034665*pt))";
5675	+	string reponse_plus_peak_up="((jzb[1]+0.034665*pt)- 3.58273 )";
5676	+	string reponse_plus_peak_down="((jzb[1]+0.034665*pt)+ 3.58273 )";
5677	+
5678	+	do_one_ttbar_test(vanilla_observed,vanilla_predicted,vanilla_ratio,vanilla_ratio2,vanilla,"vanilla");
5679	+	do_one_ttbar_test(response_only_up_observed,response_only_up_predicted,response_only_up_ratio,response_only_up_ratio2,response_only_up,"response_only_up");
5680	+	do_one_ttbar_test(response_only_down_observed,response_only_down_predicted,response_only_down_ratio,response_only_down_ratio2,response_only_down,"response_only_down");
5681	+	do_one_ttbar_test(reponse_plus_peak_up_observed,reponse_plus_peak_up_predicted,reponse_plus_peak_up_ratio,reponse_plus_peak_up_ratio2,reponse_plus_peak_up,"reponse_plus_peak_up");
5682	+	do_one_ttbar_test(reponse_plus_peak_down_observed,reponse_plus_peak_down_predicted,reponse_plus_peak_down_ratio,reponse_plus_peak_down_ratio2,reponse_plus_peak_down,"reponse_plus_peak_down");
5683	+
5684	+	ttcan->cd();
5685	+
5686	+	IllustrateVariation(vanilla_ratio,response_only_up_ratio,response_only_down_ratio,reponse_plus_peak_up_ratio,reponse_plus_peak_down_ratio,"PartialRatio");
5687	+	IllustrateVariation(vanilla_ratio2,response_only_up_ratio2,response_only_down_ratio2,reponse_plus_peak_up_ratio2,reponse_plus_peak_down_ratio2,"FullRatio");
5688	+
5689	+	delete vanilla_observed;
5690	+	delete response_only_up_observed;
5691	+	delete response_only_down_observed;
5692	+	delete reponse_plus_peak_up_observed;
5693	+	delete reponse_plus_peak_down_observed;
5694	+
5695	+	delete vanilla_predicted;
5696	+	delete response_only_up_predicted;
5697	+	delete response_only_down_predicted;
5698	+	delete reponse_plus_peak_up_predicted;
5699	+	delete reponse_plus_peak_down_predicted;
5700	+
5701	+	delete vanilla_ratio;
5702	+	delete response_only_up_ratio;
5703	+	delete response_only_down_ratio;
5704	+	delete reponse_plus_peak_up_ratio;
5705	+	delete reponse_plus_peak_down_ratio;
5706	+
5707	+	delete ttcan;
5708	+	}
5709	+
5710	+	void ttbar_region_search(string mcjzb) {
5711	+	cout << "Looking for a nice control region" << endl;
5712	+	TCanvas *can = new TCanvas("can","can",1200,1200);
5713	+	can->Divide(2,2);
5714	+	can->cd(1);
5715	+	can->cd(1)->SetLogy(1);
5716	+
5717	+	THStack JZBdistSF = allsamples.DrawStack("JZBdistSF",mcjzb,100,-150,150, "JZB [GeV]", "events", cutmass&&cutOSSF&&TCut("pfJetGoodNum40==2"),mc, luminosity);
5718	+	THStack JZBdistOF = allsamples.DrawStack("JZBdistOF",mcjzb,100,-150,150, "JZB [GeV]", "events", cutmass&&cutOSOF&&TCut("pfJetGoodNum40==2"),mc, luminosity);
5719	+
5720	+	JZBdistSF.Draw("histo");
5721	+	JZBdistSF.SetMaximum(1000);
5722	+	JZBdistSF.SetMinimum(0.1);
5723	+	JZBdistSF.Draw("histo");
5724	+	DrawMCPrelim();
5725	+
5726	+	can->cd(2);
5727	+	can->cd(2)->SetLogy(1);
5728	+	JZBdistOF.Draw("histo");
5729	+	JZBdistOF.SetMaximum(1000);
5730	+	JZBdistOF.SetMinimum(0.1);
5731	+	JZBdistOF.Draw("histo");
5732	+	DrawMCPrelim();
5733	+
5734	+	can->cd(3);
5735	+
5736	+
5737	+	write_warning(__FUNCTION__,"Missing negative JZB subtraction  ");
5738	+
5739	+	DrawMCPrelim();
5740	+
5741	+	/*
5742	+	THStack *Subtracted = new THStack();
5743	+
5744	+	TIter nextSF(JZBdistSF.GetHists());
5745	+	TIter nextOF(JZBdistOF.GetHists());
5746	+	TObject* SFobj;
5747	+	TObject* OFobj;
5748	+
5749	+	TH1* sfh = NULL;
5750	+	TH1* ofh = NULL;
5751	+
5752	+	while ( (SFobj = nextSF()) && (OFobj = nextOF()) ) {
5753	+	TH1F *sfstart = (TH1*)SFobj->Clone();
5754	+	TH1F *ofstart = (TH1*)OFobj->Clone();
5755	+	TH1F *obs, *pred;
5756	+	MakeObservationPrediction(obs,pred,sfstart,ofstart);
5757	+
5758	+	if ( !hratio ) {
5759	+	hratio = (TH1*)obj->Clone();
5760	+	hratio->SetName("hratio");
5761	+	} else hratio->Add( (TH1*)obj );
5762	+
5763	+	}
5764	+	hratio->Divide(hdata);
5765	+	*/
5766	+
5767	+
5768	+	can->cd(4);
5769	+	JZBdistOF.Draw("histo");
5770	+	DrawMCPrelim();
5771	+
5772	+
5773	+
5774	+	CompleteSave(can,"Systematics/TtbarTests/ControlRegion");
5775	+
5776	+	CleanLegends();
5777	+	delete can;
5778	+	}
5779	+
5780	+	void ttbar_tests(string mcjzb) {
5781	+	//ttbar_correction_tests();
5782	+	ttbar_region_search(mcjzb);
5783	+	}
5784
5785		void test() {
5786
5787	+
5788	+	float MCPeak, MCPeakError, DataPeak, DataPeakError;
5789	+	stringstream result, datajzb, mcjzb;
5790	+	bool doPUreweighting=true;
5791	+	bool SwitchOffNJetsCut=false;
5792	+
5793	+
5794	+	TCanvas *can = new TCanvas("can","can");
5795	+	float mcSigma,mcSigmaError,dataSigma,dataSigmaError;
5796	+
5797	+	float NumVtxBin[6] = {0,5,10,15,20,30};
5798	+
5799	+	stringstream NowCut;
5800	+
5801	+	TGraphErrors *gMCPeak = new TGraphErrors();
5802	+	gMCPeak->SetTitle("gMCPeak");
5803	+	gMCPeak->SetName("gMCPeak");
5804	+	TGraphErrors *gDataPeak = new TGraphErrors();
5805	+	gDataPeak->SetTitle("gDataPeak");
5806	+	gDataPeak->SetName("gDataPeak");
5807	+	TGraphErrors *gMCWidth = new TGraphErrors();
5808	+	gMCWidth->SetTitle("gMCWidth");
5809	+	gMCWidth->SetName("gMCWidth");
5810	+	TGraphErrors *gDataWidth = new TGraphErrors();
5811	+	gDataWidth->SetTitle("gDataWidth");
5812	+	gDataWidth->SetName("gDataWidth");
5813	+
5814	+	float AllMCPeaks[200];
5815	+	float AllMCPeaksErrors[200];
5816	+
5817	+	for(int i=10;i<180;i+=20) {
5818	+	NowCut.str("");
5819	+	NowCut << "mll>=" << i << "&&mll<" << i+20;
5820	+	find_one_peak_combination(TCut(NowCut.str().c_str()),SwitchOffNJetsCut,MCPeak,MCPeakError, DataPeak,DataPeakError,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
5821	+	AllMCPeaks[i]=MCPeak;
5822	+	AllMCPeaksErrors[i]=MCPeakError;
5823	+	}
5824	+
5825	+	dout << "Summary : " << endl;
5826	+	for(int i=10;i<180;i+=20) {
5827	+	dout << "For the slice at [" << i << " , " << i+20 << "] we get a peak at " << AllMCPeaks[i] << " +/- " << AllMCPeaksErrors[i] << endl;
5828	+	}
5829	+	/*
5830		TCanvas *testcanv = new TCanvas("testcanv","testcanv");
5831		testcanv->cd();
5832		//  switch_overunderflow(true);
5833	<	TH1F *ptdistr   = allsamples.Draw("ptdistr","pt1",100,30,200, "p_{T} [GeV]", "events", cutOSSF,data,luminosity);
5834	<	switch_overunderflow(false);
5835	<	ptdistr->Draw();
5836	<	testcanv->SaveAs("test.png");
5837	<	dout << "HELLO there!" << endl;
5833	>	TH1F *histo1 = new TH1F("histo","histo",100,0,1);
5834	>	histo1->Sumw2();
5835	>	TH1F *histo2 = new TH1F("histo2","histo2",100,0,1);
5836	>	histo2->Sumw2();
5837	>
5838	>	for(int i=0;i<100;i++) {
5839	>	histo1->Fill(1.0/i,i);
5840	>	histo2->Fill(1.0/(i*i),i);
5841	>	}
5842	>
5843	>	histo1->Draw();
5844	>	histo2->Draw("e1,same");
5845
5846	+	Save_With_Ratio( histo1, histo2, testcanv->cd(), "Bullshit", false, false, "AnyCrap" );
5847	+
5848	+
5849	+	dout << "HELLO there!" << endl;
5850	+	*/
5851		}

Diff Legend

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