From 3d649d46cb5c239a9a843ee92a2eb375b8bfd35e Mon Sep 17 00:00:00 2001 From: Vratislav Chudoba Date: Fri, 1 Jun 2018 18:45:03 +0300 Subject: [PATCH] New scripts for presentation use. --- macros/BeCorrPRC/figPresBestFitLeftES.cxx | 360 ++++++++++++++++++++++ macros/BeCorrPRC/figPresCorrGS.cxx | 134 ++++++++ macros/BeCorrPRC/figPresCorrLeftES.cxx | 290 +++++++++++++++++ macros/BeCorrPRC/figPresCorrRightES.cxx | 88 ++++++ macros/BeCorrPRC/figPresCorrRightGS.cxx | 203 ++++++++++++ macros/BeCorrPRC/figSpectraIntegral.cxx | 148 +++++++++ macros/BeCorrPRC/rootlogonPresentation.C | 110 +++++++ 7 files changed, 1333 insertions(+) create mode 100755 macros/BeCorrPRC/figPresBestFitLeftES.cxx create mode 100755 macros/BeCorrPRC/figPresCorrGS.cxx create mode 100755 macros/BeCorrPRC/figPresCorrLeftES.cxx create mode 100755 macros/BeCorrPRC/figPresCorrRightES.cxx create mode 100755 macros/BeCorrPRC/figPresCorrRightGS.cxx create mode 100644 macros/BeCorrPRC/figSpectraIntegral.cxx create mode 100644 macros/BeCorrPRC/rootlogonPresentation.C diff --git a/macros/BeCorrPRC/figPresBestFitLeftES.cxx b/macros/BeCorrPRC/figPresBestFitLeftES.cxx new file mode 100755 index 0000000..e44161a --- /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 0000000..a3cc9e4 --- /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 0000000..ddf82c3 --- /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 0000000..0e549f8 --- /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 0000000..568e630 --- /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 0000000..f25f2d3 --- /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 0000000..3702de8 --- /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"); + +} -- 2.18.1