G4NucleiProperties Class Reference

#include <G4NucleiProperties.hh>

Public Member Functions
	G4NucleiProperties ()
	~G4NucleiProperties ()

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 ( )

inline

Definition at line 48 of file G4NucleiProperties.hh.

{}

~G4NucleiProperties()

G4NucleiProperties::~G4NucleiProperties ( )

inline

Definition at line 49 of file G4NucleiProperties.hh.

{}

Member Function Documentation

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

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

Referenced by G4ParticleHPInelasticBaseFS::BaseApply(), 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(), 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);
}

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

Referenced by 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;
}

Referenced by G4EMDissociation::ApplyYourself(), G4EmCaptureCascade::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4LowEGammaNuclearModel::ApplyYourself(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4ParaFissionModel::ApplyYourself(), G4INCLXXInterface::ApplyYourself(), G4LMsdGenerator::ApplyYourself(), G4LowEIonFragmentation::ApplyYourself(), G4PreCompoundModel::ApplyYourself(), G4ParticleHPCaptureFS::ApplyYourself(), 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(), G4ParticleHPCaptureData::GetCrossSection(), G4ParticleHPElasticData::GetCrossSection(), G4ParticleHPFissionData::GetCrossSection(), G4ParticleHPInelasticData::GetCrossSection(), G4IonsKoxCrossSection::GetElementCrossSection(), G4TripathiLightCrossSection::GetElementCrossSection(), G4EvaporationChannel::GetEmissionProbability(), G4GEMChannel::GetEmissionProbability(), G4GEMChannelVI::GetEmissionProbability(), G4IonsShenCrossSection::GetIsoCrossSection(), G4MuonMinusBoundDecay::GetMuonDecayRate(), G4QMDNucleus::GetNuclearMass(), GetNuclearMass(), G4HyperNucleiProperties::GetNuclearMass(), G4InuclNuclei::getNucleiMass(), G4IonTable::GetNucleusMass(), G4ParticleHPThermalBoost::GetThermalEnergy(), 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(), and G4ScreeningMottCrossSection::SetupKinematic().

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

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

Referenced by IsInStableTable().

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

---

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

• G4NucleiProperties.hh
• G4NucleiProperties.cc