G4eIonisationParameters Class Reference

#include <G4eIonisationParameters.hh>

Public Member Functions
	G4eIonisationParameters (G4int minZ=1, G4int maxZ=99)
	~G4eIonisationParameters ()
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 61 of file G4eIonisationParameters.hh.

Constructor & Destructor Documentation

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

Definition at line 63 of file G4eIonisationParameters.cc.

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

G4eIonisationParameters::~G4eIonisationParameters

(

)

Definition at line 71 of file G4eIonisationParameters.cc.

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

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

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

  activeZ.clear();
}

Member Function Documentation

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

Definition at line 125 of file G4eIonisationParameters.cc.

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

Referenced by G4eIonisationSpectrum::Excitation().

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

  pos = excit.find(Z);
  if (pos!= excit.end()) {
    G4VEMDataSet* 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 G4eIonisationParameters::Parameter	(	G4int	Z,
		G4int	shellIndex,
		G4int	parameterIndex,
		G4double	e
	)		const

Definition at line 92 of file G4eIonisationParameters.cc.

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

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

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

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

    if(shellIndex < nShells) { 
      const G4VEMDataSet* 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 G4eIonisationParameters::PrintData

(

void

)

const

Definition at line 384 of file G4eIonisationParameters.cc.

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

Referenced by G4eIonisationSpectrum::PrintData().

{
  G4cout << G4endl;
  G4cout << "===== G4eIonisationParameters =====" << G4endl;
  G4cout << G4endl;

  size_t nZ = activeZ.size();
  std::map<G4int,G4VEMDataSet*,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()) {
        G4VEMDataSet* 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 G4VEMDataSet* comp = dataSet->GetComponent(k);
          comp->PrintData();
    }
      }
    }
  }
  G4cout << "====================================" << G4endl;
}

---

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

• G4eIonisationParameters.hh
• G4eIonisationParameters.cc