G4HyperNucleiProperties Class Reference

#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.

{}

~G4HyperNucleiProperties()

G4HyperNucleiProperties::~G4HyperNucleiProperties ( )

inline

Definition at line 50 of file G4HyperNucleiProperties.hh.

{}

Member Function Documentation

GetAtomicMass()

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

static

Definition at line 105 of file G4HyperNucleiProperties.cc.

{
  if (A < 1 || Z < 0 || Z > A-LL || LL > A || LL <0 )
  {
#ifdef G4VERBOSE
    if (G4ParticleTable::GetParticleTable()->GetVerboseLevel()>0)
    {
      G4cout << "G4HyperNucleiProperties::GetAtomicMass: " 
             << " Wrong values for A = "  << A 
             << "  Z = " << Z 
             << "  L = " << LL  << G4endl;
    }
#endif
    return 0.0;
    
  }
  else
  {
    G4double nuclearMass = GetNuclearMass(A, Z, LL);
    return nuclearMass + Z*electron_mass_c2
         - 1.433e-5*MeV*std::pow(G4double(Z),2.39);
  }
}

GetNuclearMass()

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

static

Definition at line 38 of file G4HyperNucleiProperties.cc.

{
  if (LL==0) return  G4NucleiProperties::GetNuclearMass(A, Z);
  
  G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
 
  if (A < 2 || Z < 0 || Z > A-LL || LL>A )
  {
#ifdef G4VERBOSE
    if (pTable->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 (pTable->GetVerboseLevel()>0)
    {
      G4cout << "G4HyperNucleiProperties::GetNuclearMass: "
             << " No boud state for A = " << A 
             << " Z = " << Z 
             << " L = " << LL
             << G4endl;
    }
#endif
    return 0.0;
  }
 
  
  G4ParticleDefinition* lambda = pTable->FindParticle("lambda");
  if (lambda ==0)
  {
#ifdef G4VERBOSE
    if (pTable->GetVerboseLevel()>0)
    {
      G4cout << "G4HyperNucleiProperties::GetNuclearMass: "
             << " Lambda is not defined " << G4endl;
    }
#endif
    return 0.0;
  }
  const G4double mLL= lambda->GetPDGMass(); // mLambda
 
  static const G4double b7=25.*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
  G4double bs=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*std::exp(-b8/(A-LL+1.));
  mass += LL*(mLL-bs) + eps;
 
  return mass;
}

Referenced by GetAtomicMass(), and G4IonTable::GetNucleusMass().

---

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

• G4HyperNucleiProperties.hh
• G4HyperNucleiProperties.cc