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.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 91 of file G4RDeIonisationParameters.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

void G4RDeIonisationParameters::PrintData

(

void

)

const

Definition at line 373 of file G4RDeIonisationParameters.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/eRosita/physics/include/G4RDeIonisationParameters.hh
• source/geant4.10.03.p03/examples/advanced/eRosita/physics/src/G4RDeIonisationParameters.cc