G4RIsotopeTable Class Reference

#include <G4RIsotopeTable.hh>

Inheritance diagram for G4RIsotopeTable:

Collaboration diagram for G4RIsotopeTable:

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

Public Member Functions
	G4RIsotopeTable ()
virtual	~G4RIsotopeTable ()
virtual G4bool	FindIsotope (G4IsotopeProperty *property)
virtual G4IsotopeProperty *	GetIsotope (G4int Z, G4int A, G4double E)
void	AddUserDecayDataFile (G4int Z, G4int A, G4String filename)
Public Member Functions inherited from G4VIsotopeTable
	G4VIsotopeTable ()
	G4VIsotopeTable (const G4String &)
	G4VIsotopeTable (const G4VIsotopeTable &)
G4VIsotopeTable &	operator= (const G4VIsotopeTable &)
virtual	~G4VIsotopeTable ()
virtual G4IsotopeProperty *	GetIsotopeByIsoLvl (G4int Z, G4int A, G4int level=0)
G4int	GetVerboseLevel () const
void	SetVerboseLevel (G4int level)
void	DumpTable (G4int Zmin=1, G4int Zmax=118)
const G4String &	GetName () const

Protected Member Functions
	G4RIsotopeTable (const G4RIsotopeTable &right)
G4RIsotopeTable &	operator= (const G4RIsotopeTable &right)

Private Member Functions
G4String	GetIsotopeName (G4int Z, G4int A, G4double E)
G4double	GetMeanLifeTime (G4int Z, G4int A, G4double &E)
G4IsotopeProperty *	GetIsotope (G4int index) const
G4int	GetVerboseLevel () const
G4int	Entries () const

Private Attributes
G4IsotopeList	fIsotopeList
G4IsotopeNameList	fIsotopeNameList
std::map< G4int, G4String >	theUserRadioactiveDataFiles

Static Private Attributes
static const G4double	levelTolerance = 2.0*keV

Detailed Description

Definition at line 64 of file G4RIsotopeTable.hh.

Member Typedef Documentation

G4IsotopeList

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

Definition at line 72 of file G4RIsotopeTable.hh.

G4IsotopeNameList

typedef std::vector<G4String> G4RIsotopeTable::G4IsotopeNameList

Definition at line 73 of file G4RIsotopeTable.hh.

Constructor & Destructor Documentation

G4RIsotopeTable() [1/2]

G4RIsotopeTable::G4RIsotopeTable

(

)

Definition at line 84 of file G4RIsotopeTable.cc.

                                :G4VIsotopeTable("RIsotopeTable")
{
  // Reset the list of user defined data files
  theUserRadioactiveDataFiles.clear();
}

G4RIsotopeTable() [2/2]

G4RIsotopeTable::G4RIsotopeTable ( const G4RIsotopeTable &  right )

protected

~G4RIsotopeTable()

G4RIsotopeTable::~G4RIsotopeTable ( )

virtual

Definition at line 91 of file G4RIsotopeTable.cc.

{
  for (G4int i = 0; i < G4int(fIsotopeList.size()); i++) delete fIsotopeList[i];
  fIsotopeList.clear();
  fIsotopeNameList.clear(); 
}

Member Function Documentation

AddUserDecayDataFile()

void G4RIsotopeTable::AddUserDecayDataFile	(	G4int	Z,
		G4int	A,
		G4String	filename
	)

Definition at line 281 of file G4RIsotopeTable.cc.

{ if (Z<1 || A<2) {
    G4cout<<"Z and A not valid!"<<G4endl;
  }

  std::ifstream DecaySchemeFile(filename);
  if (DecaySchemeFile){
    G4int ID_ion=A*1000+Z;
    theUserRadioactiveDataFiles[ID_ion]=filename;
  }
  else {
    G4cout<<"The file "<<filename<<" does not exist!"<<G4endl;
  }
}

Entries()

G4int G4RIsotopeTable::Entries ( ) const

inlineprivate

Definition at line 137 of file G4RIsotopeTable.hh.

  {return fIsotopeNameList.size();}

Here is the caller graph for this function:

FindIsotope()

G4bool G4RIsotopeTable::FindIsotope ( G4IsotopeProperty *  property )

virtual

