G4RDCrossSectionHandler Class Reference

#include <G4RDCrossSectionHandler.hh>

Inheritance 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

Detailed Description

Definition at line 59 of file G4RDCrossSectionHandler.hh.

Constructor & Destructor Documentation

G4RDCrossSectionHandler::G4RDCrossSectionHandler

(

)

Definition at line 55 of file G4RDCrossSectionHandler.cc.

{ }

G4RDCrossSectionHandler::~G4RDCrossSectionHandler

(

)

Definition at line 58 of file G4RDCrossSectionHandler.cc.

{ }

Member Function Documentation

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

protectedvirtual

Implements G4RDVCrossSectionHandler.

Definition at line 62 of file G4RDCrossSectionHandler.cc.

References G4RDVEMDataSet::AddComponent(), plottest35::bin, G4RDVDataSetAlgorithm::Clone(), G4RDVCrossSectionHandler::CreateInterpolation(), density, G4RDVCrossSectionHandler::FindValue(), G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4Material::GetNumberOfElements(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), python.hepunit::m, and eplot::material.

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

---

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

• G4RDCrossSectionHandler.hh
• G4RDCrossSectionHandler.cc