#include <G4HyperNucleiProperties.hh>

Public Member Functions
	G4HyperNucleiProperties ()

	~G4HyperNucleiProperties ()

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

inline

Definition at line 49 of file G4HyperNucleiProperties.hh.

49{}

◆ ~G4HyperNucleiProperties()

G4HyperNucleiProperties::~G4HyperNucleiProperties ( )

inline

Definition at line 50 of file G4HyperNucleiProperties.hh.

50{}

Member Function Documentation

◆ GetAtomicMass()

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

static

Definition at line 96 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;
  }
  else 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(), G4Fragment::ComputeGroundStateMass(), GetAtomicMass(), G4Fancy3DNucleus::GetMass(), G4IonTable::GetNucleusMass(), and G4GeneratorPrecompoundInterface::PropagateNuclNucl().

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

Public Member Functions