#include <G4INCLParticleTable.hh>

Public Types
enum	ClusterDecayType { StableCluster , NeutronDecay , ProtonDecay , AlphaDecay , TwoProtonDecay , TwoNeutronDecay , ProtonUnbound , NeutronUnbound }

typedef G4double(*	NuclearMassFn) (const G4int, const G4int)

typedef G4double(*	ParticleMassFn) (const ParticleType)

typedef G4double(*	SeparationEnergyFn) (const ParticleType, const G4int, const G4int)

Static Public Member Functions
static void	initialize (Config const *const theConfig=0)
	Initialize the particle table.

static G4int	getIsospin (const ParticleType t)
	Get the isospin of a particle.

static std::string	getName (const ParticleType t)
	Get the native INCL name of the particle.

static std::string	getShortName (const ParticleType t)
	Get the short INCL name of the particle.

static std::string	getName (const ParticleSpecies s)
	Get the native INCL name of the particle.

static std::string	getShortName (const ParticleSpecies s)
	Get the short INCL name of the particle.

static std::string	getName (const G4int A, const G4int Z)
	Get the native INCL name of the ion.

static std::string	getShortName (const G4int A, const G4int Z)
	Get the short INCL name of the ion.

static G4double	getINCLMass (const G4int A, const G4int Z)
	Get INCL nuclear mass (in MeV/c^2)

static G4double	getINCLMass (const ParticleType t)
	Get INCL particle mass (in MeV/c^2)

static G4double	getRealMass (const G4INCL::ParticleType t)
	Get particle mass (in MeV/c^2)

static G4double	getRealMass (const G4int A, const G4int Z)
	Get nuclear mass (in MeV/c^2)

static G4double	getTableQValue (const G4int A1, const G4int Z1, const G4int A2, const G4int Z2)
	Get Q-value (in MeV/c^2)

static G4double	getTableQValue (const G4int A1, const G4int Z1, const G4int A2, const G4int Z2, const G4int A3, const G4int Z3)
	Get Q-value (in MeV/c^2)

static G4double	getTableSpeciesMass (const ParticleSpecies &p)

static G4int	getMassNumber (const ParticleType t)
	Get mass number from particle type.

static G4int	getChargeNumber (const ParticleType t)
	Get charge number from particle type.

static G4double	getNuclearRadius (const G4int A, const G4int Z)

static G4double	getRadiusParameter (const G4int A, const G4int Z)

static G4double	getMaximumNuclearRadius (const G4int A, const G4int Z)

static G4double	getSurfaceDiffuseness (const G4int A, const G4int Z)

static G4double	getMomentumRMS (const G4int A, const G4int Z)
	Return the RMS of the momentum distribution (light clusters)

static G4double	getSeparationEnergyINCL (const ParticleType t, const G4int, const G4int)
	Return INCL's default separation energy.

static G4double	getSeparationEnergyReal (const ParticleType t, const G4int A, const G4int Z)
	Return the real separation energy.

static G4double	getSeparationEnergyRealForLight (const ParticleType t, const G4int A, const G4int Z)
	Return the real separation energy only for light nuclei.

static G4double	getProtonSeparationEnergy ()
	Getter for protonSeparationEnergy.

static G4double	getNeutronSeparationEnergy ()
	Getter for neutronSeparationEnergy.

static void	setProtonSeparationEnergy (const G4double s)
	Setter for protonSeparationEnergy.

static void	setNeutronSeparationEnergy (const G4double s)
	Setter for protonSeparationEnergy.

static std::string	getElementName (const G4int Z)
	Get the name of the element from the atomic number.

static std::string	getIUPACElementName (const G4int Z)
	Get the name of an unnamed element from the IUPAC convention.

static G4int	parseIUPACElement (std::string const &pS)
	Parse a IUPAC element name.

static IsotopicDistribution const &	getNaturalIsotopicDistribution (const G4int Z)

static G4int	drawRandomNaturalIsotope (const G4int Z)

Static Public Attributes
static NuclearMassFn	getTableMass
	Static pointer to the mass function for nuclei.

static ParticleMassFn	getTableParticleMass
	Static pointer to the mass function for particles.

static SeparationEnergyFn	getSeparationEnergy
	Static pointer to the separation-energy function.

static const G4int	elementTableSize = 113

static const G4double	effectiveNucleonMass = 938.2796

static const G4double	effectiveNucleonMass2 = 8.8036860777616e5

static const G4double	effectiveDeltaMass = 1232.0

static const G4double	effectivePionMass = 138.0

static const G4double	effectiveDeltaDecayThreshold

static const G4int	maxClusterMass = 12

static const G4int	maxClusterCharge = 8

static const G4int	clusterTableZSize = ParticleTable::maxClusterCharge+1

static const G4int	clusterTableASize = ParticleTable::maxClusterMass+1

static const G4double	clusterPosFact [maxClusterMass+1]

static const G4double	clusterPosFact2 [maxClusterMass+1]

static const G4int	clusterZMin [maxClusterMass+1] = {0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3}

static const G4int	clusterZMax [maxClusterMass+1] = {0, 0, 1, 2, 3, 3, 5, 5, 6, 6, 7, 7, 8}

static const G4double	clusterPhaseSpaceCut [maxClusterMass+1]

static G4IonTable *	theG4IonTable

static const ClusterDecayType	clusterDecayMode [clusterTableZSize][clusterTableASize]

static const G4double	eSquared = 1.439964
	Coulomb conversion factor, in MeV*fm.

Protected Member Functions
	ParticleTable ()

	~ParticleTable ()

Detailed Description

Definition at line 60 of file G4INCLParticleTable.hh.

Member Typedef Documentation

