G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth Class Reference

#include <G4INCLNuclearPotentialEnergyIsospinSmooth.hh>

Inheritance diagram for G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth:

Public Member Functions
	NuclearPotentialEnergyIsospinSmooth (const G4int A, const G4int Z, const G4bool pionPotential)
virtual	~NuclearPotentialEnergyIsospinSmooth ()
virtual G4double	computePotentialEnergy (const Particle *const p) const
Public Member Functions inherited from G4INCL::NuclearPotential::NuclearPotentialIsospin
	NuclearPotentialIsospin (const G4int A, const G4int Z, const G4bool pionPotential)
virtual	~NuclearPotentialIsospin ()
Public Member Functions inherited from G4INCL::NuclearPotential::INuclearPotential
	INuclearPotential (const G4int A, const G4int Z, const G4bool pionPot)
virtual	~INuclearPotential ()
G4bool	hasPionPotential () const
	Do we have a pion potential? More...
G4double	getFermiEnergy (const Particle *const p) const
	Return the Fermi energy for a particle. More...
G4double	getFermiEnergy (const ParticleType t) const
	Return the Fermi energy for a particle type. More...
G4double	getSeparationEnergy (const Particle *const p) const
	Return the separation energy for a particle. More...
G4double	getSeparationEnergy (const ParticleType t) const
	Return the separation energy for a particle type. More...
G4double	getFermiMomentum (const Particle *const p) const
	Return the Fermi momentum for a particle. More...
G4double	getFermiMomentum (const ParticleType t) const
	Return the Fermi momentum for a particle type. More...

Additional Inherited Members
Protected Member Functions inherited from G4INCL::NuclearPotential::INuclearPotential
G4double	computePionPotentialEnergy (const Particle *const p) const
	Compute the potential energy for the given pion. More...
Protected Attributes inherited from G4INCL::NuclearPotential::INuclearPotential
const G4int	theA
	The mass number of the nucleus. More...
const G4int	theZ
	The charge number of the nucleus. More...
std::map< ParticleType, G4double >	fermiEnergy
std::map< ParticleType, G4double >	fermiMomentum
std::map< ParticleType, G4double >	separationEnergy

Detailed Description

Definition at line 56 of file G4INCLNuclearPotentialEnergyIsospinSmooth.hh.

Constructor & Destructor Documentation

G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::NuclearPotentialEnergyIsospinSmooth	(	const G4int	A,
		const G4int	Z,
		const G4bool	pionPotential
	)

Definition at line 59 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

      : NuclearPotentialIsospin(A,Z,aPionPotential)
    {}

G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::~NuclearPotentialEnergyIsospinSmooth ( )

virtual

Definition at line 64 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

{}

Member Function Documentation

G4double G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy ( const Particle *const  p ) const

virtual

Reimplemented from G4INCL::NuclearPotential::NuclearPotentialIsospin.

Definition at line 66 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

References G4INCL::NuclearPotential::NuclearPotentialIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::Particle::getKineticEnergy(), G4INCL::Particle::isNucleon(), and test::v.

                                                                                                       {

      const G4double v0 = NuclearPotentialIsospin::computePotentialEnergy(particle);

      if(particle->isNucleon()) {
        const G4double t = particle->getKineticEnergy();
        const G4double tf = getFermiEnergy(particle);
        // Constant potential for T<Tf
        if(t < tf)
          return v0;

        // Linear function for Tf<T<T0, exponential function for T>T0
        const G4double t0 = tf + v0*(1.-alpha)/alpha - deltaE; // deltaE before the linear potential vanishes
        G4double v;
        if(t<t0) {
          v = v0 - alpha*(t-tf)/(1.-alpha);
        } else {
          const G4double v_at_t0 = v0 - alpha*(t0-tf)/(1.-alpha);
          const G4double kappa = alpha / (v_at_t0 * (1.-alpha));
          v = v_at_t0 * std::exp(kappa * (t0-t));
        }
        return (v>0.0) ? v : 0.0;
      } else
        return v0;
    }

---

The documentation for this class was generated from the following files:

• G4INCLNuclearPotentialEnergyIsospinSmooth.hh
• G4INCLNuclearPotentialEnergyIsospinSmooth.cc