Fixing ordering of four momentum in SecondaryParticleRecord

Harvard-Neutrino · Feb 26, 2024 · 7bb6359 · 7bb6359
1 parent 19fd4f2
commit 7bb6359
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/projects/dataclasses/private/InteractionRecord.cxx b/projects/dataclasses/private/InteractionRecord.cxx
@@ -440,7 +440,10 @@ std::array<double, 4> SecondaryParticleRecord::GetFourMomentum() const {
     if(not momentum_set) {
         UpdateMomentum();
     }
-    return {momentum.at(0), momentum.at(1), momentum.at(2), GetEnergy()};
+    if(not energy_set) {
+        UpdateEnergy();
+    }
+    return {energy, momentum.at(0), momentum.at(1), momentum.at(2)};
 }
 
 std::array<double, 3> const & SecondaryParticleRecord::GetInitialPosition() const {