Changeset 1228 for trunk/source/processes/hadronic/models/util/include/G4NuclearFermiDensity.hh

Timestamp:

Jan 8, 2010, 11:56:51 AM (14 years ago)

Author:

garnier

Message:

update geant4.9.3 tag

File:

• trunk/source/processes/hadronic/models/util/include/G4NuclearFermiDensity.hh (modified) (1 diff)

trunk/source/processes/hadronic/models/util/include/G4NuclearFermiDensity.hh

@@ -56 +56 @@
     G4double GetDeriv(const G4ThreeVector & aPosition) const
     {
-      return -std::exp((aPosition.mag()-theR)/a) * sqr(GetDensity(aPosition)) / (a*Getrho0());
+      G4double currentR=aPosition.mag();
+      if (currentR > 40*theR  ) {return 0;}
+      else return -std::exp((currentR-theR)/a) * sqr(GetDensity(aPosition)) / (a*Getrho0());
     }