G4INCL::NuclearPotential::INuclearPotential Class Reference

#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...
G4double	computePionResonancePotentialEnergy (const Particle *const p) const
	Compute the potential energy for the given pion resonances (Eta, Omega and EtaPrime and Gamma also). 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 62 of file G4INCLINuclearPotential.hh.

Constructor & Destructor Documentation

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

inline

Definition at line 64 of file G4INCLINuclearPotential.hh.

                                                                              :
          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;
          }
        }

Here is the call graph for this function:

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

inlinevirtual

Definition at line 86 of file G4INCLINuclearPotential.hh.

{}

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 166 of file G4INCLINuclearPotential.hh.

                                                                            {
// 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;
        }

Here is the call graph for this function:

Here is the caller graph for this function:

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

inlineprotected

Compute the potential energy for the given pion resonances (Eta, Omega and EtaPrime and Gamma also).

Definition at line 190 of file G4INCLINuclearPotential.hh.

                                                                                     {
// assert(p->getType()==Eta || p->getType()==Omega || p->getType()==EtaPrime || p->getType()==Photon);
          if(pionPotential && !p->isOutOfWell()) {
            switch( p->getType() ) {
              case Eta:
//jcd             return vPiZero;
//jcd              return vPiZero*1.5;
                return 0.0; // (JCD: seems to give better results)
                break;
              case Omega:
                return 15.0; // S.Friedrich et al., Physics Letters B736(2014)26-32. (V. Metag in Hyperfine Interact (2015) 234:25-31 gives 29 MeV)
                break;
              case EtaPrime:
                return 37.0; // V. Metag in Hyperfine Interact (2015) 234:25-31
                break;
              case Photon:
                return 0.0;
                break;
              default: 
                return 0.0;
                break;
            }
          }
          else
            return 0.0;
        }

Here is the call graph for this function:

Here is the caller graph for this function:

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.

Here is the caller graph for this function:

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 98 of file G4INCLINuclearPotential.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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 109 of file G4INCLINuclearPotential.hh.

                                                                   {
          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 142 of file G4INCLINuclearPotential.hh.

                                                                         {
          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;
          }
        }

Here is the call graph for this function:

Here is the caller graph for this function:

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 158 of file G4INCLINuclearPotential.hh.

                                                                     {
// 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 120 of file G4INCLINuclearPotential.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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 131 of file G4INCLINuclearPotential.hh.

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

G4bool G4INCL::NuclearPotential::INuclearPotential::hasPionPotential ( ) const

inline

Do we have a pion potential?

Definition at line 89 of file G4INCLINuclearPotential.hh.

{ return pionPotential; }

Here is the caller graph for this function:

Member Data Documentation

std::map<ParticleType,G4double> G4INCL::NuclearPotential::INuclearPotential::fermiEnergy

protected

Definition at line 228 of file G4INCLINuclearPotential.hh.

std::map<ParticleType,G4double> G4INCL::NuclearPotential::INuclearPotential::fermiMomentum

protected

Definition at line 230 of file G4INCLINuclearPotential.hh.

std::map<ParticleType,G4double> G4INCL::NuclearPotential::INuclearPotential::separationEnergy

protected

Definition at line 232 of file G4INCLINuclearPotential.hh.

const G4int G4INCL::NuclearPotential::INuclearPotential::theA

protected

The mass number of the nucleus.

Definition at line 219 of file G4INCLINuclearPotential.hh.

const G4int G4INCL::NuclearPotential::INuclearPotential::theZ

protected

The charge number of the nucleus.

Definition at line 221 of file G4INCLINuclearPotential.hh.

---

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

• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLINuclearPotential.hh
• geant4.10.03.p01/source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLINuclearPotential.cc