G4BremsstrahlungCrossSectionHandler Class Reference

#include <G4BremsstrahlungCrossSectionHandler.hh>

Inheritance diagram for G4BremsstrahlungCrossSectionHandler:

Collaboration diagram for G4BremsstrahlungCrossSectionHandler:

Public Member Functions
	G4BremsstrahlungCrossSectionHandler (const G4VEnergySpectrum *spectrum, G4VDataSetAlgorithm *interpolation)
	~G4BremsstrahlungCrossSectionHandler ()
G4double	GetCrossSectionAboveThresholdForElement (G4double energy, G4double cutEnergy, G4int Z)
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=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
G4VEMDataSet *	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	LoadNonLogData (const G4String &dataFile)
void	LoadShellData (const G4String &dataFile)
void	PrintData () const
void	Clear ()

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

Detailed Description

Definition at line 62 of file G4BremsstrahlungCrossSectionHandler.hh.

Constructor & Destructor Documentation

G4BremsstrahlungCrossSectionHandler::G4BremsstrahlungCrossSectionHandler	(	const G4VEnergySpectrum *	spectrum,
		G4VDataSetAlgorithm *	interpolation
	)

Definition at line 74 of file G4BremsstrahlungCrossSectionHandler.cc.

  : theBR(spec)
{
  interp = new G4SemiLogInterpolation();
}

G4BremsstrahlungCrossSectionHandler::~G4BremsstrahlungCrossSectionHandler

(

)

Definition at line 83 of file G4BremsstrahlungCrossSectionHandler.cc.

{
  delete interp;
}

Member Function Documentation

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

protectedvirtual

Implements G4VCrossSectionHandler.

Definition at line 91 of file G4BremsstrahlungCrossSectionHandler.cc.

{
  std::vector<G4VEMDataSet*>* set = new std::vector<G4VEMDataSet*>;

  G4DataVector* energies;
  G4DataVector* cs;

  G4DataVector* log_energies;
  G4DataVector* log_cs;

  G4int nOfBins = energyVector.size();

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

  for (size_t 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->GetVecNbOfAtomsPerVolume();
    G4int nElements = material->GetNumberOfElements();

    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());

      energies = new G4DataVector;
      cs       = new G4DataVector;

      log_energies = new G4DataVector;
      log_cs       = new G4DataVector;

      G4double density = nAtomsPerVolume[i];

      for (G4int 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 value = 0.0;

        if(e > tcut) {
          G4double elemCs = FindValue(Z, e);

          value  = theBR->Probability(Z, tcut, e, e);

          value *= elemCs*density;
    }
        cs->push_back(value);

        if (value==0.) value=1e-300;
        log_cs->push_back(std::log10(value));
      }
      G4VDataSetAlgorithm* algol = interp->Clone();

      //G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,algol,1.,1.);

      G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,log_energies,log_cs,algol,1.,1.);

      setForMat->AddComponent(elSet);
    }
    set->push_back(setForMat);
  }

  return set;
}

Here is the call graph for this function:

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

Definition at line 169 of file G4BremsstrahlungCrossSectionHandler.cc.

{
  G4double value = 0.;
  if(energy > cutEnergy)
    {
      G4double elemCs = FindValue(Z, energy);
      value  = theBR->Probability(Z,cutEnergy, energy, energy);   
      value *= elemCs;
    }
  return value;
}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/electromagnetic/lowenergy/include/G4BremsstrahlungCrossSectionHandler.hh
• source/geant4.10.03.p03/source/processes/electromagnetic/lowenergy/src/G4BremsstrahlungCrossSectionHandler.cc