Fixed apparent bug in simulating neutron-emitting reactions

6d67e755 · Greg Christian · 2dedcc41 · 6d67e755
Commit 6d67e755 authored Oct 02, 2018 by Greg Christian
Hide whitespace changes
Inline Side-by-side

Showing with 11 additions and 4 deletions

NPSimulation/Process/BeamReaction.cc NPSimulation/Process/BeamReaction.cc +11 -4

No files found.
--- a/NPSimulation/Process/BeamReaction.cc
+++ b/NPSimulation/Process/BeamReaction.cc
@@ -184,10 +184,17 @@ void NPS::BeamReaction::DoIt(const G4FastTrack& fastTrack,G4FastStep& fastStep)
   
    G4ParticleDefinition* LightName;

-    if(m_Reaction.GetUseExInGeant4())
-     LightName = IonTable->GetIon(LightZ, LightA, m_Reaction.GetExcitation3()*MeV);
-    else 
-     LightName = IonTable->GetIon(LightZ, LightA);
+		if(LightZ == 0 && LightA == 1) // neutron is special case
+		{
+			LightName = G4Neutron::Definition();
+		}
+		else
+		{	
+			if(m_Reaction.GetUseExInGeant4())
+				LightName = IonTable->GetIon(LightZ, LightA, m_Reaction.GetExcitation3()*MeV);
+			else 
+				LightName = IonTable->GetIon(LightZ, LightA);
+		}

    // Nucleus 4
    G4int HeavyZ = m_Reaction.GetNucleus4()->GetZ() ;