G4NucleiProperties Class Reference

#include <G4NucleiProperties.hh>

Public Member Functions
	G4NucleiProperties ()
	~G4NucleiProperties ()

Static Public Member Functions
static G4double	GetBindingEnergy (const G4double A, const G4double Z)
static G4double	GetBindingEnergy (const G4int A, const G4int Z)
static G4double	GetMassExcess (const G4double A, const G4double Z)
static G4double	GetMassExcess (const G4int A, const G4int Z)
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)

Private Types
enum	{ MaxZ = 120 }

Static Private Member Functions
static G4double	AtomicMass (G4double A, G4double Z)
static G4double	BindingEnergy (G4double A, G4double Z)
static G4double	GetAtomicMass (const G4double A, const G4double Z)
static G4double	MassExcess (G4double A, G4double Z)
static G4double	NuclearMass (G4double A, G4double Z)

Static Private Attributes
static G4ThreadLocal G4double	electronMass [MaxZ]
static G4ThreadLocal G4bool	isIntialized
static G4ThreadLocal G4double	mass_alpha = -1.
static G4ThreadLocal G4double	mass_deuteron = -1.
static G4ThreadLocal G4double	mass_He3 = -1.
static G4ThreadLocal G4double	mass_neutron = -1.
static G4ThreadLocal G4double	mass_proton = -1.
static G4ThreadLocal G4double	mass_triton = -1.

Detailed Description

Definition at line 44 of file G4NucleiProperties.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
MaxZ

Definition at line 85 of file G4NucleiProperties.hh.

{MaxZ = 120};

Constructor & Destructor Documentation

G4NucleiProperties()

G4NucleiProperties::G4NucleiProperties ( )

inline

Definition at line 48 of file G4NucleiProperties.hh.

{}

~G4NucleiProperties()

G4NucleiProperties::~G4NucleiProperties ( )

inline

Definition at line 49 of file G4NucleiProperties.hh.

{}

Member Function Documentation

AtomicMass()

G4double G4NucleiProperties::AtomicMass ( G4double  A, G4double  Z  )

staticprivate

Definition at line 277 of file G4NucleiProperties.cc.

{
  G4double hydrogen_mass_excess;
  G4double neutron_mass_excess;  
  hydrogen_mass_excess = G4NucleiPropertiesTableAME12::GetMassExcess(1,1);
  neutron_mass_excess = G4NucleiPropertiesTableAME12::GetMassExcess(0,1);
  G4double mass = (A-Z)*neutron_mass_excess
                + Z*hydrogen_mass_excess - BindingEnergy(A,Z) + A*amu_c2;
  return mass;
}

References A, source.hepunit::amu_c2, BindingEnergy(), G4NucleiPropertiesTableAME12::GetMassExcess(), and Z.

Referenced by GetAtomicMass(), and NuclearMass().

BindingEnergy()

G4double G4NucleiProperties::BindingEnergy ( G4double  A, G4double  Z  )

staticprivate

Definition at line 314 of file G4NucleiProperties.cc.

{ 
  //
  // Weitzsaecker's Mass formula
  //
  G4int Npairing = G4int(A-Z)%2;          // pairing
  G4int Zpairing = G4int(Z)%2;
  G4double binding =
      - 15.67*A                           // nuclear volume
      + 17.23*std::pow(A,2./3.)           // surface energy
      + 93.15*((A/2.-Z)*(A/2.-Z))/A       // asymmetry
      + 0.6984523*Z*Z*std::pow(A,-1./3.); // coulomb
  if( Npairing == Zpairing )
  {
    binding += (Npairing+Zpairing-1) * 12.0 / std::sqrt(A);  // pairing
  }
 
  return -binding*MeV;
}

References A, MeV, and Z.

Referenced by AtomicMass(), and GetBindingEnergy().

GetAtomicMass()

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

staticprivate

Definition at line 206 of file G4NucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A)
  {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>0)
    {
      G4cout << "G4NucleiProperties::GetAtomicMass: Wrong values for A = " 
             << A << " and Z = " << Z << G4endl;
    }
#endif 
    return 0.0;
 
  }
  else if (std::fabs(A - G4int(A)) > 1.e-10)
  {
    return AtomicMass(A,Z);
  }
  else
  {
    G4int iA = G4int(A);
    G4int iZ = G4int(Z);
    if ( G4NucleiPropertiesTableAME12::IsInTable(Z,A) ) {
      return G4NucleiPropertiesTableAME12::GetAtomicMass(Z,A);
    } else if (G4NucleiPropertiesTheoreticalTable::IsInTable(iZ,iA)){
      return G4NucleiPropertiesTheoreticalTable::GetAtomicMass(iZ,iA);
    } else {
      return AtomicMass(A,Z);
    }
  }
}

