G4HyperNucleiProperties Class Reference

#include <G4HyperNucleiProperties.hh>

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

Static Public Member Functions
static G4double	GetAtomicMass (G4int A, G4int Z, G4int L)
static G4double	GetNuclearMass (G4int A, G4int Z, G4int L)

Detailed Description

Definition at line 44 of file G4HyperNucleiProperties.hh.

Constructor & Destructor Documentation

G4HyperNucleiProperties()

G4HyperNucleiProperties::G4HyperNucleiProperties ( )

default

~G4HyperNucleiProperties()

G4HyperNucleiProperties::~G4HyperNucleiProperties ( )

default

Member Function Documentation

GetAtomicMass()

G4double G4HyperNucleiProperties::GetAtomicMass ( G4int A, G4int Z, G4int L )

static

Definition at line 88 of file G4HyperNucleiProperties.cc.

{
  G4double mass = GetNuclearMass(A, Z, LL);
  if (mass > 0.0) {
    mass += Z * electron_mass_c2 - 1.433e-5 * MeV * G4Pow::GetInstance()->powZ(Z, 2.39);
  }
  return mass;
}

GetNuclearMass()

G4double G4HyperNucleiProperties::GetNuclearMass ( G4int A, G4int Z, G4int L )

static

Definition at line 41 of file G4HyperNucleiProperties.cc.

{
  if (LL == 0) return G4NucleiProperties::GetNuclearMass(A, Z);
 
  if (A < 2 || Z < 0 || Z > A - LL || LL > A) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4HyperNucleiProperties::GetNuclearMass: "
             << " Wrong values for A = " << A << " Z = " << Z << " L = " << LL << G4endl;
    }
#endif
    return 0.0;
  }
  if (A == 2) {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel() > 0) {
      G4cout << "G4HyperNucleiProperties::GetNuclearMass: "
             << " No boud state for A = " << A << " Z = " << Z << " L = " << LL << G4endl;
    }
#endif
    return 0.0;
  }
 
  static const G4double b7 = 25.0 * MeV;
  static const G4double b8 = 10.5;  // Slope
  static const G4double a2 = 0.13 * MeV;  // BindingEnergy for d+Lambda(MeV)
  static const G4double a3 = 2.2 * MeV;  // BindingEnergy for (t/He3)+Lamb(MeV)
  static const G4double eps = 0.0001 * MeV;  // security value (MeV)
 
  G4double mass = G4NucleiProperties::GetNuclearMass(A - LL, Z);  // A non-"strange" nucleus
 
  // Temporarily, the mass of the Lambda below is copied from G4Lambda.cc,
  // but it should be added in CLHEP/Units/PhysicalConstants.h
  const G4double mLL = 1.115683 * CLHEP::GeV;
 
  G4double bs = 0.0;
  if (A - LL == 2)
    bs = a2;  // for nnL,npL,ppL
  else if (A - LL == 3)
    bs = a3;  // for 3nL,2npL,n2pL,3pL
  else if (A - LL > 3)
    bs = b7 * G4Exp(-b8 / (A - LL + 1.0));
  mass += LL * (mLL - bs) + eps;
 
  return mass;
}

Referenced by G4INCLXXInterface::ApplyYourself(), G4Nucleus::AtomicMass(), G4Nucleus::AtomicMass(), G4Fragment::ComputeGroundStateMass(), GetAtomicMass(), G4Fancy3DNucleus::GetMass(), G4IonTable::GetNucleusMass(), and G4GeneratorPrecompoundInterface::PropagateNuclNucl().

---

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

• G4HyperNucleiProperties.hh
• G4HyperNucleiProperties.cc