◆ NuclearMassFn

typedef G4double(* G4INCL::ParticleTable::NuclearMassFn) (const G4int, const G4int)

Definition at line 151 of file G4INCLParticleTable.hh.

◆ ParticleMassFn

typedef G4double(* G4INCL::ParticleTable::ParticleMassFn) (const ParticleType)

Definition at line 152 of file G4INCLParticleTable.hh.

◆ SeparationEnergyFn

typedef G4double(* G4INCL::ParticleTable::SeparationEnergyFn) (const ParticleType, const G4int, const G4int)

Definition at line 164 of file G4INCLParticleTable.hh.

Member Enumeration Documentation

◆ ClusterDecayType

enum G4INCL::ParticleTable::ClusterDecayType

Enumerator
StableCluster
NeutronDecay
ProtonDecay
AlphaDecay
TwoProtonDecay
TwoNeutronDecay
ProtonUnbound
NeutronUnbound

Definition at line 320 of file G4INCLParticleTable.hh.

                          {
      StableCluster,
      NeutronDecay,
      ProtonDecay,
      AlphaDecay,
      TwoProtonDecay,
      TwoNeutronDecay,
      ProtonUnbound,
      NeutronUnbound
    };

Constructor & Destructor Documentation

◆ ParticleTable()

G4INCL::ParticleTable::ParticleTable ( )

inlineprotected

Definition at line 347 of file G4INCLParticleTable.hh.

347{};

◆ ~ParticleTable()

G4INCL::ParticleTable::~ParticleTable ( )

inlineprotected

Definition at line 348 of file G4INCLParticleTable.hh.

348{};

Member Function Documentation

◆ drawRandomNaturalIsotope()

static G4int G4INCL::ParticleTable::drawRandomNaturalIsotope ( const G4int Z )

inlinestatic

Definition at line 342 of file G4INCLParticleTable.hh.

                                                         {
      return getNaturalIsotopicDistributions()->drawRandomIsotope(Z);
    }

Referenced by G4INCL::INCL::prepareReaction().

◆ getChargeNumber()

static G4int G4INCL::ParticleTable::getChargeNumber ( const ParticleType t )

inlinestatic

Get charge number from particle type.

Definition at line 192 of file G4INCLParticleTable.hh.

                                                       {
      switch(t) {
        case DeltaPlusPlus:
          return 2;
          break;
        case Proton:
        case DeltaPlus:
        case PiPlus:
          return 1;
          break;
        case Neutron:
        case DeltaZero:
        case PiZero:
          return 0;
          break;
        case DeltaMinus:
        case PiMinus:
          return -1;
          break;
        default:
/*          FATAL("Can't determine charge number for particle type " << t << std::endl);
          std::abort();*/
          return 0;
          break;
      }
    }

◆ getElementName()

std::string G4INCL::ParticleTable::getElementName ( const G4int Z )

static

Get the name of the element from the atomic number.

Definition at line 624 of file G4INCLParticleTable.cc.

                                                       {
    if(Z<1) {
      WARN("getElementName called with Z<1" << std::endl);
      return elementTable[0];
    } else if(Z<elementTableSize)
      return elementTable[Z];
    else
      return getIUPACElementName(Z);
  }

Referenced by getName(), and getShortName().

◆ getINCLMass() [1/2]

G4double G4INCL::ParticleTable::getINCLMass	(	const G4int	A,
		const G4int	Z
	)

static

Get INCL nuclear mass (in MeV/c^2)

Definition at line 518 of file G4INCLParticleTable.cc.

                                                                  {
// assert(A>=0);
    // For nuclei with Z<0 or Z>A, assume that the exotic charge state is due to pions
    if(Z<0)
      return A*neutronMass - Z*getINCLMass(PiMinus);
    else if(Z>A)
      return A*protonMass + (A-Z)*getINCLMass(PiPlus);
    else if(A>1)
      return Z*(protonMass - protonSeparationEnergy) + (A-Z)*(neutronMass - neutronSeparationEnergy);
    else if(A==1 && Z==0)
      return getINCLMass(Neutron);
    else if(A==1 && Z==1)
      return getINCLMass(Proton);
    else
      return 0.;
  }

Referenced by G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::DeltaDecayChannel::getFinalState(), G4INCL::Particle::getINCLMass(), getINCLMass(), G4INCL::Particle::getTransferQValueCorrection(), initialize(), and G4INCL::CDPP::isBlocked().

◆ getINCLMass() [2/2]

G4double G4INCL::ParticleTable::getINCLMass ( const ParticleType t )

static

Get INCL particle mass (in MeV/c^2)

Definition at line 451 of file G4INCLParticleTable.cc.

                                                           {
    if(pt == Proton) {
      return protonMass;
    } else if(pt == Neutron) {
      return neutronMass;
    } else if(pt == PiPlus) {
      return piPlusMass;
    } else if(pt == PiMinus) {
      return piMinusMass;
    } else if(pt == PiZero) {
      return piZeroMass;
    } else {
      ERROR("ParticleTable::getMass : Unknown particle type." << std::endl);
      return 0.0;
    }
  }

◆ getIsospin()

G4int G4INCL::ParticleTable::getIsospin ( const ParticleType t )

static

Get the isospin of a particle.

Definition at line 347 of file G4INCLParticleTable.cc.

                                                      {
    // Actually this is the 3rd component of isospin (I_z) multiplied by 2!
    if(t == Proton) {
      return 1;
    } else if(t == Neutron) {
      return -1;
    } else if(t == PiPlus) {
      return 2;
    } else if(t == PiMinus) {
      return -2;
    } else if(t == PiZero) {
      return 0;
    } else if(t == DeltaPlusPlus) {
      return 3;
    } else if(t == DeltaPlus) {
      return 1;
    } else if(t == DeltaZero) {
      return -1;
    } else if(t == DeltaMinus) {
      return -3;
    }
 
    ERROR("Requested isospin of an unknown particle!");
    return -10; // Unknown
  }