Definition at line 105 of file G4RIsotopeTable.cc.

{
  // do nothing, it is here just for the compiler
  // it is required by the base class
  return true;
}

GetIsotope() [1/2]

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

virtual

Implements G4VIsotopeTable.

Definition at line 113 of file G4RIsotopeTable.cc.

{
  G4String fname = GetIsotopeName(Z, A, E);  
  G4int j = -1;
  for (G4int i = 0 ; i< Entries(); i++) {
    if(fIsotopeNameList[i] == fname) j = i;}
  if (j >=0) {
    if (GetVerboseLevel() > 1) {
      G4cout <<"G4RIsotopeTable::GetIsotope No. : ";
      G4cout <<j<<G4endl;   
    }
    return  GetIsotope(j);
    // isotope property data has been loaded already - just return the pointer

  } else {
    G4double meanlife = GetMeanLifeTime(Z, A, E);
    // E is passed as a refence hence on entry E is supplied by the user and it 
    // could be slightly different from the returned value which is the one 
    // defined in the database.
    // this call is to ensure the code uses a consistent E value throughout.
    
    G4IsotopeProperty* fProperty = new G4IsotopeProperty();   
    // Set Isotope Property
    fProperty->SetLifeTime(meanlife);
    fProperty->SetAtomicNumber(Z);
    fProperty->SetAtomicMass(A);
    // Notic that the value of E may have been changed
    fProperty->SetEnergy(E);
    // The spin is not being used in the current implementation
    fProperty->SetiSpin(0);
    // the decaytable will be loaded later in G4RadioactiveDecay when it is needed
    fProperty->SetDecayTable(0);
    
    fIsotopeList.push_back(fProperty);
    fname = GetIsotopeName(Z, A, E);
    fIsotopeNameList.push_back(fname);
    if (GetVerboseLevel()>1) {
      G4cout <<"G4RIsotopeTable::GetIsotope create: ";
      G4cout <<fname <<G4endl;  
    }
    return fProperty;

  }
}

Here is the call graph for this function:

GetIsotope() [2/2]

G4IsotopeProperty * G4RIsotopeTable::GetIsotope ( G4int  index ) const

inlineprivate

Definition at line 140 of file G4RIsotopeTable.hh.

{
  if ( (index >=0 ) && (index < Entries()) ){
    return fIsotopeList[index];
  } else {
    return 0; 
  } 
}

Here is the call graph for this function:

GetIsotopeName()

G4String G4RIsotopeTable::GetIsotopeName ( G4int  Z, G4int  A, G4double  E  )

private

Definition at line 159 of file G4RIsotopeTable.cc.

