[ViewVC] Diff of: cvsroot/UserCode/cbrown/Development/Plotting/Modules/ZbTools.C

Comparing UserCode/cbrown/Development/Plotting/Modules/ZbTools.C (file contents):
Revision 1.16 by buchmann, Wed Nov 28 17:51:14 2012 UTC vs.
Revision 1.22 by buchmann, Mon Jan 21 14:43:36 2013 UTC

#	Line 1 \| Line 1
1		#include <iostream>
2		#include <vector>
3		#include <sys/stat.h>
4	+	#include <math.h>
5
6		#include <TCut.h>
7		#include <TROOT.h>
#	Line 21 \| Line 22 \| using namespace std;
22		using namespace PlottingSetup;
23
24		namespace ZbData {
24	–	vector<float> data_over_mc;
25	–
25		TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120&&id1==id2");
26	<	TCut LeadingB("Zb3010_bTagProbCSVBP[0]>0.898");
27	<	TCut EtaB("abs(Zb3010_pfJetGoodEta[0])<1.3");
28	<	TCut PhiZcut("abs(Zb3010_pfJetDphiZ[0])>2.7");
26	>	TCut LeadingB("ZbCHS3010_bTagProbCSVBP[0]>0.898");
27	>	TCut EtaB("abs(ZbCHS3010_pfJetGoodEta[0])<1.3");
28	>	TCut PhiZcut("abs(ZbCHS3010_pfJetDphiZ[0])>2.7");
29
30		}
31
32	+	TCut STDWEIGHT;
33	+
34		class ZbPointResult {
35		public:
36		ZbPointResult(string var, float ptlow, float pthigh, float alphalow, float alphahigh, Value Data, Value MC);
#	Line 102 \| Line 103 \| void SpecialCutFlow(string identifier);
103		////////////////////////////////////////* AND THIS LINE ////////////////////////////////////////////////////
104
105
106	<	void print_yield(TCut cut) {
106	>	class CutYields{
107	>	public:
108	>	float data;
109	>	float Zb;
110	>	float Zc;
111	>	float Zl;
112	>	float tt;
113	>	float SiTo;
114	>	float Dibosons;
115	>	float WJets;
116	>	float Other;
117	>
118	>	CutYields();
119	>	};
120	>
121	>	CutYields::CutYields() {
122	>	this->data=0;
123	>	this->Zb=0;
124	>	this->Zc=0;
125	>	this->Zl=0;
126	>	this->tt=0;
127	>	this->SiTo=0;
128	>	this->Dibosons=0;
129	>	this->WJets=0;
130	>	this->Other=0;
131	>	}
132	>
133	>	std::ostream &operator<<(std::ostream &ostr, CutYields y)
134	>	{
135	>	return ostr << y.data << ";" << y.Zb << ";" << y.Zc << ";" << y.Zl << ";" << y.tt << ";" << y.SiTo << ";" << y.Dibosons << ";" << y.WJets << ";" << endl;
136	>	}
137	>
138	>
139	>	CutYields print_yield(TCut cut) {
140	>
141	>	CutYields Yield;
142		TH1F *data = allsamples.Draw("data", "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,data,luminosity);
143		dout << " Yields for " << (const char*) cut << endl;
144	<	dout << " data: " << data->Integral() << endl;
144	>	// dout << " data: " << data->Integral() << endl;
145	>
146	>	Yield.data=data->Integral();
147		THStack mcm = allsamples.DrawStack("mc", "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,mc,luminosity);
148		int nhists = mcm.GetHists()->GetEntries();
111	–	dout << "Going to loop over " << nhists << " histograms " << endl;
149		TFile *test = new TFile("test.root","RECREATE");
150		mcm.Write();
151		test->Close();
152		for (int i = 0; i < nhists; i++) {
153		TH1* hist = static_cast<TH1*>(mcm.GetHists()->At(i));
154	<	cout << " " << hist->GetName() << ": " << hist->Integral() << endl;
154	>	string name=hist->GetName();
155	>	if(Contains(name,"t_bar_t")) Yield.tt+=hist->Integral();
156	>	if(Contains(name,"W_Jets")) Yield.WJets+=hist->Integral();
157	>	if(Contains(name,"Single_top")) Yield.SiTo+=hist->Integral();
158	>	if(Contains(name,"Dibosons")) Yield.Dibosons+=hist->Integral();
159	>	if(Contains(name,"Z_l_")) Yield.Zl+=hist->Integral();
160	>	if(Contains(name,"Z_c_")) Yield.Zc+=hist->Integral();
161	>	if(Contains(name,"Z_b_")) Yield.Zb+=hist->Integral();
162		}
163	<
163	>
164	>	cout << Yield << endl;
165		cout << endl << endl;
166
167		delete data;
168	+	CleanLegends();
169	+	DeleteStack(mcm);
170	+
171	+	return Yield;
172		}
173
174		void draw_kin_variable(string variable, TCut cut, int nbins, float min, float max, string xlabel, string saveas, bool dologscale, string speciallabel="") {
126	–
175		TCanvas *ckin = new TCanvas("ckin","kin variable canvas");
176	<	ckin->SetLogy(dologscale);
176	>	if(dologscale) ckin->SetLogy(1);
177		TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity);
178		datahisto->SetMarkerSize(DataMarkerSize);
179		THStack mcstack = allsamples.DrawStack("mcstack",variable,nbins,min,max, xlabel, "events",cut,mc,luminosity);
180		if(dologscale) datahisto->SetMaximum(3*datahisto->GetMaximum());
181		else datahisto->SetMaximum(1.5*datahisto->GetMaximum());
182	<	TText *speciall;
183	<	if(speciallabel!="") {
184	<	speciall = new TText(0.2,0.8,speciallabel.c_str());
185	<	speciall->Draw();
182	>
183	>	float SumMC=0.0;
184	>	TIter nextH(mcstack.GetHists());
185	>	TH1F* h;
186	>	while ( h = (TH1F*)nextH() ) {
187	>	SumMC += h->Integral();
188		}
189	+	float scale = datahisto->Integral()/SumMC;
190	+
191	+	// Re-scale stack to data integral
192	+	nextH = TIter(mcstack.GetHists());
193	+
194	+	while ( h = (TH1F*)nextH() ) {
195	+	h->Scale(scale);
196	+	}
197	+
198	+	TText *speciall;
199		datahisto->Draw("e1");
200		ckin->Update();
201	<	mcstack.Draw("same");
201	>	mcstack.Draw("histo,same");
202		datahisto->Draw("same,e1");
203		TLegend *kinleg = allsamples.allbglegend();
204	+
205		kinleg->Draw();
206
207		TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1);
208		kinpad->cd();
209	<	kinpad->SetLogy(dologscale);
209	>	if(dologscale) kinpad->SetLogy(1);
210		datahisto->Draw("e1");
211	<	mcstack.Draw("same");
211	>	mcstack.Draw("histo,same");
212		datahisto->Draw("same,e1");
213		datahisto->Draw("same,axis");
214		kinleg->Draw();
215		DrawPrelim();
216		saveas="kin/"+saveas;
217	<	save_with_ratio(datahisto,mcstack,kinpad->cd(),saveas);
157	<
158	<
217	>	kinpad->cd();
218
219	+	if(speciallabel!="") {
220	+	speciall = write_title(speciallabel);
221	+	speciall->SetX(0.18);
222	+	speciall->SetTextAlign(11);
223	+	speciall->SetTextSize(0.03);
224	+	speciall->SetY(0.85);
225	+	speciall->Draw();
226	+	}
227
228	<	ZbData::data_over_mc.push_back(datahisto->Integral()/CollapseStack(mcstack)->Integral());
228	>	save_with_ratio(datahisto,mcstack,kinpad->cd(),saveas);
229
230		datahisto->Delete();
231	+	delete kinleg;
232		delete ckin;
233	+	CleanLegends();
234	+	DeleteStack(mcstack);
235		}
236	<
236	>
237		void print_all_b_yields() {
238	+	vector<CutYields> yields;
239		cout << "Basic selection with a b jet" << endl;
240	<	print_yield(ZplusBsel&&"Zb3010_pfJetGoodNumBtag>0");
240	>	yields.push_back(print_yield(ZplusBsel&&"ZbCHS3010_pfJetGoodNumBtag>0"));
241		cout << "Leading jet is a b-jet " << endl;
242	<	print_yield(ZplusBsel&&LeadingB);
242	>	yields.push_back(print_yield(ZplusBsel&&LeadingB));
243		cout << "\|#eta_{b}\|<1.3 " << endl;
244	<	print_yield(ZplusBsel&&LeadingB&&EtaB);
244	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB));
245		cout << "\|#delta#phi(b<Z)\|>2.7" << endl;
246	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
247	<	cout << "10<ptZ<1000 GeV" << endl;
248	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000");
246	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut));
247	>	cout << "30<ptZ<1000 GeV" << endl;
248	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"));
249		cout << "#alpha < 0.35" << endl;
250	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.35");
250	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"&&"(ZbCHS3010_alpha)<0.35"));
251		cout << "#alpha < 0.3" << endl;
252	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.3");
253	<	cout << "#alpha < 0.2" << endl;
254	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.2");
252	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"&&"(ZbCHS3010_alpha)<0.3"));
253	>	/* cout << "#alpha < 0.2" << endl;
254	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"&&"(ZbCHS3010_alpha)<0.2"));
255		cout << "#alpha < 0.15" << endl;
256	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.15");
256	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"&&"(ZbCHS3010_alpha)<0.15"));
257		cout << "#alpha < 0.1" << endl;
258	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.1");
259	<	// cout << "#alpha < 0.05" << endl;
260	<	// print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.05");
258	>	yields.push_back(print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000"&&"(ZbCHS3010_alpha)<0.1"));*/
259	>
260	>	for(int iy=0;iy<yields.size();iy++) {
261	>	cout << yields[iy] << endl;
262	>	}
263		}
264
265	<	void TEST_DrawEvilCutFlow() {
265	>	void DrawEvilCutFlow() {
266		stringstream MegaCut;
267
268	<	MegaCut << " 1(" << (const char) (ZplusBsel&&TCut("Zb3010_pfJetGoodNumBtag>0")) << ")";
268	>	MegaCut << " 1(" << (const char) (ZplusBsel&&TCut("ZbCHS3010_pfJetGoodNumBtag>0")) << ")";
269		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB) << ")";
270		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB) << ")";
271		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut) << ")";
272	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000")) << ")";
273	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.3")) << ")";
274	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.2")) << ")";
275	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.15")) << ")";
276	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.1")) << ")";
277	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.05")) << ")";
272	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000")) << ")";
273	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.3")) << ")";
274	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.2")) << ")";
275	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.15")) << ")";
276	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.1")) << ")";
277	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.05")) << ")";
278		MegaCut << " ";
279
280		TCanvas *can = new TCanvas("can","can");
#	Line 214 \| Line 287 \| void TEST_DrawEvilCutFlow() {
287
288
289		TH1F *data = allsamples.Draw("data", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,data,luminosity);
290	<	THStack mcm = allsamples.DrawStack("mc", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,mc,luminosity);
290	>	THStack mcm = allsamples.DrawStack("mc", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,mc,luminosity);
291
292		float runningsum=0;
293		for(int i=data->GetNbinsX();i>0;i--) {
#	Line 245 \| Line 318 \| void TEST_DrawEvilCutFlow() {
318		data->GetXaxis()->SetBinLabel(3,"Leading jet is b-jet");
319		data->GetXaxis()->SetBinLabel(4,"\|#eta_{b}\|<1.3 ");
320		data->GetXaxis()->SetBinLabel(5,"\|#delta#phi(b,Z)\|>2.7");
321	<	data->GetXaxis()->SetBinLabel(6,"10<p_{T}^{Z}<1000 GeV");
321	>	data->GetXaxis()->SetBinLabel(6,"30<p_{T}^{Z}<1000 GeV");
322		data->GetXaxis()->SetBinLabel(7,"#alpha < 0.3");
323		data->GetXaxis()->SetBinLabel(8,"#alpha < 0.2");
324		data->GetXaxis()->SetBinLabel(9,"#alpha < 0.15");
#	Line 254 \| Line 327 \| void TEST_DrawEvilCutFlow() {
327		data->GetXaxis()->LabelsOption("v");
328
329		data->Draw("e1");
330	<	mcm.Draw("same");
330	>	mcm.Draw("histo,same");
331		data->Draw("same,e1,axis");
332		data->Draw("same,e1");
333		// data->Draw("same,TEXT");
#	Line 264 \| Line 337 \| void TEST_DrawEvilCutFlow() {
337		DrawPrelim();
338
339		CompleteSave(can,"CutFlow");
340	+	CleanLegends();
341	+	DeleteStack(mcm);
342	+
343		}
344
345		void draw_Zb_kin_vars() {
346	<	draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Z p_{T}","Official/Zpt",1);
347	<	draw_kin_variable("Zb3010_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Sub-Leading Jet Pt","Official/Jet2Pt",1);
348	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Leading Jet Pt (B)","Official/LeadingJetPt",1);
349	<	draw_kin_variable("Zb3010_pfJetGoodPt[1]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,2,"p_{T}^{J2} / p_{T}^{Z}","Official/SubLeadingJetPt_Over_ZPt",1);
350	<	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<50", 40,0,2,"#alpha","Official/alpha_pt_30_to_50",1);
351	<	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>50&&pt<75", 40,0,2,"#alpha","Official/alpha_50_to_75",1);
352	<	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>75&&pt<125", 40,0,2,"#alpha","Official/alpha_pt_75_to_125",1);
353	<	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>125&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_125_to_1000",1);
354	<	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_30_to_1000",1);
355	<
356	<	draw_kin_variable("Zb3010_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",0);
357	<
282	<	/*
283	<	draw_kin_variable("pt1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,200,"p_{T}^{l1}","Lep/pt1",1);
284	<	draw_kin_variable("pt1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id1==0",20,0,200,"p_{T}^{e1}","Lep/pt1_e",1);
285	<	draw_kin_variable("pt1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id1==1",20,0,200,"p_{T}^{#mu1}","Lep/pt1_m",1);
286	<	draw_kin_variable("pt2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,200,"p_{T}^{l2}","Lep/pt2",1);
287	<	draw_kin_variable("pt2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id2==0",20,0,200,"p_{T}^{e2}","Lep/pt2_e",1);
288	<	draw_kin_variable("pt2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id2==1",20,0,200,"p_{T}^{#mu2}","Lep/pt2_m",1);
289	<	draw_kin_variable("eta1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,-3.1,3.1,"#eta^{l1}","Lep/eta1",1);
290	<	draw_kin_variable("eta1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id1==0",20,-3.1,3.1,"#eta^{e1}","Lep/eta1_e",1);
291	<	draw_kin_variable("eta1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id1==1",20,-3.1,3.1,"#eta^{m1}","Lep/eta1_m",1);
292	<	draw_kin_variable("eta2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,-3.1,3.1,"#eta^{l1}","Lep/eta2",1);
293	<	draw_kin_variable("eta2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id2==0",20,-3.1,3.1,"#eta^{e2}","Lep/eta2_e",1);
294	<	draw_kin_variable("eta2",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000&&id2==1",20,-3.1,3.1,"#eta^{#mu2}","Lep/eta2_m",1);
295	<	draw_kin_variable("numVtx",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",30,0,30,"N(Vtx)","Lep/nVtx",1);
296	<	*/
346	>
347	>	draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000"),25,0,200,"Z p_{T}","Official/Zpt",1);
348	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000"),25,0,200,"Sub-Leading Jet Pt","Official/Jet2Pt",1);
349	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000"),25,0,200,"Leading Jet Pt (B)","Official/LeadingJetPt",1);
350	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[1]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000"),20,0,2,"p_{T}^{J2} / p_{T}^{Z}","Official/SubLeadingJetPt_Over_ZPt",1);
351	>	draw_kin_variable("ZbCHS3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<50"), 40,0,2,"#alpha","Official/alpha_pt_30_to_50",1);
352	>	draw_kin_variable("ZbCHS3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>50&&pt<75"), 40,0,2,"#alpha","Official/alpha_50_to_75",1);
353	>	draw_kin_variable("ZbCHS3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>75&&pt<125"), 40,0,2,"#alpha","Official/alpha_pt_75_to_125",1);
354	>	draw_kin_variable("ZbCHS3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>125&&pt<1000"),40,0,2,"#alpha","Official/alpha_pt_125_to_1000",1);
355	>	draw_kin_variable("ZbCHS3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000"),40,0,2,"#alpha","Official/alpha_pt_30_to_1000",1);
356	>
357	>	draw_kin_variable("ZbCHS3010_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&TCut("pt>30&&pt<1000&&pfBJetDphiZ[0]>0"),30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",0);
358		}
359
360		void draw_mpf_vars() {
361	<	ZbData::data_over_mc.clear();
362	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.3",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p3",0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
363	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.2",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p2",0,"#alpha<0.2, 10 GeV < p_{T}^{Z} < 1000 GeV");
364	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.15",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p15",0,"#alpha<0.15, 10 GeV < p_{T}^{Z} < 1000 GeV");
365	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.1",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p1",0,"#alpha<0.1, 10 GeV < p_{T}^{Z} < 1000 GeV");
366	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.05",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p05",0,"#alpha<0.05, 10 GeV < p_{T}^{Z} < 1000 GeV");
367	<
368	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.3",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p3",0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
369	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.2",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p2",0,"#alpha<0.2, 10 GeV < p_{T}^{Z} < 1000 GeV");
370	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.15",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p15",0,"#alpha<0.15, 10 GeV < p_{T}^{Z} < 1000 GeV");
371	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.1",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p1",0,"#alpha<0.1, 10 GeV < p_{T}^{Z} < 1000 GeV");
311	<	draw_kin_variable("Zb3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb3010_alpha<0.05",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz__mpf_alpha_smaller_0p05",0,"#alpha<0.05, 10 GeV < p_{T}^{Z} < 1000 GeV");
361	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.3"), 20,0.0,2.0,"Z+b MPF","mpf_alpha_smaller_0p3",0,"#alpha<0.3, 30 GeV < p_{T}^{Z} < 1000 GeV");
362	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.2"), 20,0.0,2.0,"Z+b MPF","mpf_alpha_smaller_0p2",0,"#alpha<0.2, 30 GeV < p_{T}^{Z} < 1000 GeV");
363	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.15"),20,0.0,2.0,"Z+b MPF","mpf_alpha_smaller_0p15",0,"#alpha<0.15, 30 GeV < p_{T}^{Z} < 1000 GeV");
364	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.1"), 20,0.0,2.0,"Z+b MPF","mpf_alpha_smaller_0p1",0,"#alpha<0.1, 30 GeV < p_{T}^{Z} < 1000 GeV");
365	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.05"),20,0.0,2.0,"Z+b MPF","mpf_alpha_smaller_0p05",0,"#alpha<0.05, 30 GeV < p_{T}^{Z} < 1000 GeV");
366	>
367	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.3"), 20,0.0,2.0,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p3", 0,"#alpha<0.3, 30 GeV < p_{T}^{Z} < 1000 GeV");
368	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.2"), 20,0.0,2.0,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p2", 0,"#alpha<0.2, 30 GeV < p_{T}^{Z} < 1000 GeV");
369	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.15"),20,0.0,2.0,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p15", 0,"#alpha<0.15, 30 GeV < p_{T}^{Z} < 1000 GeV");
370	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.1"), 20,0.0,2.0,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p1", 0,"#alpha<0.1, 30 GeV < p_{T}^{Z} < 1000 GeV");
371	>	draw_kin_variable("ZbCHS3010_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<1000&&ZbCHS3010_alpha<0.05"),20,0.0,2.0,"p_{T} b-jet / p_{T} Z","ptb_over_ptz__mpf_alpha_smaller_0p05",0,"#alpha<0.05, 30 GeV < p_{T}^{Z} < 1000 GeV");
372		}
373
314	–	/double GetMedian(TH1F histo) {
315	–	int numBins = histo->GetXaxis()->GetNbins();
316	–	Double_t x[numBins];
317	–	Double_t y[numBins];
318	–	for (int i = 1; i <= numBins; i++) {
319	–	x[i] = histo->GetBinCenter(i);
320	–	y[i] = histo->GetBinContent(i);
321	–	}
322	–	return TMath::Median(numBins, &x, &y);
323	–	}*/
324	–
374		Value get_Zb_data_over_mc(string ContainerName, ZbResultContainer &res, string variable, string variable2, TCut cut, string saveas, float pt1, float pt2, float a1, float a2) {
375		// write_warning(__FUNCTION__,"Debugging this function, therefore always returning 3.1415+/-1.010101"); return Value(3.1415,1.010101);
376		TCanvas *cn = new TCanvas("cn","cn");
377		string varname="MPF";
378	+
379	+	bool DoFit=false;
380	+	bool UseFit=false;
381	+
382	+
383		if(Contains(variable,"pfJetGood")) varname="R_{abs}";
384	<	TH1F *hdata = allsamples.Draw("hdata",variable,1200,0,30, varname, "events", cut,data,luminosity);
385	<	TH1F *hmc = allsamples.Draw("hmc" , variable,1200,0,30, varname, "events", cut, mc, luminosity);
384	>	// TH1F *hdata = allsamples.Draw("hdata",variable,1200,0,30, varname, "events", cut,data,luminosity);
385	>	// TH1F *hmc = allsamples.Draw("hmc" , variable,1200,0,30, varname, "events", cut, mc, luminosity);
386	>	TH1F *hdata = allsamples.Draw("hdata",variable,1200,0.5,1.5, varname, "events", cut,data,luminosity); // making it more precise
387	>	TH1F *hmc = allsamples.Draw("hmc" , variable,1200,0.5,1.5, varname, "events", cut, mc, luminosity); // making it more precise
388		hmc->SetLineColor(kBlue);
389
390		float a=hdata->GetMean();
#	Line 339 \| Line 395 \| Value get_Zb_data_over_mc(string Contain
395		float error = TMath::Sqrt( (1/(bb)) (dada) + ((aa)/(bb))db*db);
396
397		TF1 *gaus_data = new TF1("gaus_data","gaus",0.7,1.3);
342	–	gaus_data->SetParameter(0,hdata->GetMaximum());
343	–	gaus_data->SetParameter(2,hdata->GetMean());
344	–	gaus_data->SetParameter(1,hdata->GetRMS());
345	–	hdata->Fit(gaus_data);
398		TF1 *gaus_mc = new TF1("gaus_mc","gaus",0.7,1.3);
399	<	gaus_mc->SetParameter(0,hmc->GetMaximum());
400	<	gaus_mc->SetParameter(1,hmc->GetMean());
401	<	gaus_mc->SetParameter(2,hmc->GetRMS());
402	<	hmc->Fit(gaus_mc);
403	<
399	>	if(DoFit) {
400	>	gaus_data->SetParameter(0,hdata->GetMaximum());
401	>	gaus_data->SetParameter(2,hdata->GetMean());
402	>	gaus_data->SetParameter(1,hdata->GetRMS());
403	>	hdata->Fit(gaus_data);
404	>
405	>	gaus_mc->SetParameter(0,hmc->GetMaximum());
406	>	gaus_mc->SetParameter(1,hmc->GetMean());
407	>	gaus_mc->SetParameter(2,hmc->GetRMS());
408	>	hmc->Fit(gaus_mc);
409	>	}
410
411		cn->cd();
412	<	hdata->GetYaxis()->SetTitle("events / 0.1");
412	>	hdata->Rebin(0.05*1200/(1.5-0.5));
413	>	hmc->Rebin(0.05*1200/(1.5-0.5));
414	>	hdata->GetYaxis()->SetTitle("events / 0.05");
415		hdata->SetMarkerColor(kRed);
416	<	/* hdata->Rebin(4);
417	<	hmc->Rebin(4);
418	<
419	<	gaus_mc->SetParameter(0,gaus_mc->GetParameter(0)*4);//rebinning
360	<	gaus_data->SetParameter(0,gaus_data->GetParameter(0)*4);//rebinning
361	<	*/
362	<	gaus_mc->SetLineColor(kBlue);
363	<	gaus_data->SetLineColor(kRed);
416	>	if(DoFit) {
417	>	gaus_mc->SetLineColor(kBlue);
418	>	gaus_data->SetLineColor(kRed);
419	>	}
420
421		hdata->GetXaxis()->SetRangeUser(0,2);
422		hdata->GetYaxis()->SetTitleOffset(1.3);
423		hdata->Draw("e1");
424	+
425	+	hmc->Scale(hdata->Integral()/hmc->Integral());
426	+	if(hmc->GetMaximum()>hdata->GetMaximum()) hdata->SetMaximum((hmc->GetMaximum()+sqrt((double)hmc->GetMaximum()))*1.2);
427	+
428		hmc->Draw("histo,same");
429		hdata->Draw("e1,same");
430	<	gaus_data->Draw("same");
431	<	gaus_mc->Draw("same");
430	>	if(DoFit) {
431	>	gaus_data->Draw("same");
432	>	gaus_mc->Draw("same");
433	>	}
434	>
435		stringstream summary;
436		summary << "#splitline{#bf{" << varname << "}}{#splitline{" << pt1 << " < p_{T}^{Z} < " << pt2 << ", " << a1 << " #leq #alpha < " << a2 << "}{#splitline{data: " << std::setprecision(4) << a << " #pm " << std::setprecision(4) << da << "}{";
437		summary << "#splitline{mc: " << std::setprecision(4) << b << " #pm " << std::setprecision(4) << db << "}{#splitline{ratio: " << factor;
438	<	summary << " #pm " << error << "}{#splitline{}{#splitline{data fit: " << std::setprecision(4) << gaus_data->GetParameter(1) << "+/-";
439	<	summary << std::setprecision(4) << gaus_data->GetParError(1) << "}{#splitline{";
440	<	summary << "mc fit:" << std::setprecision(4) << gaus_mc->GetParameter(1) << "+/-" << std::setprecision(4) << gaus_mc->GetParError(1) << "}{ratio: ";
441	<	summary << gaus_data->GetParameter(1)/gaus_mc->GetParameter(1) << "+/-" << gaus_data->GetParameter(1)/gaus_mc->GetParameter(1) * sqrt((gaus_data->GetParError(1) * gaus_data->GetParError(1))/(gaus_data->GetParameter(1)gaus_data->GetParameter(1)) + (gaus_mc->GetParError(1) gaus_mc->GetParError(1))/(gaus_mc->GetParameter(1)*gaus_mc->GetParameter(1)) ) << "}}}}}}}}";
438	>	summary << " #pm " << error;
439	>	if(DoFit) {
440	>	summary << "}{#splitline{}{#splitline{data fit: " << std::setprecision(4) << gaus_data->GetParameter(1) << "+/-";
441	>	summary << std::setprecision(4) << gaus_data->GetParError(1) << "}{#splitline{";
442	>	summary << "mc fit:" << std::setprecision(4) << gaus_mc->GetParameter(1) << "#pm" << std::setprecision(4) << gaus_mc->GetParError(1) << "}{ratio: ";
443	>	summary << gaus_data->GetParameter(1)/gaus_mc->GetParameter(1) << "#pm" << gaus_data->GetParameter(1)/gaus_mc->GetParameter(1) * sqrt((gaus_data->GetParError(1) * gaus_data->GetParError(1))/(gaus_data->GetParameter(1)gaus_data->GetParameter(1)) + (gaus_mc->GetParError(1) gaus_mc->GetParError(1))/(gaus_mc->GetParameter(1)*gaus_mc->GetParameter(1)) ) << "}}}}";
444	>	summary << "}}}}";
445	>	} else {
446	>	summary << "}{}}}}}";
447	>	}
448	>
449	>	if(UseFit) {
450	>	write_error(__FUNCTION__,"Requested using fit result for final result - not implemented yet!");
451	>	assert(0);
452	>	}
453	>
454
455		TText *infobox = write_title(summary.str());
456		infobox->SetX(0.75);
#	Line 388 \| Line 463 \| Value get_Zb_data_over_mc(string Contain
463		CompleteSave(cn,ContainerName+"/ResponsePlots/"+saveas);
464
465		res.MPF_Rabs_location.push_back(ZbPointResult(variable,pt1,pt2,a1,a2,Value(a,da),Value(b,db)));
466	<	cout << "Data;" << a << ";"<<da<<";MC"<<b<<";"<<db<<endl;
466	>	cout << "Data;" << a << ";"<<da<<";MC;"<<b<<";"<<db<<endl;
467		delete cn;
468		delete hdata;
469		delete hmc;
#	Line 396 \| Line 471 \| Value get_Zb_data_over_mc(string Contain
471		return Value(factor,error);
472		}
473
474	<	ZbResultContainer new_data_mc_agreement_2d(string AlphaVariable, string ContainerName, TCut BTagCut, TCut BTagWgt) {
474	>	int nFakes=1;
475	>
476	>	float PurityInclusive=0.05;
477	>	float PurityLoose=0.20;
478	>	float PurityMedium=0.54;
479	>	float PurityTight=0.82;
480	>
481	>
482	>	void DeterminePurity(TCut basecut, int nptcuts, float ptcuts[], string AlphaVariable, string ContainerName, float FaceValueToBeTakenAt) {
483	>
484	>	bool ComputePurity=false;
485	>
486	>	if(ContainerName=="inclusive") ComputePurity=true;
487	>	if(ContainerName=="loose") ComputePurity=true;
488	>	if(ContainerName=="medium") ComputePurity=true;
489	>	if(ContainerName=="reference") ComputePurity=true;
490	>
491	>	if(!ComputePurity) return;
492	>
493	>	stringstream specialcut;
494	>	specialcut << "((pt>30&& pt< 1000) && (" << AlphaVariable << "<" << FaceValueToBeTakenAt << "))";
495	>
496	>	THStack mcm = allsamples.DrawStack("mcm" , "mpf",1200,0,30, "MPF", "events", TCut(basecut && specialcut.str().c_str()), mc, luminosity);
497	>	int nhists = mcm.GetHists()->GetEntries();
498	>
499	>	float total=0;
500	>	float Zb=0;
501	>
502	>	for (int i = 0; i < nhists; i++) {
503	>	TH1* hist = static_cast<TH1*>(mcm.GetHists()->At(i));
504	>	string name=hist->GetName();
505	>	if(Contains(name,"t_bar_t")) total+=hist->Integral();
506	>	if(Contains(name,"W_Jets")) total+=hist->Integral();
507	>	if(Contains(name,"Single_top")) total+=hist->Integral();
508	>	if(Contains(name,"Dibosons")) total+=hist->Integral();
509	>	if(Contains(name,"Z_l_")) total+=hist->Integral();
510	>	if(Contains(name,"Z_c_")) total+=hist->Integral();
511	>	if(Contains(name,"Z_b_")) {
512	>	total+=hist->Integral();
513	>	Zb+=hist->Integral();
514	>	}
515	>	}
516	>
517	>	float purity=Zb/total;
518	>	dout << "Determined a purity of " << 100*purity << " % for the container " << ContainerName << endl;
519	>
520	>	if(ContainerName=="inclusive") PurityInclusive=purity;
521	>	if(ContainerName=="loose") PurityLoose=purity;
522	>	if(ContainerName=="medium") PurityMedium=purity;
523	>	if(ContainerName=="reference") PurityTight=purity;
524	>
525	>	DeleteStack(mcm);
526	>	}
527	>
528	>	string NiceAlphaPrint(float alpha) {
529	>	//get in : 0.3, return 0p3 ...
530	>	stringstream res;
531	>	res << int(alpha) << "p" << int(100*alpha);
532	>	return res.str();
533	>	}
534	>
535	>	ZbResultContainer new_data_mc_agreement_2d(bool gofast, string AlphaVariable, string ContainerName, TCut BTagCut, TCut BTagWgt, float FaceValueToBeTakenAt) {
536	>
537	>
538	>	if(gofast) write_info(__FUNCTION__,"Fast mode activated for "+ContainerName);
539		ZbResultContainer results(ContainerName);
540
541	<	cutWeight=TCut("(weight(weight<1000)(is_data+(!is_data)))")*BTagWgt;
541	>	cutWeight=STDWEIGHT*BTagWgt;
542		LeadingB=BTagCut;
543
544	+
545		gStyle->SetOptFit(0);
546		TCut basecut(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
547	<	// const int nalphacuts=6;
548	<	// float alphacuts[nalphacuts] = {0.1,0.15,0.2,0.25,0.3,0.35};
547	>	const int nalphacuts=6;
548	>	float alphacuts[nalphacuts] = {0.1,0.15,0.2,0.25,0.3,0.35};
549
410	–	// const int nptcuts=5;
411	–	// float ptcuts[nptcuts]={30,50,75,125,1000};
550
551	<	write_info(__FUNCTION__,"JUST TESTING - PLEASE REMOVE THIS LINE");const int nalphacuts=2;float alphacuts[nalphacuts] = {0.1,0.3};const int nptcuts=3;float ptcuts[nptcuts]={30,75,1000};
551	>	const int nptcuts=5;
552	>	float ptcuts[nptcuts]={30,50,75,125,1000};
553
554	+	DeterminePurity(basecut, nptcuts,ptcuts, AlphaVariable, ContainerName, FaceValueToBeTakenAt);
555
556		float MPF_Results[nalphacuts][nptcuts];
557		float MPF_Errors[nalphacuts][nptcuts];
#	Line 419 \| Line 559 \| ZbResultContainer new_data_mc_agreement_
559		float RABS_Errors[nalphacuts][nptcuts];
560
561
562	<	for(int ia=0;ia<nalphacuts;ia++) {
562	>	for(int ia=0;ia<nalphacuts && !gofast;ia++) {
563		for(int ipt=0;ipt<nptcuts-1;ipt++) {
564		stringstream specialcut;
565		float alow,ahigh;
#	Line 433 \| Line 573 \| ZbResultContainer new_data_mc_agreement_
573		ahigh=alphacuts[ia];
574		}
575		cout << specialcut.str() << endl;
576	<	Value MPF_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"mpf","",TCut(basecut && specialcut.str().c_str()),"MPF___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__alpha_"+any2string(alphacuts[ia]), ptcuts[ipt], ptcuts[ipt+1],alow,ahigh);
577	<	Value RABS_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"Zb3010_pfJetGoodPt[0]/pt","",TCut(basecut && specialcut.str().c_str()),"Rabs___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__alpha_"+any2string(alphacuts[ia]), ptcuts[ipt], ptcuts[ipt+1],alow,ahigh);
576	>	Value MPF_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"mpf","",TCut(basecut && specialcut.str().c_str()),"MPF___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__alpha_"+any2string(NiceAlphaPrint(alphacuts[ia])), ptcuts[ipt], ptcuts[ipt+1],alow,ahigh);
577	>	Value RABS_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"ZbCHS3010_pfJetGoodPt[0]/pt","",TCut(basecut && specialcut.str().c_str()),"Rabs___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__alpha_"+any2string(NiceAlphaPrint(alphacuts[ia])), ptcuts[ipt], ptcuts[ipt+1],alow,ahigh);
578
579		MPF_Results[ia][ipt]=MPF_data_over_mc.getValue();
580		MPF_Errors[ia][ipt]=MPF_data_over_mc.getError();
#	Line 468 \| Line 608 \| ZbResultContainer new_data_mc_agreement_
608		int rabs_counter=0;
609		int mpf_counter=0;
610
611	<	for(int ia=0;ia<nalphacuts;ia++) {
611	>	for(int ia=0;ia<nalphacuts && !gofast;ia++) {
612		cout << "Setting a point for an alpha cut at " << alphacuts[ia] << " with value " << MPF_Results[ia][ipt] << " (MPF) and " << RABS_Results[ia][ipt] << " (RABS)" << endl;
613		if(MPF_Results[ia][ipt]==MPF_Results[ia][ipt] && MPF_Errors[ia][ipt]==MPF_Errors[ia][ipt]) { // remove nan's
614		mpf_gr->SetPoint(mpf_counter,alphacuts[ia],MPF_Results[ia][ipt]);
#	Line 491 \| Line 631 \| ZbResultContainer new_data_mc_agreement_
631		}
632
633		stringstream specialcut;
634	<	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (" << AlphaVariable << "<0.3))";
635	<	Value MPF_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"mpf","",TCut(basecut && specialcut.str().c_str()),"MPF___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__INCLUSIVEalpha_0.3",ptcuts[ipt],ptcuts[ipt+1],0,0.3);
636	<	Value RABS_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"Zb3010_pfJetGoodPt[0]/pt","",TCut(basecut && specialcut.str().c_str()),"Rabs___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__INCLUSIVEalpha_0.3",ptcuts[ipt],ptcuts[ipt+1],0,0.3);
634	>	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (" << AlphaVariable << "<" << FaceValueToBeTakenAt << "))";
635	>	Value MPF_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"mpf","",TCut(basecut && specialcut.str().c_str()),"MPF___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__INCLUSIVEalpha_"+any2string(FaceValueToBeTakenAt),ptcuts[ipt],ptcuts[ipt+1],0,FaceValueToBeTakenAt);
636	>	Value RABS_data_over_mc = get_Zb_data_over_mc(ContainerName,results,"ZbCHS3010_pfJetGoodPt[0]/pt","",TCut(basecut && specialcut.str().c_str()),"Rabs___pt_"+any2string(ptcuts[ipt])+"_"+any2string(ptcuts[ipt+1])+"__INCLUSIVEalpha_"+any2string(FaceValueToBeTakenAt),ptcuts[ipt],ptcuts[ipt+1],0,FaceValueToBeTakenAt);
637		cout << "About to add a point in pt (" << (ptcuts[ipt]+ptcuts[ipt+1])/2 << " ) to face value : " << MPF_data_over_mc << " (MPF) , " << RABS_data_over_mc << "(RABS)" << endl;
638		MPF_FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1])/2,MPF_data_over_mc.getValue());
639		MPF_FaceValueAtPoint3->SetPointError(ipt,0,MPF_data_over_mc.getError());
640		RABS_FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1])/2,RABS_data_over_mc.getValue());
641		RABS_FaceValueAtPoint3->SetPointError(ipt,0,RABS_data_over_mc.getError());
642
503	–	can->cd();
504	–	mpf_gr->SetMarkerStyle(21);
505	–	mpf_gr->SetMarkerColor(kBlue);
506	–	mpf_gr->Draw("AP");
507	–	mpf_gr->Fit("pol1","");
508	–	mpf_gr->GetXaxis()->SetTitle("#alpha");
509	–	mpf_gr->GetXaxis()->CenterTitle();
510	–	mpf_gr->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
511	–	mpf_gr->GetYaxis()->CenterTitle();
512	–	TF1 mpf_pol=(TF1)mpf_gr->GetFunction("pol1");
513	–	mpf_pol->SetLineWidth(0);
514	–	TF1 mpf_pol_complete = new TF1("mpf_pol_complete","[0]+[1]x");
515	–	mpf_pol_complete->SetParameters(mpf_pol->GetParameters());
516	–	mpf_pol_complete->SetLineWidth(1);
517	–	mpf_pol_complete->SetLineColor(kBlue);
518	–	float min,max;
519	–	FindMinMax(mpf_gr,min,max);
520	–	TH1F *histo = new TH1F("histo","histo",1,0,0.4);
521	–	histo->GetXaxis()->SetTitle("#alpha");
522	–	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
523	–	histo->GetXaxis()->CenterTitle();
524	–	histo->GetYaxis()->CenterTitle();
525	–	histo->GetYaxis()->SetTitleOffset(1.3);
526	–	histo->SetStats(0);
527	–	histo->GetXaxis()->SetRangeUser(0,0.4);
528	–	histo->GetYaxis()->SetRangeUser(min,max);
529	–	mpf_gr->GetYaxis()->SetRangeUser(min,max);
530	–
531	–	TPolyLine *mpf_fit_uncert = GetFitUncertaintyShape(mpf_gr, "pol1", min, max,0.0,0.4);
532	–	mpf_pol->SetLineColor(TColor::GetColor("#0101DF"));
533	–	mpf_fit_uncert->SetFillColor(TColor::GetColor("#5353FC"));
534	–	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
535	–	histo->Draw();
536	–	mpf_fit_uncert->Draw("F");
537	–	mpf_fit_uncert->Draw("");
538	–	mpf_gr->Draw("P");
539	–	mpf_pol_complete->Draw("same");
540	–
541	–	float mpf_a=mpf_pol->GetParameter(0);
542	–	float mpf_b=mpf_pol->GetParameter(1);
543	–	MPF_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),mpf_a);
544	–	MPF_FinalGraph->SetPointError(ipt,0,mpf_pol->GetParError(0));
545	–
546	–	MPF_ExtraPolatedResults[ipt]=mpf_a;
547	–
548	–	stringstream mpf_info;
549	–	mpf_info << "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (MPF) }}{Extrapolated value: " << std::setprecision(4) << mpf_a << "}}{#chi^{2}/NDF : " << mpf_pol->GetChisquare() << " / " << mpf_pol->GetNDF() << "}";
550	–
551	–	TText *mpf_mark = write_title(mpf_info.str());
552	–	mpf_mark->SetX(0.75);
553	–	mpf_mark->SetTextSize(0.03);
554	–	mpf_mark->SetY(0.75);
555	–	mpf_mark->Draw();
643
557	–	string filenamebkp=filename.str();
558	–	filename << "__MPF";
559	–	DrawPrelim();
560	–	CompleteSave(can,ContainerName+"/"+filename.str());
561	–	cout << "MPF : " << mpf_a << " + " << mpf_b << " * alpha " << endl;
644
645	<	filename.str("");
645	>	if(!gofast) {
646	>	can->cd();
647	>	mpf_gr->SetMarkerStyle(21);
648	>	mpf_gr->SetMarkerColor(kBlue);
649	>	mpf_gr->Draw("AP");
650	>	mpf_gr->Fit("pol1","");
651	>	mpf_gr->GetXaxis()->SetTitle("#alpha");
652	>	mpf_gr->GetXaxis()->CenterTitle();
653	>	mpf_gr->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
654	>	mpf_gr->GetYaxis()->CenterTitle();
655	>	TF1 mpf_pol=(TF1)mpf_gr->GetFunction("pol1");
656	>	mpf_pol->SetLineWidth(0);
657	>	TF1 mpf_pol_complete = new TF1("mpf_pol_complete","[0]+[1]x");
658	>	mpf_pol_complete->SetParameters(mpf_pol->GetParameters());
659	>	mpf_pol_complete->SetLineWidth(1);
660	>	mpf_pol_complete->SetLineColor(kBlue);
661	>
662	>	float min,max;
663	>	FindMinMax(mpf_gr,min,max);
664	>	TH1F *histo = new TH1F("histo","histo",1,0,0.4);
665	>	histo->GetXaxis()->SetTitle("#alpha");
666	>	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
667	>	histo->GetXaxis()->CenterTitle();
668	>	histo->GetYaxis()->CenterTitle();
669	>	histo->GetYaxis()->SetTitleOffset(1.3);
670	>	histo->SetStats(0);
671	>	histo->GetXaxis()->SetRangeUser(0,0.4);
672	>	histo->GetYaxis()->SetRangeUser(min,max);
673	>	mpf_gr->GetYaxis()->SetRangeUser(min,max);
674	>
675	>	TPolyLine *mpf_fit_uncert = GetFitUncertaintyShape(mpf_gr, "pol1", min, max,0.0,0.4);
676	>	mpf_pol->SetLineColor(TColor::GetColor("#0101DF"));
677	>	mpf_fit_uncert->SetFillColor(TColor::GetColor("#5353FC"));
678	>	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
679	>	histo->Draw();
680	>	mpf_fit_uncert->Draw("F");
681	>	mpf_fit_uncert->Draw("");
682	>	mpf_gr->Draw("P");
683	>	mpf_pol_complete->Draw("same");
684	>
685	>	float mpf_a=mpf_pol->GetParameter(0);
686	>	float mpf_b=mpf_pol->GetParameter(1);
687	>	MPF_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),mpf_a);
688	>	MPF_FinalGraph->SetPointError(ipt,0,mpf_pol->GetParError(0));
689	>
690	>	MPF_ExtraPolatedResults[ipt]=mpf_a;
691	>
692	>	stringstream mpf_info;
693	>	mpf_info << "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (MPF) }}{Extrapolated value: " << std::setprecision(4) << mpf_a << "}}{#chi^{2}/NDF : " << mpf_pol->GetChisquare() << " / " << mpf_pol->GetNDF() << "}";
694	>
695	>	TText *mpf_mark = write_title(mpf_info.str());
696	>	mpf_mark->SetX(0.75);
697	>	mpf_mark->SetTextSize(0.03);
698	>	mpf_mark->SetY(0.75);
699	>	mpf_mark->Draw();
700	>
701	>	string filenamebkp=filename.str();
702	>	filename << "__MPF";
703	>	DrawPrelim();
704	>	CompleteSave(can,ContainerName+"/"+filename.str());
705	>	cout << "MPF : " << mpf_a << " + " << mpf_b << " * alpha " << endl;
706	>
707	>	filename.str("");
708	>
709	>	rabs_gr->SetMarkerStyle(21);
710	>	rabs_gr->SetMarkerColor(kRed);
711	>	rabs_gr->Draw("AP");
712	>	rabs_gr->Fit("pol1","");
713	>	rabs_gr->GetXaxis()->SetTitle("#alpha");
714	>	rabs_gr->GetXaxis()->CenterTitle();
715	>	rabs_gr->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
716	>	rabs_gr->GetYaxis()->CenterTitle();
717	>	TF1 rabs_pol=(TF1)rabs_gr->GetFunction("pol1");
718	>	rabs_pol->SetLineWidth(0);
719	>	TF1 rabs_pol_complete = new TF1("rabs_pol_complete","[0]+[1]x");
720	>	rabs_pol_complete->SetParameters(rabs_pol->GetParameters());
721	>	rabs_pol_complete->SetLineColor(kRed);
722	>	rabs_pol_complete->SetLineWidth(1);
723	>	FindMinMax(rabs_gr,min,max);
724	>	rabs_gr->GetYaxis()->SetRangeUser(min,max);
725	>	histo->GetYaxis()->SetRangeUser(min,max);
726	>	TPolyLine *rabs_fit_uncert = GetFitUncertaintyShape(rabs_gr, "pol1", min, max,0.0,0.4);
727	>	rabs_pol->SetLineColor(TColor::GetColor("#FA5858"));
728	>	rabs_pol->SetFillColor(TColor::GetColor("#FA5858"));
729	>	rabs_fit_uncert->SetLineColor(TColor::GetColor("#FA5858"));
730	>	rabs_fit_uncert->SetFillColor(TColor::GetColor("#FA5858"));
731	>	histo->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
732	>	histo->Draw();
733	>	rabs_fit_uncert->Draw("F");
734	>	rabs_fit_uncert->Draw("");
735	>	rabs_gr->Draw("P");
736	>	rabs_pol_complete->Draw("same");
737	>
738	>	float rabs_a=rabs_pol->GetParameter(0);
739	>	float rabs_b=rabs_pol->GetParameter(1);
740	>	RABS_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),rabs_a);
741	>	RABS_FinalGraph->SetPointError(ipt,0,rabs_pol->GetParError(0));
742	>
743	>	cout << "!+!+!+!+!+!+!+!+!+ Just added a point to the final plot : " << rabs_a << " +/- " << rabs_pol->GetParError(0) << endl;
744	>
745	>	stringstream rabs_info;
746	>	rabs_info << "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (R_{abs}) }}{Extrapolated value: " << std::setprecision(4) << rabs_a << "}}{#chi^{2}/NDF : " << rabs_pol->GetChisquare() << " / " << rabs_pol->GetNDF() << "}";
747	>	TText *rabs_mark = write_title(rabs_info.str());
748	>	rabs_mark->SetX(0.75);
749	>	rabs_mark->SetTextSize(0.03);
750	>	rabs_mark->SetY(0.75);
751	>	rabs_mark->Draw();
752	>
753	>	RABS_ExtraPolatedResults[ipt]=rabs_a;
754	>	filename << filenamebkp << "__RABS";
755	>	DrawPrelim();
756	>	CompleteSave(can,ContainerName+"/"+filename.str());
757	>	cout << "RABS : " << rabs_a << " + " << rabs_b << " * alpha " << endl;
758	>	}
759	>	}
760	>
761	>	float MPFResult;
762	>	float MPFResultError;
763	>
764	>	float RabsResult;
765	>	float RabsResultError;
766	>
767	>	if(!gofast) {
768	>	MPF_FinalGraph->SetMarkerStyle(21);
769	>	MPF_FinalGraph->SetMarkerColor(kBlue);
770	>	MPF_FinalGraph->Fit("pol0",""); // quiet, use minos
771	>
772	>	TF1 mpf_pol0=(TF1)MPF_FinalGraph->GetFunction("pol0");
773	>	TF1 mpf_pol0_complete = new TF1("mpf_pol0_complete","[0]+[1]x");
774	>	mpf_pol0_complete->SetLineWidth(1);
775	>	mpf_pol0_complete->SetLineColor(kBlue);
776	>	mpf_pol0->SetLineColor(kBlue);
777	>	MPF_FinalGraph->Draw("AP");
778	>	MPF_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
779	>	MPF_FinalGraph->GetXaxis()->CenterTitle();
780	>	MPF_FinalGraph->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
781	>	MPF_FinalGraph->GetYaxis()->CenterTitle();
782	>	MPF_FinalGraph->Draw("AP");
783	>	mpf_pol0_complete->Draw("same");
784	>
785	>	stringstream mpf_result;
786	>	mpf_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/mpf_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << mpf_pol0->GetParError(0)/(mpf_pol0->GetParameter(0)*mpf_pol0->GetParameter(0)) << "}{ (MPF)}";
787	>	TText *rmpf_final_mark = write_title(mpf_result.str());
788	>	rmpf_final_mark->SetX(0.75);
789	>	rmpf_final_mark->SetY(0.75);
790	>	rmpf_final_mark->SetTextSize(0.03);
791	>	rmpf_final_mark->Draw();
792	>	DrawPrelim();
793	>	MPFResult=1/mpf_pol0->GetParameter(0);
794	>	MPFResultError=mpf_pol0->GetParError(0)/(mpf_pol0->GetParameter(0)*mpf_pol0->GetParameter(0));
795
796	<	rabs_gr->SetMarkerStyle(21);
797	<	rabs_gr->SetMarkerColor(kRed);
798	<	rabs_gr->Draw("AP");
799	<	rabs_gr->Fit("pol1","");
800	<	rabs_gr->GetXaxis()->SetTitle("#alpha");
801	<	rabs_gr->GetXaxis()->CenterTitle();
802	<	rabs_gr->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
803	<	rabs_gr->GetYaxis()->CenterTitle();
804	<	TF1 rabs_pol=(TF1)rabs_gr->GetFunction("pol1");
805	<	rabs_pol->SetLineWidth(0);
806	<	TF1 rabs_pol_complete = new TF1("rabs_pol_complete","[0]+[1]x");
807	<	rabs_pol_complete->SetParameters(rabs_pol->GetParameters());
808	<	rabs_pol_complete->SetLineColor(kRed);
809	<	rabs_pol_complete->SetLineWidth(1);
810	<	FindMinMax(rabs_gr,min,max);
811	<	rabs_gr->GetYaxis()->SetRangeUser(min,max);
812	<	histo->GetYaxis()->SetRangeUser(min,max);
813	<	TPolyLine *rabs_fit_uncert = GetFitUncertaintyShape(rabs_gr, "pol1", min, max,0.0,0.4);
814	<	rabs_pol->SetLineColor(TColor::GetColor("#FA5858"));
815	<	rabs_pol->SetFillColor(TColor::GetColor("#FA5858"));
816	<	rabs_fit_uncert->SetLineColor(TColor::GetColor("#FA5858"));
817	<	rabs_fit_uncert->SetFillColor(TColor::GetColor("#FA5858"));
818	<	histo->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
588	<	histo->Draw();
589	<	rabs_fit_uncert->Draw("F");
590	<	rabs_fit_uncert->Draw("");
591	<	rabs_gr->Draw("P");
592	<	rabs_pol_complete->Draw("same");
593	<
594	<
595	<	float rabs_a=rabs_pol->GetParameter(0);
596	<	float rabs_b=rabs_pol->GetParameter(1);
597	<	RABS_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),rabs_a);
598	<	RABS_FinalGraph->SetPointError(ipt,0,rabs_pol->GetParError(0));
599	<
600	<	cout << "!+!+!+!+!+!+!+!+!+ Just added a point to the final plot : " << rabs_a << " +/- " << rabs_pol->GetParError(0) << endl;
601	<
602	<	stringstream rabs_info;
603	<	rabs_info << "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (R_{abs}) }}{Extrapolated value: " << std::setprecision(4) << rabs_a << "}}{#chi^{2}/NDF : " << rabs_pol->GetChisquare() << " / " << rabs_pol->GetNDF() << "}";
604	<	TText *rabs_mark = write_title(rabs_info.str());
605	<	rabs_mark->SetX(0.75);
606	<	rabs_mark->SetTextSize(0.03);
607	<	rabs_mark->SetY(0.75);
608	<	rabs_mark->Draw();
796	>	stringstream filename2;
797	>	filename2 << "Extrapolation/FINAL_RESULT_MPF";
798	>	CompleteSave(can,ContainerName+"/"+filename2.str());
799	>
800	>	RABS_FinalGraph->SetMarkerStyle(21);
801	>	RABS_FinalGraph->SetMarkerStyle(21);
802	>	RABS_FinalGraph->SetMarkerColor(kRed);
803	>	RABS_FinalGraph->Fit("pol0","QE"); // quiet, use minos
804	>	RABS_FinalGraph->Draw("AP");
805	>	RABS_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
806	>	RABS_FinalGraph->GetXaxis()->CenterTitle();
807	>	RABS_FinalGraph->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
808	>	RABS_FinalGraph->GetYaxis()->CenterTitle();
809	>	RABS_FinalGraph->Draw("AP");
810	>	TF1 rabs_pol0=(TF1)RABS_FinalGraph->GetFunction("pol0");
811	>	stringstream rabs_result;
812	>	rabs_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/rabs_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << rabs_pol0->GetParError(0)/(rabs_pol0->GetParameter(0)*rabs_pol0->GetParameter(0)) << "}{ (R_{abs})}";;
813	>	TText *rabs_final_mark = write_title(rabs_result.str());
814	>	rabs_final_mark->SetX(0.75);
815	>	rabs_final_mark->SetY(0.75);
816	>	rabs_final_mark->SetTextSize(0.03);
817	>	rabs_final_mark->Draw();
818	>	DrawPrelim();
819
820	+	RabsResult=1/rabs_pol0->GetParameter(0);
821	+	RabsResultError=rabs_pol0->GetParError(0)/(rabs_pol0->GetParameter(0)*rabs_pol0->GetParameter(0));
822
823	<	RABS_ExtraPolatedResults[ipt]=rabs_a;
824	<	filename << filenamebkp << "__RABS";
825	<	DrawPrelim();
614	<	CompleteSave(can,ContainerName+"/"+filename.str());
615	<	cout << "RABS : " << rabs_a << " + " << rabs_b << " * alpha " << endl;
823	>	filename2.str("");
824	>	filename2 << "Extrapolation/FINAL_RESULT_RABS";
825	>	CompleteSave(can,ContainerName+"/"+filename2.str());
826		}
827	<
618	<	MPF_FinalGraph->SetMarkerStyle(21);
619	<	MPF_FinalGraph->SetMarkerColor(kBlue);
620	<	MPF_FinalGraph->Fit("pol0",""); // quiet, use minos
621	<
622	<	TF1 mpf_pol0=(TF1)MPF_FinalGraph->GetFunction("pol0");
623	<	TF1 mpf_pol0_complete = new TF1("mpf_pol0_complete","[0]+[1]x");
624	<	mpf_pol0_complete->SetLineWidth(1);
625	<	mpf_pol0_complete->SetLineColor(kBlue);
626	<	mpf_pol0->SetLineColor(kBlue);
627	<	MPF_FinalGraph->Draw("AP");
628	<	MPF_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
629	<	MPF_FinalGraph->GetXaxis()->CenterTitle();
630	<	MPF_FinalGraph->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
631	<	MPF_FinalGraph->GetYaxis()->CenterTitle();
632	<	MPF_FinalGraph->Draw("AP");
633	<	mpf_pol0_complete->Draw("same");
634	<
635	<	stringstream mpf_result;
636	<	mpf_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/mpf_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << mpf_pol0->GetParError(0)/(mpf_pol0->GetParameter(0)*mpf_pol0->GetParameter(0)) << "}{ (MPF)}";
637	<	TText *rmpf_final_mark = write_title(mpf_result.str());
638	<	rmpf_final_mark->SetX(0.75);
639	<	rmpf_final_mark->SetY(0.75);
640	<	rmpf_final_mark->SetTextSize(0.03);
641	<	rmpf_final_mark->Draw();
642	<	DrawPrelim();
643	<	float MPFResult=1/mpf_pol0->GetParameter(0);
644	<	float MPFResultError=mpf_pol0->GetParError(0)/(mpf_pol0->GetParameter(0)*mpf_pol0->GetParameter(0));
645	<
646	<	stringstream filename2;
647	<	filename2 << "Extrapolation/FINAL_RESULT_MPF";
648	<	CompleteSave(can,ContainerName+"/"+filename2.str());
649	<
650	<
651	<	RABS_FinalGraph->SetMarkerStyle(21);
652	<	RABS_FinalGraph->SetMarkerStyle(21);
653	<	RABS_FinalGraph->SetMarkerColor(kRed);
654	<	RABS_FinalGraph->Fit("pol0","QE"); // quiet, use minos
655	<	RABS_FinalGraph->Draw("AP");
656	<	RABS_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
657	<	RABS_FinalGraph->GetXaxis()->CenterTitle();
658	<	RABS_FinalGraph->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
659	<	RABS_FinalGraph->GetYaxis()->CenterTitle();
660	<	RABS_FinalGraph->Draw("AP");
661	<	TF1 rabs_pol0=(TF1)RABS_FinalGraph->GetFunction("pol0");
662	<	stringstream rabs_result;
663	<	rabs_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/rabs_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << rabs_pol0->GetParError(0)/(rabs_pol0->GetParameter(0)*rabs_pol0->GetParameter(0)) << "}{ (R_{abs})}";;
664	<	TText *rabs_final_mark = write_title(rabs_result.str());
665	<	rabs_final_mark->SetX(0.75);
666	<	rabs_final_mark->SetY(0.75);
667	<	rabs_final_mark->SetTextSize(0.03);
668	<	rabs_final_mark->Draw();
669	<	DrawPrelim();
670	<
671	<	float RabsResult=1/rabs_pol0->GetParameter(0);
672	<	float RabsResultError=rabs_pol0->GetParError(0)/(rabs_pol0->GetParameter(0)*rabs_pol0->GetParameter(0));
673	<
674	<	filename2.str("");
675	<	filename2 << "Extrapolation/FINAL_RESULT_RABS";
676	<	CompleteSave(can,ContainerName+"/"+filename2.str());
677	<
827	>
828		can->cd();
829		MPF_FaceValueAtPoint3->SetMarkerStyle(21);
830		MPF_FaceValueAtPoint3->SetMarkerColor(kBlue);
#	Line 682 \| Line 832 \| ZbResultContainer new_data_mc_agreement_
832		MPF_FaceValueAtPoint3->Draw("AP");
833		MPF_FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
834		MPF_FaceValueAtPoint3->GetXaxis()->CenterTitle();
835	<	MPF_FaceValueAtPoint3->GetYaxis()->SetTitle("<MPF>__{data}/<MPF>_{mc} for #alpha<0.3");
835	>	MPF_FaceValueAtPoint3->GetYaxis()->SetTitle(((string)"<MPF>__{data}/<MPF>_{mc} for #alpha<"+any2string(FaceValueToBeTakenAt)).c_str());
836		MPF_FaceValueAtPoint3->GetYaxis()->CenterTitle();
837		MPF_FaceValueAtPoint3->Draw("AP");
838		TF1 faceval_pol0=(TF1)MPF_FaceValueAtPoint3->GetFunction("pol0");
#	Line 692 \| Line 842 \| ZbResultContainer new_data_mc_agreement_
842		TF1 faceval_pol0_complete = new TF1("faceval_pol0_complete","[0]+[1]x");
843		faceval_pol0_complete->SetParameters(faceval_pol0->GetParameters());
844		stringstream faceval_result;
845	<	faceval_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/faceval_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << faceval_pol0->GetParError(0)/(faceval_pol0->GetParameter(0)*faceval_pol0->GetParameter(0)) << "}{ (MPF at #alpha<0.3)}";;
845	>	faceval_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/faceval_pol0->GetParameter(0) << " #pm " << std::setprecision(2) << faceval_pol0->GetParError(0)/(faceval_pol0->GetParameter(0)*faceval_pol0->GetParameter(0)) << "}{ (MPF at #alpha<" << FaceValueToBeTakenAt << ")}";;
846		TText *faceval_final_mark = write_title(faceval_result.str());
847		faceval_final_mark->SetX(0.75);
848		faceval_final_mark->SetY(0.75);
#	Line 703 \| Line 853 \| ZbResultContainer new_data_mc_agreement_
853		float FaceValResult=1/faceval_pol0->GetParameter(0);
854		float FaceValResultError=faceval_pol0->GetParError(0)/(faceval_pol0->GetParameter(0)*faceval_pol0->GetParameter(0));
855
856	+	stringstream filename2;
857		filename2.str("");
858		filename2 << "Extrapolation/FINAL_RESULT_MPF_FaceValp3";
859		CompleteSave(can,ContainerName+"/"+filename2.str());
#	Line 714 \| Line 865 \| ZbResultContainer new_data_mc_agreement_
865		RABS_FaceValueAtPoint3->Draw("AP");
866		RABS_FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
867		RABS_FaceValueAtPoint3->GetXaxis()->CenterTitle();
868	<	RABS_FaceValueAtPoint3->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc} for #alpha<0.3");
868	>	RABS_FaceValueAtPoint3->GetYaxis()->SetTitle(((string)"<R_{abs}>_{data}/<R_{abs}>_{mc} for #alpha<"+any2string(FaceValueToBeTakenAt)).c_str());
869		RABS_FaceValueAtPoint3->GetYaxis()->CenterTitle();
870		RABS_FaceValueAtPoint3->Draw("AP");
871		TF1 faceval_pol0RABS=(TF1)RABS_FaceValueAtPoint3->GetFunction("pol0");
#	Line 724 \| Line 875 \| ZbResultContainer new_data_mc_agreement_
875		faceval_pol0RABS->SetLineColor(kRed);
876		faceval_pol0RABS_complete->Draw("same");
877		stringstream RABSfaceval_result;
878	<	RABSfaceval_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/faceval_pol0RABS->GetParameter(0) << " #pm " << std::setprecision(2) << faceval_pol0RABS->GetParError(0)/(faceval_pol0RABS->GetParameter(0)*faceval_pol0RABS->GetParameter(0)) << "}{ (R_{abs} at #alpha<0.3)}";;
878	>	RABSfaceval_result << "#splitline{C_{abs}= " << std::setprecision(5) << 1/faceval_pol0RABS->GetParameter(0) << " #pm " << std::setprecision(2) << faceval_pol0RABS->GetParError(0)/(faceval_pol0RABS->GetParameter(0)*faceval_pol0RABS->GetParameter(0)) << "}{ (R_{abs} at #alpha<" << FaceValueToBeTakenAt << ")}";;
879		TText *RABSfaceval_final_mark = write_title(RABSfaceval_result.str());
880		RABSfaceval_final_mark->SetX(0.75);
881		RABSfaceval_final_mark->SetY(0.75);
#	Line 740 \| Line 891 \| ZbResultContainer new_data_mc_agreement_
891		CompleteSave(can,ContainerName+"/"+filename2.str());
892
893		cout << "FINAL RESULTS: " << endl;
894	<	cout << "MPF : " << MPFResult << " +/- " << MPFResultError << endl;
895	<	cout << "Rabs: " << RabsResult << " +/- " << RabsResultError << endl;
896	<	cout << "Face value at 0.3 :" << FaceValResult << " +/- " << FaceValResultError << " (MPF) " << endl;
897	<	cout << "Face value at 0.3 :" << RABSFaceValResult << " +/- " << RABSFaceValResultError << " (RABS) " << endl;
894	>	if(!gofast) cout << "MPF : " << MPFResult << " +/- " << MPFResultError << endl;
895	>	if(!gofast) cout << "Rabs: " << RabsResult << " +/- " << RabsResultError << endl;
896	>	cout << "Face value at " << FaceValueToBeTakenAt << " :" << FaceValResult << " +/- " << FaceValResultError << " (MPF) " << endl;
897	>	cout << "Face value at " << FaceValueToBeTakenAt << " :" << RABSFaceValResult << " +/- " << RABSFaceValResultError << " (RABS) " << endl;
898
899		delete can;
900
901		results.MPF_Result_FaceValue=Value(FaceValResult,FaceValResultError);
902		results.Rabs_Result_FaceValue=Value(RABSFaceValResult,RABSFaceValResultError);
903	<	results.MPF_Result_Extrapolated=Value(MPFResult,MPFResultError);
904	<	results.Rabs_Result_Extrapolated=Value(RabsResult,RabsResultError);
903	>	if(!gofast) results.MPF_Result_Extrapolated=Value(MPFResult,MPFResultError);
904	>	if(!gofast) results.Rabs_Result_Extrapolated=Value(RabsResult,RabsResultError);
905		cout << results << endl;
906		return results;
907		}
908
909
910	<	void TEST_DoPUStudy(string identifier) {
910	>	void DoPUStudy(string identifier) {
911
912		float numVtxcuts[7]={0,10,15,20};
913
#	Line 775 \| Line 926 \| void TEST_DoPUStudy(string identifier) {
926		cout << numVtxcuts[i-1] << ";";
927		}
928
929	+	sSpecialCut << " && " << identifier << "_alpha<0.3";
930		TCut SpecialCut(sSpecialCut.str().c_str());
931
932
933	<	malpha[i] = allsamples.Draw("malpha"+any2string(i),"mpf",1000,0,10,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,mc,luminosity);
934	<	dalpha[i] = allsamples.Draw("dalpha"+any2string(i),"mpf",1000,0,10,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,data,luminosity);
933	>	malpha[i] = allsamples.Draw("malpha"+any2string(i),"mpf",1000,0,2,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,mc,luminosity);
934	>	dalpha[i] = allsamples.Draw("dalpha"+any2string(i),"mpf",1000,0,2,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,data,luminosity);
935
936
937		float a=dalpha[i]->GetMean();
#	Line 790 \| Line 942 \| void TEST_DoPUStudy(string identifier) {
942		float error = TMath::Sqrt( (1/(bb)) (dada) + ((aa)/(bb))db*db);
943
944		cout << malpha[i]->GetMean() << ";" << malpha[i]->GetMeanError() << ";" << dalpha[i]->GetMean() << ";" << dalpha[i]->GetMeanError() << ";" << malpha[i]->Integral()<<";"<<dalpha[i]->Integral() << ";"<<factor << ";" << error << endl;
793	–	malpha[i]->SetLineColor(i+1);
794	–	if(i==0) malpha[i]->Draw("histo");
795	–	else malpha[i]->Draw("histo,same");
796	–	}
797	–
798	–	CompleteSave(can,"AlphaOverview_"+identifier);
799	–	delete can;
800	–
801	–	/* for(int i=1;i<8;i++) {
945		delete malpha[i];
946		delete dalpha[i];
947		}
948	<	*/
948	>
949	>	delete can;
950
951		}
952
953	<	void TEST_ScenarioComparison() {
953	>	void ScenarioComparison() {
954		//very dumb way to do this.
955	+
956		TGraphAsymmErrors *gr[5];
957		TF1 *fit[5];
958	<	bool dofit=false;
958	>	bool dofit=true;
959
960		TCanvas *can = new TCanvas("can","can",500,500);
961
#	Line 818 \| Line 963 \| void TEST_ScenarioComparison() {
963		leg->SetBorderSize(0);
964		leg->SetFillColor(kWhite);
965
966	<	for(int i=0;i<5;i++) {
966	>	for(int i=0;i<6;i++) {
967		gr[i] = new TGraphAsymmErrors();
968	<	double x[4] = {5,12.5,17.5,20};
969	<	double dxl[4] = {5,2.5,2.5,2.5};
968	>	double x[4] = {5,12.5,17.5,25};
969	>	double dxl[4] = {5,2.5,2.5,5};
970		double dxh[4] = {5,2.5,2.5,10};
971
972		string name="";
973		string color="";
974
975		if(i==0) {
976	<	double y[4] = {1.008,1.013,1.003,1.007};
977	<	double dyl[4] = {0.011,0.009,0.013,0.020};
978	<	double dyh[4] = {0.011,0.009,0.013,0.020};
976	>	double y[4] = {1.01935,1.01558,0.996438,1.01901};
977	>	double dyl[4] = {0.0132973,0.0131661,0.017629,0.0319767};
978	>	double dyh[4] = {0.0132973,0.0131661,0.017629,0.0319767};
979		name="30/30";
980		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
981		color="#a01d99";
982		}
983		if(i==1) {
984	<	double y[4] = {1.008+0.001,1.013+0.001,1.003+0.001,1.007+0.001};
985	<	double dyl[4] = {0.011,0.009,0.013,0.020};
986	<	double dyh[4] = {0.011,0.009,0.013,0.020};
984	>	double y[4] = {1.00577,1.00706,0.991029,0.954098};
985	>	double dyl[4] = {0.0135575,0.0136548,0.0170778,0.0263919};
986	>	double dyh[4] = {0.0135575,0.0136548,0.0170778,0.0263919};
987		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
988		name="30/15";
989		color="#2464bc";
990		}
991		if(i==2) {
992	<	double y[4] = {1.014,1.026,1.014,1.059};
993	<	double dyl[4] = {0.011,0.010,0.014,0.023};
994	<	double dyh[4] = {0.011,0.010,0.014,0.023};
992	>	double y[4] = {1.023,1.00648,1.02217,1.03181};
993	>	double dyl[4] = {0.0157372,0.0187778,0.0261512,0.0422295};
994	>	double dyh[4] = {0.0157372,0.0187778,0.0261512,0.0422295};
995		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
996		name="15/15";
997		color="#f49230";
998		}
999		if(i==3) {
1000	<	double y[4] = {1.015,1.020,1.008,1.050};
1001	<	double dyl[4] = {0.011,0.010,0.014,0.022};
1002	<	double dyh[4] = {0.011,0.010,0.014,0.022};
1000	>	double y[4] = {1.02029,1.01338,0.999769,1.00357};
1001	>	double dyl[4] = {0.0129509,0.0120127,0.0160537,0.0268972};
1002	>	double dyh[4] = {0.0129509,0.0120127,0.0160537,0.0268972};
1003		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1004	<	name="20/15";
1004	>	name="CHS 30/30";
1005		color="#c72f3c";
1006		}
1007		if(i==4) {
1008	<	double y[4] = {1.015+0.001,1.020+0.001,1.008+0.001,1.050+0.001};
1009	<	double dyl[4] = {0.011,0.010,0.014,0.022};
1010	<	double dyh[4] = {0.011,0.010,0.014,0.022};
1008	>	double y[4] = {1.00274,1.01056,0.993214,0.997963};
1009	>	double dyl[4] = {0.0127179,0.0140344,0.0162814,0.0278914};
1010	>	double dyh[4] = {0.0127179,0.0140344,0.0162814,0.0278914};
1011		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1012	<	name="20/20";
1012	>	name="CHS 30/15";
1013		color="#7cb433";
1014		}
1015	<
1016	<
1015	>	if(i==5) {
1016	>	double y[4] = {1.04103,1.03726,1.04649,1.05383};
1017	>	double dyl[4] = {0.0156465,0.0168194,0.0231982,0.0376114};
1018	>	double dyh[4] = {0.0156465,0.0168194,0.0231982,0.0376114};
1019	>	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1020	>	name="CHS 15/15";
1021	>	color="#DF013A";
1022	>	}
1023	>
1024	>
1025		gr[i]->SetTitle();
1026		gr[i]->GetXaxis()->SetTitle("N(Vtx)");
1027		gr[i]->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
#	Line 896 \| Line 1049 \| void TEST_ScenarioComparison() {
1049
1050		}
1051
1052	<	void TEST_ScenarioComparisonInclusive() {
1052	>	void ScenarioComparisonInclusive() {
1053		//dumbest way ever to do this. but ok we only need to do it once.
1054		TGraphAsymmErrors *gr[5];
1055		TF1 *fit[5];
#	Line 904 \| Line 1057 \| void TEST_ScenarioComparisonInclusive()
1057
1058		TCanvas *can = new TCanvas("can","can",500,500);
1059
1060	<	TLegend *leg = new TLegend(0.2,0.2,0.5,0.45);
1060	>	TLegend *leg = new TLegend(0.2,0.2,0.65,0.45);
1061		leg->SetBorderSize(0);
1062		leg->SetFillColor(kWhite);
1063
1064	<	for(int i=0;i<5;i++) {
1064	>	for(int i=0;i<6;i++) {
1065		gr[i] = new TGraphAsymmErrors();
1066	<	double x[4] = {5,12.5,17.5,20};
1067	<	double dxl[4] = {5,2.5,2.5,2.5};
1068	<	double dxh[4] = {5,2.5,2.5,10};
1066	>	double x[4] = {5,12.5,17.5,25};
1067	>	double dxl[4] = {5,2.5,2.5,5};
1068	>	double dxh[4] = {5,2.5,2.5,5};
1069
1070		string name="";
1071		string color="";
1072
1073		if(i==0) {
1074	<	double y[4] = {1.01103,1.01496,1.01486,1.02785};
1075	<	double dyl[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
1076	<	double dyh[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
1074	>	double y[4] = {1.01093,1.01379,1.0147,1.02766};
1075	>	double dyl[4] = {0.00243596,0.00242534,0.00341941,0.00600663};
1076	>	double dyh[4] = {0.00243596,0.00242534,0.00341941,0.00600663};
1077		name="30/30";
1078		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1079		color="#a01d99";
1080		}
1081		if(i==1) {
1082	<	double y[4] = {1.01103+0.001,1.01496+0.001,1.01486+0.001,1.02785+0.001};
1083	<	double dyl[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
1084	<	double dyh[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
1082	>	double y[4] = {1.0068,1.00779,0.999692,0.994484};
1083	>	double dyl[4] = {0.00249233,0.00257169,0.00371399,0.00631138};
1084	>	double dyh[4] = {0.00249233,0.00257169,0.00371399,0.00631138};
1085		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1086		name="30/15";
1087		color="#2464bc";
1088		}
1089		if(i==2) {
1090	<	double y[4] = {1.01575,1.02879,1.0325,1.05456};
1091	<	double dyl[4] = {0.00288281,0.00278536,0.00367451,0.00590421};
1092	<	double dyh[4] = {0.00288281,0.00278536,0.00367451,0.00590421};
1090	>	double y[4] = {1.01208,1.01991,1.02185,1.03479};
1091	>	double dyl[4] = {0.00366745,0.00400001,0.00584054,0.0104661};
1092	>	double dyh[4] = {0.00366745,0.00400001,0.00584054,0.0104661};
1093		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1094		name="15/15";
1095		color="#f49230";
1096		}
1097		if(i==3) {
1098	<	double y[4] = {1.01221,1.01922,1.02396,1.04677};
1099	<	double dyl[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
1100	<	double dyh[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
1098	>	double y[4] = {1.01177,1.01568,1.01616,1.03211};
1099	>	double dyl[4] = {0.00217317,0.00203014,0.00276838,0.00465643};
1100	>	double dyh[4] = {0.00217317,0.00203014,0.00276838,0.00465643};
1101		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1102	<	name="20/15";
1102	>	name="CHS 30/30";
1103		color="#c72f3c";
1104		}
1105		if(i==4) {
1106	<	double y[4] = {1.01221+0.001,1.01922+0.001,1.02396+0.001,1.04677+0.001};
1107	<	double dyl[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
1108	<	double dyh[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
1106	>	double y[4] = {1.00751,1.00964,1.00804,1.01073};
1107	>	double dyl[4] = {0.00224851,0.00215563,0.0029592,0.00499805};
1108	>	double dyh[4] = {0.00224851,0.00215563,0.0029592,0.00499805};
1109		gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1110	<	name="20/20";
1110	>	name="CHS 30/15";
1111		color="#7cb433";
1112		}
1113	+	if(i==5) {
1114	+	double y[4] = {1.01557,1.0318,1.03535,1.06331};
1115	+	double dyl[4] = {0.00300783,0.00302486,0.00421064,0.00739416};
1116	+	double dyh[4] = {0.00300783,0.00302486,0.00421064,0.00739416};
1117	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
1118	+	name="CHS 15/15";
1119	+	color="#DF013A";
1120	+	}
1121
1122	<
1122	>	gStyle->SetOptFit(0);
1123		gr[i]->SetTitle();
1124		gr[i]->GetXaxis()->SetTitle("N(Vtx)");
1125		gr[i]->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
#	Line 974 \| Line 1135 \| void TEST_ScenarioComparisonInclusive()
1135		fit[i] = (TF1*)gr[i]->GetFunction("pol1");
1136		fit[i]->SetLineColor(TColor::GetColor(color.c_str()));
1137		}
1138	<	leg->AddEntry(gr[i],name.c_str(),"l");
1138	>
1139	>	stringstream nameplus;
1140	>	nameplus << name << " (slope: " << std::setprecision(3) << fit[i]->GetParameter(1) << ")";
1141	>
1142	>	leg->AddEntry(gr[i],nameplus.str().c_str(),"l");
1143		if(i==0) gr[i]->Draw("AP*");
1144		else gr[i]->Draw("P");
1145		}
#	Line 986 \| Line 1151 \| void TEST_ScenarioComparisonInclusive()
1151
1152		}
1153
1154	<	void TEST_compare_selection(string identifier) {
1154	>	void compare_selection(string identifier) {
1155		bool recognized_scenario=false;
1156
1157		cout << "Running with identifier " << identifier << endl;
1158	+	if(identifier=="Zb30") {
1159	+	LeadingB=TCut ("Zb30_bTagProbCSVBP[0]>0.898");
1160	+	EtaB=TCut("abs(Zb30_pfJetGoodEta[0])<1.3");
1161	+	PhiZcut=TCut("abs(Zb30_pfJetDphiZ[0])>2.7");
1162	+	recognized_scenario=true;
1163	+	}
1164		if(identifier=="Zb3010") {
1165		LeadingB=TCut ("Zb3010_bTagProbCSVBP[0]>0.898");
1166		EtaB=TCut("abs(Zb3010_pfJetGoodEta[0])<1.3");
#	Line 997 \| Line 1168 \| void TEST_compare_selection(string ident
1168		recognized_scenario=true;
1169		}
1170
1171	<	if(identifier=="Zb3010") {
1172	<	LeadingB=TCut ("Zb3010_bTagProbCSVBP[0]>0.898");
1173	<	EtaB=TCut("abs(Zb3010_pfJetGoodEta[0])<1.3");
1174	<	PhiZcut=TCut("abs(Zb3010_pfJetDphiZ[0])>2.7");
1171	>	if(identifier=="ZbCHS3010") {
1172	>	LeadingB=TCut ("ZbCHS1510_bTagProbCSVBP[0]>0.898");
1173	>	EtaB=TCut("abs(ZbCHS1510_pfJetGoodEta[0])<1.3");
1174	>	PhiZcut=TCut("abs(ZbCHS1510_pfJetDphiZ[0])>2.7");
1175		recognized_scenario=true;
1176		}
1177
#	Line 1011 \| Line 1182 \| void TEST_compare_selection(string ident
1182		recognized_scenario=true;
1183		}
1184
1185	<	if(identifier=="Zb301030") {
1186	<	LeadingB=TCut ("Zb301030_bTagProbCSVBP[0]>0.898");
1187	<	EtaB=TCut("abs(Zb301030_pfJetGoodEta[0])<1.3");
1188	<	PhiZcut=TCut("abs(Zb301030_pfJetDphiZ[0])>2.7");
1185	>	if(identifier=="ZbCHS301030") {
1186	>	LeadingB=TCut ("ZbCHS301030_bTagProbCSVBP[0]>0.898");
1187	>	EtaB=TCut("abs(ZbCHS301030_pfJetGoodEta[0])<1.3");
1188	>	PhiZcut=TCut("abs(ZbCHS301030_pfJetDphiZ[0])>2.7");
1189		recognized_scenario=true;
1190		}
1191
#	Line 1025 \| Line 1196 \| void TEST_compare_selection(string ident
1196		recognized_scenario=true;
1197		}
1198
1199	<	if(identifier=="Zb301010") {
1200	<	LeadingB=TCut ("Zb301010_bTagProbCSVBP[0]>0.898");
1201	<	EtaB=TCut("abs(Zb301010_pfJetGoodEta[0])<1.3");
1202	<	PhiZcut=TCut("abs(Zb301010_pfJetDphiZ[0])>2.7");
1199	>	if(identifier=="ZbCHS301010") {
1200	>	LeadingB=TCut ("ZbCHS301010_bTagProbCSVBP[0]>0.898");
1201	>	EtaB=TCut("abs(ZbCHS301010_pfJetGoodEta[0])<1.3");
1202	>	PhiZcut=TCut("abs(ZbCHS301010_pfJetDphiZ[0])>2.7");
1203		recognized_scenario=true;
1204		}
1205
#	Line 1039 \| Line 1210 \| void TEST_compare_selection(string ident
1210		recognized_scenario=true;
1211		}
1212
1213	<	if(identifier=="Zb301010") {
1214	<	LeadingB=TCut ("Zb301010_bTagProbCSVBP[0]>0.898");
1215	<	EtaB=TCut("abs(Zb301010_pfJetGoodEta[0])<1.3");
1216	<	PhiZcut=TCut("abs(Zb301010_pfJetDphiZ[0])>2.7");
1213	>	if(identifier=="ZbCHS301010") {
1214	>	LeadingB=TCut ("ZbCHS301010_bTagProbCSVBP[0]>0.898");
1215	>	EtaB=TCut("abs(ZbCHS301010_pfJetGoodEta[0])<1.3");
1216	>	PhiZcut=TCut("abs(ZbCHS301010_pfJetDphiZ[0])>2.7");
1217		recognized_scenario=true;
1218		}
1219
#	Line 1053 \| Line 1224 \| void TEST_compare_selection(string ident
1224		recognized_scenario=true;
1225		}
1226
1227	<	if(identifier=="Zb3010_SecEta3") {
1228	<	LeadingB=TCut ("Zb3010_SecEta3_bTagProbCSVBP[0]>0.898");
1229	<	EtaB=TCut("abs(Zb3010_SecEta3_pfJetGoodEta[0])<1.3");
1230	<	PhiZcut=TCut("abs(Zb3010_SecEta3_pfJetDphiZ[0])>2.7");
1227	>	if(identifier=="ZbCHS3010_SecEta3") {
1228	>	LeadingB=TCut ("ZbCHS3010_SecEta3_bTagProbCSVBP[0]>0.898");
1229	>	EtaB=TCut("abs(ZbCHS3010_SecEta3_pfJetGoodEta[0])<1.3");
1230	>	PhiZcut=TCut("abs(ZbCHS3010_SecEta3_pfJetDphiZ[0])>2.7");
1231		recognized_scenario=true;
1232		}
1233
1234	<	if(identifier=="Zb3010_SecEta5") {
1235	<	LeadingB=TCut ("Zb3010_SecEta5_bTagProbCSVBP[0]>0.898");
1236	<	EtaB=TCut("abs(Zb3010_SecEta5_pfJetGoodEta[0])<1.3");
1237	<	PhiZcut=TCut("abs(Zb3010_SecEta5_pfJetDphiZ[0])>2.7");
1234	>	if(identifier=="ZbCHS3010_SecEta5") {
1235	>	LeadingB=TCut ("ZbCHS3010_SecEta5_bTagProbCSVBP[0]>0.898");
1236	>	EtaB=TCut("abs(ZbCHS3010_SecEta5_pfJetGoodEta[0])<1.3");
1237	>	PhiZcut=TCut("abs(ZbCHS3010_SecEta5_pfJetDphiZ[0])>2.7");
1238		recognized_scenario=true;
1239		}
1240
1241	<	if(identifier=="Zb3010_p5Clean") {
1242	<	LeadingB=TCut ("Zb3010_p5Clean_bTagProbCSVBP[0]>0.898 && Zb3010_p5Clean_IsClean");
1243	<	EtaB=TCut("abs(Zb3010_p5Clean_pfJetGoodEta[0])<1.3 && Zb3010_p5Clean_IsClean");
1244	<	PhiZcut=TCut("abs(Zb3010_p5Clean_pfJetDphiZ[0])>2.7 && Zb3010_p5Clean_IsClean");
1241	>	if(identifier=="ZbCHS3010_p5Clean") {
1242	>	LeadingB=TCut ("ZbCHS3010_p5Clean_bTagProbCSVBP[0]>0.898 && ZbCHS3010_p5Clean_IsClean");
1243	>	EtaB=TCut("abs(ZbCHS3010_p5Clean_pfJetGoodEta[0])<1.3 && ZbCHS3010_p5Clean_IsClean");
1244	>	PhiZcut=TCut("abs(ZbCHS3010_p5Clean_pfJetDphiZ[0])>2.7 && ZbCHS3010_p5Clean_IsClean");
1245		recognized_scenario=true;
1246		}
1247
1248	<	if(identifier=="Zb3010CHS") {
1249	<	LeadingB=TCut ("Zb3010CHS_bTagProbCSVBP[0]>0.898");
1250	<	EtaB=TCut("abs(Zb3010CHS_pfJetGoodEta[0])<1.3");
1251	<	PhiZcut=TCut("abs(Zb3010CHS_pfJetDphiZ[0])>2.7");
1248	>	if(identifier=="ZbCHS3010") {
1249	>	LeadingB=TCut ("ZbCHS3010_bTagProbCSVBP[0]>0.898");
1250	>	EtaB=TCut("abs(ZbCHS3010_pfJetGoodEta[0])<1.3");
1251	>	PhiZcut=TCut("abs(ZbCHS3010_pfJetDphiZ[0])>2.7");
1252	>	recognized_scenario=true;
1253	>	}
1254	>	if(identifier=="ZbCHS3010") {
1255	>	LeadingB=TCut ("ZbCHS3010_bTagProbCSVBP[0]>0.898");
1256	>	EtaB=TCut("abs(ZbCHS3010_pfJetGoodEta[0])<1.3");
1257	>	PhiZcut=TCut("abs(ZbCHS3010_pfJetDphiZ[0])>2.7");
1258	>	recognized_scenario=true;
1259	>	}
1260	>	if(identifier=="ZbCHS30") {
1261	>	LeadingB=TCut ("ZbCHS30_bTagProbCSVBP[0]>0.898");
1262	>	EtaB=TCut("abs(ZbCHS30_pfJetGoodEta[0])<1.3");
1263	>	PhiZcut=TCut("abs(ZbCHS30_pfJetDphiZ[0])>2.7");
1264		recognized_scenario=true;
1265		}
1266
#	Line 1089 \| Line 1272 \| void TEST_compare_selection(string ident
1272		cout << " " << (const char*) PhiZcut << endl;
1273
1274		cout << endl << endl;
1275	+
1276	+	// ScenarioComparison();
1277	+	// ScenarioComparisonInclusive();
1278	+
1279	+	// DoPUStudy(identifier);
1280	+	// SpecialCutFlow(identifier);
1281
1282	<	// ScenarioComparison();
1283	<	// ScenarioComparisonInclusive();
1284	<
1285	<
1097	<	// DoPUStudy(identifier);
1098	<
1099	<	// SpecialCutFlow(identifier);
1100	<	// draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",17000,30,200,"Z p_{T}","Official/"+identifier+"/Zpt__"+identifier,1);
1101	<	// draw_kin_variable(identifier+"_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",40,0,200,"Leading Jet Pt (B)","Official/"+identifier+"/LeadingJetPt__"+identifier,1);
1102	<	// draw_kin_variable(identifier+"_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",40,0,200,"Sub-Leading Jet Pt [GeV]","Official/"+identifier+"/Jet2Pt__"+identifier,1);
1103	<	// draw_kin_variable(identifier+"_pfJetGoodPt[1]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,2,"p_{T}^{J2} / p_{T}^{Z}","Official/"+identifier+"/SubLeadingJetPt_Over_ZPt__"+identifier,1);
1282	>	// draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",17000,30,200,"Z p_{T}","Official/"+identifier+"/Zpt__"+identifier,1);
1283	>	// draw_kin_variable(identifier+"_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",40,0,200,"Leading Jet Pt (B)","Official/"+identifier+"/LeadingJetPt__"+identifier,1);
1284	>	// draw_kin_variable(identifier+"_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",40,0,200,"Sub-Leading Jet Pt [GeV]","Official/"+identifier+"/Jet2Pt__"+identifier,1);
1285	>	// draw_kin_variable(identifier+"_pfJetGoodPt[1]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",20,0,2,"p_{T}^{J2} / p_{T}^{Z}","Official/"+identifier+"/SubLeadingJetPt_Over_ZPt__"+identifier,1);
1286		draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>30&&pt<50"),20,0,2,"#alpha","Official/"+identifier+"/alpha_pt_30_to_50__"+identifier,1);
1287		draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>50&&pt<75"),20,0,2,"#alpha","Official/"+identifier+"/alpha_50_to_75__"+identifier,1);
1288		draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>75&&pt<125"),20,0,2,"#alpha","Official/"+identifier+"/alpha_pt_75_to_125__"+identifier,1);
1289		draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>125&&pt<1000"),20,0,2,"#alpha","Official/"+identifier+"/alpha_pt_125_to_1000__"+identifier,1);
1290		draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut,20,0,2,"#alpha","Official/"+identifier+"/alpha_full__"+identifier,1);
1291	<	// draw_kin_variable(identifier+"_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","Official/"+identifier+"/DeltaPhiZBlead__"+identifier,0);
1291	>	draw_kin_variable(identifier+"_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>30&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","Official/"+identifier+"/DeltaPhiZBlead__"+identifier,0);
1292
1293	<	// draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.3").c_str()),20,0,2,"Z+b MPF","Official/"+identifier+"/mpf_alpha_smaller_0p3___"+identifier,0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
1294	<	// draw_kin_variable(identifier+"_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.3").c_str()),20,0,2,"p_{T} b-jet / p_{T} Z","Official/"+identifier+"/ptb_over_ptz___alpha_smaller_0p3______"+identifier,0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
1295	<
1114	<
1115	<	// */
1293	>	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.3").c_str()),20,0,2,"Z+b MPF","Official/"+identifier+"/mpf_alpha_smaller_0p3___"+identifier,0,"#alpha<0.3, 30 GeV < p_{T}^{Z} < 1000 GeV");
1294	>	draw_kin_variable(identifier+"_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.3").c_str()),20,0,2,"p_{T} b-jet / p_{T} Z","Official/"+identifier+"/ptb_over_ptz___alpha_smaller_0p3______"+identifier,0,"#alpha<0.3, 30 GeV < p_{T}^{Z} < 1000 GeV");
1295	>	draw_kin_variable(identifier+"_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<50&&"+identifier+"_alpha<0.3").c_str()),20,0,2,"p_{T} b-jet / p_{T} Z","Official/"+identifier+"/ptb_over_ptz___alpha_smaller_0p3___pt3050___"+identifier,0,"#alpha<0.3, 30 GeV < p_{T}^{Z} < 1000 GeV}{30<p_{T}^{Z}<50 GeV}");
1296
1297		}
1298
1299		void GetNumberEventsInsideOutsideAlphaWindow(TCut specialcut, string title) {
1300
1301	<	// TCut cut = ZplusBsel&&LeadingB&&PhiZcut&&specialcut&&TCut("Zb3010_alpha<0.3");
1301	>	// TCut cut = ZplusBsel&&LeadingB&&PhiZcut&&specialcut&&TCut("ZbCHS3010_alpha<0.3");
1302		TCut cut = specialcut;
1303		TCut in = EtaB;
1304	<	TCut out = TCut("abs(Zb3010_pfJetGoodEta[0])>1.3");
1304	>	TCut out = TCut("abs(ZbCHS3010_pfJetGoodEta[0])>1.3");
1305
1306		TH1F *data_IN = allsamples.Draw("data_IN", "mll",1,0,150, "nothing [GeV]", "events", cut&&in, data,luminosity);
1307		TH1F *data_OUT = allsamples.Draw("data_OUT", "mll",1,0,150, "nothing [GeV]", "events", cut&&out,data,luminosity);
#	Line 1167 \| Line 1347 \| void TEST_GetNumberEventsInsideOutsideAl
1347		GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2"), "ee/mm");
1348
1349		cout << "alpha cut" << endl;
1350	<	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==0")&&TCut("Zb3010_alpha<0.3"), "ee");
1351	<	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==1")&&TCut("Zb3010_alpha<0.3"), "mm");
1352	<	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2")&&TCut("Zb3010_alpha<0.3"), "ee/mm");
1350	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==0")&&TCut("ZbCHS3010_alpha<0.3"), "ee");
1351	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==1")&&TCut("ZbCHS3010_alpha<0.3"), "mm");
1352	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2")&&TCut("ZbCHS3010_alpha<0.3"), "ee/mm");
1353
1354		delete ca;
1355
#	Line 1180 \| Line 1360 \| void TEST_GetNumberEventsInsideOutsideAl
1360
1361
1362		void do_basic_ZB_analysis(bool DoCompleteAnalysis) {
1363	<
1364	<	//https://twiki.cern.ch/twiki/bin/viewauth/CMS/BtagPOG#2011_Data_and_MC
1185	<	//https://twiki.cern.ch/twiki/pub/CMS/BtagPOG/SFb-ttbar_payload.txt
1186	<
1363	>	STDWEIGHT=TCut(cutWeight);
1364	>	cutWeight=TCut(STDWEIGHT*TCut("ZbCHS3010_BTagWgtT"));
1365		cout << "The lepton requirement is " << (const char*) leptoncut << endl;
1366	<	bool doquick=true;
1366	>	bool doquick=false;
1367
1368	<	TCanvas *zbcanvas = new TCanvas("zbcanvas","zbcanvas");
1368	>	TCanvas *zbcanvas = new TCanvas("zbcanvas","zbcanvas");
1369
1370	<	// compare_selection("Zb3010_p5Clean");
1371	<	// compare_selection("Zb3010CHS");
1372	<	// compare_selection("Zb3010");
1373	<	// compare_selection("Zb301030");
1374	<	// compare_selection("Zb1530");
1375	<	// compare_selection("Zb301010");
1376	<	// compare_selection("Zb3010_SecEta3");
1377	<	// compare_selection("Zb3010_SecEta5");
1378	<	// compare_selection("Zb1510");
1379	<	// compare_selection("Zb301010");
1380	<	// compare_selection("ZbMikko");
1381	<	// compare_selection("Zb3010");
1382	<	// compare_selection("Zb40");
1370	>	// ScenarioComparison();
1371	>	// ScenarioComparisonInclusive();
1372	>
1373	>	// compare_selection("ZbCHS3010_p5Clean");
1374	>	// compare_selection("ZbCHS3010");
1375	>	// compare_selection("Zb30");
1376	>	// compare_selection("Zb3010");
1377	>	// compare_selection("Zb1510");
1378	>	// compare_selection("ZbCHS30");
1379	>	// compare_selection("ZbCHS3010");
1380	>	// compare_selection("ZbCHS3010");
1381	>	// compare_selection("ZbCHS301010");
1382	>	// compare_selection("ZbCHS3010_SecEta3");
1383	>	// compare_selection("ZbCHS3010_SecEta5");
1384	>	// compare_selection("Zb1510");
1385	>	// compare_selection("ZbCHS301010");
1386	>	// compare_selection("ZbMikko");
1387	>	// compare_selection("ZbCHS3010");
1388	>	// compare_selection("Zb40");
1389
1390		if(!doquick) {
1391	<	print_all_b_yields();
1392	<	draw_mpf_vars();
1393	<	// DrawEvilCutFlow();
1391	>	//draw_kin_variable("mll","",28,60.0,200.0,"m_{ll}","mll",1,"");//nice and quick test to see if the normalization is ok (i.e. all data is processed etc.)
1392	>	//print_all_b_yields();
1393	>	// draw_mpf_vars();
1394	>	//DrawEvilCutFlow();
1395
1396	<	draw_Zb_kin_vars();
1396	>	// draw_Zb_kin_vars();
1397
1398		}
1399
1400		// GetNumberEventsInsideOutsideAlphaWindow();
1401
1217	–	dout <<"THIS IS WHERE IT GETS REALLY INTERESTING!!!!" << endl;
1402		if(DoCompleteAnalysis) {
1403		/* need to run the following scenarios:
1404		* - normal
#	Line 1224 \| Line 1408 \| void do_basic_ZB_analysis(bool DoComplet
1408		* - alpha down
1409		*/
1410
1411	<	ZbResultContainer inclusive = new_data_mc_agreement_2d("Zb3010_alpha","inclusive",TCut("Zb3010_bTagProbCSVBP[0]>-10000"),TCut("1.0"));
1412	<	ZbResultContainer Reference = new_data_mc_agreement_2d("Zb3010_alpha","reference",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"));
1413	<	ZbResultContainer effmistagUp = new_data_mc_agreement_2d("Zb3010_alpha","effmistagUp",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTUp"));
1414	<	ZbResultContainer effmistagDown = new_data_mc_agreement_2d("Zb3010_alpha","effmistagDown",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTDown"));
1231	<	ZbResultContainer medium = new_data_mc_agreement_2d("Zb3010_alpha","medium",TCut("Zb3010_bTagProbCSVBP[0]>0.679"),TCut("ZbCHS3010_BTagWgtM"));
1232	<	ZbResultContainer loose = new_data_mc_agreement_2d("Zb3010_alpha","loose",TCut("Zb3010_bTagProbCSVBP[0]>0.244"),TCut("ZbCHS3010_BTagWgtL"));
1233	<	ZbResultContainer AlphaUp = new_data_mc_agreement_2d("Zb3010_alphaUp","AlphaUp",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"));
1234	<	ZbResultContainer AlphaDown = new_data_mc_agreement_2d("Zb3010_alphaDown","AlphaDown",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"));
1411	>	bool fast=true;bool slow=false;//don't change these
1412	>
1413	>	ZbResultContainer inclusive = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","inclusive",TCut("mll>0"),TCut("1.0"),0.3);
1414	>	ZbResultContainer Reference = new_data_mc_agreement_2d(slow,"ZbCHS3010_alpha","reference",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"),0.3);
1415
1416	+	ZbResultContainer NeutrinoQ = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","NeutrinoQ",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTUp(1.0softMuon)"),0.3);
1417	+	ZbResultContainer ANeutrinoQ = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","ANeutrinoQ",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTUp(1.0(1.0-softMuon))"),0.3);
1418	+
1419	+	ZbResultContainer effmistagUp = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","effmistagUp",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTUp"),0.3);
1420	+	ZbResultContainer effmistagDown = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","effmistagDown",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtTDown"),0.3);
1421	+	ZbResultContainer medium = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","medium",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.679"),TCut("ZbCHS3010_BTagWgtM"),0.3);
1422	+	ZbResultContainer loose = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","loose",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.244"),TCut("ZbCHS3010_BTagWgtL"),0.3);
1423	+	ZbResultContainer AlphaUp = new_data_mc_agreement_2d(fast,"ZbCHS3010_alphaUp","AlphaUp",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"),0.3);
1424	+	ZbResultContainer AlphaDown = new_data_mc_agreement_2d(fast,"ZbCHS3010_alphaDown","AlphaDown",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"),0.3);
1425	+
1426	+	ZbResultContainer L5inclusive = new_data_mc_agreement_2d(fast,"ZbCHS3010_alphaL5","L5inclusive",TCut("mll>0"),TCut("1.0"),0.3);
1427	+	ZbResultContainer L5Reference = new_data_mc_agreement_2d(fast,"ZbCHS3010_alphaL5","L5reference",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"),0.3);
1428	+
1429	+	ZbResultContainer AlphaThresholdVariationDown10 = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationDown10", TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"), TCut("ZbCHS3010_BTagWgtT"),0.3-0.1*0.3);
1430	+	ZbResultContainer AlphaThresholdVariationUp10 = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationUp10", TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"), TCut("ZbCHS3010_BTagWgtT"),0.3+0.1*0.3);
1431	+	ZbResultContainer AlphaThresholdVariationDown30 = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationDown30", TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"), TCut("ZbCHS3010_BTagWgtT"),0.3-0.3*0.3);
1432	+	ZbResultContainer AlphaThresholdVariationUp30 = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationUp30", TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"), TCut("ZbCHS3010_BTagWgtT"),0.3+0.3*0.3);
1433	+	ZbResultContainer AlphaThresholdVariationDown10inc = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationDown10inc",TCut("mll>0"), TCut("1.0"),0.3-0.1*0.3);
1434	+	ZbResultContainer AlphaThresholdVariationUp10inc = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationUp10inc", TCut("mll>0"), TCut("1.0"),0.3+0.1*0.3);
1435	+	ZbResultContainer AlphaThresholdVariationDown30inc = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationDown30inc",TCut("mll>0"), TCut("1.0"),0.3-0.3*0.3);
1436	+	ZbResultContainer AlphaThresholdVariationUp30inc = new_data_mc_agreement_2d(fast,"ZbCHS3010_alpha","AlphaThresholdVariationUp30inc", TCut("mll>0"), TCut("1.0"),0.3+0.3*0.3);
1437	+
1438		//and now let's compute some systematics!
1439		//1 alpha variation
1440		float SysMPF_Alpha_down = Reference.MPF_Result_FaceValue.getValue()-AlphaDown.MPF_Result_FaceValue.getValue();
#	Line 1242 \| Line 1444 \| void do_basic_ZB_analysis(bool DoComplet
1444		float SysMPF_EFF_down = Reference.MPF_Result_FaceValue.getValue()-effmistagDown.MPF_Result_FaceValue.getValue();
1445		float SysMPF_EFF_up = Reference.MPF_Result_FaceValue.getValue()-effmistagUp.MPF_Result_FaceValue.getValue();
1446
1447	<	//3 purity
1447	>	//3 Presence of soft muons (->neutrino question)
1448	>	float SysMPF_Neutrino_up = Reference.MPF_Result_FaceValue.getValue()-NeutrinoQ.MPF_Result_FaceValue.getValue();
1449	>	float SysMPF_Neutrino_down = Reference.MPF_Result_FaceValue.getValue()-ANeutrinoQ.MPF_Result_FaceValue.getValue();
1450	>
1451	>	//4 purity
1452		zbcanvas->cd();
1453	+	gStyle->SetOptFit(0);
1454		write_info(__FUNCTION__,"Don't know the exact purity at the moment, still needs to be determined!");
1455	+	TH1F *purityhisto = new TH1F("purityhisto","purityhisto",1,0,1);
1456	+	purityhisto->GetXaxis()->SetTitle("Purity");
1457	+	purityhisto->GetYaxis()->SetTitle("C_{abs}");
1458	+	purityhisto->GetXaxis()->CenterTitle();
1459	+	purityhisto->GetYaxis()->CenterTitle();
1460	+
1461		TGraphErrors *gSysMPF_Purity = new TGraphErrors();
1462	<	gSysMPF_Purity->SetPoint(0,0.05,inclusive.MPF_Result_FaceValue.getValue());
1462	>
1463	>	gSysMPF_Purity->SetPoint(0,PurityInclusive,inclusive.MPF_Result_FaceValue.getValue());
1464		gSysMPF_Purity->SetPointError(0,0.0,inclusive.MPF_Result_FaceValue.getError());
1465	<	gSysMPF_Purity->SetPointError(1,0.20,loose.MPF_Result_FaceValue.getError());//x value?
1465	>
1466	>	gSysMPF_Purity->SetPoint(1,PurityLoose,loose.MPF_Result_FaceValue.getValue());//x value?
1467		gSysMPF_Purity->SetPointError(1,0.0,loose.MPF_Result_FaceValue.getError());
1468	<	gSysMPF_Purity->SetPointError(2,0.54,medium.MPF_Result_FaceValue.getError());//x value?
1468	>	gSysMPF_Purity->SetPoint(2,PurityMedium,medium.MPF_Result_FaceValue.getValue());//x value?
1469		gSysMPF_Purity->SetPointError(2,0.0,medium.MPF_Result_FaceValue.getError());
1470	<	gSysMPF_Purity->SetPointError(3,0.82,Reference.MPF_Result_FaceValue.getError());//x value?
1470	>	gSysMPF_Purity->SetPoint(3,PurityTight,Reference.MPF_Result_FaceValue.getValue());//x value?
1471		gSysMPF_Purity->SetPointError(3,0.0,Reference.MPF_Result_FaceValue.getError());
1472
1473	<	gSysMPF_Purity->GetXaxis()->SetTitle("Purity");
1474	<	gSysMPF_Purity->GetXaxis()->CenterTitle();
1475	<	gSysMPF_Purity->GetYaxis()->SetTitle("C_{abs}");
1476	<	gSysMPF_Purity->GetYaxis()->CenterTitle();
1473	>
1474	>	float min,max;
1475	>	FindMinMax(gSysMPF_Purity,min,max);
1476	>	purityhisto->SetStats(0);
1477	>	purityhisto->GetYaxis()->SetRangeUser(min,max);
1478
1479		gSysMPF_Purity->Fit("pol1");
1480	<	gSysMPF_Purity->Draw();
1480	>	TPolyLine *purity_fit_uncert = GetFitUncertaintyShape(gSysMPF_Purity, "pol1", min, max,0.0,1.0);
1481	>	TF1* pufit = (TF1*)gSysMPF_Purity->GetFunction("pol1");
1482	>
1483	>	TF1* pufit_c = new TF1("pufit_c","[0]+[1]*x");
1484	>	pufit_c->SetParameters(pufit->GetParameters());
1485	>	pufit_c->SetLineColor(TColor::GetColor("#0101DF"));
1486	>	purity_fit_uncert->SetFillColor(TColor::GetColor("#5353FC"));
1487	>	pufit->SetLineColor(TColor::GetColor("#5353FC"));
1488	>
1489	>	purityhisto->Draw();
1490	>	purity_fit_uncert->Draw("f");
1491	>	gSysMPF_Purity->Draw("P*");
1492	>	pufit_c->Draw("same");
1493	>	DrawPrelim();
1494	>
1495	>	float ExtrapolatedResult=pufit->GetParameter(0)+1.0*pufit->GetParameter(1);
1496	>	float dExtrapolatedResult=sqrt(pufit->GetParError(0)pufit->GetParError(0)+1.01.0pufit->GetParError(1)pufit->GetParError(1));
1497	>
1498	>	stringstream summary;
1499	>	summary << "#splitline{Reference: " << std::setprecision(4) << Reference.MPF_Result_FaceValue.getValue() << "+/-" << std::setprecision(4) << Reference.MPF_Result_FaceValue.getError() << "}{";
1500	>	summary << "Extrapolation: " << std::setprecision(4) << ExtrapolatedResult << " +/- " << std::setprecision(4) << dExtrapolatedResult << "}";
1501	>
1502	>	dout << "Have found an extrapolated result at 1.0 of " << ExtrapolatedResult << " ; will use the difference between this result and the face value as ";
1503	>	dout << "purity-related systematic uncertainty" << endl;
1504	>
1505	>	TText *infobox = write_title(summary.str());
1506	>	infobox->SetX(0.75);
1507	>	infobox->SetTextSize(0.03);
1508	>	infobox->SetY(0.75);
1509	>	infobox->Draw();
1510	>
1511		CompleteSave(zbcanvas,"Systematics/Purity");
1512
1267	–	TF1 PurityFit = (TF1)gSysMPF_Purity->GetFunction("pol1");
1268	–	float ExtrapolatedResult=PurityFit->GetParameter(0)+1.0*PurityFit->GetParameter(1);
1513
1514	+
1515		float SysMPF_Purity = abs(ExtrapolatedResult-Reference.MPF_Result_FaceValue.getValue());
1516
1517	<	float sys_alpha = abs(SysMPF_Alpha_down)>abs(SysMPF_Alpha_up)?abs(SysMPF_Alpha_down):abs(SysMPF_Alpha_up);
1518	<	float sys_eff = abs(SysMPF_EFF_down)>abs(SysMPF_EFF_up)?abs(SysMPF_EFF_down):abs(SysMPF_EFF_up);
1519	<	float sys_purity = abs(SysMPF_Purity);
1520	<	float sys = sqrt(sys_alphasys_alpha+sys_effsys_eff+sys_purity*sys_purity);
1517	>	float sys_alpha = abs(SysMPF_Alpha_down)>abs(SysMPF_Alpha_up)?abs(SysMPF_Alpha_down):abs(SysMPF_Alpha_up);
1518	>	float sys_eff = abs(SysMPF_EFF_down)>abs(SysMPF_EFF_up)?abs(SysMPF_EFF_down):abs(SysMPF_EFF_up);
1519	>	float sys_purity = abs(SysMPF_Purity);
1520	>	float sys_neutrino = abs(SysMPF_Neutrino_down)>abs(SysMPF_Neutrino_up)?abs(SysMPF_Neutrino_down):abs(SysMPF_Neutrino_up);
1521	>	float sys = sqrt(sys_alphasys_alpha+sys_effsys_eff+sys_puritysys_purity+sys_neutrinosys_neutrino);
1522
1523		dout << "MPF METHOD:" << endl;
1524	<	dout << " Systematics : " << endl;
1525	<	dout << " Alpha variation : " << SysMPF_Alpha_down << " , " << SysMPF_Alpha_up << " ( " << sys_alpha << ") " << endl;
1524	>	dout << " Systematics : down up (selected)" << endl;
1525	>	dout << " Oversmearing : " << SysMPF_Alpha_down << " , " << SysMPF_Alpha_up << " ( " << sys_alpha << ") " << endl;
1526		dout << " eff/mistag variation : " << SysMPF_EFF_down << " , " << SysMPF_EFF_up << " ( " << sys_eff << ") " << endl;
1527		dout << " purity : " << SysMPF_Purity << " ( " << sys_purity << " )" << endl;
1528	<	dout << " total : " << sys << endl;
1528	>	dout << " neutrino : " << SysMPF_Neutrino_down << " , " << SysMPF_Neutrino_up << " ( " << sys_neutrino << ") " << endl;
1529	>	dout << " total : " << sys << " (note that this is the sys on C_{abs}^{b}, not yet on C_{corr} - there will also be an additional source!)" << endl;
1530		dout << endl << endl;
1531	<	dout << " FINAL RESULT : " << endl;
1532	<	dout << " " << Reference.MPF_Result_FaceValue/inclusive.MPF_Result_FaceValue << " (stat) +/- " << sys << " (sys) " << endl;
1531	>	dout << " FINAL RESULTS : " << endl;
1532	>	dout << " C_{abs}^{b} = " << Reference.MPF_Result_FaceValue << " (stat) +/- " << sys << " (sys) " << endl;
1533	>	dout << " C_{abs}^{inc} = " << inclusive.MPF_Result_FaceValue << " (stat) " << endl;
1534	>
1535	>	float sysfinal = Reference.MPF_Result_FaceValue.getValue()/inclusive.MPF_Result_FaceValue.getValue();
1536	>	sysfinal = sqrt(syssys/(Reference.MPF_Result_FaceValue.getValue()*Reference.MPF_Result_FaceValue.getValue()));
1537	>	//yes, this could be simplified but it would look like a bug.
1538
1539	<	LeadingB=TCut("Zb3010_bTagProbCSVBP[0]>0.898");
1539	>
1540	>
1541	>
1542	>	dout << "Note: Varying the alpha cut (->face value) we get the following results: " << endl;
1543	>	dout << "-30%: " << AlphaThresholdVariationDown30.MPF_Result_FaceValue/AlphaThresholdVariationDown30inc.MPF_Result_FaceValue << " (stat)" << endl;
1544	>	dout << "-10%: " << AlphaThresholdVariationDown10.MPF_Result_FaceValue/AlphaThresholdVariationDown10inc.MPF_Result_FaceValue << " (stat)" << endl;
1545	>	dout << " 0%: " << Reference.MPF_Result_FaceValue/inclusive.MPF_Result_FaceValue << " (stat) +/- " << endl;
1546	>	dout << "+10%: " << AlphaThresholdVariationUp10.MPF_Result_FaceValue/AlphaThresholdVariationUp10inc.MPF_Result_FaceValue << " (stat)" << endl;
1547	>	dout << "+30%: " << AlphaThresholdVariationUp30.MPF_Result_FaceValue/AlphaThresholdVariationUp30inc.MPF_Result_FaceValue << " (stat)" << endl;
1548	>
1549	>	dout << " XCheck for Rabs, varying the alpha cut (->face value) we get the following results: " << endl;
1550	>	dout << " -30%: " << AlphaThresholdVariationDown30.Rabs_Result_FaceValue/AlphaThresholdVariationDown30inc.Rabs_Result_FaceValue << " (stat)" << endl;
1551	>	dout << " -10%: " << AlphaThresholdVariationDown10.Rabs_Result_FaceValue/AlphaThresholdVariationDown10inc.Rabs_Result_FaceValue << " (stat)" << endl;
1552	>	dout << " 0%: " << Reference.Rabs_Result_FaceValue/inclusive.Rabs_Result_FaceValue << " (stat) +/- " << endl;
1553	>	dout << " +10%: " << AlphaThresholdVariationUp10.Rabs_Result_FaceValue/AlphaThresholdVariationUp10inc.Rabs_Result_FaceValue << " (stat)" << endl;
1554	>	dout << " +30%: " << AlphaThresholdVariationUp30.Rabs_Result_FaceValue/AlphaThresholdVariationUp30inc.Rabs_Result_FaceValue << " (stat)" << endl;
1555	>
1556	>
1557	>	float sys_alphavar_up = abs(AlphaThresholdVariationUp30.MPF_Result_FaceValue.getValue()/AlphaThresholdVariationUp30inc.MPF_Result_FaceValue.getValue() - Reference.MPF_Result_FaceValue.getValue()/inclusive.MPF_Result_FaceValue.getValue());
1558	>	float sys_alphavar_down = abs(AlphaThresholdVariationDown30.MPF_Result_FaceValue.getValue()/AlphaThresholdVariationDown30inc.MPF_Result_FaceValue.getValue() - Reference.MPF_Result_FaceValue.getValue()/inclusive.MPF_Result_FaceValue.getValue());
1559	>	float sys_alphavar = sys_alphavar_up>sys_alphavar_down?sys_alphavar_up:sys_alphavar_down;
1560	>
1561	>	dout << " -> Additional systematic (from 30%) : " << sys_alphavar << endl;
1562	>	sysfinal=sqrt(sysfinalsysfinal+sys_alphavarsys_alphavar);
1563	>	dout << "Final systematic is : " << sysfinal << endl;
1564	>
1565	>	dout << " ****************************************** < FINAL RESULT > ****************************************** " << endl;
1566	>	dout << " C_{corr} = " << Reference.MPF_Result_FaceValue/inclusive.MPF_Result_FaceValue << " (stat) +/- " << sysfinal << " (sys) " << endl;
1567	>	dout << " ****************************************** < /FINAL RESULT > ****************************************** " << endl;
1568	>
1569	>	dout << "Note that if L5 corrections are applied, we get " << L5Reference.MPF_Result_FaceValue/L5inclusive.MPF_Result_FaceValue << " (stat) " << endl;
1570	>
1571	>	LeadingB=TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898");
1572
1573		} else {
1574	<	ZbResultContainer Reference = new_data_mc_agreement_2d("Zb3010_alpha","reference",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgt"));
1574	>	ZbResultContainer Reference = new_data_mc_agreement_2d(false,"ZbCHS3010_alpha","reference",TCut("ZbCHS3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtT"),0.3);
1575		}
1576
1577
#	Line 1315 \| Line 1599 \| void SpecialCutFlow(string identifier) {
1599		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB) << ")";
1600		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB) << ")";
1601		MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut) << ")";
1602	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000").c_str())) << ")";
1603	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.3").c_str())) << ")";
1604	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.2").c_str())) << ")";
1605	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.15").c_str())) << ")";
1606	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.1").c_str())) << ")";
1607	<	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.05").c_str())) << ")";
1602	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000").c_str())) << ")";
1603	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.3").c_str())) << ")";
1604	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.2").c_str())) << ")";
1605	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.15").c_str())) << ")";
1606	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.1").c_str())) << ")";
1607	>	MegaCut << " + 1(" << (const char) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>30&&pt<1000&&"+identifier+"_alpha<0.05").c_str())) << ")";
1608		MegaCut << " ";
1609
1610		TCanvas *can = new TCanvas("can","can");
#	Line 1333 \| Line 1617 \| void SpecialCutFlow(string identifier) {
1617
1618
1619		TH1F *data = allsamples.Draw("data", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,data,luminosity);
1620	<	THStack mcm = allsamples.DrawStack("mc", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,mc,luminosity);
1620	>	THStack mcm = allsamples.DrawStack("mc", MegaCut.str(),11,-0.5,10.5, "", "events", basecut,mc,luminosity);
1621
1622		float runningsum=0;
1623		for(int i=data->GetNbinsX();i>0;i--) {
#	Line 1371 \| Line 1655 \| void SpecialCutFlow(string identifier) {
1655		data->GetXaxis()->SetBinLabel(3,"Leading jet is b-jet");
1656		data->GetXaxis()->SetBinLabel(4,"\|#eta_{b}\|<1.3 ");
1657		data->GetXaxis()->SetBinLabel(5,"\|#delta#phi(b,Z)\|>2.7");
1658	<	data->GetXaxis()->SetBinLabel(6,"10<p_{T}^{Z}<1000 GeV");
1658	>	data->GetXaxis()->SetBinLabel(6,"30<p_{T}^{Z}<1000 GeV");
1659		data->GetXaxis()->SetBinLabel(7,"#alpha < 0.3");
1660		data->GetXaxis()->SetBinLabel(8,"#alpha < 0.2");
1661		data->GetXaxis()->SetBinLabel(9,"#alpha < 0.15");
#	Line 1389 \| Line 1673 \| void SpecialCutFlow(string identifier) {
1673		leg->AddEntry(crap,neventsinfo.str().c_str(),"l");
1674
1675		data->Draw("e1");
1676	<	mcm.Draw("same");
1676	>	mcm.Draw("histo,same");
1677		data->Draw("same,e1,axis");
1678		data->Draw("same,e1");
1679		// data->Draw("same,TEXT");
#	Line 1402 \| Line 1686 \| void SpecialCutFlow(string identifier) {
1686
1687		delete crap;
1688		delete data;
1689	<
1689	>	CleanLegends();
1690	>	DeleteStack(mcm);
1691	>
1692
1693		}

Diff Legend

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

Comparing UserCode/cbrown/Development/Plotting/Modules/ZbTools.C (file contents): Revision 1.16 by buchmann, Wed Nov 28 17:51:14 2012 UTC vs. Revision 1.22 by buchmann, Mon Jan 21 14:43:36 2013 UTC

Diff Legend

Comparing UserCode/cbrown/Development/Plotting/Modules/ZbTools.C (file contents):
Revision 1.16 by buchmann, Wed Nov 28 17:51:14 2012 UTC vs.
Revision 1.22 by buchmann, Mon Jan 21 14:43:36 2013 UTC