diff --git a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
index b7fadcdd975093f5b7daba0aeef747cf4aa37330..b14c07996e17b9f8887243b8e9a7853493ef2ec4 100644
--- a/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
+++ b/benchmarks/barrel_ecal/scripts/emcal_barrel_pions_analysis.cxx
@@ -128,6 +128,13 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
     return result;
   };
 
+  // Thrown Energy [GeV]
+  auto Ethr = [](std::vector<dd4pod::Geant4ParticleData> const& input) {
+    std::vector<double> result;
+    result.push_back(input[2].charge);
+  return result;
+  };
+
   // Define variables
   auto d1 = d0.Define("Ethr",  Ethr,  {"mcparticles"})
 	      .Define("nhits",  nhits,      {"EcalBarrelHits"})
@@ -136,7 +143,7 @@ void emcal_barrel_pions_analysis(const char* input_fname = "sim_output/sim_emcal
         .Define("fsam",   fsam2,      {"Esim","Ethr"})
         .Define("dE",     eResol,     {"Esim", "Ethr"})
         .Define("dE_rel", eResol_rel, {"Esim", "Ethr"})
-        .Define("charge", {"charge"})
+        .Define("charge", {"mcparticles"})
 	      ;
 
   // Define Histograms