G4ShellEMDataSet Class Reference

#include <G4ShellEMDataSet.hh>

Inheritance diagram for G4ShellEMDataSet:

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

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

Detailed Description

Definition at line 60 of file G4ShellEMDataSet.hh.

Constructor & Destructor Documentation

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

Definition at line 64 of file G4ShellEMDataSet.cc.

References FatalErrorInArgument, G4Exception(), and G4ShellEMDataSet().

Referenced by G4ShellEMDataSet().

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

G4ShellEMDataSet::~G4ShellEMDataSet ( )

virtual

Definition at line 77 of file G4ShellEMDataSet.cc.

References CleanUpComponents().

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

Member Function Documentation

virtual void G4ShellEMDataSet::AddComponent ( G4VEMDataSet *  dataSet )

inlinevirtual

Implements G4VEMDataSet.

Definition at line 75 of file G4ShellEMDataSet.hh.

Referenced by LoadData(), and LoadNonLogData().

{ components.push_back(dataSet); }

void G4ShellEMDataSet::CleanUpComponents ( void  )

protected

Definition at line 389 of file G4ShellEMDataSet.cc.

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

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

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

virtual

Implements G4VEMDataSet.

Definition at line 84 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;
}

const G4VDataSetAlgorithm* G4ShellEMDataSet::GetAlgorithm ( ) const

inlineprotected

Definition at line 101 of file G4ShellEMDataSet.hh.

{ return algorithm; }

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

inlinevirtual

Implements G4VEMDataSet.

Definition at line 74 of file G4ShellEMDataSet.hh.

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

{ return components[componentId]; }

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

inlinevirtual

Implements G4VEMDataSet.

Definition at line 79 of file G4ShellEMDataSet.hh.

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

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

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

inlinevirtual

Implements G4VEMDataSet.

Definition at line 78 of file G4ShellEMDataSet.hh.

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

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

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

inlinevirtual

Implements G4VEMDataSet.

Definition at line 81 of file G4ShellEMDataSet.hh.

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

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

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

inlinevirtual

Implements G4VEMDataSet.

Definition at line 80 of file G4ShellEMDataSet.hh.

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

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

G4double G4ShellEMDataSet::GetUnitData ( ) const

inlineprotected

Definition at line 100 of file G4ShellEMDataSet.hh.

{ return unitData; }

G4double G4ShellEMDataSet::GetUnitEnergies ( ) const

inlineprotected

Definition at line 99 of file G4ShellEMDataSet.hh.

{ return unitEnergies; }

G4bool G4ShellEMDataSet::LoadData ( const G4String &  fileName )

virtual

Implements G4VEMDataSet.

Definition at line 160 of file G4ShellEMDataSet.cc.

References test::a, AddComponent(), CleanUpComponents(), G4VDataSetAlgorithm::Clone(), FatalException, and G4Exception().

{
  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;
}

G4bool G4ShellEMDataSet::LoadNonLogData ( const G4String &  fileName )

virtual

Implements G4VEMDataSet.

Definition at line 244 of file G4ShellEMDataSet.cc.

References test::a, AddComponent(), CleanUpComponents(), G4VDataSetAlgorithm::Clone(), FatalException, and G4Exception().

{
  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;
}

virtual size_t G4ShellEMDataSet::NumberOfComponents ( void  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 76 of file G4ShellEMDataSet.hh.

Referenced by PrintData(), and SaveData().

{ return components.size(); }

void G4ShellEMDataSet::PrintData ( void  ) const

virtual

Implements G4VEMDataSet.

Definition at line 102 of file G4ShellEMDataSet.cc.

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

{
  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++;
    }
}

virtual G4double G4ShellEMDataSet::RandomSelect ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 95 of file G4ShellEMDataSet.hh.

{ return -1.; };

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

virtual

Implements G4VEMDataSet.

Definition at line 317 of file G4ShellEMDataSet.cc.

References FatalException, G4Exception(), GetComponent(), G4VEMDataSet::GetData(), G4VEMDataSet::GetEnergies(), left, n, and NumberOfComponents().

{
  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;
}

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

virtual

Implements G4VEMDataSet.

Definition at line 120 of file G4ShellEMDataSet.cc.

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

{
  G4VEMDataSet* component = components[componentId];
 
  if (component)
    {
      component->SetEnergiesData(energies, data, 0);
      return;
    }

  G4String msg = "component " + (G4String)componentId + " not found";
 
  G4Exception("G4ShellEMDataSet::SetEnergiesData()","em0008", FatalErrorInArgument ,msg);
}

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

virtual

Implements G4VEMDataSet.

Definition at line 138 of file G4ShellEMDataSet.cc.

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

{
  G4VEMDataSet* component = components[componentId];
 
  if (component)
    {
      component->SetLogEnergiesData(energies, data, log_energies, log_data, 0);
      return;
    }

  G4String msg = "component " + (G4String)componentId + " not found";
 
  G4Exception("G4ShellEMDataSet::SetLogEnergiesData()","em0008", FatalErrorInArgument ,msg);

}

---

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

• G4ShellEMDataSet.hh
• G4ShellEMDataSet.cc