From 191caabb7537a1d293a8c431382b1b52a9f8fd5e Mon Sep 17 00:00:00 2001
From: morfouace <pierre.morfouace@cea.fr>
Date: Wed, 6 Nov 2024 16:34:17 +0100
Subject: [PATCH] Debugging NPReaction.cxx

---
 NPLib/Physics/NPReaction.cxx | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/NPLib/Physics/NPReaction.cxx b/NPLib/Physics/NPReaction.cxx
index ecf1c4d35..6f890404c 100644
--- a/NPLib/Physics/NPReaction.cxx
+++ b/NPLib/Physics/NPReaction.cxx
@@ -301,7 +301,7 @@ double Reaction::ReconstructRelativistic(double EnergyLab, double ThetaLab, doub
   double E3 = m3 + EnergyLab;
   double p_Lab_3 = sqrt(E3 * E3 - m3 * m3);
   fEnergyImpulsionLab_3 =
-      TLorentzVector(p_Lab_3 * sin(ThetaLab) * cos(PhiLab), p_Lab_3 * sin(PhiLab), p_Lab_3 * cos(ThetaLab), E3);
+      TLorentzVector(p_Lab_3 * sin(ThetaLab) * cos(PhiLab), p_Lab_3 * sin(ThetaLab) * sin(PhiLab), p_Lab_3 * cos(ThetaLab), E3);
   fEnergyImpulsionLab_4 = fTotalEnergyImpulsionLab - fEnergyImpulsionLab_3;
 
   double Eex = fEnergyImpulsionLab_4.Mag() - fParticle4.Mass();
@@ -313,7 +313,7 @@ TLorentzVector Reaction::LorentzAfterReaction(double EnergyLab, double ThetaLab,
   double E3 = m3 + EnergyLab;
   double p_Lab_3 = sqrt(E3 * E3 - m3 * m3);
   fEnergyImpulsionLab_3 =
-      TLorentzVector(p_Lab_3 * sin(ThetaLab) * cos(PhiLab), p_Lab_3 * sin(PhiLab), p_Lab_3 * cos(ThetaLab), E3);
+      TLorentzVector(p_Lab_3 * sin(ThetaLab) * cos(PhiLab), p_Lab_3 * sin(ThetaLab) * sin(PhiLab), p_Lab_3 * cos(ThetaLab), E3);
   fEnergyImpulsionLab_4 = fTotalEnergyImpulsionLab - fEnergyImpulsionLab_3;
   return fEnergyImpulsionLab_4;
 }
-- 
GitLab