References A, AtomicMass(), G4cout, G4endl, G4NucleiPropertiesTableAME12::GetAtomicMass(), G4NucleiPropertiesTheoreticalTable::GetAtomicMass(), G4ParticleTable::GetParticleTable(), G4NucleiPropertiesTableAME12::IsInTable(), G4NucleiPropertiesTheoreticalTable::IsInTable(), and Z.

Referenced by MassExcess().

GetBindingEnergy() [1/2]

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

static

Definition at line 239 of file G4NucleiProperties.cc.

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

References A, GetBindingEnergy(), and Z.

GetBindingEnergy() [2/2]

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

static

Definition at line 246 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;
 
  }
  else
  {
    if ( G4NucleiPropertiesTableAME12::IsInTable(Z,A) ) {
      return G4NucleiPropertiesTableAME12::GetBindingEnergy(Z,A);
    } else if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z,A)) {
      return G4NucleiPropertiesTheoreticalTable::GetBindingEnergy(Z,A);
    }else {
      return BindingEnergy(A,Z);
    }
  }
}

References A, BindingEnergy(), G4cout, G4endl, G4NucleiPropertiesTableAME12::GetBindingEnergy(), G4NucleiPropertiesTheoreticalTable::GetBindingEnergy(), G4ParticleTable::GetParticleTable(), G4NucleiPropertiesTableAME12::IsInTable(), G4NucleiPropertiesTheoreticalTable::IsInTable(), and Z.

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), G4Fancy3DNucleus::BindingEnergy(), G4InuclSpecialFunctions::bindingEnergy(), G4StatMFMacroBiNucleon::CalcEnergy(), G4StatMFMacroTetraNucleon::CalcEnergy(), G4StatMFMacroTriNucleon::CalcEnergy(), G4StatMFMacroBiNucleon::CalcMeanMultiplicity(), G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(), G4StatMFMacroTriNucleon::CalcMeanMultiplicity(), 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(), G4QMDGroundStateNucleus::packNucleons(), and G4KM_OpticalEqRhs::SetFactor().

GetMassExcess() [1/2]

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

static

Definition at line 170 of file G4NucleiProperties.cc.

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

References A, GetMassExcess(), and Z.

GetMassExcess() [2/2]

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

static

Definition at line 177 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;
    
  }
  else
  {
 
    if ( G4NucleiPropertiesTableAME12::IsInTable(Z,A) ){
      // AME table
      return G4NucleiPropertiesTableAME12::GetMassExcess(Z,A);
    } else if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z,A)){
      return G4NucleiPropertiesTheoreticalTable::GetMassExcess(Z,A);
    } else {
      return MassExcess(A,Z);
    }
  }
}

References A, G4cout, G4endl, G4NucleiPropertiesTableAME12::GetMassExcess(), G4NucleiPropertiesTheoreticalTable::GetMassExcess(), G4ParticleTable::GetParticleTable(), G4NucleiPropertiesTableAME12::IsInTable(), G4NucleiPropertiesTheoreticalTable::IsInTable(), MassExcess(), and Z.

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

GetNuclearMass() [1/2]

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

static

Definition at line 50 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
    G4int iZ = G4int(Z);
    G4int iA = G4int(A);
    mass =GetNuclearMass(iA,iZ);
  }
  return mass;
}

References A, GetNuclearMass(), NuclearMass(), and Z.

Referenced by G4EMDissociation::ApplyYourself(), G4EmCaptureCascade::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4INCLXXInterface::ApplyYourself(), G4LMsdGenerator::ApplyYourself(), G4PreCompoundModel::ApplyYourself(), G4LowEIonFragmentation::ApplyYourself(), G4ParticleHPCaptureFS::ApplyYourself(), G4Nucleus::AtomicMass(), G4ParticleHPInelasticBaseFS::BaseApply(), G4InuclEvaporation::BreakItUp(), G4UnstableFragmentBreakUp::BreakUpChain(), G4Fragment::CalculateGroundStateMass(), 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(), G4ParticleHPCaptureData::GetCrossSection(), G4ParticleHPElasticData::GetCrossSection(), G4ParticleHPFissionData::GetCrossSection(), G4ParticleHPInelasticData::GetCrossSection(), G4EvaporationChannel::GetEmissionProbability(), G4GEMChannel::GetEmissionProbability(), G4GEMChannelVI::GetEmissionProbability(), G4IonsShenCrossSection::GetIsoCrossSection(), G4MuonMinusBoundDecay::GetMuonDecayRate(), G4QMDNucleus::GetNuclearMass(), GetNuclearMass(), G4HyperNucleiProperties::GetNuclearMass(), G4InuclNuclei::getNucleiMass(), G4IonTable::GetNucleusMass(), G4ParticleHPThermalBoost::GetThermalEnergy(), G4FermiFragmentsPoolVI::Initialise(), G4VPreCompoundFragment::Initialize(), G4eCoulombScatteringModel::MinPrimaryEnergy(), G4hCoulombScatteringModel::MinPrimaryEnergy(), G4GeneratorPrecompoundInterface::Propagate(), G4GeneratorPrecompoundInterface::PropagateNuclNucl(), G4ParticleHPContAngularPar::Sample(), G4DiffuseElastic::SampleInvariantT(), G4DiffuseElasticV2::SampleInvariantT(), G4NuclNuclDiffuseElastic::SampleInvariantT(), G4AntiNuclElastic::SampleInvariantT(), G4eSingleCoulombScatteringModel::SampleSecondaries(), G4IonCoulombScatteringModel::SampleSecondaries(), G4eCoulombScatteringModel::SampleSecondaries(), G4hCoulombScatteringModel::SampleSecondaries(), G4BetheHeitler5DModel::SampleSecondaries(), G4ScreeningMottCrossSection::SetupKinematic(), G4ElectroNuclearCrossSection::ThresholdEnergy(), and G4PhotoNuclearCrossSection::ThresholdEnergy().