Referenced by G4INCL::CrossSections::deltaProduction(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), G4INCL::Nucleus::insertParticle(), G4INCL::CrossSections::pionNucleon(), and G4INCL::CrossSections::recombination().

◆ getIUPACElementName()

std::string G4INCL::ParticleTable::getIUPACElementName ( const G4int Z )

static

Get the name of an unnamed element from the IUPAC convention.

Definition at line 678 of file G4INCLParticleTable.cc.

                                                            {
    std::stringstream elementStream;
    elementStream << Z;
    std::string elementName = elementStream.str();
    std::transform(elementName.begin(), elementName.end(), elementName.begin(), ParticleTable::intToIUPAC);
    elementName[0] = std::toupper(elementName.at(0));
    return elementName;
  }

Referenced by getElementName().

◆ getMassNumber()

static G4int G4INCL::ParticleTable::getMassNumber ( const ParticleType t )

inlinestatic

Get mass number from particle type.

Definition at line 168 of file G4INCLParticleTable.hh.

                                                     {
      switch(t) {
        case Proton:
        case Neutron:
        case DeltaPlusPlus:
        case DeltaPlus:
        case DeltaZero:
        case DeltaMinus:
          return 1;
          break;
        case PiPlus:
        case PiMinus:
        case PiZero:
          return 0;
          break;
        default:
/*          FATAL("Can't determine mass number for particle type " << t << std::endl);
          std::abort();*/
          return 0;
          break;
      }
    }

◆ getMaximumNuclearRadius()

G4double G4INCL::ParticleTable::getMaximumNuclearRadius	(	const G4int	A,
		const G4int	Z
	)

static

Definition at line 589 of file G4INCLParticleTable.cc.

                                                                              {
    const G4double XFOISA = 8.0;
    if(A >= 19) {
      return getNuclearRadius(A,Z) + XFOISA * getSurfaceDiffuseness(A,Z);
    } else if(A < 19 && A >= 6) {
      return 5.5 + 0.3 * (G4double(A) - 6.0)/12.0;
    } else if(A >= 2) {
      return ParticleTable::getNuclearRadius(A, Z) + 4.5;
    } else {
      ERROR("ParticleTable::getMaximumNuclearRadius : No maximum radius for nucleus A = " << A << " Z = " << Z << std::endl);
      return 0.0;
    }
  }

Referenced by G4INCL::NuclearDensityFactory::createRCDFTable(), and G4INCL::NuclearDensityFactory::createRPCorrelationTable().

◆ getMomentumRMS()

static G4double G4INCL::ParticleTable::getMomentumRMS	(	const G4int	A,
		const G4int	Z
	)

inlinestatic

Return the RMS of the momentum distribution (light clusters)

Definition at line 225 of file G4INCLParticleTable.hh.

                                                                 {
// assert(Z>=0 && A>=0 && Z<=A);
      if(Z<clusterTableZSize && A<clusterTableASize)
        return momentumRMS[Z][A];
      else
        return Math::sqrtThreeFifths * PhysicalConstants::Pf;
    }

Referenced by G4INCL::NuclearDensityFactory::createPCDFTable().

◆ getName() [1/3]

std::string G4INCL::ParticleTable::getName	(	const G4int	A,
		const G4int	Z
	)

static

Get the native INCL name of the ion.

Definition at line 387 of file G4INCLParticleTable.cc.

                                                               {
    std::stringstream stream;
    stream << getElementName(Z) << "-" << A;
    return stream.str();
  }

◆ getName() [2/3]

std::string G4INCL::ParticleTable::getName ( const ParticleSpecies s )

static

Get the native INCL name of the particle.

Definition at line 380 of file G4INCLParticleTable.cc.

                                                          {
    if(s.theType==Composite)
      return getName(s.theA,s.theZ);
    else
      return getName(s.theType);
  }

◆ getName() [3/3]

std::string G4INCL::ParticleTable::getName ( const ParticleType t )

static

Get the native INCL name of the particle.

Definition at line 401 of file G4INCLParticleTable.cc.

                                                       {
    if(p == G4INCL::Proton) {
      return std::string("proton");
    } else if(p == G4INCL::Neutron) {
      return std::string("neutron");
    } else if(p == G4INCL::DeltaPlusPlus) {
      return std::string("delta++");
    } else if(p == G4INCL::DeltaPlus) {
      return std::string("delta+");
    } else if(p == G4INCL::DeltaZero) {
      return std::string("delta0");
    } else if(p == G4INCL::DeltaMinus) {
      return std::string("delta-");
    } else if(p == G4INCL::PiPlus) {
      return std::string("pi+");
    } else if(p == G4INCL::PiZero) {
      return std::string("pi0");
    } else if(p == G4INCL::PiMinus) {
      return std::string("pi-");
    } else if(p == G4INCL::Composite) {
      return std::string("composite");
    }
    return std::string("unknown");
  }

Referenced by G4INCL::Particle::dump(), getName(), G4INCL::Cluster::print(), G4INCL::Particle::print(), and G4INCL::Config::summary().

◆ getNaturalIsotopicDistribution()

static IsotopicDistribution const & G4INCL::ParticleTable::getNaturalIsotopicDistribution ( const G4int Z )

inlinestatic

Definition at line 338 of file G4INCLParticleTable.hh.

                                                                                     {
      return getNaturalIsotopicDistributions()->getIsotopicDistribution(Z);
    }

