#include <G4RDCrossSectionHandler.hh>

Inheritance diagram for G4RDCrossSectionHandler:

Collaboration diagram for G4RDCrossSectionHandler:

Public Member Functions
	G4RDCrossSectionHandler ()

	~G4RDCrossSectionHandler ()

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
virtual std::vector< G4RDVEMDataSet * > *	BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts=0)

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

void	ActiveElements ()

virtual G4RDVDataSetAlgorithm *	CreateInterpolation ()

const G4RDVDataSetAlgorithm *	GetInterpolation () const

Private Member Functions
	G4RDCrossSectionHandler (const G4RDCrossSectionHandler &)

G4RDCrossSectionHandler &	operator= (const G4RDCrossSectionHandler &right)

Detailed Description

Definition at line 59 of file G4RDCrossSectionHandler.hh.

Constructor & Destructor Documentation

◆ G4RDCrossSectionHandler() [1/2]

G4RDCrossSectionHandler::G4RDCrossSectionHandler ( )

Definition at line 55 of file G4RDCrossSectionHandler.cc.

56 { }

◆ ~G4RDCrossSectionHandler()

G4RDCrossSectionHandler::~G4RDCrossSectionHandler ( )

Definition at line 58 of file G4RDCrossSectionHandler.cc.

59 { }

◆ G4RDCrossSectionHandler() [2/2]

G4RDCrossSectionHandler::G4RDCrossSectionHandler ( const G4RDCrossSectionHandler & )

private

Member Function Documentation

◆ BuildCrossSectionsForMaterials()

std::vector< G4RDVEMDataSet * > * G4RDCrossSectionHandler::BuildCrossSectionsForMaterials	(	const G4DataVector &	energyVector,
		const G4DataVector *	energyCuts = `0`
	)

protectedvirtual

Implements G4RDVCrossSectionHandler.

Definition at line 62 of file G4RDCrossSectionHandler.cc.

 {
   G4DataVector* energies;
   G4DataVector* data;
 
   std::vector<G4RDVEMDataSet*>* matCrossSections = new std::vector<G4RDVEMDataSet*>;
 
   const G4ProductionCutsTable* theCoupleTable=
         G4ProductionCutsTable::GetProductionCutsTable();
   size_t numOfCouples = theCoupleTable->GetTableSize();
 
   size_t nOfBins = energyVector.size();
   const G4RDVDataSetAlgorithm* interpolationAlgo = CreateInterpolation();
 
   for (size_t m=0; m<numOfCouples; m++)
     {
       const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(m);
       const G4Material* material= couple->GetMaterial();
       G4int nElements = material->GetNumberOfElements();
       const G4ElementVector* elementVector = material->GetElementVector();
       const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
 
       G4RDVDataSetAlgorithm* algo = interpolationAlgo->Clone();
 
       G4RDVEMDataSet* setForMat = new G4RDCompositeEMDataSet(algo,1.,1.);
 
       for (G4int i=0; i<nElements; i++) {
  
         G4int Z = (G4int) (*elementVector)[i]->GetZ();
         G4double density = nAtomsPerVolume[i];
 
         energies = new G4DataVector;
         data = new G4DataVector;
 
 
         for (size_t bin=0; bin<nOfBins; bin++)
       {
         G4double e = energyVector[bin];
         energies->push_back(e);
         G4double cross = density*FindValue(Z,e);
         data->push_back(cross);
       }
 
         G4RDVDataSetAlgorithm* algo1 = interpolationAlgo->Clone();
         G4RDVEMDataSet* elSet = new G4RDEMDataSet(i,energies,data,algo1,1.,1.);
         setForMat->AddComponent(elSet);
       }
 
       matCrossSections->push_back(setForMat);
     }
   return matCrossSections;
 }

Here is the call graph for this function:

◆ operator=()

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

private

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

Public Member Functions

Protected Member Functions

Private Member Functions