G4HyperNucleiProperties Class Reference

#include <G4HyperNucleiProperties.hh>

Collaboration diagram for G4HyperNucleiProperties:

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::~G4HyperNucleiProperties ( )

inline

Definition at line 56 of file G4HyperNucleiProperties.hh.

{ };

G4HyperNucleiProperties()

G4HyperNucleiProperties::G4HyperNucleiProperties ( )

inline

Definition at line 59 of file G4HyperNucleiProperties.hh.

{};

Here is the call graph for this function:

Member Function Documentation

GetAtomicMass()

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

static

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

Here is the call graph for this function:

Here is the caller graph for this function:

GetNuclearMass()

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

static

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4HyperNucleiProperties.hh
• G4HyperNucleiProperties.cc