G4LevelManager Class Reference

#include <G4LevelManager.hh>

Public Member Functions
	G4LevelManager (G4int Z, G4int A, size_t nlev, const std::vector< G4double > &energies, const std::vector< G4int > &spin, const std::vector< const G4NucLevel * > &levels)
	~G4LevelManager ()
size_t	NumberOfTransitions () const
const G4NucLevel *	GetLevel (size_t i) const
G4double	LevelEnergy (size_t i) const
G4double	MaxLevelEnergy () const
size_t	NearestLevelIndex (G4double energy, size_t index=0) const
size_t	NearestLowEdgeLevelIndex (G4double energy) const
const G4NucLevel *	NearestLevel (G4double energy, size_t index=0) const
G4double	NearestLevelEnergy (G4double energy, size_t index=0) const
G4double	NearestLowEdgeLevelEnergy (G4double energy) const
G4double	LifeTime (size_t i) const
G4int	SpinTwo (size_t i) const
G4int	Parity (size_t i) const
G4int	FloatingLevel (size_t i) const
G4double	ShellCorrection () const
G4double	LevelDensity (G4double U) const
const G4String &	FloatingType (size_t i) const
void	StreamInfo (std::ostream &os) const

Detailed Description

Definition at line 53 of file G4LevelManager.hh.

Constructor & Destructor Documentation

G4LevelManager()

G4LevelManager::G4LevelManager ( G4int  Z, G4int  A, size_t  nlev, const std::vector< G4double > &  energies, const std::vector< G4int > &  spin, const std::vector< const G4NucLevel * > &  levels  )

explicit

Definition at line 52 of file G4LevelManager.cc.

  : nTransitions(0)
{ 
  if(0 < ntrans) {
    nTransitions = ntrans - 1;
    fLevelEnergy.reserve(ntrans);
    fSpin.reserve(ntrans);
    fLevels.reserve(ntrans);
    for(size_t i=0; i<ntrans; ++i) {
      fLevelEnergy.push_back(energies[i]);
      fSpin.push_back(spin[i]);
      fLevels.push_back(levels[i]);
    }
    //G4cout << "New G4LevelManager N= " << nTransitions << " " 
    //<< fLevelEnergy.size() << " <" << this << ">" << G4endl;
  }
  // J. Nucl. Sci. Tech. 31(2): 151-162 (1994)
  fShellCorrection = G4NuclearLevelData::GetInstance()->
    GetShellCorrection()->GetShellCorrection(A,Z);
  G4int N = A - Z;
  G4int In = N - (N/2)*2; 
  G4int Iz = Z - (Z/2)*2;
  G4double a13 = 1.0/G4Pow::GetInstance()->Z13(A);
  if(In == 0 && Iz == 0) {
    fLevelDensity = 0.067946*A*(1.0 + 4.1277*a13);
  } else if(In == 0 && Iz == 1) {
    fLevelDensity = 0.053061*A*(1.0 + 7.1862*a13);
  } else if(In == 1 && Iz == 0) {
    fLevelDensity = 0.060920*A*(1.0 + 3.8767*a13);
  } else {
    fLevelDensity = 0.065291*A*(1.0 + 4.4505*a13);
  }
}

~G4LevelManager()

G4LevelManager::~G4LevelManager

(

)

Definition at line 89 of file G4LevelManager.cc.

{
  for(size_t i=0; i<=nTransitions; ++i) { delete fLevels[i]; }
}

Member Function Documentation

FloatingLevel()

G4int G4LevelManager::FloatingLevel ( size_t  i ) const

inline

Definition at line 210 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "Floating"); }
#endif
  return fSpin[i]/100000; 
}

Referenced by StreamInfo().

FloatingType()

const G4String & G4LevelManager::FloatingType

(

size_t

i

)

const

Definition at line 122 of file G4LevelManager.cc.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "FloatingType(idx)"); }
#endif
  return fFloatingLevels[fSpin[i]/100000]; 
}

GetLevel()

const G4NucLevel * G4LevelManager::GetLevel ( size_t  i ) const

inline

Definition at line 138 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "GetLevel(idx)"); }
#endif
  return fLevels[i]; 
}

Referenced by NearestLevel().

LevelDensity()

G4double G4LevelManager::LevelDensity ( G4double  U ) const

inline

Definition at line 223 of file G4LevelManager.hh.

{
  return fLevelDensity;
}

Referenced by G4NuclearLevelData::GetLevelDensity().

LevelEnergy()

G4double G4LevelManager::LevelEnergy ( size_t  i ) const

inline

Definition at line 146 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "LevelEnergy(idx)"); }
#endif
  return fLevelEnergy[i]; 
}

Referenced by G4NuclearLevelData::FindLevel(), NearestLevelEnergy(), and NearestLowEdgeLevelEnergy().

LifeTime()

G4double G4LevelManager::LifeTime ( size_t  i ) const

inline

Definition at line 186 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "LifeTime"); }
#endif
  return (fLevels[i]) ? fLevels[i]->GetTimeGamma() : 0.0;
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), and G4RadioactiveDecay::CalculateChainsFromParent().

MaxLevelEnergy()

G4double G4LevelManager::MaxLevelEnergy ( ) const

inline

Definition at line 154 of file G4LevelManager.hh.

{
  return fLevelEnergy[nTransitions]; 
}

Referenced by G4GEMProbabilityVI::G4GEMProbabilityVI().

