G4IsotopeMagneticMomentTable Class Reference

#include <G4IsotopeMagneticMomentTable.hh>

Inheritance diagram for G4IsotopeMagneticMomentTable:

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

Public Member Functions
	G4IsotopeMagneticMomentTable ()
virtual	~G4IsotopeMagneticMomentTable ()
virtual G4bool	FindIsotope (G4IsotopeProperty *property)
virtual G4IsotopeProperty *	GetIsotope (G4int Z, G4int A, G4double E)
virtual G4IsotopeProperty *	GetIsotopeByIsoLvl (G4int Z, G4int A, G4int lvl=0)
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
	G4IsotopeMagneticMomentTable (const G4IsotopeMagneticMomentTable &right)
G4IsotopeMagneticMomentTable &	operator= (const G4IsotopeMagneticMomentTable &right)

Detailed Description

Definition at line 53 of file G4IsotopeMagneticMomentTable.hh.

Member Typedef Documentation

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

Definition at line 62 of file G4IsotopeMagneticMomentTable.hh.

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

Definition at line 63 of file G4IsotopeMagneticMomentTable.hh.

Constructor & Destructor Documentation

G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable

(

)

Definition at line 56 of file G4IsotopeMagneticMomentTable.cc.

References G4cout, G4endl, G4Exception(), JustWarning, python.hepunit::MeV, G4IsotopeProperty::SetAtomicMass(), G4IsotopeProperty::SetAtomicNumber(), G4IsotopeProperty::SetEnergy(), G4IsotopeProperty::SetIsomerLevel(), G4IsotopeProperty::SetiSpin(), and G4IsotopeProperty::SetMagneticMoment().

  :G4VIsotopeTable("MagneticMoment")
{
  if ( !getenv("G4IONMAGNETICMOMENT")) {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>1) {
      G4cout << "G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable():  " 
         <<  "Please setenv G4IONMAGNETICMOMENT for the magnetic moment data." 
         << G4endl;
      G4Exception( "G4IsotopeMagneticMomentTable",
           "File Not Found",
           JustWarning, 
           "Please setenv G4IONMAGNETICMOMENT");
    }
#endif
    G4Exception( "G4IsotopeMagneticMomentTable",
         "File Not Found",
         JustWarning, 
         "Please setenv G4IONMAGNETICMOMENT");
    return;
  }
  
  G4String file = getenv("G4IONMAGNETICMOMENT");
  std::ifstream DataFile(file);

  if (!DataFile ) {
#ifdef G4VERBOSE
    if (GetVerboseLevel()>0) {
      G4cout << "G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable():  " 
        << file << " is not found " << G4endl;
    }
#endif
    G4Exception( "G4IsotopeMagneticMomentTable",
         "File Not Found",
         JustWarning, 
         "Can not open G4IONMAGNETICMOMENT file");
    return;
  }
  
  char inputChars[80]={' '};
   
  while ( !DataFile.eof() ) {
    DataFile.getline(inputChars, 80);
    G4String inputLine = inputChars;
    G4int ionA, ionZ, ionJ, isomer;
    G4double ionE, ionMu, ionLife;
    G4String ionName, ionLifeUnit;
 
    if (inputChars[0] != '#' && inputLine.length() != 0) {
      std::istringstream tmpstream(inputLine);
      tmpstream >> ionZ >>  ionName >> ionA 
        >> isomer >>  ionE 
        >> ionLife >> ionLifeUnit
        >> ionJ >> ionMu;
          
      G4IsotopeProperty* fProperty = new G4IsotopeProperty();   
      // Set Isotope Property
      fProperty->SetAtomicNumber(ionZ);
      fProperty->SetAtomicMass(ionA);
      fProperty->SetIsomerLevel(isomer);
      fProperty->SetEnergy(ionE * MeV);
      fProperty->SetiSpin(ionJ);
      fProperty->SetMagneticMoment(ionMu*nuclearMagneton);
      
      fIsotopeList.push_back(fProperty);

      //if (GetVerboseLevel()>2) {
      // fProperty->DumpInfo();
      //}
    
    }
  }

  DataFile.close();
}

G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable ( const G4IsotopeMagneticMomentTable &  right )

protected

Definition at line 142 of file G4IsotopeMagneticMomentTable.cc.

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

G4IsotopeMagneticMomentTable::~G4IsotopeMagneticMomentTable ( )

virtual

Definition at line 133 of file G4IsotopeMagneticMomentTable.cc.

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

Member Function Documentation

G4bool G4IsotopeMagneticMomentTable::FindIsotope ( G4IsotopeProperty *  property )

virtual

Definition at line 155 of file G4IsotopeMagneticMomentTable.cc.

References G4IsotopeProperty::GetAtomicMass(), G4IsotopeProperty::GetAtomicNumber(), G4IsotopeProperty::GetEnergy(), and G4IsotopeProperty::GetIsomerLevel().

{
  for (size_t i = 0 ; i< fIsotopeList.size(); ++i) {
    G4IsotopeProperty*  fP = fIsotopeList[i];
    
    // check Z
    if ( fP->GetAtomicNumber() > pP->GetAtomicNumber()) {
      // Not Found
      break;
    }
    if ( fP->GetAtomicNumber() < pP->GetAtomicNumber()) {
      // next
      continue;
    }
    
    // check A
    if ( fP->GetAtomicMass() != pP->GetAtomicMass()) {
      // next
      continue;
    }
    
    //check isomerLevel
    if (fP->GetIsomerLevel() != pP->GetIsomerLevel()) {
      // next
      continue;     
    }

    //check E
    if (std::fabs(fP->GetEnergy() - pP->GetEnergy()) < levelTolerance) {
      // Found
      return true;     
    }
    
  }
  return false;
}

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

virtual

Implements G4VIsotopeTable.

Definition at line 194 of file G4IsotopeMagneticMomentTable.cc.

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

{
  G4IsotopeProperty* fProperty = 0;
  for (size_t i = 0 ; i< fIsotopeList.size(); ++i) {
    G4IsotopeProperty*  fP = fIsotopeList[i];
 
     // check Z
    if ( fP->GetAtomicNumber() > Z) {
      // Not Found
      break;
    }
    if ( fP->GetAtomicNumber() < Z) {
      // next
      continue;
    }
    
    // check A
    if ( fP->GetAtomicMass() != A ) {
      // next
      continue;
    }
    
    //check E
    if (std::fabs(fP->GetEnergy() - E) < levelTolerance) {
      // Found
      fProperty = fP;
      // fP->DumpInfo();
      break;     
    }
    
  }

  return fProperty;
  
}

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

virtual

Reimplemented from G4VIsotopeTable.

Definition at line 232 of file G4IsotopeMagneticMomentTable.cc.

References G4IsotopeProperty::GetAtomicMass(), G4IsotopeProperty::GetAtomicNumber(), and G4IsotopeProperty::GetIsomerLevel().

{
  G4IsotopeProperty* fProperty = 0;
  for (size_t i = 0 ; i< fIsotopeList.size(); ++i) {
    G4IsotopeProperty*  fP = fIsotopeList[i];
 
     // check Z
    if ( fP->GetAtomicNumber() > Z) {
      // Not Found
      break;
    }
    if ( fP->GetAtomicNumber() < Z) {
      // next
      continue;
    }
    // check A
    if ( fP->GetAtomicMass() != A ) {
      // next
      continue;
    }
    
    
    //check isomerLevel
    if (fP->GetIsomerLevel() == lvl) {
      // Found
      fProperty = fP;
      //fP->DumpInfo();
      break;     
    }
    
  }

  return fProperty;
  
}

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

protected

Definition at line 149 of file G4IsotopeMagneticMomentTable.cc.

{
  return *this;
}

---

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

• G4IsotopeMagneticMomentTable.hh
• G4IsotopeMagneticMomentTable.cc