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 45 of file G4HyperNucleiProperties.hh.

Constructor & Destructor Documentation

G4HyperNucleiProperties::~G4HyperNucleiProperties ( )

inline

Definition at line 56 of file G4HyperNucleiProperties.hh.

{ };

G4HyperNucleiProperties::G4HyperNucleiProperties ( )

inline

Definition at line 59 of file G4HyperNucleiProperties.hh.

{};

Member Function Documentation

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

static

Definition at line 102 of file G4HyperNucleiProperties.cc.

References python.hepunit::electron_mass_c2, G4cout, G4endl, GetNuclearMass(), G4ParticleTable::GetParticleTable(), and python.hepunit::MeV.

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

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

static

Definition at line 42 of file G4HyperNucleiProperties.cc.

References G4ParticleTable::FindParticle(), G4cout, G4endl, G4NucleiProperties::GetNuclearMass(), G4ParticleTable::GetParticleTable(), G4ParticleDefinition::GetPDGMass(), G4ParticleTable::GetVerboseLevel(), G4InuclParticleNames::lambda, and python.hepunit::MeV.

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

{
  if (L==0) return  G4NucleiProperties::GetNuclearMass(A, Z);
  
  G4ParticleTable* pTable = G4ParticleTable::GetParticleTable();
  if (A < 2 || Z < 0 || Z > A-L || L>A ) {
#ifdef G4VERBOSE
    if (pTable->GetVerboseLevel()>0) {
      G4cout << "G4HyperNucleiProperties::GetNuclearMass: "
         << " Wrong values for A = " << A 
         << " Z = " << Z 
         << " L = " << L
         << 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 = " << L
         << 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 mL= 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-L, Z); 
  // A non-"strange" nucleus
  G4double bs=0.;
  if     (A-L ==2) bs=a2;         // for nnL,npL,ppL
  else if(A-L ==3) bs=a3;         // for 3nL,2npL,n2pL,3pL
  else if(A-L >3)  bs=b7*std::exp(-b8/(A-L+1.));
  mass += L*(mL-bs) + eps;

  return mass;
}

---

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

• G4HyperNucleiProperties.hh
• G4HyperNucleiProperties.cc