G4IsomerTable Class Reference

#include <G4IsomerTable.hh>

Inheritance diagram for G4IsomerTable:

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

Public Member Functions
	G4IsomerTable ()
virtual	~G4IsomerTable ()
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

Protected Member Functions
	G4IsomerTable (const G4IsomerTable &right)
G4IsomerTable &	operator= (const G4IsomerTable &right)
void	FillIsotopeList ()

Detailed Description

Definition at line 55 of file G4IsomerTable.hh.

Member Typedef Documentation

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

Definition at line 64 of file G4IsomerTable.hh.

Constructor & Destructor Documentation

G4IsomerTable::G4IsomerTable

(

)

Definition at line 53 of file G4IsomerTable.cc.

References FillIsotopeList(), G4ParticleTable::GetParticleTable(), G4VIsotopeTable::GetVerboseLevel(), and G4VIsotopeTable::SetVerboseLevel().

  :G4VIsotopeTable("Isomer"),
   fIsotopeList(0)
{
  SetVerboseLevel(G4ParticleTable::GetParticleTable()->GetVerboseLevel());
  FillIsotopeList();
  return;
}

G4IsomerTable::G4IsomerTable ( const G4IsomerTable &  right )

protected

Definition at line 109 of file G4IsomerTable.cc.

References FillIsotopeList().

  :G4VIsotopeTable(right),
   fIsotopeList(0)
{
  FillIsotopeList();
}

G4IsomerTable::~G4IsomerTable ( )

virtual

Definition at line 96 of file G4IsomerTable.cc.

{
  if (fIsotopeList!=0) {
    for (size_t i = 0 ; i<fIsotopeList->size(); i++) {
      delete (*fIsotopeList)[i];
    }
    fIsotopeList->clear();
    delete fIsotopeList;
    fIsotopeList = 0;
  }
}

Member Function Documentation

size_t G4IsomerTable::entries ( ) const

inline

Definition at line 118 of file G4IsomerTable.hh.

Referenced by GetIsotope(), and GetIsotopeByIsoLvl().

{
  return fIsotopeList->size();
}

void G4IsomerTable::FillIsotopeList ( )

protected

Definition at line 63 of file G4IsomerTable.cc.

References python.hepunit::joule, python.hepunit::keV, ns, and python.hepunit::tesla.

Referenced by G4IsomerTable(), and operator=().

{
  if(fIsotopeList !=0) return;
  
  fIsotopeList = new G4IsotopeList();
  fIsotopeList->reserve(nEntries);

  for (size_t i=0; i<nEntries; i++) {
    G4int    pid      = (G4int)(isomerTable[i][idxPID]);
    G4int    ionZ     = pid/10000; 
    G4int    ionA     = (pid-ionZ*10000)/10;
    G4int    lvl      = pid%10;
    G4double ionE     = isomerTable[i][idxEnergy]*keV;
    G4double ionLife  = isomerTable[i][idxLife]*ns;
    G4int    ionJ     = (G4int)(isomerTable[i][idxSpin]);
    G4double ionMu    = isomerTable[i][idxMu]*(joule/tesla);

    G4IsotopeProperty* fProperty = new G4IsotopeProperty(); 
    // Set Isotope Property
    fProperty->SetAtomicNumber(ionZ);
    fProperty->SetAtomicMass(ionA);
    fProperty->SetIsomerLevel(lvl);
    fProperty->SetEnergy(ionE);
    fProperty->SetiSpin(ionJ);
    fProperty->SetLifeTime(ionLife);
    fProperty->SetDecayTable(0);
    fProperty->SetMagneticMoment(ionMu);
    
    fIsotopeList->push_back(fProperty);
  }
}

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

virtual

Implements G4VIsotopeTable.

Definition at line 125 of file G4IsomerTable.cc.

References entries(), G4IsotopeProperty::GetAtomicMass(), and G4IsotopeProperty::GetEnergy().

