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 60 of file G4eIonisationParameters.hh.

Constructor & Destructor Documentation

G4eIonisationParameters()

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

Definition at line 62 of file G4eIonisationParameters.cc.

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

~G4eIonisationParameters()

G4eIonisationParameters::~G4eIonisationParameters

(

)

Definition at line 70 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

Excitation()

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

Definition at line 124 of file G4eIonisationParameters.cc.

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

Referenced by G4eIonisationSpectrum::Excitation().

Parameter()

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

Definition at line 91 of file G4eIonisationParameters.cc.

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

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

PrintData()

void G4eIonisationParameters::PrintData

(

)

const

Definition at line 383 of file G4eIonisationParameters.cc.

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

Referenced by G4eIonisationSpectrum::PrintData().

---

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

• G4eIonisationParameters.hh
• G4eIonisationParameters.cc