#include <G4INCLNuclearPotentialIsospin.hh>

Inheritance diagram for G4INCL::NuclearPotential::NuclearPotentialIsospin:

Public Member Functions
virtual G4double	computePotentialEnergy (const Particle *const p) const

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

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

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

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

virtual	~NuclearPotentialIsospin ()

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

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
std::map< ParticleType, G4double >	fermiEnergy

std::map< ParticleType, G4double >	fermiMomentum

std::map< ParticleType, G4double >	separationEnergy

const G4int	theA
	The mass number of the nucleus. More...

const G4int	theZ
	The charge number of the nucleus. More...

Private Member Functions
void	initialize ()

Private Attributes
const G4bool	pionPotential

G4double	vDeltaMinus

G4double	vDeltaPlus

G4double	vDeltaPlusPlus

G4double	vDeltaZero

G4double	vKMinus

G4double	vKPlus

G4double	vKZero

G4double	vKZeroBar

G4double	vLambda

G4double	vNeutron

G4double	vPiMinus

G4double	vPiPlus

G4double	vPiZero

G4double	vProton

G4double	vSigmaMinus

G4double	vSigmaPlus

G4double	vSigmaZero

Static Private Attributes
static const G4double	vKMinusDefault = 60.

static const G4double	vKPlusDefault = -25.

static const G4double	vPionDefault = 30.6

Detailed Description

Definition at line 56 of file G4INCLNuclearPotentialIsospin.hh.

Constructor & Destructor Documentation

◆ NuclearPotentialIsospin()

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

Definition at line 57 of file G4INCLNuclearPotentialIsospin.cc.

      : INuclearPotential(A, Z, aPionPotential)
    {
      initialize();
    }

References initialize().

◆ ~NuclearPotentialIsospin()

G4INCL::NuclearPotential::NuclearPotentialIsospin::~NuclearPotentialIsospin ( )

virtual

Definition at line 64 of file G4INCLNuclearPotentialIsospin.cc.

64{}

Member Function Documentation

◆ computeKaonPotentialEnergy()

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

inlineprotectedinherited

Compute the potential energy for the given kaon.

Definition at line 197 of file G4INCLINuclearPotential.hh.

                                                                            {
// assert(p->getType()==KPlus || p->getType()==KZero || p->getType()==KZeroBar || p->getType()==KMinus|| p->getType()==KShort|| p->getType()==KLong);
          if(pionPotential && !p->isOutOfWell()) { // if pionPotental false -> kaonPotential false
            switch( p->getType() ) {
              case KPlus:
                return vKPlus;
                break;
              case KZero:
                return vKZero;
                break;
              case KZeroBar:
                return vKZeroBar;
                break;
              case KShort:
              case KLong:
                 return 0.0; // Should never be in the nucleus
                 break;
               case KMinus:
                 return vKMinus;
                 break;
               default:
                 return 0.0;
                 break;
               }
            }
        else
          return 0.0;
        }

References G4INCL::Particle::getType(), G4INCL::Particle::isOutOfWell(), G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::NuclearPotential::INuclearPotential::pionPotential, G4INCL::NuclearPotential::INuclearPotential::vKMinus, G4INCL::NuclearPotential::INuclearPotential::vKPlus, G4INCL::NuclearPotential::INuclearPotential::vKZero, and G4INCL::NuclearPotential::INuclearPotential::vKZeroBar.

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

◆ computePionPotentialEnergy()

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

inlineprotectedinherited

Compute the potential energy for the given pion.

Definition at line 173 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;
        }

References G4INCL::Particle::getType(), G4INCL::Particle::isOutOfWell(), G4INCL::PiMinus, G4INCL::NuclearPotential::INuclearPotential::pionPotential, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::NuclearPotential::INuclearPotential::vPiMinus, G4INCL::NuclearPotential::INuclearPotential::vPiPlus, and G4INCL::NuclearPotential::INuclearPotential::vPiZero.

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

◆ computePionResonancePotentialEnergy()

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

inlineprotectedinherited

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

Definition at line 228 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;
        }

References G4INCL::Eta, G4INCL::EtaPrime, G4INCL::Particle::getType(), G4INCL::Particle::isOutOfWell(), G4INCL::Omega, G4INCL::Photon, and G4INCL::NuclearPotential::INuclearPotential::pionPotential.

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

