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)
void	AddState (G4int, G4int, G4double, G4int, G4double, G4int ionJ=0, G4double ionMu=0.0)
void	AddState (G4int, G4int, G4double, G4Ions::G4FloatLevelBase, G4double, G4int ionJ=0, G4double ionMu=0.0)
size_t	GetSizeOfIsotopeList ()
virtual G4IsotopeProperty *	GetIsotope (G4int Z, G4int A, G4double E, G4Ions::G4FloatLevelBase flb=G4Ions::G4FloatLevelBase::no_Float)
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 ()

Detailed Description

Definition at line 58 of file G4NuclideTable.hh.

Member Typedef Documentation

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

Definition at line 73 of file G4NuclideTable.hh.

Constructor & Destructor Documentation

G4NuclideTable::~G4NuclideTable ( )

virtual

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

}

Member Function Documentation

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

Definition at line 311 of file G4NuclideTable.cc.

{
   if ( G4Threading::IsMasterThread() ) {
      G4int flbIndex = 0;
      ionE = StripFloatLevelBase( ionE, flbIndex );
      AddState(ionZ,ionA,ionE,flbIndex,ionLife,ionJ,ionMu);
   }
}

Here is the call graph for this function:

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

Definition at line 320 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);
   fProperty->SetFloatLevelBase(flbIndex);

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

   }
}

Here is the call graph for this function:

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

Definition at line 346 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);
      fProperty->SetFloatLevelBase( flb );

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

   }
}

Here is the call graph for this function:

size_t G4NuclideTable::entries ( ) const

inline

Definition at line 173 of file G4NuclideTable.hh.

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

Here is the caller graph for this function:

void G4NuclideTable::FillHardCodeList ( )

protected

Definition at line 170 of file G4NuclideTable.cc.

{
   ;
}

void G4NuclideTable::GenerateNuclide

(

)

Definition at line 176 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;
      G4String ionFL;
      G4double ionLife;
      G4int ionJ;
      G4double ionMu;

      //ifs >> ionZ >> ionA >> ionE >> ionLife >> ionJ >> ionMu;
      ifs >> ionZ >> ionA >> ionE >> ionFL >> 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;
         //G4int flbIndex = 0;
         //ionE = StripFloatLevelBase( ionE, flbIndex );
         G4Ions::G4FloatLevelBase flb = StripFloatLevelBase( ionFL );
         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 && flb == G4Ions::G4FloatLevelBase::no_Float ) //
           || ( threshold_of_half_life <= ionLife*std::log(2.0) && ionLife*std::log(2.0) < minimum_threshold_of_half_life ) ) {

            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);
            fProperty->SetFloatLevelBase( flb );

            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 >> ionFL >> 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:

G4NuclideTable * G4NuclideTable::GetInstance ( )

static

Definition at line 56 of file G4NuclideTable.cc.

                                            {
   static G4NuclideTable instance;
   return &instance;
}

Here is the caller graph for this function:

G4IsotopeProperty * G4NuclideTable::GetIsotope ( G4int  Z, G4int  A, G4double  E, G4Ions::G4FloatLevelBase  flb = G4Ions::G4FloatLevelBase::no_Float  )

virtual

Implements G4VIsotopeTable.

Definition at line 108 of file G4NuclideTable.cc.

{

   G4IsotopeProperty* fProperty = nullptr;

   // At first searching UserDefined
   if ( fUserDefinedList ) {
      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 ) {
               if( flb == (*it)->GetFloatLevelBase() )
               { 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 ) {
            if( flb == (lower_bound_itr->second)->GetFloatLevelBase() )
            { return lower_bound_itr->second; } // found 
         }
      }
*/
      while ( lower_bound_itr != itf -> second.end() ) {
         levelE = lower_bound_itr->first;
         if ( levelE - flevelTolerance/2 <= E && E < levelE + flevelTolerance/2 ) {
            if ( flb == (lower_bound_itr->second)->GetFloatLevelBase() ) return lower_bound_itr->second; // found 
         } else {
            break;
         } 
         lower_bound_itr++;    
      }
        
   }

   return fProperty; // not found
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4IsotopeProperty * G4NuclideTable::GetIsotopeByIndex ( size_t  idx ) const

inline

Definition at line 179 of file G4NuclideTable.hh.

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

Here is the caller graph for this function:

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

virtual

Reimplemented from G4VIsotopeTable.

Definition at line 163 of file G4NuclideTable.cc.

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

Here is the call graph for this function:

G4double G4NuclideTable::GetLevelTolerance ( )

inline

Definition at line 93 of file G4NuclideTable.hh.

{ return flevelTolerance; };

Here is the caller graph for this function:

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:

size_t G4NuclideTable::GetSizeOfIsotopeList ( )

inline

Definition at line 104 of file G4NuclideTable.hh.

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

G4double G4NuclideTable::GetThresholdOfHalfLife ( )

inline

Definition at line 90 of file G4NuclideTable.hh.

{ return threshold_of_half_life; };

static G4double G4NuclideTable::GetTrancationError ( G4double  eex )

inlinestatic

Definition at line 129 of file G4NuclideTable.hh.

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

Here is the call graph for this function:

static G4double G4NuclideTable::Round ( G4double  eex )

inlinestatic

Definition at line 132 of file G4NuclideTable.hh.

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

Here is the call graph for this function:

void G4NuclideTable::SetLevelTolerance ( G4double  x )

inline

Definition at line 92 of file G4NuclideTable.hh.

{ flevelTolerance=x;};

Here is the caller graph for this function:

void G4NuclideTable::SetThresholdOfHalfLife

(

G4double

t

)

Definition at line 373 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:

static G4double G4NuclideTable::Tolerance ( )

inlinestatic

Definition at line 138 of file G4NuclideTable.hh.

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

Here is the call graph for this function:

static G4long G4NuclideTable::Trancate ( G4double  eex )

inlinestatic

Definition at line 135 of file G4NuclideTable.hh.

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

Here is the call graph for this function:

---

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

• geant4.10.03.p01/source/particles/management/include/G4NuclideTable.hh
• geant4.10.03.p01/source/particles/management/src/G4NuclideTable.cc