daughters

benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx

@@ -41,7 +41,8 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
   std::string test_tag = "Barrel_emcal_pi0";
   //TODO: Change test_tag to something else
   std:string detector = "Barrel_emcal";
-  //Energy resolution in teh barrel region(-1 < eta < 1)
+  //Energy resolution in the barrel region(-1 < eta < 1)
+  // Taken from : Initial considerations for EMCal of the EIC detector by A. Bazilevsky
   // sigma_E / E = 12% / E^0.5 convoluted with 2%
   // sigma_E / E = [ (0.12/E^0.5)^2 + 0.02^2]^0.5, with E in [GeV]
   double thrown_energy = 5; // Current thrown energy, will need to grab from json file
@@ -122,12 +123,11 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
 
   // Returns the pdgID of the particle daughters(hopefully?)
   auto getdau = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
-    std::vector<double> result;
-    /*
+    std::vector<int> result;
+    
     for (auto part : input[2].daughters_begin)
-      result.push_back(*part)
-      */
-    return *input[2].daughters_begin;
+      result.push_back(part)
+    return result;
   };
 
   // Define variables