◆ getNeutronSeparationEnergy()

static G4double G4INCL::ParticleTable::getNeutronSeparationEnergy ( )

inlinestatic

Getter for neutronSeparationEnergy.

Definition at line 271 of file G4INCLParticleTable.hh.

271{ return neutronSeparationEnergy; }

◆ getNuclearRadius()

G4double G4INCL::ParticleTable::getNuclearRadius	(	const G4int	A,
		const G4int	Z
	)

static

Definition at line 535 of file G4INCLParticleTable.cc.

                                                                       {
// assert(A>0 && Z>=0 && Z<=A);
    if(A >= 19 || (A < 6 && A >= 2)) {
      // For large (Woods-Saxon or Modified Harmonic Oscillator) or small
      // (Gaussian) nuclei, the radius parameter is just the nuclear radius
      return getRadiusParameter(A,Z);
    } else if(A < clusterTableASize && Z < clusterTableZSize && A >= 6) {
      const G4double thisRMS = positionRMS[Z][A];
      if(thisRMS>0.0)
        return thisRMS;
      else {
        ERROR("ParticleTable::getRadiusParameter : Radius for nucleus A = " << A << " Z = " << Z << " is ==0.0" << std::endl);
        return 0.0;
      }
    } else if(A < 19) {
      const G4double theRadiusParameter = getRadiusParameter(A, Z);
      const G4double theDiffusenessParameter = getSurfaceDiffuseness(A, Z);
      // The formula yields the nuclear RMS radius based on the parameters of
      // the nuclear-density function
      return 1.581*theDiffusenessParameter*
        (2.+5.*theRadiusParameter)/(2.+3.*theRadiusParameter);
    } else {
      ERROR("ParticleTable::getNuclearRadius: No radius for nucleus A = " << A << " Z = " << Z << std::endl);
      return 0.0;
    }
  }

Referenced by getMaximumNuclearRadius(), G4INCL::NuclearDensity::getTransmissionRadius(), G4INCL::CoulombNonRelativistic::maxImpactParameter(), G4INCL::CoulombNone::maxImpactParameter(), and G4INCL::StandardPropagationModel::shootComposite().

◆ getProtonSeparationEnergy()

static G4double G4INCL::ParticleTable::getProtonSeparationEnergy ( )

inlinestatic

Getter for protonSeparationEnergy.

Definition at line 268 of file G4INCLParticleTable.hh.

268{ return protonSeparationEnergy; }

◆ getRadiusParameter()

G4double G4INCL::ParticleTable::getRadiusParameter	(	const G4int	A,
		const G4int	Z
	)

static

Definition at line 562 of file G4INCLParticleTable.cc.

                                                                         {
// assert(A>0 && Z>=0 && Z<=A);
    if(A >= 28) {
      // phenomenological radius fit
      return (2.745e-4 * A + 1.063) * std::pow(A, 1.0/3.0);
    } else if(A < 6 && A >= 2) {
      if(Z<clusterTableZSize) {
        const G4double thisRMS = positionRMS[Z][A];
        if(thisRMS>0.0)
          return thisRMS;
        else {
          ERROR("ParticleTable::getRadiusParameter : Radius for nucleus A = " << A << " Z = " << Z << " is ==0.0" << std::endl);
          return 0.0;
        }
      } else {
        ERROR("ParticleTable::getRadiusParameter : No radius for nucleus A = " << A << " Z = " << Z << std::endl);
        return 0.0;
      }
    } else if(A < 28 && A >= 6) {
      return mediumRadius[A-1];
      //      return 1.581*mediumDiffuseness[A-1]*(2.+5.*mediumRadius[A-1])/(2.+3.*mediumRadius[A-1]);
    } else {
      ERROR("ParticleTable::getRadiusParameter: No radius for nucleus A = " << A << " Z = " << Z << std::endl);
      return 0.0;
    }
  }

Referenced by G4INCL::NuclearDensityFactory::createRCDFTable(), G4INCL::NuclearDensityFactory::createRPCorrelationTable(), and getNuclearRadius().

◆ getRealMass() [1/2]

G4double G4INCL::ParticleTable::getRealMass ( const G4INCL::ParticleType t )

static

Get particle mass (in MeV/c^2)

Definition at line 468 of file G4INCLParticleTable.cc.

                                                          {
    switch(t) {
      case Proton:
        return theRealProtonMass;
        break;
      case Neutron:
        return theRealNeutronMass;
        break;
      case PiPlus:
      case PiMinus:
        return theRealChargedPiMass;
        break;
      case PiZero:
        return theRealPiZeroMass;
        break;
      default:
        ERROR("Particle::getRealMass : Unknown particle type." << std::endl);
        return 0.0;
        break;
    }
  }

Referenced by G4INCL::Particle::getRealMass(), getRealMass(), and initialize().

◆ getRealMass() [2/2]

G4double G4INCL::ParticleTable::getRealMass	(	const G4int	A,
		const G4int	Z
	)

static

Get nuclear mass (in MeV/c^2)

Definition at line 490 of file G4INCLParticleTable.cc.

                                                                  {
// assert(A>=0);
    // For nuclei with Z<0 or Z>A, assume that the exotic charge state is due to pions
    if(Z<0)
      return A*neutronMass - Z*getRealMass(PiMinus);
    else if(Z>A)
      return A*protonMass + (A-Z)*getRealMass(PiPlus);
    else if(Z==0)
      return A*getRealMass(Neutron);
    else if(A==Z)
      return A*getRealMass(Proton);
    else if(A>1) {
#ifndef INCLXX_IN_GEANT4_MODE
      if(ParticleTable::hasMassTable(A,Z))
        return ParticleTable::massTable.at(Z).at(A);
      else {
        DEBUG("Real mass unavailable for isotope A=" << A << ", Z=" << Z
            << ", using Weizsaecker's formula"
            << std::endl);
        return getWeizsaeckerMass(A,Z);
      }
#else
      return theG4IonTable->GetNucleusMass(Z,A) / MeV;
#endif
    } else
      return 0.;
  }

