G4RDeIonisationParameters Class Reference

#include <G4RDeIonisationParameters.hh>

Public Member Functions
	G4RDeIonisationParameters (G4int minZ=1, G4int maxZ=99)
	~G4RDeIonisationParameters ()
G4double	Parameter (G4int Z, G4int shellIndex, G4int parameterIndex, G4double e) const
G4double	Excitation (G4int Z, G4double e) const
void	PrintData () const

Detailed Description

Definition at line 62 of file G4RDeIonisationParameters.hh.

Constructor & Destructor Documentation

G4RDeIonisationParameters::G4RDeIonisationParameters	(	G4int	minZ = 1,
		G4int	maxZ = 99
	)

Definition at line 62 of file G4RDeIonisationParameters.cc.

  : zMin(minZ), zMax(maxZ),
  length(24)
{
  LoadData();
}

G4RDeIonisationParameters::~G4RDeIonisationParameters

(

)

Definition at line 70 of file G4RDeIonisationParameters.cc.

{ 
  // Reset the map of data sets: remove the data sets from the map 
  std::map<G4int,G4RDVEMDataSet*,std::less<G4int> >::iterator pos;

  for (pos = param.begin(); pos != param.end(); ++pos)
    {
      G4RDVEMDataSet* dataSet = (*pos).second;
      delete dataSet;
    }

  for (pos = excit.begin(); pos != excit.end(); ++pos)
    {
      G4RDVEMDataSet* dataSet = (*pos).second;
      delete dataSet;
    }

  activeZ.clear();
}

Member Function Documentation

G4double G4RDeIonisationParameters::Excitation	(	G4int	Z,
		G4double	e
	)		const

Definition at line 124 of file G4RDeIonisationParameters.cc.

References G4RDVEMDataSet::FindValue(), G4cout, G4endl, G4RDVEMDataSet::GetEnergies(), G4INCL::Math::max(), and G4INCL::Math::min().

Referenced by G4RDeIonisationSpectrum::Excitation().

{
  G4double value = 0.;
  std::map<G4int,G4RDVEMDataSet*,std::less<G4int> >::const_iterator pos;

  pos = excit.find(Z);
  if (pos!= excit.end()) {
    G4RDVEMDataSet* dataSet = (*pos).second;

    const G4DataVector ener = dataSet->GetEnergies(0);
    G4double ee = std::max(ener.front(),std::min(ener.back(),e));
    value = dataSet->FindValue(ee);
  } else {
    G4cout << "WARNING: G4IonisationParameters::Excitation "
           << "did not find ID = "
           << Z << G4endl;
  }

  return value;
}

G4double G4RDeIonisationParameters::Parameter	(	G4int	Z,
		G4int	shellIndex,
		G4int	parameterIndex,
		G4double	e
	)		const

Definition at line 91 of file G4RDeIonisationParameters.cc.

References G4RDVEMDataSet::FindValue(), G4cout, G4endl, G4RDVEMDataSet::GetComponent(), G4RDVEMDataSet::GetEnergies(), G4INCL::Math::max(), G4INCL::Math::min(), and G4RDVEMDataSet::NumberOfComponents().

Referenced by G4RDeIonisationSpectrum::AverageEnergy(), G4RDeIonisationSpectrum::Probability(), and G4RDeIonisationSpectrum::SampleEnergy().

{
  G4double value = 0.;
  G4int id = Z*100 + parameterIndex;
  std::map<G4int,G4RDVEMDataSet*,std::less<G4int> >::const_iterator pos;

  pos = param.find(id);
  if (pos!= param.end()) {
    G4RDVEMDataSet* dataSet = (*pos).second;
    G4int nShells = dataSet->NumberOfComponents();

    if(shellIndex < nShells) { 
      const G4RDVEMDataSet* component = dataSet->GetComponent(shellIndex);
      const G4DataVector ener = component->GetEnergies(0);
      G4double ee = std::max(ener.front(),std::min(ener.back(),e));
      value = component->FindValue(ee);
    } else {
      G4cout << "WARNING: G4IonisationParameters::FindParameter "
             << "has no parameters for shell= " << shellIndex 
             << "; Z= " << Z
             << G4endl;
    }
  } else {
    G4cout << "WARNING: G4IonisationParameters::Parameter "
           << "did not find ID = "
           << shellIndex << G4endl;
  }

  return value;
}

void G4RDeIonisationParameters::PrintData

(

void

)

const

Definition at line 373 of file G4RDeIonisationParameters.cc.

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

Referenced by G4RDeIonisationSpectrum::PrintData().

{
  G4cout << G4endl;
  G4cout << "===== G4RDeIonisationParameters =====" << G4endl;
  G4cout << G4endl;

  size_t nZ = activeZ.size();
  std::map<G4int,G4RDVEMDataSet*,std::less<G4int> >::const_iterator pos;

  for (size_t i=0; i<nZ; i++) {
    G4int Z = (G4int)activeZ[i];      

    for (size_t j=0; j<length; j++) {
    
      G4int index = Z*100 + j;

      pos = param.find(index);
      if (pos!= param.end()) {
        G4RDVEMDataSet* dataSet = (*pos).second;
        size_t nShells = dataSet->NumberOfComponents();

        for (size_t k=0; k<nShells; k++) {

          G4cout << "===== Z= " << Z << " shell= " << k 
                 << " parameter[" << j << "]  =====" 
                 << G4endl;
          const G4RDVEMDataSet* comp = dataSet->GetComponent(k);
          comp->PrintData();
    }
      }
    }
  }
  G4cout << "====================================" << G4endl;
}

---

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

• G4RDeIonisationParameters.hh
• G4RDeIonisationParameters.cc