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.2 by buchmann, Tue Jul 10 16:00:42 2012 UTC vs.
Revision 1.5 by buchmann, Wed Aug 15 13:52:13 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 159 | 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 169 | 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");
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;
#	Line 179 | Line 180 | void print_all_b_yields() {
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,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.2, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.15, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.1, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.05, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.3, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.2, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.15, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.1, 10 GeV < p_{T}^{Z} < 1000 GeV");
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,"#alpha<0.05, 10 GeV < p_{T}^{Z} < 1000 GeV");
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		Value get_Zb_data_over_mc(string variable, TCut cut) {
#	Line 232 | Line 233 | Value get_Zb_data_over_mc(string variabl
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	+	int PointThree=-1;
242	+	for(int i=0;i<nalphacuts;i++) {
243	+	if((alphacuts[i]-0.3)<0.01) PointThree=i;
244	+	}
245	+	if(PointThree<0) {
246	+	write_warning(__FUNCTION__,"alpha=0.3 MUST BE ONE OF THE ALPHA VALUES!");
247	+	assert(PointThree<0);
248	+	}
249	+
250		const int nptcuts=5;
251		float ptcuts[nptcuts]={10,50,100,200,1000};
252
#	Line 250 | Line 261 | void new_data_mc_agreement_2d() {
261		stringstream specialcut;
262		specialcut << "((pt>" << ptcuts[ipt] << " && pt< " << ptcuts[ipt+1] << ") && (alpha<" << alphacuts[ia] << "))";
263		Value MPF_data_over_mc = get_Zb_data_over_mc("mpf",TCut(basecut && specialcut.str().c_str()));
264	<	Value RABS_data_over_mc = get_Zb_data_over_mc("pfJetGoodPt[0]/pt",TCut(basecut && specialcut.str().c_str()));
264	>	Value RABS_data_over_mc = get_Zb_data_over_mc("Zb30_pfJetGoodPt[0]/pt",TCut(basecut && specialcut.str().c_str()));
265
266		MPF_Results[ia][ipt]=MPF_data_over_mc.getValue();
267		MPF_Errors[ia][ipt]=MPF_data_over_mc.getError();
#	Line 269 | Line 280 | void new_data_mc_agreement_2d() {
280
281		TGraphErrors *MPF_FinalGraph = new TGraphErrors();
282		TGraphErrors *RABS_FinalGraph = new TGraphErrors();
283	+	TGraphErrors *FaceValueAtPoint3 = new TGraphErrors();
284
285		for(int ipt=0;ipt<nptcuts-1;ipt++) {
286
#	Line 283 | Line 295 | void new_data_mc_agreement_2d() {
295		rabs_gr->SetPoint(ia,alphacuts[ia],RABS_Results[ia][ipt]);
296		rabs_gr->SetPointError(ia,0,RABS_Errors[ia][ipt]);
297		}
298	+	FaceValueAtPoint3->SetPoint(ipt,(ptcuts[ipt]+ptcuts[ipt+1]),MPF_Results[PointThree][ipt]);
299	+	FaceValueAtPoint3->SetPointError(ipt,0,MPF_Errors[PointThree][ipt]);
300
301		mpf_gr->SetMarkerStyle(21);
302		mpf_gr->SetMarkerColor(kBlue);
303		mpf_gr->Draw("AP");
304	<	mpf_gr->Fit("pol1");
304	>	mpf_gr->Fit("pol1","QE"); // quiet, use minos
305	>	mpf_gr->GetXaxis()->SetTitle("#alpha");
306	>	mpf_gr->GetXaxis()->CenterTitle();
307	>	mpf_gr->GetYaxis()->SetTitle("data/mc");
308	>	mpf_gr->GetYaxis()->CenterTitle();
309		TF1 *mpf_pol=(TF1*)mpf_gr->GetFunction("pol1");
310		float mpf_a=mpf_pol->GetParameter(0);
311		float mpf_b=mpf_pol->GetParameter(1);
#	Line 297 | Line 315 | void new_data_mc_agreement_2d() {
315		MPF_ExtraPolatedResults[ipt]=mpf_a;
316
317		stringstream mpf_info;
318	<	mpf_info << "#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (MPF) }}{Extrapolated value: " << std::setprecision(3) << mpf_a << "}";
318	>	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() << "}";
319	>
320		TText *mpf_mark = write_title(mpf_info.str());
321		mpf_mark->SetX(0.75);
322		mpf_mark->SetTextSize(0.03);
#	Line 306 | Line 325 | void new_data_mc_agreement_2d() {
325
326		string filenamebkp=filename.str();
327		filename << "__MPF";
328	+	DrawPrelim();
329		CompleteSave(can,filename.str());
330		cout << "MPF : " << mpf_a << " + " << mpf_b << " * alpha " << endl;
331
#	Line 313 | Line 333 | void new_data_mc_agreement_2d() {
333
334		rabs_gr->SetMarkerStyle(21);
335		rabs_gr->SetMarkerColor(kBlue);
336	+	rabs_gr->Fit("pol1","QE"); // quiet, use minos
337		rabs_gr->Draw("AP");
338	<	rabs_gr->Fit("pol1");
338	>	rabs_gr->GetXaxis()->SetTitle("#alpha");
339	>	rabs_gr->GetXaxis()->CenterTitle();
340	>	rabs_gr->GetYaxis()->SetTitle("data/mc");
341	>	rabs_gr->GetYaxis()->CenterTitle();
342	>
343		TF1 *rabs_pol=(TF1*)rabs_gr->GetFunction("pol1");
344		float rabs_a=rabs_pol->GetParameter(0);
345		float rabs_b=rabs_pol->GetParameter(1);
346		RABS_FinalGraph->SetPoint(ipt,0.5*(ptcuts[ipt]+ptcuts[ipt+1]),rabs_a);
347		RABS_FinalGraph->SetPointError(ipt,0,rabs_pol->GetParError(0));
348	+
349		cout << "!+!+!+!+!+!+!+!+!+ Just added a point to the final plot : " << rabs_a << " +/- " << rabs_pol->GetParError(0) << endl;
350
351
352		stringstream rabs_info;
353	<	rabs_info << "#splitline{#splitline{" << ptcuts[ipt] << " GeV < p_{T}^{Z} < " << ptcuts[ipt+1] << " GeV }{ (R_{abs}) }}{Extrapolated value: " << std::setprecision(3) << rabs_a << "}";
353	>	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() << "}";
354		TText *rabs_mark = write_title(rabs_info.str());
355		rabs_mark->SetX(0.75);
356		rabs_mark->SetTextSize(0.03);
357		rabs_mark->SetY(0.75);
358		rabs_mark->Draw();
359	+
360	+
361
362		RABS_ExtraPolatedResults[ipt]=rabs_a;
363		filename << filenamebkp << "__RABS";
364	+	DrawPrelim();
365		CompleteSave(can,filename.str());
366		cout << "RABS : " << rabs_a << " + " << rabs_b << " * alpha " << endl;
367		}
368
369		MPF_FinalGraph->SetMarkerStyle(21);
370		MPF_FinalGraph->SetMarkerColor(kRed);
371	<	MPF_FinalGraph->Fit("pol0");
371	>	MPF_FinalGraph->Fit("pol0","QE"); // quiet, use minos
372	>
373		TF1 *mpf_pol0=(TF1*)MPF_FinalGraph->GetFunction("pol0");
374		MPF_FinalGraph->Draw("AP");
375	+	MPF_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
376	+	MPF_FinalGraph->GetXaxis()->CenterTitle();
377	+	MPF_FinalGraph->GetYaxis()->SetTitle("data/mc");
378	+	MPF_FinalGraph->GetYaxis()->CenterTitle();
379	+	MPF_FinalGraph->Draw("AP");
380	+
381		stringstream mpf_result;
382	<	mpf_result << "C_{abs}= " << mpf_pol0->GetParameter(0) << " #pm " << mpf_pol0->GetParError(0);
382	>	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)}";
383		TText *rmpf_final_mark = write_title(mpf_result.str());
384		rmpf_final_mark->SetX(0.75);
385		rmpf_final_mark->SetY(0.75);
386		rmpf_final_mark->SetTextSize(0.03);
387		rmpf_final_mark->Draw();
388	+	DrawPrelim();
389	+	float MPFResult=1/mpf_pol0->GetParameter(0);
390	+	float MPFResultError=mpf_pol0->GetParError(0)/(mpf_pol0->GetParameter(0)*mpf_pol0->GetParameter(0));
391
392		stringstream filename2;
393		filename2 << "Extrapolation/FINAL_RESULT_MPF";
#	Line 358 | Line 397 | void new_data_mc_agreement_2d() {
397		RABS_FinalGraph->SetMarkerStyle(21);
398		RABS_FinalGraph->SetMarkerStyle(21);
399		RABS_FinalGraph->SetMarkerColor(kRed);
400	<	RABS_FinalGraph->Fit("pol0");
400	>	RABS_FinalGraph->Fit("pol0","QE"); // quiet, use minos
401	>	RABS_FinalGraph->Draw("AP");
402	>	RABS_FinalGraph->GetXaxis()->SetTitle("p_{T}^{Z}");
403	>	RABS_FinalGraph->GetXaxis()->CenterTitle();
404	>	RABS_FinalGraph->GetYaxis()->SetTitle("data/mc");
405	>	RABS_FinalGraph->GetYaxis()->CenterTitle();
406		RABS_FinalGraph->Draw("AP");
407		TF1 *rabs_pol0=(TF1*)RABS_FinalGraph->GetFunction("pol0");
408		stringstream rabs_result;
409	<	rabs_result << "C_{abs}= " << rabs_pol0->GetParameter(0) << " #pm " << rabs_pol0->GetParError(0);
409	>	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})}";;
410		TText *rabs_final_mark = write_title(rabs_result.str());
411		rabs_final_mark->SetX(0.75);
412		rabs_final_mark->SetY(0.75);
413		rabs_final_mark->SetTextSize(0.03);
414		rabs_final_mark->Draw();
415	+	DrawPrelim();
416	+
417	+	float RabsResult=1/rabs_pol0->GetParameter(0);
418	+	float RabsResultError=rabs_pol0->GetParError(0)/(rabs_pol0->GetParameter(0)*rabs_pol0->GetParameter(0));
419	+
420	+	filename2.str("");
421	+	filename2 << "Extrapolation/FINAL_RESULT_RABS";
422	+	CompleteSave(can,filename2.str());
423	+
424	+
425	+	TCanvas *FaceValueCanvas = new TCanvas("FaceValueCanvas","FaceValueCanvas");
426	+	FaceValueAtPoint3->SetMarkerStyle(21);
427	+	FaceValueAtPoint3->SetMarkerStyle(21);
428	+	FaceValueAtPoint3->Fit("pol0","QE"); // quiet, use minos
429	+	FaceValueAtPoint3->Draw("AP");
430	+	FaceValueAtPoint3->GetXaxis()->SetTitle("p_{T}^{Z}");
431	+	FaceValueAtPoint3->GetXaxis()->CenterTitle();
432	+	FaceValueAtPoint3->GetYaxis()->SetTitle("data/mc for #alpha<0.3");
433	+	FaceValueAtPoint3->GetYaxis()->CenterTitle();
434	+	FaceValueAtPoint3->Draw("AP");
435	+	TF1 *faceval_pol0=(TF1*)FaceValueAtPoint3->GetFunction("pol0");
436	+	stringstream faceval_result;
437	+	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)}";;
438	+	TText *faceval_final_mark = write_title(faceval_result.str());
439	+	faceval_final_mark->SetX(0.75);
440	+	faceval_final_mark->SetY(0.75);
441	+	faceval_final_mark->SetTextSize(0.03);
442	+	faceval_final_mark->Draw();
443	+	DrawPrelim();
444	+
445	+	float FaceValResult=1/faceval_pol0->GetParameter(0);
446	+	float FaceValResultError=faceval_pol0->GetParError(0)/(faceval_pol0->GetParameter(0)*faceval_pol0->GetParameter(0));
447
448		filename2.str("");
449		filename2 << "Extrapolation/FINAL_RESULT_RABS";
450		CompleteSave(can,filename2.str());
451
452	+	cout << "FINAL RESULTS: " << endl;
453	+	cout << "MPF : " << MPFResult  << " +/- " << MPFResultError  << endl;
454	+	cout << "Rabs: " << RabsResult << " +/- " << RabsResultError << endl;
455	+	cout << "Face value at 0.3 :" << FaceValResult << " +/- " << FaceValResultError << endl;
456
457		delete can;
458
#	Line 386 | Line 466 | void do_basic_ZB_analysis() {
466		write_warning(__FUNCTION__,"TESTING NEW DATA MC AGREEMENT");
467
468		/*
389	–	write_warning(__FUNCTION__,"Not doing b yields (already implemented)");
469		//  print_all_b_yields();
391	–	write_warning(__FUNCTION__,"Not drawing mpf variables (already implemented)");
470		//  draw_mpf_vars();
393	–	write_warning(__FUNCTION__,"Not drawing mpf variables (already implemented)");
471		draw_Zb_kin_vars();
472
473		data_mc_agreement_2d();

Diff Legend

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