◆ computePotentialEnergy()

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

virtual

Implements G4INCL::NuclearPotential::INuclearPotential.

Reimplemented in G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin, and G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth.

Definition at line 198 of file G4INCLNuclearPotentialIsospin.cc.

                                                                                           {
 
      switch( particle->getType() )
      {
        case Proton:
          return vProton;
          break;
        case Neutron:
          return vNeutron;
          break;
 
        case PiPlus:
        case PiZero:
        case PiMinus:
          return computePionPotentialEnergy(particle);
          break;
        
        case SigmaPlus:
          return vSigmaPlus;
          break;
        case SigmaZero:
          return vSigmaZero;
          break;
        case Lambda:
          return vLambda;
          break;
        case SigmaMinus:
          return vSigmaMinus;
          break;
 
        case Eta:
        case Omega:
        case EtaPrime:
          return computePionResonancePotentialEnergy(particle);
          break;
 
        case KPlus:
        case KZero:
        case KZeroBar:
        case KMinus:
        case KShort:
        case KLong:
          return computeKaonPotentialEnergy(particle);
          break;
 
        case Photon:
          return 0.0;
          break;
 
        case DeltaPlusPlus:
          return vDeltaPlusPlus;
          break;
        case DeltaPlus:
          return vDeltaPlus;
          break;
        case DeltaZero:
          return vDeltaZero;
          break;
        case DeltaMinus:
          return vDeltaMinus;
          break;
      case Composite:
    INCL_ERROR("No potential computed for particle of type Cluster.");
    return 0.0;
    break;
      case UnknownParticle:
    INCL_ERROR("Trying to compute potential energy for an unknown particle.");
    return 0.0;
    break;
      }
 
      INCL_ERROR("There is no potential for this type of particle.");
      return 0.0;
    }

References G4INCL::Composite, G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::computePionResonancePotentialEnergy(), G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Eta, G4INCL::EtaPrime, G4INCL::Particle::getType(), INCL_ERROR, G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Neutron, G4INCL::Omega, G4INCL::Photon, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Proton, G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, G4INCL::UnknownParticle, vDeltaMinus, vDeltaPlus, vDeltaPlusPlus, vDeltaZero, vLambda, vNeutron, vProton, vSigmaMinus, vSigmaPlus, and vSigmaZero.

Referenced by G4INCL::NuclearPotential::NuclearPotentialEnergyIsospin::computePotentialEnergy(), and G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy().

◆ getFermiEnergy() [1/2]

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

inlineinherited

Return the Fermi energy for a particle.

Parameters

p	pointer to a Particle

Returns: Fermi energy for that particle type

Definition at line 105 of file G4INCLINuclearPotential.hh.

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

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

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

◆ getFermiEnergy() [2/2]

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

inlineinherited

Return the Fermi energy for a particle type.

Parameters

t	particle type

Returns: Fermi energy for that particle type

Definition at line 116 of file G4INCLINuclearPotential.hh.

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

References G4INCL::NuclearPotential::INuclearPotential::fermiEnergy.

◆ getFermiMomentum() [1/2]

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

inlineinherited

Return the Fermi momentum for a particle.

Parameters

p	pointer to a Particle

Returns: Fermi momentum for that particle type

Definition at line 149 of file G4INCLINuclearPotential.hh.

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

References G4INCL::NuclearPotential::INuclearPotential::fermiMomentum, G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::Particle::getMass(), G4INCL::Particle::getType(), and G4INCL::Particle::isDelta().

Referenced by G4INCL::PauliStandard::getBlockingProbability(), G4INCL::Nucleus::getSurfaceRadius(), G4INCL::ParticleSampler::sampleOneParticleWithFuzzyRPCorrelation(), and G4INCL::ParticleSampler::sampleOneParticleWithRPCorrelation().

◆ getFermiMomentum() [2/2]

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

inlineinherited

Return the Fermi momentum for a particle type.

Parameters

t	particle type

Returns: Fermi momentum for that particle type

Definition at line 165 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;
        }

References G4INCL::NuclearPotential::INuclearPotential::fermiMomentum.

◆ getSeparationEnergy() [1/2]

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

inlineinherited

Return the separation energy for a particle.

Parameters

