#include <G4eIonisationCrossSectionHandler.hh>

Inheritance diagram for G4eIonisationCrossSectionHandler:

Public Member Functions
	G4eIonisationCrossSectionHandler (const G4VEnergySpectrum spec, G4VDataSetAlgorithm alg, G4double emin, G4double emax, G4int nbin)

	~G4eIonisationCrossSectionHandler ()

G4double	GetCrossSectionAboveThresholdForElement (G4double energy, G4double cutEnergy, G4int Z)

G4eIonisationCrossSectionHandler &	operator= (const G4eIonisationCrossSectionHandler &right)=delete

	G4eIonisationCrossSectionHandler (const G4eIonisationCrossSectionHandler &)=delete

Public Member Functions inherited from G4VCrossSectionHandler
	G4VCrossSectionHandler ()

	G4VCrossSectionHandler (G4VDataSetAlgorithm interpolation, G4double minE=250 CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int nBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)

virtual	~G4VCrossSectionHandler ()

void	Initialise (G4VDataSetAlgorithm interpolation=nullptr, G4double minE=250 CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int numberOfBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)

G4int	SelectRandomAtom (const G4MaterialCutsCouple *couple, G4double e) const

const G4Element *	SelectRandomElement (const G4MaterialCutsCouple *material, G4double e) const

G4int	SelectRandomShell (G4int Z, G4double e) const

G4VEMDataSet *	BuildMeanFreePathForMaterials (const G4DataVector *energyCuts=nullptr)

G4double	FindValue (G4int Z, G4double e) const

G4double	FindValue (G4int Z, G4double e, G4int shellIndex) const

G4double	ValueForMaterial (const G4Material *material, G4double e) const

void	LoadData (const G4String &dataFile)

void	LoadNonLogData (const G4String &dataFile)

void	LoadShellData (const G4String &dataFile)

void	PrintData () const

void	Clear ()

	G4VCrossSectionHandler (const G4VCrossSectionHandler &)=delete

G4VCrossSectionHandler &	operator= (const G4VCrossSectionHandler &right)=delete

Protected Member Functions
std::vector< G4VEMDataSet * > *	BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts) override

Protected Member Functions inherited from G4VCrossSectionHandler
G4int	NumberOfComponents (G4int Z) const

void	ActiveElements ()

virtual std::vector< G4VEMDataSet * > *	BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts=nullptr)=0

virtual G4VDataSetAlgorithm *	CreateInterpolation ()

const G4VDataSetAlgorithm *	GetInterpolation () const

Detailed Description

Definition at line 59 of file G4eIonisationCrossSectionHandler.hh.

Constructor & Destructor Documentation

◆ G4eIonisationCrossSectionHandler() [1/2]

G4eIonisationCrossSectionHandler::G4eIonisationCrossSectionHandler	(	const G4VEnergySpectrum *	spec,
		G4VDataSetAlgorithm *	alg,
		G4double	emin,
		G4double	emax,
		G4int	nbin
	)

explicit

Definition at line 72 of file G4eIonisationCrossSectionHandler.cc.

 :  G4VCrossSectionHandler(),
    theParam(spec),verbose(0)
{
  G4VCrossSectionHandler::Initialise(alg, emin, emax, nbin);
  interp = new G4LinLogLogInterpolation();
}

◆ ~G4eIonisationCrossSectionHandler()

G4eIonisationCrossSectionHandler::~G4eIonisationCrossSectionHandler ( )

Definition at line 84 of file G4eIonisationCrossSectionHandler.cc.

{
  delete interp;
}

◆ G4eIonisationCrossSectionHandler() [2/2]

G4eIonisationCrossSectionHandler::G4eIonisationCrossSectionHandler ( const G4eIonisationCrossSectionHandler & )

delete

Member Function Documentation

◆ BuildCrossSectionsForMaterials()

std::vector< G4VEMDataSet * > * G4eIonisationCrossSectionHandler::BuildCrossSectionsForMaterials	(	const G4DataVector &	energyVector,
		const G4DataVector *	energyCuts
	)

