From 5e1fe07e695b5aaa9b8070373bed2a791e5ee3a9 Mon Sep 17 00:00:00 2001 From: Vratislav Chudoba Date: Fri, 4 Aug 2017 14:15:38 +0300 Subject: [PATCH] Figures of 2+ left slope --- macros/BeCorrPRC/figCorrLeftES.cxx | 212 +++++++++++++++++++++++++++++ 1 file changed, 212 insertions(+) create mode 100755 macros/BeCorrPRC/figCorrLeftES.cxx diff --git a/macros/BeCorrPRC/figCorrLeftES.cxx b/macros/BeCorrPRC/figCorrLeftES.cxx new file mode 100755 index 0000000..26f9c43 --- /dev/null +++ b/macros/BeCorrPRC/figCorrLeftES.cxx @@ -0,0 +1,212 @@ +void figCorrLeftES() +{ + + 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 *c3 = new TCanvas(); + c3->ToggleEditor(); + c3->SetWindowSize(800, 1227); + c3->Divide(2, 3, .001, .001); + + TFile *fr2 = new TFile("figures/PRCthetaA.root"); + + TH1F *hAlphaExp[6]; + TH1F *hAlphaSim[6]; + TH1F *hAlphaInput[6]; + + TString hName; + + for (Int_t i = 0; i < 6; i++) { + c3->cd(i+1); + + 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); + + 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(""); + hAlphaExp[i]->Draw("E same"); + hAlphaInput[i]->Draw("same"); + } + + + c3->SaveAs("figures/ESThetaAleft7590.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"); + +} -- 2.18.1