p	pointer to a Particle

Returns: separation energy for that particle type

Definition at line 127 of file G4INCLINuclearPotential.hh.

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

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

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

◆ getSeparationEnergy() [2/2]

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

inlineinherited

Return the separation energy for a particle type.

Parameters

t	particle type

Returns: separation energy for that particle type

Definition at line 138 of file G4INCLINuclearPotential.hh.

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

References G4INCL::NuclearPotential::INuclearPotential::separationEnergy.

◆ hasPionPotential()

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

inlineinherited

Do we have a pion potential?

Definition at line 96 of file G4INCLINuclearPotential.hh.

96{ return pionPotential; }

References G4INCL::NuclearPotential::INuclearPotential::pionPotential.

Referenced by G4INCL::Nucleus::decayInsideDeltas().

◆ initialize()

void G4INCL::NuclearPotential::NuclearPotentialIsospin::initialize ( )

private

Definition at line 66 of file G4INCLNuclearPotentialIsospin.cc.

                                             {
      const G4double ZOverA = ((G4double) theZ) / ((G4double) theA);
 
      const G4double mp = ParticleTable::getINCLMass(Proton);
      const G4double mn = ParticleTable::getINCLMass(Neutron);
      const G4double ml = ParticleTable::getINCLMass(Lambda);
 
      const G4double theFermiMomentum = ParticleTable::getFermiMomentum(theA,theZ);
 
      fermiMomentum[Proton] = theFermiMomentum * Math::pow13(2.*ZOverA);
      const G4double theProtonFermiEnergy = std::sqrt(fermiMomentum[Proton]*fermiMomentum[Proton] + mp*mp) - mp;
      fermiEnergy[Proton] = theProtonFermiEnergy;
      // Use separation energies from the ParticleTable
      const G4double theProtonSeparationEnergy = ParticleTable::getSeparationEnergy(Proton,theA,theZ);
      separationEnergy[Proton] = theProtonSeparationEnergy;
      vProton = theProtonFermiEnergy + theProtonSeparationEnergy;
 
      fermiMomentum[Neutron] = theFermiMomentum * Math::pow13(2.*(1.-ZOverA));
      const G4double theNeutronFermiEnergy = std::sqrt(fermiMomentum[Neutron]*fermiMomentum[Neutron] + mn*mn) - mn;
      fermiEnergy[Neutron] = theNeutronFermiEnergy;
      // Use separation energies from the ParticleTable
      const G4double theNeutronSeparationEnergy = ParticleTable::getSeparationEnergy(Neutron,theA,theZ);
      separationEnergy[Neutron] = theNeutronSeparationEnergy;
      vNeutron = theNeutronFermiEnergy + theNeutronSeparationEnergy;
 
      const G4double separationEnergyDeltaPlusPlus = 2.*theProtonSeparationEnergy - theNeutronSeparationEnergy;
      separationEnergy[DeltaPlusPlus] = separationEnergyDeltaPlusPlus;
      separationEnergy[DeltaPlus] = theProtonSeparationEnergy;
      separationEnergy[DeltaZero] = theNeutronSeparationEnergy;
      const G4double separationEnergyDeltaMinus = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy;
      separationEnergy[DeltaMinus] = separationEnergyDeltaMinus;
 
      const G4double tinyMargin = 1E-7;
      vDeltaPlus = vProton;
      vDeltaZero = vNeutron;
      vDeltaPlusPlus = std::max(separationEnergyDeltaPlusPlus + tinyMargin, 2.*vDeltaPlus - vDeltaZero);
      vDeltaMinus = std::max(separationEnergyDeltaMinus + tinyMargin, 2.*vDeltaZero - vDeltaPlus);
      
      vSigmaMinus = -16.; // Repulsive potential, from Eur. Phys.J.A. (2016) 52:21
      vSigmaZero = -16.;  // hypothesis: same potential for each sigma
      vSigmaPlus = -16.;
 
      vLambda = 30.;
      const G4double asy = (theA - 2.*theZ)/theA;
      // Jose Luis Rodriguez-Sanchez et al., Rapid Communication PRC 98, 021602 (2018)
      if      (asy > 0.236) vLambda = 40.91;
      else if (asy > 0.133) vLambda = 56.549 - 678.73*asy + 4905.35*asy*asy - 9789.1*asy*asy*asy;
 
      const G4double theLambdaSeparationEnergy = ParticleTable::getSeparationEnergy(Lambda,theA,theZ);
 
      separationEnergy[PiPlus] = theProtonSeparationEnergy - theNeutronSeparationEnergy;
      separationEnergy[PiZero] = 0.;
      separationEnergy[PiMinus] = theNeutronSeparationEnergy - theProtonSeparationEnergy;
 
      separationEnergy[Eta]      = 0.;
      separationEnergy[Omega]    = 0.;
      separationEnergy[EtaPrime] = 0.;
      separationEnergy[Photon]   = 0.;
      
      separationEnergy[Lambda]      = theLambdaSeparationEnergy;
      separationEnergy[SigmaPlus]   = theProtonSeparationEnergy + theLambdaSeparationEnergy - theNeutronSeparationEnergy;
      separationEnergy[SigmaZero]   = theLambdaSeparationEnergy;
      separationEnergy[SigmaMinus]  = theNeutronSeparationEnergy + theLambdaSeparationEnergy - theProtonSeparationEnergy;
 
      separationEnergy[KPlus]       = theProtonSeparationEnergy - theLambdaSeparationEnergy;
      separationEnergy[KZero]       = (theNeutronSeparationEnergy - theLambdaSeparationEnergy);
      separationEnergy[KZeroBar]    = (theLambdaSeparationEnergy - theNeutronSeparationEnergy);
      separationEnergy[KMinus]      = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy-theLambdaSeparationEnergy;
 
      separationEnergy[KShort]      = (theNeutronSeparationEnergy - theLambdaSeparationEnergy);
      separationEnergy[KLong]       = (theNeutronSeparationEnergy - theLambdaSeparationEnergy);
 
      fermiEnergy[DeltaPlusPlus] = vDeltaPlusPlus - separationEnergy[DeltaPlusPlus];
      fermiEnergy[DeltaPlus] = vDeltaPlus - separationEnergy[DeltaPlus];
      fermiEnergy[DeltaZero] = vDeltaZero - separationEnergy[DeltaZero];
      fermiEnergy[DeltaMinus] = vDeltaMinus - separationEnergy[DeltaMinus];
      
      fermiEnergy[Lambda] = vLambda - separationEnergy[Lambda];
      if (fermiEnergy[Lambda] <= 0.)
         fermiMomentum[Lambda]=0.;
      else
         fermiMomentum[Lambda]=std::sqrt(std::pow(fermiEnergy[Lambda]+ml,2.0)-ml*ml);
 
      fermiEnergy[SigmaPlus] = vSigmaPlus - separationEnergy[SigmaPlus];
      fermiEnergy[SigmaZero] = vSigmaZero - separationEnergy[SigmaZero];
      fermiEnergy[SigmaMinus] = vSigmaMinus - separationEnergy[SigmaMinus];
 
      INCL_DEBUG("Table of separation energies [MeV] for A=" << theA << ", Z=" << theZ << ":" << '\n'
            << "  proton:  " << separationEnergy[Proton] << '\n'
            << "  neutron: " << separationEnergy[Neutron] << '\n'
            << "  delta++: " << separationEnergy[DeltaPlusPlus] << '\n'
            << "  delta+:  " << separationEnergy[DeltaPlus] << '\n'
            << "  delta0:  " << separationEnergy[DeltaZero] << '\n'
            << "  delta-:  " << separationEnergy[DeltaMinus] << '\n'
            << "  pi+:     " << separationEnergy[PiPlus] << '\n'
            << "  pi0:     " << separationEnergy[PiZero] << '\n'
            << "  pi-:     " << separationEnergy[PiMinus] << '\n'
            << "  eta:     " << separationEnergy[Eta] << '\n'
            << "  omega:   " << separationEnergy[Omega] << '\n'
            << "  etaprime:" << separationEnergy[EtaPrime] << '\n'
            << "  photon:  " << separationEnergy[Photon] << '\n'
            << "  lambda:  " << separationEnergy[Lambda] << '\n'
            << "  sigmaplus:  " << separationEnergy[SigmaPlus] << '\n'
            << "  sigmazero:  " << separationEnergy[SigmaZero] << '\n'
            << "  sigmaminus:  " << separationEnergy[SigmaMinus] << '\n'
            << "  kplus:  " << separationEnergy[KPlus] << '\n'
            << "  kzero:  " << separationEnergy[KZero] << '\n'
            << "  kzerobar:  " << separationEnergy[KZeroBar] << '\n'
            << "  kminus:  " << separationEnergy[KMinus] << '\n'
            << "  kshort:  " << separationEnergy[KShort] << '\n'
            << "  klong:  " << separationEnergy[KLong] << '\n'
            );
 
      INCL_DEBUG("Table of Fermi energies [MeV] for A=" << theA << ", Z=" << theZ << ":" << '\n'
            << "  proton:  " << fermiEnergy[Proton] << '\n'
            << "  neutron: " << fermiEnergy[Neutron] << '\n'
            << "  delta++: " << fermiEnergy[DeltaPlusPlus] << '\n'
            << "  delta+:  " << fermiEnergy[DeltaPlus] << '\n'
            << "  delta0:  " << fermiEnergy[DeltaZero] << '\n'
            << "  delta-:  " << fermiEnergy[DeltaMinus] << '\n'
            << "  lambda:  " << fermiEnergy[Lambda] << '\n'
            << "  sigma+:  " << fermiEnergy[SigmaPlus] << '\n'
            << "  sigma0:  " << fermiEnergy[SigmaZero] << '\n'
            << "  sigma-:  " << fermiEnergy[SigmaMinus] << '\n'
            );
 
      INCL_DEBUG("Table of Fermi momenta [MeV/c] for A=" << theA << ", Z=" << theZ << ":" << '\n'
            << "  proton:  " << fermiMomentum[Proton] << '\n'
            << "  neutron: " << fermiMomentum[Neutron] << '\n'
            );
    }

