G4RDShellEMDataSet Class Reference

#include <G4RDShellEMDataSet.hh>

Inheritance diagram for G4RDShellEMDataSet:

Public Member Functions
	G4RDShellEMDataSet (G4int Z, G4RDVDataSetAlgorithm *algo, G4double eUnit=CLHEP::MeV, G4double dataUnit=CLHEP::barn)
virtual	~G4RDShellEMDataSet ()
virtual G4double	FindValue (G4double energy, G4int componentId=0) const
virtual void	PrintData (void) const
virtual const G4RDVEMDataSet *	GetComponent (G4int componentId) const
virtual void	AddComponent (G4RDVEMDataSet *dataSet)
virtual size_t	NumberOfComponents (void) const
virtual const G4DataVector &	GetEnergies (G4int componentId) const
virtual const G4DataVector &	GetData (G4int componentId) const
virtual void	SetEnergiesData (G4DataVector *energies, G4DataVector *data, G4int componentId)
virtual G4bool	LoadData (const G4String &fileName)
virtual G4bool	SaveData (const G4String &fileName) const
virtual G4double	RandomSelect (G4int) const
Public Member Functions inherited from G4RDVEMDataSet
	G4RDVEMDataSet ()
virtual	~G4RDVEMDataSet ()

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

Detailed Description

Definition at line 59 of file G4RDShellEMDataSet.hh.

Constructor & Destructor Documentation

G4RDShellEMDataSet::G4RDShellEMDataSet	(	G4int	Z,
		G4RDVDataSetAlgorithm *	algo,
		G4double	eUnit = CLHEP::MeV,
		G4double	dataUnit = CLHEP::barn
	)

Definition at line 48 of file G4RDShellEMDataSet.cc.

References FatalException, G4Exception(), and G4RDShellEMDataSet().

Referenced by G4RDShellEMDataSet().

  :
  z(zeta),
  algorithm(algo),
  unitEnergies(eUnit),
  unitData(dataUnit)
{
  if (algorithm == 0)
    G4Exception("G4RDShellEMDataSet::G4RDShellEMDataSet()", "InvalidSetup",
                FatalException, "Interpolation == 0!");
}

G4RDShellEMDataSet::~G4RDShellEMDataSet ( )

virtual

Definition at line 63 of file G4RDShellEMDataSet.cc.

References CleanUpComponents().

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

Member Function Documentation

virtual void G4RDShellEMDataSet::AddComponent ( G4RDVEMDataSet *  dataSet )

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 74 of file G4RDShellEMDataSet.hh.

Referenced by LoadData().

{ components.push_back(dataSet); }

void G4RDShellEMDataSet::CleanUpComponents ( void  )

protected

Definition at line 259 of file G4RDShellEMDataSet.cc.

Referenced by LoadData(), and ~G4RDShellEMDataSet().

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

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

virtual

Implements G4RDVEMDataSet.

Definition at line 70 of file G4RDShellEMDataSet.cc.

{
  // Returns the sum over the shells corresponding to e
  G4double value = 0.;

  std::vector<G4RDVEMDataSet *>::const_iterator i(components.begin());
  std::vector<G4RDVEMDataSet *>::const_iterator end(components.end());

  while (i != end)
    {
      value += (*i)->FindValue(energy);
      i++;
    }

  return value;
}

const G4RDVDataSetAlgorithm* G4RDShellEMDataSet::GetAlgorithm ( ) const

inlineprotected

Definition at line 90 of file G4RDShellEMDataSet.hh.

{ return algorithm; }

virtual const G4RDVEMDataSet* G4RDShellEMDataSet::GetComponent ( G4int  componentId ) const

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 73 of file G4RDShellEMDataSet.hh.

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

{ return components[componentId]; }

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

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 78 of file G4RDShellEMDataSet.hh.

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

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

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

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 77 of file G4RDShellEMDataSet.hh.

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

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

G4double G4RDShellEMDataSet::GetUnitData ( ) const

inlineprotected

Definition at line 89 of file G4RDShellEMDataSet.hh.

{ return unitData; }

G4double G4RDShellEMDataSet::GetUnitEnergies ( ) const

inlineprotected

Definition at line 88 of file G4RDShellEMDataSet.hh.

{ return unitEnergies; }

G4bool G4RDShellEMDataSet::LoadData ( const G4String &  fileName )

virtual

Implements G4RDVEMDataSet.

Definition at line 126 of file G4RDShellEMDataSet.cc.

References test::a, AddComponent(), CleanUpComponents(), G4RDVDataSetAlgorithm::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("G4RDShellEMDataSet::LoadData()", "DataNotFound",
                  FatalException, message);
    }

  G4DataVector* energies = 0;
  G4DataVector* data = 0;

  G4double a = 0.;
  G4int shellIndex = 0;
  bool energyColumn = true;

  do
    {
      in >> a;
  
      if (a == -1)
    {
      if (energyColumn && energies!=0)
        {
          AddComponent(new G4RDEMDataSet(shellIndex, energies, data, algorithm->Clone(), unitEnergies, unitData));
          energies = 0;
          data = 0;
        }
   
      energyColumn = (!energyColumn);
    }
      else if (a != -2)
    {
      if (energies == 0)
        {
          energies = new G4DataVector;
          data = new G4DataVector;
        }
  
      if (energyColumn)
        energies->push_back(a * unitEnergies);
      else
        data->push_back(a * unitData);

      energyColumn = (!energyColumn);
    }
    }
  while (a != -2);

  return true;
}

virtual size_t G4RDShellEMDataSet::NumberOfComponents ( void  ) const

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 75 of file G4RDShellEMDataSet.hh.

Referenced by PrintData(), and SaveData().

{ return components.size(); }

void G4RDShellEMDataSet::PrintData ( void  ) const

virtual

Implements G4RDVEMDataSet.

Definition at line 88 of file G4RDShellEMDataSet.cc.

References G4cout, G4endl, GetComponent(), n, NumberOfComponents(), and G4RDVEMDataSet::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 G4RDShellEMDataSet::RandomSelect ( G4int  ) const

inlinevirtual

Implements G4RDVEMDataSet.

Definition at line 84 of file G4RDShellEMDataSet.hh.

{ return -1.; };

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

virtual

Implements G4RDVEMDataSet.

Definition at line 186 of file G4RDShellEMDataSet.cc.

References FatalException, G4Exception(), GetComponent(), G4RDVEMDataSet::GetData(), G4RDVEMDataSet::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("G4RDEMDataSet::SaveData()", "CannotOpenFile",
                  FatalException, message);
    }
 
  const size_t n = NumberOfComponents();
  size_t k = 0;
 
  while (k < n)
    {
      const G4RDVEMDataSet* 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 G4RDShellEMDataSet::SetEnergiesData ( G4DataVector *  energies, G4DataVector *  data, G4int  componentId  )

virtual

Implements G4RDVEMDataSet.

Definition at line 106 of file G4RDShellEMDataSet.cc.

References FatalException, G4Exception(), and G4RDVEMDataSet::SetEnergiesData().

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

  std::ostringstream message;
  message << "Component " << componentId << " not found";
 
  G4Exception("G4RDShellEMDataSet::SetEnergiesData()", "DataNotFound",
              FatalException, message.str().c_str());
}

---

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

• G4RDShellEMDataSet.hh
• G4RDShellEMDataSet.cc