#include <G4INCLINuclearPotential.hh>

Inheritance diagram for G4INCL::NuclearPotential::INuclearPotential:

Public Member Functions
	INuclearPotential (const G4int A, const G4int Z, const G4bool pionPot)

virtual	~INuclearPotential ()

G4bool	hasPionPotential () const
	Do we have a pion potential? More...

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

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...

Protected Member Functions
G4double	computePionPotentialEnergy (const Particle *const p) const
	Compute the potential energy for the given pion. More...

Protected Attributes
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 61 of file G4INCLINuclearPotential.hh.

Constructor & Destructor Documentation

G4INCL::NuclearPotential::INuclearPotential::INuclearPotential	(	const G4int	A,
		const G4int	Z,
		const G4bool	pionPot
	)

inline

Definition at line 63 of file G4INCLINuclearPotential.hh.

References G4INCL::PhysicalConstants::eSquared, G4INCL::Math::pow13(), theA, and theZ.

                                                                               :
           theA(A),
           theZ(Z),
           pionPotential(pionPot)
         {
           if(pionPotential) {
             const G4double ZOverA = ((G4double) theZ) / ((G4double) theA);
             // As in INCL4.6, use the r0*A^(1/3) formula to estimate vc
             const G4double r = 1.12*Math::pow13((G4double)theA);
 
             const G4double xsi = 1. - 2.*ZOverA;
             const G4double vc = 1.25*PhysicalConstants::eSquared*theZ/r;
             vPiPlus = vPionDefault + 71.*xsi - vc;
             vPiZero = vPionDefault;
             vPiMinus = vPionDefault - 71.*xsi + vc;
           } else {
             vPiPlus = 0.0;
             vPiZero = 0.0;
             vPiMinus = 0.0;
           }
         }

virtual G4INCL::NuclearPotential::INuclearPotential::~INuclearPotential ( )

inlinevirtual

Definition at line 85 of file G4INCLINuclearPotential.hh.

85 {}

Member Function Documentation

G4double G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy ( const Particle *const p ) const

inlineprotected

Compute the potential energy for the given pion.

Definition at line 165 of file G4INCLINuclearPotential.hh.

References G4INCL::Particle::getType(), G4INCL::Particle::isOutOfWell(), G4INCL::PiMinus, G4INCL::PiPlus, and G4INCL::PiZero.

Referenced by G4INCL::NuclearPotential::NuclearPotentialIsospin::computePotentialEnergy(), and G4INCL::NuclearPotential::NuclearPotentialConstant::computePotentialEnergy().

                                                                             {
 // assert(p->getType()==PiPlus || p->getType()==PiZero || p->getType()==PiMinus);
           if(pionPotential && !p->isOutOfWell()) {
             switch( p->getType() ) {
               case PiPlus:
                 return vPiPlus;
                 break;
               case PiZero:
                 return vPiZero;
                 break;
               case PiMinus:
                 return vPiMinus;
                 break;
               default: // Pion potential is defined and non-zero only for pions
                 return 0.0;
                 break;
             }
           }
           else
             return 0.0;
         }

virtual G4double G4INCL::NuclearPotential::INuclearPotential::computePotentialEnergy ( const Particle *const p ) const

pure virtual

Implemented in G4INCL::NuclearPotential::NuclearPotentialConstant, G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth, G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin, and G4INCL::NuclearPotential::NuclearPotentialIsospin.

Referenced by G4INCL::PionNucleonChannel::getFinalState(), and G4INCL::Nucleus::updatePotentialEnergy().

G4double G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy ( const Particle *const p ) const

inline

Return the Fermi energy for a particle.

Parameters

p	pointer to a Particle

Returns: Fermi energy for that particle type

Definition at line 97 of file G4INCLINuclearPotential.hh.

References fermiEnergy, and G4INCL::Particle::getType().

Referenced by G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy(), G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy(), G4INCL::SurfaceAvatar::getChannel(), getFermiMomentum(), G4INCL::ParticleEntryChannel::getFinalState(), and G4INCL::CDPP::processOneParticle().

                                                                        {
           std::map<ParticleType, G4double>::const_iterator i = fermiEnergy.find(p->getType());
 // assert(i!=fermiEnergy.end());
           return i->second;
         }

G4double G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy ( const ParticleType t ) const

inline

Return the Fermi energy for a particle type.

Parameters

t	particle type

Returns: Fermi energy for that particle type

Definition at line 108 of file G4INCLINuclearPotential.hh.

References fermiEnergy.

                                                                    {
           std::map<ParticleType, G4double>::const_iterator i = fermiEnergy.find(t);
 // assert(i!=fermiEnergy.end());
           return i->second;
         }

G4double G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum ( const Particle *const p ) const

inline

Return the Fermi momentum for a particle.

Parameters

p	pointer to a Particle

Returns: Fermi momentum for that particle type

Definition at line 141 of file G4INCLINuclearPotential.hh.

References fermiMomentum, getFermiEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getType(), G4INCL::Particle::isDelta(), and python.hepunit::m.

Referenced by G4INCL::PauliStandard::getBlockingProbability(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::Nucleus::getSurfaceRadius(), G4INCL::PauliGlobal::isBlocked(), and G4INCL::PauliStrict::isBlocked().

                                                                          {
           if(p->isDelta()) {
             const G4double Tf = getFermiEnergy(p), m = p->getMass();
             return std::sqrt(Tf*(Tf+2.*m));
           } else {
             std::map<ParticleType, G4double>::const_iterator i = fermiMomentum.find(p->getType());
 // assert(i!=fermiMomentum.end());
             return i->second;
           }
         }

G4double G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum ( const ParticleType t ) const

inline

Return the Fermi momentum for a particle type.

Parameters

t	particle type

Returns: Fermi momentum for that particle type

Definition at line 157 of file G4INCLINuclearPotential.hh.

References fermiMomentum.

                                                                      {
 // assert(t!=DeltaPlusPlus && t!=DeltaPlus && t!=DeltaZero && t!=DeltaMinus);
           std::map<ParticleType, G4double>::const_iterator i = fermiMomentum.find(t);
           return i->second;
         }

G4double G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy ( const Particle *const p ) const

inline

Return the separation energy for a particle.

Parameters

p	pointer to a Particle

Returns: separation energy for that particle type

Definition at line 119 of file G4INCLINuclearPotential.hh.

References G4INCL::Particle::getType(), and separationEnergy.

Referenced by G4INCL::Nucleus::computeSeparationEnergyBalance(), G4INCL::KinematicsUtils::getLocalEnergy(), G4INCL::Nucleus::Nucleus(), and G4INCL::CDPP::processOneParticle().

                                                                             {
           std::map<ParticleType, G4double>::const_iterator i = separationEnergy.find(p->getType());
 // assert(i!=separationEnergy.end());
           return i->second;
         }

G4double G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy ( const ParticleType t ) const

inline

Return the separation energy for a particle type.

Parameters

t	particle type

Returns: separation energy for that particle type

Definition at line 130 of file G4INCLINuclearPotential.hh.

References separationEnergy.

                                                                         {
           std::map<ParticleType, G4double>::const_iterator i = separationEnergy.find(t);
 // assert(i!=separationEnergy.end());
           return i->second;
         }