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)
Protected Member Functions
	G4IsotopeMagneticMomentTable (const G4IsotopeMagneticMomentTable &right)
G4IsotopeMagneticMomentTable &	operator= (const G4IsotopeMagneticMomentTable &right)

---

Detailed Description

Definition at line 52 of file G4IsotopeMagneticMomentTable.hh.

---

Member Typedef Documentation

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

Definition at line 61 of file G4IsotopeMagneticMomentTable.hh.

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

Definition at line 62 of file G4IsotopeMagneticMomentTable.hh.

---

Constructor & Destructor Documentation

G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable

(

)

Definition at line 55 of file G4IsotopeMagneticMomentTable.cc.

References G4cout, G4endl, G4Exception(), JustWarning, G4IsotopeProperty::SetAtomicMass(), G4IsotopeProperty::SetAtomicNumber(), G4IsotopeProperty::SetEnergy(), 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;
    G4double ionE, ionMu, ionLife;
    G4String ionName, ionLifeUnit;
 
    if (inputChars[0] != '#' && inputLine.length() != 0) {
      std::istringstream tmpstream(inputLine);
      tmpstream >> ionZ >>  ionName >> ionA >>  ionE 
                >> ionLife >> ionLifeUnit
                >> ionJ >> ionMu;
          
      G4IsotopeProperty* fProperty = new G4IsotopeProperty();   
      // Set Isotope Property
      fProperty->SetAtomicNumber(ionZ);
      fProperty->SetAtomicMass(ionA);
      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 139 of file G4IsotopeMagneticMomentTable.cc.

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

G4IsotopeMagneticMomentTable::~G4IsotopeMagneticMomentTable

(

)

[virtual]

Definition at line 130 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]

Implements G4VIsotopeTable.

Definition at line 152 of file G4IsotopeMagneticMomentTable.cc.

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

{
  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 E
    if (std::fabs(fP->GetEnergy() - pP->GetEnergy()) <= fP->GetEnergy()*levelTolerance) {
      // Found
      return true;     
    }
    
  }
  return false;
}

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

Implements G4VIsotopeTable.

Definition at line 184 of file G4IsotopeMagneticMomentTable.cc.

References G4IsotopeProperty::DumpInfo().

{
  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) <= fP->GetEnergy()*levelTolerance) {
      // Found
      fProperty = fP;
      fP->DumpInfo();
      break;     
    }
    
  }

  return fProperty;
  
}

G4IsotopeMagneticMomentTable & G4IsotopeMagneticMomentTable::operator=

(

const G4IsotopeMagneticMomentTable &

right

)

[protected]

Definition at line 146 of file G4IsotopeMagneticMomentTable.cc.

{
  return *this;
}

---

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

• G4IsotopeMagneticMomentTable.hh
• G4IsotopeMagneticMomentTable.cc

---