From 5cf020599a96934fb6afade2c0fae98f51078f20 Mon Sep 17 00:00:00 2001 From: Vratislav Chudoba Date: Fri, 1 Jun 2018 18:44:05 +0300 Subject: [PATCH] New scripts for PRC figures. --- macros/BeCorrPRC/figBestFitLeftES.cxx | 342 ++++++++++++++++++++++++++ macros/BeCorrPRC/figCorrRightES.cxx | 242 ++++++++++++++++++ macros/BeCorrPRC/figKinematics.cxx | 90 +++++++ 3 files changed, 674 insertions(+) create mode 100755 macros/BeCorrPRC/figBestFitLeftES.cxx create mode 100755 macros/BeCorrPRC/figCorrRightES.cxx create mode 100644 macros/BeCorrPRC/figKinematics.cxx diff --git a/macros/BeCorrPRC/figBestFitLeftES.cxx b/macros/BeCorrPRC/figBestFitLeftES.cxx new file mode 100755 index 0000000..ac18642 --- /dev/null +++ b/macros/BeCorrPRC/figBestFitLeftES.cxx @@ -0,0 +1,342 @@ +//canvas c1: new Figure for revised version + +using std::cout; +using std::endl; + +void figBestFitLeftES(Bool_t saveFigs = 0) +{ + + 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", 800, 800); + 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->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->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->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->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", 800, 800); + 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->GetYaxis()->SetRangeUser(0., 650.); +// hAlphaComb4560low->GetYaxis()->SetTitle(""); + + hAlphaComb4560low->Draw(""); + + + + + + + + + +// hAlphaComb4560low->Draw("same"); + 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->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->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->Draw("same"); + hAlphaExp90120low[4]->Draw("E same"); + + hAlphaExp90120low[4]->Chi2Test(hAlphaComb90120low, "P"); + + c1->Update(); + + + if (saveFigs) { + } + +} diff --git a/macros/BeCorrPRC/figCorrRightES.cxx b/macros/BeCorrPRC/figCorrRightES.cxx new file mode 100755 index 0000000..34fd3d6 --- /dev/null +++ b/macros/BeCorrPRC/figCorrRightES.cxx @@ -0,0 +1,242 @@ +//canvas c3: Fig. 16 in first submitted version + +void figCorrRightES(Bool_t saveFigs = 0) +{ + + /*TFile *fr = new TFile("figures/PRCepsilonT.root", "READ"); + + 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(800, 800); + 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); + + hSim[i]->Draw(""); + hExp[i]->Draw("E same"); + hInput[i]->Draw("same"); + } + + c1->SaveAs("figures/ESepsilonTleft.eps"); + +// return; + + TCanvas *c2 = new TCanvas(); + c2->ToggleEditor(); + c2->SetWindowSize(800, 1227); + c2->Divide(2, 3, .001, .001); + +// TFile *fr3 = new TFile("figures/PRCepsilonY.root"); + + TH1F *hCorrExp[6]; + TH1F *hCorrSim[6]; + TH1F *hCorrInput[6]; + + TString hName; + + for (Int_t i = 0; i < 6; i++) { + c2->cd(i+1); + + 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]->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]->Draw(""); + hCorrExp[i]->Draw("E same"); + hCorrInput[i]->Draw("same"); + } + + c2->SaveAs("figures/ESepsilonT7590left.eps");*/ + +// return; + + TCanvas *c1 = new TCanvas(); + c1->Update(); + + TH1F *hCurr = 0; + + TCanvas *c3 = new TCanvas(); + c3->ToggleEditor(); + c3->SetWindowSize(800, 1227); + c3->Divide(2, 3, .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(i+1); + + 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]->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.eps"); + } + + + /*TCanvas *c4 = new TCanvas(); + c4->ToggleEditor(); + c4->SetWindowSize(800, 1227); + c4->Divide(2, 3, .001, .001); + +// TFile *fr2 = new TFile("figures/PRCthetaA.root"); + + TH1F *hAlphaExp90120[6]; + TH1F *hAlphaSim90120[6]; + TH1F *hAlphaInput90120[6]; + + TString hName; + + for (Int_t i = 0; i < 6; i++) { + c4->cd(i+1); + + 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); + + 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]->Draw(""); + hAlphaExp90120[i]->Draw("E same"); + hAlphaInput90120[i]->Draw("same"); + } + + + c4->SaveAs("figures/ESThetaAleft90120.eps");*/ + +} diff --git a/macros/BeCorrPRC/figKinematics.cxx b/macros/BeCorrPRC/figKinematics.cxx new file mode 100644 index 0000000..d8c9b03 --- /dev/null +++ b/macros/BeCorrPRC/figKinematics.cxx @@ -0,0 +1,90 @@ +void figKinematics() +{ + +// gSystem->Load("/home/vratik/workspace/daf/trunk/libUtilities.so"); +// gSystem->Load("/home/vratik/workspace/daf/trunk/libEvent.so"); + + + TStopwatch sw; + sw.Start(); + +// const Bool_t savePictures = 1; + const Bool_t savePictures = 0; + TString pictureName; + TString picturePath = "figures/"; +// TString pictFormat = ".jpg"; +// TString pictFormat = ".eps"; + TString pictFormat = ".gif"; + + const Bool_t kKinematics = 1; + +// if (kKinematics) { + TCanvas *cKin = new TCanvas("kinematicsBe", "kinematics of (p,n) reaction"); + cKin->ToggleEditor(); + cKin->SetWindowSize(1500, 1000); + + +// TFile fr("../simulace/kin.root"); + TFile *fr = new TFile("./kin.root"); + TTree *tr = (TTree*)fr->Get("kin"); + + + tr->SetMarkerStyle(20); + tr->SetMarkerSize(0.4); + + TH2F *hKin = new TH2F("hKin", "", 1000, 0., 30., 1000, 4., 56.); +// TCut cP = "p.fImpulse.Theta()*TMath::RadToDeg()>100"; +// TCut cP = "p.fImpulse.Theta()*TMath::RadToDeg()>0"; + +// tr->Draw("(p1.fImpulse.T()-p1.fMass)/p1.fA:p1.fImpulse.Theta()*TMath::RadToDeg()>>hKin", "", "col"); +// hKin->Print(); +// return; + +// tr->SetMarkerColor(kRed); +// tr->Draw("(p1.fImpulse.T()-p1.fMass)/p1.fA:p1.fImpulse.Theta()*TMath::RadToDeg()>>hKin", "", "scat"); + tr->Draw("(p1.fImpulse.T()-p1.fMass)/p1.fA:p1.fImpulse.Theta()*TMath::RadToDeg()>>hKin", "", ""); + hKin->SetMarkerColor(kRed); + hKin->SetMarkerStyle(20); + hKin->SetMarkerSize(0.4); + + tr->SetMarkerColor(kRed); + tr->Draw("(p2.fImpulse.T()-p2.fMass)/p2.fA:p2.fImpulse.Theta()*TMath::RadToDeg()", "", "same"); + tr->SetMarkerColor(kGreen); + tr->Draw("(4_He.fImpulse.T()-4_He.fMass)/4_He.fA:4_He.fImpulse.Theta()*TMath::RadToDeg()", "", "same"); + tr->SetMarkerColor(kBlue); + tr->Draw("(6_Be.fImpulse.T()-6_Be.fMass)/6_Be.fA:6_Be.fImpulse.Theta()*TMath::RadToDeg()", "", "same"); + + + + hKin->GetXaxis()->SetTitle("#it{#theta}_{lab} (deg)"); + hKin->GetXaxis()->CenterTitle(); + hKin->GetXaxis()->SetTitleSize(0.07); + hKin->GetXaxis()->SetTitleOffset(0.79); + hKin->GetXaxis()->SetLabelSize(0.05); + + hKin->GetYaxis()->SetTitle("#it{E}_{lab} (AMeV)"); + hKin->GetYaxis()->CenterTitle(); + hKin->GetYaxis()->SetTitleSize(0.07); + hKin->GetYaxis()->SetTitleOffset(0.79); + hKin->GetYaxis()->SetLabelSize(0.05); + + cKin->cd(); + + cKin->Update(); + TCanvas *c1 = new TCanvas(); + c1->Update(); + delete c1; + cKin->Update(); + + if (savePictures) { +// gStyle->SetCanvasDefW(649); +// gStyle->SetCanvasDefH(727); + pictureName.Form("%s%s%s", picturePath.Data(), cKin->GetName(), pictFormat.Data()); +// cKin->SaveAs(pictureName.Data()); + cKin->SaveAs("figures/PresKinematics.png"); + cKin->cd(); +// cKin->Close(); + } +// } + +} -- 2.18.1