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

Fixed according to local branch

parent faf58e2f
No related branches found
No related tags found
No related merge requests found
This commit is part of merge request !51. Comments created here will be created in the context of that merge request.
...@@ -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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment