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

commented out fitting resolution parameters

parent 8cf8a0f4
No related branches found
No related tags found
1 merge request!24Resolve "pi0 resolution for ECal barrel"
This commit is part of merge request !24. Comments created here will be created in the context of that merge request.
...@@ -181,7 +181,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -181,7 +181,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
auto hdE = d1.Histo1D({"hdE", "dE; dE[GeV]; Events", 20, -7.5, 7.5}, "dE");//changed from 100 auto hdE = d1.Histo1D({"hdE", "dE; dE[GeV]; Events", 20, -7.5, 7.5}, "dE");//changed from 100
auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 20, -7.5, 7.5}, "dE_rel");//changed from 100 auto hdE_rel = d1.Histo1D({"hdE_rel", "dE Relative; dE Relative; Events", 20, -7.5, 7.5}, "dE_rel");//changed from 100
hdE_rel->Fit("gaus", "", "", -7.5, 7.5); hdE_rel->Fit("gaus", "", "", -7.5, 7.5);
double* res = hdE_rel->GetFunction("gaus")->GetParameters(); //double* res = hdE_rel->GetFunction("gaus")->GetParameters();
//Pass/Fail //Pass/Fail
/* /*
...@@ -191,7 +191,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -191,7 +191,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
pion0_energy_resolution.fail(res[2]); pion0_energy_resolution.fail(res[2]);
} }
*/ */
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);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment