From 3edc704872dafc95e0b021422f07e5aa79d2b792 Mon Sep 17 00:00:00 2001
From: fysquare <freddy.flavigny@laposte.net>
Date: Fri, 3 Jul 2020 09:55:38 +0200
Subject: [PATCH] Fix inversion in NPQFS between frag. and removed cluster
---
NPLib/Physics/NPQFS.cxx | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/NPLib/Physics/NPQFS.cxx b/NPLib/Physics/NPQFS.cxx
index d3156d6d9..663cc53e0 100644
--- a/NPLib/Physics/NPQFS.cxx
+++ b/NPLib/Physics/NPQFS.cxx
@@ -227,13 +227,13 @@ void QFS::CalculateVariables(){
double EA = sqrt(mA*mA + PA*PA); // Beam total energy
fEnergyImpulsionLab_A = TLorentzVector(0.,0.,PA,EA);
- //Internal momentum of removed cluster/nucleon
- Pa.SetX(fInternalMomentum.X());
- Pa.SetY(fInternalMomentum.Y());
- Pa.SetZ(fInternalMomentum.Z());
-
- //Internal momentum of heavy recoil after removal
- PB.SetXYZ( (-Pa.X()) , (-Pa.Y()) , (-Pa.Z()) );
+ // Internal momentum of removed cluster/nucleon (Pa) and recoil (PB)
+ // here fInternalMomentum contains PB (recoil fragment momentum)
+ // readout from the input file (theoretical)
+ PB.SetX(fInternalMomentum.X());
+ PB.SetY(fInternalMomentum.Y());
+ PB.SetZ(fInternalMomentum.Z());
+ Pa.SetXYZ( (-PB.X()) , (-PB.Y()) , (-PB.Z()) );
// Off-shell mass of the bound nucleon from E conservation
// in virtual dissociation of A -> B + a
--
GitLab