Skip to content
Snippets Groups Projects

Resolve "pi0 resolution for ECal barrel"

Merged Marshall Scott requested to merge 19-pi0-resolution-for-ecal-barrel into master
Compare and Show latest version
1 file
+ 6
2
Compare changes
  • Side-by-side
  • Inline
@@ -41,6 +41,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
@@ -41,6 +41,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
//TODO: Change test_tag to something else
//TODO: Change test_tag to something else
std:string detector = "Barrel_emcal";
std:string detector = "Barrel_emcal";
double resolutionTarget = 0.1;
double resolutionTarget = 0.1;
 
/*
eic::util::Test pion0_Energy_resolution{
eic::util::Test pion0_Energy_resolution{
{{"name", fmt::format("{}_energy_resolution", test_tag)},
{{"name", fmt::format("{}_energy_resolution", test_tag)},
{"title", "Pion0 Energy resolution"},
{"title", "Pion0 Energy resolution"},
@@ -48,7 +49,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
@@ -48,7 +49,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
fmt::format("Pion0 energy resolution with {}, estimated using a Gaussian fit.", detector)},
fmt::format("Pion0 energy resolution with {}, estimated using a Gaussian fit.", detector)},
{"quantity", "resolution (in %)"},
{"quantity", "resolution (in %)"},
{"target", std::to_string(resolutionTarget)}}};
{"target", std::to_string(resolutionTarget)}}};
*/
ROOT::EnableImplicitMT();
ROOT::EnableImplicitMT();
ROOT::RDataFrame d0("events", input_fname);
ROOT::RDataFrame d0("events", input_fname);
@@ -148,11 +149,14 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
@@ -148,11 +149,14 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
const double* res = f1->GetParams();
const double* res = f1->GetParams();
//Pass/Fail
//Pass/Fail
 
/*
if (res[2] <= resolutionTarget) {
if (res[2] <= resolutionTarget) {
pion0_energy_resolution.pass(res[2]);
pion0_energy_resolution.pass(res[2]);
} else {
} else {
pion0_energy_resolution.fail(res[2]);
pion0_energy_resolution.fail(res[2]);
}
}
 
*/
 
std::printf("Energy Resolution is %f\n", res[2]);
//Energy Resolution Histogram Plotting
//Energy Resolution Histogram Plotting
auto *cdE = new TCanvas("cdE", "cdE", 700, 500);
auto *cdE = new TCanvas("cdE", "cdE", 700, 500);
@@ -175,7 +179,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
@@ -175,7 +179,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.png");
cdE_rel->SaveAs("results/emcal_barrel_pi0_dE_rel.pdf");
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));
}
}
Loading