◆ getSeparationEnergyINCL()

static G4double G4INCL::ParticleTable::getSeparationEnergyINCL	(	const ParticleType	t,
		const	G4int,
		const	G4int
	)

inlinestatic

Return INCL's default separation energy.

Definition at line 234 of file G4INCLParticleTable.hh.

                                                                                                {
      if(t==Proton)
        return theINCLProtonSeparationEnergy;
      else if(t==Neutron)
        return theINCLNeutronSeparationEnergy;
      else {
        ERROR("ParticleTable::getSeparationEnergyINCL : Unknown particle type." << std::endl);
        return 0.0;
      }
    }

Referenced by getSeparationEnergyRealForLight(), and initialize().

◆ getSeparationEnergyReal()

static G4double G4INCL::ParticleTable::getSeparationEnergyReal	(	const ParticleType	t,
		const G4int	A,
		const G4int	Z
	)

inlinestatic

Return the real separation energy.

Definition at line 246 of file G4INCLParticleTable.hh.

                                                                                                {
      // Real separation energies for all nuclei
      if(t==Proton)
        return (*getTableParticleMass)(Proton) + (*getTableMass)(A-1,Z-1) - (*getTableMass)(A,Z);
      else if(t==Neutron)
        return (*getTableParticleMass)(Neutron) + (*getTableMass)(A-1,Z) - (*getTableMass)(A,Z);
      else {
        ERROR("ParticleTable::getSeparationEnergyReal : Unknown particle type." << std::endl);
        return 0.0;
      }
    }

Referenced by getSeparationEnergyRealForLight(), and initialize().

◆ getSeparationEnergyRealForLight()

static G4double G4INCL::ParticleTable::getSeparationEnergyRealForLight	(	const ParticleType	t,
		const G4int	A,
		const G4int	Z
	)

inlinestatic

Return the real separation energy only for light nuclei.

Definition at line 259 of file G4INCLParticleTable.hh.

                                                                                                        {
      // Real separation energies for light nuclei, fixed values for heavy nuclei
      if(Z<clusterTableZSize && A<clusterTableASize)
        return getSeparationEnergyReal(t, A, Z);
      else
        return getSeparationEnergyINCL(t, A, Z);
    }

Referenced by initialize().

◆ getShortName() [1/3]

std::string G4INCL::ParticleTable::getShortName	(	const G4int	A,
		const G4int	Z
	)

static

Get the short INCL name of the ion.

Definition at line 393 of file G4INCLParticleTable.cc.

                                                                    {
    std::stringstream stream;
    stream << getElementName(Z);
    if(A>0)
      stream << A;
    return stream.str();
  }

◆ getShortName() [2/3]

std::string G4INCL::ParticleTable::getShortName ( const ParticleSpecies s )

static

Get the short INCL name of the particle.

Definition at line 373 of file G4INCLParticleTable.cc.

                                                               {
    if(s.theType==Composite)
      return getShortName(s.theA,s.theZ);
    else
      return getShortName(s.theType);
  }

◆ getShortName() [3/3]

std::string G4INCL::ParticleTable::getShortName ( const ParticleType t )

static

Get the short INCL name of the particle.

Definition at line 426 of file G4INCLParticleTable.cc.

                                                            {
    if(p == G4INCL::Proton) {
      return std::string("p");
    } else if(p == G4INCL::Neutron) {
      return std::string("n");
    } else if(p == G4INCL::DeltaPlusPlus) {
      return std::string("d++");
    } else if(p == G4INCL::DeltaPlus) {
      return std::string("d+");
    } else if(p == G4INCL::DeltaZero) {
      return std::string("d0");
    } else if(p == G4INCL::DeltaMinus) {
      return std::string("d-");
    } else if(p == G4INCL::PiPlus) {
      return std::string("pi+");
    } else if(p == G4INCL::PiZero) {
      return std::string("pi0");
    } else if(p == G4INCL::PiMinus) {
      return std::string("pi-");
    } else if(p == G4INCL::Composite) {
      return std::string("comp");
    }
    return std::string("unknown");
  }

Referenced by getShortName().

◆ getSurfaceDiffuseness()

G4double G4INCL::ParticleTable::getSurfaceDiffuseness	(	const G4int	A,
		const G4int	Z
	)

static

