Skip to content
Snippets Groups Projects

Resolve "pi0 resolution for ECal barrel"

Merged Marshall Scott requested to merge 19-pi0-resolution-for-ecal-barrel into master
Compare and Show latest version
1 file
+ 10
1
Compare changes
  • Side-by-side
  • Inline
@@ -20,6 +20,8 @@
#include "TF1.h"
#include "TH1D.h"
#include <TFitResult.h>
#include "TSystem.h"
#include <fstream>
using ROOT::RDataFrame;
using namespace ROOT::VecOps;
@@ -51,10 +53,17 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
{"quantity", "resolution (in %)"},
{"target", std::to_string(resolutionTarget)}}};
*/
gSystem->Exec("cp sim_output/sim_emcal_barrel_uniform_pions.root results/sim_emcal_barrel_uniform_pions.root");
//gSystem->Exec("cp sim_output/sim_emcal_barrel_uniform_pions.root results/sim_emcal_barrel_uniform_pions.root");
ROOT::EnableImplicitMT();
ROOT::RDataFrame d0("events", input_fname);
auto colNames = d0.GetColumnNames();
ofstream out;
out.open("results/column_list.txt");
for (auto &&col : colNames){out << col << endl; cout << col << endl;}
out.close();
// Thrown Energy [GeV]
auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
std::vector<double> result;
Loading