G4IsotopeMagneticMomentTable Class Reference

#include <G4IsotopeMagneticMomentTable.hh>

Inheritance diagram for G4IsotopeMagneticMomentTable:

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

Public Member Functions
	G4IsotopeMagneticMomentTable ()
	~G4IsotopeMagneticMomentTable () override
virtual G4bool	FindIsotope (G4IsotopeProperty *property)
G4IsotopeProperty *	GetIsotope (G4int Z, G4int A, G4double E, G4Ions::G4FloatLevelBase flb=G4Ions::G4FloatLevelBase::no_Float) override
G4IsotopeProperty *	GetIsotopeByIsoLvl (G4int Z, G4int A, G4int lvl=0) override
Public Member Functions inherited from G4VIsotopeTable
	G4VIsotopeTable ()=default
	G4VIsotopeTable (const G4String &)
	G4VIsotopeTable (const G4VIsotopeTable &)=default
G4VIsotopeTable &	operator= (const G4VIsotopeTable &)
virtual	~G4VIsotopeTable ()=default
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

IsomerLevel is added 30 Apr. 2013 H.Kurashige

Definition at line 49 of file G4IsotopeMagneticMomentTable.hh.

Member Typedef Documentation

G4IsotopeList

using G4IsotopeMagneticMomentTable::G4IsotopeList = std::vector<G4IsotopeProperty*>

Definition at line 57 of file G4IsotopeMagneticMomentTable.hh.

G4IsotopeNameList

using G4IsotopeMagneticMomentTable::G4IsotopeNameList = std::vector<G4String>

Definition at line 58 of file G4IsotopeMagneticMomentTable.hh.

Constructor & Destructor Documentation

G4IsotopeMagneticMomentTable() [1/2]

G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable

(

)

Definition at line 55 of file G4IsotopeMagneticMomentTable.cc.

                                                           : G4VIsotopeTable("MagneticMoment")
{
  if (std::getenv("G4IONMAGNETICMOMENT") == nullptr) {
#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 = std::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()) {  // Loop checking, 09.08.2015, K.Kurashige
    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;
 
      auto 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::~G4IsotopeMagneticMomentTable ( )

override

Definition at line 120 of file G4IsotopeMagneticMomentTable.cc.

{
  for (const auto& i : fIsotopeList) {
    delete i;
  }
  fIsotopeList.clear();
}

G4IsotopeMagneticMomentTable() [2/2]

G4IsotopeMagneticMomentTable::G4IsotopeMagneticMomentTable ( const G4IsotopeMagneticMomentTable & right )

protected

Definition at line 128 of file G4IsotopeMagneticMomentTable.cc.

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

Member Function Documentation

FindIsotope()

G4bool G4IsotopeMagneticMomentTable::FindIsotope ( G4IsotopeProperty * property )

virtual

Definition at line 139 of file G4IsotopeMagneticMomentTable.cc.

{
  for (const auto fP : fIsotopeList) {
    // 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;
}

GetIsotope()

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

overridevirtual

Implements G4VIsotopeTable.

Definition at line 173 of file G4IsotopeMagneticMomentTable.cc.

{
  G4IsotopeProperty* fProperty = nullptr;
  for (const auto fP : fIsotopeList) {
    // 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;
}

GetIsotopeByIsoLvl()

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

overridevirtual

Reimplemented from G4VIsotopeTable.

Definition at line 206 of file G4IsotopeMagneticMomentTable.cc.

{
  G4IsotopeProperty* fProperty = nullptr;
  for (const auto fP : fIsotopeList) {
    // 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;
}

operator=()

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

protected

Definition at line 134 of file G4IsotopeMagneticMomentTable.cc.

{
  return *this;
}

---

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

• G4IsotopeMagneticMomentTable.hh
• G4IsotopeMagneticMomentTable.cc