{
  std::ostringstream os;
  os.setf(std::ios::fixed);
  os <<"A"<< A << "Z" << Z <<'[' << std::setprecision(1) << E/keV << ']';
  G4String name = os.str();
  if (GetVerboseLevel()>1) {
    G4cout <<"G4RIsotopeTable::GetIsotope Name: ";
    G4cout <<name <<G4endl;   
  }
  return name;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetMeanLifeTime()

G4double G4RIsotopeTable::GetMeanLifeTime ( G4int  Z, G4int  A, G4double &  E  )

private

Definition at line 173 of file G4RIsotopeTable.cc.

{

  G4double lifetime = -1.0;

  //Check if data have been provided by the user
  std::map<int,G4String>::iterator it = theUserRadioactiveDataFiles.find(1000*A+Z);
  G4String file="";
  if ( it != theUserRadioactiveDataFiles.end() ){
      file=it->second;
  }
  if (file ==""){
    if (!getenv("G4RADIOACTIVEDATA")) {
        G4cout << "Please setenv G4RADIOACTIVEDATA to point to the radioactive decay data files." << G4endl;
        throw G4HadronicException(__FILE__, __LINE__,
                  "Please setenv G4RADIOACTIVEDATA to point to the radioactive decay data files.");
    }
    G4String dirName = getenv("G4RADIOACTIVEDATA");

    std::ostringstream os;
    os <<dirName <<"/z" <<Z <<".a" <<A ;
    file = os.str();
  }
  std::ifstream DecaySchemeFile(file);

  G4bool found_in_raddecay_data(false);
  if (!DecaySchemeFile) {
    if (GetVerboseLevel()>1) {
      G4cout <<"G4RIsotopeTable::GetMeanLife() : "
         <<"cannot find ion radioactive decay file: " 
         <<file <<G4endl;
    }
  } else {
    char inputChars[100]={' '};
    G4String inputLine;
    G4String recordType("");
    G4double a(0.0);
    G4double b(0.0);

    G4int loop = 0;
    G4ExceptionDescription ed;
    ed << " While count exceeded " << G4endl;
    while (!found_in_raddecay_data && !DecaySchemeFile.getline(inputChars, 100).eof()) { /* Loop checking, 01.09.2015, D.Wright */
      loop++;
      if (loop > 100000) {
        G4Exception("G4RIsotopeTable::GetMeanLifeTime()", "HAD_RDM_100", JustWarning, ed);
        break;
      }
 
      inputLine = inputChars;
      inputLine = inputLine.strip(1);

      if (inputChars[0] != '#' && inputLine.length() != 0) {
        std::istringstream tmpstream(inputLine);
        tmpstream >> recordType >> a >> b;
        if (recordType == "P") {
          if (std::abs(a*keV-aE) < levelTolerance) {
            found_in_raddecay_data    = true;
            lifetime = b/0.693147*s ;
          }
        }
      }
    }
    DecaySchemeFile.close();
  }

    if (!found_in_raddecay_data && aE) {
      G4double half_life=-1.;
      lifetime = 1.0E-20*s;


      //added by L.Desorgher If the life time is not found in  raddecay database
      // then it is deduced from photo-evaporation level
      const G4NuclearLevel* aLevel =
            G4NuclearLevelStore::GetInstance()->GetManager(Z, A)
                                        ->NearestLevel(aE,levelTolerance);
      if (aLevel) {
          half_life = aLevel->HalfLife();
          lifetime = half_life/0.693147;
      }

      if (GetVerboseLevel()>1 && half_life<0) {
        G4cout << "G4RIsotopeTable::GetMeanLife() : ";
        G4cout << "cannot find ion of required excitation E = " << aE << G4endl;
        G4cout << "state in radioactive or photoevaporation data file " << G4endl;
        G4cout <<"The nucleus is assumed to be IT decayed with life = 1E-20 s" << G4endl;
        G4cout <<" -----------* THIS MAY CAUSE PROBLEM IN ITS DECAY-----------" << G4endl;
      }
    }

    if (!found_in_raddecay_data && !aE) {
      if (GetVerboseLevel()>1) {
        G4cout <<"G4RIsotopeTable::GetMeanLife() : ";
        G4cout <<"cannot find ion of required excitation E = " << aE << G4endl;
        G4cout <<"state in radioactive or photoevaporation data file" <<G4endl;
        G4cout <<"The nucleus is assumed to be stable" <<G4endl;
        lifetime = -1.0;
      }
    }

    if (GetVerboseLevel()>1) {
       G4cout <<"G4RIsotopeTable::GetMeanLifeTime: ";
       G4cout <<lifetime << " for " << GetIsotopeName(Z, A, aE) <<G4endl;
    }
  return lifetime;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetVerboseLevel()

G4int G4RIsotopeTable::GetVerboseLevel ( ) const

private

Definition at line 99 of file G4RIsotopeTable.cc.

{
  return G4ParticleTable::GetParticleTable()->GetVerboseLevel();
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

protected

Member Data Documentation

fIsotopeList

G4IsotopeList G4RIsotopeTable::fIsotopeList

private

Definition at line 128 of file G4RIsotopeTable.hh.

fIsotopeNameList

G4IsotopeNameList G4RIsotopeTable::fIsotopeNameList

private

Definition at line 129 of file G4RIsotopeTable.hh.

levelTolerance

const G4double G4RIsotopeTable::levelTolerance = 2.0*keV

staticprivate

Definition at line 130 of file G4RIsotopeTable.hh.

theUserRadioactiveDataFiles

std::map<G4int, G4String> G4RIsotopeTable::theUserRadioactiveDataFiles

private

Definition at line 133 of file G4RIsotopeTable.hh.

---

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

• G4RIsotopeTable.hh
• G4RIsotopeTable.cc