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 bool	IsInStableTable (const G4double A, const G4double Z)
static bool	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 50 of file G4NucleiProperties.hh.

Constructor & Destructor Documentation

~G4NucleiProperties()

G4NucleiProperties::~G4NucleiProperties ( )

inline

Definition at line 59 of file G4NucleiProperties.hh.

{ };

G4NucleiProperties()

G4NucleiProperties::G4NucleiProperties ( )

inline

Definition at line 62 of file G4NucleiProperties.hh.

{};

Member Function Documentation

GetBindingEnergy() [1/2]

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

static

Definition at line 216 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 223 of file G4NucleiProperties.cc.

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

Referenced by G4KaonMinusAbsorptionAtRest::AtRestDoIt(), G4NeutronHPInelasticBaseFS::BaseApply(), G4InuclSpecialFunctions::bindingEnergy(), G4StatMFMacroBiNucleon::CalcEnergy(), G4StatMFMacroTetraNucleon::CalcEnergy(), G4StatMFMacroTriNucleon::CalcEnergy(), G4StatMFMacroBiNucleon::CalcMeanMultiplicity(), G4StatMFMacroTetraNucleon::CalcMeanMultiplicity(), G4StatMFMacroTriNucleon::CalcMeanMultiplicity(), G4QMDNucleus::CalEnergyAndAngularMomentumInCM(), G4PiMinusStopAbsorption::DoAbsorption(), G4PiMinusStopAbsorption::Energy(), GetBindingEnergy(), G4AntiProtonField::GetField(), G4KaonMinusField::GetField(), G4KaonPlusField::GetField(), G4KaonZeroField::GetField(), G4PionMinusField::GetField(), G4PionPlusField::GetField(), G4PionZeroField::GetField(), G4SigmaMinusField::GetField(), G4SigmaPlusField::GetField(), G4SigmaZeroField::GetField(), G4NeutronHPInelasticBaseFS::InitGammas(), G4NeutronHPContAngularPar::Sample(), and G4KM_OpticalEqRhs::SetFactor().

GetMassExcess() [1/2]

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

static

Definition at line 157 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 164 of file G4NucleiProperties.cc.

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

Referenced by G4BEDeuteronChannel::G4BEDeuteronChannel(), G4StatMFFragment::GetEnergy(), GetMassExcess(), and G4BertiniEvaporationChannel::getQ().

GetNuclearMass() [1/2]

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

static

