From 67e7be4360ad8cf1fa63d09e7e7b5fcc86b2595e Mon Sep 17 00:00:00 2001 From: Morfouace <pierre.morfouace@gmail.com> Date: Thu, 3 Feb 2022 17:37:29 +0100 Subject: [PATCH] Changing ParentID in NPFissionDecau so that the beam can fission --- NPLib/Physics/NPFissionDecay.cxx | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/NPLib/Physics/NPFissionDecay.cxx b/NPLib/Physics/NPFissionDecay.cxx index 5b7cb1f7a..7a0713866 100644 --- a/NPLib/Physics/NPFissionDecay.cxx +++ b/NPLib/Physics/NPFissionDecay.cxx @@ -146,16 +146,16 @@ bool NPL::FissionDecay::GenerateEvent(string CompoundName, double MEx,double MEK double Phil = m_FissionModel->GetPhffl(); double Phih = m_FissionModel->GetPhffh(); - TVector3 uxy = TVector3(cos(TMath::Pi()/2-PhiCN), -sin(TMath::Pi()/2-PhiCN), 0); + TVector3 uxy = TVector3(cos(TMath::Pi()/2-PhiCN), sin(TMath::Pi()/2-PhiCN), 0); TVector3 Momentuml = Pl * TVector3(sin(Thetal)*cos(Phil), sin(Thetal)*sin(Phil), cos(Thetal)); - //Momentuml.Rotate(-ThetaCN, uxy); + Momentuml.Rotate(-ThetaCN, uxy); TVector3 Momentumh = Ph * TVector3(sin(Thetah)*cos(Phih), sin(Thetah)*sin(Phih), cos(Thetah)); - //Momentumh.Rotate(-ThetaCN, uxy); + Momentumh.Rotate(-ThetaCN, uxy); DPx.push_back(Momentuml.X()); DPx.push_back(Momentumh.X()); -- GitLab