References G4INCL::DeltaMinus, G4INCL::DeltaPlus, G4INCL::DeltaPlusPlus, G4INCL::DeltaZero, G4INCL::Eta, G4INCL::EtaPrime, G4INCL::NuclearPotential::INuclearPotential::fermiEnergy, G4INCL::NuclearPotential::INuclearPotential::fermiMomentum, G4INCL::ParticleTable::getFermiMomentum, G4INCL::ParticleTable::getINCLMass(), G4INCL::ParticleTable::getSeparationEnergy, INCL_DEBUG, G4INCL::KLong, G4INCL::KMinus, G4INCL::KPlus, G4INCL::KShort, G4INCL::KZero, G4INCL::KZeroBar, G4INCL::Lambda, G4INCL::Math::max(), G4INCL::Neutron, G4INCL::Omega, G4INCL::Photon, G4INCL::PiMinus, G4INCL::PiPlus, G4INCL::PiZero, G4INCL::Math::pow13(), G4INCL::Proton, G4INCL::NuclearPotential::INuclearPotential::separationEnergy, G4INCL::SigmaMinus, G4INCL::SigmaPlus, G4INCL::SigmaZero, G4INCL::NuclearPotential::INuclearPotential::theA, G4INCL::NuclearPotential::INuclearPotential::theZ, vDeltaMinus, vDeltaPlus, vDeltaPlusPlus, vDeltaZero, vLambda, vNeutron, vProton, vSigmaMinus, vSigmaPlus, and vSigmaZero.

