ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/cbrown/Development/Plotting/Modules/ZbTools.C (file contents):
Revision 1.7 by buchmann, Wed Nov 7 10:30:51 2012 UTC vs.
Revision 1.15 by buchmann, Tue Nov 27 07:25:55 2012 UTC

#	Line 23 | Line 23 | using namespace PlottingSetup;
23		namespace ZbData {
24		vector<float> data_over_mc;
25
26	<	TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120");
27	<	TCut LeadingB("Zb30_bTagProbCSVBP[0]>0.679");
28	<	TCut EtaB("abs(Zb30_pfJetGoodEta[0])<1.3");
29	<	TCut PhiZcut("abs(pfJetDphiZ[0])>2.7");
26	>	TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120&&id1==id2");
27	>	TCut LeadingB("Zb3010_bTagProbCSVBP[0]>0.898");
28	>	TCut EtaB("abs(Zb3010_pfJetGoodEta[0])<1.3");
29	>	TCut PhiZcut("abs(Zb3010_pfJetDphiZ[0])>2.7");
30
31		}
32
33	+	class ZbPointResult {
34	+	public:
35	+	ZbPointResult(string var, float ptlow, float pthigh, float alphalow, float alphahigh, Value Data, Value MC);
36	+	float PtLow;
37	+	float PtHigh;
38	+	float AlphaLow;
39	+	float AlphaHigh;
40	+	string variable;
41	+	Value Data;
42	+	Value MC;
43	+	Value Ratio;
44	+	};
45	+
46	+	ZbPointResult::ZbPointResult(string var, float ptlow, float pthigh, float alphalow, float alphahigh, Value data, Value mc) :
47	+	variable(var), PtLow(ptlow), PtHigh(pthigh), AlphaLow(alphalow), AlphaHigh(alphahigh), Data(data), MC(mc) {
48	+	Ratio = data/mc;
49	+	}
50	+
51	+	std::ostream &operator<<(std::ostream &ostr, ZbPointResult &b)
52	+	{
53	+	ostr << b.PtLow << ";" << b.PtHigh << ";" << b.AlphaLow << ";" << b.AlphaHigh << ";" << b.Data.getValue() << ";" << b.Data.getError() << ";" << b.MC.getValue() << ";" << b.MC.getError() << ";" << b.Ratio.getValue() << ";" << b.Ratio.getError();
54	+	return ostr;
55	+	}
56	+
57	+
58	+	class ZbResultContainer {
59	+	public:
60	+	string name;
61	+	vector<ZbPointResult> MPF_Rabs_location;
62	+	Value MPF_Result_FaceValue;
63	+	Value Rabs_Result_FaceValue;
64	+	Value MPF_Result_Extrapolated;
65	+	Value Rabs_Result_Extrapolated;
66	+
67	+	ZbResultContainer(string name);
68	+	ZbResultContainer(const ZbResultContainer &old);
69	+	};
70	+
71	+
72	+	ZbResultContainer::ZbResultContainer(string _name) {
73	+	name=_name;
74	+	cout<< "Result Container for \"" << name << "\" has been initialized" << endl;
75	+
76	+	}
77	+
78	+	std::ostream &operator<<(std::ostream &ostr, ZbResultContainer &b)
79	+	{
80	+	ostr << "Results for " << b.name << endl;
81	+	ostr << " MPF and Rabs locations: " << endl;
82	+	for(int i=0;i<b.MPF_Rabs_location.size();i++) ostr << "   " << b.MPF_Rabs_location[i] << endl;
83	+	ostr << " Results: " << endl;
84	+	ostr << "   Face Value: " << endl;
85	+	ostr << "     MPF " << b.MPF_Result_FaceValue << endl;
86	+	ostr << "     Rabs " << b.Rabs_Result_FaceValue << endl;
87	+	ostr << "   Extrapolated: " << endl;
88	+	ostr << "     MPF " << b.MPF_Result_Extrapolated << endl;
89	+	ostr << "             Rabs " << b.Rabs_Result_Extrapolated << endl;
90	+
91	+	return ostr;
92	+	}
93	+
94	+
95	+
96		using namespace ZbData;
97
98
#	Line 43 | Line 106 | void print_yield(TCut cut) {
106		TH1F *data  = allsamples.Draw("data",    "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,data,luminosity);
107		dout << "  Yields for " << (const char*) cut << endl;
108		dout << "     data: " << data->Integral() << endl;
109	<	/*  THStack mcm = allsamples.DrawStack("mc", "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,mc,luminosity);
109	>	THStack mcm = allsamples.DrawStack("mc", "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,mc,luminosity);
110		int nhists = mcm.GetHists()->GetEntries();
111		dout << "Going to loop over " << nhists << " histograms " << endl;
112		TFile *test = new TFile("test.root","RECREATE");
#	Line 53 | Line 116 | void print_yield(TCut cut) {
116		TH1* hist = static_cast<TH1*>(mcm.GetHists()->At(i));
117		cout << "     " << hist->GetName() << ": " << hist->Integral() << endl;
118		}
119	+
120		cout << endl << endl;
121	<	*/
121	>
122		delete data;
123		}
124
#	Line 100 | Line 164 | void draw_kin_variable(string variable,
164		delete ckin;
165		}
166
103	–	Value getfrom2Dmap(TH2F *map, int ixbin, int iybin) {
104	–	Value sum(0,0);
105	–	for(int iy=1;iy<=iybin;iy++) {
106	–	sum=sum+Value(map->GetBinContent(ixbin,iy),map->GetBinError(ixbin,iy));
107	–	}
108	–	return sum;
109	–	}
110	–
167		void print_all_b_yields() {
168		cout << "Basic selection with a b jet" << endl;
169	<	print_yield(ZplusBsel&&"Zb30_pfJetGoodNumBtag>0");
169	>	print_yield(ZplusBsel&&"Zb3010_pfJetGoodNumBtag>0");
170		cout << "Leading jet is a b-jet " << endl;
171		print_yield(ZplusBsel&&LeadingB);
172		cout << "|#eta_{b}|<1.3 " << endl;
#	Line 119 | Line 175 | void print_all_b_yields() {
175		print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
176		cout << "10<ptZ<1000 GeV" << endl;
177		print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000");
178	+	cout << "#alpha < 0.35" << endl;
179	+	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.35");
180		cout << "#alpha < 0.3" << endl;
181	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb30_pfJetGoodPt[1]/pt)<0.3");
181	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.3");
182		cout << "#alpha < 0.2" << endl;
183	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb30_pfJetGoodPt[1]/pt)<0.2");
183	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.2");
184		cout << "#alpha < 0.15" << endl;
185	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb30_pfJetGoodPt[1]/pt)<0.15");
185	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.15");
186		cout << "#alpha < 0.1" << endl;
187	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb30_pfJetGoodPt[1]/pt)<0.1");
188	<	cout << "#alpha < 0.05" << endl;
189	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb30_pfJetGoodPt[1]/pt)<0.05");
187	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.1");
188	>	//  cout << "#alpha < 0.05" << endl;
189	>	//  print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"(Zb3010_alpha)<0.05");
190		}
191
192	<	void DrawEvilCutFlow() {
192	>	void TEST_DrawEvilCutFlow() {
193		stringstream MegaCut;
194
195	<	MegaCut << "   1*(" << (const char*) (ZplusBsel&&TCut("Zb30_pfJetGoodNumBtag>0")) << ")";
195	>	MegaCut << "   1*(" << (const char*) (ZplusBsel&&TCut("Zb3010_pfJetGoodNumBtag>0")) << ")";
196		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB) << ")";
197		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB) << ")";
198		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut) << ")";
199		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000")) << ")";
200	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.3")) << ")";
201	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.2")) << ")";
202	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.15")) << ")";
203	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.1")) << ")";
204	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.05")) << ")";
200	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.3")) << ")";
201	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.2")) << ")";
202	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.15")) << ")";
203	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.1")) << ")";
204	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb3010_alpha<0.05")) << ")";
205		MegaCut << " ";
206
207		TCanvas *can = new TCanvas("can","can");
#	Line 210 | Line 268 | void DrawEvilCutFlow() {
268
269		void draw_Zb_kin_vars() {
270		draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Z p_{T}","Official/Zpt",1);
271	<	draw_kin_variable("Zb30_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Sub-Leading Jet Pt","Official/Jet2Pt",1);
272	<	draw_kin_variable("Zb30_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Leading Jet Pt (B)","Official/LeadingJetPt",1);
273	<	draw_kin_variable("Zb30_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);
274	<	draw_kin_variable("Zb30_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<50",   20,0,2,"#alpha","Official/alpha_pt_10_to_50",1);
275	<	draw_kin_variable("Zb30_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>50&&pt<100",  20,0,2,"#alpha","Official/alpha_50_to_100",1);
276	<	draw_kin_variable("Zb30_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>100&&pt<200", 20,0,2,"#alpha","Official/alpha_pt_100_to_200",1);
277	<	draw_kin_variable("Zb30_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>200&&pt<1000",20,0,2,"#alpha","Official/alpha_pt_200_to_1000",1);
271	>	draw_kin_variable("Zb3010_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Sub-Leading Jet Pt","Official/Jet2Pt",1);
272	>	draw_kin_variable("Zb3010_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Leading Jet Pt (B)","Official/LeadingJetPt",1);
273	>	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);
274	>	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<50",   40,0,2,"#alpha","Official/alpha_pt_30_to_50",1);
275	>	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>50&&pt<75",  40,0,2,"#alpha","Official/alpha_50_to_75",1);
276	>	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>75&&pt<125", 40,0,2,"#alpha","Official/alpha_pt_75_to_125",1);
277	>	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>125&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_125_to_1000",1);
278	>	draw_kin_variable("Zb3010_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>30&&pt<1000",40,0,2,"#alpha","Official/alpha_pt_30_to_1000",1);
279
280	<	draw_kin_variable("pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",0);
280	>	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);
281
282		/*
283		draw_kin_variable("pt1",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,200,"p_{T}^{l1}","Lep/pt1",1);
#	Line 239 | Line 298 | void draw_Zb_kin_vars() {
298
299		void draw_mpf_vars() {
300		ZbData::data_over_mc.clear();
301	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.3",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p3",0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
302	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.2",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p2",0,"#alpha<0.2, 10 GeV < p_{T}^{Z} < 1000 GeV");
303	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.15",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p15",0,"#alpha<0.15, 10 GeV < p_{T}^{Z} < 1000 GeV");
304	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.1",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p1",0,"#alpha<0.1, 10 GeV < p_{T}^{Z} < 1000 GeV");
305	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.05",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p05",0,"#alpha<0.05, 10 GeV < p_{T}^{Z} < 1000 GeV");
306	<
307	<	draw_kin_variable("Zb30_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_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");
308	<	draw_kin_variable("Zb30_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_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");
309	<	draw_kin_variable("Zb30_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_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");
310	<	draw_kin_variable("Zb30_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_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("Zb30_pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_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");
301	>	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");
302	>	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");
303	>	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");
304	>	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");
305	>	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");
306	>
307	>	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");
308	>	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");
309	>	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");
310	>	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");
312		}
313
314	<	Value get_Zb_data_over_mc(string variable, TCut cut) {
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	>
325	>	Value get_Zb_data_over_mc(ZbResultContainer &res, string variable, string variable2, TCut cut, string saveas, float pt1, float pt2, float a1, float a2) {
326		//  write_warning(__FUNCTION__,"Debugging this function, therefore always returning 3.1415+/-1.010101"); return Value(3.1415,1.010101);
327	<	TH1F *hdata  = allsamples.Draw("hdata",variable,1000,0,30, "#alpha [GeV]", "events", cut,data,luminosity);
328	<	TH1F *hmc    = allsamples.Draw("hmc" , variable,1000,0,30, "#alpha [GeV]", "events", cut,mc,  luminosity);
327	>	TCanvas *cn = new TCanvas("cn","cn");
328	>	string varname="MPF";
329	>	if(Contains(variable,"pfJetGood")) varname="R_{abs}";
330	>	TH1F *hdata  = allsamples.Draw("hdata",variable,1200,0,30, varname, "events",  cut,data,luminosity);
331	>	TH1F *hmc    = allsamples.Draw("hmc" , variable,1200,0,30, varname, "events", cut, mc,  luminosity);
332	>	hmc->SetLineColor(kBlue);
333
334		float a=hdata->GetMean();
335		float b=hmc->GetMean();
#	Line 263 | Line 337 | Value get_Zb_data_over_mc(string variabl
337		float db=hmc->GetMeanError();
338		float factor = a / b;
339		float error = TMath::Sqrt(  (1/(b*b)) * (da*da) + ((a*a)/(b*b))*db*db);
340	+
341	+	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);
346	+	TF1 *gaus_mc = new TF1("gaus_mc","gaus",0.7,1.3);
347	+	gaus_mc->SetParameter(0,hmc->GetMaximum());
348	+	gaus_mc->SetParameter(1,hmc->GetMean());
349	+	gaus_mc->SetParameter(2,hmc->GetRMS());
350	+	hmc->Fit(gaus_mc);
351	+
352	+
353	+	cn->cd();
354	+	hdata->GetYaxis()->SetTitle("events / 0.1");
355	+	hdata->SetMarkerColor(kRed);
356	+	/*  hdata->Rebin(4);
357	+	hmc->Rebin(4);
358	+
359	+	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);
364	+
365	+	hdata->GetXaxis()->SetRangeUser(0,2);
366	+	hdata->Draw("e1");
367	+	hmc->Draw("histo,same");
368	+	hdata->Draw("e1,same");
369	+	gaus_data->Draw("same");
370	+	gaus_mc->Draw("same");
371	+	stringstream summary;
372	+	summary << "#splitline{" << varname << "}{#splitline{data: " << std::setprecision(4) << a << " #pm " << std::setprecision(4) << da << "}{";
373	+	summary << "#splitline{mc:   " << std::setprecision(4) << b << " #pm " << std::setprecision(4) << db << "}{#splitline{ratio: " << factor;
374	+	summary << " #pm " << error << "}{#splitline{}{#splitline{data fit: " << std::setprecision(4) << gaus_data->GetParameter(1) << "+/-";
375	+	summary << std::setprecision(4) << gaus_data->GetParError(1) << "}{#splitline{";
376	+	summary << "mc fit:" << std::setprecision(4) << gaus_mc->GetParameter(1) << "+/-" << std::setprecision(4) << gaus_mc->GetParError(1) << "}{ratio: ";
377	+	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)) ) << "}}}}}}}";
378	+	TText *infobox = write_title(summary.str());
379	+	infobox->SetX(0.75);
380	+	infobox->SetTextSize(0.03);
381	+	infobox->SetY(0.75);
382	+	infobox->Draw();
383	+
384	+
385	+	DrawPrelim();
386	+	CompleteSave(cn,"ResponsePlots/"+saveas);
387	+
388	+	res.MPF_Rabs_location.push_back(ZbPointResult(variable,pt1,pt2,a1,a2,Value(a,da),Value(b,db)));
389	+	cout << "Data;" << a << ";"<<da<<";MC"<<b<<";"<<db<<endl;
390	+	delete cn;
391		delete hdata;
392		delete hmc;
393	+
394		return Value(factor,error);
395		}
396
397	+	ZbResultContainer new_data_mc_agreement_2d(string AlphaVariable, string ContainerName, TCut BTagCut, TCut BTagWgt) {
398	+	ZbResultContainer results(ContainerName);
399	+
400	+	cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)))")*BTagWgt;
401	+	LeadingB=BTagCut;
402
272	–
273	–	void new_data_mc_agreement_2d() {
403		gStyle->SetOptFit(0);
404		TCut basecut(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
405	<	const int nalphacuts=6;
406	<	float alphacuts[nalphacuts] = {0.0,0.1,0.15,0.2,0.25,0.3};
407	<	int PointThree=-1;
408	<	for(int i=0;i<nalphacuts;i++) {
409	<	if((alphacuts[i]-0.3)<0.01) PointThree=i;
410	<	}
411	<	if(PointThree<0) {
283	<	write_warning(__FUNCTION__,"alpha=0.3 MUST BE ONE OF THE ALPHA VALUES!");
284	<	assert(PointThree<0);
285	<	}
405	>	//  const int nalphacuts=6;
406	>	//  float alphacuts[nalphacuts] = {0.1,0.15,0.2,0.25,0.3,0.35};
407	>
408	>	//  const int nptcuts=5;
409	>	//  float ptcuts[nptcuts]={30,50,75,125,1000};
410	>
411	>	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};
412
287	–	const int nptcuts=5;
288	–	float ptcuts[nptcuts]={10,50,100,200,1000};
413
414		float MPF_Results[nalphacuts][nptcuts];
415		float MPF_Errors[nalphacuts][nptcuts];
#	Line 293 | Line 417 | void new_data_mc_agreement_2d() {
417		float RABS_Errors[nalphacuts][nptcuts];
418
419
420	<	for(int ia=0;ia<nalphacuts-1;ia++) {
420	>	for(int ia=0;ia<nalphacuts;ia++) {
421		for(int ipt=0;ipt<nptcuts-1;ipt++) {
422		stringstream specialcut;
423	<	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (Zb30_alpha<" << alphacuts[ia+1] << "&&Zb30_alpha>" << alphacuts[ia] << "))";
424	<	Value MPF_data_over_mc = get_Zb_data_over_mc("mpf",TCut(basecut && specialcut.str().c_str()));
425	<	Value RABS_data_over_mc = get_Zb_data_over_mc("Zb30_pfJetGoodPt[0]/pt",TCut(basecut && specialcut.str().c_str()));
423	>	float alow,ahigh;
424	>	if(ia>0) {
425	>	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (" << AlphaVariable << "<" << alphacuts[ia] << " && " << AlphaVariable << ">" << alphacuts[ia-1] << "))";
426	>	alow=alphacuts[ia-1];
427	>	ahigh=alphacuts[ia];
428	>	} else {
429	>	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (" << AlphaVariable << "<" << alphacuts[ia] << "))";
430	>	alow=0;
431	>	ahigh=alphacuts[ia];
432	>	}
433	>	cout << specialcut.str() << endl;
434	>	Value MPF_data_over_mc = get_Zb_data_over_mc(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);
435	>	Value RABS_data_over_mc = get_Zb_data_over_mc(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);
436
437		MPF_Results[ia][ipt]=MPF_data_over_mc.getValue();
438		MPF_Errors[ia][ipt]=MPF_data_over_mc.getError();
#	Line 306 | Line 440 | void new_data_mc_agreement_2d() {
440		RABS_Results[ia][ipt]=RABS_data_over_mc.getValue();
441		RABS_Errors[ia][ipt]=RABS_data_over_mc.getError();
442
443	<	cout << "alpha cut at " << alphacuts[ia] << " and pt range " << ptcuts[ipt] << " , " << ptcuts[ipt+1] << " \t : \t" << specialcut.str() << " \t \t --> " << MPF_data_over_mc << " (MPF) " << RABS_data_over_mc << " (RABS)" << endl;
443	>	cout << "alpha cut at " << alphacuts[ia+1] << " and pt range " << ptcuts[ipt] << " , " << ptcuts[ipt+1] << " \t : \t" << specialcut.str() << " \t \t --> " << MPF_data_over_mc << " (MPF) " << RABS_data_over_mc << " (RABS)" << endl;
444		}
445		}
446
#	Line 314 | Line 448 | void new_data_mc_agreement_2d() {
448		float RABS_ExtraPolatedResults[nptcuts];
449
450		TCanvas *can = new TCanvas("can");
451	+	can->cd(1)->SetLeftMargin(0.15);
452
453		TGraphErrors *MPF_FinalGraph = new TGraphErrors();
454		TGraphErrors *RABS_FinalGraph = new TGraphErrors();
455	<	TGraphErrors *FaceValueAtPoint3 = new TGraphErrors();
455	>	TGraphErrors *MPF_FaceValueAtPoint3 = new TGraphErrors();
456	>	TGraphErrors *RABS_FaceValueAtPoint3 = new TGraphErrors();
457
458		for(int ipt=0;ipt<nptcuts-1;ipt++) {
459
#	Line 326 | Line 462 | void new_data_mc_agreement_2d() {
462		cout << "Going to create histo for " << filename.str() << endl;
463		TGraphErrors *mpf_gr =new TGraphErrors();
464		TGraphErrors *rabs_gr =new TGraphErrors();
465	+
466	+	int rabs_counter=0;
467	+	int mpf_counter=0;
468	+
469		for(int ia=0;ia<nalphacuts;ia++) {
470	<	mpf_gr->SetPoint(ia,alphacuts[ia],MPF_Results[ia][ipt]);
471	<	mpf_gr->SetPointError(ia,0,MPF_Errors[ia][ipt]);
472	<	rabs_gr->SetPoint(ia,alphacuts[ia],RABS_Results[ia][ipt]);
473	<	rabs_gr->SetPointError(ia,0,RABS_Errors[ia][ipt]);
470	>	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;
471	>	if(MPF_Results[ia][ipt]==MPF_Results[ia][ipt] && MPF_Errors[ia][ipt]==MPF_Errors[ia][ipt]) { // remove nan's
472	>	mpf_gr->SetPoint(mpf_counter,alphacuts[ia],MPF_Results[ia][ipt]);
473	>	mpf_gr->SetPointError(mpf_counter,0,MPF_Errors[ia][ipt]);
474	>	mpf_counter++;
475	>	} else {
476	>	dout << "Detected BS for MPF method (bin ignored) : " << endl;
477	>	dout << "       alpha: " << alphacuts[ia] << " , pt range: " << ptcuts[ipt] << " < pt < " << ptcuts[ipt+1] << endl;
478	>	cout << "       value (MPF) : " << MPF_Results[ia][ipt] << " +/- " << MPF_Errors[ia][ipt] << endl;
479	>	}
480	>	if(RABS_Results[ia][ipt]==RABS_Results[ia][ipt] && RABS_Errors[ia][ipt]==RABS_Errors[ia][ipt]) {
481	>	rabs_gr->SetPoint(rabs_counter,alphacuts[ia],RABS_Results[ia][ipt]);
482	>	rabs_gr->SetPointError(rabs_counter,0,RABS_Errors[ia][ipt]);
483	>	rabs_counter++;
484	>	} else {
485	>	dout << "Detected BS for RABS method (bin ignored) : " << endl;
486	>	dout << "       alpha: " << alphacuts[ia] << " , pt range: " << ptcuts[ipt] << " < pt < " << ptcuts[ipt+1] << endl;
487	>	cout << "       value (RABS) : " << RABS_Results[ia][ipt] << " +/- " << RABS_Errors[ia][ipt] << endl;
488	>	}
489		}
335	–	FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1]),MPF_Results[PointThree][ipt]);
336	–	FaceValueAtPoint3->SetPointError(ipt,0,MPF_Errors[PointThree][ipt]);
490
491	+	stringstream specialcut;
492	+	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (" << AlphaVariable << "<0.3))";
493	+	Value MPF_data_over_mc = get_Zb_data_over_mc(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);
494	+	Value RABS_data_over_mc = get_Zb_data_over_mc(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);
495	+	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;
496	+	MPF_FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1])/2,MPF_data_over_mc.getValue());
497	+	MPF_FaceValueAtPoint3->SetPointError(ipt,0,MPF_data_over_mc.getError());
498	+	RABS_FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1])/2,RABS_data_over_mc.getValue());
499	+	RABS_FaceValueAtPoint3->SetPointError(ipt,0,RABS_data_over_mc.getError());
500	+
501	+	can->cd();
502		mpf_gr->SetMarkerStyle(21);
503		mpf_gr->SetMarkerColor(kBlue);
504		mpf_gr->Draw("AP");
505	<	mpf_gr->Fit("pol1","QE"); // quiet, use minos
505	>	mpf_gr->Fit("pol1","");
506		mpf_gr->GetXaxis()->SetTitle("#alpha");
507		mpf_gr->GetXaxis()->CenterTitle();
508	<	mpf_gr->GetYaxis()->SetTitle("data/mc");
508	>	mpf_gr->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
509		mpf_gr->GetYaxis()->CenterTitle();
510		TF1 *mpf_pol=(TF1*)mpf_gr->GetFunction("pol1");
511	+	mpf_pol->SetLineWidth(0);
512	+	TF1 *mpf_pol_complete = new TF1("mpf_pol_complete","[0]+[1]*x");
513	+	mpf_pol_complete->SetParameters(mpf_pol->GetParameters());
514	+	mpf_pol_complete->SetLineWidth(1);
515	+	mpf_pol_complete->SetLineColor(kBlue);
516	+	float min,max;
517	+	FindMinMax(mpf_gr,min,max);
518	+	TH1F *histo = new TH1F("histo","histo",1,0,0.4);
519	+	histo->GetXaxis()->SetTitle("#alpha");
520	+	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
521	+	histo->GetXaxis()->CenterTitle();
522	+	histo->GetYaxis()->CenterTitle();
523	+	histo->SetStats(0);
524	+	histo->GetXaxis()->SetRangeUser(0,0.4);
525	+	histo->GetYaxis()->SetRangeUser(min,max);
526	+	mpf_gr->GetYaxis()->SetRangeUser(min,max);
527	+
528	+	TPolyLine *mpf_fit_uncert = GetFitUncertaintyShape(mpf_gr, "pol1", min, max,0.0,0.4);
529	+	mpf_pol->SetLineColor(TColor::GetColor("#0101DF"));
530	+	mpf_fit_uncert->SetFillColor(TColor::GetColor("#5353FC"));
531	+	histo->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
532	+	histo->Draw();
533	+	mpf_fit_uncert->Draw("F");
534	+	mpf_fit_uncert->Draw("");
535	+	mpf_gr->Draw("P");
536	+	mpf_pol_complete->Draw("same");
537	+
538		float mpf_a=mpf_pol->GetParameter(0);
539		float mpf_b=mpf_pol->GetParameter(1);
540		MPF_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),mpf_a);
#	Line 352 | Line 543 | void new_data_mc_agreement_2d() {
543		MPF_ExtraPolatedResults[ipt]=mpf_a;
544
545		stringstream mpf_info;
546	<	mpf_info <<  "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (MPF) }}{Extrapolated value: " << std::setprecision(3) << mpf_a << "}}{#chi^{2}/NDF : " << mpf_pol->GetChisquare() << " / " << mpf_pol->GetNDF() << "}";
546	>	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() << "}";
547
548		TText *mpf_mark = write_title(mpf_info.str());
549		mpf_mark->SetX(0.75);
#	Line 363 | Line 554 | void new_data_mc_agreement_2d() {
554		string filenamebkp=filename.str();
555		filename << "__MPF";
556		DrawPrelim();
557	<	CompleteSave(can,filename.str());
557	>	CompleteSave(can,ContainerName+"/"+filename.str());
558		cout << "MPF : " << mpf_a << " + " << mpf_b << " * alpha " << endl;
559
560		filename.str("");
561
562		rabs_gr->SetMarkerStyle(21);
563	<	rabs_gr->SetMarkerColor(kBlue);
373	<	rabs_gr->Fit("pol1","QE"); // quiet, use minos
563	>	rabs_gr->SetMarkerColor(kRed);
564		rabs_gr->Draw("AP");
565	+	rabs_gr->Fit("pol1","");
566		rabs_gr->GetXaxis()->SetTitle("#alpha");
567		rabs_gr->GetXaxis()->CenterTitle();
568	<	rabs_gr->GetYaxis()->SetTitle("data/mc");
568	>	rabs_gr->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
569		rabs_gr->GetYaxis()->CenterTitle();
379	–
570		TF1 *rabs_pol=(TF1*)rabs_gr->GetFunction("pol1");
571	+	rabs_pol->SetLineWidth(0);
572	+	TF1 *rabs_pol_complete = new TF1("rabs_pol_complete","[0]+[1]*x");
573	+	rabs_pol_complete->SetParameters(rabs_pol->GetParameters());
574	+	rabs_pol_complete->SetLineColor(kRed);
575	+	rabs_pol_complete->SetLineWidth(1);
576	+	FindMinMax(rabs_gr,min,max);
577	+	rabs_gr->GetYaxis()->SetRangeUser(min,max);
578	+	histo->GetYaxis()->SetRangeUser(min,max);
579	+	TPolyLine *rabs_fit_uncert = GetFitUncertaintyShape(rabs_gr, "pol1", min, max,0.0,0.4);
580	+	rabs_pol->SetLineColor(TColor::GetColor("#FA5858"));
581	+	rabs_pol->SetFillColor(TColor::GetColor("#FA5858"));
582	+	rabs_fit_uncert->SetLineColor(TColor::GetColor("#FA5858"));
583	+	rabs_fit_uncert->SetFillColor(TColor::GetColor("#FA5858"));
584	+	histo->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
585	+	histo->Draw();
586	+	rabs_fit_uncert->Draw("F");
587	+	rabs_fit_uncert->Draw("");
588	+	rabs_gr->Draw("P");
589	+	rabs_pol_complete->Draw("same");
590	+
591	+
592		float rabs_a=rabs_pol->GetParameter(0);
593		float rabs_b=rabs_pol->GetParameter(1);
594		RABS_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),rabs_a);
#	Line 385 | Line 596 | void new_data_mc_agreement_2d() {
596
597		cout << "!+!+!+!+!+!+!+!+!+ Just added a point to the final plot : " << rabs_a << " +/- " << rabs_pol->GetParError(0) << endl;
598
388	–
599		stringstream rabs_info;
600	<	rabs_info << "#splitline{#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (R_{abs}) }}{Extrapolated value: " << std::setprecision(3) << rabs_a << "}}{#chi^{2}/NDF : " << rabs_pol->GetChisquare() << " / " << rabs_pol->GetNDF() << "}";
600	>	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() << "}";
601		TText *rabs_mark = write_title(rabs_info.str());
602		rabs_mark->SetX(0.75);
603		rabs_mark->SetTextSize(0.03);
#	Line 395 | Line 605 | void new_data_mc_agreement_2d() {
605		rabs_mark->Draw();
606
607
398	–
608		RABS_ExtraPolatedResults[ipt]=rabs_a;
609		filename << filenamebkp << "__RABS";
610		DrawPrelim();
611	<	CompleteSave(can,filename.str());
611	>	CompleteSave(can,ContainerName+"/"+filename.str());
612		cout << "RABS : " << rabs_a << " + " << rabs_b << " * alpha " << endl;
613		}
614
615		MPF_FinalGraph->SetMarkerStyle(21);
616	<	MPF_FinalGraph->SetMarkerColor(kRed);
617	<	MPF_FinalGraph->Fit("pol0","QE"); // quiet, use minos
616	>	MPF_FinalGraph->SetMarkerColor(kBlue);
617	>	MPF_FinalGraph->Fit("pol0",""); // quiet, use minos
618
619		TF1 *mpf_pol0=(TF1*)MPF_FinalGraph->GetFunction("pol0");
620	+	TF1 *mpf_pol0_complete = new TF1("mpf_pol0_complete","[0]+[1]*x");
621	+	mpf_pol0_complete->SetLineWidth(1);
622	+	mpf_pol0_complete->SetLineColor(kBlue);
623	+	mpf_pol0->SetLineColor(kBlue);
624		MPF_FinalGraph->Draw("AP");
625		MPF_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
626		MPF_FinalGraph->GetXaxis()->CenterTitle();
627	<	MPF_FinalGraph->GetYaxis()->SetTitle("data/mc");
627	>	MPF_FinalGraph->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
628		MPF_FinalGraph->GetYaxis()->CenterTitle();
629		MPF_FinalGraph->Draw("AP");
630	+	mpf_pol0_complete->Draw("same");
631
632		stringstream mpf_result;
633		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)}";
#	Line 428 | Line 642 | void new_data_mc_agreement_2d() {
642
643		stringstream filename2;
644		filename2 << "Extrapolation/FINAL_RESULT_MPF";
645	<	CompleteSave(can,filename2.str());
645	>	CompleteSave(can,ContainerName+"/"+filename2.str());
646
647
648		RABS_FinalGraph->SetMarkerStyle(21);
#	Line 438 | Line 652 | void new_data_mc_agreement_2d() {
652		RABS_FinalGraph->Draw("AP");
653		RABS_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
654		RABS_FinalGraph->GetXaxis()->CenterTitle();
655	<	RABS_FinalGraph->GetYaxis()->SetTitle("data/mc");
655	>	RABS_FinalGraph->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc}");
656		RABS_FinalGraph->GetYaxis()->CenterTitle();
657		RABS_FinalGraph->Draw("AP");
658		TF1 *rabs_pol0=(TF1*)RABS_FinalGraph->GetFunction("pol0");
#	Line 456 | Line 670 | void new_data_mc_agreement_2d() {
670
671		filename2.str("");
672		filename2 << "Extrapolation/FINAL_RESULT_RABS";
673	<	CompleteSave(can,filename2.str());
460	<
673	>	CompleteSave(can,ContainerName+"/"+filename2.str());
674
675	<	TCanvas *FaceValueCanvas = new TCanvas("FaceValueCanvas","FaceValueCanvas");
676	<	FaceValueAtPoint3->SetMarkerStyle(21);
677	<	FaceValueAtPoint3->SetMarkerStyle(21);
678	<	FaceValueAtPoint3->Fit("pol0","QE"); // quiet, use minos
679	<	FaceValueAtPoint3->Draw("AP");
680	<	FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
681	<	FaceValueAtPoint3->GetXaxis()->CenterTitle();
682	<	FaceValueAtPoint3->GetYaxis()->SetTitle("data/mc for #alpha<0.3");
683	<	FaceValueAtPoint3->GetYaxis()->CenterTitle();
684	<	FaceValueAtPoint3->Draw("AP");
685	<	TF1 *faceval_pol0=(TF1*)FaceValueAtPoint3->GetFunction("pol0");
675	>	can->cd();
676	>	MPF_FaceValueAtPoint3->SetMarkerStyle(21);
677	>	MPF_FaceValueAtPoint3->SetMarkerColor(kBlue);
678	>	MPF_FaceValueAtPoint3->Fit("pol0","QE"); // quiet, use minos
679	>	MPF_FaceValueAtPoint3->Draw("AP");
680	>	MPF_FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
681	>	MPF_FaceValueAtPoint3->GetXaxis()->CenterTitle();
682	>	MPF_FaceValueAtPoint3->GetYaxis()->SetTitle("<MPF>__{data}/<MPF>_{mc} for #alpha<0.3");
683	>	MPF_FaceValueAtPoint3->GetYaxis()->CenterTitle();
684	>	MPF_FaceValueAtPoint3->Draw("AP");
685	>	TF1 *faceval_pol0=(TF1*)MPF_FaceValueAtPoint3->GetFunction("pol0");
686	>	faceval_pol0->SetLineColor(kBlue);
687	>	faceval_pol0->SetLineWidth(1);
688	>	faceval_pol0->Draw("same");
689	>	TF1 *faceval_pol0_complete = new TF1("faceval_pol0_complete","[0]+[1]*x");
690	>	faceval_pol0_complete->SetParameters(faceval_pol0->GetParameters());
691		stringstream faceval_result;
692		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)}";;
693		TText *faceval_final_mark = write_title(faceval_result.str());
#	Line 483 | Line 701 | void new_data_mc_agreement_2d() {
701		float FaceValResultError=faceval_pol0->GetParError(0)/(faceval_pol0->GetParameter(0)*faceval_pol0->GetParameter(0));
702
703		filename2.str("");
704	<	filename2 << "Extrapolation/FINAL_RESULT_RABS";
705	<	CompleteSave(can,filename2.str());
704	>	filename2 << "Extrapolation/FINAL_RESULT_MPF_FaceValp3";
705	>	CompleteSave(can,ContainerName+"/"+filename2.str());
706	>
707	>	can->cd();
708	>	RABS_FaceValueAtPoint3->SetMarkerStyle(21);
709	>	RABS_FaceValueAtPoint3->SetMarkerColor(kRed);
710	>	RABS_FaceValueAtPoint3->Fit("pol0","QE"); // quiet, use minos
711	>	RABS_FaceValueAtPoint3->Draw("AP");
712	>	RABS_FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
713	>	RABS_FaceValueAtPoint3->GetXaxis()->CenterTitle();
714	>	RABS_FaceValueAtPoint3->GetYaxis()->SetTitle("<R_{abs}>_{data}/<R_{abs}>_{mc} for #alpha<0.3");
715	>	RABS_FaceValueAtPoint3->GetYaxis()->CenterTitle();
716	>	RABS_FaceValueAtPoint3->Draw("AP");
717	>	TF1 *faceval_pol0RABS=(TF1*)RABS_FaceValueAtPoint3->GetFunction("pol0");
718	>	TF1 *faceval_pol0RABS_complete = new TF1("faceval_pol0_complete","[0]+[1]*x");
719	>	faceval_pol0RABS_complete->SetParameters(faceval_pol0RABS->GetParameters());
720	>	faceval_pol0RABS->SetLineWidth(1);
721	>	faceval_pol0RABS->SetLineColor(kRed);
722	>	faceval_pol0RABS_complete->Draw("same");
723	>	stringstream RABSfaceval_result;
724	>	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)}";;
725	>	TText *RABSfaceval_final_mark = write_title(RABSfaceval_result.str());
726	>	RABSfaceval_final_mark->SetX(0.75);
727	>	RABSfaceval_final_mark->SetY(0.75);
728	>	RABSfaceval_final_mark->SetTextSize(0.03);
729	>	RABSfaceval_final_mark->Draw();
730	>	DrawPrelim();
731	>
732	>	float RABSFaceValResult=1/faceval_pol0RABS->GetParameter(0);
733	>	float RABSFaceValResultError=faceval_pol0RABS->GetParError(0)/(faceval_pol0RABS->GetParameter(0)*faceval_pol0RABS->GetParameter(0));
734	>
735	>	filename2.str("");
736	>	filename2 << "Extrapolation/FINAL_RESULT_RABS_FaceValp3";
737	>	CompleteSave(can,ContainerName+"/"+filename2.str());
738
739		cout << "FINAL RESULTS: " << endl;
740		cout << "MPF : " << MPFResult  << " +/- " << MPFResultError  << endl;
741		cout << "Rabs: " << RabsResult << " +/- " << RabsResultError << endl;
742	<	cout << "Face value at 0.3 :" << FaceValResult << " +/- " << FaceValResultError << endl;
742	>	cout << "Face value at 0.3 :" << FaceValResult << " +/- " << FaceValResultError << "  (MPF) " << endl;
743	>	cout << "Face value at 0.3 :" << RABSFaceValResult << " +/- " << RABSFaceValResultError << "  (RABS) " << endl;
744
745		delete can;
746
747	+	write_info(__FUNCTION__,"Need to fill result capsule");
748	+	results.MPF_Result_FaceValue=Value(FaceValResult,FaceValResultError);
749	+	results.Rabs_Result_FaceValue=Value(RABSFaceValResult,RABSFaceValResultError);
750	+	results.MPF_Result_Extrapolated=Value(MPFResult,MPFResultError);
751	+	results.Rabs_Result_Extrapolated=Value(RabsResult,RabsResultError);
752	+	cout << results << endl;
753	+	return results;
754	+	}
755	+
756	+
757	+	void TEST_DoPUStudy(string identifier) {
758	+
759	+	float numVtxcuts[7]={0,10,15,20};
760	+
761	+	TCanvas *can = new TCanvas("can","can");
762	+	TH1F *malpha[8];
763	+	TH1F *dalpha[8];
764	+	cout << identifier << ";";
765	+	cout << endl;
766	+	for(int i=1;i<5;i++) {
767	+	stringstream sSpecialCut;
768	+	if(i<4) {
769	+	sSpecialCut << "numVtx>" << numVtxcuts[i-1] << " && numVtx<" << numVtxcuts[i];
770	+	cout << numVtxcuts[i-1] << "  < nVtx < " << numVtxcuts[i] << ";";
771	+	} else {
772	+	sSpecialCut << "numVtx>" << numVtxcuts[i-1];
773	+	cout << numVtxcuts[i-1] << ";";
774	+	}
775	+
776	+	TCut SpecialCut(sSpecialCut.str().c_str());
777	+
778	+
779	+	malpha[i] = allsamples.Draw("malpha"+any2string(i),"mpf",1000,0,10,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,mc,luminosity);
780	+	dalpha[i] = allsamples.Draw("dalpha"+any2string(i),"mpf",1000,0,10,"MPF","events",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&SpecialCut,data,luminosity);
781	+
782	+
783	+	float a=dalpha[i]->GetMean();
784	+	float b=malpha[i]->GetMean();
785	+	float da=dalpha[i]->GetMeanError();
786	+	float db=malpha[i]->GetMeanError();
787	+	float factor = a / b;
788	+	float error = TMath::Sqrt(  (1/(b*b)) * (da*da) + ((a*a)/(b*b))*db*db);
789	+
790	+	cout << malpha[i]->GetMean() << ";" << malpha[i]->GetMeanError() << ";" << dalpha[i]->GetMean() << ";" << dalpha[i]->GetMeanError() << ";" << malpha[i]->Integral()<<";"<<dalpha[i]->Integral() << ";"<<factor << ";" << error <<  endl;
791	+	malpha[i]->SetLineColor(i+1);
792	+	if(i==0) malpha[i]->Draw("histo");
793	+	else malpha[i]->Draw("histo,same");
794	+	}
795	+
796	+	CompleteSave(can,"AlphaOverview_"+identifier);
797	+	delete can;
798	+
799	+	/*  for(int i=1;i<8;i++) {
800	+	delete malpha[i];
801	+	delete dalpha[i];
802	+	}
803	+	*/
804
805		}
806
807	+	void TEST_ScenarioComparison() {
808	+	//very dumb way to do this.
809	+	TGraphAsymmErrors *gr[5];
810	+	TF1 *fit[5];
811	+	bool dofit=false;
812	+
813	+	TCanvas *can = new TCanvas("can","can",500,500);
814	+
815	+	TLegend *leg = new TLegend(0.2,0.2,0.5,0.45);
816	+	leg->SetBorderSize(0);
817	+	leg->SetFillColor(kWhite);
818	+
819	+	for(int i=0;i<5;i++) {
820	+	gr[i] = new TGraphAsymmErrors();
821	+	double x[4]   = {5,12.5,17.5,20};
822	+	double dxl[4] = {5,2.5,2.5,2.5};
823	+	double dxh[4] = {5,2.5,2.5,10};
824	+
825	+	string name="";
826	+	string color="";
827	+
828	+	if(i==0) {
829	+	double y[4]   = {1.008,1.013,1.003,1.007};
830	+	double dyl[4] = {0.011,0.009,0.013,0.020};
831	+	double dyh[4] = {0.011,0.009,0.013,0.020};
832	+	name="30/30";
833	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
834	+	color="#a01d99";
835	+	}
836	+	if(i==1) {
837	+	double y[4]   = {1.008+0.001,1.013+0.001,1.003+0.001,1.007+0.001};
838	+	double dyl[4] = {0.011,0.009,0.013,0.020};
839	+	double dyh[4] = {0.011,0.009,0.013,0.020};
840	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
841	+	name="30/15";
842	+	color="#2464bc";
843	+	}
844	+	if(i==2) {
845	+	double y[4]   = {1.014,1.026,1.014,1.059};
846	+	double dyl[4] = {0.011,0.010,0.014,0.023};
847	+	double dyh[4] = {0.011,0.010,0.014,0.023};
848	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
849	+	name="15/15";
850	+	color="#f49230";
851	+	}
852	+	if(i==3) {
853	+	double y[4]   = {1.015,1.020,1.008,1.050};
854	+	double dyl[4] = {0.011,0.010,0.014,0.022};
855	+	double dyh[4] = {0.011,0.010,0.014,0.022};
856	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
857	+	name="20/15";
858	+	color="#c72f3c";
859	+	}
860	+	if(i==4) {
861	+	double y[4]   = {1.015+0.001,1.020+0.001,1.008+0.001,1.050+0.001};
862	+	double dyl[4] = {0.011,0.010,0.014,0.022};
863	+	double dyh[4] = {0.011,0.010,0.014,0.022};
864	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
865	+	name="20/20";
866	+	color="#7cb433";
867	+	}
868	+
869	+
870	+	gr[i]->SetTitle();
871	+	gr[i]->GetXaxis()->SetTitle("N(Vtx)");
872	+	gr[i]->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
873	+	gr[i]->GetXaxis()->CenterTitle();
874	+	gr[i]->GetYaxis()->CenterTitle();
875	+	gr[i]->SetName(name.c_str());
876	+	gr[i]->SetMarkerSize(0.00001);
877	+	gr[i]->GetYaxis()->SetRangeUser(0.90,1.10);
878	+	gr[i]->SetLineColor(TColor::GetColor(color.c_str()));
879	+	gr[i]->SetMarkerColor(TColor::GetColor(color.c_str()));
880	+	if(dofit) {
881	+	gr[i]->Fit("pol1");
882	+	fit[i] = (TF1*)gr[i]->GetFunction("pol1");
883	+	fit[i]->SetLineColor(TColor::GetColor(color.c_str()));
884	+	}
885	+	leg->AddEntry(gr[i],name.c_str(),"l");
886	+	if(i==0) gr[i]->Draw("AP*");
887	+	else gr[i]->Draw("P");
888	+	}
889	+	leg->Draw();
890	+	DrawPrelim();
891	+
892	+	CompleteSave(can,"ScenarioComparison");
893	+	delete can;
894	+
895	+	}
896
897	+	void TEST_ScenarioComparisonInclusive() {
898	+	//dumbest way ever to do this. but ok we only need to do it once.
899	+	TGraphAsymmErrors *gr[5];
900	+	TF1 *fit[5];
901	+	bool dofit=true;
902	+
903	+	TCanvas *can = new TCanvas("can","can",500,500);
904	+
905	+	TLegend *leg = new TLegend(0.2,0.2,0.5,0.45);
906	+	leg->SetBorderSize(0);
907	+	leg->SetFillColor(kWhite);
908	+
909	+	for(int i=0;i<5;i++) {
910	+	gr[i] = new TGraphAsymmErrors();
911	+	double x[4]   = {5,12.5,17.5,20};
912	+	double dxl[4] = {5,2.5,2.5,2.5};
913	+	double dxh[4] = {5,2.5,2.5,10};
914	+
915	+	string name="";
916	+	string color="";
917	+
918	+	if(i==0) {
919	+	double y[4]   = {1.01103,1.01496,1.01486,1.02785};
920	+	double dyl[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
921	+	double dyh[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
922	+	name="30/30";
923	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
924	+	color="#a01d99";
925	+	}
926	+	if(i==1) {
927	+	double y[4]   = {1.01103+0.001,1.01496+0.001,1.01486+0.001,1.02785+0.001};
928	+	double dyl[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
929	+	double dyh[4] = {0.00235044,0.00225264,0.00313574,0.00541047};
930	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
931	+	name="30/15";
932	+	color="#2464bc";
933	+	}
934	+	if(i==2) {
935	+	double y[4]   = {1.01575,1.02879,1.0325,1.05456};
936	+	double dyl[4] = {0.00288281,0.00278536,0.00367451,0.00590421};
937	+	double dyh[4] = {0.00288281,0.00278536,0.00367451,0.00590421};
938	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
939	+	name="15/15";
940	+	color="#f49230";
941	+	}
942	+	if(i==3) {
943	+	double y[4]   = {1.01221,1.01922,1.02396,1.04677};
944	+	double dyl[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
945	+	double dyh[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
946	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
947	+	name="20/15";
948	+	color="#c72f3c";
949	+	}
950	+	if(i==4) {
951	+	double y[4]   = {1.01221+0.001,1.01922+0.001,1.02396+0.001,1.04677+0.001};
952	+	double dyl[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
953	+	double dyh[4] = {0.00263821,0.00260853,0.00353673,0.00584728};
954	+	gr[i] = new TGraphAsymmErrors(4,x,y,dxl,dxh,dyl,dyh);
955	+	name="20/20";
956	+	color="#7cb433";
957	+	}
958	+
959	+
960	+	gr[i]->SetTitle();
961	+	gr[i]->GetXaxis()->SetTitle("N(Vtx)");
962	+	gr[i]->GetYaxis()->SetTitle("<MPF>_{data}/<MPF>_{mc}");
963	+	gr[i]->GetXaxis()->CenterTitle();
964	+	gr[i]->GetYaxis()->CenterTitle();
965	+	gr[i]->SetName(name.c_str());
966	+	gr[i]->SetMarkerSize(0.00001);
967	+	gr[i]->GetYaxis()->SetRangeUser(0.90,1.10);
968	+	gr[i]->SetLineColor(TColor::GetColor(color.c_str()));
969	+	gr[i]->SetMarkerColor(TColor::GetColor(color.c_str()));
970	+	if(dofit) {
971	+	gr[i]->Fit("pol1");
972	+	fit[i] = (TF1*)gr[i]->GetFunction("pol1");
973	+	fit[i]->SetLineColor(TColor::GetColor(color.c_str()));
974	+	}
975	+	leg->AddEntry(gr[i],name.c_str(),"l");
976	+	if(i==0) gr[i]->Draw("AP*");
977	+	else gr[i]->Draw("P");
978	+	}
979	+	leg->Draw();
980	+	DrawPrelim();
981	+
982	+	CompleteSave(can,"ScenarioComparisonInclusive");
983	+	delete can;
984	+
985	+	}
986
987	<	void compare_selection(string identifier) {
987	>	void TEST_compare_selection(string identifier) {
988		bool recognized_scenario=false;
989
990		cout << "Running with identifier " << identifier << endl;
991	<	if(identifier=="Zb30") {
992	<	LeadingB=TCut ("Zb30_bTagProbCSVBP[0]>0.679");
993	<	EtaB=TCut("abs(Zb30_pfJetGoodEta[0])<1.3");
994	<	PhiZcut=TCut("abs(Zb30_pfJetDphiZ[0])>2.7");
991	>	if(identifier=="Zb3010") {
992	>	LeadingB=TCut ("Zb3010_bTagProbCSVBP[0]>0.898");
993	>	EtaB=TCut("abs(Zb3010_pfJetGoodEta[0])<1.3");
994	>	PhiZcut=TCut("abs(Zb3010_pfJetDphiZ[0])>2.7");
995		recognized_scenario=true;
996		}
997
998	<	if(identifier=="Zb2030") {
999	<	LeadingB=TCut ("Zb2030_bTagProbCSVBP[0]>0.679");
1000	<	EtaB=TCut("abs(Zb2030_pfJetGoodEta[0])<1.3");
1001	<	PhiZcut=TCut("abs(Zb2030_pfJetDphiZ[0])>2.7");
998	>	if(identifier=="Zb3010") {
999	>	LeadingB=TCut ("Zb3010_bTagProbCSVBP[0]>0.898");
1000	>	EtaB=TCut("abs(Zb3010_pfJetGoodEta[0])<1.3");
1001	>	PhiZcut=TCut("abs(Zb3010_pfJetDphiZ[0])>2.7");
1002	>	recognized_scenario=true;
1003	>	}
1004	>
1005	>	if(identifier=="Zb40") {
1006	>	LeadingB=TCut ("Zb40_bTagProbCSVBP[0]>0.898");
1007	>	EtaB=TCut("abs(Zb40_pfJetGoodEta[0])<1.3");
1008	>	PhiZcut=TCut("abs(Zb40_pfJetDphiZ[0])>2.7");
1009	>	recognized_scenario=true;
1010	>	}
1011	>
1012	>	if(identifier=="Zb301030") {
1013	>	LeadingB=TCut ("Zb301030_bTagProbCSVBP[0]>0.898");
1014	>	EtaB=TCut("abs(Zb301030_pfJetGoodEta[0])<1.3");
1015	>	PhiZcut=TCut("abs(Zb301030_pfJetDphiZ[0])>2.7");
1016		recognized_scenario=true;
1017		}
1018
1019		if(identifier=="Zb1530") {
1020	<	LeadingB=TCut ("Zb1530_bTagProbCSVBP[0]>0.679");
1020	>	LeadingB=TCut ("Zb1530_bTagProbCSVBP[0]>0.898");
1021		EtaB=TCut("abs(Zb1530_pfJetGoodEta[0])<1.3");
1022		PhiZcut=TCut("abs(Zb1530_pfJetDphiZ[0])>2.7");
1023		recognized_scenario=true;
1024		}
1025
1026	<	if(identifier=="Zb3010") {
1027	<	LeadingB=TCut ("Zb3010_bTagProbCSVBP[0]>0.679");
1028	<	EtaB=TCut("abs(Zb3010_pfJetGoodEta[0])<1.3");
1029	<	PhiZcut=TCut("abs(Zb3010_pfJetDphiZ[0])>2.7");
1026	>	if(identifier=="Zb301010") {
1027	>	LeadingB=TCut ("Zb301010_bTagProbCSVBP[0]>0.898");
1028	>	EtaB=TCut("abs(Zb301010_pfJetGoodEta[0])<1.3");
1029	>	PhiZcut=TCut("abs(Zb301010_pfJetDphiZ[0])>2.7");
1030	>	recognized_scenario=true;
1031	>	}
1032	>
1033	>	if(identifier=="Zb1510") {
1034	>	LeadingB=TCut ("Zb1510_bTagProbCSVBP[0]>0.898");
1035	>	EtaB=TCut("abs(Zb1510_pfJetGoodEta[0])<1.3");
1036	>	PhiZcut=TCut("abs(Zb1510_pfJetDphiZ[0])>2.7");
1037	>	recognized_scenario=true;
1038	>	}
1039	>
1040	>	if(identifier=="Zb301010") {
1041	>	LeadingB=TCut ("Zb301010_bTagProbCSVBP[0]>0.898");
1042	>	EtaB=TCut("abs(Zb301010_pfJetGoodEta[0])<1.3");
1043	>	PhiZcut=TCut("abs(Zb301010_pfJetDphiZ[0])>2.7");
1044	>	recognized_scenario=true;
1045	>	}
1046	>
1047	>	if(identifier=="ZbMikko") {
1048	>	LeadingB=TCut ("ZbMikko_bTagProbCSVBP[0]>0.898 && ZbMikko__IsClean");
1049	>	EtaB=TCut("abs(ZbMikko_pfJetGoodEta[0])<1.3 && ZbMikko__IsClean");
1050	>	PhiZcut=TCut("abs(ZbMikko_pfJetDphiZ[0])>2.7 && ZbMikko__IsClean");
1051		recognized_scenario=true;
1052		}
1053
1054	<	if(identifier=="Zb30_SecEta3") {
1055	<	LeadingB=TCut ("Zb30_SecEta3_bTagProbCSVBP[0]>0.679");
1056	<	EtaB=TCut("abs(Zb30_SecEta3_pfJetGoodEta[0])<1.3");
1057	<	PhiZcut=TCut("abs(Zb30_SecEta3_pfJetDphiZ[0])>2.7");
1054	>	if(identifier=="Zb3010_SecEta3") {
1055	>	LeadingB=TCut ("Zb3010_SecEta3_bTagProbCSVBP[0]>0.898");
1056	>	EtaB=TCut("abs(Zb3010_SecEta3_pfJetGoodEta[0])<1.3");
1057	>	PhiZcut=TCut("abs(Zb3010_SecEta3_pfJetDphiZ[0])>2.7");
1058		recognized_scenario=true;
1059		}
1060
1061	<	if(identifier=="Zb30_SecEta5") {
1062	<	LeadingB=TCut ("Zb30_SecEta5_bTagProbCSVBP[0]>0.679");
1063	<	EtaB=TCut("abs(Zb30_SecEta5_pfJetGoodEta[0])<1.3");
1064	<	PhiZcut=TCut("abs(Zb30_SecEta5_pfJetDphiZ[0])>2.7");
1061	>	if(identifier=="Zb3010_SecEta5") {
1062	>	LeadingB=TCut ("Zb3010_SecEta5_bTagProbCSVBP[0]>0.898");
1063	>	EtaB=TCut("abs(Zb3010_SecEta5_pfJetGoodEta[0])<1.3");
1064	>	PhiZcut=TCut("abs(Zb3010_SecEta5_pfJetDphiZ[0])>2.7");
1065		recognized_scenario=true;
1066		}
1067
1068	<	if(identifier=="Zb30_p5Clean") {
1069	<	LeadingB=TCut ("Zb30_p5Clean_bTagProbCSVBP[0]>0.679");
1070	<	EtaB=TCut("abs(Zb30_p5Clean_pfJetGoodEta[0])<1.3");
1071	<	PhiZcut=TCut("abs(Zb30_p5Clean_pfJetDphiZ[0])>2.7");
1068	>	if(identifier=="Zb3010_p5Clean") {
1069	>	LeadingB=TCut ("Zb3010_p5Clean_bTagProbCSVBP[0]>0.898 && Zb3010_p5Clean_IsClean");
1070	>	EtaB=TCut("abs(Zb3010_p5Clean_pfJetGoodEta[0])<1.3 && Zb3010_p5Clean_IsClean");
1071	>	PhiZcut=TCut("abs(Zb3010_p5Clean_pfJetDphiZ[0])>2.7 && Zb3010_p5Clean_IsClean");
1072		recognized_scenario=true;
1073		}
1074
1075	<	if(identifier=="Zb30CHS") {
1076	<	LeadingB=TCut ("Zb30CHS_bTagProbCSVBP[0]>0.679");
1077	<	EtaB=TCut("abs(Zb30CHS_pfJetGoodEta[0])<1.3");
1078	<	PhiZcut=TCut("abs(Zb30CHS_pfJetDphiZ[0])>2.7");
1075	>	if(identifier=="Zb3010CHS") {
1076	>	LeadingB=TCut ("Zb3010CHS_bTagProbCSVBP[0]>0.898");
1077	>	EtaB=TCut("abs(Zb3010CHS_pfJetGoodEta[0])<1.3");
1078	>	PhiZcut=TCut("abs(Zb3010CHS_pfJetDphiZ[0])>2.7");
1079		recognized_scenario=true;
1080		}
1081
#	Line 567 | Line 1088 | void compare_selection(string identifier
1088
1089		cout << endl << endl;
1090
1091	<	SpecialCutFlow(identifier);
1092	<	draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",25,0,200,"Z p_{T}","Official/Zpt__"+identifier,1);
572	<	draw_kin_variable(identifier+"_pfJetGoodPt[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",40,0,200,"Leading Jet Pt (B)","Official/LeadingJetPt__"+identifier,1);
573	<	draw_kin_variable(identifier+"_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",40,0,200,"Sub-Leading Jet Pt [GeV]","Official/Jet2Pt__"+identifier,1);
574	<	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/SubLeadingJetPt_Over_ZPt__"+identifier,1);
575	<	draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<50",20,0,2,"#alpha","Official/alpha_pt_10_to_50__"+identifier,1);
576	<	draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>50&&pt<100",20,0,2,"#alpha","Official/alpha_50_to_100__"+identifier,1);
577	<	draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>100&&pt<200",20,0,2,"#alpha","Official/alpha_pt_100_to_200__"+identifier,1);
578	<	draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>200&&pt<1000",20,0,2,"#alpha","Official/alpha_pt_200_to_1000__"+identifier,1);
579	<	draw_kin_variable(identifier+"_pfBJetDphiZ[0]",ZplusBsel&&LeadingB&&"pt>10&&pt<1000&&pfBJetDphiZ[0]>0",30,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead__"+identifier,0);
1091	>	//    ScenarioComparison();
1092	>	//    ScenarioComparisonInclusive();
1093
1094	<	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","mpf_alpha_smaller_0p3___"+identifier,0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
1095	<	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","ptb_over_ptz___alpha_smaller_0p3______"+identifier,0,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
1094	>
1095	>	//     DoPUStudy(identifier);
1096	>
1097	>	//   SpecialCutFlow(identifier);
1098	>	//    draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",17000,30,200,"Z p_{T}","Official/"+identifier+"/Zpt__"+identifier,1);
1099	>	//   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);
1100	>	//   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);
1101	>	//   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);
1102	>	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);
1103	>	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);
1104	>	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);
1105	>	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);
1106	>	draw_kin_variable(identifier+"_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut,20,0,2,"#alpha","Official/"+identifier+"/alpha_full__"+identifier,1);
1107	>	//   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);
1108	>
1109	>	//   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");
1110	>	//   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");
1111
1112
1113	+	//   */
1114
1115	+	}
1116	+
1117	+	void GetNumberEventsInsideOutsideAlphaWindow(TCut specialcut, string title) {
1118
1119	+	//  TCut cut =  ZplusBsel&&LeadingB&&PhiZcut&&specialcut&&TCut("Zb3010_alpha<0.3");
1120	+	TCut cut = specialcut;
1121	+	TCut in  =  EtaB;
1122	+	TCut out =  TCut("abs(Zb3010_pfJetGoodEta[0])>1.3");
1123	+
1124	+	TH1F *data_IN  = allsamples.Draw("data_IN",    "mll",1,0,150, "nothing [GeV]", "events", cut&&in, data,luminosity);
1125	+	TH1F *data_OUT = allsamples.Draw("data_OUT",   "mll",1,0,150, "nothing [GeV]", "events", cut&&out,data,luminosity);
1126	+	TH1F *mc_IN    = allsamples.Draw("mc_IN",      "mll",1,0,150, "nothing [GeV]", "events", cut&&in, mc  ,luminosity);
1127	+	TH1F *mc_OUT   = allsamples.Draw("mc_OUT",     "mll",1,0,150, "nothing [GeV]", "events", cut&&out,mc  ,luminosity);
1128	+
1129	+	dout << title << " : " << endl;
1130	+	dout << "  Data: In " << data_IN->Integral() << " vs out " << data_OUT->Integral() << endl;
1131	+	dout << "  MC  : In " << mc_IN->Integral()   << " vs out " << mc_OUT->Integral()   << endl;
1132	+	dout << "       ratio in  : " << data_IN->Integral()/mc_IN->Integral() << " +/- " << (data_IN->Integral()/mc_IN->Integral()) * sqrt(1.0/mc_IN->Integral()+ 1.0/data_IN->Integral()) << endl;
1133	+	dout << "       ratio out : " << data_OUT->Integral()/mc_OUT->Integral() << " +/- " << (data_OUT->Integral()/mc_OUT->Integral()) * sqrt(1.0/mc_OUT->Integral()+ 1.0/data_OUT->Integral()) << endl;
1134	+
1135	+	delete data_IN;
1136	+	delete data_OUT;
1137	+	delete mc_IN;
1138	+	delete mc_OUT;
1139		}
1140
1141	<	void do_basic_ZB_analysis() {
1142	<	TCanvas *zbcanvas = new TCanvas("zbcanvas","zbcanvas");
1141	>	void TEST_GetNumberEventsInsideOutsideAlphaWindow() {
1142	>	TCanvas *ca = new TCanvas("ca","ca");
1143	>	//   GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2&&id1==0"), "ee");
1144	>	//   GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2&&id1==1"), "mm");
1145	>	//   GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2"), "ee/mm");
1146	>
1147	>	cout << "BASE SELECTION" << endl;
1148	>	GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2&&id1==0"), "ee");
1149	>	GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2&&id1==1"), "mm");
1150	>	GetNumberEventsInsideOutsideAlphaWindow(TCut("id1==id2"), "ee/mm");
1151	>
1152	>	cout << "BASE SELECTION: Z+b" << endl;
1153	>	GetNumberEventsInsideOutsideAlphaWindow(ZplusBsel&&TCut("id1==id2&&id1==0"), "ee");
1154	>	GetNumberEventsInsideOutsideAlphaWindow(ZplusBsel&&TCut("id1==id2&&id1==1"), "mm");
1155	>	GetNumberEventsInsideOutsideAlphaWindow(ZplusBsel&&TCut("id1==id2"), "ee/mm");
1156	>
1157	>	cout << "Leading B" << endl;
1158	>	GetNumberEventsInsideOutsideAlphaWindow(LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==0"), "ee");
1159	>	GetNumberEventsInsideOutsideAlphaWindow(LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==1"), "mm");
1160	>	GetNumberEventsInsideOutsideAlphaWindow(LeadingB&&ZplusBsel&&TCut("id1==id2"), "ee/mm");
1161	>
1162	>	cout << "PhiZcut" << endl;
1163	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==0"), "ee");
1164	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==1"), "mm");
1165	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2"), "ee/mm");
1166	>
1167	>	cout << "alpha cut" << endl;
1168	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==0")&&TCut("Zb3010_alpha<0.3"), "ee");
1169	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2&&id1==1")&&TCut("Zb3010_alpha<0.3"), "mm");
1170	>	GetNumberEventsInsideOutsideAlphaWindow(PhiZcut&&LeadingB&&ZplusBsel&&TCut("id1==id2")&&TCut("Zb3010_alpha<0.3"), "ee/mm");
1171
1172	<	//    compare_selection("Zb30_p5Clean");
593	<	compare_selection("Zb30CHS");
594	<	compare_selection("Zb30");
595	<	compare_selection("Zb2030");
596	<	compare_selection("Zb1530");
597	<	compare_selection("Zb3010");
598	<	compare_selection("Zb30_SecEta3");
599	<	compare_selection("Zb30_SecEta5"); // buggy, ntuples need to be updated
1172	>	delete ca;
1173
1174
1175	+	}
1176
1177	<	//    print_all_b_yields();
1178	<	//  draw_mpf_vars();
1179	<	//  DrawEvilCutFlow();
1177	>
1178	>
1179	>
1180	>	void do_basic_ZB_analysis(bool DoCompleteAnalysis) {
1181	>
1182	>	//https://twiki.cern.ch/twiki/bin/viewauth/CMS/BtagPOG#2011_Data_and_MC
1183	>	//https://twiki.cern.ch/twiki/pub/CMS/BtagPOG/SFb-ttbar_payload.txt
1184
1185	<	//  draw_Zb_kin_vars();
1185	>	cout << "The lepton requirement is " << (const char*) leptoncut << endl;
1186	>	bool doquick=true;
1187	>
1188	>	TCanvas *zbcanvas = new TCanvas("zbcanvas","zbcanvas");
1189
1190	<	//  new_data_mc_agreement_2d();
1190	>	//     compare_selection("Zb3010_p5Clean");
1191	>	//     compare_selection("Zb3010CHS");
1192	>	//     compare_selection("Zb3010");
1193	>	//     compare_selection("Zb301030");
1194	>	//     compare_selection("Zb1530");
1195	>	//     compare_selection("Zb301010");
1196	>	//     compare_selection("Zb3010_SecEta3");
1197	>	//     compare_selection("Zb3010_SecEta5");
1198	>	//     compare_selection("Zb1510");
1199	>	//     compare_selection("Zb301010");
1200	>	//     compare_selection("ZbMikko");
1201	>	//     compare_selection("Zb3010");
1202	>	//     compare_selection("Zb40");
1203	>
1204	>	if(!doquick) {
1205	>	print_all_b_yields();
1206	>	draw_mpf_vars();
1207	>	//  DrawEvilCutFlow();
1208	>
1209	>	draw_Zb_kin_vars();
1210	>
1211	>	}
1212	>
1213	>	//  GetNumberEventsInsideOutsideAlphaWindow();
1214	>
1215	>	dout <<"THIS IS WHERE IT GETS REALLY INTERESTING!!!!" << endl;
1216	>	if(DoCompleteAnalysis) {
1217	>	/* need to run the following scenarios:
1218	>	* - normal
1219	>	* - inclusive
1220	>	* - medium WP
1221	>	* - alpha up
1222	>	* - alpha down
1223	>	*/
1224	>
1225	>	ZbResultContainer inclusive = new_data_mc_agreement_2d("Zb3010_alpha","inclusive",TCut("Zb3010_bTagProbCSVBP[0]>-10000"),TCut("1.0"));
1226	>	ZbResultContainer Reference = new_data_mc_agreement_2d("Zb3010_alpha","reference",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgt"));
1227	>	ZbResultContainer effmistagUp = new_data_mc_agreement_2d("Zb3010_alpha","effmistagUp",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtUp"));
1228	>	ZbResultContainer effmistagDown = new_data_mc_agreement_2d("Zb3010_alpha","effmistagDown",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgtDown"));
1229	>	ZbResultContainer medium = new_data_mc_agreement_2d("Zb3010_alpha","medium",TCut("Zb3010_bTagProbCSVBP[0]>0.679"),TCut("ZbCHS3010_BTagWgtMedium"));
1230	>	ZbResultContainer AlphaUp = new_data_mc_agreement_2d("Zb3010_alphaUP","AlphaUp",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgt"));
1231	>	ZbResultContainer AlphaDown = new_data_mc_agreement_2d("Zb3010_alphaDOWN","AlphaDown",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgt"));
1232	>
1233	>	LeadingB=TCut("Zb3010_bTagProbCSVBP[0]>0.898");
1234	>	} else {
1235	>	ZbResultContainer Reference = new_data_mc_agreement_2d("Zb3010_alpha","reference",TCut("Zb3010_bTagProbCSVBP[0]>0.898"),TCut("ZbCHS3010_BTagWgt"));
1236	>	}
1237
1238
1239		delete zbcanvas;
#	Line 633 | Line 1260 | void SpecialCutFlow(string identifier) {
1260		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB) << ")";
1261		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB) << ")";
1262		MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut) << ")";
1263	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000")) << ")";
1264	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.3")) << ")";
1265	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.2")) << ")";
1266	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.15")) << ")";
1267	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.1")) << ")";
1268	<	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut("pt>10&&pt<1000&&Zb30_alpha<0.05")) << ")";
1263	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000").c_str())) << ")";
1264	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.3").c_str())) << ")";
1265	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.2").c_str())) << ")";
1266	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.15").c_str())) << ")";
1267	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.1").c_str())) << ")";
1268	>	MegaCut << " + 1*(" << (const char*) (ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&TCut(((string)"pt>10&&pt<1000&&"+identifier+"_alpha<0.05").c_str())) << ")";
1269		MegaCut << " ";
1270
1271		TCanvas *can = new TCanvas("can","can");
1272		can->SetLogy(1);
1273	<	TCut basecut("mll>6||mll<10");
1273	>	TCut basecut(ZplusBsel);
1274
1275
1276

Diff Legend

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