#include <G4INCLNuclearPotentialEnergyIsospin.hh>
Inheritance diagram for G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin:
Public Member Functions | |
NuclearPotentialEnergyIsospin (const G4int A, const G4int Z, const G4bool pionPotential) | |
virtual | ~NuclearPotentialEnergyIsospin () |
virtual G4double | computePotentialEnergy (const Particle *const p) const |
Definition at line 55 of file G4INCLNuclearPotentialEnergyIsospin.hh.
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::NuclearPotentialEnergyIsospin | ( | const G4int | A, | |
const G4int | Z, | |||
const G4bool | pionPotential | |||
) |
Definition at line 57 of file G4INCLNuclearPotentialEnergyIsospin.cc.
00058 : NuclearPotentialIsospin(A,Z,aPionPotential) 00059 {}
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::~NuclearPotentialEnergyIsospin | ( | ) | [virtual] |
G4double G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy | ( | const Particle *const | p | ) | const [virtual] |
Reimplemented from G4INCL::NuclearPotential::NuclearPotentialIsospin.
Definition at line 64 of file G4INCLNuclearPotentialEnergyIsospin.cc.
References G4INCL::NuclearPotential::NuclearPotentialIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::Particle::getKineticEnergy(), and G4INCL::Particle::isNucleon().
00064 { 00065 00066 const G4double v0 = NuclearPotentialIsospin::computePotentialEnergy(particle); 00067 00068 if(particle->isNucleon()) { 00069 const G4double t = particle->getKineticEnergy(); 00070 const G4double tf = getFermiEnergy(particle); 00071 // Constant potential for T<Tf 00072 if(t < tf) 00073 return v0; 00074 00075 // Linear function for T>Tf 00076 const G4double v = v0 - alpha*(t-tf)/(1-alpha); 00077 return (v>0.0) ? v : 0.0; // return 0.0 if v is negative 00078 } else 00079 return v0; 00080 }