33 #define INCLXX_IN_GEANT4_MODE 1
53 namespace NuclearPotential {
55 const G4double NuclearPotentialEnergyIsospinSmooth::alpha= 0.223;
56 const G4double NuclearPotentialEnergyIsospinSmooth::deltaE= 25.;
78 const G4double t0 = tf + v0*(1.-alpha)/alpha - deltaE;
81 v = v0 - alpha*(t-tf)/(1.-alpha);
83 const G4double v_at_t0 = v0 - alpha*(t0-tf)/(1.-alpha);
84 const G4double kappa = alpha / (v_at_t0 * (1.-alpha));
85 v = v_at_t0 * std::exp(kappa * (t0-t));
87 return (v>0.0) ? v : 0.0;