genie::NaturalIsotopes Class Reference

Singleton class to load & serve tables of natural occurring isotopes. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/genie/v3_00_06_p01/Linux64bit+3.10-2.17-e17-debug/GENIE-Generator/src/Framework/ParticleData/NaturalIsotopes.h"

Classes
struct	Cleaner

Public Member Functions
int	NElements (int Z) const
const NaturalIsotopeElementData *	ElementData (int Z, int ielement) const

Static Public Member Functions
static NaturalIsotopes *	Instance (void)

Private Member Functions
	NaturalIsotopes ()
	NaturalIsotopes (const NaturalIsotopes &)
virtual	~NaturalIsotopes ()
bool	LoadTable (void)

Private Attributes
map< int, vector< NaturalIsotopeElementData * > >	fNaturalIsotopesTable

Static Private Attributes
static NaturalIsotopes *	fInstance = 0

Friends
struct	Cleaner

Detailed Description

Singleton class to load & serve tables of natural occurring isotopes.

Author

Costas Andreopoulos <costas.andreopoulos stfc.ac.uk> University of Liverpool & STFC Rutherford Appleton Lab

Jim Dobson j.dob.nosp@m.son0.nosp@m.7@imp.nosp@m.eria.nosp@m.l.ac..nosp@m.uk Imperial College London

May 30, 2008

Definition at line 31 of file NaturalIsotopes.h.

Constructor & Destructor Documentation

NaturalIsotopes::NaturalIsotopes ( )

private

Definition at line 40 of file NaturalIsotopes.cxx.

References fInstance, LoadTable(), LOG, and pERROR.

Referenced by Instance().

{
  if( ! this->LoadTable() ) {
     LOG("NatIsotop", pERROR) << "NaturalIsotopes initialization failed!";
  }
  fInstance =  0;
}

genie::NaturalIsotopes::NaturalIsotopes ( const NaturalIsotopes &  )

private

NaturalIsotopes::~NaturalIsotopes ( )

privatevirtual

Definition at line 48 of file NaturalIsotopes.cxx.

References om::cout, allTimeWatchdog::endl, fInstance, and fNaturalIsotopesTable.

{
  cout << "NaturalIsotopes singleton dtor: "
       << "Deleting natural isotope data tables" << endl;

  map<int, vector<NaturalIsotopeElementData*> >::iterator miter;
  vector<NaturalIsotopeElementData*>::iterator viter;

  for(miter = fNaturalIsotopesTable.begin(); 
                 miter != fNaturalIsotopesTable.end(); ++miter) {
     vector<NaturalIsotopeElementData*> vec = miter->second;
     for(viter = vec.begin(); viter != vec.end(); ++viter) {
        NaturalIsotopeElementData * element_data = *viter;
        if(element_data) {
            delete element_data;
            element_data = 0;
        }
     }
     vec.clear();
  }
  fNaturalIsotopesTable.clear();
  fInstance = 0;
}

Member Function Documentation

const NaturalIsotopeElementData * NaturalIsotopes::ElementData	(	int	Z,
		int	ielement
	)		const

Definition at line 99 of file NaturalIsotopes.cxx.

References fNaturalIsotopesTable, LOG, pWARN, and Z.

Referenced by main().

{
  map<int, vector<NaturalIsotopeElementData*> >::const_iterator miter;

  if( (miter=fNaturalIsotopesTable.find(Z)) == fNaturalIsotopesTable.end()) {
    LOG("NatIsotop", pWARN)  
       << "Table has no elements for natural isotope  Z = " << Z;
    return 0;
  }
  vector<NaturalIsotopeElementData*> vec = miter->second;
  if(ielement >= (int)vec.size() || ielement < 0) {
    LOG("NatIsotop", pWARN)  
      << "Natural isotope Z = " << Z << " has " << vec.size() << " elements"
      << " (element = " << ielement << " was requested)";
    return 0;
  }
  return vec[ielement];
}

NaturalIsotopes * NaturalIsotopes::Instance ( void  )

static

Definition at line 72 of file NaturalIsotopes.cxx.

References genie::NaturalIsotopes::Cleaner::DummyMethodAndSilentCompiler(), fInstance, LOG, NaturalIsotopes(), and pINFO.

Referenced by main().

{
  if(fInstance == 0) {
    LOG("NatIsotop", pINFO) << "NaturalIsotopes late initialization";

    static NaturalIsotopes::Cleaner cleaner;
    cleaner.DummyMethodAndSilentCompiler();

    fInstance = new NaturalIsotopes;
  }
  return fInstance;
}

bool NaturalIsotopes::LoadTable ( void  )

private

eof

Definition at line 118 of file NaturalIsotopes.cxx.

References shutoffs::filename, fNaturalIsotopesTable, MECModelEnuComparisons::i, LOG, getGoodRuns4SAM::n, pDEBUG, pINFO, pWARN, string, and Z.

Referenced by NaturalIsotopes().

{
  // get the natural isotopes table filename
  string filename = string(gSystem->Getenv("GENIE")) + 
                    string("/data/evgen/catalogues/iso/natural-isotopes.data");

  LOG("NatIsotop", pINFO)  
    << "Loading natural occurring isotope table from file: " << filename;

  bool is_accessible = ! (gSystem->AccessPathName( filename.c_str() ));
  if (!is_accessible) {
     LOG("NatIsotop", pWARN) << "Can not read file: " << filename;
     return false;
  }

  // load the natural isotopes .txt file
  string input_buf;
  std::ifstream input(filename.c_str());  
  if (input.is_open()){
    
    //skip first 8 lines (comments)
    for(int i=0; i<8; i++){ 
      string buffer;
      getline(input, buffer);
    } 

    int Z = -1, Z_previous = -1, nelements = 0, pdgcode = 0;
    double atomicmass = 0, abundance = 0;
    string elementname, subelementname;
    
    while( !input.eof() ) {
      
      //read in naturally occuring element info
      input >> Z;
      input >> elementname; 
      input >> nelements; 
      
      vector<NaturalIsotopeElementData *> vec;
      NaturalIsotopeElementData * data = 0;
      
      // check not re-reading same element
      if(Z!=Z_previous){    
        LOG("NatIsotop", pDEBUG) << "Reading entry for Z = " << Z;
        for(int n=0 ; n < nelements; n++){
          input >> subelementname;
          input >> pdgcode;
          input >> atomicmass;
          input >> abundance;
          LOG("NatIsotop", pDEBUG) 
            << " - Element: " << n << ", pdg = " << pdgcode 
            << ", A = " << atomicmass << ", abundance = " << abundance;
          data = new NaturalIsotopeElementData(pdgcode, abundance);
          vec.push_back(data);
        }
        fNaturalIsotopesTable.insert(
           map<int,vector<NaturalIsotopeElementData*> >::value_type(Z,vec));
      }      
      Z_previous = Z;
    } //!eof

  } else { 
    return false;       
  } //open?
  
  return true;
}

int NaturalIsotopes::NElements

(

int

Z

)

const

Definition at line 85 of file NaturalIsotopes.cxx.

References fNaturalIsotopesTable, LOG, pWARN, and Z.

Referenced by main().

{
  map<int, vector<NaturalIsotopeElementData*> >::const_iterator miter;

  if( (miter=fNaturalIsotopesTable.find(Z)) == fNaturalIsotopesTable.end()) {
    LOG("NatIsotop", pWARN)  
       << "Table has no elements for natural isotope  Z = " << Z;
    return 0;
  } 
  vector<NaturalIsotopeElementData*> vec = miter->second;
  return vec.size();
}

Friends And Related Function Documentation

friend struct Cleaner

friend

Definition at line 59 of file NaturalIsotopes.h.

Member Data Documentation

NaturalIsotopes * NaturalIsotopes::fInstance = 0

staticprivate

Definition at line 46 of file NaturalIsotopes.h.

Referenced by Instance(), NaturalIsotopes(), genie::NaturalIsotopes::Cleaner::~Cleaner(), and ~NaturalIsotopes().

map<int, vector<NaturalIsotopeElementData*> > genie::NaturalIsotopes::fNaturalIsotopesTable

private

Definition at line 48 of file NaturalIsotopes.h.

Referenced by ElementData(), LoadTable(), NElements(), and ~NaturalIsotopes().

---

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

• Framework/ParticleData/NaturalIsotopes.h
• Framework/ParticleData/NaturalIsotopes.cxx