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.35 by buchmann, Mon Jul 9 16:26:00 2012 UTC vs.
Revision 1.47 by buchmann, Tue Aug 14 09:32:53 2012 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 &DataSigma, stringstream &result, bool doPUreweighting = true, string saveas="")
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="")
37		{
38		// Temporarily switch off PU reweighting, if asked
39		TCut weightbackup=cutWeight;
#	Line 52 | Line 52 | void find_one_peak_combination(TCut spec
52		if(method==doKM) {
53		//we only need this histo for the KM fitting...
54		rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&specialcut,mc,luminosity,allsamples.FindSample("TTJet"));
55	<	MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method,saveas);
56	<	DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method,saveas);
55	>	MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas);
56	>	DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas);
57		delete rawttbarjzbeemmMC;
58		}
59		else {
#	Line 62 | Line 62 | void find_one_peak_combination(TCut spec
62		reducedMC->Add(rawJZBemMC,-1);
63		reducedData->Add(rawJZBemData,-1);
64		//this is Kostas' way of doing it - we subtract em to get rid of some of the ttbar contribution (in reality, of flavor-symmetric contribution)
65	<	MCPeak=find_peak(reducedMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method,saveas);
66	<	DataPeak=find_peak(reducedData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method,saveas);
65	>	MCPeak=find_peak(reducedMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas);
66	>	DataPeak=find_peak(reducedData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas);
67		delete reducedMC;
68		delete reducedData;
69		}
#	Line 73 | Line 73 | void find_one_peak_combination(TCut spec
73
74		//  MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method);
75		//  DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method);
76	<	dout   << "   We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- ?? (not impl.)" << endl;
77	<	result << "   We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- ?? (not impl.)" << endl;
78	<	dout   << "   We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- ?? (not impl.)" << endl;
79	<	result << "   We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- ?? (not impl.)" << endl;
76	>	dout   << "   We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- " << DataSigmaError << endl;
77	>	result << "   We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- " << DataSigmaError << endl;
78	>	dout   << "   We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- " << MCSigmaError << endl;
79	>	result << "   We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- " << MCSigmaError << endl;
80		delete rawJZBeemmData;
81		delete rawJZBeemmMC;
82		delete rawJZBemData;
#	Line 98 | Line 98 | void find_peaks(float &MCPeak,float &MCP
98		float datammpeak,dataeepeak;
99		float datammpeakerr,dataeepeakerr;
100
101	<	float mcSigma, dataSigma;
101	>	float mcSigma,mcSigmaError, dataSigma, dataSigmaError;
102
103		dout << "Finding global peak : " << endl;
104	<	find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma, dataSigma,result,doPUreweighting,"");
104	>	find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
105
106		if(DoInvidualeemmPeaks) {
107		dout << "Finding peak for electrons : " << endl;
108	<	find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr, dataeepeak,dataeepeakerr,mcSigma, dataSigma,result,doPUreweighting,"_ele");
108	>	find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele");
109		dout << "Finding peak for muons : " << endl;
110	<	find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr, datammpeak,datammpeakerr,mcSigma, dataSigma,result,doPUreweighting,"_mu");
110	>	find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu");
111
112		datajzb << "(" << jzbvariabledata;
113		mcjzb << "(" << jzbvariablemc;
#	Line 147 | Line 147 | void find_peaks(float &MCPeak,float &MCP
147
148		mcjzb << ")";
149		}
150	<
150	>
151
152		}
153
154	<	void make_special_obs_pred_mll_plot(string datajzb, string mcjzb, float jzbthreshold) {
154	>	void make_special_obs_pred_mll_plot(string datajzb, string mcjzb, float jzbthreshold, float binWidth = 5.0) {
155		float min=70.0;
156		float max=115.0;
157		if(!PlottingSetup::RestrictToMassPeak) {
158		min=20;
159	<	max=200;
159	>	max=300;
160		}
161	<	int nbins=int((max-min)/5);
161	>	int nbins=int((max-min)/binWidth);
162
163		TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
164
#	Line 195 | Line 195 | void make_special_obs_pred_mll_plot(stri
195		LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm","mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets&&smallerzeroD,data, luminosity);
196		}
197
198	<	TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
199	<	if(PlottingSetup::RestrictToMassPeak) {
200	<	Bpred->Add(RcorrJZBem,1.0/3.);
201	<	Bpred->Add(LcorrJZBem,-1.0/3.);
202	<	Bpred->Add(RcorrJZBSBem,1.0/3.);
203	<	Bpred->Add(LcorrJZBSBem,-1.0/3.);
204	<	Bpred->Add(RcorrJZBSBeemm,1.0/3.);
205	<	Bpred->Add(LcorrJZBSBeemm,-1.0/3.);
206	<	} else {
207	<	Bpred->Add(RcorrJZBem,1.0);
208	<	Bpred->Add(LcorrJZBem,-1.0);
209	<	}
210	<
198	>	// Separate predictions
199	>	TH1F* SFN = (TH1F*)LcorrJZBeemm->Clone("SFN");
200	>	TH1F* OFP = (TH1F*)RcorrJZBem->Clone("OFP");
201	>	TH1F* OFN = (TH1F*)LcorrJZBem->Clone("OFN");
202	>	if(PlottingSetup::RestrictToMassPeak) {
203	>	OFP->Scale(1.0/3.0);
204	>	OFP->Add(RcorrJZBSBem,1.0/3.);
205	>	OFP->Add(RcorrJZBSBeemm,1.0/3.);
206	>	OFN->Scale(1.0/3.0);
207	>	OFN->Add(LcorrJZBSBem,1.0/3.);
208	>	OFN->Add(LcorrJZBSBeemm,1.0/3.);
209	>	}
210	>
211	>	TH1F* Bpred = (TH1F*)SFN->Clone("Bpred");
212	>	Bpred->Add(OFP);
213	>	Bpred->Add(OFN,-1);
214		Bpred->SetLineColor(kRed);
215
216	+	RcorrJZBeemm->SetTitleOffset(1.3,"y");
217		RcorrJZBeemm->Draw();
218		mcRcorrJZBeemm.Draw("same");
219		Bpred->Draw("histo,same");
#	Line 225 | Line 229 | void make_special_obs_pred_mll_plot(stri
229		saveas << "kin/Mll_After_Cut/Cut_At" << jzbthreshold;
230		CompleteSave(ckin,saveas.str());
231
232	<
232	>	// Draw all predictions overlayed
233	>	unsigned int w = gStyle->GetHistLineWidth()+1; // Make line a bit wider, since we dash it
234	>	SFN->SetLineColor(kGreen+2);
235	>	SFN->SetLineStyle(2);
236	>	SFN->SetLineWidth(w);
237	>	OFP->SetLineColor(kBlue+2);
238	>	OFP->SetLineStyle(2);
239	>	OFP->SetLineWidth(w);
240	>	OFN->SetLineColor(kMagenta+2);
241	>	OFN->SetLineStyle(3);
242	>	OFN->SetLineWidth(w);
243
244		RcorrJZBeemm->Draw();
245	+	SFN->Draw("histo,same");
246	+	OFP->Draw("histo,same");
247	+	OFN->Draw("histo,same");
248		Bpred->Draw("histo,same");
249		RcorrJZBeemm->Draw("same");
250
251	<	TLegend *leg2 = make_legend();
252	<	leg2->SetX1(0.58);
236	<	leg2->AddEntry(RcorrJZBeemm,"observed (data)","l");
251	>	TLegend *leg2 = make_legend("",0.52,0.7);
252	>	leg2->AddEntry(RcorrJZBeemm,"observed (data)","lp");
253		leg2->AddEntry(Bpred,"predicted (data)","l");
254	+	leg2->AddEntry(SFN,  "  SF JZB<0","l");
255	+	leg2->AddEntry(OFN,  "  OF JZB<0","l");
256	+	leg2->AddEntry(OFP,  "  OF JZB>0","l");
257		leg2->Draw("same");
258
259		saveas.str("");
260		saveas << "kin/Mll_After_Cut/Cut_At" << jzbthreshold << "_nomc";
261		CompleteSave(ckin,saveas.str());
262
244	–
245	–
263		delete RcorrJZBeemm;
264		delete LcorrJZBeemm;
265		delete RcorrJZBem;
#	Line 341 | Line 358 | void draw_ratio_plot(TH1* hdata, THStack
358		oneline->Draw("same");
359		}
360
361	+	void make_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, bool logscale,
362	+	string xlabel, string filename, bool isPF=true, bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.3) {
363	+
364	+	TCut ibasiccut=basiccut;
365	+	bool draw_separation_lines=false;
366	+
367	+	if(isPF) ibasiccut=basiccut&&"pfjzb[0]>-998";
368	+
369	+	if(addcut != "") ibasiccut = ibasiccut && addcut.c_str();
370	+
371	+	TCut cutSF;
372	+	TCut cutOF;
373	+
374	+	cutOF=cutOSOF&&cutnJets&&ibasiccut;
375	+	cutSF=cutOSSF&&cutnJets&&ibasiccut;
376	+
377	+
378	+	TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
379	+	ckin->SetLogy(logscale);
380	+	TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cutSF,data,luminosity);
381	+	TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cutOF,data,luminosity);
382	+	string signal("LM3");
383	+	TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
384	+	int idx = signalsamples.FindSample(signal)[0];
385	+	(signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF);
386	+	signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
387	+	signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor);
388	+	datahistoSF->SetMarkerSize(DataMarkerSize);
389	+	datahistoOF->SetLineColor(kRed);
390	+
391	+	if ( !logscale ) {
392	+	datahistoSF->SetMinimum(0); // Defaults
393	+	datahistoOF->SetMinimum(0); // Defaults
394	+	} else {
395	+	datahistoSF->SetMinimum(0.5);
396	+	datahistoOF->SetMinimum(0.5);
397	+	}
398	+	if (logscale) {
399	+	datahistoSF->SetMaximum(5.3*datahistoSF->GetMaximum());
400	+	datahistoOF->SetMaximum(5.3*datahistoSF->GetMaximum());
401	+	} else {
402	+	datahistoSF->SetMaximum(2.0*datahistoSF->GetMaximum());
403	+	datahistoOF->SetMaximum(2.0*datahistoOF->GetMaximum());
404	+	}
405	+
406	+	datahistoSF->GetXaxis()->SetTitle(xlabel.c_str());
407	+	datahistoOF->GetXaxis()->SetTitle(xlabel.c_str());
408	+	signalhisto->GetXaxis()->SetTitle(xlabel.c_str());
409	+	datahistoSF->GetYaxis()->SetTitle("N. Events");
410	+	datahistoOF->GetYaxis()->SetTitle("N. Events");
411	+	signalhisto->GetYaxis()->SetTitle("N. Events");
412	+	datahistoSF->GetXaxis()->CenterTitle();
413	+	datahistoOF->GetXaxis()->CenterTitle();
414	+	signalhisto->GetXaxis()->CenterTitle();
415	+	datahistoSF->GetYaxis()->CenterTitle();
416	+	datahistoOF->GetYaxis()->CenterTitle();
417	+	signalhisto->GetYaxis()->CenterTitle();
418	+
419	+	TLegend *mleg = new TLegend(0.2+legendPosition, 0.72, 0.55+legendPosition, 0.87, legendTitle.c_str());
420	+	mleg->SetFillColor(kWhite);
421	+	mleg->SetTextFont(42);
422	+	mleg->SetTextSize(0.048);
423	+	mleg->SetLineWidth(0);
424	+	mleg->SetBorderSize(0);
425	+	mleg->AddEntry(datahistoSF, "S. Flavor", "PL");
426	+	mleg->AddEntry(datahistoOF, "O. Flavor", "L");
427	+	mleg->AddEntry(signalhisto, "LM3", "L");
428	+
429	+	datahistoSF->Draw("E1");
430	+	datahistoOF->Draw("HIST,SAMES");
431	+	signalhisto->Draw("HIST,SAMES");
432	+	mleg->Draw();
433	+	DrawPrelim();
434	+	CompleteSave(ckin, "./SFOF/" + filename);
435	+
436	+	datahistoSF->Delete();
437	+	datahistoOF->Delete();
438	+	signalhisto->Delete();
439	+	delete mleg;
440	+	delete ckin;
441	+
442	+	}
443	+
444	+
445	+
446		float lastrange_min=0;
447		float lastrange_max=0;
448
#	Line 398 | Line 500 | void make_kin_plot(string variable, stri
500		ckin->SetLogy(logscale);
501		TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity);
502		datahisto->SetMarkerSize(DataMarkerSize);
503	<	THStack mcstack = allsamples.DrawStack("mcstack",variable,nbins,min,max, xlabel, "events",cut,mc,luminosity);
503	>	if ( !logscale ) datahisto->SetMinimum(0); // Defaults
504	>	else datahisto->SetMinimum(0.5);
505	>	// Custom max.
506	>	if(variable=="TMath::Abs(pfJetDphiMet[0])") datahisto->SetMaximum(1.5*datahisto->GetMaximum());
507		if(variable=="pfJetGoodPt[0]") datahisto->SetMaximum(10*datahisto->GetMaximum());
508		if(variable=="pt") datahisto->SetMaximum(10*datahisto->GetMaximum());
509		if(filename=="mll_inclusive"||filename=="mll_inclusive_mm"||filename=="mll_inclusive_ee") datahisto->SetMinimum(1);
510		if(filename=="mll_osof") datahisto->SetMaximum(10*datahisto->GetMaximum());
511		if(filename=="mll_osof") datahisto->SetMinimum(9);
512	<	datahisto->SetMaximum(5.3*datahisto->GetMaximum());
513	<	datahisto->Draw("e1");
514	<	ckin->Update();
512	>	if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum());
513	>	else datahisto->SetMaximum(1.3*datahisto->GetMaximum());
514	>
515		if(loadlastminmax) {
516		datahisto->SetMinimum(lastrange_min);
517		datahisto->SetMaximum(lastrange_max);
#	Line 416 | Line 521 | void make_kin_plot(string variable, stri
521		}
522		}
523
524	+	// Draw signal by hand (for some reason I don't manage to use the sample class: it adds it to the stack!)
525	+	string signal("LM3");
526	+	TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max);
527	+	int idx = signalsamples.FindSample(signal)[0];
528	+	(signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut);
529	+	signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity);
530	+	signalhisto->SetLineColor(kOrange);
531	+
532	+	THStack mcstack   = allsamples.DrawStack("mcstack",  variable,nbins,min,max,xlabel,"events",cut,mc,luminosity);
533	+	datahisto->Draw("e1");
534	+	ckin->Update();
535		mcstack.Draw("same");
536	+
537		datahisto->Draw("same,e1");
538		TLegend *kinleg = allsamples.allbglegend();
539		kinleg->Draw();
#	Line 459 | Line 576 | void make_kin_plot(string variable, stri
576		kinpad->SetLogy(logscale);
577		datahisto->Draw("e1");
578		mcstack.Draw("same");
579	+	signalhisto->Draw("same");
580		datahisto->Draw("same,e1");
581		datahisto->Draw("same,axis");
582		if(RestrictToMassPeak&&draw_separation_lines) {
#	Line 466 | Line 584 | void make_kin_plot(string variable, stri
584		upperboundary->Draw("same");
585		}
586
587	+	kinleg->AddEntry("signalhisto",signal.c_str(),"l");
588		kinleg->Draw();
589		write_cut->Draw();
590		DrawPrelim();
#	Line 483 | Line 602 | void make_kin_plot(string variable, stri
602		delete ckin;
603		}
604
605	<	void make_JES_plot() {
605	>	void make_JES_plot(TCut cut, string name) {
606
607		int nbins=10;
608		float min=-0.5;
609		float max = 9.5;
610		bool logscale=true;
611		string xlabel="nJets";
493	–	TCut cut=cutmass&&cutOSSF&&basiccut;
612
613		TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600);
614		ckin->SetLogy(logscale);
615	<	TH1F *datahisto = allsamples.Draw("datahisto","pfJetGoodNum",nbins,min,max, xlabel, "events",cut,data,luminosity);
615	>	TH1F *datahisto = allsamples.Draw("datahisto","pfJetGoodNum40",nbins,min,max, xlabel, "events",cut,data,luminosity);
616		datahisto->SetMarkerSize(DataMarkerSize);
617	<	THStack mcstack = allsamples.DrawStack("mcstack","pfJetGoodNum",nbins,min,max, xlabel, "events",cut,mc,luminosity);
618	<	TH1F *JESup = allsamples.Draw("JESup","pfJetGoodNump1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity);
619	<	TH1F *JESdn = allsamples.Draw("JESdn","pfJetGoodNumn1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity);
617	>	THStack mcstack = allsamples.DrawStack("mcstack","pfJetGoodNum40",nbins,min,max, xlabel, "events",cut,mc,luminosity);
618	>	TH1F *JESup = allsamples.Draw("JESup","pfJetGoodNum40p1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity);
619	>	TH1F *JESdn = allsamples.Draw("JESdn","pfJetGoodNum40n1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity);
620
621		datahisto->SetMinimum(1);
622		datahisto->SetMaximum(5.3*datahisto->GetMaximum()); // in line with kinematic plots style
#	Line 523 | Line 641 | void make_JES_plot() {
641		TLegend *kinleg = allsamples.allbglegend();
642		kinleg->AddEntry(JESunc,"JES uncertainty","f");
643		kinleg->Draw();
644	<	CompleteSave(ckin,"Systematics/JES");
644	>	CompleteSave(ckin,"Systematics/JES"+name);
645		datahisto->Delete();
646		delete ckin;
647
#	Line 538 | Line 656 | void do_kinematic_plots(string mcjzb, st
656		float mll_hi=160;
657		if(!PlottingSetup::RestrictToMassPeak) {
658		mll_low=20;
659	<	mll_hi=210;
659	>	mll_hi=300;
660		}
661	<
662	<	make_kin_plot("met[4]","",40,0,200,dolog,"MET [GeV]","met",doPF,true);
661	>
662	>	make_OFSF_plot("pfJetGoodNum40", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 5, 3, 8, true, "N. jets", "njetsLow", false, false, false, 0.2);
663	>	make_OFSF_plot("pfJetGoodNum40", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 5, 3, 8, true, "N. jets", "njetsMed", false, false, false, 0.2);
664	>	make_OFSF_plot("pfJetGoodPt[0]/pfJetGoodPt[1]", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2Low", false, false, false, 0.2);
665	>	make_OFSF_plot("pfJetGoodPt[0]/pfJetGoodPt[1]", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2Med", false, false, false, 0.2);
666	>	make_OFSF_plot("TMath::Abs(pfJetDphiMet[0])", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphijetmetLow", false, false, false, 0.0);
667	>	make_OFSF_plot("TMath::Abs(pfJetDphiMet[0])", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphijetmetMed", false, false, false, 0.0);
668	>	make_OFSF_plot("TMath::Abs(dphi)", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphiLow", false, false, false, 0.2);
669	>	make_OFSF_plot("TMath::Abs(dphi)", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphiMed", false, false, false, 0.2);
670	>	make_OFSF_plot("TMath::Abs(dphiMet1)", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1Low", false, false, false, 0.2);
671	>	make_OFSF_plot("TMath::Abs(dphiMet1)", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1Med", false, false, false, 0.2);
672	>	make_OFSF_plot("TMath::Abs(dphiMet2)", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2Low", false, false, false, 0.2);
673	>	make_OFSF_plot("TMath::Abs(dphiMet2)", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2Med", false, false, false, 0.2);
674	>	make_OFSF_plot("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilcLow", false, false, false, 0.2);
675	>	make_OFSF_plot("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "mll>40&&mll<70&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilcMed", false, false, false, 0.2);
676	>	make_OFSF_plot("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijcLow", false, false, false, 0.2);
677	>	make_OFSF_plot("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "mll>20&&mll<40&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijcMed", false, false, false, 0.2);
678	>	make_OFSF_plot("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "Min(#Pi - |#Delta#phi(l,MET)|)", "dphilcoLow", false, false, false, 0.2);
679	>	make_OFSF_plot("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "mll>20&&mll<40&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "Min(#Pi - |#Delta#phi(l,MET)|)", "dphilcoMed", false, false, false, 0.2);
680	>	make_OFSF_plot("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "mll>20&&mll<40&&met[4]>100", "20 GeV < mll < 40 GeV", 16, 0, 3.2, false, "Min(#Pi - |#Delta#phi(jet,MET)|)", "dphijcoLow", false, false, false, 0.2);
681	>	make_OFSF_plot("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "mll>20&&mll<40&&met[4]>100", "40 GeV < mll < 70 GeV", 16, 0, 3.2, false, "Min(#Pi - |#Delta#phi(jet,MET)|)", "dphijcoMed", false, false, false, 0.2);
682	>
683	>
684	>
685	>	make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true);
686	>
687		make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true);
688		make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true);
689		make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true);
690		make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true);
691	<	make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
692	<	make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true);
693	<	make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true);
694	<	make_kin_plot("mll","",(int)((mll_hi-mll_low))/2,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
695	<	make_kin_plot("mll","",(int)((350-mll_low)),mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
696	<	make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF);
697	<	//  make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF);
698	<	//  make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF);
699	<	make_kin_plot("pt","",50,0,400,dolog,"Z p_{T} [GeV]","Zpt",doPF);
700	<	make_kin_plot("pt1","",50,0,100,nolog,"p_{T} [GeV]","pt1",doPF);
559	<	make_kin_plot("pt2","",50,0,100,nolog,"p_{T} [GeV]","pt2",doPF);
560	<	make_kin_plot("eta1","",40,-5,5,nolog,"#eta_{l}","eta",doPF);
561	<	make_kin_plot("jzb[1]","",100,-150,150,dolog,"JZB [GeV]","jzb_ossf",doPF);
691	>	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true);
692	>	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true);
693	>	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true);
694	>	make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true);
695	>	make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF);
696	>	if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]");
697	>
698	>	make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true);
699	>	make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true);
700	>
701		make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets",doPF);
702		make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF);
703		make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF);
704	<	if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]");
704	>
705	>	make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF);
706	>	//  make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF);
707	>	//  make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF);
708	>	make_kin_plot("pt","",50,0,500,dolog,"Z p_{T} [GeV]","Zpt",doPF);
709	>	make_kin_plot("pt1","",50,0,200,nolog,"p_{T} [GeV]","pt1",doPF);
710	>	make_kin_plot("pt2","",50,0,200,nolog,"p_{T} [GeV]","pt2",doPF);
711	>	make_kin_plot("eta1","",40,-3,3,nolog,"#eta_{l}","eta",doPF);
712	>	make_kin_plot("jzb[1]","",100,-150,200,dolog,"JZB [GeV]","jzb_ossf",doPF);
713		stringstream jzbcut;
714		jzbcut << "((is_data&&("<<datajzb<<")>100)||(!is_data&&("<<mcjzb<<")>100))";
715		make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true);
#	Line 581 | Line 728 | void do_kinematic_plots(string mcjzb, st
728		make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB50",doPF,true,true);
729		make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB50",doPF,true);
730		make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB50",doPF,true);
731	<
732	<
731	>
732	>	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",doPF,true);
733	>
734		make_special_obs_pred_mll_plot(datajzb,mcjzb,0);
735		make_special_obs_pred_mll_plot(datajzb,mcjzb,50);
736	+	make_special_obs_pred_mll_plot(datajzb,mcjzb,80);
737		make_special_obs_pred_mll_plot(datajzb,mcjzb,100);
738	<	make_special_obs_pred_mll_plot(datajzb,mcjzb,150);
739	<	make_special_obs_pred_mll_plot(datajzb,mcjzb,200);
740	<	make_special_obs_pred_mll_plot(datajzb,mcjzb,250);
741	<
742	<	make_JES_plot();
738	>	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,150);
739	>	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,200);
740	>	//   make_special_obs_pred_mll_plot(datajzb,mcjzb,250);
741	>
742	>	make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf");
743	>	make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof");
744	>
745	>
746		}
747
748		void make_comp_plot( string var, string xlabel, string filename, float jzbcut, string mcjzb, string datajzb,
749		int nbins, float xmin, float xmax, bool log,
750		float ymin=0, float ymax=0, bool leftJustified=false ) {
751	<	flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- WATCH OUT: the argument in the function changed!
751	>	flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- WATCH OUT: the argument in the function changed!
752
753		TCut weightbackup=cutWeight;//backing up the correct weight (restoring below!)
754		if(weightbackup==TCut("1.0")||weightbackup==TCut("1")) write_warning(__FUNCTION__,"WATCH OUT THE WEIGHT HAS POSSIBLY NOT BEEN RESET!!!! PLEASE CHANGE LINE "+any2string(__LINE__));
755	<	if(var=="numVtx") cutWeight=TCut("1.0");
755	>	//if(var=="numVtx") cutWeight=TCut("1.0");
756		TCut jzbData[]= { TCut(TString(datajzb+">"+any2string(jzbcut))),TCut(TString(datajzb+"<-"+any2string(jzbcut))) };
757		TCut jzbMC[]  = { TCut(TString(mcjzb+">"+any2string(jzbcut))),TCut(TString(mcjzb+"<-"+any2string(jzbcut))) };
758
#	Line 624 | Line 776 | void make_comp_plot( string var, string
776		TCanvas *ccomp = new TCanvas("ccomp","Comparison plot",600,400);
777		ccomp->SetLogy(log);
778		TH1F *datahisto = allsamples.Draw("datahisto",   var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbData[ijzb],data,luminosity);
779	<	TH1F *lm4histo = signalsamples.Draw("lm4histo",   var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbData[ijzb],data,luminosity,signalsamples.FindSample("LM"));
779	>	TH1F *lm3histo  = signalsamples.Draw("lm3histo",   var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbMC[ijzb],data,luminosity,signalsamples.FindSample("LM3"));
780		THStack mcstack = allsamples.DrawStack("mcstack",var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbMC[ijzb],  mc,  luminosity);
781		datahisto->SetMarkerSize(DataMarkerSize);
782		if (ymax>ymin) datahisto->SetMaximum(ymax);
783	<	lm4histo->SetLineStyle(2);
783	>	lm3histo->SetLineStyle(2);
784		datahisto->Draw("e1");
785		mcstack.Draw("same");
786		datahisto->Draw("same,e1");
787	<	//      lm4histo->Draw("hist,same");
787	>	lm3histo->Draw("hist,same");
788		TLegend *kinleg = allsamples.allbglegend((sRegions[iregion]+(ijzb?"neg":"pos")).c_str());
789		if ( leftJustified ) {
790		Float_t w = kinleg->GetX2()-kinleg->GetX1();
791		kinleg->SetX1(0.2);
792		kinleg->SetX2(0.2+w);
793		}
794	<	//      kinleg->AddEntry(lm4histo,"LM4","l");
794	>	kinleg->AddEntry(lm3histo,"LM3","l");
795		kinleg->Draw();
796		TText* write_variable = write_text(0.99,0.01,var);
797		write_variable->SetTextAlign(31);
#	Line 648 | Line 800 | void make_comp_plot( string var, string
800		CompleteSave(ccomp,"compare/JZBcut_at_"+any2string(jzbcut)+"/"+filename+"/"+filename+sRegions[iregion]+(ijzb?"neg":"pos"));
801		delete datahisto;
802		delete ccomp;
803	<	delete lm4histo;
803	>	delete lm3histo;
804		}
805		cutWeight=weightbackup;
806		}
#	Line 666 | Line 818 | void region_comparison_plots(string mcjz
818		float jzbcut=jzb_cuts[ijzb]; // Comparison plots are done for this JZB cut
819		float mll_low=50;float mll_high=170;
820		if(!PlottingSetup::RestrictToMassPeak) {
821	<	mll_high=200;
822	<	mll_low=50;
821	>	mll_high=300;
822	>	mll_low=20;
823		}
824	<	make_comp_plot("mll","m_{ll} [GeV]","mll",jzbcut,mcjzb,datajzb,30,mll_low,mll_high,false,0,16.);
824	>	make_comp_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",jzbcut,mcjzb,datajzb,100,0,10,true);
825	>	make_comp_plot("TMath::Abs(pfJetDphiMet[0])","|#Delta#phi(jet1,MET)|","dphiJ1MET",jzbcut,mcjzb,datajzb,32,0,3.2,false,0,0,true);
826	>
827	>	make_comp_plot("mll","m_{ll} [GeV]","mll",jzbcut,mcjzb,datajzb,56,mll_low,mll_high,false,0,16.);
828		make_comp_plot("met[4]","pfMET [GeV]","pfmet",jzbcut,mcjzb,datajzb,18,0,360,false,0,16.);
829	<	make_comp_plot("pfJetGoodNum","#(jets)","njets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.);
829	>	make_comp_plot("pfJetGoodNum40","#(jets)","njets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.);
830	>	make_comp_plot("pfJetGoodNumBtag","#(b-jets)","nBjets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.);
831		make_comp_plot("pt","Z p_{T} [GeV]","Zpt",jzbcut,mcjzb,datajzb,26,0,525,false,0.,21.);
832	<	make_comp_plot("numVtx","#(prim. vertices)","nvtx",jzbcut,mcjzb,datajzb,20,0.,20.,false,0,16.);
832	>	make_comp_plot("numVtx","#(prim. vertices)","nvtx",jzbcut,mcjzb,datajzb,40,0.,40.,false,0,16.);
833		make_comp_plot("TMath::Abs(dphi)","#Delta#phi(leptons)","dphilep",jzbcut,mcjzb,datajzb,10,0.,3.1415,false,0,16.,true);
834		make_comp_plot("TMath::Abs(dphi_sumJetVSZ[1])","#Delta#phi(Z,jets)","dphiZjets",jzbcut,mcjzb,datajzb,10,0.,3.1415,false,0,16.,true);
835	+	make_comp_plot("eta1","#eta_1","eta1",jzbcut,mcjzb,datajzb,10,0.,2.5,false,0,16.);
836	+	make_comp_plot("eta2","#eta_2","eta2",jzbcut,mcjzb,datajzb,10,0.,2.5,false,0,16.);
837		}
838
839		switch_overunderflow(false); // switching overflow/underflow bins off
#	Line 1278 | Line 1436 | void do_prediction_plot(string jzb, TCan
1436		string ytitle("ratio");
1437		if ( use_data==1 ) ytitle = "data/pred";
1438		//save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,Bpredsaveas,true,use_data!=1,ytitle);
1439	<	save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,true,ytitle,BpredSys);//not extending the y range anymore up to 4
1439	>	save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,false,ytitle,BpredSys);//not extending the y range anymore up to 4
1440
1441
1442		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
#	Line 1722 | Line 1880 | void draw_pure_jzb_histo(TCut cut,string
1880		datahisto->Draw("same,e1");
1881
1882		TLegend *leg;
1883	<	if(is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("OFZP",datahisto);
1884	<	else if(!is_OF(cut) &&  is_ZP(cut)) leg = allsamples.allbglegend("SFZP",datahisto);
1885	<	else if( is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("OFSB",datahisto);
1886	<	else if(!is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("SFSB",datahisto);
1887	<	else {
1888	<	std::cerr << "Unable to decode cut: " << cut.GetTitle() << std::endl;
1889	<	exit(-1);
1890	<	}
1883	>	if (!PlottingSetup::RestrictToMassPeak) {
1884	>	if(is_OF(cut)) leg = allsamples.allbglegend("Opposite flavor",datahisto);
1885	>	else leg = allsamples.allbglegend("Same flavor",datahisto);
1886	>	} else {
1887	>	if(is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("OFZP",datahisto);
1888	>	else if(!is_OF(cut) &&  is_ZP(cut)) leg = allsamples.allbglegend("SFZP",datahisto);
1889	>	else if( is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("OFSB",datahisto);
1890	>	else if(!is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("SFSB",datahisto);
1891	>	else {
1892	>	std::cerr << "Unable to decode cut: " << cut.GetTitle() << std::endl;
1893	>	exit(-1);
1894	>	}
1895	>	}
1896		leg->Draw();
1897		string write_cut = decipher_cut(cut,"");
1898		TText *writeline1 = write_cut_on_canvas(write_cut.c_str());
#	Line 1889 | Line 2052 | void jzb_plots(string mcjzb, string data
2052		flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
2053		if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_SFSB",can,binning);
2054		if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB",can,binning);
1892	–
2055		draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm",can,binning);
2056		draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm_coarse",can,coarse_binning);
2057	+	draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm_coarsest",can,coarsest_binning);
2058
2059		}
2060
#	Line 2243 | Line 2406 | float find_one_correction_factor(string
2406		TCanvas *cancorr = new TCanvas("cancorr","Canvas for Response Correction");
2407		cancorr->SetLogz();
2408		cancorr->SetRightMargin(0.13);
2409	<	TCut zptforresponsepresentation("pt<600"&&cutmass&&cutOSSF&&"((sumJetPt[1]/pt)<5.0)"&&SpecialCut&&passtrig);
2409	>	TCut zptforresponsepresentation(Restrmasscut&&SpecialCut&&passtrig);
2410	>
2411		if(PlottingSetup::DoBTag) zptforresponsepresentation=zptforresponsepresentation&&PlottingSetup::bTagRequirement;
2412		TH2F *niceresponseplotd = new TH2F("niceresponseplotd","",100,0,600,100,0,5);
2413	+	niceresponseplotd->Sumw2();
2414	+	TH2F* emuResponse = (TH2F*)niceresponseplotd->Clone("emuResponse");
2415		vector<int> SampleIndices=allsamples.FindSample(FindKeyword);
2416		for(int iSample=0;iSample<(int)SampleIndices.size();iSample++) {
2417		if((allsamples.collection)[SampleIndices[iSample]].is_data && !dodata) continue;
2418		if((allsamples.collection)[SampleIndices[iSample]].is_data ==false && dodata) continue;
2419
2420		dout << "   Response correction : Using sample " << (allsamples.collection)[SampleIndices[iSample]].filename << " for " << FindKeyword << endl;
2421	<	(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+niceresponseplotd",zptforresponsepresentation*cutWeight);
2421	>	(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+niceresponseplotd",(zptforresponsepresentation&&cutOSSF)*cutWeight);
2422	>	(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+emuResponse",(zptforresponsepresentation&&cutOSOF)*cutWeight);
2423		}
2424	<
2424	>	niceresponseplotd->Add(emuResponse,-1);
2425	>
2426		niceresponseplotd->SetStats(0);
2427		niceresponseplotd->GetXaxis()->SetTitle("Z p_{T} [GeV]");
2428		niceresponseplotd->GetYaxis()->SetTitle("Response");
#	Line 2262 | Line 2430 | float find_one_correction_factor(string
2430		niceresponseplotd->GetYaxis()->CenterTitle();
2431		niceresponseplotd->Draw("COLZ");
2432		TProfile * profd = (TProfile*)niceresponseplotd->ProfileX();
2433	+	profd->Rebin(2);
2434		profd->SetMarkerSize(DataMarkerSize);
2435		profd->Fit("pol0","","same,e1",100,400);
2436		DrawPrelim();
#	Line 2281 | Line 2450 | float find_one_correction_factor(string
2450		CompleteSave(cancorr,"ResponseCorrection/Response_Correction_Illustration_New_"+SaveAs);
2451		delete cancorr;
2452		delete niceresponseplotd;
2453	+	delete profd;
2454		return correction;
2455		}
2456
#	Line 2440 | Line 2610 | void ttbar_sidebands_comparison(string m
2610
2611
2612		if(!doPURW) {
2613	+	dout << "Not doing PU reweighting for ttbar closure test" << endl;
2614		cutWeight="1.0";
2615		// Do it without PU re-weighting
2616		float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0;
#	Line 2449 | Line 2620 | void ttbar_sidebands_comparison(string m
2620
2621		stringstream mcjzbnoPU;
2622		find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,resultsNoPU,true,noPUdatajzb,noPUmcjzb);
2623	<	if(MCPeakNoPU>0) mcjzbnoPU<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeakNoPU)<<")";
2453	<	else mcjzbnoPU<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeakNoPU)<<")";
2454	<	dout << "The peak corrected JZB expression for MC without pileup is : " <<  mcjzbnoPU.str() << endl;
2455	<
2456	<	mcjzb = mcjzbnoPU.str();
2457	<
2623	>	mcjzb = noPUmcjzb.str();
2624		}
2625
2626	<
2461	<
2462	<
2463	<
2464	<
2626	>
2627		float simulatedlumi = luminosity; //in pb please - adjust to your likings
2628
2629	<	TH1F *TZem = systsamples.Draw("TZem",mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2629	>	TH1F *TZem  = systsamples.Draw("TZem",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2630		TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2631		TH1F *TSem;
2632		TH1F *nTSem;
2633	<	TH1F *TZeemm = systsamples.Draw("TZeemm",mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2633	>	TH1F *TZeemm  = systsamples.Draw("TZeemm",     mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2634		TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2635		TH1F *TSeemm;
2636		TH1F *nTSeemm;
2637
2638		if(PlottingSetup::RestrictToMassPeak) {
2639	<	TSem = systsamples.Draw("TSem",mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2640	<	nTSem = systsamples.Draw("nTSem","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2641	<	TSeemm = systsamples.Draw("TSeemm",mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2639	>	TSem    = systsamples.Draw("TSem",       mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2640	>	nTSem   = systsamples.Draw("nTSem",  "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2641	>	TSeemm  = systsamples.Draw("TSeemm",     mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2642		nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets"));
2643		}
2644
2645		TCanvas *tcan = new TCanvas("tcan","tcan");
2484	–
2646		tcan->SetLogy(1);
2647
2648		TZeemm->SetLineColor(kBlack);
#	Line 2510 | Line 2671 | void ttbar_sidebands_comparison(string m
2671		histos.push_back(TSem);
2672		histos.push_back(TSeemm);
2673		}
2674	<	draw_all_ttbar_histos(tcan,histos,"histo",0.04);
2674	>	draw_all_ttbar_histos(tcan,histos,"histo",1);
2675
2676		TLegend *leg = make_legend("MC t#bar{t}",0.6,0.65,false);
2677		leg->AddEntry(TZeemm,"SFZP","l");
#	Line 2538 | Line 2699 | void ttbar_sidebands_comparison(string m
2699		TSeemm->Divide(TZeemm);
2700		TSem->SetMarkerStyle(24);
2701		TSeemm->SetMarkerStyle(32);
2702	<	}
2702	>	}
2703
2704		tcan->SetLogy(0);
2705		TZem->GetYaxis()->SetRangeUser(0,2.5);
#	Line 2587 | Line 2748 | void ttbar_sidebands_comparison(string m
2748		DrawMCPrelim(simulatedlumi);
2749		CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison_ratio");
2750
2751	+	if (0) { // Turn this off: we don't need this
2752
2753	<	///-------------- second part: only look at the quantity we actually care about!
2754	<	TH1F *leftsfzp  = (TH1F*)nTZeemm->Clone("leftsfzp");
2755	<	TH1F *rightsfzp = (TH1F*)TZeemmcopy->Clone("rightsfzp");
2756	<	rightsfzp->Add(leftsfzp,-1);
2757	<	TH1F *leftofzp  = (TH1F*)nTZem->Clone("leftofzp");
2758	<	TH1F *rightofzp = (TH1F*)TZemcopy->Clone("rightofzp");
2759	<	rightofzp->Add(leftofzp,-1);
2760	<	TH1F *leftofsb;
2761	<	TH1F *rightofsb;
2762	<	TH1F *leftsfsb;
2763	<	TH1F *rightsfsb;
2764	<	if(PlottingSetup::RestrictToMassPeak) {
2765	<	leftofsb  = (TH1F*)nTSem->Clone("leftofsb");
2766	<	rightofsb = (TH1F*)TSemcopy->Clone("rightofsb");
2767	<	rightofsb->Add(leftofsb,-1);
2768	<	leftsfsb  = (TH1F*)nTSeemm->Clone("leftsfsb");
2769	<	rightsfsb = (TH1F*)TSeemmcopy->Clone("rightsfsb");
2770	<	rightsfsb->Add(leftsfsb,-1);
2771	<	}
2753	>	///-------------- second part: only look at the quantity we actually care about!
2754	>	TH1F *leftsfzp  = (TH1F*)nTZeemm->Clone("leftsfzp");
2755	>	TH1F *rightsfzp = (TH1F*)TZeemmcopy->Clone("rightsfzp");
2756	>	rightsfzp->Add(leftsfzp,-1);
2757	>	TH1F *leftofzp  = (TH1F*)nTZem->Clone("leftofzp");
2758	>	TH1F *rightofzp = (TH1F*)TZemcopy->Clone("rightofzp");
2759	>	rightofzp->Add(leftofzp,-1);
2760	>	TH1F *leftofsb;
2761	>	TH1F *rightofsb;
2762	>	TH1F *leftsfsb;
2763	>	TH1F *rightsfsb;
2764	>	if(PlottingSetup::RestrictToMassPeak) {
2765	>	leftofsb  = (TH1F*)nTSem->Clone("leftofsb");
2766	>	rightofsb = (TH1F*)TSemcopy->Clone("rightofsb");
2767	>	rightofsb->Add(leftofsb,-1);
2768	>	leftsfsb  = (TH1F*)nTSeemm->Clone("leftsfsb");
2769	>	rightsfsb = (TH1F*)TSeemmcopy->Clone("rightsfsb");
2770	>	rightsfsb->Add(leftsfsb,-1);
2771	>	}
2772
2773	<	tcan->SetLogy(1);
2774	<	rightsfzp->GetXaxis()->SetRangeUser(0,binning[binning.size()-1]);
2775	<	rightsfzp->GetYaxis()->SetTitle("#deltaJZB / 25 GeV");
2776	<	rightsfzp->Draw("histo");
2777	<	rightofzp->Draw("histo,same");
2778	<	if(PlottingSetup::RestrictToMassPeak) {
2779	<	rightofsb->Draw("histo,same");
2780	<	rightsfsb->Draw("histo,same");
2781	<	}
2782	<	DrawMCPrelim(simulatedlumi);
2773	>	tcan->SetLogy(1);
2774	>	rightsfzp->GetXaxis()->SetRangeUser(0,binning[binning.size()-1]);
2775	>	rightsfzp->GetYaxis()->SetTitle("#deltaJZB / 25 GeV");
2776	>	rightsfzp->Draw("histo");
2777	>	rightofzp->Draw("histo,same");
2778	>	if(PlottingSetup::RestrictToMassPeak) {
2779	>	rightofsb->Draw("histo,same");
2780	>	rightsfsb->Draw("histo,same");
2781	>	}
2782	>	DrawMCPrelim(simulatedlumi);
2783
2784	<	TLegend *legA = make_legend("MC t#bar{t}",0.6,0.65,false);
2785	<	legA->AddEntry(rightsfzp,"SFZP","l");
2786	<	legA->AddEntry(rightofzp,"OFZP","l");
2787	<	if(PlottingSetup::RestrictToMassPeak) {
2788	<	legA->AddEntry(rightofsb,"SFSB","l");
2789	<	legA->AddEntry(rightsfsb,"OFSB","l");
2790	<	}
2791	<	legA->Draw();
2784	>	TLegend *legA = make_legend("MC t#bar{t}",0.6,0.65,false);
2785	>	legA->AddEntry(rightsfzp,"SFZP","l");
2786	>	legA->AddEntry(rightofzp,"OFZP","l");
2787	>	if(PlottingSetup::RestrictToMassPeak) {
2788	>	legA->AddEntry(rightofsb,"SFSB","l");
2789	>	legA->AddEntry(rightsfsb,"OFSB","l");
2790	>	}
2791	>	legA->Draw();
2792
2793	<	CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison");
2793	>	CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison");
2794
2795	<	tcan->SetLogy(0);
2796	<	rightofzp->Divide(rightsfzp);
2797	<	rightofzp->GetXaxis()->SetRangeUser(0.0,binning[binning.size()-1]);
2798	<	rightofzp->GetYaxis()->SetRangeUser(0.0,2.5);
2799	<	rightofzp->GetYaxis()->SetTitle("#deltaJZB ratio");
2800	<	rightofzp->Draw();
2801	<	if(PlottingSetup::RestrictToMassPeak) {
2802	<	rightofsb->Divide(rightsfzp);
2803	<	rightsfsb->Divide(rightsfzp);
2804	<	rightofsb->Draw("same");
2805	<	rightsfsb->Draw("same");
2806	<	}
2795	>	tcan->SetLogy(0);
2796	>	rightofzp->Divide(rightsfzp);
2797	>	rightofzp->GetXaxis()->SetRangeUser(0.0,binning[binning.size()-1]);
2798	>	rightofzp->GetYaxis()->SetRangeUser(0.0,2.5);
2799	>	rightofzp->GetYaxis()->SetTitle("#deltaJZB ratio");
2800	>	rightofzp->Draw();
2801	>	if(PlottingSetup::RestrictToMassPeak) {
2802	>	rightofsb->Divide(rightsfzp);
2803	>	rightsfsb->Divide(rightsfzp);
2804	>	rightofsb->Draw("same");
2805	>	rightsfsb->Draw("same");
2806	>	}
2807
2808	<	TLine *top2 = new TLine(0.0,1.0+linepos,binning[binning.size()-1],1.0+linepos);
2809	<	TLine *center2 = new TLine(0.0,1.0,binning[binning.size()-1],1.0);
2810	<	TLine *bottom2 = new TLine(0.0,1.0-linepos,binning[binning.size()-1],1.0-linepos);
2811	<
2812	<	top2->SetLineColor(kBlue);top2->SetLineStyle(2);
2813	<	bottom2->SetLineColor(kBlue);bottom2->SetLineStyle(2);
2814	<	center2->SetLineColor(kBlue);
2815	<
2816	<	top2->Draw("same");
2817	<	center2->Draw("same");
2818	<	bottom2->Draw("same");
2808	>	TLine *top2 = new TLine(0.0,1.0+linepos,binning[binning.size()-1],1.0+linepos);
2809	>	TLine *center2 = new TLine(0.0,1.0,binning[binning.size()-1],1.0);
2810	>	TLine *bottom2 = new TLine(0.0,1.0-linepos,binning[binning.size()-1],1.0-linepos);
2811	>
2812	>	top2->SetLineColor(kBlue);top2->SetLineStyle(2);
2813	>	bottom2->SetLineColor(kBlue);bottom2->SetLineStyle(2);
2814	>	center2->SetLineColor(kBlue);
2815	>
2816	>	top2->Draw("same");
2817	>	center2->Draw("same");
2818	>	bottom2->Draw("same");
2819
2820	<	rightofzp->SetMarkerStyle(21);
2821	<	if(PlottingSetup::RestrictToMassPeak) {
2822	<	rightofsb->SetMarkerStyle(24);
2823	<	rightsfsb->SetMarkerStyle(32);
2824	<	}
2820	>	rightofzp->SetMarkerStyle(21);
2821	>	if(PlottingSetup::RestrictToMassPeak) {
2822	>	rightofsb->SetMarkerStyle(24);
2823	>	rightsfsb->SetMarkerStyle(32);
2824	>	}
2825
2826	<	TLegend *leg3 = make_legend("MC t#bar{t}",0.55,0.75,false);
2827	<	leg3->AddEntry(rightofzp,"OFZP / SFZP","ple");
2828	<	if(PlottingSetup::RestrictToMassPeak) {
2829	<	leg3->AddEntry(rightsfsb,"SFSB / SFZP","ple");
2830	<	leg3->AddEntry(rightofsb,"OFSB / SFZP","ple");
2831	<	}
2832	<	leg3->AddEntry(bottom,"syst. envelope","l");
2833	<	leg3->SetX1(0.25);leg3->SetX2(0.6);
2834	<	leg3->SetY1(0.65);
2826	>	TLegend *leg3 = make_legend("MC t#bar{t}",0.55,0.75,false);
2827	>	leg3->AddEntry(rightofzp,"OFZP / SFZP","ple");
2828	>	if(PlottingSetup::RestrictToMassPeak) {
2829	>	leg3->AddEntry(rightsfsb,"SFSB / SFZP","ple");
2830	>	leg3->AddEntry(rightofsb,"OFSB / SFZP","ple");
2831	>	}
2832	>	leg3->AddEntry(bottom,"syst. envelope","l");
2833	>	leg3->SetX1(0.25);leg3->SetX2(0.6);
2834	>	leg3->SetY1(0.65);
2835
2836	<	leg3->Draw("same");
2836	>	leg3->Draw("same");
2837
2838	<	DrawMCPrelim(simulatedlumi);
2839	<	CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison_ratio");
2838	>	DrawMCPrelim(simulatedlumi);
2839	>	CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison_ratio");
2840	>
2841	>	}
2842
2843		delete TZem;
2844		delete nTZem;
#	Line 2721 | Line 2885 | void ttbar_sidebands_comparison(string m
2885		nicer_binning.push_back(100);
2886		nicer_binning.push_back(125);
2887		nicer_binning.push_back(150);
2888	<	nicer_binning.push_back(175);
2888	>	//nicer_binning.push_back(175);
2889		nicer_binning.push_back(200);
2890	<	nicer_binning.push_back(230);
2891	<	nicer_binning.push_back(280);
2890	>	nicer_binning.push_back(250);
2891	>	nicer_binning.push_back(300);
2892		nicer_binning.push_back(400);
2893
2894		ttbar_sidebands_comparison(mcjzb,nicer_binning, "ttbar/");
#	Line 2733 | Line 2897 | void ttbar_sidebands_comparison(string m
2897
2898		void zjets_prediction_comparison(string mcjzbWithPUa) {
2899		TCanvas *zcan = new TCanvas("zcan","zcan");
2900	<	zcan->SetLogy(1);
2900	>	//  zcan->SetLogy(1);
2901		TCut weightbackup=cutWeight;
2902	<
2903	<
2904	<	// Do it without PU re-weighting
2905	<	float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0;
2906	<	stringstream resultsNoPU;
2907	<	stringstream noPUdatajzb;
2908	<	stringstream noPUmcjzb;
2909	<
2910	<	stringstream mcjzbnoPU;
2911	<	//  void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb)
2912	<	find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,resultsNoPU,true,noPUdatajzb,noPUmcjzb);
2913	<	if(MCPeakNoPU>0) mcjzbnoPU<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeakNoPU)<<")";
2914	<	else mcjzbnoPU<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeakNoPU)<<")";
2915	<	dout << "The peak corrected JZB expression for MC without pileup is : " <<  mcjzbnoPU.str() << endl;
2916	<
2917	<	string mcjzb = mcjzbnoPU.str();
2918	<
2919	<	cutWeight="1.0";
2902	>
2903	>	bool UsePURW=false;
2904	>
2905	>
2906	>	string mcjzb;
2907	>	if(UsePURW) {
2908	>	mcjzb=mcjzbWithPUa;
2909	>	} else {
2910	>	// Do it without PU re-weighting
2911	>	float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0;
2912	>	stringstream resultsNoPU;
2913	>	stringstream noPUdatajzb;
2914	>	stringstream noPUmcjzb;
2915	>
2916	>	find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,resultsNoPU,false,noPUdatajzb,noPUmcjzb);
2917	>	dout << "The peak corrected JZB expression for MC without pileup is : " <<  noPUmcjzb.str() << endl;
2918	>
2919	>	mcjzb = noPUmcjzb.str();
2920	>
2921	>	cutWeight="1.0";
2922	>	}
2923
2924
2925		vector<float> binning;
#	Line 2760 | Line 2927 | void zjets_prediction_comparison(string
2927		binning.push_back(10);
2928		binning.push_back(20);
2929		binning.push_back(40);
2930	<	//  binning.push_back(80);
2930	>	binning.push_back(60);
2931	>	//   binning.push_back(50);
2932	>	//   binning.push_back(60);
2933	>	//   binning.push_back(70);
2934	>	//   binning.push_back(80);
2935	>	//   binning.push_back(90);
2936		binning.push_back(100);
2937	<	//  float sbg_min=0.;
2766	<	//  float sbg_max=100.;
2767	<	//  int sbg_nbins=5;
2937	>
2938		float simulatedlumi = luminosity;//in pb please - adjust to your likings
2939
2940		TCut kPos((mcjzb+">0").c_str());
2941		TCut kNeg((mcjzb+"<0").c_str());
2942		string var( "abs("+mcjzb+")" );
2943	+
2944	+	TCut notTau("abs(genMID1)!=15");
2945	+	TCut ONLYTau("mll>20");
2946	+
2947
2948	<	TCut kcut(cutmass&&cutOSSF&&cutnJets);
2949	<	TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",kcut&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
2950	<	TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",kcut&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
2948	>	TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
2949	>	TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&notTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
2950	>
2951		hJZBpos->SetLineColor(kBlack);
2952		hJZBneg->SetLineColor(kRed);
2953
#	Line 2782 | Line 2956 | void zjets_prediction_comparison(string
2956		hJZBneg->Draw("same,hist");
2957		hJZBpos->Draw("same,e1"); // So it's on top...
2958
2959	<	TLegend *leg = make_legend("MC Z+jets",0.55,0.75,false);
2959	>	TLegend *leg = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.55,0.75,false);
2960		leg->AddEntry(hJZBpos,"Observed","pe");
2961		leg->AddEntry(hJZBneg,"Predicted","l");
2962		leg->Draw("same");
2963		DrawMCPrelim(simulatedlumi);
2964	<	CompleteSave(zcan,"Systematics/zjets_prediction");
2964	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction");
2965
2966		TH1F* hratio = (TH1F*)hJZBpos->Clone("hratio");
2967		hratio->Divide(hJZBneg);
2968
2969		for(int i=1;i<=hJZBpos->GetNbinsX();i++) {
2970	<	cout << "Positive: " << hJZBpos->GetBinContent(i) << "   vs    Negative : " << hJZBneg->GetBinContent(i) << endl;
2970	>	cout << "Positive: " << hJZBpos->GetBinContent(i) << "   vs    Negative : " << hJZBneg->GetBinContent(i) << "   (ratio : " << hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i) << endl;
2971		}
2972
2973		zcan->SetLogy(0);
2974	<	hratio->GetYaxis()->SetRangeUser(0,2.5);
2974	>	hratio->GetYaxis()->SetRangeUser(0,2.0);
2975		hratio->GetYaxis()->SetTitle("Observed/Predicted");
2976		hratio->Draw("e1");
2977
#	Line 2814 | Line 2988 | void zjets_prediction_comparison(string
2988		center->Draw("same");
2989		bottom->Draw("same");
2990
2991	<	TLegend *leg2 = make_legend("MC Z+jets",0.25,0.75,false);
2991	>	TLegend *leg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false);
2992		leg2->AddEntry(hratio,"obs / pred","pe");
2993		leg2->AddEntry(bottom,"syst. envelope","l");
2994		leg2->Draw("same");
2995		DrawMCPrelim(simulatedlumi);
2996	<	CompleteSave(zcan,"Systematics/zjets_prediction_ratio");
2996	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio");
2997	>
2998	>	TCut reducedNJets(cutnJets);
2999	>
3000	>	TH1F *TAUhJZBpos = systsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3001	>	TH1F *LcorrJZBeemm     = systsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3002	>	TH1F *RcorrJZBem       = systsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3003	>	TH1F *LcorrJZBem       = systsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3004	>
3005	>	TH1F *RcorrJZBSBem;
3006	>	TH1F *LcorrJZBSBem;
3007	>	TH1F *RcorrJZBSBeemm;
3008	>	TH1F *LcorrJZBSBeemm;
3009	>
3010	>	if(PlottingSetup::RestrictToMassPeak) {
3011	>	RcorrJZBSBem   = systsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3012	>	LcorrJZBSBem   = systsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3013	>	RcorrJZBSBeemm = systsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY"));
3014	>	LcorrJZBSBeemm = systsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ONLYTau&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY"));
3015	>	}
3016	>
3017	>	TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred");
3018	>	if(PlottingSetup::RestrictToMassPeak) {
3019	>	Bpred->Add(RcorrJZBem,1.0/3.);
3020	>	Bpred->Add(LcorrJZBem,-1.0/3.);
3021	>	Bpred->Add(RcorrJZBSBem,1.0/3.);
3022	>	Bpred->Add(LcorrJZBSBem,-1.0/3.);
3023	>	Bpred->Add(RcorrJZBSBeemm,1.0/3.);
3024	>	Bpred->Add(LcorrJZBSBeemm,-1.0/3.);
3025	>	} else {
3026	>	Bpred->Add(RcorrJZBem,1.0);
3027	>	Bpred->Add(LcorrJZBem,-1.0);
3028	>	}
3029	>
3030	>	Bpred->SetLineColor(kRed);
3031	>
3032	>	TAUhJZBpos->SetLineColor(kBlack);
3033	>	Bpred->SetLineColor(kRed);
3034	>
3035	>	TAUhJZBpos->SetMinimum(1.0);
3036	>	TAUhJZBpos->Draw("e1");
3037	>	Bpred->Draw("same,hist");
3038	>	TAUhJZBpos->Draw("same,e1");
3039	>
3040	>	TLegend *TAUleg = make_legend("MC Z+jets #rightarrow ee,#mu#mu,#tau#tau",0.55,0.75,false);
3041	>	TAUleg->AddEntry(TAUhJZBpos,"Observed","pe");
3042	>	TAUleg->AddEntry(Bpred,"Predicted","l");
3043	>	TAUleg->Draw("same");
3044	>	DrawMCPrelim(simulatedlumi);
3045	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction");
3046	>
3047	>	TH1F* TAUhratio = (TH1F*)TAUhJZBpos->Clone("TAUhratio");
3048	>	TAUhratio->Divide(Bpred);
3049	>
3050	>	for(int i=1;i<=TAUhJZBpos->GetNbinsX();i++) {
3051	>	cout << "ee/mm/tautau observed: " << TAUhJZBpos->GetBinContent(i) << "   vs    predicted : " << Bpred->GetBinContent(i) << "   (ratio : " << TAUhJZBpos->GetBinContent(i) / Bpred->GetBinContent(i) << endl;
3052	>	}
3053	>
3054	>	zcan->SetLogy(0);
3055	>	TAUhratio->GetYaxis()->SetRangeUser(0,2.5);
3056	>	TAUhratio->GetYaxis()->SetTitle("Observed/Predicted");
3057	>	TAUhratio->Draw("e1");
3058	>
3059	>	top->Draw("same");
3060	>	center->Draw("same");
3061	>	bottom->Draw("same");
3062	>
3063	>	TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow #tau#tau",0.25,0.75,false);
3064	>	TAUleg2->AddEntry(TAUhratio,"obs / pred","pe");
3065	>	TAUleg2->AddEntry(bottom,"syst. envelope","l");
3066	>	TAUleg2->Draw("same");
3067	>	DrawMCPrelim(simulatedlumi);
3068	>	CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio");
3069	>
3070	>	delete Bpred;
3071	>	delete TAUhJZBpos;
3072	>	delete LcorrJZBeemm;
3073	>	delete RcorrJZBem;
3074	>	delete LcorrJZBem;
3075	>
3076	>	if(PlottingSetup::RestrictToMassPeak) {
3077	>	delete RcorrJZBSBem;
3078	>	delete LcorrJZBSBem;
3079	>	delete RcorrJZBSBeemm;
3080	>	delete LcorrJZBSBeemm;
3081	>	}
3082	>
3083
3084		delete zcan;
3085		cutWeight=weightbackup;
#	Line 3332 | Line 3592 | void do_mlls_plot(string mcjzb) {
3592		}
3593		}
3594
3595	<	void met_vs_jzb_plots() {
3595	>	void met_vs_jzb_plots(string datajzb, string mcjzb) {
3596
3597		TCanvas *canmetjzb = new TCanvas("canmet","MET vs JZB canvas");
3598		canmetjzb->SetRightMargin(0.16);
#	Line 3346 | Line 3606 | void met_vs_jzb_plots() {
3606		vector<int> selsamples = allsamples.FindSample(findme[ifind]);
3607		TH2F *metvsjzb = new TH2F("metvsjzb","metvsjzb",200,0,100,400,-100,100);
3608		for(int isel=0;isel<(int)selsamples.size();isel++) {
3609	<	cout << "Producing MET:JZB plot ... working on sample: " << allsamples.collection[selsamples[isel]].filename << endl;
3610	<	allsamples.collection[selsamples[isel]].events->Draw("jzb[1]:met[4]>>+metvsjzb",cutmass&&cutOSSF);
3609	>	dout << "Producing MET:JZB plot ... working on sample: " << allsamples.collection[selsamples[isel]].filename << endl;
3610	>	allsamples.collection[selsamples[isel]].events->Draw("jzb[1]:met[4]>>+metvsjzb",cutmass&&cutOSSF&&cutnJets);
3611		}
3612		metvsjzb->Scale(allsamples.collection[selsamples[0]].weight);
3613		metvsjzb->SetStats(0);
#	Line 3361 | Line 3621 | void met_vs_jzb_plots() {
3621		title->Draw();
3622		CompleteSave(canmetjzb,(string)"METvsJZBplots/"+findme[ifind]);
3623		}
3624	+
3625	+
3626	+	dout << "About to produce MET plot for DY split up by JZB" << endl;
3627	+
3628	+	int nbins=14;
3629	+	float low=0;
3630	+	float high=140;
3631	+
3632	+	stringstream sLEFT;
3633	+	sLEFT << "((" << mcjzb << ")<0)";
3634	+	TCut LEFT(sLEFT.str().c_str());
3635	+	stringstream sRIGHT;
3636	+	sRIGHT << "((" << mcjzb << ")>0)";
3637	+	TCut RIGHT(sRIGHT.str().c_str());
3638	+
3639	+	TH1F *metleft   = allsamples.Draw("metleft","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
3640	+	TH1F *metleftO  = allsamples.Draw("metleftO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
3641	+	TH1F *metright  = allsamples.Draw("metright","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
3642	+	TH1F *metrightO = allsamples.Draw("metrightO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
3643	+
3644	+	metleft->Add(metleftO,-1);
3645	+	metright->Add(metrightO,-1);
3646	+
3647	+	metleft->SetLineColor(kRed);
3648	+	metright->SetLineColor(kBlack);
3649	+	TPad *metpad = new TPad("metpad","metpad",0,0,1,1);
3650	+	metpad->cd();
3651	+	metpad->SetLogy(1);
3652	+	metleft->Draw("histo");
3653	+	metright->Draw("same");
3654	+	TLegend *lg = make_legend();
3655	+	lg->SetX1(0.5);
3656	+	lg->AddEntry(metright,"JZB>0 (OSOF corrected)","P");
3657	+	lg->AddEntry(metleft,"JZB<0 (OSOF corrected)","L");
3658	+	lg->Draw();
3659	+	save_with_ratio(metright,metleft,metpad->cd(),"METvsJZBplots/ComparingLeftToRightinMETspectrum");
3660	+	TPad *metpad2 = new TPad("metpad2","metpad2",0,0,1,1);
3661	+
3662	+	TH1F *metlefta   = allsamples.Draw("metlefta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
3663	+	TH1F *metleftOa  = allsamples.Draw("metleftOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets"));
3664	+	TH1F *metrighta  = allsamples.Draw("metrighta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
3665	+	TH1F *metrightOa = allsamples.Draw("metrightOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets"));
3666	+
3667	+	TH1F *Bpred =  (TH1F*)metlefta->Clone("Bpred");
3668	+	Bpred->Add(metleftOa,-1);
3669	+	Bpred->Add(metrightOa);
3670	+	TH1F *obs = (TH1F*)metrighta->Clone("obs");
3671	+
3672	+	metlefta->Add(metleftOa,-1);
3673	+	metrighta->Add(metrightOa,-1);
3674	+
3675	+	metlefta->SetLineColor(kRed);
3676	+	metpad2->cd();
3677	+	metlefta->Draw("histo");
3678	+	metrighta->Draw("same");
3679	+	lg->Draw();
3680	+	save_with_ratio(metrighta,metlefta,metpad2->cd(),"METvsJZBplots/ComparingLeftToRightinMET_type1_spectrum");
3681	+
3682	+	TPad *metpad3 = new TPad("metpad3","metpad3",0,0,1,1);
3683	+	metpad3->cd();
3684	+	Bpred->SetLineColor(kRed);
3685	+	Bpred->Draw("histo");
3686	+	obs->SetLineColor(kBlack);
3687	+	obs->Draw("same");
3688	+	TLegend *lg2 = make_legend();
3689	+	lg2->SetX1(0.5);
3690	+	lg2->AddEntry(obs,"observed","P");
3691	+	lg2->AddEntry(Bpred,"predicted","L");
3692	+	lg2->Draw();
3693	+
3694	+	save_with_ratio(obs,Bpred,metpad3->cd(),"METvsJZBplots/ComparingPredObsinMET_type1_spectrum");
3695	+	delete lg;
3696	+	delete canmetjzb;
3697	+	delete metleft;
3698	+	delete metleftO;
3699	+	delete metright;
3700	+	delete metrightO;
3701		}
3702
3703

Diff Legend

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