Changeset 962 for trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPKallbachMannSyst.cc

Timestamp:

Apr 6, 2009, 12:30:29 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPKallbachMannSyst.cc (modified) (2 diffs)

trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPKallbachMannSyst.cc

@@ -27 +27 @@
 // J.P. Wellisch, Nov-1996
 // A prototype of the low energy neutron transport model.
+//
+// 080801 Protect div0 error, when theCompundFraction is 1 by T. Koi
 //
 #include "G4NeutronHPKallbachMannSyst.hh" 
@@ -69 +71 @@
 {
   G4double result;
+  if ( theCompoundFraction == 1 ) 
+  { 
+     //G4cout << "080730b Adjust theCompoundFraction " << G4endl;
+     theCompoundFraction *= (1-1.0e-15);   
+  } 
   result = 0.5 * (1./A(anEnergy)) * std::log((1-theCompoundFraction)/(1+theCompoundFraction));
   return result;