Skip to content
Snippets Groups Projects
Commit 71af800e authored by Marshall Scott's avatar Marshall Scott
Browse files

Checking if benchmark testing is the issue

parent 772ec88c
No related branches found
No related tags found
No related merge requests found
...@@ -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));
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment