G4NucleiProperties Class Reference

#include <G4NucleiProperties.hh>

Public Member Functions
	G4NucleiProperties ()=default
	~G4NucleiProperties ()=default

Static Public Member Functions
static G4double	GetNuclearMass (const G4double A, const G4double Z)
static G4double	GetNuclearMass (const G4int A, const G4int Z)
static G4bool	IsInStableTable (const G4double A, const G4double Z)
static G4bool	IsInStableTable (const G4int A, const G4int Z)
static G4double	GetBindingEnergy (const G4int A, const G4int Z)
static G4double	GetBindingEnergy (const G4double A, const G4double Z)
static G4double	GetMassExcess (const G4int A, const G4int Z)
static G4double	GetMassExcess (const G4double A, const G4double Z)

Detailed Description

Definition at line 44 of file G4NucleiProperties.hh.

Constructor & Destructor Documentation

G4NucleiProperties()

G4NucleiProperties::G4NucleiProperties ( )

default

~G4NucleiProperties()

G4NucleiProperties::~G4NucleiProperties ( )

default

Member Function Documentation

GetBindingEnergy() [1/2]

G4double G4NucleiProperties::GetBindingEnergy ( const G4double A, const G4double Z )

static

Definition at line 222 of file G4NucleiProperties.cc.

{
  auto iA = G4int(A);
  auto iZ = G4int(Z);
  return GetBindingEnergy(iA, iZ);
}

GetBindingEnergy() [2/2]

G4double G4NucleiProperties::GetBindingEnergy ( const G4int A, const G4int Z )

static

Definition at line 229 of file G4NucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4NucleiProperties::GetMassExccess: Wrong values for A = " << A << " and Z = " << Z
             << G4endl;
    }
#endif
    return 0.0;
  }
 
  if (G4NucleiPropertiesTableAME12::IsInTable(Z, A)) {
    return G4NucleiPropertiesTableAME12::GetBindingEnergy(Z, A);
  }
  if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z, A)) {
    return G4NucleiPropertiesTheoreticalTable::GetBindingEnergy(Z, A);
  }
  return BindingEnergy(A, Z);
}

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), G4InuclSpecialFunctions::bindingEnergy(), G4StatMFMacroBiNucleon::CalcEnergy(), G4StatMFMacroTetraNucleon::CalcEnergy(), G4StatMFMacroTriNucleon::CalcEnergy(), G4StatMFMacroBiNucleon::CalcMeanMultiplicity(), G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(), G4StatMFMacroTriNucleon::CalcMeanMultiplicity(), G4LightIonQMDNucleus::CalEnergyAndAngularMomentumInCM(), G4QMDNucleus::CalEnergyAndAngularMomentumInCM(), GetBindingEnergy(), G4AntiProtonField::GetField(), G4KaonMinusField::GetField(), G4KaonPlusField::GetField(), G4KaonZeroField::GetField(), G4PionMinusField::GetField(), G4PionPlusField::GetField(), G4PionZeroField::GetField(), G4SigmaMinusField::GetField(), G4SigmaPlusField::GetField(), G4SigmaZeroField::GetField(), G4ParticleHPInelasticBaseFS::InitGammas(), and G4KM_OpticalEqRhs::SetFactor().

GetMassExcess() [1/2]

G4double G4NucleiProperties::GetMassExcess ( const G4double A, const G4double Z )

static

Definition at line 167 of file G4NucleiProperties.cc.

{
  auto iA = G4int(A);
  auto iZ = G4int(Z);
  return GetMassExcess(iA, iZ);
}

GetMassExcess() [2/2]

G4double G4NucleiProperties::GetMassExcess ( const G4int A, const G4int Z )

static

Definition at line 174 of file G4NucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4NucleiProperties::GetMassExccess: Wrong values for A = " << A << " and Z = " << Z
             << G4endl;
    }
#endif
    return 0.0;
  }
 
  if (G4NucleiPropertiesTableAME12::IsInTable(Z, A)) {
    // AME table
    return G4NucleiPropertiesTableAME12::GetMassExcess(Z, A);
  }
  if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z, A)) {
    return G4NucleiPropertiesTheoreticalTable::GetMassExcess(Z, A);
  }
  return MassExcess(A, Z);
}

Referenced by G4StatMFFragment::GetEnergy(), and GetMassExcess().

GetNuclearMass() [1/2]

G4double G4NucleiProperties::GetNuclearMass ( const G4double A, const G4double Z )

static

Definition at line 49 of file G4NucleiProperties.cc.

{
  G4double mass = 0.0;
 
  if (std::fabs(A - G4int(A)) > 1.e-10) {
    mass = NuclearMass(A, Z);
  }
  else {
    // use mass table
    auto iZ = G4int(Z);
    auto iA = G4int(A);
    mass = GetNuclearMass(iA, iZ);
  }
  return mass;
}

Referenced by G4AblaInterface::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4EmCaptureCascade::ApplyYourself(), G4EMDissociation::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4INCLXXInterface::ApplyYourself(), G4LMsdGenerator::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4LowEIonFragmentation::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4NeutronFissionVI::ApplyYourself(), G4NeutronHPCaptureFS::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4NeutronRadCaptureHP::ApplyYourself(), G4PreCompoundModel::ApplyYourself(), G4Nucleus::AtomicMass(), G4Nucleus::AtomicMass(), G4ParticleHPInelasticBaseFS::BaseApply(), G4InuclEvaporation::BreakItUp(), G4UnstableFragmentBreakUp::BreakUpChain(), G4HadronicProcess::CheckEnergyMomentumConservation(), G4HadronicProcess::CheckResult(), G4NeutrinoNucleusModel::CoherentPion(), G4ParticleHPInelasticCompFS::CompositeApply(), G4ComponentGGNuclNuclXsc::ComputeCoulombBarier(), G4Fragment::ComputeGroundStateMass(), G4HadronElastic::ComputeMomentumCMS(), G4NuclearRadii::CoulombFactor(), G4GEMProbability::Dump(), G4CompetitiveFission::EmittedFragment(), G4NeutrinoNucleusModel::FinalBarion(), G4ElasticData::G4ElasticData(), G4EvaporationChannel::G4EvaporationChannel(), G4FermiFragment::G4FermiFragment(), G4FermiPair::G4FermiPair(), G4GEMChannel::G4GEMChannel(), G4GEMChannelVI::G4GEMChannelVI(), G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp(), G4VEmissionProbability::G4VEmissionProbability(), G4NeutronHPCaptureData::GetCrossSection(), G4ParticleHPElasticData::GetCrossSection(), G4ParticleHPFissionData::GetCrossSection(), G4ParticleHPInelasticData::GetCrossSection(), G4EvaporationChannel::GetEmissionProbability(), G4GEMChannel::GetEmissionProbability(), G4GEMChannelVI::GetEmissionProbability(), G4ChargeExchangeXS::GetIsoCrossSection(), G4IonsShenCrossSection::GetIsoCrossSection(), G4MuonMinusBoundDecay::GetMuonDecayRate(), G4HyperNucleiProperties::GetNuclearMass(), G4LightIonQMDNucleus::GetNuclearMass(), GetNuclearMass(), G4QMDNucleus::GetNuclearMass(), G4InuclNuclei::getNucleiMass(), G4IonTable::GetNucleusMass(), G4ParticleHPThermalBoost::GetThermalEnergy(), G4FermiFragmentsPoolVI::Initialise(), G4VPreCompoundFragment::Initialize(), G4eCoulombScatteringModel::MinPrimaryEnergy(), G4hCoulombScatteringModel::MinPrimaryEnergy(), G4GeneratorPrecompoundInterface::Propagate(), G4GeneratorPrecompoundInterface::PropagateNuclNucl(), G4ParticleHPContAngularPar::Sample(), G4AntiNuclElastic::SampleInvariantT(), G4DiffuseElastic::SampleInvariantT(), G4DiffuseElasticV2::SampleInvariantT(), G4NuclNuclDiffuseElastic::SampleInvariantT(), G4BetheHeitler5DModel::SampleSecondaries(), G4eCoulombScatteringModel::SampleSecondaries(), G4eSingleCoulombScatteringModel::SampleSecondaries(), G4hCoulombScatteringModel::SampleSecondaries(), G4IonCoulombScatteringModel::SampleSecondaries(), and G4ScreeningMottCrossSection::SetupKinematic().

