From 772ec88cd0f77c8455dd327f49b0d81922d35a0b Mon Sep 17 00:00:00 2001
From: Marshall Scott <mbscott@anl.gov>
Date: Mon, 26 Apr 2021 15:27:45 -0400
Subject: [PATCH] Fixed naming
---
.../scripts/emcal_barrel_pions_analysis.cxx | 13 +++++++------
1 file changed, 7 insertions(+), 6 deletions(-)
diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
index 925141df..bd7b15be 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
@@ -96,7 +96,6 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
auto hNhits = d1.Histo1D({"hNhits", "Number of hits per events; Number of hits; Events", 100, 0.0, 2000.0}, "nhits");
auto hEsim = d1.Histo1D({"hEsim", "Energy Deposit; Energy Deposit [GeV]; Events", 100, 0.0, 1.0}, "Esim");
auto hfsam = d1.Histo1D({"hfsam", "Sampling Fraction; Sampling Fraction; Events", 100, 0.0, 0.1}, "fsam");
-
// Event Counts
auto nevents_thrown = d1.Count();
@@ -143,16 +142,19 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
c4->SaveAs("results/emcal_barrel_pions_fsam.pdf");
//Energy Resolution Work
- auto hdE = d1.Histo1D({"hdE", "dE; dE[GeV]; Events", 100, -7.5, 7.5}, "dE");
- auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -7.5, 7.5}, "dE_rel");
- auto f1 = hdE_rel->Fit("gaus", "S");
+ auto hdE = d1.Histo1D({"hdE", "dE; dE[GeV]; Events", 100, -7.5, 7.5}, "dE");
+ auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -7.5, 7.5}, "dE_rel");
+ auto f1 = hdE_rel->Fit("gaus", "S");
const double* res = f1->GetParams();
+
+ //Pass/Fail
if (res[2] <= resolutionTarget) {
pion0_energy_resolution.pass(res[2]);
} else {
pion0_energy_resolution.fail(res[2]);
}
+ //Energy Resolution Histogram Plotting
auto *cdE = new TCanvas("cdE", "cdE", 700, 500);
cDE->SetLogy(1);
hdE->GetYaxis()->SetTitleOffset(1.4);
@@ -163,7 +165,6 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
cdE->SaveAs("results/emcal_barrel_pi0_dE.pdf");
auto *cdE_rel = new TCanvas("cdE_rel", "cdE_rel", 700, 500);
- cdE_rel->SetLogy(1);
hdE_rel->GetYaxis()->SetTitleOffset(1.4);
hdE_rel->SetLineWidth(2);
hdE_rel->SetLineColor(kBlue);
@@ -174,7 +175,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.png");
cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.pdf");
- eic::util::write_test({pion0_Energy_resolution}, fmt::format("{}_pions.json", detector));
+ eic::util::write_test({pion0_energy_resolution}, fmt::format("{}_pions.json", detector));
}
--
GitLab