G4RDShellData Class Reference

#include <G4RDShellData.hh>

Public Member Functions
	G4RDShellData (G4int minZ=1, G4int maxZ=100, G4bool isOccupancy=false)
	~G4RDShellData ()
size_t	NumberOfShells (G4int Z) const
G4int	ShellId (G4int Z, G4int shellIndex) const
G4double	ShellOccupancyProbability (G4int Z, G4int shellIndex) const
const std::vector< G4double > &	ShellIdVector (G4int Z) const
G4double	BindingEnergy (G4int Z, G4int shellIndex) const
void	SetOccupancyData ()
void	LoadData (const G4String &fileName)
void	PrintData () const
G4int	SelectRandomShell (G4int Z) const

Detailed Description

Definition at line 53 of file G4RDShellData.hh.

Constructor & Destructor Documentation

G4RDShellData::G4RDShellData	(	G4int	minZ = 1,
		G4int	maxZ = 100,
		G4bool	isOccupancy = false
	)

Definition at line 54 of file G4RDShellData.cc.

  : zMin(minZ), zMax(maxZ), occupancyData(isOccupancy)
{  }

G4RDShellData::~G4RDShellData

(

)

Definition at line 59 of file G4RDShellData.cc.

{
  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos;
  for (pos = idMap.begin(); pos != idMap.end(); ++pos)
    {
      std::vector<G4double>* dataSet = (*pos).second;
      delete dataSet;
    }

  std::map<G4int,G4DataVector*,std::less<G4int> >::iterator pos2;
  for (pos2 = bindingMap.begin(); pos2 != bindingMap.end(); ++pos2)
    {
      G4DataVector* dataSet = (*pos2).second;
      delete dataSet;
    }

  if (occupancyData)
    {
      std::map<G4int,std::vector<G4double>*,std::less<G4int> >::iterator pos3;
      for (pos3 = occupancyPdfMap.begin(); pos3 != occupancyPdfMap.end(); ++pos3)
    {
      std::vector<G4double>* dataSet = (*pos3).second;
      delete dataSet;
    }
    }
}

Member Function Documentation

G4double G4RDShellData::BindingEnergy	(	G4int	Z,
		G4int	shellIndex
	)		const

Definition at line 169 of file G4RDShellData.cc.

Referenced by G4RDAtomicTransitionManager::G4RDAtomicTransitionManager(), G4LowEnergyCompton::PostStepDoIt(), and G4LowEnergyPolarizedCompton::PostStepDoIt().

{
  G4double value = 0.;

  if (Z >= zMin && Z <= zMax)
    {
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator pos;
      pos = bindingMap.find(Z);
      if (pos!= bindingMap.end())
    {
      G4DataVector dataSet = *((*pos).second);
      G4int nData = dataSet.size();
      if (shellIndex >= 0 && shellIndex < nData)
        {
          value = dataSet[shellIndex];
        }
    }
    }
  return value;
}

void G4RDShellData::LoadData

(

const G4String &

fileName

)

Definition at line 237 of file G4RDShellData.cc.

References test::a, FatalException, G4Exception(), python.hepunit::MeV, n, and ShellIdVector().

Referenced by G4LowEnergyCompton::BuildPhysicsTable(), G4LowEnergyPolarizedCompton::BuildPhysicsTable(), and G4RDAtomicTransitionManager::G4RDAtomicTransitionManager().

{ 
  // Build the complete string identifying the file with the data set
  
  std::ostringstream ost;
  
  ost << fileName << ".dat";
  
  G4String name(ost.str());
  
  char* path = getenv("G4LEDATA");
  if (!path)
    { 
      G4String excep("G4LEDATA environment variable not set!");
      G4Exception("G4RDShellData::LoadData()", "InvalidSetup",
                  FatalException, excep);
    }
  
  G4String pathString(path);
  G4String dirFile = pathString + name;
  std::ifstream file(dirFile);
  std::filebuf* lsdp = file.rdbuf();

  if (! (lsdp->is_open()) )
    {
      G4String s1("Data file: ");
      G4String s2(" not found");
      G4String excep = s1 + dirFile + s2;
      G4Exception("G4RDShellData::LoadData()", "DataNotFound",
                  FatalException, excep);
    }

  G4double a = 0;
  G4int k = 1;
  G4int s = 0;
  
  G4int Z = 1;
  G4DataVector* energies = new G4DataVector;
  std::vector<G4double>* ids = new std::vector<G4double>;

  do {
    file >> a;
    G4int nColumns = 2;
    if (a == -1)
      {
    if (s == 0)
      {
        // End of a shell data set
        idMap[Z] = ids;
            bindingMap[Z] = energies;
            G4int n = ids->size();
        nShells.push_back(n);
        // Start of new shell data set
        ids = new std::vector<G4double>;
            energies = new G4DataVector;
            Z++;        
      }      
    s++;
    if (s == nColumns)
    {
      s = 0;
    }
      }
    else if (a == -2)
      {
    // End of file; delete the empty vectors created when encountering the last -1 -1 row
    delete energies;
    delete ids;
    //nComponents = components.size();
      }
    else
      {
    // 1st column is shell id
    if(k%nColumns != 0)
      {     
        ids->push_back(a);
        k++;
      }
    else if (k%nColumns == 0)
      {
        // 2nd column is binding energy
        G4double e = a * MeV;
        energies->push_back(e);
        k = 1;
      }
      }
  } while (a != -2); // end of file
  file.close();    

  // For Doppler broadening: the data set contains shell occupancy and binding energy for each shell
  // Build additional map with probability for each shell based on its occupancy

  if (occupancyData)
    {
      // Build cumulative from raw shell occupancy

      for (G4int Z=zMin; Z <= zMax; Z++)
    {
      std::vector<G4double> occupancy = ShellIdVector(Z);

      std::vector<G4double>* prob = new std::vector<G4double>;
      G4double scale = 1. / G4double(Z);

      prob->push_back(occupancy[0] * scale);
      for (size_t i=1; i<occupancy.size(); i++)
        {
          prob->push_back(occupancy[i]*scale + (*prob)[i-1]);
        }
      occupancyPdfMap[Z] = prob;

      /*
        G4double scale = 1. / G4double(Z);
        //      transform((*prob).begin(),(*prob).end(),(*prob).begin(),bind2nd(multiplies<G4double>(),scale));

        for (size_t i=0; i<occupancy.size(); i++)
        {
        (*prob)[i] *= scale;
        }
      */
    }
    }
}

size_t G4RDShellData::NumberOfShells

(

G4int

Z

)

const

Definition at line 87 of file G4RDShellData.cc.

References n, and z.

Referenced by G4RDAtomicTransitionManager::G4RDAtomicTransitionManager(), and SelectRandomShell().

{
  G4int z = Z - 1;
  G4int n = 0;

  if (Z>= zMin && Z <= zMax)
    {
      n = nShells[z];
    }
  return n;
}

void G4RDShellData::PrintData

(

void

)

const

Definition at line 190 of file G4RDShellData.cc.

References G4cout, G4endl, and python.hepunit::keV.

{
  for (G4int Z = zMin; Z <= zMax; Z++)
    {
      G4cout << "---- Shell data for Z = "
         << Z
         << " ---- "
         << G4endl;
      G4int nSh = nShells[Z-1];
      std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator posId;
      posId = idMap.find(Z);
      std::vector<G4double>* ids = (*posId).second;
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator posE;
      posE = bindingMap.find(Z);
      G4DataVector* energies = (*posE).second;
      for (G4int i=0; i<nSh; i++)
    {
      G4int id = (G4int) (*ids)[i];
      G4double e = (*energies)[i] / keV;
      G4cout << i << ") ";

      if (occupancyData) 
        {
          G4cout << " Occupancy: ";
        }
      else 
        {
          G4cout << " Shell id: ";
        }
      G4cout << id << " - Binding energy = "
         << e << " keV ";
        if (occupancyData)
          {
        std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator posOcc;
        posOcc = occupancyPdfMap.find(Z);
                std::vector<G4double> probs = *((*posOcc).second);
                G4double prob = probs[i];
        G4cout << "- Probability = " << prob;
          }
        G4cout << G4endl;
    }
      G4cout << "-------------------------------------------------" 
         << G4endl;
    }
}

G4int G4RDShellData::SelectRandomShell

(

G4int

Z

)

const

Definition at line 361 of file G4RDShellData.cc.

References FatalException, G4Exception(), G4UniformRand, and NumberOfShells().

Referenced by G4LowEnergyCompton::PostStepDoIt(), and G4LowEnergyPolarizedCompton::PostStepDoIt().

{
  if (Z < zMin || Z > zMax)
    G4Exception("G4RDShellData::SelectRandomShell()", "OutOfRange",
                FatalException, "Z outside boundaries!");

  G4int shellIndex = 0;    
  std::vector<G4double> prob = ShellVector(Z);
  G4double random = G4UniformRand();

  // std::vector<G4double>::const_iterator pos;
  // pos = lower_bound(prob.begin(),prob.end(),random);

  // Binary search the shell with probability less or equal random

  G4int nShells = NumberOfShells(Z);
  G4int upperBound = nShells;

  while (shellIndex <= upperBound) 
    {
      G4int midShell = (shellIndex + upperBound) / 2;
      if ( random < prob[midShell] ) 
    upperBound = midShell - 1;
      else 
    shellIndex = midShell + 1;
    }  
  if (shellIndex >= nShells) shellIndex = nShells - 1;

  return shellIndex;
}

void G4RDShellData::SetOccupancyData ( )

inline

Definition at line 71 of file G4RDShellData.hh.

Referenced by G4LowEnergyCompton::G4LowEnergyCompton(), and G4LowEnergyPolarizedCompton::G4LowEnergyPolarizedCompton().

{ occupancyData = true; }

G4int G4RDShellData::ShellId	(	G4int	Z,
		G4int	shellIndex
	)		const

Definition at line 124 of file G4RDShellData.cc.

References n.

Referenced by G4RDAtomicTransitionManager::G4RDAtomicTransitionManager().

{
  G4int n = -1;

  if (Z >= zMin && Z <= zMax)
    {
      std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
      pos = idMap.find(Z);
      if (pos!= idMap.end())
    {
      std::vector<G4double> dataSet = *((*pos).second);
      G4int nData = dataSet.size();
      if (shellIndex >= 0 && shellIndex < nData)
        {
          n = (G4int) dataSet[shellIndex];
        }
    }
    }
  return n;
}

const std::vector< G4double > & G4RDShellData::ShellIdVector

(

G4int

Z

)

const

Definition at line 100 of file G4RDShellData.cc.

References FatalException, and G4Exception().

Referenced by LoadData().

{
  std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
  if (Z < zMin || Z > zMax)
    G4Exception("G4RDShellData::ShellIdVector()", "OutOfRange",
                FatalException, "Z outside boundaries!");
  pos = idMap.find(Z);
  std::vector<G4double>* dataSet = (*pos).second;
  return *dataSet;
}

G4double G4RDShellData::ShellOccupancyProbability	(	G4int	Z,
		G4int	shellIndex
	)		const

Definition at line 146 of file G4RDShellData.cc.

{
  G4double prob = -1.;

  if (Z >= zMin && Z <= zMax)
    {
      std::map<G4int,std::vector<G4double>*,std::less<G4int> >::const_iterator pos;
      pos = idMap.find(Z);
      if (pos!= idMap.end())
    {
      std::vector<G4double> dataSet = *((*pos).second);
      G4int nData = dataSet.size();
      if (shellIndex >= 0 && shellIndex < nData)
        {
          prob = dataSet[shellIndex];
        }
    }
    }
  return prob;
}

---

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

• G4RDShellData.hh
• G4RDShellData.cc