#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.

73{ components.push_back(dataSet); }

G4ShellEMDataSet::components

std::vector< G4VEMDataSet * > components

Definition: G4ShellEMDataSet.hh:106

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.

100{ return algorithm; }

References algorithm.

◆ GetComponent()

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

inlineoverridevirtual

Implements G4VEMDataSet.

Definition at line 72 of file G4ShellEMDataSet.hh.

72{ 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.

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

G4ShellEMDataSet::GetComponent

const G4VEMDataSet * GetComponent(G4int componentId) const override

Definition: G4ShellEMDataSet.hh:72

G4VEMDataSet::GetData

virtual const G4DataVector & GetData(G4int componentId) const =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.

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

G4VEMDataSet::GetEnergies

virtual const G4DataVector & GetEnergies(G4int componentId) const =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.

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

G4VEMDataSet::GetLogData

virtual const G4DataVector & GetLogData(G4int componentId) const =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.

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

G4VEMDataSet::GetLogEnergies

virtual const G4DataVector & GetLogEnergies(G4int componentId) const =0

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

◆ GetUnitData()

G4double G4ShellEMDataSet::GetUnitData ( ) const

inlineprotected

Definition at line 99 of file G4ShellEMDataSet.hh.

99{ return unitData; }

References unitData.

◆ GetUnitEnergies()

G4double G4ShellEMDataSet::GetUnitEnergies ( ) const

inlineprotected

Definition at line 98 of file G4ShellEMDataSet.hh.

98{ 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.

74{ 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.

92{ 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

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ShellEMDataSet() [1/2]

◆ ~G4ShellEMDataSet()

◆ G4ShellEMDataSet() [2/2]

Member Function Documentation

◆ AddComponent()

◆ CleanUpComponents()

◆ FindValue()

◆ FullFileName()

◆ GetAlgorithm()

◆ GetComponent()

◆ GetData()

◆ GetEnergies()

◆ GetLogData()

◆ GetLogEnergies()

◆ GetUnitData()

◆ GetUnitEnergies()

◆ LoadData()

◆ LoadNonLogData()

◆ NumberOfComponents()

◆ operator=()

◆ PrintData()

◆ RandomSelect()

◆ SaveData()

◆ SetEnergiesData()

◆ SetLogEnergiesData()

Field Documentation

◆ algorithm

◆ components

◆ unitData

◆ unitEnergies

◆ z