G4NucleiProperties Class Reference

#include <G4NucleiProperties.hh>

Collaboration diagram for G4NucleiProperties:

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)

Private Types
enum	{ MaxZ = 120 }

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

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

Detailed Description

Definition at line 50 of file G4NucleiProperties.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
MaxZ

Definition at line 101 of file G4NucleiProperties.hh.

{MaxZ = 120};

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.

{};

Here is the call graph for this function:

Member Function Documentation

AtomicMass()

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

staticprivate

Definition at line 252 of file G4NucleiProperties.cc.

{
  const G4double hydrogen_mass_excess = G4NucleiPropertiesTableAME03::GetMassExcess(1,1);
  const G4double neutron_mass_excess =  G4NucleiPropertiesTableAME03::GetMassExcess(0,1);

  G4double mass =
      (A-Z)*neutron_mass_excess + Z*hydrogen_mass_excess - BindingEnergy(A,Z) + A*amu_c2;

  return mass;
}

Here is the call graph for this function:

Here is the caller graph for this function:

BindingEnergy()

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

staticprivate

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

Here is the caller graph for this function:

GetAtomicMass()

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

staticprivate

Definition at line 189 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 (G4NucleiPropertiesTableAME03::IsInTable(iZ,iA)) {
      return G4NucleiPropertiesTableAME03::GetAtomicMass(iZ,iA);
    } else if (G4NucleiPropertiesTheoreticalTable::IsInTable(iZ,iA)){
      return G4NucleiPropertiesTheoreticalTable::GetAtomicMass(iZ,iA);
    } else {
      return AtomicMass(A,Z);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetBindingEnergy() [1/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) {
      G4cout << "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);
    }

  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetBindingEnergy() [2/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);
}

Here is the call graph for this function:

GetMassExcess() [1/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) {
      G4cout << "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);
    }
  }

}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMassExcess() [2/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);
}

Here is the call graph for this function:

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

Here is the call graph for this function:

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

Here is the call graph for this function:

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

Here is the caller graph for this function:

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) {
      G4cout << "G4NucleiProperties::IsInStableTable: Wrong values for A = " 
         << A << " and Z = " << Z << G4endl;    
    }
#endif 
    return false;
  } 

  return G4NucleiPropertiesTableAME03::IsInTable(Z,A);

}

Here is the call graph for this function:

MassExcess()

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

staticprivate

Definition at line 247 of file G4NucleiProperties.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

NuclearMass()

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

staticprivate

Definition at line 263 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 formula in  AME03
  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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

electronMass

G4ThreadLocal G4double G4NucleiProperties::electronMass[MaxZ]

staticprivate

Definition at line 102 of file G4NucleiProperties.hh.

isIntialized

G4ThreadLocal G4bool G4NucleiProperties::isIntialized

staticprivate

Definition at line 105 of file G4NucleiProperties.hh.

mass_alpha

G4ThreadLocal G4double G4NucleiProperties::mass_alpha = -1.

staticprivate

Definition at line 110 of file G4NucleiProperties.hh.

mass_deuteron

G4ThreadLocal G4double G4NucleiProperties::mass_deuteron = -1.

staticprivate

Definition at line 108 of file G4NucleiProperties.hh.

mass_He3

G4ThreadLocal G4double G4NucleiProperties::mass_He3 = -1.

staticprivate

Definition at line 111 of file G4NucleiProperties.hh.

mass_neutron

G4ThreadLocal G4double G4NucleiProperties::mass_neutron = -1.

staticprivate

Definition at line 107 of file G4NucleiProperties.hh.

mass_proton

G4ThreadLocal G4double G4NucleiProperties::mass_proton = -1.

staticprivate

Definition at line 106 of file G4NucleiProperties.hh.

mass_triton

G4ThreadLocal G4double G4NucleiProperties::mass_triton = -1.

staticprivate

Definition at line 109 of file G4NucleiProperties.hh.

---

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

• G4NucleiProperties.hh
• G4NucleiProperties.cc