Commit 3d649d46 authored by Vratislav Chudoba's avatar Vratislav Chudoba

New scripts for presentation use.

parent 5cf02059
This diff is collapsed.
//#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");
}
}
//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");
}
}
//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");
}
}
//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");
}
}
#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();
// }
// }
}
//#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");
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment