#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 ()

virtual G4double	computePotentialEnergy (const Particle *const p) const

Public Member Functions inherited from G4INCL::NuclearPotential::INuclearPotential
	INuclearPotential (const G4int A, const G4int Z, const G4bool pionPot)

virtual	~INuclearPotential ()

G4bool	hasPionPotential ()
	Do we have a pion potential?

virtual G4double	computePotentialEnergy (const Particle *const p) const =0

G4double	getFermiEnergy (const Particle *const p) const
	Return the Fermi energy for a particle.

G4double	getFermiEnergy (const ParticleType t) const
	Return the Fermi energy for a particle type.

G4double	getSeparationEnergy (const Particle *const p) const
	Return the separation energy for a particle.

G4double	getSeparationEnergy (const ParticleType t) const
	Return the separation energy for a particle type.

G4double	getFermiMomentum (const Particle *const p) const
	Return the Fermi momentum for a particle.

G4double	getFermiMomentum (const ParticleType t) const
	Return the Fermi momentum for a particle type.

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.

Protected Attributes inherited from G4INCL::NuclearPotential::INuclearPotential
const G4int	theA
	The mass number of the nucleus.

const G4int	theZ
	The charge number of the nucleus.

std::map< ParticleType, G4double >	fermiEnergy

std::map< ParticleType, G4double >	fermiMomentum

std::map< ParticleType, G4double >	separationEnergy

Detailed Description

Definition at line 58 of file G4INCLNuclearPotentialEnergyIsospinSmooth.hh.

Constructor & Destructor Documentation

◆ NuclearPotentialEnergyIsospinSmooth()

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

Definition at line 61 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

62 : NuclearPotentialIsospin(A,Z,aPionPotential)

63 {}

G4INCL::NuclearPotential::NuclearPotentialIsospin::NuclearPotentialIsospin

NuclearPotentialIsospin(const G4int A, const G4int Z, const G4bool pionPotential)

Definition: G4INCLNuclearPotentialIsospin.cc:58

◆ ~NuclearPotentialEnergyIsospinSmooth()

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

virtual

Definition at line 66 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

66{}

Member Function Documentation

◆ computePotentialEnergy()

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

virtual

Reimplemented from G4INCL::NuclearPotential::NuclearPotentialIsospin.

Definition at line 68 of file G4INCLNuclearPotentialEnergyIsospinSmooth.cc.

                                                                                                       {
 
      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: