Resolve "Add layer energy-deposit profile to energy resolution benchmark"
Merged
requested to merge 26-add-layer-energy-deposit-profile-to-energy-resolution-benchmark into master
Compare changes
@@ -18,6 +18,10 @@
@@ -18,6 +18,10 @@
@@ -44,7 +48,7 @@ void save_canvas(TCanvas* c, std::string var_label, std::string E_label, std::st
@@ -44,7 +48,7 @@ void save_canvas(TCanvas* c, std::string var_label, std::string E_label, std::st
std::tuple <double, double, double, double> extract_sampling_fraction_parameters(std::string particle_label, std::string E_label)
std::string input_fname = fmt::format("sim_output/energy_scan/{}/sim_emcal_barrel_{}.root", E_label, particle_label);
@@ -73,6 +77,25 @@ std::tuple <double, double, double, double> extract_sampling_fraction_parameters
@@ -73,6 +77,25 @@ std::tuple <double, double, double, double> extract_sampling_fraction_parameters
@@ -91,8 +114,35 @@ std::tuple <double, double, double, double> extract_sampling_fraction_parameters
@@ -91,8 +114,35 @@ std::tuple <double, double, double, double> extract_sampling_fraction_parameters
@@ -190,38 +240,51 @@ void emcal_barrel_energy_scan_analysis(std::string particle_label = "electron")
@@ -190,38 +240,51 @@ void emcal_barrel_energy_scan_analysis(std::string particle_label = "electron")
auto scanned_energies = read_scanned_energies(fmt::format("sim_output/emcal_barrel_energy_scan_points_{}.txt", particle_label));
auto [fsam, fsam_res, fsam_err, fsam_res_err] = extract_sampling_fraction_parameters(particle_label, E_label);
auto fsam_res_rel_err = 100.0*(sqrt(pow((fsam_res_err/fsam),2)+pow((fsam_err*fsam_res)/(fsam*fsam),2)));