G4NuclearLevelManager Class Reference

#include <G4NuclearLevelManager.hh>

Collaboration diagram for G4NuclearLevelManager:

Public Member Functions
	G4NuclearLevelManager (G4int Z, G4int A, const G4String &filename)
	G4NuclearLevelManager (const G4NuclearLevelManager &right)
	~G4NuclearLevelManager ()
void	SetNucleus (G4int Z, G4int A, const G4String &filename)
G4bool	IsValid () const
G4int	NumberOfLevels () const
const G4NuclearLevel *	GetLevel (G4int i) const
const G4NuclearLevel *	NearestLevel (G4double energy, G4double eDiffMax=1.e+8) const
const G4NuclearLevel *	LowestLevel () const
const G4NuclearLevel *	HighestLevel () const
G4double	MinLevelEnergy () const
G4double	MaxLevelEnergy () const
void	PrintLevels ()
void	PrintAll ()

Private Member Functions
const G4NuclearLevelManager &	operator= (const G4NuclearLevelManager &right)
G4bool	operator== (const G4NuclearLevelManager &right) const
G4bool	operator!= (const G4NuclearLevelManager &right) const
G4bool	Read (std::ifstream &aDataFile)
G4bool	ReadDataLine (std::ifstream &dataFile)
G4bool	ReadDataItem (std::istream &dataFile, G4double &x)
void	ProcessDataLine ()
void	MakeLevels (const G4String &filename)
void	ClearLevels ()
G4NuclearLevel *	UseLevelOrMakeNew (G4NuclearLevel *level)
void	AddDataToLevel (G4NuclearLevel *level)
void	FinishLevel (G4NuclearLevel *level)

Private Attributes
G4int	_nucleusA
G4int	_nucleusZ
G4bool	_validity
G4PtrLevelVector *	_levels

Static Private Attributes
static G4double	_levelEnergy =0.
static G4double	_gammaEnergy =0.
static G4double	_probability =0.
static G4double	_polarity =0.
static G4double	_halfLife =0.
static G4double	_angularMomentum =0.
static G4double	_kCC =0.
static G4double	_l1CC =0.
static G4double	_l2CC =0.
static G4double	_l3CC =0.
static G4double	_m1CC =0.
static G4double	_m2CC =0.
static G4double	_m3CC =0.
static G4double	_m4CC =0.
static G4double	_m5CC =0.
static G4double	_nPlusCC =0.
static G4double	_totalCC =0.

Detailed Description

Definition at line 65 of file G4NuclearLevelManager.hh.

Constructor & Destructor Documentation

G4NuclearLevelManager() [1/2]

G4NuclearLevelManager::G4NuclearLevelManager	(	G4int	Z,
		G4int	A,
		const G4String &	filename
	)

Definition at line 82 of file G4NuclearLevelManager.cc.

                                                   :
    _nucleusA(A), _nucleusZ(Z), _validity(false), 
    _levels(0)
{ 
  if (A <= 0 || Z <= 0 || Z > A ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                  "==== G4NuclearLevelManager ==== (Z,A) <0, or Z>A");
  }
  MakeLevels(filename);
  /*
  G4cout << "+++++ New G4NuclearLevelManager Z= " << Z << " A= " << A 
     << "  " << _levels;
  if(_levels) { G4cout << " Nlevels= " << _levels->size(); }
  G4cout << G4endl;
  */
}

Here is the call graph for this function:

G4NuclearLevelManager() [2/2]

G4NuclearLevelManager::G4NuclearLevelManager

(

const G4NuclearLevelManager &

right

)

Definition at line 406 of file G4NuclearLevelManager.cc.

{
  _nucleusA = right._nucleusA;
  _nucleusZ = right._nucleusZ;
  _validity = right._validity;

  if (right._levels != 0)   
    {
      _levels = new G4PtrLevelVector;
      G4int n = right._levels->size();
      G4int i;
      for (i=0; i<n; ++i)
    {
      _levels->push_back(new G4NuclearLevel(*(right.GetLevel(i))));
    }
      G4PtrSort<G4NuclearLevel>(_levels);
    }
  else 
    {
      _levels = 0;
    }
}

Here is the call graph for this function:

~G4NuclearLevelManager()

G4NuclearLevelManager::~G4NuclearLevelManager

(

)

Definition at line 100 of file G4NuclearLevelManager.cc.

{
  ClearLevels();
}

Here is the call graph for this function:

Member Function Documentation

AddDataToLevel()

void G4NuclearLevelManager::AddDataToLevel ( G4NuclearLevel *  level )

private

Definition at line 346 of file G4NuclearLevelManager.cc.

{
  if (!level) return;       // Sanity check

  level->_energies.push_back(_gammaEnergy);
  level->_weights.push_back(_probability);
  level->_polarities.push_back(_polarity);
  level->_kCC.push_back(_kCC);
  level->_l1CC.push_back(_l1CC);
  level->_l2CC.push_back(_l2CC);
  level->_l3CC.push_back(_l3CC);
  level->_m1CC.push_back(_m1CC);
  level->_m2CC.push_back(_m2CC);
  level->_m3CC.push_back(_m3CC);
  level->_m4CC.push_back(_m4CC);
  level->_m5CC.push_back(_m5CC);
  level->_nPlusCC.push_back(_nPlusCC);
  level->_totalCC.push_back(_totalCC);
}

Here is the caller graph for this function:

ClearLevels()

void G4NuclearLevelManager::ClearLevels ( )

private

Definition at line 278 of file G4NuclearLevelManager.cc.

{
  _validity = false;

  if (NumberOfLevels() > 0) {
    std::for_each(_levels->begin(), _levels->end(), DeleteLevel());
    _levels->clear();
  }

  delete _levels;
  _levels = 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

FinishLevel()

void G4NuclearLevelManager::FinishLevel ( G4NuclearLevel *  level )

private

Definition at line 366 of file G4NuclearLevelManager.cc.

{
  if (!level || !_levels) return;       // Sanity check

  level->Finalize();
  _levels->push_back(level);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetLevel()

const G4NuclearLevel * G4NuclearLevelManager::GetLevel

(

G4int

i

)

const

Definition at line 119 of file G4NuclearLevelManager.cc.

                                                                   {
  return (i>=0 && i<NumberOfLevels()) ? (*_levels)[i] : 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

HighestLevel()

const G4NuclearLevel * G4NuclearLevelManager::HighestLevel

(

)

const

Definition at line 159 of file G4NuclearLevelManager.cc.

{
  return (NumberOfLevels() > 0) ? _levels->front() : 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

IsValid()

G4bool G4NuclearLevelManager::IsValid ( void  ) const

inline

Definition at line 77 of file G4NuclearLevelManager.hh.

{ return _validity; }

LowestLevel()

const G4NuclearLevel * G4NuclearLevelManager::LowestLevel

(

)

const

Definition at line 164 of file G4NuclearLevelManager.cc.

{
  return (NumberOfLevels() > 0) ? _levels->back() : 0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

MakeLevels()

void G4NuclearLevelManager::MakeLevels ( const G4String &  filename )

private

Definition at line 291 of file G4NuclearLevelManager.cc.

{
  _validity = false;
  if (NumberOfLevels() > 0) ClearLevels();  // Discard existing data

  std::ifstream inFile(filename, std::ios::in);
  if (! inFile) 
    {
#ifdef GAMMAFILEWARNING
      if (_nucleusZ > 10) G4cout << " G4NuclearLevelManager: nuclide (" 
                                 << _nucleusZ << "," << _nucleusA 
                                 << ") does not have a gamma levels file" << G4endl;
#endif
      return;
    }

  _levels = new G4PtrLevelVector;

  // Read individual gamma data and fill levels for this nucleus
 
  G4NuclearLevel* thisLevel = 0;
  G4int nData = 0;

  // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  while (Read(inFile)) {
    thisLevel = UseLevelOrMakeNew(thisLevel);   // May create new pointer
    AddDataToLevel(thisLevel);
    nData++;                    // For debugging purposes
  }

  FinishLevel(thisLevel);       // Final  must be completed by hand
  
  // ---- MGP ---- Don't forget to close the file 
  inFile.close();
    
  //  G4cout << " ==== MakeLevels ===== " << nData << " data read " << G4endl;

  G4PtrSort<G4NuclearLevel>(_levels);
  
  _validity = true;
  
  return;
}

Here is the call graph for this function:

Here is the caller graph for this function:

MaxLevelEnergy()

G4double G4NuclearLevelManager::MaxLevelEnergy

(

)

const

Definition at line 154 of file G4NuclearLevelManager.cc.

{
  return (NumberOfLevels() > 0) ? _levels->back()->Energy() : 0.*GeV;
}

Here is the call graph for this function:

Here is the caller graph for this function:

MinLevelEnergy()

G4double G4NuclearLevelManager::MinLevelEnergy

(

)

const

Definition at line 149 of file G4NuclearLevelManager.cc.

{
  return (NumberOfLevels() > 0) ? _levels->front()->Energy() : 9999.*GeV;
}

Here is the call graph for this function:

Here is the caller graph for this function:

NearestLevel()

const G4NuclearLevel * G4NuclearLevelManager::NearestLevel	(	G4double	energy,
		G4double	eDiffMax = 1.e+8
	)		const

Definition at line 124 of file G4NuclearLevelManager.cc.

{
  if (NumberOfLevels() <= 0) return 0;

  G4int iNear = -1;

  //G4cout << "G4NuclearLevelManager::NearestLevel E(MeV)= " 
  //     << energy/MeV << " dEmax(MeV)= " << eDiffMax/MeV << G4endl;
  
  G4double diff = 1.e+10;
  for (unsigned int i=0; i<_levels->size(); ++i)
    {
      G4double e = GetLevel(i)->Energy();
      G4double eDiff = std::fabs(e - energy);
      //G4cout << i << ".   eDiff(MeV)= " << eDiff/MeV << G4endl;
      if (eDiff <= diff)
    { 
      diff = eDiff; 
      iNear = i;
    }
    }
  
  return GetLevel(iNear);   // Includes range checking on iNear
}

Here is the call graph for this function:

Here is the caller graph for this function:

NumberOfLevels()

G4int G4NuclearLevelManager::NumberOfLevels ( ) const

inline

Definition at line 79 of file G4NuclearLevelManager.hh.

{ return (_levels ? _levels->size() : 0); }

Here is the call graph for this function:

Here is the caller graph for this function:

operator!=()

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

private

Here is the caller graph for this function:

operator=()

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

private

Here is the caller graph for this function:

operator==()

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

private

Here is the caller graph for this function:

PrintAll()

void G4NuclearLevelManager::PrintAll

(

void

)

Definition at line 375 of file G4NuclearLevelManager.cc.

{
  G4int nLevels = NumberOfLevels();
    
  G4cout << " ==== G4NuclearLevelManager ==== (" << _nucleusZ << ", " << _nucleusA
     << ") has " << nLevels << " levels" << G4endl
     << "Highest level is at energy " << MaxLevelEnergy() << " MeV "
     << G4endl << "Lowest level is at energy " << MinLevelEnergy()
     << " MeV " << G4endl;
    
  for (G4int i=0; i<nLevels; ++i) {
    GetLevel(i)->PrintAll();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

PrintLevels()

void G4NuclearLevelManager::PrintLevels

(

)

Definition at line 390 of file G4NuclearLevelManager.cc.

{
  G4int nLevels = NumberOfLevels();
  G4double efermi = G4NucleiProperties::GetNuclearMass(_nucleusA-1, _nucleusZ)
    + neutron_mass_c2 -
    G4NucleiProperties::GetNuclearMass(_nucleusA, _nucleusZ);

  G4cout << "Z= " << _nucleusZ << " A= " << _nucleusA
     << "  " << nLevels << " levels" 
     << "  Efermi(MeV)= " << efermi << G4endl;
    
  for (G4int i=0; i<nLevels; ++i) {
    GetLevel(i)->PrintLevels();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

ProcessDataLine()

void G4NuclearLevelManager::ProcessDataLine ( )

private

Definition at line 227 of file G4NuclearLevelManager.cc.

{
  const G4double minProbability = 1e-8;
  
  // Assign units for dimensional quantities
  _levelEnergy *= keV;
  _gammaEnergy *= keV;
  _halfLife *= second;
  
  // The following adjustment is needed to take care of anomalies in 
  // data files, where some transitions show up with relative probability
  // zero
  if (_probability < minProbability) _probability = minProbability;
  // the folowwing is to convert icc probability to accumulative ones
  _l1CC += _kCC;
  _l2CC += _l1CC;
  _l3CC += _l2CC;
  _m1CC += _l3CC;
  _m2CC += _m1CC;
  _m3CC += _m2CC;
  _m4CC += _m3CC;
  _m5CC += _m4CC;
  _nPlusCC += _m5CC;

  if (_nPlusCC!=0) {    // Normalize to probabilities
    _kCC /= _nPlusCC;
    _l1CC /= _nPlusCC;
    _l2CC /= _nPlusCC;
    _l3CC /= _nPlusCC;
    _m1CC /= _nPlusCC;
    _m2CC /= _nPlusCC;
    _m3CC /= _nPlusCC;
    _m4CC /= _nPlusCC;
    _m5CC /= _nPlusCC;
    _nPlusCC /= _nPlusCC;  
  } else {      // Total was zero, reset to unity
    _kCC = 1;
    _l1CC = 1;
    _l2CC = 1;
    _l3CC = 1;
    _m1CC = 1;
    _m2CC = 1;
    _m3CC = 1;
    _m4CC = 1;
    _m5CC = 1;
    _nPlusCC = 1;
  }
    
  // G4cout << "Read " << _levelEnergy << " " << _gammaEnergy << " " << _probability << G4endl;
}

Here is the caller graph for this function:

Read()

G4bool G4NuclearLevelManager::Read ( std::ifstream &  aDataFile )

private

Definition at line 169 of file G4NuclearLevelManager.cc.

{
  G4bool goodRead = ReadDataLine(dataFile);
  
  if (goodRead) ProcessDataLine();
  return goodRead;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReadDataItem()

G4bool G4NuclearLevelManager::ReadDataItem ( std::istream &  dataFile, G4double &  x  )

private

Definition at line 213 of file G4NuclearLevelManager.cc.

{
  // G4bool okay = (dataFile >> buffer) != 0;       // Get next token
  // if (okay) x = strtod(buffer, NULL);
  char buffer[30];
  for(G4int i=0; i<30; ++i) { buffer[i] = 0; }
  G4bool okay = true;
  dataFile >> buffer;
  if(dataFile.fail()) { okay = false; }
  else { x = strtod(buffer, NULL); }

  return okay;
}

Here is the caller graph for this function:

ReadDataLine()

G4bool G4NuclearLevelManager::ReadDataLine ( std::ifstream &  dataFile )

private

Definition at line 179 of file G4NuclearLevelManager.cc.

                                                                {
  /***** DO NOT USE REGULAR STREAM I/O
  G4bool result = true;
  if (dataFile >> _levelEnergy)
    {
      dataFile >> _gammaEnergy >> _probability >> _polarity >> _halfLife
           >> _angularMomentum  >> _totalCC >> _kCC >> _l1CC >> _l2CC 
           >> _l3CC >> _m1CC >> _m2CC >> _m3CC >> _m4CC >> _m5CC
           >> _nPlusCC;
    }
  else result = false;
  *****/

  // Each item will return iostream status
  return (ReadDataItem(dataFile, _levelEnergy) &&
      ReadDataItem(dataFile, _gammaEnergy) &&
      ReadDataItem(dataFile, _probability) &&
      ReadDataItem(dataFile, _polarity) &&
      ReadDataItem(dataFile, _halfLife) &&
      ReadDataItem(dataFile, _angularMomentum) &&
      ReadDataItem(dataFile, _totalCC) &&
      ReadDataItem(dataFile, _kCC) &&
      ReadDataItem(dataFile, _l1CC) &&
      ReadDataItem(dataFile, _l2CC) &&
      ReadDataItem(dataFile, _l3CC) &&
      ReadDataItem(dataFile, _m1CC) &&
      ReadDataItem(dataFile, _m2CC) &&
      ReadDataItem(dataFile, _m3CC) &&
      ReadDataItem(dataFile, _m4CC) &&
      ReadDataItem(dataFile, _m5CC) &&
      ReadDataItem(dataFile, _nPlusCC) );
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetNucleus()

void G4NuclearLevelManager::SetNucleus	(	G4int	Z,
		G4int	A,
		const G4String &	filename
	)

Definition at line 105 of file G4NuclearLevelManager.cc.

{
  if (A <= 0 || Z <= 0 || Z > A ) {
    throw G4HadronicException(__FILE__, __LINE__, 
                  "==== G4NuclearLevelManager ==== (Z,A) <0, or Z>A");
  }
  if (_nucleusZ != Z || _nucleusA != A)
    {
      _nucleusA = A;
      _nucleusZ = Z;
      MakeLevels(filename);
    }
}

Here is the call graph for this function:

UseLevelOrMakeNew()

G4NuclearLevel * G4NuclearLevelManager::UseLevelOrMakeNew ( G4NuclearLevel *  level )

private

Definition at line 336 of file G4NuclearLevelManager.cc.

{
  if (level && _levelEnergy == level->Energy()) return level;   // No change

  if (level) FinishLevel(level);    // Save what we have up to now

  //  G4cout << "Making a new level... " << _levelEnergy << G4endl;
  return new G4NuclearLevel(_levelEnergy, _halfLife, _angularMomentum);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

_angularMomentum

G4double G4NuclearLevelManager::_angularMomentum =0.

staticprivate

Definition at line 126 of file G4NuclearLevelManager.hh.

_gammaEnergy

G4double G4NuclearLevelManager::_gammaEnergy =0.

staticprivate

Definition at line 122 of file G4NuclearLevelManager.hh.

_halfLife

G4double G4NuclearLevelManager::_halfLife =0.

staticprivate

Definition at line 125 of file G4NuclearLevelManager.hh.

_kCC

G4double G4NuclearLevelManager::_kCC =0.

staticprivate

Definition at line 127 of file G4NuclearLevelManager.hh.

_l1CC

G4double G4NuclearLevelManager::_l1CC =0.

staticprivate

Definition at line 128 of file G4NuclearLevelManager.hh.

_l2CC

G4double G4NuclearLevelManager::_l2CC =0.

staticprivate

Definition at line 129 of file G4NuclearLevelManager.hh.

_l3CC

G4double G4NuclearLevelManager::_l3CC =0.

staticprivate

Definition at line 130 of file G4NuclearLevelManager.hh.

_levelEnergy

G4double G4NuclearLevelManager::_levelEnergy =0.

staticprivate

Definition at line 121 of file G4NuclearLevelManager.hh.

_levels

G4PtrLevelVector* G4NuclearLevelManager::_levels

private

Definition at line 117 of file G4NuclearLevelManager.hh.

_m1CC

G4double G4NuclearLevelManager::_m1CC =0.

staticprivate

Definition at line 131 of file G4NuclearLevelManager.hh.

_m2CC

G4double G4NuclearLevelManager::_m2CC =0.

staticprivate

Definition at line 132 of file G4NuclearLevelManager.hh.

_m3CC

G4double G4NuclearLevelManager::_m3CC =0.

staticprivate

Definition at line 133 of file G4NuclearLevelManager.hh.

_m4CC

G4double G4NuclearLevelManager::_m4CC =0.

staticprivate

Definition at line 134 of file G4NuclearLevelManager.hh.

_m5CC

G4double G4NuclearLevelManager::_m5CC =0.

staticprivate

Definition at line 135 of file G4NuclearLevelManager.hh.

_nPlusCC

G4double G4NuclearLevelManager::_nPlusCC =0.

staticprivate

Definition at line 136 of file G4NuclearLevelManager.hh.

_nucleusA

G4int G4NuclearLevelManager::_nucleusA

private

Definition at line 113 of file G4NuclearLevelManager.hh.

_nucleusZ

G4int G4NuclearLevelManager::_nucleusZ

private

Definition at line 114 of file G4NuclearLevelManager.hh.

_polarity

G4double G4NuclearLevelManager::_polarity =0.

staticprivate

Definition at line 124 of file G4NuclearLevelManager.hh.

_probability

G4double G4NuclearLevelManager::_probability =0.

staticprivate

Definition at line 123 of file G4NuclearLevelManager.hh.

_totalCC

G4double G4NuclearLevelManager::_totalCC =0.

staticprivate

Definition at line 137 of file G4NuclearLevelManager.hh.

_validity

G4bool G4NuclearLevelManager::_validity

private

Definition at line 116 of file G4NuclearLevelManager.hh.

---

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

• G4NuclearLevelManager.hh
• G4NuclearLevelManager.cc