G4RDeIonisationCrossSectionHandler Class Reference

#include <G4RDeIonisationCrossSectionHandler.hh>

Inheritance diagram for G4RDeIonisationCrossSectionHandler:

Collaboration diagram for G4RDeIonisationCrossSectionHandler:

Public Member Functions
	G4RDeIonisationCrossSectionHandler (const G4RDVEnergySpectrum *spec, G4RDVDataSetAlgorithm *alg, G4double emin, G4double emax, G4int nbin)
	~G4RDeIonisationCrossSectionHandler ()
Public Member Functions inherited from G4RDVCrossSectionHandler
	G4RDVCrossSectionHandler ()
	G4RDVCrossSectionHandler (G4RDVDataSetAlgorithm *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	~G4RDVCrossSectionHandler ()
void	Initialise (G4RDVDataSetAlgorithm *interpolation=0, 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
G4RDVEMDataSet *	BuildMeanFreePathForMaterials (const G4DataVector *energyCuts=0)
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	LoadShellData (const G4String &dataFile)
void	PrintData () const
void	Clear ()

Protected Member Functions
std::vector< G4RDVEMDataSet * > *	BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts)
Protected Member Functions inherited from G4RDVCrossSectionHandler
G4int	NumberOfComponents (G4int Z) const
void	ActiveElements ()
virtual G4RDVDataSetAlgorithm *	CreateInterpolation ()
const G4RDVDataSetAlgorithm *	GetInterpolation () const

Detailed Description

Definition at line 63 of file G4RDeIonisationCrossSectionHandler.hh.

Constructor & Destructor Documentation

G4RDeIonisationCrossSectionHandler::G4RDeIonisationCrossSectionHandler	(	const G4RDVEnergySpectrum *	spec,
		G4RDVDataSetAlgorithm *	alg,
		G4double	emin,
		G4double	emax,
		G4int	nbin
	)

Definition at line 62 of file G4RDeIonisationCrossSectionHandler.cc.

 :  G4RDVCrossSectionHandler(),
    theParam(spec)
{
  G4RDVCrossSectionHandler::Initialise(alg, emin, emax, nbin);
  interp = new G4RDLinLogLogInterpolation();
}

Here is the call graph for this function:

G4RDeIonisationCrossSectionHandler::~G4RDeIonisationCrossSectionHandler

(

)

Definition at line 73 of file G4RDeIonisationCrossSectionHandler.cc.

{
  delete interp;
}

Member Function Documentation

std::vector< G4RDVEMDataSet * > * G4RDeIonisationCrossSectionHandler::BuildCrossSectionsForMaterials ( const G4DataVector &  energyVector, const G4DataVector *  energyCuts  )

protectedvirtual

Implements G4RDVCrossSectionHandler.

Definition at line 79 of file G4RDeIonisationCrossSectionHandler.cc.

{
  G4int verbose = 0;
  std::vector<G4RDVEMDataSet*>* set = new std::vector<G4RDVEMDataSet*>;

  G4DataVector* energies;
  G4DataVector* cs;
  G4int nOfBins = energyVector.size();

  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();

  for (size_t m=0; m<numOfCouples; m++) {

    const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(m);
    const G4Material* material= couple->GetMaterial();
    const G4ElementVector* elementVector = material->GetElementVector();
    const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
    G4int nElements = material->GetNumberOfElements();

    if(verbose > 0) {
      G4cout << "eIonisation CS for " << m << "th material "
             << material->GetName()
             << "  eEl= " << nElements << G4endl;
    }

    G4double tcut  = (*energyCuts)[m];

    G4RDVDataSetAlgorithm* algo = interp->Clone();
    G4RDVEMDataSet* setForMat = new G4RDCompositeEMDataSet(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;
      G4double density = nAtomsPerVolume[i];

      for (G4int bin=0; bin<nOfBins; bin++) {

        G4double e = energyVector[bin];
        energies->push_back(e);
        G4double value = 0.0;

        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 && m == 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;
        }

      }
    }
        cs->push_back(value);
      }
      G4RDVDataSetAlgorithm* algo = interp->Clone();
      G4RDVEMDataSet* elSet = new G4RDEMDataSet(i,energies,cs,algo,1.,1.);
      setForMat->AddComponent(elSet);
    }
    set->push_back(setForMat);
  }

  return set;
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/eRosita/physics/include/G4RDeIonisationCrossSectionHandler.hh
• source/geant4.10.03.p02/examples/advanced/eRosita/physics/src/G4RDeIonisationCrossSectionHandler.cc