Referenced by NuclearPotentialIsospin().

Field Documentation

◆ fermiEnergy

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

protectedinherited

Definition at line 269 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy(), G4INCL::NuclearPotential::NuclearPotentialConstant::initialize(), and initialize().

◆ fermiMomentum

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

protectedinherited

Definition at line 271 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::getFermiMomentum(), G4INCL::NuclearPotential::NuclearPotentialConstant::initialize(), and initialize().

◆ pionPotential

const G4bool G4INCL::NuclearPotential::INuclearPotential::pionPotential

privateinherited

Definition at line 261 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::computePionResonancePotentialEnergy(), G4INCL::NuclearPotential::INuclearPotential::hasPionPotential(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ separationEnergy

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

protectedinherited

Definition at line 273 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::getSeparationEnergy(), G4INCL::NuclearPotential::NuclearPotentialConstant::initialize(), and initialize().

◆ theA

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

protectedinherited

The mass number of the nucleus.

Definition at line 257 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::NuclearPotentialConstant::initialize(), initialize(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ theZ

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

protectedinherited

The charge number of the nucleus.

Definition at line 259 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::NuclearPotentialConstant::initialize(), initialize(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vDeltaMinus

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vDeltaMinus

private

Definition at line 66 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vDeltaPlus

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vDeltaPlus

private

Definition at line 66 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vDeltaPlusPlus

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vDeltaPlusPlus

private

Definition at line 66 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vDeltaZero

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vDeltaZero

private

Definition at line 66 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vKMinus

G4double G4INCL::NuclearPotential::INuclearPotential::vKMinus

privateinherited

Definition at line 264 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vKMinusDefault

const G4double G4INCL::NuclearPotential::INuclearPotential::vKMinusDefault = 60.

staticprivateinherited

Definition at line 266 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vKPlus

G4double G4INCL::NuclearPotential::INuclearPotential::vKPlus

privateinherited

Definition at line 264 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vKPlusDefault

const G4double G4INCL::NuclearPotential::INuclearPotential::vKPlusDefault = -25.

staticprivateinherited

Definition at line 265 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vKZero

G4double G4INCL::NuclearPotential::INuclearPotential::vKZero

privateinherited

Definition at line 264 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vKZeroBar

G4double G4INCL::NuclearPotential::INuclearPotential::vKZeroBar

privateinherited

Definition at line 264 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computeKaonPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vLambda

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vLambda

private

Definition at line 67 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vNeutron

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vNeutron

private

Definition at line 65 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vPiMinus

G4double G4INCL::NuclearPotential::INuclearPotential::vPiMinus

privateinherited

Definition at line 262 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vPionDefault

const G4double G4INCL::NuclearPotential::INuclearPotential::vPionDefault = 30.6

staticprivateinherited

Definition at line 263 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vPiPlus

G4double G4INCL::NuclearPotential::INuclearPotential::vPiPlus

privateinherited

Definition at line 262 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vPiZero

G4double G4INCL::NuclearPotential::INuclearPotential::vPiZero

privateinherited

Definition at line 262 of file G4INCLINuclearPotential.hh.

Referenced by G4INCL::NuclearPotential::INuclearPotential::computePionPotentialEnergy(), and G4INCL::NuclearPotential::INuclearPotential::INuclearPotential().

◆ vProton

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vProton

private

Definition at line 65 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vSigmaMinus

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vSigmaMinus

private

Definition at line 67 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vSigmaPlus

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vSigmaPlus

private

Definition at line 67 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

◆ vSigmaZero

G4double G4INCL::NuclearPotential::NuclearPotentialIsospin::vSigmaZero

private

Definition at line 67 of file G4INCLNuclearPotentialIsospin.hh.

Referenced by computePotentialEnergy(), and initialize().

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

source/processes/hadronic/models/inclxx/incl_physics/include/G4INCLNuclearPotentialIsospin.hh
source/processes/hadronic/models/inclxx/incl_physics/src/G4INCLNuclearPotentialIsospin.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ NuclearPotentialIsospin()

◆ ~NuclearPotentialIsospin()

Member Function Documentation

◆ computeKaonPotentialEnergy()

◆ computePionPotentialEnergy()

◆ computePionResonancePotentialEnergy()

◆ computePotentialEnergy()

◆ getFermiEnergy() [1/2]

◆ getFermiEnergy() [2/2]

◆ getFermiMomentum() [1/2]

◆ getFermiMomentum() [2/2]

◆ getSeparationEnergy() [1/2]

◆ getSeparationEnergy() [2/2]

◆ hasPionPotential()

◆ initialize()

Field Documentation

◆ fermiEnergy

◆ fermiMomentum

◆ pionPotential

◆ separationEnergy

◆ theA

◆ theZ

◆ vDeltaMinus

◆ vDeltaPlus

◆ vDeltaPlusPlus

◆ vDeltaZero

◆ vKMinus

◆ vKMinusDefault

◆ vKPlus

◆ vKPlusDefault

◆ vKZero

◆ vKZeroBar

◆ vLambda

◆ vNeutron

◆ vPiMinus

◆ vPionDefault

◆ vPiPlus

◆ vPiZero

◆ vProton

◆ vSigmaMinus

◆ vSigmaPlus

◆ vSigmaZero