GetNuclearMass() [2/2]

G4double G4NucleiProperties::GetNuclearMass ( const G4int A, const G4int Z )

static

Definition at line 65 of file G4NucleiProperties.cc.

{
  if (mass_proton <= 0.0) {
    const G4ParticleDefinition* nucleus = nullptr;
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("neutron");
    if (nucleus != nullptr) mass_neutron = nucleus->GetPDGMass();
 
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("deuteron");
    if (nucleus != nullptr) mass_deuteron = nucleus->GetPDGMass();
 
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("triton");
    if (nucleus != nullptr) mass_triton = nucleus->GetPDGMass();
 
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("alpha");
    if (nucleus != nullptr) mass_alpha = nucleus->GetPDGMass();
 
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("He3");
    if (nucleus != nullptr) mass_He3 = nucleus->GetPDGMass();
 
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("proton");
    if (nucleus != nullptr) mass_proton = nucleus->GetPDGMass();
  }
 
  if (A < 1 || Z < 0 || Z > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4NucleiProperties::GetNuclearMass: Wrong values for A = " << A << " and Z = " << Z
             << G4endl;
    }
#endif
    return 0.0;
  }
 
  G4double mass = -1.;
  if ((Z <= 2)) {
    // light nuclei
    if ((Z == 1) && (A == 1)) {
      mass = mass_proton;
    }
    else if ((Z == 0) && (A == 1)) {
      mass = mass_neutron;
    }
    else if ((Z == 1) && (A == 2)) {
      mass = mass_deuteron;
    }
    else if ((Z == 1) && (A == 3)) {
      mass = mass_triton;
    }
    else if ((Z == 2) && (A == 4)) {
      mass = mass_alpha;
    }
    else if ((Z == 2) && (A == 3)) {
      mass = mass_He3;
    }
  }
 
  if (mass < 0.) {
    if (G4NucleiPropertiesTableAME12::IsInTable(Z, A)) {
      // AME table
      mass = G4NucleiPropertiesTableAME12::GetNuclearMass(Z, A);
    }
    else if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z, A)) {
      // Theoretical table
      mass = G4NucleiPropertiesTheoreticalTable::GetNuclearMass(Z, A);
    }
    else if (Z == A) {
      mass = A * mass_proton;
    }
    else if (0 == Z) {
      mass = A * mass_neutron;
    }
    else {
      mass = NuclearMass(G4double(A), G4double(Z));
    }
  }
 
  if (mass < 0.) mass = 0.0;
  return mass;
}

IsInStableTable() [1/2]

G4bool G4NucleiProperties::IsInStableTable ( const G4double A, const G4double Z )

static

Definition at line 145 of file G4NucleiProperties.cc.

{
  auto iA = G4int(A);
  auto iZ = G4int(Z);
  return IsInStableTable(iA, iZ);
}

Referenced by IsInStableTable().

IsInStableTable() [2/2]

G4bool G4NucleiProperties::IsInStableTable ( const G4int A, const G4int Z )

static

Definition at line 152 of file G4NucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4NucleiProperties::IsInStableTable: Wrong values for A = " << A
             << " and Z = " << Z << G4endl;
    }
#endif
    return false;
  }
 
  return G4NucleiPropertiesTableAME12::IsInTable(Z, A);
}

---

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

• G4NucleiProperties.hh
• G4NucleiProperties.cc