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.1 by buchmann, Mon Jun 25 13:34:21 2012 UTC vs.
Revision 1.4 by buchmann, Fri Aug 10 11:47:53 2012 UTC

#	Line 24 | Line 24 | namespace ZbData {
24		vector<float> data_over_mc;
25
26		TCut ZplusBsel("pt1>20&&pt2>20&&mll>60&&mll<120");
27	<	TCut LeadingB("bTagProbCSVBP[0]>0.679");
28	<	TCut EtaB("abs(pfJetGoodEta[0])<1.3");
29	<	TCut PhiZcut("pfJetDphiZ[0]>2.7");
27	>	TCut LeadingB("Zb30_bTagProbCSVBP[0]>0.679");
28	>	TCut EtaB("abs(Zb30_pfJetGoodEta[0])<1.3");
29	>	TCut PhiZcut("Zb30_pfJetDphiZ[0]>2.7");
30
31		}
32
#	Line 36 | Line 36 | void print_yield(TCut cut) {
36		TH1F *data  = allsamples.Draw("data",    "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,data,luminosity);
37		cout << "  Yields for " << (const char*) cut << endl;
38		cout << "     data: " << data->Integral() << endl;
39	<	THStack mcm = allsamples.DrawStack("mc", "mll",1,0,10000000, "m_{ll} [GeV}", "events", cut,mc,luminosity);
39	>	THStack mcm = allsamples.DrawStack("mc", "mll",1,0,10000000, "m_{ll} [GeV]", "events", cut,mc,luminosity);
40		int nhists = mcm.GetHists()->GetEntries();
41		cout << "Going to loop over " << nhists << " histograms " << endl;
42		TFile *test = new TFile("test.root","RECREATE");
#	Line 50 | Line 50 | void print_yield(TCut cut) {
50		delete data;
51		}
52
53	<	void draw_kin_variable(string variable, TCut cut, int nbins, float min, float max, string xlabel, string saveas, bool dologscale) {
53	>	void draw_kin_variable(string variable, TCut cut, int nbins, float min, float max, string xlabel, string saveas, bool dologscale, string speciallabel="") {
54
55		TCanvas *ckin = new TCanvas("ckin","kin variable canvas");
56		ckin->SetLogy(dologscale);
57		TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity);
58		datahisto->SetMarkerSize(DataMarkerSize);
59		THStack mcstack = allsamples.DrawStack("mcstack",variable,nbins,min,max, xlabel, "events",cut,mc,luminosity);
60	<	datahisto->SetMaximum(5.3*datahisto->GetMaximum());
60	>	if(dologscale) datahisto->SetMaximum(3*datahisto->GetMaximum());
61	>	else datahisto->SetMaximum(1.5*datahisto->GetMaximum());
62	>	TText *speciall;
63	>	if(speciallabel!="") {
64	>	speciall = new TText(0.2,0.8,speciallabel.c_str());
65	>	speciall->Draw();
66	>	}
67		datahisto->Draw("e1");
68		ckin->Update();
69		mcstack.Draw("same");
#	Line 84 | Line 90 | void draw_kin_variable(string variable,
90
91		datahisto->Delete();
92		delete ckin;
87	–
88	–
89	–
90	–	write_warning(__FUNCTION__,"Not implemented yet");
91	–	write_warning(__FUNCTION__,"Need to store overall data_over_mc");
93		}
94
95		Value getfrom2Dmap(TH2F *map, int ixbin, int iybin) {
#	Line 158 | Line 159 | void data_mc_agreement_2d() {
159
160		stringstream saveas;
161		saveas << "Zplusb_data_over_mc___" << binningx[ix-1] << "to" << binningx[ix];
162	+	DrawPrelim();
163		CompleteSave(agrcan,saveas.str());
164		//    agrcan->SaveAs("tester.png");
165		}
#	Line 168 | Line 170 | void data_mc_agreement_2d() {
170
171		void print_all_b_yields() {
172		cout << "Basic selection with a b jet" << endl;
173	<	print_yield(ZplusBsel&&"pfJetGoodNumBtag>0");
174	<	cout << "Leading jet is a b " << endl;
173	>	print_yield(ZplusBsel&&"Zb30_pfJetGoodNumBtag>0");
174	>	cout << "Leading jet is a b-jet " << endl;
175		print_yield(ZplusBsel&&LeadingB);
176		cout << "|#eta_{b}|<1.3 " << endl;
177		print_yield(ZplusBsel&&LeadingB&&EtaB);
178		cout << "|#delta#phi(b<Z)|>2.7" << endl;
179		print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
180	<	cout << "10<ptZ<100 GeV" << endl;
180	>	cout << "10<ptZ<1000 GeV" << endl;
181		print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000");
182		cout << "#alpha < 0.3" << endl;
183	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.3");
183	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.3");
184		cout << "#alpha < 0.2" << endl;
185	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.2");
185	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.2");
186		cout << "#alpha < 0.15" << endl;
187	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.15");
187	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.15");
188		cout << "#alpha < 0.1" << endl;
189	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.1");
189	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.1");
190		cout << "#alpha < 0.05" << endl;
191	<	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.05");
191	>	print_yield(ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"Zb30_alpha<0.05");
192		}
193
194		void draw_Zb_kin_vars() {
195	<	draw_kin_variable("pfBJetDphiZ",ZplusBsel&&LeadingB,100,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",1);
195	>	draw_kin_variable("Zb30_pfBJetDphiZ",ZplusBsel&&LeadingB,100,0,3.2,"#delta#phi (Z,b lead)","DeltaPhiZBlead",1);
196		draw_kin_variable("pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",50,0,200,"Z p_{T}","Zpt",1);
197	<	draw_kin_variable("pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",50,0,200,"Sub-Leading Jet Pt","Jet2Pt",1);
198	<	draw_kin_variable("pfJetGoodPtlBtag[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",50,0,200,"Leading Jet Pt (B)","LeadingJetPt",1);
199	<	draw_kin_variable("alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,2,"#alpha","alpha",1);
197	>	draw_kin_variable("Zb30_pfJetGoodPt[1]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",50,0,200,"Sub-Leading Jet Pt","Jet2Pt",1);
198	>	draw_kin_variable("Zb30_pfJetGoodPtlBtag[0]",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",50,0,200,"Leading Jet Pt (B)","LeadingJetPt",1);
199	>	draw_kin_variable("Zb30_alpha",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000",20,0,2,"#alpha","alpha",1);
200		}
201
202		void draw_mpf_vars() {
203		ZbData::data_over_mc.clear();
204	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.3",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p3",0);
205	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.2",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p2",0);
206	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.15",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p15",0);
207	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.1",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p1",0);
208	<	draw_kin_variable("mpf",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.05",20,0,2,"Z+b MPF","mpf_alpha_smaller_0p05",0);
209	<
210	<	draw_kin_variable("pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.3",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p3",0);
211	<	draw_kin_variable("pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.2",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p2",0);
212	<	draw_kin_variable("pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.15",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p15",0);
213	<	draw_kin_variable("pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.1",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz___alpha_smaller_0p1",0);
214	<	draw_kin_variable("pfJetGoodPt[0]/pt",ZplusBsel&&LeadingB&&EtaB&&PhiZcut&&"pt>10&&pt<1000"&&"alpha<0.05",20,0,2,"p_{T} b-jet / p_{T} Z","ptb_over_ptz__mpf_alpha_smaller_0p05",0);
204	>	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");
205	>	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");
206	>	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");
207	>	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");
208	>	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");
209	>
210	>	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");
211	>	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");
212	>	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");
213	>	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");
214	>	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");
215		}
216
217	<	void do_basic_ZB_analysis() {
217	>	Value get_Zb_data_over_mc(string variable, TCut cut) {
218	>	//  write_warning(__FUNCTION__,"Debugging this function, therefore always returning 3.1415+/-1.010101"); return Value(3.1415,1.010101);
219	>	TH1F *hdata  = allsamples.Draw("hdata",variable,1000,0,30, "#alpha [GeV]", "events", cut,data,luminosity);
220	>	TH1F *hmc    = allsamples.Draw("hmc" , variable,1000,0,30, "#alpha [GeV]", "events", cut,mc,  luminosity);
221	>
222	>	float a=hdata->GetMean();
223	>	float b=hmc->GetMean();
224	>	float da=hdata->GetMeanError();
225	>	float db=hmc->GetMeanError();
226	>	float factor = a / b;
227	>	float error = TMath::Sqrt(  (1/(b*b)) * (da*da) + ((a*a)/(b*b))*db*db);
228	>	delete hdata;
229	>	delete hmc;
230	>	return Value(factor,error);
231	>	}
232	>
233	>
234	>
235	>	void new_data_mc_agreement_2d() {
236	>	gStyle->SetOptFit(0);
237	>	TCut basecut(ZplusBsel&&LeadingB&&EtaB&&PhiZcut);
238	>
239	>	const int nalphacuts=5;
240	>	float alphacuts[nalphacuts] = {0.1,0.15,0.2,0.25,0.3};
241	>	const int nptcuts=5;
242	>	float ptcuts[nptcuts]={10,50,100,200,1000};
243	>
244	>	float MPF_Results[nalphacuts][nptcuts];
245	>	float MPF_Errors[nalphacuts][nptcuts];
246	>	float RABS_Results[nalphacuts][nptcuts];
247	>	float RABS_Errors[nalphacuts][nptcuts];
248	>
249	>
250	>	for(int ia=0;ia<nalphacuts;ia++) {
251	>	for(int ipt=0;ipt<nptcuts-1;ipt++) {
252	>	stringstream specialcut;
253	>	specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (alpha<" << alphacuts[ia] << "))";
254	>	Value MPF_data_over_mc = get_Zb_data_over_mc("mpf",TCut(basecut && specialcut.str().c_str()));
255	>	Value RABS_data_over_mc = get_Zb_data_over_mc("Zb30_pfJetGoodPt[0]/pt",TCut(basecut && specialcut.str().c_str()));
256	>
257	>	MPF_Results[ia][ipt]=MPF_data_over_mc.getValue();
258	>	MPF_Errors[ia][ipt]=MPF_data_over_mc.getError();
259	>
260	>	RABS_Results[ia][ipt]=RABS_data_over_mc.getValue();
261	>	RABS_Errors[ia][ipt]=RABS_data_over_mc.getError();
262	>
263	>	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;
264	>	}
265	>	}
266	>
267	>	float MPF_ExtraPolatedResults[nptcuts];
268	>	float RABS_ExtraPolatedResults[nptcuts];
269	>
270	>	TCanvas *can = new TCanvas("can");
271	>
272	>	TGraphErrors *MPF_FinalGraph = new TGraphErrors();
273	>	TGraphErrors *RABS_FinalGraph = new TGraphErrors();
274	>
275	>	for(int ipt=0;ipt<nptcuts-1;ipt++) {
276	>
277	>	stringstream filename;
278	>	filename << "Extrapolation/Zplusb_data_over_mc___" << ptcuts[ipt] << "_to_" << ptcuts[ipt+1];
279	>	cout << "Going to create histo for " << filename.str() << endl;
280	>	TGraphErrors *mpf_gr =new TGraphErrors();
281	>	TGraphErrors *rabs_gr =new TGraphErrors();
282	>	for(int ia=0;ia<nalphacuts;ia++) {
283	>	mpf_gr->SetPoint(ia,alphacuts[ia],MPF_Results[ia][ipt]);
284	>	mpf_gr->SetPointError(ia,0,MPF_Errors[ia][ipt]);
285	>	rabs_gr->SetPoint(ia,alphacuts[ia],RABS_Results[ia][ipt]);
286	>	rabs_gr->SetPointError(ia,0,RABS_Errors[ia][ipt]);
287	>	}
288	>
289	>	mpf_gr->SetMarkerStyle(21);
290	>	mpf_gr->SetMarkerColor(kBlue);
291	>	mpf_gr->Draw("AP");
292	>	mpf_gr->Fit("pol1","QE"); // quiet, use minos
293	>	mpf_gr->GetXaxis()->SetTitle("#alpha");
294	>	mpf_gr->GetXaxis()->CenterTitle();
295	>	mpf_gr->GetYaxis()->SetTitle("data/mc");
296	>	mpf_gr->GetYaxis()->CenterTitle();
297	>	TF1 *mpf_pol=(TF1*)mpf_gr->GetFunction("pol1");
298	>	float mpf_a=mpf_pol->GetParameter(0);
299	>	float mpf_b=mpf_pol->GetParameter(1);
300	>	MPF_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),mpf_a);
301	>	MPF_FinalGraph->SetPointError(ipt,0,mpf_pol->GetParError(0));
302	>
303	>	MPF_ExtraPolatedResults[ipt]=mpf_a;
304	>
305	>	stringstream mpf_info;
306	>	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() << "}";
307	>
308	>	TText *mpf_mark = write_title(mpf_info.str());
309	>	mpf_mark->SetX(0.75);
310	>	mpf_mark->SetTextSize(0.03);
311	>	mpf_mark->SetY(0.75);
312	>	mpf_mark->Draw();
313	>
314	>	string filenamebkp=filename.str();
315	>	filename << "__MPF";
316	>	DrawPrelim();
317	>	CompleteSave(can,filename.str());
318	>	cout << "MPF : " << mpf_a << " + " << mpf_b << " * alpha " << endl;
319	>
320	>	filename.str("");
321	>
322	>	rabs_gr->SetMarkerStyle(21);
323	>	rabs_gr->SetMarkerColor(kBlue);
324	>	rabs_gr->Fit("pol1","QE"); // quiet, use minos
325	>	rabs_gr->Draw("AP");
326	>	rabs_gr->GetXaxis()->SetTitle("#alpha");
327	>	rabs_gr->GetXaxis()->CenterTitle();
328	>	rabs_gr->GetYaxis()->SetTitle("data/mc");
329	>	rabs_gr->GetYaxis()->CenterTitle();
330	>
331	>	TF1 *rabs_pol=(TF1*)rabs_gr->GetFunction("pol1");
332	>	float rabs_a=rabs_pol->GetParameter(0);
333	>	float rabs_b=rabs_pol->GetParameter(1);
334	>	RABS_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),rabs_a);
335	>	RABS_FinalGraph->SetPointError(ipt,0,rabs_pol->GetParError(0));
336	>
337	>	cout << "!+!+!+!+!+!+!+!+!+ Just added a point to the final plot : " << rabs_a << " +/- " << rabs_pol->GetParError(0) << endl;
338	>
339	>
340	>	stringstream rabs_info;
341	>	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() << "}";
342	>	TText *rabs_mark = write_title(rabs_info.str());
343	>	rabs_mark->SetX(0.75);
344	>	rabs_mark->SetTextSize(0.03);
345	>	rabs_mark->SetY(0.75);
346	>	rabs_mark->Draw();
347	>
348	>
349	>
350	>	RABS_ExtraPolatedResults[ipt]=rabs_a;
351	>	filename << filenamebkp << "__RABS";
352	>	DrawPrelim();
353	>	CompleteSave(can,filename.str());
354	>	cout << "RABS : " << rabs_a << " + " << rabs_b << " * alpha " << endl;
355	>	}
356	>
357	>	MPF_FinalGraph->SetMarkerStyle(21);
358	>	MPF_FinalGraph->SetMarkerColor(kRed);
359	>	MPF_FinalGraph->Fit("pol0","QE"); // quiet, use minos
360	>
361	>	TF1 *mpf_pol0=(TF1*)MPF_FinalGraph->GetFunction("pol0");
362	>	MPF_FinalGraph->Draw("AP");
363	>	MPF_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
364	>	MPF_FinalGraph->GetXaxis()->CenterTitle();
365	>	MPF_FinalGraph->GetYaxis()->SetTitle("data/mc");
366	>	MPF_FinalGraph->GetYaxis()->CenterTitle();
367	>	MPF_FinalGraph->Draw("AP");
368	>
369	>	stringstream mpf_result;
370	>	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)}";
371	>	TText *rmpf_final_mark = write_title(mpf_result.str());
372	>	rmpf_final_mark->SetX(0.75);
373	>	rmpf_final_mark->SetY(0.75);
374	>	rmpf_final_mark->SetTextSize(0.03);
375	>	rmpf_final_mark->Draw();
376	>	DrawPrelim();
377	>
378	>	stringstream filename2;
379	>	filename2 << "Extrapolation/FINAL_RESULT_MPF";
380	>	CompleteSave(can,filename2.str());
381	>
382	>
383	>	RABS_FinalGraph->SetMarkerStyle(21);
384	>	RABS_FinalGraph->SetMarkerStyle(21);
385	>	RABS_FinalGraph->SetMarkerColor(kRed);
386	>	RABS_FinalGraph->Fit("pol0","QE"); // quiet, use minos
387	>	RABS_FinalGraph->Draw("AP");
388	>	RABS_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
389	>	RABS_FinalGraph->GetXaxis()->CenterTitle();
390	>	RABS_FinalGraph->GetYaxis()->SetTitle("data/mc");
391	>	RABS_FinalGraph->GetYaxis()->CenterTitle();
392	>	RABS_FinalGraph->Draw("AP");
393	>	TF1 *rabs_pol0=(TF1*)RABS_FinalGraph->GetFunction("pol0");
394	>	stringstream rabs_result;
395	>	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})}";;
396	>	TText *rabs_final_mark = write_title(rabs_result.str());
397	>	rabs_final_mark->SetX(0.75);
398	>	rabs_final_mark->SetY(0.75);
399	>	rabs_final_mark->SetTextSize(0.03);
400	>	rabs_final_mark->Draw();
401	>	DrawPrelim();
402	>
403	>	filename2.str("");
404	>	filename2 << "Extrapolation/FINAL_RESULT_RABS";
405	>	CompleteSave(can,filename2.str());
406	>
407	>
408	>	delete can;
409
410
411	+	}
412	+
413	+
414	+
415	+	void do_basic_ZB_analysis() {
416	+	TCanvas *zbcanvas = new TCanvas("zbcanvas","zbcanvas");
417	+	write_warning(__FUNCTION__,"TESTING NEW DATA MC AGREEMENT");
418	+
419	+	/*
420	+	//  print_all_b_yields();
421	+	//  draw_mpf_vars();
422	+	draw_Zb_kin_vars();
423
424		data_mc_agreement_2d();
425	+	*/
426	+	new_data_mc_agreement_2d();
427
428	<	write_warning(__FUNCTION__,"Replaced B jet dphi cut with jet dphi cut - need to put this back!");
428	>	delete zbcanvas;
429		}

Diff Legend

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