G4LevelManager Class Reference

#include <G4LevelManager.hh>

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

Private Member Functions
	G4LevelManager (const G4LevelManager &right)=delete
G4bool	operator!= (const G4LevelManager &right) const =delete
const G4LevelManager &	operator= (const G4LevelManager &right)=delete
G4bool	operator== (const G4LevelManager &right) const =delete

Private Attributes
G4double	fLevelDensity
std::vector< G4double >	fLevelEnergy
std::vector< const G4NucLevel * >	fLevels
G4double	fShellCorrection
std::vector< G4int >	fSpin
size_t	nTransitions

Static Private Attributes
static G4String	fFloatingLevels [nfloting]
static const G4int	nfloting = 13

Detailed Description

Definition at line 53 of file G4LevelManager.hh.

Constructor & Destructor Documentation

G4LevelManager() [1/2]

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);
  }
}

References A, fLevelDensity, fLevelEnergy, fLevels, fShellCorrection, fSpin, G4Pow::GetInstance(), G4NuclearLevelData::GetInstance(), G4NuclearLevelData::GetShellCorrection(), nTransitions, Z, and G4Pow::Z13().

~G4LevelManager()

G4LevelManager::~G4LevelManager

(

)

Definition at line 89 of file G4LevelManager.cc.

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

References fLevels, and nTransitions.

G4LevelManager() [2/2]

G4LevelManager::G4LevelManager ( const G4LevelManager &  right )

privatedelete

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; 
}

References fSpin, and nTransitions.

Referenced by G4PhotonEvaporation::GenerateGamma(), and 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]; 
}

References fFloatingLevels, fSpin, and nTransitions.

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]; 
}

References fLevels, and nTransitions.

Referenced by G4PhotonEvaporation::GenerateGamma(), and NearestLevel().

LevelDensity()

G4double G4LevelManager::LevelDensity ( G4double  U ) const

inline

Definition at line 223 of file G4LevelManager.hh.

{
  return fLevelDensity;
}

References 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]; 
}

References fLevelEnergy, and nTransitions.

Referenced by G4NuclearLevelData::FindLevel(), G4PhotonEvaporation::GenerateGamma(), G4FermiFragmentsPoolVI::Initialise(), 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;
}

References fLevels, and nTransitions.

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4PhotonEvaporation::GenerateGamma(), and G4FermiFragmentsPoolVI::Initialise().

MaxLevelEnergy()

G4double G4LevelManager::MaxLevelEnergy ( ) const

inline

Definition at line 154 of file G4LevelManager.hh.

{
  return fLevelEnergy[nTransitions]; 
}

References fLevelEnergy, and nTransitions.

Referenced by G4GEMProbabilityVI::G4GEMProbabilityVI(), and G4PhotonEvaporation::InitialiseLevelManager().

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));
}

References G4INCL::KinematicsUtils::energy(), GetLevel(), and NearestLevelIndex().

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));
}

References G4INCL::KinematicsUtils::energy(), LevelEnergy(), and NearestLevelIndex().

Referenced by G4Radioactivation::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;
}

References G4INCL::KinematicsUtils::energy(), fLevelEnergy, G4INCL::Math::min(), and nTransitions.

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4NuclearLevelData::FindLevel(), G4VEmissionProbability::FindRecoilExcitation(), G4PhotonEvaporation::GenerateGamma(), NearestLevel(), and NearestLevelEnergy().

NearestLowEdgeLevelEnergy()

G4double G4LevelManager::NearestLowEdgeLevelEnergy ( G4double  energy ) const

inline

Definition at line 181 of file G4LevelManager.hh.

{   
  return LevelEnergy(NearestLowEdgeLevelIndex(energy));
}

References G4INCL::KinematicsUtils::energy(), LevelEnergy(), and NearestLowEdgeLevelIndex().

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;
}

References G4INCL::KinematicsUtils::energy(), fLevelEnergy, and nTransitions.

Referenced by NearestLowEdgeLevelEnergy().

NumberOfTransitions()

size_t G4LevelManager::NumberOfTransitions ( ) const

inline

Definition at line 133 of file G4LevelManager.hh.

{
  return nTransitions;
}

References nTransitions.

Referenced by G4Radioactivation::CalculateChainsFromParent(), G4NuclearLevelData::FindLevel(), G4PhotonEvaporation::GenerateGamma(), G4NuclearLevelData::GetLevelEnergy(), G4FermiFragmentsPoolVI::Initialise(), and G4NuclearLevelData::StreamLevels().

operator!=()

G4bool G4LevelManager::operator!= ( const G4LevelManager &  right ) const

privatedelete

operator=()

const G4LevelManager & G4LevelManager::operator= ( const G4LevelManager &  right )

privatedelete

operator==()

G4bool G4LevelManager::operator== ( const G4LevelManager &  right ) const

privatedelete

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; 
}

References fSpin, and nTransitions.

Referenced by StreamInfo().

ShellCorrection()

G4double G4LevelManager::ShellCorrection ( ) const

inline

Definition at line 218 of file G4LevelManager.hh.

{
  return fShellCorrection;
}

References 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); 
}

References fSpin, and nTransitions.

Referenced by G4PhotonEvaporation::GenerateGamma(), G4FermiFragmentsPoolVI::Initialise(), and 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); }
  }
}

References fLevelEnergy, fLevels, FloatingLevel(), nTransitions, Parity(), CLHEP::prec, and SpinTwo().

Referenced by G4NuclearLevelData::StreamLevels().

Field Documentation

fFloatingLevels

G4String G4LevelManager::fFloatingLevels

staticprivate

Initial value:

= {
  "-", "+X", "+Y", "+Z", "+U", "+V", "+W", "+R", "+S", "+T", "+A", "+B", "+C"}

Definition at line 129 of file G4LevelManager.hh.

Referenced by FloatingType().

fLevelDensity

G4double G4LevelManager::fLevelDensity

private

Definition at line 124 of file G4LevelManager.hh.

Referenced by G4LevelManager(), and LevelDensity().

fLevelEnergy

std::vector<G4double> G4LevelManager::fLevelEnergy

private

Definition at line 119 of file G4LevelManager.hh.

Referenced by G4LevelManager(), LevelEnergy(), MaxLevelEnergy(), NearestLevelIndex(), NearestLowEdgeLevelIndex(), and StreamInfo().

fLevels

std::vector<const G4NucLevel*> G4LevelManager::fLevels

private

Definition at line 121 of file G4LevelManager.hh.

Referenced by G4LevelManager(), GetLevel(), LifeTime(), StreamInfo(), and ~G4LevelManager().

fShellCorrection

G4double G4LevelManager::fShellCorrection

private

Definition at line 123 of file G4LevelManager.hh.

Referenced by G4LevelManager(), and ShellCorrection().

fSpin

std::vector<G4int> G4LevelManager::fSpin

private

Definition at line 120 of file G4LevelManager.hh.

Referenced by FloatingLevel(), FloatingType(), G4LevelManager(), Parity(), and SpinTwo().

nfloting

const G4int G4LevelManager::nfloting = 13

staticprivate

Definition at line 128 of file G4LevelManager.hh.

nTransitions

size_t G4LevelManager::nTransitions

private

Definition at line 126 of file G4LevelManager.hh.

Referenced by FloatingLevel(), FloatingType(), G4LevelManager(), GetLevel(), LevelEnergy(), LifeTime(), MaxLevelEnergy(), NearestLevelIndex(), NearestLowEdgeLevelIndex(), NumberOfTransitions(), Parity(), SpinTwo(), StreamInfo(), and ~G4LevelManager().

---

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

• source/processes/hadronic/models/de_excitation/management/include/G4LevelManager.hh
• source/processes/hadronic/models/de_excitation/management/src/G4LevelManager.cc