Definition at line 604 of file G4INCLParticleTable.cc.

                                                                             {
#else
  G4double ParticleTable::getSurfaceDiffuseness(const G4int A, const G4int Z) {
#endif // INCLXX_IN_GEANT4_MODE
 
    if(A >= 28) {
      return 1.63e-4 * A + 0.510;
    } else if(A < 28 && A >= 19) {
      return mediumDiffuseness[A-1];
    } else if(A < 19 && A >= 6) {
      return mediumDiffuseness[A-1];
    } else if(A < 6 && A >= 2) {
      ERROR("ParticleTable::getSurfaceDiffuseness: was called for A = " << A << " Z = " << Z << std::endl);
      return 0.0;
    } else {
      ERROR("ParticleTable::getSurfaceDiffuseness: No diffuseness for nucleus A = " << A << " Z = " << Z << std::endl);
      return 0.0;
    }
  }

Referenced by G4INCL::NuclearDensityFactory::createRCDFTable(), G4INCL::NuclearDensityFactory::createRPCorrelationTable(), getMaximumNuclearRadius(), getNuclearRadius(), and getSurfaceDiffuseness().

◆ getTableQValue() [1/2]

static G4double G4INCL::ParticleTable::getTableQValue	(	const G4int	A1,
		const G4int	Z1,
		const G4int	A2,
		const G4int	Z2
	)

inlinestatic

Get Q-value (in MeV/c^2)

Uses the getTableMass function to compute the Q-value for the following reaction:

$(A_1,Z_1) + (A_2, Z_2) --> (A_1+A_2,Z_1+Z_2)$

Definition at line 136 of file G4INCLParticleTable.hh.

                                                                                                   {
      return getTableMass(A1,Z1) + getTableMass(A2,Z2) - getTableMass(A1+A2,Z1+Z2);
    }

Referenced by G4INCL::Particle::getEmissionQValueCorrection(), and G4INCL::Particle::getTransferQValueCorrection().

◆ getTableQValue() [2/2]

static G4double G4INCL::ParticleTable::getTableQValue	(	const G4int	A1,
		const G4int	Z1,
		const G4int	A2,
		const G4int	Z2,
		const G4int	A3,
		const G4int	Z3
	)

inlinestatic

Get Q-value (in MeV/c^2)

Uses the getTableMass function to compute the Q-value for the following reaction:

$(A_1,Z_1) + (A_2, Z_2) --> (A_3,Z_3) + (A1+A2-A3,Z1+Z2-Z3)$

Definition at line 146 of file G4INCLParticleTable.hh.

                                                                                                                                   {
      return getTableMass(A1,Z1) + getTableMass(A2,Z2) - getTableMass(A3,Z3) - getTableMass(A1+A2-A3,Z1+Z2-Z3);
    }

◆ getTableSpeciesMass()

static G4double G4INCL::ParticleTable::getTableSpeciesMass ( const ParticleSpecies & p )

inlinestatic

Definition at line 156 of file G4INCLParticleTable.hh.

                                                                  {
      if(p.theType == Composite)
        return (*getTableMass)(p.theA, p.theZ);
      else
        return (*getTableParticleMass)(p.theType);
    }

◆ initialize()

void G4INCL::ParticleTable::initialize ( Config const *const theConfig = 0 )

static

Initialize the particle table.

Definition at line 303 of file G4INCLParticleTable.cc.

                                                                 {
    protonMass = theINCLNucleonMass;
    neutronMass = theINCLNucleonMass;
    piPlusMass = theINCLPionMass;
    piMinusMass = theINCLPionMass;
    piZeroMass = theINCLPionMass;
 
#ifndef INCLXX_IN_GEANT4_MODE
    std::string dataFilePath;
    if(theConfig)
      dataFilePath = theConfig->getINCLXXDataFilePath();
    readRealMasses(dataFilePath);
#endif
 
    if(theConfig && theConfig->getUseRealMasses()) {
      getTableMass = ParticleTable::getRealMass;
      getTableParticleMass = ParticleTable::getRealMass;
    } else {
      getTableMass = ParticleTable::getINCLMass;
      getTableParticleMass = ParticleTable::getINCLMass;
    }
#ifdef INCLXX_IN_GEANT4_MODE
    G4ParticleTable *theG4ParticleTable = G4ParticleTable::GetParticleTable();
    theG4IonTable = theG4ParticleTable->GetIonTable();
    theRealProtonMass = theG4ParticleTable->FindParticle("proton")->GetPDGMass() / MeV;
    theRealNeutronMass = theG4ParticleTable->FindParticle("neutron")->GetPDGMass() / MeV;
    theRealChargedPiMass = theG4ParticleTable->FindParticle("pi+")->GetPDGMass() / MeV;
    theRealPiZeroMass = theG4ParticleTable->FindParticle("pi0")->GetPDGMass() / MeV;
#endif
 
    // Initialise the separation-energy function
    if(!theConfig || theConfig->getSeparationEnergyType()==INCLSeparationEnergy)
      getSeparationEnergy = ParticleTable::getSeparationEnergyINCL;
    else if(theConfig->getSeparationEnergyType()==RealSeparationEnergy)
      getSeparationEnergy = ParticleTable::getSeparationEnergyReal;
    else if(theConfig->getSeparationEnergyType()==RealForLightSeparationEnergy)
      getSeparationEnergy = ParticleTable::getSeparationEnergyRealForLight;
    else {
      FATAL("Unrecognized separation-energy type in ParticleTable initialization: " << theConfig->getSeparationEnergyType() << std::endl);
      std::abort();
    }
 
  }

Referenced by G4INCL::INCL::INCL().

◆ parseIUPACElement()

G4int G4INCL::ParticleTable::parseIUPACElement ( std::string const & pS )

static

Parse a IUPAC element name.

Note: this function is UGLY. Look at it at your own peril.

Parameters

pS	a normalised string (lowercase)

Returns: the charge number of the nuclide, or zero on fail

Definition at line 687 of file G4INCLParticleTable.cc.

                                                           {
    // Normalise to lower case
    std::string elementName(s);
    std::transform(elementName.begin(), elementName.end(), elementName.begin(), ::tolower);
    // Return 0 if the element name contains anything but IUPAC digits
    if(elementName.find_first_not_of(elementIUPACDigits)!=std::string::npos)
      return 0;
    std::transform(elementName.begin(), elementName.end(), elementName.begin(), ParticleTable::iupacToInt);
    std::stringstream elementStream(elementName);
    G4int Z;
    elementStream >> Z;
    return Z;
  }

◆ setNeutronSeparationEnergy()

static void G4INCL::ParticleTable::setNeutronSeparationEnergy ( const G4double s )

inlinestatic

Setter for protonSeparationEnergy.

Definition at line 277 of file G4INCLParticleTable.hh.

277{ neutronSeparationEnergy = s; }

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

◆ setProtonSeparationEnergy()

static void G4INCL::ParticleTable::setProtonSeparationEnergy ( const G4double s )

inlinestatic

Setter for protonSeparationEnergy.

Definition at line 274 of file G4INCLParticleTable.hh.

274{ protonSeparationEnergy = s; }

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

Member Data Documentation

◆ clusterDecayMode

const ParticleTable::ClusterDecayType G4INCL::ParticleTable::clusterDecayMode

static

Initial value:

=
  {
 
    {StableCluster, StableCluster,   NeutronDecay, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,   NeutronDecay, TwoNeutronDecay,   NeutronDecay, TwoNeutronDecay, NeutronUnbound, NeutronUnbound,  NeutronUnbound, NeutronUnbound, NeutronUnbound},
    {StableCluster, StableCluster,    ProtonDecay,  StableCluster,  StableCluster,    NeutronDecay,  StableCluster,    NeutronDecay,  StableCluster,   NeutronDecay, TwoNeutronDecay, NeutronUnbound, NeutronUnbound},
    {StableCluster, StableCluster,  StableCluster,  ProtonUnbound,    ProtonDecay,     ProtonDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster,    NeutronDecay,  StableCluster,   NeutronDecay},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,  ProtonUnbound,     ProtonDecay, TwoProtonDecay,   StableCluster,     AlphaDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   ProtonUnbound, TwoProtonDecay,     ProtonDecay,  StableCluster,    ProtonDecay,   StableCluster,  StableCluster,  StableCluster},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  ProtonUnbound,   ProtonUnbound, TwoProtonDecay,  StableCluster,   StableCluster,  StableCluster,  StableCluster},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  StableCluster,   ProtonUnbound,  ProtonUnbound,  ProtonUnbound,     ProtonDecay,    ProtonDecay,  StableCluster},
    {StableCluster, StableCluster,  StableCluster,  StableCluster,  StableCluster,   StableCluster,  StableCluster,   StableCluster,  ProtonUnbound,  ProtonUnbound,   ProtonUnbound,  ProtonUnbound,    ProtonDecay}
  }