GetNuclearMass() [2/2]

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

static

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

References A, G4ParticleTable::FindParticle(), G4cout, G4endl, G4NucleiPropertiesTableAME12::GetNuclearMass(), G4NucleiPropertiesTheoreticalTable::GetNuclearMass(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4NucleiPropertiesTableAME12::IsInTable(), G4NucleiPropertiesTheoreticalTable::IsInTable(), mass_alpha, mass_deuteron, mass_He3, mass_neutron, mass_proton, mass_triton, NuclearMass(), and Z.

IsInStableTable() [1/2]

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

static

Definition at line 146 of file G4NucleiProperties.cc.

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

References A, IsInStableTable(), and Z.

Referenced by IsInStableTable(), G4ElectroNuclearCrossSection::ThresholdEnergy(), and G4PhotoNuclearCrossSection::ThresholdEnergy().

IsInStableTable() [2/2]

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

static

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

References A, G4cout, G4endl, G4ParticleTable::GetParticleTable(), G4NucleiPropertiesTableAME12::IsInTable(), and Z.

MassExcess()

G4double G4NucleiProperties::MassExcess ( G4double  A, G4double  Z  )

staticprivate

Definition at line 272 of file G4NucleiProperties.cc.

{
  return GetAtomicMass(A,Z) - A*amu_c2;
}

References A, source.hepunit::amu_c2, GetAtomicMass(), and Z.

Referenced by GetMassExcess().

NuclearMass()

G4double G4NucleiProperties::NuclearMass ( G4double  A, G4double  Z  )

staticprivate

Definition at line 288 of file G4NucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A)
  {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>0)
    {
      G4cout << "G4NucleiProperties::NuclearMass: Wrong values for A = " 
             << A << " and Z = " << Z << G4endl;
    }
#endif 
    return 0.0;
  }
 
  G4double mass = AtomicMass(A,Z);
 
  // atomic mass is converted to nuclear mass according to
  // formula in  AME03 and 12
  //
  mass -= Z*electron_mass_c2;
  mass += ( 14.4381*std::pow ( Z , 2.39 )
        + 1.55468*1e-6*std::pow ( Z , 5.35 ) )*eV;      
 
  return mass;
}

References A, AtomicMass(), source.hepunit::electron_mass_c2, eV, G4cout, G4endl, G4ParticleTable::GetParticleTable(), and Z.

Referenced by GetNuclearMass().

Field Documentation

electronMass

G4ThreadLocal G4double G4NucleiProperties::electronMass[MaxZ]

staticprivate

Definition at line 86 of file G4NucleiProperties.hh.

isIntialized

G4ThreadLocal G4bool G4NucleiProperties::isIntialized

staticprivate

Definition at line 89 of file G4NucleiProperties.hh.

mass_alpha

G4ThreadLocal G4double G4NucleiProperties::mass_alpha = -1.

staticprivate

Definition at line 94 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

mass_deuteron

G4ThreadLocal G4double G4NucleiProperties::mass_deuteron = -1.

staticprivate

Definition at line 92 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

mass_He3

G4ThreadLocal G4double G4NucleiProperties::mass_He3 = -1.

staticprivate

Definition at line 95 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

mass_neutron

G4ThreadLocal G4double G4NucleiProperties::mass_neutron = -1.

staticprivate

Definition at line 91 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

mass_proton

G4ThreadLocal G4double G4NucleiProperties::mass_proton = -1.

staticprivate

Definition at line 90 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

mass_triton

G4ThreadLocal G4double G4NucleiProperties::mass_triton = -1.

staticprivate

Definition at line 93 of file G4NucleiProperties.hh.

Referenced by GetNuclearMass().

---

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

• source/particles/management/include/G4NucleiProperties.hh
• source/particles/management/src/G4NucleiProperties.cc