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

Changed binning

parent 50d85eef
No related branches found
No related tags found
1 merge request!24Resolve "pi0 resolution for ECal barrel"
Show whitespace changes
Inline Side-by-side
benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
+ 3
− 4
View file @ c88e7ac2
...@@ -213,14 +213,13 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -213,14 +213,13 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
c6->SaveAs("results/emcal_barrel_pions_dau.pdf"); c6->SaveAs("results/emcal_barrel_pions_dau.pdf");
//Energy Resolution Calculation //Energy Resolution Calculation
auto hdE = d2.Histo1D({"hdE", "dE; dE[GeV]; Events", 20, -2.5, 2.5}, "dE");//changed from 100 auto hdE = d2.Histo1D({"hdE", "dE; dE[GeV]; Events", 100, -3.0, 3.0}, "dE");
auto hdE_rel = d2.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 20, -2.5, 2.5}, "dE_rel");//changed from 100 auto hdE_rel = d2.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 100, -3.0, 3.0}, "dE_rel");
hdE->Fit("gaus", "", "", -2.5, 2.5); hdE->Fit("gaus", "", "", -3.0, 3.0);
double* res = hdE->GetFunction("gaus")->GetParameters(); double* res = hdE->GetFunction("gaus")->GetParameters();
double sigmaOverE = res[2] / thrown_energy; double sigmaOverE = res[2] / thrown_energy;
//Pass/Fail //Pass/Fail
if (sigmaOverE <= resolutionTarget) { if (sigmaOverE <= resolutionTarget) {
pi0_energy_resolution.pass(sigmaOverE); pi0_energy_resolution.pass(sigmaOverE);
} else { } else {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment