G4NuclideTable Class Reference

#include <G4NuclideTable.hh>

Inheritance diagram for G4NuclideTable:

Collaboration diagram for G4NuclideTable:

Public Types
typedef std::vector< G4IsotopeProperty * >	G4IsotopeList

Public Member Functions
virtual	~G4NuclideTable ()
void	GenerateNuclide ()
void	SetThresholdOfHalfLife (G4double)
G4double	GetThresholdOfHalfLife ()
void	SetLevelTolerance (G4double x)
G4double	GetLevelTolerance ()
void	AddState (G4int, G4int, G4double, G4double, G4int ionJ=0, G4double ionMu=0.0)
size_t	GetSizeOfIsotopeList ()
virtual G4IsotopeProperty *	GetIsotope (G4int Z, G4int A, G4double E)
virtual G4IsotopeProperty *	GetIsotopeByIsoLvl (G4int Z, G4int A, G4int lvl=0)
size_t	entries () const
G4IsotopeProperty *	GetIsotopeByIndex (size_t idx) const
Public Member Functions inherited from G4VIsotopeTable
	G4VIsotopeTable ()
	G4VIsotopeTable (const G4String &)
	G4VIsotopeTable (const G4VIsotopeTable &)
G4VIsotopeTable &	operator= (const G4VIsotopeTable &)
virtual	~G4VIsotopeTable ()
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int level)
void	DumpTable (G4int Zmin=1, G4int Zmax=118)
const G4String &	GetName () const

Static Public Member Functions
static G4NuclideTable *	GetInstance ()
static G4NuclideTable *	GetNuclideTable ()
static G4double	GetTrancationError (G4double eex)
static G4double	Round (G4double eex)
static G4long	Trancate (G4double eex)
static G4double	Tolerance ()

Protected Member Functions
void	FillHardCodeList ()

Private Types
enum	{   idxZ =0, idxA, idxEnergy, idxLife,   idxSpin, idxMu }

Private Member Functions
	G4NuclideTable ()

Private Attributes
G4double	threshold_of_half_life
G4double	minimum_threshold_of_half_life
G4IsotopeList *	fUserDefinedList
std::map< G4int, std::multimap< G4double, G4IsotopeProperty *> >	map_pre_load_list
std::map< G4int, std::multimap< G4double, G4IsotopeProperty *> >	map_full_list
G4IsotopeList *	fIsotopeList
G4double	flevelTolerance
G4NuclideTableMessenger *	fMessenger

Detailed Description

Definition at line 58 of file G4NuclideTable.hh.

Member Typedef Documentation

G4IsotopeList

typedef std::vector<G4IsotopeProperty*> G4NuclideTable::G4IsotopeList

Definition at line 73 of file G4NuclideTable.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
idxZ
idxA
idxEnergy
idxLife
idxSpin
idxMu

Definition at line 157 of file G4NuclideTable.hh.

{idxZ=0, idxA,idxEnergy, idxLife, idxSpin, idxMu };

Constructor & Destructor Documentation

G4NuclideTable()

G4NuclideTable::G4NuclideTable ( )

private

Definition at line 61 of file G4NuclideTable.cc.

  :G4VIsotopeTable("Isomer"),
   threshold_of_half_life(1000.0*ns),
   minimum_threshold_of_half_life(DBL_MAX),
   fUserDefinedList(NULL), 
   fIsotopeList(NULL),
   flevelTolerance(1.0*eV)
{
  //SetVerboseLevel(G4ParticleTable::GetParticleTable()->GetVerboseLevel());
  //FillHardCodeList();
  fMessenger = new G4NuclideTableMessenger(this);
  fIsotopeList = new G4IsotopeList();
  GenerateNuclide();
}

Here is the call graph for this function:

~G4NuclideTable()

G4NuclideTable::~G4NuclideTable ( )

virtual

Definition at line 77 of file G4NuclideTable.cc.

{

  for ( std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator 
     it = map_pre_load_list.begin(); it != map_pre_load_list.end(); it++ ) {
     it->second.clear();
  }
  map_pre_load_list.clear();


  for ( std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator 
     it = map_full_list.begin(); it != map_full_list.end(); it++ ) {
     it->second.clear();
  }
  map_full_list.clear();

  if (fIsotopeList!=0) {
    for (size_t i = 0 ; i<fIsotopeList->size(); i++) {
       //G4IsotopeProperty* fProperty = (*fIsotopeList)[i]; std::cout << fProperty->GetAtomicNumber() << " " << fProperty->GetAtomicMass() << " " << fProperty->GetEnergy() << std::endl;
      delete (*fIsotopeList)[i];
    }
    fIsotopeList->clear();
    delete fIsotopeList;
    fIsotopeList = 0;
  }

}

Here is the caller graph for this function:

Member Function Documentation

AddState()

void G4NuclideTable::AddState	(	G4int	ionZ,
		G4int	ionA,
		G4double	ionE,
		G4double	ionLife,
		G4int	ionJ = 0,
		G4double	ionMu = 0.0
	)

Definition at line 287 of file G4NuclideTable.cc.

{
   if ( G4Threading::IsMasterThread() ) {

   if ( fUserDefinedList == NULL ) fUserDefinedList = new G4IsotopeList();

   G4IsotopeProperty* fProperty = new G4IsotopeProperty(); 

   // Set Isotope Property
   fProperty->SetAtomicNumber(ionZ);
   fProperty->SetAtomicMass(ionA);
   fProperty->SetIsomerLevel(9);
   fProperty->SetEnergy(ionE);
   fProperty->SetiSpin(ionJ);
   fProperty->SetLifeTime(ionLife);
   fProperty->SetDecayTable(0);
   fProperty->SetMagneticMoment(ionMu);

   fUserDefinedList->push_back(fProperty);
   fIsotopeList->push_back(fProperty);

   }
}

Here is the call graph for this function:

Here is the caller graph for this function:

entries()

size_t G4NuclideTable::entries ( ) const

inline

Definition at line 166 of file G4NuclideTable.hh.

{
  return (fIsotopeList ? fIsotopeList->size() : static_cast<size_t>(0) );
}

Here is the caller graph for this function:

FillHardCodeList()

void G4NuclideTable::FillHardCodeList ( )

protected

Definition at line 155 of file G4NuclideTable.cc.

{
   ;
}

Here is the caller graph for this function:

GenerateNuclide()

void G4NuclideTable::GenerateNuclide

(

)

Definition at line 161 of file G4NuclideTable.cc.

{

   if ( threshold_of_half_life < minimum_threshold_of_half_life ) {

      // Need to update full list

      char* path = getenv("G4ENSDFSTATEDATA");

      if ( !path ) {
         G4Exception("G4NuclideTable", "PART70000",
                     FatalException, "G4ENSDFSTATEDATA environment variable must be set");
         return;
      }

      std::ifstream ifs;
      G4String filename(path);
      filename += "/ENSDFSTATE.dat";

      ifs.open( filename.c_str() );

      if ( !ifs.good() ) {
         G4Exception("G4NuclideTable", "PART70001",
                     FatalException, "ENSDFSTATE.dat is not found.");
         return;
      }


      G4int ionCode=0;
      G4int iLevel=0;

      G4int ionZ;
      G4int ionA;
      G4double ionE;
      G4double ionLife;
      G4int ionJ;
      G4double ionMu;

      ifs >> ionZ >> ionA >> ionE >> ionLife >> ionJ >> ionMu;

      while ( ifs.good() ) {// Loop checking, 09.08.2015, K.Kurashige

         if ( ionCode != 1000*ionZ + ionA ) {
              iLevel = 0;
              ionCode = 1000*ionZ + ionA;
         }

         ionE *= keV;
         ionLife *= ns;
         ionMu *= (joule/tesla);

         if ( ( ionE == 0 && minimum_threshold_of_half_life == DBL_MAX ) // ground state is alwyas build in very first attempt
           || ( threshold_of_half_life <= ionLife*std::log(2.0) && ionLife*std::log(2.0) < minimum_threshold_of_half_life && ionE != 0 ) ) {

            if ( ionE > 0 ) iLevel++;
            if ( iLevel > 9 ) iLevel=9;

            G4IsotopeProperty* fProperty = new G4IsotopeProperty();

            // Set Isotope Property
            fProperty->SetAtomicNumber(ionZ);
            fProperty->SetAtomicMass(ionA);
            fProperty->SetIsomerLevel(iLevel);
            fProperty->SetEnergy(ionE);
            fProperty->SetiSpin(ionJ);
            fProperty->SetLifeTime(ionLife);
            fProperty->SetDecayTable(0);
            fProperty->SetMagneticMoment(ionMu);

            fIsotopeList->push_back(fProperty);

            std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator itf = map_full_list.find( ionCode );
            if ( itf == map_full_list.end() ) {
               std::multimap<G4double, G4IsotopeProperty*> aMultiMap;
               //itf = map_full_list.insert( std::pair< G4int , std::multimap< G4double , G4IsotopeProperty* > > ( ionCode , aMultiMap ) );
               itf = ( map_full_list.insert( std::pair< G4int , std::multimap< G4double , G4IsotopeProperty* > > ( ionCode , aMultiMap ) ) ).first;
            } 
            itf -> second.insert( std::pair< G4double, G4IsotopeProperty* >( ionE , fProperty ) );
         }

         ifs >> ionZ >> ionA >> ionE >> ionLife >> ionJ >> ionMu;
      }

      minimum_threshold_of_half_life = threshold_of_half_life;

   }


   // Clear current map
   for ( std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator
      it = map_pre_load_list.begin(); it != map_pre_load_list.end(); it++ ) {
      it->second.clear();
   }
   map_pre_load_list.clear();

   // Build map based on current threshold value 
   for ( std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator
      it = map_full_list.begin(); it != map_full_list.end(); it++ ) {

      G4int ionCode = it->first;
      std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator itf = map_pre_load_list.find( ionCode );
      if ( itf == map_pre_load_list.end() ) {
         std::multimap<G4double, G4IsotopeProperty*> aMultiMap;
         itf = ( map_pre_load_list.insert( std::pair< G4int , std::multimap< G4double , G4IsotopeProperty* > > ( ionCode , aMultiMap ) ) ).first;
      }
      G4int iLevel = 0;
      for ( std::multimap< G4double , G4IsotopeProperty* >::iterator
         itt = it->second.begin(); itt != it->second.end(); itt++ ) {

         G4double exEnergy = itt->first;
         G4double meanLife = itt->second->GetLifeTime();

         if ( exEnergy == 0.0
           || meanLife*std::log(2.0) > threshold_of_half_life ) {

            if ( itt->first != 0.0 ) iLevel++;
            if ( iLevel > 9 ) iLevel=9;
            itt->second->SetIsomerLevel( iLevel );

            itf -> second.insert( std::pair< G4double, G4IsotopeProperty* >( exEnergy , itt->second ) );
         }
      }
   }

}

Here is the call graph for this function:

Here is the caller graph for this function:

GetInstance()

G4NuclideTable * G4NuclideTable::GetInstance ( void  )

static

Definition at line 55 of file G4NuclideTable.cc.

                                            {
   static G4NuclideTable instance;
   return &instance;
}

Here is the caller graph for this function:

GetIsotope()

G4IsotopeProperty * G4NuclideTable::GetIsotope ( G4int  Z, G4int  A, G4double  E  )

virtual

Implements G4VIsotopeTable.

Definition at line 107 of file G4NuclideTable.cc.

{

   G4IsotopeProperty* fProperty = NULL;

   // At first searching UserDefined
   if ( fUserDefinedList != NULL ) {
      for ( G4IsotopeList::iterator it = fUserDefinedList->begin() ; it != fUserDefinedList->end() ; it++ ) {
        
         if ( Z == (*it)->GetAtomicNumber() && A == (*it)->GetAtomicMass() ) {
            G4double levelE = (*it)->GetEnergy();         
            if ( levelE - flevelTolerance/2 <= E && E < levelE + flevelTolerance/2 ) {
               return *it; //found
            }
         }

      }
   } 

   //Serching pre-load
   //Note: isomer level is properly set only for pre_load_list.
   //
   G4int ionCode = 1000*Z + A;
   std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > >::iterator itf = map_pre_load_list.find( ionCode );

   if ( itf !=  map_pre_load_list.end() ) {
      std::multimap< G4double , G4IsotopeProperty* >::iterator lower_bound_itr = itf -> second.lower_bound ( E - flevelTolerance/2 );
      G4double levelE = DBL_MAX;
      if ( lower_bound_itr !=  itf -> second.end() ) {
         levelE = lower_bound_itr->first;
         if ( levelE - flevelTolerance/2 <= E && E < levelE + flevelTolerance/2 ) {
            return lower_bound_itr->second; // found
         }
      }
   }

   return fProperty; // not found;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetIsotopeByIndex()

G4IsotopeProperty * G4NuclideTable::GetIsotopeByIndex ( size_t  idx ) const

inline

Definition at line 172 of file G4NuclideTable.hh.

{
  if ( fIsotopeList && idx<fIsotopeList->size()) return (*fIsotopeList)[idx];
  else                          return 0;
}

Here is the caller graph for this function:

GetIsotopeByIsoLvl()

G4IsotopeProperty * G4NuclideTable::GetIsotopeByIsoLvl ( G4int  Z, G4int  A, G4int  lvl = 0  )

virtual

Reimplemented from G4VIsotopeTable.

Definition at line 148 of file G4NuclideTable.cc.

{
  if(lvl==0) return GetIsotope(Z,A,0.0);
  return (G4IsotopeProperty*)0;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetLevelTolerance()

G4double G4NuclideTable::GetLevelTolerance ( )

inline

Definition at line 93 of file G4NuclideTable.hh.

{ return flevelTolerance; };

Here is the call graph for this function:

Here is the caller graph for this function:

GetNuclideTable()

static G4NuclideTable* G4NuclideTable::GetNuclideTable ( )

inlinestatic

Definition at line 73 of file G4NuclideTable.hh.

{ return GetInstance(); };

Here is the call graph for this function:

Here is the caller graph for this function:

GetSizeOfIsotopeList()

size_t G4NuclideTable::GetSizeOfIsotopeList ( )

inline

Definition at line 103 of file G4NuclideTable.hh.

{ return ( fIsotopeList ? fIsotopeList->size() : static_cast<size_t>(0)) ; };

Here is the call graph for this function:

GetThresholdOfHalfLife()

G4double G4NuclideTable::GetThresholdOfHalfLife ( )

inline

Definition at line 90 of file G4NuclideTable.hh.

{ return threshold_of_half_life; };

GetTrancationError()

static G4double G4NuclideTable::GetTrancationError ( G4double  eex )

inlinestatic

Definition at line 125 of file G4NuclideTable.hh.

  { G4double tolerance= G4NuclideTable::GetInstance()->GetLevelTolerance();
  return eex - (G4long)(eex/tolerance)*tolerance; }

Here is the call graph for this function:

Round()

static G4double G4NuclideTable::Round ( G4double  eex )

inlinestatic

Definition at line 128 of file G4NuclideTable.hh.

  { G4double tolerance= G4NuclideTable::GetInstance()->GetLevelTolerance();
    return round(eex/tolerance)*tolerance; }

Here is the call graph for this function:

SetLevelTolerance()

void G4NuclideTable::SetLevelTolerance ( G4double  x )

inline

Definition at line 92 of file G4NuclideTable.hh.

{ flevelTolerance=x;};

Here is the caller graph for this function:

SetThresholdOfHalfLife()

void G4NuclideTable::SetThresholdOfHalfLife

(

G4double

t

)

Definition at line 312 of file G4NuclideTable.cc.

{
   if ( G4Threading::IsMasterThread() ) {
      threshold_of_half_life=t; 
      GenerateNuclide();
   }
}

Here is the call graph for this function:

Here is the caller graph for this function:

Tolerance()

static G4double G4NuclideTable::Tolerance ( )

inlinestatic

Definition at line 134 of file G4NuclideTable.hh.

  { return G4NuclideTable::GetInstance()->GetLevelTolerance(); }

Here is the call graph for this function:

Trancate()

static G4long G4NuclideTable::Trancate ( G4double  eex )

inlinestatic

Definition at line 131 of file G4NuclideTable.hh.

  { G4double tolerance= G4NuclideTable::GetInstance()->GetLevelTolerance();
    return (G4long)(eex/tolerance); }

Here is the call graph for this function:

Member Data Documentation

fIsotopeList

G4IsotopeList* G4NuclideTable::fIsotopeList

private

Definition at line 159 of file G4NuclideTable.hh.

flevelTolerance

G4double G4NuclideTable::flevelTolerance

private

Definition at line 160 of file G4NuclideTable.hh.

fMessenger

G4NuclideTableMessenger* G4NuclideTable::fMessenger

private

Definition at line 161 of file G4NuclideTable.hh.

fUserDefinedList

G4IsotopeList* G4NuclideTable::fUserDefinedList

private

Definition at line 141 of file G4NuclideTable.hh.

map_full_list

std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > > G4NuclideTable::map_full_list

private

Definition at line 148 of file G4NuclideTable.hh.

map_pre_load_list

std::map< G4int , std::multimap< G4double , G4IsotopeProperty* > > G4NuclideTable::map_pre_load_list

private

Definition at line 147 of file G4NuclideTable.hh.

minimum_threshold_of_half_life

G4double G4NuclideTable::minimum_threshold_of_half_life

private

Definition at line 140 of file G4NuclideTable.hh.

threshold_of_half_life

G4double G4NuclideTable::threshold_of_half_life

private

Definition at line 139 of file G4NuclideTable.hh.

---

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

• G4NuclideTable.hh
• G4NuclideTable.cc