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

checking include issues 2

parent a931ef85
Branches
Tags
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.
...@@ -8,11 +8,11 @@ ...@@ -8,11 +8,11 @@
#include "dd4pod/Geant4ParticleCollection.h" #include "dd4pod/Geant4ParticleCollection.h"
#include "dd4pod/CalorimeterHitCollection.h" #include "dd4pod/CalorimeterHitCollection.h"
/*
#include "benchmark.h" #include "benchmark.h"
#include "mt.h" #include "mt.h"
#include "util.h" #include "util.h"
*/
#include "TCanvas.h" #include "TCanvas.h"
#include "TStyle.h" #include "TStyle.h"
#include "TMath.h" #include "TMath.h"
...@@ -42,7 +42,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -42,7 +42,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
//TODO: Change test_tag to something else //TODO: Change test_tag to something else
std:string detector = "Barrel_emcal"; std:string detector = "Barrel_emcal";
double resolutionTarget = 0.1; double resolutionTarget = 0.1;
/*
eic::util::Test pion0_Energy_resolution{ eic::util::Test pion0_Energy_resolution{
{{"name", fmt::format("{}_energy_resolution", test_tag)}, {{"name", fmt::format("{}_energy_resolution", test_tag)},
{"title", "Pion0 Energy resolution"}, {"title", "Pion0 Energy resolution"},
...@@ -50,7 +50,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -50,7 +50,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
fmt::format("Pion0 energy resolution with {}, estimated using a Gaussian fit.", detector)}, fmt::format("Pion0 energy resolution with {}, estimated using a Gaussian fit.", detector)},
{"quantity", "resolution (in %)"}, {"quantity", "resolution (in %)"},
{"target", std::to_string(resolutionTarget)}}}; {"target", std::to_string(resolutionTarget)}}};
*/
ROOT::EnableImplicitMT(); ROOT::EnableImplicitMT();
ROOT::RDataFrame d0("events", input_fname); ROOT::RDataFrame d0("events", input_fname);
...@@ -105,7 +105,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal ...@@ -105,7 +105,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
return result; return result;
}; };
// Relative Energy Resolution = (Esampling - Ethrown)/Ethrownnnn // Relative Energy Resolution = (Esampling - Ethrown)/Ethrown
auto eResol_rel = [](const std::vector<double>& sampled, const std::vector<double>& thrown) { auto eResol_rel = [](const std::vector<double>& sampled, const std::vector<double>& thrown) {
std::vector<double> result; std::vector<double> result;
auto it_sam = sampled.cbegin(); auto it_sam = sampled.cbegin();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment