G4ShellEMDataSet Class Reference

#include <G4ShellEMDataSet.hh>

Inheritance diagram for G4ShellEMDataSet:

Public Member Functions
void	AddComponent (G4VEMDataSet *dataSet) override
G4double	FindValue (G4double energy, G4int componentId=0) const override
	G4ShellEMDataSet (const G4ShellEMDataSet &copy)=delete
	G4ShellEMDataSet (G4int Z, G4VDataSetAlgorithm *algo, G4double eUnit=CLHEP::MeV, G4double dataUnit=CLHEP::barn)
const G4VEMDataSet *	GetComponent (G4int componentId) const override
const G4DataVector &	GetData (G4int componentId) const override
const G4DataVector &	GetEnergies (G4int componentId) const override
const G4DataVector &	GetLogData (G4int componentId) const override
const G4DataVector &	GetLogEnergies (G4int componentId) const override
G4bool	LoadData (const G4String &fileName) override
G4bool	LoadNonLogData (const G4String &fileName) override
size_t	NumberOfComponents (void) const override
G4ShellEMDataSet &	operator= (const G4ShellEMDataSet &right)=delete
void	PrintData (void) const override
G4double	RandomSelect (G4int) const override
G4bool	SaveData (const G4String &fileName) const override
void	SetEnergiesData (G4DataVector *energies, G4DataVector *data, G4int componentId) override
void	SetLogEnergiesData (G4DataVector *energies, G4DataVector *data, G4DataVector *log_energies, G4DataVector *log_data, G4int componentId) override
virtual	~G4ShellEMDataSet ()

Protected Member Functions
void	CleanUpComponents ()
const G4VDataSetAlgorithm *	GetAlgorithm () const
G4double	GetUnitData () const
G4double	GetUnitEnergies () const

Private Member Functions
G4String	FullFileName (const G4String &fileName) const

Private Attributes
G4VDataSetAlgorithm *	algorithm
std::vector< G4VEMDataSet * >	components
G4double	unitData
G4double	unitEnergies
G4int	z

Detailed Description

Definition at line 59 of file G4ShellEMDataSet.hh.

Constructor & Destructor Documentation

G4ShellEMDataSet() [1/2]

G4ShellEMDataSet::G4ShellEMDataSet ( G4int  Z, G4VDataSetAlgorithm *  algo, G4double  eUnit = CLHEP::MeV, G4double  dataUnit = CLHEP::barn  )

explicit

Definition at line 63 of file G4ShellEMDataSet.cc.

  : 
  algorithm(algo),
  unitEnergies(eUnit),
  unitData(dataUnit),
  z(zeta)
{
  if (algorithm == nullptr) 
    G4Exception("G4ShellEMDataSet::G4ShellEMDataSet()","em0007",
        FatalErrorInArgument, "Interpolation == 0");
}

References algorithm, FatalErrorInArgument, and G4Exception().

~G4ShellEMDataSet()

G4ShellEMDataSet::~G4ShellEMDataSet ( )

virtual

Definition at line 78 of file G4ShellEMDataSet.cc.

{
  CleanUpComponents();
  if (algorithm) delete algorithm;
}

References algorithm, and CleanUpComponents().

G4ShellEMDataSet() [2/2]

G4ShellEMDataSet::G4ShellEMDataSet ( const G4ShellEMDataSet &  copy )

delete

Member Function Documentation

AddComponent()

void G4ShellEMDataSet::AddComponent ( G4VEMDataSet *  dataSet )

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 73 of file G4ShellEMDataSet.hh.

{ components.push_back(dataSet); }

References components.

Referenced by LoadData(), and LoadNonLogData().

CleanUpComponents()

void G4ShellEMDataSet::CleanUpComponents ( void  )

protected

Definition at line 399 of file G4ShellEMDataSet.cc.

{
  while (!components.empty())
    {
      if (components.back()) delete components.back();
      components.pop_back();
    }
}

References components.

Referenced by LoadData(), LoadNonLogData(), and ~G4ShellEMDataSet().

FindValue()

G4double G4ShellEMDataSet::FindValue ( G4double  energy, G4int  componentId = 0  ) const

overridevirtual

Implements G4VEMDataSet.

Definition at line 85 of file G4ShellEMDataSet.cc.

{
  // Returns the sum over the shells corresponding to e
  G4double value = 0.;
 
  std::vector<G4VEMDataSet *>::const_iterator i(components.begin());
  std::vector<G4VEMDataSet *>::const_iterator end(components.end());
 
  while (i != end)
    {
      value += (*i)->FindValue(energy);
      i++;
    }
 
  return value;
}

References components, and G4INCL::KinematicsUtils::energy().

FullFileName()

G4String G4ShellEMDataSet::FullFileName ( const G4String &  fileName ) const

private

Definition at line 410 of file G4ShellEMDataSet.cc.

{
  char* path = std::getenv("G4LEDATA");
  
  if (!path)
  {
    G4Exception("G4ShellEMDataSet::FullFileName()","em0006",JustWarning,"Please set G4LEDATA");
    return "";
  }
  
  std::ostringstream fullFileName;
 
  fullFileName << path << '/' << fileName << z << ".dat";
                      
  return G4String(fullFileName.str().c_str());
}

References G4Exception(), JustWarning, and z.

Referenced by LoadData(), LoadNonLogData(), and SaveData().

GetAlgorithm()

const G4VDataSetAlgorithm * G4ShellEMDataSet::GetAlgorithm ( ) const

inlineprotected

Definition at line 100 of file G4ShellEMDataSet.hh.

{ return algorithm; }   

References algorithm.

GetComponent()

const G4VEMDataSet * G4ShellEMDataSet::GetComponent ( G4int  componentId ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 72 of file G4ShellEMDataSet.hh.

{ return components[componentId]; }

References components.

Referenced by GetData(), GetEnergies(), GetLogData(), GetLogEnergies(), PrintData(), and SaveData().

GetData()

const G4DataVector & G4ShellEMDataSet::GetData ( G4int  componentId ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 77 of file G4ShellEMDataSet.hh.

{ return GetComponent(componentId)->GetData(0); }

References GetComponent(), and G4VEMDataSet::GetData().

GetEnergies()

const G4DataVector & G4ShellEMDataSet::GetEnergies ( G4int  componentId ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 76 of file G4ShellEMDataSet.hh.

{ return GetComponent(componentId)->GetEnergies(0); }

References GetComponent(), and G4VEMDataSet::GetEnergies().

GetLogData()

const G4DataVector & G4ShellEMDataSet::GetLogData ( G4int  componentId ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 80 of file G4ShellEMDataSet.hh.

{ return GetComponent(componentId)->GetLogData(0); }

References GetComponent(), and G4VEMDataSet::GetLogData().

GetLogEnergies()

const G4DataVector & G4ShellEMDataSet::GetLogEnergies ( G4int  componentId ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 78 of file G4ShellEMDataSet.hh.

  { return GetComponent(componentId)->GetLogEnergies(0); }

References GetComponent(), and G4VEMDataSet::GetLogEnergies().

GetUnitData()

G4double G4ShellEMDataSet::GetUnitData ( ) const

inlineprotected

Definition at line 99 of file G4ShellEMDataSet.hh.

{ return unitData; }

References unitData.

GetUnitEnergies()

G4double G4ShellEMDataSet::GetUnitEnergies ( ) const

inlineprotected

Definition at line 98 of file G4ShellEMDataSet.hh.

{ return unitEnergies; }

References unitEnergies.

LoadData()

G4bool G4ShellEMDataSet::LoadData ( const G4String &  fileName )

overridevirtual

Implements G4VEMDataSet.

Definition at line 165 of file G4ShellEMDataSet.cc.

{
  CleanUpComponents();
 
  G4String fullFileName = FullFileName(file);
  std::ifstream in(fullFileName);
 
  if (!in.is_open())
    {
      G4String message("Data file \"");
      message += fullFileName;
      message += "\" not found";
      G4Exception("G4ShellEMDataSet::LoadData()", "em0003",FatalException, message);
      return 0;
    }
  
  G4DataVector* orig_shell_energies = 0;
  G4DataVector* orig_shell_data = 0;
  G4DataVector* log_shell_energies = 0;
  G4DataVector* log_shell_data = 0;
 
  G4double a = 0.;
  G4int shellIndex = 0;
  G4int k = 0;
  G4int nColumns = 2;
 
  do
    {
      in >> a;
 
      if (a==0.) a=1e-300;
 
      // The file is organized into four columns:
      // 1st column contains the values of energy
      // 2nd column contains the corresponding data value
      // The file terminates with the pattern: -1   -1
      //                                       -2   -2
      //
      if (a == -1)
    {
      if ((k%nColumns == 0) && (orig_shell_energies != 0) )
        {
         AddComponent(new G4EMDataSet(shellIndex, orig_shell_energies, 
                      orig_shell_data, log_shell_energies, 
                      log_shell_data, algorithm->Clone(), 
                      unitEnergies, unitData));
          orig_shell_energies = 0;
          orig_shell_data = 0;
              log_shell_energies = 0;
              log_shell_data = 0;
        }
    }
      else if (a != -2)
    {
      if (orig_shell_energies == 0)
        {
         orig_shell_energies = new G4DataVector;
         orig_shell_data = new G4DataVector;
             log_shell_energies = new G4DataVector;
             log_shell_data = new G4DataVector;
        }
      if (k%nColumns == 0)
            {
         orig_shell_energies->push_back(a*unitEnergies);
         log_shell_energies->push_back(std::log10(a) + std::log10(unitEnergies));
            }
      else if (k%nColumns == 1)
            {
         orig_shell_data->push_back(a*unitData);
             log_shell_data->push_back(std::log10(a) + std::log10(unitData));
            }
          k++;
    }
      else k = 1;
    }
  while (a != -2);  // End of file
 
  delete orig_shell_energies;
  delete orig_shell_data;
  delete log_shell_energies;
  delete log_shell_data;
 
  return true;
}

References AddComponent(), algorithm, CleanUpComponents(), G4VDataSetAlgorithm::Clone(), FatalException, geant4_check_module_cycles::file, FullFileName(), G4Exception(), unitData, and unitEnergies.

LoadNonLogData()

G4bool G4ShellEMDataSet::LoadNonLogData ( const G4String &  fileName )

overridevirtual

Implements G4VEMDataSet.

Definition at line 252 of file G4ShellEMDataSet.cc.

{
  CleanUpComponents();
 
  G4String fullFileName = FullFileName(file);
  std::ifstream in(fullFileName);
 
  if (!in.is_open())
    {
      G4String message("G4ShellEMDataSet::LoadData - data file \"");
      message += fullFileName;
      message += "\" not found";
      G4Exception("G4ShellEMDataSet::LoadNonLogData()", "em0003",FatalException, message);
      return 0;
    }
 
  G4DataVector* orig_shell_energies = 0;
  G4DataVector* orig_shell_data = 0;
 
  G4double a = 0.;
  G4int shellIndex = 0;
  G4int k = 0;
  G4int nColumns = 2;
 
  do
    {
      in >> a;
      // The file is organized into four columns:
      // 1st column contains the values of energy
      // 2nd column contains the corresponding data value
      // The file terminates with the pattern: -1   -1
      //                                       -2   -2
      //
      if (a == -1)
    {
      if ((k%nColumns == 0) && (orig_shell_energies != 0) )
        {
         AddComponent(new G4EMDataSet(shellIndex, orig_shell_energies, 
                      orig_shell_data, algorithm->Clone(), 
                      unitEnergies, unitData));
          orig_shell_energies = 0;
          orig_shell_data = 0;
        }
    }
      else if (a != -2)
    {
      if (orig_shell_energies == 0)
        {
         orig_shell_energies = new G4DataVector;
         orig_shell_data = new G4DataVector;
        }
      if (k%nColumns == 0)
            {
         orig_shell_energies->push_back(a*unitEnergies);
            }
      else if (k%nColumns == 1)
            {
         orig_shell_data->push_back(a*unitData);
            }
          k++;
    }
      else k = 1;
    }
  while (a != -2);  // End of file
 
 
  delete orig_shell_energies;
  delete orig_shell_data;
 
  return true;
}

References AddComponent(), algorithm, CleanUpComponents(), G4VDataSetAlgorithm::Clone(), FatalException, geant4_check_module_cycles::file, FullFileName(), G4Exception(), unitData, and unitEnergies.

NumberOfComponents()

size_t G4ShellEMDataSet::NumberOfComponents ( void  ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 74 of file G4ShellEMDataSet.hh.

{ return components.size(); }

References components.

Referenced by PrintData(), and SaveData().

operator=()

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

delete

PrintData()

void G4ShellEMDataSet::PrintData ( void  ) const

overridevirtual

Implements G4VEMDataSet.

Definition at line 104 of file G4ShellEMDataSet.cc.

{
  const size_t n = NumberOfComponents();
 
  G4cout << "The data set has " << n << " components" << G4endl;
  G4cout << G4endl;
 
  size_t i = 0;
 
  while (i < n)
    {
      G4cout << "--- Component " << i << " ---" << G4endl;
      GetComponent(i)->PrintData();
      i++;
    }
}

References G4cout, G4endl, GetComponent(), CLHEP::detail::n, NumberOfComponents(), and G4VEMDataSet::PrintData().

RandomSelect()

G4double G4ShellEMDataSet::RandomSelect ( G4int  ) const

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 92 of file G4ShellEMDataSet.hh.

{ return -1.; };

SaveData()

G4bool G4ShellEMDataSet::SaveData ( const G4String &  fileName ) const

overridevirtual

Implements G4VEMDataSet.

Definition at line 326 of file G4ShellEMDataSet.cc.

{
  G4String fullFileName = FullFileName(file);
  std::ofstream out(fullFileName);
 
  if (!out.is_open())
    {
      G4String message("Cannot open \"");
      message += fullFileName;
      message += "\"";
      G4Exception("G4EMDataSet::SaveData()","em0005",FatalException,message);
    }
 
  const size_t n = NumberOfComponents();
  size_t k = 0;
 
  while (k < n)
    {
      const G4VEMDataSet* component = GetComponent(k);
  
      if (component)
    {
      const G4DataVector& energies = component->GetEnergies(0);
      const G4DataVector& data = component->GetData(0);
 
      G4DataVector::const_iterator i = energies.begin();
      G4DataVector::const_iterator endI = energies.end();
      G4DataVector::const_iterator j = data.begin();
  
      while (i != endI)
        {
          out.precision(10);
          out.width(15);
          out.setf(std::ofstream::left);
          out << ((*i)/unitEnergies) << ' ';
 
          out.precision(10);
          out.width(15);
          out.setf(std::ofstream::left);
          out << ((*j)/unitData) << std::endl;
          i++;
          j++;
        }
    }
  
      out.precision(10);
      out.width(15);
      out.setf(std::ofstream::left);
      out << -1.f << ' ';
 
      out.precision(10);
      out.width(15);
      out.setf(std::ofstream::left);
      out << -1.f << std::endl;
  
      k++;
    }
 
  out.precision(10);
  out.width(15);
  out.setf(std::ofstream::left);
  out << -2.f << ' ';
 
  out.precision(10);
  out.width(15);
  out.setf(std::ofstream::left);
  out << -2.f << std::endl;
 
  return true;
}

References FatalException, geant4_check_module_cycles::file, FullFileName(), G4Exception(), GetComponent(), G4VEMDataSet::GetData(), G4VEMDataSet::GetEnergies(), CLHEP::detail::n, NumberOfComponents(), unitData, and unitEnergies.

SetEnergiesData()

void G4ShellEMDataSet::SetEnergiesData ( G4DataVector *  energies, G4DataVector *  data, G4int  componentId  )

overridevirtual

Implements G4VEMDataSet.

Definition at line 123 of file G4ShellEMDataSet.cc.

{
  G4VEMDataSet* component = components[componentId]; 
  if (component)
    {
      component->SetEnergiesData(energies, data, 0);
      return;
    }
 
  G4String msg = "component ";
  msg += static_cast<char>(componentId);
  msg += " not found";
 
  G4Exception("G4ShellEMDataSet::SetEnergiesData()","em0008", FatalErrorInArgument ,msg);
}

References components, FatalErrorInArgument, G4Exception(), and G4VEMDataSet::SetEnergiesData().

SetLogEnergiesData()

void G4ShellEMDataSet::SetLogEnergiesData ( G4DataVector *  energies, G4DataVector *  data, G4DataVector *  log_energies, G4DataVector *  log_data, G4int  componentId  )

overridevirtual

Implements G4VEMDataSet.

Definition at line 143 of file G4ShellEMDataSet.cc.

{
  G4VEMDataSet* component = components[componentId]; 
  if (component)
    {
      component->SetLogEnergiesData(energies, data, log_energies, log_data, 0);
      return;
    }
 
  G4String msg = "component ";
  msg += static_cast<char>(componentId);
  msg += " not found";
  G4Exception("G4ShellEMDataSet::SetLogEnergiesData()","em0008", 
          FatalErrorInArgument ,msg);
}

References components, FatalErrorInArgument, G4Exception(), and G4VEMDataSet::SetLogEnergiesData().

Field Documentation

algorithm

G4VDataSetAlgorithm* G4ShellEMDataSet::algorithm

private

Definition at line 107 of file G4ShellEMDataSet.hh.

Referenced by G4ShellEMDataSet(), GetAlgorithm(), LoadData(), LoadNonLogData(), and ~G4ShellEMDataSet().

components

std::vector<G4VEMDataSet*> G4ShellEMDataSet::components

private

Definition at line 106 of file G4ShellEMDataSet.hh.

Referenced by AddComponent(), CleanUpComponents(), FindValue(), GetComponent(), NumberOfComponents(), SetEnergiesData(), and SetLogEnergiesData().

unitData

G4double G4ShellEMDataSet::unitData

private

Definition at line 110 of file G4ShellEMDataSet.hh.

Referenced by GetUnitData(), LoadData(), LoadNonLogData(), and SaveData().

unitEnergies

G4double G4ShellEMDataSet::unitEnergies

private

Definition at line 109 of file G4ShellEMDataSet.hh.

Referenced by GetUnitEnergies(), LoadData(), LoadNonLogData(), and SaveData().

z

G4int G4ShellEMDataSet::z

private

Definition at line 111 of file G4ShellEMDataSet.hh.

Referenced by mcscore.MCVertex::dump_vertex(), FullFileName(), and mcscore.MCVertex::printout().

---

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

• source/processes/electromagnetic/lowenergy/include/G4ShellEMDataSet.hh
• source/processes/electromagnetic/lowenergy/src/G4ShellEMDataSet.cc