#include <G4CrossSectionHandler.hh>

Inheritance diagram for G4CrossSectionHandler:

Public Member Functions
	G4CrossSectionHandler ()

	~G4CrossSectionHandler ()

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=0) 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 56 of file G4CrossSectionHandler.hh.

Constructor & Destructor Documentation

◆ G4CrossSectionHandler()

G4CrossSectionHandler::G4CrossSectionHandler ( )

explicit

Definition at line 65 of file G4CrossSectionHandler.cc.

66{ }

◆ ~G4CrossSectionHandler()

G4CrossSectionHandler::~G4CrossSectionHandler ( )

Definition at line 70 of file G4CrossSectionHandler.cc.

71{ }

Member Function Documentation

◆ BuildCrossSectionsForMaterials()

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

overrideprotectedvirtual

Implements G4VCrossSectionHandler.

Definition at line 76 of file G4CrossSectionHandler.cc.

{
  G4DataVector* energies;
  G4DataVector* data;
 
  G4DataVector* log_energies;
  G4DataVector* log_data;
 
  std::vector<G4VEMDataSet*>* matCrossSections = new std::vector<G4VEMDataSet*>;
 
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  G4int numOfCouples = (G4int)theCoupleTable->GetTableSize();
 
  std::size_t nOfBins = energyVector.size();
  const G4VDataSetAlgorithm* interpolationAlgo = CreateInterpolation();
 
  for (G4int mLocal=0; mLocal<numOfCouples; ++mLocal)
    {
      const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(mLocal);
      const G4Material* material= couple->GetMaterial();
      G4int nElements = (G4int)material->GetNumberOfElements();
      const G4ElementVector* elementVector = material->GetElementVector();
      const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
 
      G4VDataSetAlgorithm* algo = interpolationAlgo->Clone();
 
      G4VEMDataSet* setForMat = new G4CompositeEMDataSet(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;
 
        log_energies = new G4DataVector;
        log_data = new G4DataVector;
 
        for (std::size_t bin=0; bin<nOfBins; ++bin)
      {
        G4double e = energyVector[bin];
        energies->push_back(e);
            if (e==0.) e=1e-300;
            log_energies->push_back(std::log10(e));
        G4double cross = density*FindValue(Z,e);
        data->push_back(cross);
            if (cross==0.) cross=1e-300;
            log_data->push_back(std::log10(cross));
      }
 
        G4VDataSetAlgorithm* algo1 = interpolationAlgo->Clone();
        G4VEMDataSet* elSet = new G4EMDataSet(i,energies,data,log_energies,log_data,algo1,1.,1.);
 
        setForMat->AddComponent(elSet);
      }
 
      matCrossSections->push_back(setForMat);
    }
  delete interpolationAlgo;
  return matCrossSections;
}

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