Definition at line 330 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusterDecay::isStable().

◆ clusterPhaseSpaceCut

const G4double G4INCL::ParticleTable::clusterPhaseSpaceCut

static

Initial value:

= {0.0, 70000.0, 180000.0,
                                90000.0, 90000.0,
                                128941.0 ,145607.0,
                                161365.0, 176389.0,
                                190798.0, 204681.0,
                                218109.0, 231135.0}

Definition at line 310 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::getCluster().

◆ clusterPosFact

const G4double G4INCL::ParticleTable::clusterPosFact

static

Initial value:

= {0.0, 1.0, 0.5,
                              0.33333, 0.25,
                              0.2, 0.16667,
                              0.14286, 0.125,
                              0.11111, 0.1,
                              0.09091, 0.083333}

Precomputed factor 1.0/A

Definition at line 306 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusterUtils::getNewPositionVector().

◆ clusterPosFact2

const G4double G4INCL::ParticleTable::clusterPosFact2

static

Initial value:

= {0.0, 1.0, 0.25,
                               0.11111, 0.0625,
                               0.04, 0.0277778,
                               0.020408, 0.015625,
                               0.012346, 0.01,
                               0.0082645, 0.0069444}

Precomputed factor (1.0/A)^2

Definition at line 307 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::getCluster(), and G4INCL::ClusterUtils::getPhaseSpace().

◆ clusterTableASize

const G4int G4INCL::ParticleTable::clusterTableASize = ParticleTable::maxClusterMass+1

static

Definition at line 305 of file G4INCLParticleTable.hh.

Referenced by getMomentumRMS(), getNuclearRadius(), and getSeparationEnergyRealForLight().

◆ clusterTableZSize

const G4int G4INCL::ParticleTable::clusterTableZSize = ParticleTable::maxClusterCharge+1

static

Definition at line 304 of file G4INCLParticleTable.hh.

Referenced by getMomentumRMS(), getRadiusParameter(), and getSeparationEnergyRealForLight().

◆ clusterZMax

const G4int G4INCL::ParticleTable::clusterZMax = {0, 0, 1, 2, 3, 3, 5, 5, 6, 6, 7, 7, 8}

static

Definition at line 309 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::ClusteringModelIntercomparison().

◆ clusterZMin

const G4int G4INCL::ParticleTable::clusterZMin = {0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3}

static

Definition at line 308 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::ClusteringModelIntercomparison().

◆ effectiveDeltaDecayThreshold

const G4double G4INCL::ParticleTable::effectiveDeltaDecayThreshold

static

Initial value:

=
    ParticleTable::theRealNeutronMass + ParticleTable::theRealChargedPiMass
    + 0.5

Definition at line 299 of file G4INCLParticleTable.hh.

Referenced by G4INCL::InteractionAvatar::enforceEnergyConservation(), and G4INCL::InteractionAvatar::postInteraction().

◆ effectiveDeltaMass

const G4double G4INCL::ParticleTable::effectiveDeltaMass = 1232.0

static

Definition at line 297 of file G4INCLParticleTable.hh.

◆ effectiveNucleonMass

const G4double G4INCL::ParticleTable::effectiveNucleonMass = 938.2796

static

Definition at line 295 of file G4INCLParticleTable.hh.

Referenced by G4INCL::DeltaDecayChannel::computeDecayTime(), G4INCL::Nucleus::computeTotalEnergy(), G4INCL::CrossSections::deltaProduction(), G4INCL::DeltaProductionChannel::getFinalState(), G4INCL::ElasticChannel::getFinalState(), and G4INCL::CrossSections::recombination().

