diff --git a/macros/BeCorrPRC/figPresBestFitLeftES.cxx b/macros/BeCorrPRC/figPresBestFitLeftES.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..e44161a61314599f7b13fd8035a141a529e7546a
--- /dev/null
+++ b/macros/BeCorrPRC/figPresBestFitLeftES.cxx
@@ -0,0 +1,360 @@
+//canvas c1: new Figure for revised version
+
+using std::cout;
+using std::endl;
+
+void figPresBestFitLeftES(Bool_t saveFigs = 0)
+{
+
+	const Int_t canMap[6] = {1,4,3,6,2,5};
+
+
+	TH1F *hCurr = 0;
+	TString hName;
+
+	TFile *fr2 = new TFile("figures/PRCthetaA.root");
+
+	TH1F *hAlphaExp7590[6];
+	TH1F *hAlphaSim7590[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT3_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp7590[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim7590[i] = new TH1F(*hCurr);
+	}
+
+	TH1F *hAlphaExp90120[6];
+	TH1F *hAlphaSim90120[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT3_%d_90120", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp90120[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_90120", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim90120[i] = new TH1F(*hCurr);
+	}
+
+	TFile *fr4560 = new TFile("figures/BestFit4560.root");
+
+	TH1F *hAlphaExp4560[6];
+	TH1F *hAlphaSim4560[6];
+	TH1F *hAlphaExp4560low[6];
+	TH1F *hAlphaSim4560low[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT3_%d_4560", i);
+		hCurr = (TH1F*)fr4560->Get(hName.Data());
+		hAlphaExp4560[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_4560", i);
+		hCurr = (TH1F*)fr4560->Get(hName.Data());
+		hAlphaSim4560[i] = new TH1F(*hCurr);
+
+		hName.Form("hethetaAT2_%d_4560", i);
+		hCurr = (TH1F*)fr4560->Get(hName.Data());
+		hAlphaExp4560low[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT2_%d_4560", i);
+		hCurr = (TH1F*)fr4560->Get(hName.Data());
+		hAlphaSim4560low[i] = new TH1F(*hCurr);
+	}
+
+	TFile *fr90120 = new TFile("figures/BestFit90120.root");
+
+	TH1F *hAlphaExp90120low[6];
+	TH1F *hAlphaSim90120low[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT2_%d_90120", i);
+		hCurr = (TH1F*)fr90120->Get(hName.Data());
+		hAlphaExp90120low[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT2_%d_90120", i);
+		hCurr = (TH1F*)fr90120->Get(hName.Data());
+		hAlphaSim90120low[i] = new TH1F(*hCurr);
+	}
+
+	TFile *fr7590 = new TFile("figures/BestFit7590.root");
+
+	TH1F *hAlphaExp7590low[6];
+	TH1F *hAlphaSim7590low[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT2_%d_7590", i);
+		hCurr = (TH1F*)fr7590->Get(hName.Data());
+		hAlphaExp7590low[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT2_%d_7590", i);
+		hCurr = (TH1F*)fr7590->Get(hName.Data());
+		hAlphaSim7590low[i] = new TH1F(*hCurr);
+	}
+
+
+	TFile *fr6075 = new TFile("figures/BestFit6075.root");
+
+	TH1F *hAlphaExp6075[6];
+	TH1F *hAlphaSim6075[6];
+	TH1F *hAlphaExp6075low[6];
+	TH1F *hAlphaSim6075low[6];
+
+	for (Int_t i = 0; i < 6; i++) {
+		hName.Form("hethetaAT3_%d_6075", i);
+		hCurr = (TH1F*)fr6075->Get(hName.Data());
+		hAlphaExp6075[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_6075", i);
+		hCurr = (TH1F*)fr6075->Get(hName.Data());
+		hAlphaSim6075[i] = new TH1F(*hCurr);
+
+		hName.Form("hethetaAT2_%d_6075", i);
+		hCurr = (TH1F*)fr6075->Get(hName.Data());
+		hAlphaExp6075low[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT2_%d_6075", i);
+		hCurr = (TH1F*)fr6075->Get(hName.Data());
+		hAlphaSim6075low[i] = new TH1F(*hCurr);
+	}
+
+////////////////////////
+//	2.5 - 3.1 MeV
+////////////////////////
+
+
+	TCanvas *c1 = new TCanvas("c1", "Combined fits  2.5 - 3.1 MeV", 1000, 1000);
+	c1->ToggleEditor();
+	c1->Divide(2, 2, .001, .001);
+
+	cout << "new Figure" << endl;
+
+	c1->cd(1);
+	TH1F *hAlphaComb4560 = new TH1F(*hAlphaSim4560[3]);
+	//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb4560->Add(hAlphaSim4560[2], 0.1);
+	hAlphaComb4560->Scale(0.9);
+
+	hAlphaComb4560->SetTitle("");
+	hAlphaComb4560->GetYaxis()->SetRangeUser(0., 480.);
+	hAlphaComb4560->GetYaxis()->SetTitle("");
+	hAlphaComb4560->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb4560->GetXaxis()->SetNdivisions(506);
+	hAlphaComb4560->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb4560->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb4560->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb4560->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb4560->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb4560->SetFillColor(kGray+1);
+	hAlphaComb4560->SetLineColor(kGray+1);
+
+	hAlphaComb4560->Draw("same");
+	hAlphaExp4560[4]->Draw("E same");
+
+	hAlphaExp4560[4]->Chi2Test(hAlphaComb4560, "P");
+
+	c1->cd(2);
+	TH1F *hAlphaComb6075 = new TH1F(*hAlphaSim6075[3]);
+//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb6075->Add(hAlphaSim6075[2], 0.3);
+	hAlphaComb6075->Scale(0.8);
+
+	hAlphaComb6075->SetTitle("");
+	hAlphaComb6075->GetYaxis()->SetRangeUser(0., 650.);
+	hAlphaComb6075->GetYaxis()->SetTitle("");
+	hAlphaComb6075->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb6075->GetXaxis()->SetNdivisions(506);
+	hAlphaComb6075->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb6075->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb6075->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb6075->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb6075->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb6075->SetFillColor(kGray+1);
+	hAlphaComb6075->SetLineColor(kGray+1);
+
+	hAlphaComb6075->Draw("same");
+	hAlphaExp6075[4]->Draw("E same");
+
+	hAlphaExp6075[4]->Chi2Test(hAlphaComb6075, "P");
+
+	c1->cd(3);
+
+	TH1F *hAlphaComb7590 = new TH1F(*hAlphaSim7590[3]);
+//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb7590->Add(hAlphaSim7590[2], 0.2);
+	hAlphaComb7590->Scale(0.9);
+
+	hAlphaComb7590->SetTitle("");
+	hAlphaComb7590->GetYaxis()->SetRangeUser(0., 499.);
+	hAlphaComb7590->GetYaxis()->SetTitle("");
+	hAlphaComb7590->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb7590->GetXaxis()->SetNdivisions(506);
+	hAlphaComb7590->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb7590->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb7590->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb7590->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb7590->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb7590->SetFillColor(kGray+1);
+	hAlphaComb7590->SetLineColor(kGray+1);
+
+	hAlphaComb7590->Draw("same");
+	hAlphaExp7590[4]->Draw("E same");
+
+	hAlphaExp7590[4]->Chi2Test(hAlphaComb7590, "P");
+
+	c1->cd(4);
+
+	TH1F *hAlphaComb90120 = new TH1F(*hAlphaSim90120[5]);
+//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb90120->Add(hAlphaSim90120[3], 1);
+	hAlphaComb90120->Scale(0.5);
+
+	hAlphaComb90120->SetTitle("");
+	hAlphaComb90120->GetYaxis()->SetRangeUser(0., 399.);
+	hAlphaComb90120->GetYaxis()->SetTitle("");
+	hAlphaComb90120->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb90120->GetXaxis()->SetNdivisions(506);
+	hAlphaComb90120->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb90120->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb90120->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb90120->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb90120->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb90120->SetFillColor(kGray+1);
+	hAlphaComb90120->SetLineColor(kGray+1);
+
+	hAlphaComb90120->Draw("same");
+	hAlphaExp90120[4]->Draw("E same");
+
+	hAlphaExp90120[4]->Chi2Test(hAlphaComb90120, "P");
+
+//	return;
+
+////////////////////////
+//	1.9 - 2.5 MeV
+////////////////////////
+
+	TCanvas *c2 = new TCanvas("c2", "Combined fits 1.9 - 2.5 MeV", 1000, 1000);
+	c2->ToggleEditor();
+	c2->Divide(2, 2, .001, .001);
+
+	cout << endl << endl << endl << "new Figure for lower energy interval" << endl;
+
+	c2->cd(1);
+	TH1F *hAlphaComb4560low = new TH1F(*hAlphaSim4560low[5]);
+	//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb4560low->Add(hAlphaSim4560low[4], 0.1);
+	hAlphaComb4560low->Add(hAlphaSim4560low[2], 0.1);
+	hAlphaComb4560low->Add(hAlphaSim4560low[3], 1);
+	hAlphaComb4560low->Scale(0.5);
+
+	hAlphaComb4560low->SetTitle("");
+	hAlphaComb4560low->GetYaxis()->SetRangeUser(0., 300.);
+	hAlphaComb4560low->GetYaxis()->SetTitle("");
+	hAlphaComb4560low->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb4560low->GetXaxis()->SetNdivisions(506);
+	hAlphaComb4560low->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb4560low->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb4560low->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb4560low->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb4560low->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb4560low->SetLineColor(kGray+1);
+	hAlphaComb4560low->SetFillColor(kGray+1);
+
+	hAlphaComb4560low->Draw("");
+	hAlphaExp4560low[4]->Draw("E same");
+
+	hAlphaExp4560low[4]->Chi2Test(hAlphaComb4560low, "P");
+
+
+	c2->cd(2);
+	TH1F *hAlphaComb6075low = new TH1F(*hAlphaSim6075low[3]);
+	//	hAlphaComb7590->Add(hAlphaSim[5], 1);
+	hAlphaComb6075low->Add(hAlphaSim6075low[2], 0.1);
+	hAlphaComb6075low->Scale(0.9);
+
+	hAlphaComb6075low->SetTitle("");
+	hAlphaComb6075low->GetYaxis()->SetRangeUser(0., 400.);
+	hAlphaComb6075low->GetYaxis()->SetTitle("");
+	hAlphaComb6075low->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb6075low->GetXaxis()->SetNdivisions(506);
+	hAlphaComb6075low->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb6075low->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb6075low->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb6075low->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb6075low->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb6075low->SetLineColor(kGray+1);
+	hAlphaComb6075low->SetFillColor(kGray+1);
+
+	hAlphaComb6075low->Draw("same");
+	hAlphaExp6075low[4]->Draw("E same");
+
+	hAlphaExp6075low[4]->Chi2Test(hAlphaComb6075low, "P");
+
+	c2->cd(3);
+	TH1F *hAlphaComb7590low = new TH1F(*hAlphaSim7590low[3]);
+//	hAlphaComb7590low->Add(hAlphaSim7590low[2], 0.2);
+//	hAlphaComb7590low->Scale(0.8);
+
+	hAlphaComb7590low->SetTitle("");
+	hAlphaComb7590low->GetYaxis()->SetRangeUser(0., 300.);
+	hAlphaComb7590low->GetYaxis()->SetTitle("");
+	hAlphaComb7590low->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb7590low->GetXaxis()->SetNdivisions(506);
+	hAlphaComb7590low->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb7590low->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb7590low->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb7590low->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb7590low->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb7590low->SetLineColor(kGray+1);
+	hAlphaComb7590low->SetFillColor(kGray+1);
+
+	hAlphaComb7590low->Draw("same");
+	hAlphaExp7590low[4]->Draw("E same");
+
+	hAlphaExp7590low[4]->Chi2Test(hAlphaComb7590low, "P");
+
+	c2->cd(4);
+	TH1F *hAlphaComb90120low = new TH1F(*hAlphaSim90120low[5]);
+	hAlphaComb90120low->Add(hAlphaSim90120low[4], 0.1);
+	hAlphaComb90120low->Scale(0.8);
+
+	hAlphaComb90120low->SetTitle("");
+	hAlphaComb90120low->GetYaxis()->SetRangeUser(0., 300.);
+	hAlphaComb90120low->GetYaxis()->SetTitle("");
+	hAlphaComb90120low->GetYaxis()->SetNdivisions(30505);
+	hAlphaComb90120low->GetXaxis()->SetNdivisions(506);
+	hAlphaComb90120low->GetXaxis()->SetLabelSize(0.07);
+	hAlphaComb90120low->GetXaxis()->SetLabelOffset(0.002);
+	hAlphaComb90120low->GetXaxis()->SetTitleSize(0.07);
+	hAlphaComb90120low->GetXaxis()->SetTitleOffset(0.74);
+	hAlphaComb90120low->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+
+	hAlphaComb90120low->SetLineColor(kGray+1);
+	hAlphaComb90120low->SetFillColor(kGray+1);
+
+	hAlphaComb90120low->Draw("same");
+	hAlphaExp90120low[4]->Draw("E same");
+
+	hAlphaExp90120low[4]->Chi2Test(hAlphaComb90120low, "P");
+
+	c1->Update();
+	c2->Update();
+	c1->Update();
+
+
+	if (saveFigs) {
+		c1->SaveAs("figures/PresBestFit2531.png");
+		c2->SaveAs("figures/PresBestFit1925.png");
+	}
+
+}
diff --git a/macros/BeCorrPRC/figPresCorrGS.cxx b/macros/BeCorrPRC/figPresCorrGS.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..a3cc9e4b6f195455fff2f79f48bcb4e3e2cbff44
--- /dev/null
+++ b/macros/BeCorrPRC/figPresCorrGS.cxx
@@ -0,0 +1,134 @@
+//#include "TFile.h"
+//#include "TCanvas.h"
+//#include "TH1F.h"
+
+void figPresCorrGS(const Bool_t saveFigs = 0)
+{
+
+	TFile *fr1 = new TFile("figures/PRCepsilonT.root", "READ");
+//	TFile *fr2 = new TFile("figures/PRCepsilonT.root", "READ");
+
+	TH1F *hExp[4];
+	TH1F *hSim[4];
+	TH1F *hInput[4];
+
+	TH1F *hCurr = 0;
+
+	hCurr = (TH1F*)fr1->Get("heeT0_0_4560");
+	hExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT0_5_6075");
+	hExp[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT0_5_7590");
+	hExp[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT0_0_90120");
+	hExp[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hseT0_0_4560");
+	hSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT0_5_6075");
+	hSim[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT0_5_7590");
+	hSim[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT0_0_90120");
+	hSim[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hieT0_0_4560");
+	hInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT0_5_6075");
+	hInput[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT0_5_7590");
+	hInput[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT0_0_90120");
+	hInput[3] = new TH1F(*hCurr);
+
+
+
+	TCanvas *c1 = new TCanvas();
+	c1->ToggleEditor();
+//	c1->SetWindowSize(800, 1200);
+	c1->SetWindowSize(1000, 1000);
+	c1->Divide(2, 2, .001, .001);
+
+
+	for (Int_t i = 0; i < 4; i++) {
+		c1->cd(i+1);
+
+//		cout << i << endl;
+
+		hSim[i]->SetTitle("");
+		hSim[i]->SetLineColor(kGray+1);
+		hSim[i]->SetFillColor(kGray+1);
+		hSim[i]->GetXaxis()->SetNdivisions(405);
+		hSim[i]->GetYaxis()->SetTitle("");
+		if (i == 0 || i == 1) {
+			hSim[i]->GetYaxis()->SetNdivisions(1005);
+		}
+//		hInput[i]->SetLineWidth(2);
+		hInput[i]->SetLineColor(kRed);
+
+		hSim[i]->Draw("");
+		hExp[i]->Draw("E same");
+		hInput[i]->Draw("same");
+
+		hExp[i]->Chi2Test(hSim[i], "P");
+	}
+	c1->Update();
+
+
+	TCanvas *c2 = new TCanvas();
+	c2->ToggleEditor();
+	c2->SetWindowSize(1000, 500);
+	c2->Divide(2, 1, .001, .001);
+
+	TFile *fr3 = new TFile("figures/PRCepsilonY.root");
+	TFile *fr4 = new TFile("figures/PRCcosThetaT.root");
+
+	TH1F *hCorrExp[2];
+	TH1F *hCorrSim[2];
+	TH1F *hCorrInput[2];
+
+	hCurr = (TH1F*)fr3->Get("heeY0_0");
+	hCorrExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hecoskT0_0_7590");
+	hCorrExp[1] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr3->Get("hseY0_0");
+	hCorrSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hscoskT0_0_7590");
+	hCorrSim[1] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr3->Get("hieY0_0");
+	hCorrInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr4->Get("hicoskT0_0_7590");
+	hCorrInput[1] = new TH1F(*hCurr);
+
+	for (Int_t i = 0; i < 2; i++) {
+		c2->cd(i+1);
+
+		hCorrSim[i]->SetTitle("");
+		hCorrSim[i]->SetLineColor(kGray+1);
+		hCorrSim[i]->SetFillColor(kGray+1);
+		hCorrSim[i]->GetXaxis()->SetNdivisions(505);
+		hCorrSim[i]->GetYaxis()->SetTitle("");
+		if (i == 1) {
+			hCorrSim[i]->GetXaxis()->SetTitleSize(0.08);
+			hCorrSim[i]->GetXaxis()->SetTitleOffset(0.73);
+		}
+//		hCorrInput[i]->SetLineWidth(2);
+		hCorrInput[i]->SetLineColor(kRed);
+
+		hCorrSim[i]->Draw("");
+		hCorrExp[i]->Draw("E same");
+		hCorrInput[i]->Draw("same");
+
+		hCorrExp[i]->Chi2Test(hCorrSim[i], "P");
+	}
+	c2->Update();
+	c1->Update();
+
+	if (saveFigs) {
+		c1->SaveAs("figures/PresGSepsilonT.png");
+		c2->SaveAs("figures/PresGSotherCorrs.png");
+	}
+
+}
diff --git a/macros/BeCorrPRC/figPresCorrLeftES.cxx b/macros/BeCorrPRC/figPresCorrLeftES.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..ddf82c370c6dc288e6963d53d07c23c306132f0f
--- /dev/null
+++ b/macros/BeCorrPRC/figPresCorrLeftES.cxx
@@ -0,0 +1,290 @@
+//canvas c1: Fig. 12 in first submitted version
+//canvas c2: Fig. 13 in first submitted version
+//canvas c3: Fig. 14 in first submitted version
+//canvas c4: Fig. 15 in first submitted version
+
+using std::cout;
+using std::endl;
+
+void figPresCorrLeftES(const Bool_t saveFigs = 0)
+{
+
+//	gStyle->SetLabelFont();
+
+	TFile *fr = new TFile("figures/PRCepsilonT.root", "READ");
+
+	const Int_t canMap[6] = {1,4,3,6,2,5};
+
+	TH1F *hExp[4];
+	TH1F *hSim[4];
+	TH1F *hInput[4];
+
+	TH1F *hCurr = 0;
+
+	hCurr = (TH1F*)fr->Get("heeT3_5_4560");
+	hExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("heeT3_5_6075");
+	hExp[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("heeT3_5_7590");
+	hExp[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("heeT3_5_90120");
+	hExp[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr->Get("hseT3_5_4560");
+	hSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hseT3_5_6075");
+	hSim[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hseT3_5_7590");
+	hSim[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hseT3_5_90120");
+	hSim[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr->Get("hieT3_5_4560");
+	hInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hieT3_5_6075");
+	hInput[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hieT3_5_7590");
+	hInput[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr->Get("hieT3_5_90120");
+	hInput[3] = new TH1F(*hCurr);
+
+
+
+	TCanvas *c1 = new TCanvas();
+	c1->ToggleEditor();
+//	c1->SetWindowSize(800, 1200);
+	c1->SetWindowSize(1000, 1000);
+	c1->Divide(2, 2, .001, .001);
+
+	cout << "Fig. 12" << endl;
+
+	for (Int_t i = 0; i < 4; i++) {
+		c1->cd(i+1);
+
+//		cout << i << endl;
+
+		hSim[i]->SetTitle("");
+		hSim[i]->SetLineColor(kGray+1);
+		hSim[i]->SetFillColor(kGray+1);
+		hSim[i]->GetXaxis()->SetNdivisions(405);
+		hSim[i]->GetYaxis()->SetTitle("");
+		if (i == 0 || i == 1) {
+//			hSim[i]->GetYaxis()->SetNdivisions(1005);
+		}
+		hInput[i]->SetLineWidth(2);
+		hInput[i]->SetLineColor(kRed);
+
+		hSim[i]->Draw("");
+		hExp[i]->Draw("E same");
+		hInput[i]->Draw("same");
+
+		hExp[i]->Chi2Test(hSim[i], "P");
+	}
+	c1->Update();
+
+	TCanvas *c2 = new TCanvas();
+	c2->ToggleEditor();
+	c2->SetWindowSize(1500, 1000);
+	c2->Divide(3, 2, .001, .001);
+
+//	TFile *fr3 = new TFile("figures/PRCepsilonY.root");
+
+	TH1F *hCorrExp[6];
+	TH1F *hCorrSim[6];
+	TH1F *hCorrInput[6];
+
+	TString hName;
+
+	cout << "Fig. 13" << endl;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c2->cd(canMap[i]);
+
+		hName.Form("heeT3_%d_7590", i);
+		hCurr = (TH1F*)fr->Get(hName.Data());
+		hCorrExp[i] = new TH1F(*hCurr);
+
+		hName.Form("hseT3_%d_7590", i);
+		hCurr = (TH1F*)fr->Get(hName.Data());
+		hCorrSim[i] = new TH1F(*hCurr);
+
+		hName.Form("hieT3_%d_7590", i);
+		hCurr = (TH1F*)fr->Get(hName.Data());
+		hCorrInput[i] = new TH1F(*hCurr);
+		hCorrInput[i]->SetLineColor(kRed);
+
+		hCorrInput[i]->SetLineWidth(2);
+		hCorrSim[i]->SetTitle("");
+		hCorrSim[i]->GetXaxis()->SetNdivisions(505);
+
+		hCorrSim[i]->GetYaxis()->SetRangeUser(0., 1390.);
+		hCorrSim[i]->GetYaxis()->SetNdivisions(406);
+		hCorrSim[i]->GetYaxis()->SetTitle("");
+		hCorrSim[i]->SetLineColor(kGray+1);
+		hCorrSim[i]->SetFillColor(kGray+1);
+
+		hCorrSim[i]->Draw("");
+		hCorrExp[i]->Draw("E same");
+		hCorrInput[i]->Draw("same");
+
+		hCorrExp[i]->Chi2Test(hCorrSim[i], "P");
+	}
+	c2->Update();
+
+
+	TCanvas *c3 = new TCanvas();
+	c3->ToggleEditor();
+//	c3->SetWindowSize(800, 1227);
+	c3->SetWindowSize(1500, 1000);
+	c3->Divide(3, 2, .001, .001);
+
+	TFile *fr2 = new TFile("figures/PRCthetaA.root");
+
+	TH1F *hAlphaExp[6];
+	TH1F *hAlphaSim[6];
+	TH1F *hAlphaInput[6];
+	TH1F *hAlphaInputCorr[6];
+
+
+	cout << "Fig. 14" << endl;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c3->cd(canMap[i]);
+
+		hName.Form("hethetaAT3_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim[i] = new TH1F(*hCurr);
+
+		hName.Form("hithetaAT3_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaInput[i] = new TH1F(*hCurr);
+
+		hAlphaInputCorr[i] = new TH1F(*hAlphaInput[i]);
+		Double_t binContent = 0;
+		for (Int_t binN = 1; binN <= 15; binN++) {
+			binContent = hAlphaInputCorr[i]->GetBinContent(binN);
+			hAlphaInputCorr[i]->SetBinContent(binN, binContent + hAlphaInputCorr[i]->GetBinContent(31-binN));
+			hAlphaInputCorr[i]->SetBinContent(31-binN, hAlphaInputCorr[i]->GetBinContent(binN));
+		}
+
+		hAlphaInput[i]->SetLineWidth(2);
+		hAlphaSim[i]->SetTitle("");
+		hAlphaSim[i]->GetXaxis()->SetNdivisions(506);
+		hAlphaSim[i]->GetXaxis()->SetLabelSize(0.07);
+		hAlphaSim[i]->GetXaxis()->SetLabelOffset(0.002);
+		hAlphaSim[i]->GetXaxis()->SetTitleSize(0.08);
+		hAlphaSim[i]->GetXaxis()->SetTitleOffset(0.74);
+
+		hAlphaSim[i]->GetYaxis()->SetRangeUser(0., 690.);
+		hAlphaSim[i]->GetYaxis()->SetTitle("");
+
+		hAlphaSim[i]->SetLineColor(kGray+1);
+		hAlphaSim[i]->SetFillColor(kGray+1);
+
+		hAlphaSim[i]->Draw("");
+		hAlphaSim[i]->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+		hAlphaExp[i]->Draw("E same");
+//		hAlphaInput[i]->Draw("same");
+		hAlphaInputCorr[i]->SetLineColor(kRed);
+		hAlphaInputCorr[i]->SetLineWidth(2);
+		hAlphaInputCorr[i]->Scale(0.5);
+		hAlphaInputCorr[i]->Smooth();
+		hAlphaInputCorr[i]->Draw("same");
+
+//		cout << "Chi2: " << hAlphaSim[i]->Chi2Test()endl;
+		hAlphaExp[i]->Chi2Test(hAlphaSim[i], "P");
+//		hAlphaSim[i]->Chi2Test(hAlphaInputCorr[i], "P");
+
+	}
+	c3->Update();
+
+
+
+
+	TCanvas *c4 = new TCanvas();
+	c4->ToggleEditor();
+//	c4->SetWindowSize(800, 1227);
+//	c4->Divide(2, 3, .001, .001);
+	c4->SetWindowSize(1500, 1000);
+	c4->Divide(3, 2, .001, .001);
+
+//	TFile *fr2 = new TFile("figures/PRCthetaA.root");
+
+	TH1F *hAlphaExp90120[6];
+	TH1F *hAlphaSim90120[6];
+	TH1F *hAlphaInput90120[6];
+	TH1F *hAlphaInput90120Corr[6];
+
+
+	cout << "Fig. 15" << endl;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c4->cd(canMap[i]);
+
+		hName.Form("hethetaAT3_%d_90120", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp90120[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT3_%d_90120", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim90120[i] = new TH1F(*hCurr);
+
+		hName.Form("hithetaAT3_%d_90120", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaInput90120[i] = new TH1F(*hCurr);
+
+		hAlphaInput90120Corr[i] = new TH1F(*hAlphaInput90120[i]);
+		Double_t binContent = 0;
+		for (Int_t binN = 1; binN <= 15; binN++) {
+			binContent = hAlphaInput90120Corr[i]->GetBinContent(binN);
+			hAlphaInput90120Corr[i]->SetBinContent(binN, binContent + hAlphaInput90120Corr[i]->GetBinContent(31-binN));
+			hAlphaInput90120Corr[i]->SetBinContent(31-binN, hAlphaInput90120Corr[i]->GetBinContent(binN));
+		}
+
+		hAlphaInput90120[i]->SetLineWidth(2);
+		hAlphaSim90120[i]->SetTitle("");
+		hAlphaSim90120[i]->GetXaxis()->SetNdivisions(506);
+		hAlphaSim90120[i]->GetXaxis()->SetLabelSize(0.07);
+		hAlphaSim90120[i]->GetXaxis()->SetLabelOffset(0.002);
+		hAlphaSim90120[i]->GetXaxis()->SetTitleSize(0.08);
+		hAlphaSim90120[i]->GetXaxis()->SetTitleOffset(0.74);
+
+		hAlphaSim90120[i]->GetYaxis()->SetRangeUser(0., 650.);
+		hAlphaSim90120[i]->GetYaxis()->SetTitle("");
+
+		hAlphaSim90120[i]->SetLineColor(kGray+1);
+		hAlphaSim90120[i]->SetFillColor(kGray+1);
+
+		hAlphaSim90120[i]->Draw("");
+		hAlphaSim90120[i]->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+		hAlphaExp90120[i]->Draw("E same");
+//		hAlphaInput90120[i]->Draw("same");
+		hAlphaInput90120Corr[i]->SetLineColor(kRed);
+		hAlphaInput90120Corr[i]->SetLineWidth(2);
+		hAlphaInput90120Corr[i]->Scale(0.5);
+		hAlphaInput90120Corr[i]->Smooth();
+		hAlphaInput90120Corr[i]->Draw("same");
+
+		hAlphaExp90120[i]->Chi2Test(hAlphaSim90120[i], "P");
+
+	}
+	c4->Update();
+	c1->Update();
+
+	if (saveFigs) {
+		c1->SaveAs("figures/PresESepsilonTleft.png");
+		c2->SaveAs("figures/PresESepsilonT7590left.png");
+		c3->SaveAs("figures/PresESThetaAleft7590.png");
+		c4->SaveAs("figures/PresESThetaAleft90120.png");
+
+//		c1->SaveAs("figures/PresESepsilonTleft.gif");
+//		c2->SaveAs("figures/PresESepsilonT7590left.gif");
+//		c3->SaveAs("figures/PresESThetaAleft7590.gif");
+//		c4->SaveAs("figures/PresESThetaAleft90120.gif");
+	}
+
+}
diff --git a/macros/BeCorrPRC/figPresCorrRightES.cxx b/macros/BeCorrPRC/figPresCorrRightES.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..0e549f86064d0e2963237331d08fe6973eb498b0
--- /dev/null
+++ b/macros/BeCorrPRC/figPresCorrRightES.cxx
@@ -0,0 +1,88 @@
+//canvas c3: Fig. 16 in first submitted version
+
+void figPresCorrRightES(const Bool_t saveFigs = 0)
+{
+
+	const Int_t canMap[6] = {1,4,3,6,2,5};
+
+
+	TCanvas *c1 = new TCanvas();
+	c1->Update();
+
+	TH1F *hCurr = 0;
+
+	TCanvas *c3 = new TCanvas();
+	c3->ToggleEditor();
+	c3->SetWindowSize(1500, 1000);
+	c3->Divide(3, 2, .001, .001);
+
+	TFile *fr2 = new TFile("figures/PRCthetaARightES.root");
+
+	TH1F *hAlphaExp[6];
+	TH1F *hAlphaSim[6];
+	TH1F *hAlphaInput[6];
+	TH1F *hAlphaInputCorr[6];
+
+	TString hName;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c3->cd(canMap[i]);
+
+		hName.Form("hethetaAT4_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT4_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim[i] = new TH1F(*hCurr);
+
+		hName.Form("hithetaAT4_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaInput[i] = new TH1F(*hCurr);
+
+		hAlphaInputCorr[i] = new TH1F(*hAlphaInput[i]);
+		Double_t binContent = 0;
+		for (Int_t binN = 1; binN <= 15; binN++) {
+			binContent = hAlphaInputCorr[i]->GetBinContent(binN);
+			hAlphaInputCorr[i]->SetBinContent(binN, binContent + hAlphaInputCorr[i]->GetBinContent(31-binN));
+			hAlphaInputCorr[i]->SetBinContent(31-binN, hAlphaInputCorr[i]->GetBinContent(binN));
+		}
+
+		hAlphaInput[i]->SetLineWidth(2);
+		hAlphaSim[i]->SetTitle("");
+		hAlphaSim[i]->GetXaxis()->SetNdivisions(506);
+		hAlphaSim[i]->GetXaxis()->SetLabelSize(0.07);
+		hAlphaSim[i]->GetXaxis()->SetLabelOffset(0.002);
+		hAlphaSim[i]->GetXaxis()->SetTitleSize(0.08);
+		hAlphaSim[i]->GetXaxis()->SetTitleOffset(0.74);
+
+		hAlphaSim[i]->GetYaxis()->SetRangeUser(0., 690.);
+		hAlphaSim[i]->GetYaxis()->SetTitle("");
+
+		hAlphaSim[i]->SetFillColor(kGray+1);
+		hAlphaSim[i]->SetLineColor(kGray+1);
+
+		hAlphaSim[i]->Draw("");
+		hAlphaSim[i]->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+		hAlphaExp[i]->Draw("E same");
+//		hAlphaInput[i]->Draw("same");
+		hAlphaInputCorr[i]->SetLineColor(kRed);
+		hAlphaInputCorr[i]->SetLineWidth(2);
+		hAlphaInputCorr[i]->Scale(0.5);
+		hAlphaInputCorr[i]->Smooth();
+		hAlphaInputCorr[i]->Draw("same");
+
+		hAlphaExp[i]->Chi2Test(hAlphaSim[i], "P");
+
+	}
+	c3->Update();
+	c1->Update();
+
+	delete c1;
+
+	if (saveFigs) {
+		c3->SaveAs("figures/ESThetaAright7590.png");
+	}
+
+
+}
diff --git a/macros/BeCorrPRC/figPresCorrRightGS.cxx b/macros/BeCorrPRC/figPresCorrRightGS.cxx
new file mode 100755
index 0000000000000000000000000000000000000000..568e63029b3e30a12b9fc3572ee6fe9e98e3a76f
--- /dev/null
+++ b/macros/BeCorrPRC/figPresCorrRightGS.cxx
@@ -0,0 +1,203 @@
+//canvas c1: Fig. 9 in first submitted version
+//canvas c2: Fig. 10 in first submitted version
+//canvas c3: Fig. 11 in first submitted version
+
+//#include "TFile.h"
+//#include "TCanvas.h"
+//#include "TH1F.h"
+
+using std::cout;
+using std::endl;
+
+void figPresCorrRightGS(const Bool_t saveFigs = 0)
+{
+
+	const Int_t canMap[6] = {1,4,3,6,2,5};
+
+	TFile *fr1 = new TFile("figures/PRCepsilonT.root", "READ");
+//	TFile *fr2 = new TFile("figures/PRCepsilonT.root", "READ");
+
+	TH1F *hExp[4];
+	TH1F *hSim[4];
+	TH1F *hInput[4];
+
+	TH1F *hCurr = 0;
+
+	hCurr = (TH1F*)fr1->Get("heeT1_2_4560");
+	hExp[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT1_2_6075");
+	hExp[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT1_2_7590");
+	hExp[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("heeT1_2_90120");
+	hExp[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hseT1_2_4560");
+	hSim[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT1_2_6075");
+	hSim[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT1_2_7590");
+	hSim[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hseT1_2_90120");
+	hSim[3] = new TH1F(*hCurr);
+
+	hCurr = (TH1F*)fr1->Get("hieT1_2_4560");
+	hInput[0] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT1_2_6075");
+	hInput[1] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT1_2_7590");
+	hInput[2] = new TH1F(*hCurr);
+	hCurr = (TH1F*)fr1->Get("hieT1_2_90120");
+	hInput[3] = new TH1F(*hCurr);
+
+
+
+	TCanvas *c1 = new TCanvas();
+	c1->ToggleEditor();
+//	c1->SetWindowSize(800, 1200);
+	c1->SetWindowSize(1000, 1000);
+	c1->Divide(2, 2, .001, .001);
+
+	cout << "Fig. 9" << endl;
+
+	for (Int_t i = 0; i < 4; i++) {
+		c1->cd(i+1);
+
+//		cout << i << endl;
+
+		hSim[i]->SetTitle("");
+		hSim[i]->SetLineColor(kGray+1);
+		hSim[i]->SetFillColor(kGray+1);
+		hSim[i]->GetXaxis()->SetNdivisions(405);
+		hSim[i]->GetYaxis()->SetTitle("");
+		if (i == 0 || i == 1) {
+//			hSim[i]->GetYaxis()->SetNdivisions(1005);
+		}
+		hInput[i]->SetLineColor(kRed);
+
+		hSim[i]->Draw("");
+		hExp[i]->Draw("E same");
+		hInput[i]->Draw("same");
+		hExp[i]->Chi2Test(hSim[i], "P");
+	}
+	c1->Update();
+
+	TCanvas *c2 = new TCanvas();
+	c2->ToggleEditor();
+	c2->SetWindowSize(1500, 1000);
+	c2->Divide(3, 2, .001, .001);
+
+	TH1F *hCorrExp[6];
+	TH1F *hCorrSim[6];
+	TH1F *hCorrInput[6];
+
+	TString hName;
+
+	cout << "Fig. 10" << endl;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c2->cd(canMap[i]);
+
+		hName.Form("heeT1_%d_7590", i);
+		hCurr = (TH1F*)fr1->Get(hName.Data());
+		hCorrExp[i] = new TH1F(*hCurr);
+
+		hName.Form("hseT1_%d_7590", i);
+		hCurr = (TH1F*)fr1->Get(hName.Data());
+		hCorrSim[i] = new TH1F(*hCurr);
+
+		hName.Form("hieT1_%d_7590", i);
+		hCurr = (TH1F*)fr1->Get(hName.Data());
+		hCorrInput[i] = new TH1F(*hCurr);
+
+		hCorrInput[i]->SetLineColor(kRed);
+		hCorrSim[i]->SetTitle("");
+		hCorrSim[i]->GetXaxis()->SetNdivisions(505);
+
+		hCorrSim[i]->GetYaxis()->SetRangeUser(0., 850.);
+		hCorrSim[i]->GetYaxis()->SetNdivisions(505);
+		hCorrSim[i]->GetYaxis()->SetTitle("");
+
+		hCorrSim[i]->SetLineColor(kGray+1);
+		hCorrSim[i]->SetFillColor(kGray+1);
+
+		hCorrSim[i]->Draw("");
+		hCorrExp[i]->Draw("E same");
+		hCorrInput[i]->Draw("same");
+
+		hCorrExp[i]->Chi2Test(hCorrSim[i], "P");
+	}
+	c2->Update();
+
+	TCanvas *c3 = new TCanvas();
+	c3->ToggleEditor();
+	c3->SetWindowSize(1500, 1000);
+	c3->Divide(3, 2, .001, .001);
+
+	TFile *fr2 = new TFile("figures/PRCthetaA.root");
+
+	TH1F *hAlphaExp[6];
+	TH1F *hAlphaSim[6];
+	TH1F *hAlphaInput[6];
+	TH1F *hAlphaInputCorr[6];
+
+	cout << "Fig. 11" << endl;
+
+	for (Int_t i = 0; i < 6; i++) {
+		c3->cd(canMap[i]);
+
+		hName.Form("hethetaAT1_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaExp[i] = new TH1F(*hCurr);
+
+		hName.Form("hsthetaAT1_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaSim[i] = new TH1F(*hCurr);
+
+		hName.Form("hithetaAT1_%d_7590", i);
+		hCurr = (TH1F*)fr2->Get(hName.Data());
+		hAlphaInput[i] = new TH1F(*hCurr);
+
+		hAlphaInputCorr[i] = new TH1F(*hAlphaInput[i]);
+		Double_t binContent = 0;
+		for (Int_t binN = 1; binN <= 15; binN++) {
+			binContent = hAlphaInputCorr[i]->GetBinContent(binN);
+			hAlphaInputCorr[i]->SetBinContent(binN, binContent + hAlphaInputCorr[i]->GetBinContent(31-binN));
+			hAlphaInputCorr[i]->SetBinContent(31-binN, hAlphaInputCorr[i]->GetBinContent(binN));
+		}
+
+		hAlphaInput[i]->SetLineWidth(2);
+		hAlphaSim[i]->SetTitle("");
+		hAlphaSim[i]->GetXaxis()->SetNdivisions(506);
+		hAlphaSim[i]->GetXaxis()->SetLabelSize(0.07);
+		hAlphaSim[i]->GetXaxis()->SetLabelOffset(0.002);
+		hAlphaSim[i]->GetXaxis()->SetTitleSize(0.08);
+		hAlphaSim[i]->GetXaxis()->SetTitleOffset(0.74);
+
+		hAlphaSim[i]->GetYaxis()->SetRangeUser(0., 290.);
+		hAlphaSim[i]->GetYaxis()->SetTitle("");
+
+		hAlphaSim[i]->SetLineColor(kGray+1);
+		hAlphaSim[i]->SetFillColor(kGray+1);
+
+		hAlphaSim[i]->Draw("");
+		hAlphaSim[i]->GetXaxis()->SetTitle("#it{#theta}_{#alpha} (rad)");
+		hAlphaExp[i]->Draw("E same");
+//		hAlphaInput[i]->Draw("same");
+		hAlphaInputCorr[i]->SetLineColor(kRed);
+		hAlphaInputCorr[i]->SetLineWidth(2);
+		hAlphaInputCorr[i]->Scale(0.5);
+		hAlphaInputCorr[i]->Smooth();
+		hAlphaInputCorr[i]->Draw("same");
+
+		hAlphaExp[i]->Chi2Test(hAlphaSim[i], "P");
+
+	}
+	c3->Update();
+
+	if (saveFigs) {
+		c1->SaveAs("figures/PresGSepsilonTright.png");
+		c2->SaveAs("figures/PresGSepsilonTright7590.png");
+		c3->SaveAs("figures/PresGSThetaAright7590.png");
+	}
+}
diff --git a/macros/BeCorrPRC/figSpectraIntegral.cxx b/macros/BeCorrPRC/figSpectraIntegral.cxx
new file mode 100644
index 0000000000000000000000000000000000000000..f25f2d3412c92a003fffe58d68794920d2738d5a
--- /dev/null
+++ b/macros/BeCorrPRC/figSpectraIntegral.cxx
@@ -0,0 +1,148 @@
+#include "TCanvas.h"
+
+void figSpectraIntegral()
+{
+	TStopwatch sw;
+	sw.Start();
+
+	const Int_t lowfile = 0;
+	const Int_t upfile = 40;
+
+	const Int_t lowrfile = 0;
+	const Int_t uprfile = 0;
+
+//	TChain *che = BeWork::OpenChain("/data2/be/rootdata/correlations/Be.",
+	TChain *che = BeWork::OpenChain("../../../be/rootdata/correlations/v5_exp/Be.",
+			lowfile, upfile, "beonly");										//original file
+	Info("figSpectraIntegral", "%lld events in chain \"%s\" containing experimental data",
+			che->GetEntries(), che->GetName());
+//
+//	TChain *cheRaw = BeWork::OpenChain("/data2/be/rootdata/correlations/Exp.", lowfile, upfile, "be", kBlack);
+//
+//	TChain *chraw = BeWork::OpenChain("/data2/be/rootdata/correlations/Exp.",
+//			lowrfile, uprfile, "be");										//original file
+//	Info("statesRatioFitting.cxx", "%d events in chain \"%s\" containing experimental data",
+//			chraw->GetEntries(), chraw->GetName());
+
+	const Bool_t savePictures = 0;
+	TString pictureName;
+	TString picturePath = "figures/";
+//	TString pictFormat = ".jpg";
+	TString pictFormat = ".eps";
+//	TString pictFormat = ".gif";
+
+	Int_t wtopx = 0, wtopy = 0;
+	UInt_t ww = 0, wh = 0;
+
+
+//	if (kspectrum1D) {
+	TCanvas *cSpectrum = new TCanvas("spectrumIntegral", "spectrum up to 16 MeV");
+	cSpectrum->Divide(1,2, 0., 0.);
+//	cSpectrum->SetBottomMargin(0.15);
+
+	cSpectrum->cd(1);
+	che->Draw("fBeIM>>hBe(200,0,7.9)", "fBeIM<7.9");
+	TH1F *hBe = (TH1F*)gPad->FindObject("hBe");
+	hBe->SetTitle("");
+	hBe->SetLineColor(kGray+2);
+	hBe->SetFillColor(kGray+2);
+	hBe->SetXTitle("E_{T} (MeV)");
+	hBe->GetXaxis()->SetTitleOffset(0.71);
+	hBe->GetXaxis()->SetLabelSize(0.05);
+	hBe->GetXaxis()->SetTitleSize(0.09);
+	hBe->GetXaxis()->CenterTitle();
+	hBe->SetYTitle("Counts");
+	hBe->GetYaxis()->SetTitleOffset(0.75);
+	hBe->GetYaxis()->SetLabelSize(0.05);
+	hBe->GetYaxis()->SetTitleSize(0.09);
+	hBe->GetYaxis()->CenterTitle();
+	hBe->GetYaxis()->SetRangeUser(0,19000);
+	hBe->GetXaxis()->SetNdivisions(516);
+
+//	TLine *lineE[4];
+//
+//	lineE[0] = new TLine(3, 45, 3, 120);
+//	lineE[1] = new TLine(0, 45, 0, 120);
+//	lineE[2] = new TLine(0, 120, 3, 120);
+//	lineE[3] = new TLine(0, 45, 3, 45);
+
+
+	cSpectrum->cd(2);
+	che->Draw("fBeThetaCM1*TMath::RadToDeg():fBeIM>>hBe2D(100,0,7.9, 100,0,179)", "fBeIM<7.9", "cont");
+	TH2F *hBe2D = (TH2F*)gPad->FindObject("hBe2D");
+	hBe2D->SetTitle("");
+	hBe2D->SetLineColor(kGray+2);
+	hBe2D->SetFillColor(kGray+2);
+	hBe2D->SetXTitle("E_{T} (MeV)");
+	hBe2D->GetXaxis()->SetTitleOffset(0.71);
+	hBe2D->GetXaxis()->SetLabelSize(0.05);
+	hBe2D->GetXaxis()->SetTitleSize(0.09);
+	hBe2D->GetXaxis()->CenterTitle();
+	hBe2D->SetYTitle("\\theta_{Be} (deg.)");
+	hBe2D->GetYaxis()->SetTitleOffset(0.60);
+	hBe2D->GetYaxis()->SetLabelSize(0.05);
+	hBe2D->GetYaxis()->SetTitleSize(0.09);
+	hBe2D->GetYaxis()->CenterTitle();
+	hBe2D->GetXaxis()->SetNdivisions(516);
+
+	TLine *line[4];
+
+	line[0] = new TLine(3, 45, 3, 120);
+	line[1] = new TLine(0, 45, 0, 120);
+	line[2] = new TLine(0, 120, 3, 120);
+	line[3] = new TLine(0, 45, 3, 45);
+
+	for (Int_t i = 0; i < 4; i++) {
+		line[i]->SetLineWidth(3);
+		line[i]->SetLineColor(kRed);
+		line[i]->Draw("same");
+	}
+
+	if (savePictures) {
+		pictureName.Form("%s%s%s", picturePath.Data(), cSpectrum->GetName(), pictFormat.Data());
+		cSpectrum->SaveAs(pictureName.Data());
+		cSpectrum->cd();
+		cSpectrum->Close();
+	}
+
+
+
+
+
+	//	}
+
+	//	if (kThetaLab) {
+	//		TCanvas *cThetas = new TCanvas("thetaLab", "really measured thetas in lab frame");
+	//		che->Draw("fAThetaLab*TMath::RadToDeg()>>hATheta(500,0,27)", "");
+	//		che->Draw("fP1ThetaLab*TMath::RadToDeg()>>hAP1(100,0,27)", "", "same goff");
+	//		che->Draw("fP2ThetaLab*TMath::RadToDeg()>>hAP2(100,0,27)", "", "same goff");
+	//		TH1F *hATheta = (TH1F*)gPad->FindObject("hATheta");
+	//		hATheta->SetTitle("");
+	////		hATheta->SetLineColor(kGray+2);
+	////		hATheta->SetFillColor(kGray+2);
+//		hATheta->SetXTitle("\\theta_{lab} [\\circ]");
+//		hATheta->GetXaxis()->SetTitleOffset(0.95);
+//		hATheta->GetXaxis()->CenterTitle();
+//		hATheta->SetYTitle("Events");
+//		hATheta->GetYaxis()->SetTitleOffset(1.1);
+//		hATheta->GetYaxis()->CenterTitle();
+//		hATheta->SetLineColor(kRed);
+//		hATheta->SetFillColor(kRed);
+//		hATheta->SetFillStyle(3005);
+//
+//		TH1F hP(*hAP2);
+//		hP.Add(hAP1);
+//		hP.SetLineColor(kBlue);
+//		hP.SetFillColor(kBlue);
+//		hP.SetFillStyle(3004);
+//		hP.Draw("same");
+//
+//		if (savePictures) {
+//			pictureName.Form("%s%s%s", picturePath.Data(), cThetas->GetName(), pictFormat.Data());
+//			cThetas->SaveAs(pictureName.Data());
+//			cThetas->cd();
+//			cThetas->Close();
+//		}
+//	}
+
+}
diff --git a/macros/BeCorrPRC/rootlogonPresentation.C b/macros/BeCorrPRC/rootlogonPresentation.C
new file mode 100644
index 0000000000000000000000000000000000000000..3702de8d7045ff2ed3eaef3544ef59af6403e9de
--- /dev/null
+++ b/macros/BeCorrPRC/rootlogonPresentation.C
@@ -0,0 +1,110 @@
+//#include "TStyle.h"
+
+void rootlogonPresentation()
+{
+
+	Info("rootlogonPresentation.C", "This settings are used to make a nice pictures for presentation");
+
+	gROOT->ProcessLine(".x load_libs.cxx");
+
+	//gStyle choice
+//	gROOT->SetStyle("Default");
+//	gROOT->SetStyle("Classic");
+	gROOT->SetStyle("Plain");
+//	gROOT->SetStyle("Bold");
+//	gROOT->SetStyle("Video");
+
+//	return;
+
+//	gStyle->SetOptStat("n");
+//	gStyle->SetOptStat("ne");
+	gStyle->SetOptStat("");
+
+	//canvas size and position
+//	gStyle->SetCanvasDefX(85);
+//	gStyle->SetCanvasDefY(224);
+
+//	gStyle->SetCanvasDefW(649);
+//	gStyle->SetCanvasDefH(727);
+
+	//marker attributes
+	gStyle->SetMarkerStyle(20);
+	gStyle->SetMarkerSize(0.1);
+	gStyle->SetMarkerColor(1);
+
+//	gStyle->SetLabelOffset(0.55, "X");
+
+//	gStyle->SetLabelSize(0.06, "X");
+	gStyle->SetLabelSize(0.07, "X");
+	gStyle->SetTitleSize(0.11, "X");
+	gStyle->SetTitleOffset(0.55, "X");
+	gStyle->SetTitleXSize(0.11);
+
+//	gStyle->SetLabelSize(0.06, "Y");
+	gStyle->SetLabelSize(0.07, "Y");
+	gStyle->SetTitleSize(0.11, "Y");
+	gStyle->SetTitleOffset(0.61, "Y");
+	gStyle->SetTitleYSize(0.11);
+
+
+
+	//line attributes
+//	gStyle->SetHistLineWidth(1);
+//	gStyle->SetHistLineWidth(2);
+	gStyle->SetHistLineWidth(3);
+
+	gStyle->SetTitleBorderSize(0);
+
+	gStyle->SetFrameLineWidth(2);
+
+
+	gStyle->SetTitleFont(62, "X");
+	gStyle->SetTitleFont(62, "Y");
+	gStyle->SetLabelFont(62, "X");
+	gStyle->SetLabelFont(62, "Y");
+
+	gStyle->SetPadTopMargin(0.04);
+	gStyle->SetPadBottomMargin(0.13);
+	gStyle->SetPadLeftMargin(0.13);
+	gStyle->SetPadRightMargin(0.03);
+//	gStyle->SetTitleOffset(1.3, "Y");
+//	gStyle->SetLabelSize(0.05, "X");
+//	gStyle->SetLabelSize(0.05, "Y");
+//	gStyle->SetTitleSize(0.06, "X");
+//	gStyle->SetTitleSize(0.06, "Y");
+
+//
+
+	gROOT->ForceStyle();
+
+	return;
+
+	//parameters used for pictures in EXON paper
+//	gStyle->SetStatBorderSize(0);
+//	gStyle->SetTitleBorderSize(0);
+//	gStyle->SetTitleAlign(12);
+//	gStyle->SetTitleX(0.83);
+//	gStyle->SetTitleY(0.85);
+	gStyle->SetTitleX(0.83);	//for Fig2
+	gStyle->SetTitleY(0.9);	//for Fig2
+	gStyle->SetTitleFontSize(0.08);
+//	gStyle->SetTitleFillColor(0);
+//	gStyle->SetLineWidth(1);
+//	gStyle->SetHistLineWidth(1);
+	gStyle->SetHistLineWidth(3);
+	gStyle->SetPadTopMargin(0.05);
+	gStyle->SetPadBottomMargin(0.13);
+	gStyle->SetPadLeftMargin(0.14);
+	gStyle->SetPadRightMargin(0.05);
+	gStyle->SetTitleOffset(1.3, "Y");
+	gStyle->SetLabelSize(0.05, "X");
+	gStyle->SetLabelSize(0.05, "Y");
+	gStyle->SetTitleSize(0.06, "X");
+	gStyle->SetTitleSize(0.06, "Y");
+
+
+//	gStyle->SetAxisColor(1);
+	gStyle->SetLabelColor(1, "X");
+	gStyle->SetLabelColor(1, "Y");
+
+}