overrideprotectedvirtual

Implements G4VCrossSectionHandler.

Definition at line 91 of file G4eIonisationCrossSectionHandler.cc.

{
  std::vector<G4VEMDataSet*>* set = new std::vector<G4VEMDataSet*>;
 
  G4DataVector* energies;
  G4DataVector* cs;
 
  G4DataVector* log_energies;
  G4DataVector* log_cs;
 
  std::size_t nOfBins = energyVector.size();
 
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  G4int numOfCouples = (G4int)theCoupleTable->GetTableSize();
 
  for (G4int mLocal=0; mLocal<numOfCouples; ++mLocal) {
 
    const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(mLocal);
    const G4Material* material= couple->GetMaterial();
    const G4ElementVector* elementVector = material->GetElementVector();
    const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
    G4int nElements = (G4int)material->GetNumberOfElements();
 
    if(verbose > 0) 
      {
    G4cout << "eIonisation CS for " << mLocal << "th material "
           << material->GetName()
           << "  eEl= " << nElements << G4endl;
      }
 
    G4double tcut  = (*energyCuts)[mLocal];
 
    G4VDataSetAlgorithm* algo = interp->Clone();
    G4VEMDataSet* setForMat = new G4CompositeEMDataSet(algo,1.,1.);
 
    for (G4int i=0; i<nElements; ++i) {
 
      G4int Z = (G4int) (*elementVector)[i]->GetZ();
      G4int nShells = NumberOfComponents(Z);
 
      energies = new G4DataVector;
      cs       = new G4DataVector;
 
      log_energies = new G4DataVector;
      log_cs       = new G4DataVector;
 
      G4double density = nAtomsPerVolume[i];
 
      for (std::size_t bin=0; bin<nOfBins; ++bin) {
 
        G4double e = energyVector[bin];
        energies->push_back(e);
        log_energies->push_back(std::log10(e));
        G4double value = 0.0;
        G4double log_value = -300;
 
        if(e > tcut) {
          for (G4int n=0; n<nShells; n++) {
            G4double cross = FindValue(Z, e, n);
            G4double p = theParam->Probability(Z, tcut, e, e, n);
            value += cross * p * density;
 
        if(verbose>0 && mLocal == 0 && e>=1. && e<=0.) 
        {
              G4cout << "G4eIonCrossSH: e(MeV)= " << e/MeV
                     << " n= " << n
                     << " cross= " << cross
                     << " p= " << p
                     << " value= " << value
                     << " tcut(MeV)= " << tcut/MeV
                     << " rho= " << density
                     << " Z= " << Z
                     << G4endl;
          }
      }
          if (value == 0.) value = 1e-300;
          log_value = std::log10(value);
    }
        cs->push_back(value);
        log_cs->push_back(log_value);
      }
      G4VDataSetAlgorithm* algoLocal = interp->Clone();
      G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,log_energies,log_cs,algoLocal,1.,1.);
 
      setForMat->AddComponent(elSet);
    }
    set->push_back(setForMat);
  }
 
  return set;
}

◆ GetCrossSectionAboveThresholdForElement()

G4double G4eIonisationCrossSectionHandler::GetCrossSectionAboveThresholdForElement	(	G4double	energy,
		G4double	cutEnergy,
		G4int	Z
	)

Definition at line 188 of file G4eIonisationCrossSectionHandler.cc.

{
  G4int nShells = NumberOfComponents(Z);
  G4double value = 0.;
  if(energy > cutEnergy) 
    {
      for (G4int n=0; n<nShells; ++n) {
    G4double cross = FindValue(Z, energy, n);
    G4double p = theParam->Probability(Z, cutEnergy, energy, energy, n);
    value += cross * p;
      }
    }
  return value;
}

Referenced by G4LivermoreIonisationModel::ComputeCrossSectionPerAtom().

◆ operator=()

G4eIonisationCrossSectionHandler & G4eIonisationCrossSectionHandler::operator= ( const G4eIonisationCrossSectionHandler & right )