{
  G4IsotopeProperty* fProperty = 0;
  if ((Z<MinZ)||(Z>MaxZ)||(A>MaxA)) return fProperty; // not found
  
  G4int low  = 0;
  G4int high = entries() -1;
  G4int ptr = (low+high)/2;
 
  // find Z
  G4int zptr   = (*fIsotopeList)[ptr]->GetAtomicNumber();
  while ( high-low > 1){
    if (zptr == Z){
      while ((zptr == Z)&&(ptr>0)) {
    ptr -= 1;
    zptr = (*fIsotopeList)[ptr]->GetAtomicNumber();
      }
      if (ptr!=0) ptr += 1;
      break;
    } else if (zptr >Z) {
      high = ptr;
      ptr = (low+high)/2;
    } else {
      low = ptr;
      ptr = (low+high)/2;
    }
    zptr   = (*fIsotopeList)[ptr]->GetAtomicNumber();
  }
  if ( Z == (*fIsotopeList)[low]->GetAtomicNumber()) ptr =low;
  else if ( Z != zptr ) return fProperty; // not found

  // find A
  G4int aptr   = (*fIsotopeList)[ptr]->GetAtomicMass();
  while ( (aptr<A) && (ptr+1<(G4int)(entries())) ) {
    ptr +=1;
    aptr   = (*fIsotopeList)[ptr]->GetAtomicMass();
    if (Z != (*fIsotopeList)[ptr]->GetAtomicNumber()) break;
  }
  if ( aptr != A)  return fProperty;  // not found

  // find E 
  G4double ptrE = (*fIsotopeList)[ptr]->GetEnergy();
  while ( (ptrE < E-levelTolerance ) && (ptr+1<(G4int)(entries())) ){
    ptr +=1;
    ptrE = (*fIsotopeList)[ptr]->GetEnergy();
    if (A != (*fIsotopeList)[ptr]->GetAtomicMass())  return fProperty; // not found
  }
  if (ptrE > E+levelTolerance) return fProperty; // not found
  
  // FOUND !!
  fProperty =  (*fIsotopeList)[ptr];

  return fProperty;
  
}

G4IsotopeProperty * G4IsomerTable::GetIsotopeByIndex ( size_t  idx ) const

inline

Definition at line 124 of file G4IsomerTable.hh.

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

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

virtual

Reimplemented from G4VIsotopeTable.

Definition at line 183 of file G4IsomerTable.cc.

References entries().

{
  G4IsotopeProperty* fProperty = 0;
  if ((Z<MinZ)||(Z>MaxZ)||(A>MaxA)) return fProperty; // not found

  G4int low  = 0;
  G4int high = entries() -1;
  G4int ptr = (low+high)/2;


  // find PID
  const G4int PID = 10000*Z +10*A + lvl; 
  G4int id   =  (fIsotopeList->at(ptr))->GetAtomicNumber() *10000
              + (fIsotopeList->at(ptr))->GetAtomicMass() *10
              + (fIsotopeList->at(ptr))->GetIsomerLevel();
  while ( high-low > 1){
    if (id == PID){
      break;
    } else if (id >PID) {
      high = ptr;
      ptr = (low+high)/2;
    } else {
      low = ptr;
      ptr = (low+high)/2;
    }
    id   =  (fIsotopeList->at(ptr))->GetAtomicNumber() *10000
          + (fIsotopeList->at(ptr))->GetAtomicMass() *10
          + (fIsotopeList->at(ptr))->GetIsomerLevel();
  }

  if (id == PID) {
    fProperty =  (*fIsotopeList)[ptr];
  } else {
    id   =  (fIsotopeList->at(low))->GetAtomicNumber() *10000
          + (fIsotopeList->at(low))->GetAtomicMass() *10
          + (fIsotopeList->at(low))->GetIsomerLevel();
    if (id == PID) {
      fProperty =  (*fIsotopeList)[low];
    } else {
      id   =  (fIsotopeList->at(high))->GetAtomicNumber() *10000
            + (fIsotopeList->at(high))->GetAtomicMass() *10
            + (fIsotopeList->at(high))->GetIsomerLevel();
      if (id == PID) fProperty =  (*fIsotopeList)[high];
    }
  }
  return fProperty;
  
}

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

protected

Definition at line 117 of file G4IsomerTable.cc.

References FillIsotopeList().

{
  FillIsotopeList();
  return *this;
}

---

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

• G4IsomerTable.hh
• G4IsomerTable.cc
• G4IsomerTable2.cc