Definition at line 53 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(), G4ChargeExchange::ApplyYourself(), G4HadronElastic::ApplyYourself(), G4NeutronRadCapture::ApplyYourself(), G4EmCaptureCascade::ApplyYourself(), G4MuMinusCapturePrecompound::ApplyYourself(), G4ParaFissionModel::ApplyYourself(), G4INCLXXInterface::ApplyYourself(), G4LowEIonFragmentation::ApplyYourself(), G4PreCompoundModel::ApplyYourself(), G4NeutronHPCaptureFS::ApplyYourself(), G4Nucleus::AtomicMass(), G4MuonMinusCaptureAtRest::AtRestDoIt(), G4NeutronHPInelasticBaseFS::BaseApply(), G4InuclEvaporation::BreakItUp(), G4BertiniEvaporation::BreakItUp(), G4StopTheoDeexcitation::BreakUp(), G4UnstableFragmentBreakUp::BreakUpFragment(), G4HadronicProcess::CheckEnergyMomentumConservation(), G4HadronicProcess::CheckResult(), G4NeutronHPInelasticCompFS::CompositeApply(), G4Fragment::ComputeGroundStateMass(), G4HadronElastic::ComputeMomentumCMS(), G4PiMinusStopAbsorption::DoAbsorption(), G4PiMinusStopAbsorption::Energy(), G4B9FermiFragment::G4B9FermiFragment(), G4Be8FermiFragment::G4Be8FermiFragment(), G4EvaporationChannel::G4EvaporationChannel(), G4GEMChannel::G4GEMChannel(), G4He5FermiFragment::G4He5FermiFragment(), G4Li5FermiFragment::G4Li5FermiFragment(), G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp(), G4VFermiFragment::G4VFermiFragment(), G4FermiFragmentsPool::GetConfigurationList(), G4NeutronHPCaptureData::GetCrossSection(), G4NeutronHPElasticData::GetCrossSection(), G4NeutronHPFissionData::GetCrossSection(), G4NeutronHPInelasticData::GetCrossSection(), G4NeutronHPorLCaptureData::GetCrossSection(), G4NeutronHPorLEInelasticData::GetCrossSection(), G4NeutronHPorLElasticData::GetCrossSection(), G4NeutronHPorLFissionData::GetCrossSection(), G4IonsKoxCrossSection::GetElementCrossSection(), G4TripathiLightCrossSection::GetElementCrossSection(), G4EvaporationChannel::GetEmissionProbability(), G4GEMChannel::GetEmissionProbability(), G4PromptPhotonEvaporation::GetEmissionProbability(), G4IonsShenCrossSection::GetIsoCrossSection(), G4QMDNucleus::GetNuclearMass(), GetNuclearMass(), G4HyperNucleiProperties::GetNuclearMass(), G4InuclNuclei::getNucleiMass(), G4IonTable::GetNucleusMass(), G4NeutronHPThermalBoost::GetThermalEnergy(), G4VPreCompoundFragment::Initialize(), G4GeneratorPrecompoundInterface::Propagate(), G4AntiNuclElastic::SampleInvariantT(), G4DiffuseElastic::SampleInvariantT(), G4NuclNuclDiffuseElastic::SampleInvariantT(), G4eCoulombScatteringModel::SampleSecondaries(), G4hCoulombScatteringModel::SampleSecondaries(), G4IonCoulombScatteringModel::SampleSecondaries(), G4ScreeningMottCrossSection::SetupKinematic(), G4QElectronNuclearCrossSection::ThresholdEnergy(), G4QMuonNuclearCrossSection::ThresholdEnergy(), G4QTauNuclearCrossSection::ThresholdEnergy(), and G4QPhotonNuclearCrossSection::ThresholdEnergy().

GetNuclearMass() [2/2]

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

static

Definition at line 71 of file G4NucleiProperties.cc.

{
  if (mass_proton  <= 0.0 ) {
    const G4ParticleDefinition * nucleus = 0;
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("proton"); // proton 
    if (nucleus!=0) mass_proton = nucleus->GetPDGMass();
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("neutron"); // neutron 
    if (nucleus!=0) mass_neutron = nucleus->GetPDGMass();
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("deuteron"); // deuteron 
    if (nucleus!=0) mass_deuteron = nucleus->GetPDGMass();
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("triton"); // triton 
    if (nucleus!=0) mass_triton = nucleus->GetPDGMass();
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("alpha"); // alpha 
    if (nucleus!=0) mass_alpha = nucleus->GetPDGMass();
    nucleus = G4ParticleTable::GetParticleTable()->FindParticle("He3"); // He3 
    if (nucleus!=0) mass_He3 = nucleus->GetPDGMass();
 
  }
 
  if (A < 1 || Z < 0 || Z > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>0) {
      G4cerr << "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 (G4NucleiPropertiesTableAME03::IsInTable(Z,A)) {
      // AME 03 table
      mass = G4NucleiPropertiesTableAME03::GetNuclearMass(Z,A);
    } else if (G4NucleiPropertiesTheoreticalTable::IsInTable(Z,A)){
      // Theoretical table
      mass = G4NucleiPropertiesTheoreticalTable::GetNuclearMass(Z,A);
    } 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 134 of file G4NucleiProperties.cc.

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

Referenced by G4IonTable::CreateAllIon(), IsInStableTable(), G4QElectronNuclearCrossSection::ThresholdEnergy(), G4QMuonNuclearCrossSection::ThresholdEnergy(), G4QTauNuclearCrossSection::ThresholdEnergy(), and G4QPhotonNuclearCrossSection::ThresholdEnergy().

IsInStableTable() [2/2]

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

static

Definition at line 141 of file G4NucleiProperties.cc.

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

---

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

• G4NucleiProperties.hh
• G4NucleiProperties.cc