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

---

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(), and G4INCL::Particle::isNucleon().

                                                                                                       {

      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

---