◆ effectiveNucleonMass2

const G4double G4INCL::ParticleTable::effectiveNucleonMass2 = 8.8036860777616e5

static

Definition at line 296 of file G4INCLParticleTable.hh.

Referenced by G4INCL::DeltaProductionChannel::getFinalState(), and G4INCL::CrossSections::recombination().

◆ effectivePionMass

const G4double G4INCL::ParticleTable::effectivePionMass = 138.0

static

Definition at line 298 of file G4INCLParticleTable.hh.

Referenced by G4INCL::DeltaDecayChannel::computeDecayTime(), and G4INCL::CrossSections::deltaProduction().

◆ elementTableSize

const G4int G4INCL::ParticleTable::elementTableSize = 113

static

Definition at line 293 of file G4INCLParticleTable.hh.

Referenced by getElementName().

◆ eSquared

const G4double G4INCL::ParticleTable::eSquared = 1.439964

static

Coulomb conversion factor, in MeV*fm.

$e^2/(4 pi epsilon_0)$

Definition at line 336 of file G4INCLParticleTable.hh.

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

◆ getSeparationEnergy

ParticleTable::SeparationEnergyFn G4INCL::ParticleTable::getSeparationEnergy

static

Static pointer to the separation-energy function.

Definition at line 165 of file G4INCLParticleTable.hh.

Referenced by initialize().

◆ getTableMass

ParticleTable::NuclearMassFn G4INCL::ParticleTable::getTableMass

static

Static pointer to the mass function for nuclei.

Definition at line 153 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ProjectileRemnant::addMostDynamicalSpectators(), G4INCL::Nucleus::computeRecoilKinematics(), G4INCL::Nucleus::finalizeProjectileRemnant(), G4INCL::KinematicsUtils::gammaFromKineticEnergy(), G4INCL::Nucleus::getConservationBalance(), G4INCL::Particle::getEmissionQValueCorrection(), G4INCL::ParticleEntryChannel::getFinalState(), getSeparationEnergyReal(), G4INCL::Particle::getTableMass(), getTableQValue(), getTableSpeciesMass(), initialize(), G4INCL::StandardPropagationModel::shootComposite(), and G4INCL::StandardPropagationModel::shootParticle().

◆ getTableParticleMass

ParticleTable::ParticleMassFn G4INCL::ParticleTable::getTableParticleMass

static

Static pointer to the mass function for particles.

Definition at line 154 of file G4INCLParticleTable.hh.

Referenced by G4INCL::KinematicsUtils::gammaFromKineticEnergy(), getSeparationEnergyReal(), G4INCL::Particle::getTableMass(), getTableSpeciesMass(), initialize(), and G4INCL::StandardPropagationModel::shootParticle().

◆ maxClusterCharge

const G4int G4INCL::ParticleTable::maxClusterCharge = 8

static

Definition at line 302 of file G4INCLParticleTable.hh.

◆ maxClusterMass

const G4int G4INCL::ParticleTable::maxClusterMass = 12

static

Definition at line 301 of file G4INCLParticleTable.hh.

Referenced by G4INCL::ClusteringModelIntercomparison::ClusteringModelIntercomparison().

◆ theG4IonTable

G4IonTable * G4INCL::ParticleTable::theG4IonTable

static

Definition at line 313 of file G4INCLParticleTable.hh.

Referenced by getRealMass(), and initialize().

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

Public Types

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Detailed Description

Member Typedef Documentation

◆ NuclearMassFn

◆ ParticleMassFn

◆ SeparationEnergyFn

Member Enumeration Documentation

◆ ClusterDecayType

Constructor & Destructor Documentation

◆ ParticleTable()

◆ ~ParticleTable()

Member Function Documentation

◆ drawRandomNaturalIsotope()

◆ getChargeNumber()

◆ getElementName()

◆ getINCLMass() [1/2]

◆ getINCLMass() [2/2]

◆ getIsospin()

◆ getIUPACElementName()

◆ getMassNumber()

◆ getMaximumNuclearRadius()

◆ getMomentumRMS()

◆ getName() [1/3]

◆ getName() [2/3]

◆ getName() [3/3]

◆ getNaturalIsotopicDistribution()

◆ getNeutronSeparationEnergy()

◆ getNuclearRadius()

◆ getProtonSeparationEnergy()

◆ getRadiusParameter()

◆ getRealMass() [1/2]

◆ getRealMass() [2/2]

◆ getSeparationEnergyINCL()

◆ getSeparationEnergyReal()

◆ getSeparationEnergyRealForLight()

◆ getShortName() [1/3]

◆ getShortName() [2/3]

◆ getShortName() [3/3]

◆ getSurfaceDiffuseness()

◆ getTableQValue() [1/2]

◆ getTableQValue() [2/2]

◆ getTableSpeciesMass()

◆ initialize()

◆ parseIUPACElement()

◆ setNeutronSeparationEnergy()

◆ setProtonSeparationEnergy()

Member Data Documentation

◆ clusterDecayMode

◆ clusterPhaseSpaceCut

◆ clusterPosFact

◆ clusterPosFact2

◆ clusterTableASize

◆ clusterTableZSize

◆ clusterZMax

◆ clusterZMin

◆ effectiveDeltaDecayThreshold

◆ effectiveDeltaMass

◆ effectiveNucleonMass

◆ effectiveNucleonMass2

◆ effectivePionMass

◆ elementTableSize

◆ eSquared

◆ getSeparationEnergy

◆ getTableMass

◆ getTableParticleMass

◆ maxClusterCharge

◆ maxClusterMass

◆ theG4IonTable