NearestLevel()

const G4NucLevel * G4LevelManager::NearestLevel ( G4double  energy, size_t  index = 0  ) const

inline

Definition at line 170 of file G4LevelManager.hh.

{   
  return GetLevel(NearestLevelIndex(energy, index));
}

NearestLevelEnergy()

G4double G4LevelManager::NearestLevelEnergy ( G4double  energy, size_t  index = 0  ) const

inline

Definition at line 176 of file G4LevelManager.hh.

{   
  return LevelEnergy(NearestLevelIndex(energy, index));
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4RadioactiveDecay::CalculateChainsFromParent(), G4PhotonEvaporation::GetFinalLevelEnergy(), and G4NuclearLevelData::GetLevelEnergy().

NearestLevelIndex()

size_t G4LevelManager::NearestLevelIndex	(	G4double	energy,
		size_t	index = 0
	)		const

Definition at line 95 of file G4LevelManager.cc.

{
  //G4cout<< "index= " << index << " max= " << nTransitions << " exc= " << ener 
  //     << " Emax= " << fLevelEnergy[nTransitions] << G4endl;
  size_t idx = std::min(index, nTransitions);
  static const G4double tolerance = 1.0f-6;
  if(0 == nTransitions || std::abs(energy - fLevelEnergy[idx]) <= tolerance) {
    return idx;
  }
  // ground state
  if(energy <= fLevelEnergy[1]*0.5)  
    { idx = 0; }
  // take top level
  else if((fLevelEnergy[nTransitions] + fLevelEnergy[nTransitions-1])*0.5 <= energy) 
    { idx = nTransitions; }
 
  // if shortcuts are not working, make binary search
  else {
    idx = std::lower_bound(fLevelEnergy.begin(), fLevelEnergy.end(), energy)
      - fLevelEnergy.begin() - 1;
    if(energy - fLevelEnergy[idx] > fLevelEnergy[idx+1] - energy) { ++idx; }
    //G4cout << "E= " << energy << " " << fLevelEnergy[idx-1] 
    //<< " " << fLevelEnergy[idx] << G4endl;
  }
  return idx;
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4RadioactiveDecay::CalculateChainsFromParent(), G4NuclearLevelData::FindLevel(), NearestLevel(), and NearestLevelEnergy().

NearestLowEdgeLevelEnergy()

G4double G4LevelManager::NearestLowEdgeLevelEnergy ( G4double  energy ) const

inline

Definition at line 181 of file G4LevelManager.hh.

{   
  return LevelEnergy(NearestLowEdgeLevelIndex(energy));
}

Referenced by G4NuclearLevelData::GetLowEdgeLevelEnergy().

NearestLowEdgeLevelIndex()

size_t G4LevelManager::NearestLowEdgeLevelIndex ( G4double  energy ) const

inline

Definition at line 159 of file G4LevelManager.hh.

{
  size_t idx = nTransitions;
  if(energy < fLevelEnergy[nTransitions]) {
    idx = std::lower_bound(fLevelEnergy.begin(), fLevelEnergy.end(), energy)
      - fLevelEnergy.begin() - 1;
  }
  return idx;
}

Referenced by NearestLowEdgeLevelEnergy().

NumberOfTransitions()

size_t G4LevelManager::NumberOfTransitions ( ) const

inline

Definition at line 133 of file G4LevelManager.hh.

{
  return nTransitions;
}

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4RadioactiveDecay::CalculateChainsFromParent(), G4NuclearLevelData::FindLevel(), G4NuclearLevelData::GetLevelEnergy(), and G4NuclearLevelData::StreamLevels().

Parity()

G4int G4LevelManager::Parity ( size_t  i ) const

inline

Definition at line 202 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "SpinTwo"); }
#endif
  return (fSpin[i]%100000 - 100 > 0) ? 1 : -1; 
}

Referenced by StreamInfo().

ShellCorrection()

G4double G4LevelManager::ShellCorrection ( ) const

inline

Definition at line 218 of file G4LevelManager.hh.

{
  return fShellCorrection;
}

SpinTwo()

G4int G4LevelManager::SpinTwo ( size_t  i ) const

inline

Definition at line 194 of file G4LevelManager.hh.

{
#ifdef G4VERBOSE
  if(i > nTransitions) { PrintError(i, "SpinTwo"); }
#endif
  return std::abs(fSpin[i]%100000 - 100); 
}

Referenced by StreamInfo().

StreamInfo()

void G4LevelManager::StreamInfo

(

std::ostream &

os

)

const

Definition at line 141 of file G4LevelManager.cc.

{
  for(size_t i=0; i<=nTransitions; ++i) {
    G4int prec = out.precision(6);
    out << std::setw(6) << i << ". " 
    << std::setw(8) << fLevelEnergy[i];
    if(fLevels[i]) {
    out << std::setw(8) << fLevels[i]->GetTimeGamma()
        << std::setw(4) << fLevels[i]->NumberOfTransitions()
        << std::setw(4) << SpinTwo(i)
        << std::setw(4) << Parity(i)
        << std::setw(4) << FloatingLevel(i);
    }
    out << "\n";
    out.precision(prec);
    if(fLevels[i]) { fLevels[i]->StreamInfo(out); }
  }
}

Referenced by G4NuclearLevelData::StreamLevels().

---

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

• G4LevelManager.hh
• G4LevelManager.cc