Skip to content
Snippets Groups Projects

Resolve "Fix Pi0 benchmark"

Merged Marshall Scott requested to merge 37-fix-pi0-benchmark into master
1 file
+ 3
3
Compare changes
  • Side-by-side
  • Inline
@@ -177,7 +177,7 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
@@ -177,7 +177,7 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
c6->SaveAs("results/emcal_barrel_pi0_dau.pdf");
c6->SaveAs("results/emcal_barrel_pi0_dau.pdf");
// Energy Resolution Calculation
// Energy Resolution Calculation
std::string test_tag = "Barrel_emcal_pi0";//TODO: Change test_tag to something else
std::string test_tag = "Barrel_emcal_pi0";// TODO: Change test_tag to something else
std:string detEle = "Barrel_emcal";
std:string detEle = "Barrel_emcal";
// Energy resolution in the barrel region (-1 < eta < 1)
// Energy resolution in the barrel region (-1 < eta < 1)
@@ -203,11 +203,11 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
@@ -203,11 +203,11 @@ void emcal_barrel_pi0_analysis(const char* input_fname = "sim_output/sim_emcal_b
double* res = hdE->GetFunction("gaus")->GetParameters();
double* res = hdE->GetFunction("gaus")->GetParameters();
double sigmaOverE = res[2] / meanE;
double sigmaOverE = res[2] / meanE;
//Pass/Fail
// Pass/Fail
sigmaOverE <= resolutionTarget ? pi0_energy_resolution.pass(sigmaOverE) : pi0_energy_resolution.fail(sigmaOverE);
sigmaOverE <= resolutionTarget ? pi0_energy_resolution.pass(sigmaOverE) : pi0_energy_resolution.fail(sigmaOverE);
std::printf("Energy Resolution is %f\n", 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);
cdE->SetLogy(1);
cdE->SetLogy(1);
hdE->GetYaxis()->SetTitleOffset(1.4);
hdE->GetYaxis()->SetTitleOffset(1.4);
Loading