G4HadElementSelector Class Reference

#include <G4HadronXSDataTable.hh>

Public Member Functions
	G4HadElementSelector (G4DynamicParticle *, G4CrossSectionDataStore *, const G4Material *, G4int bins, G4double emin, G4double emax, G4bool spline)
	~G4HadElementSelector ()
void	Dump ()
const G4Element *	SelectRandomAtom (G4double e) const

Detailed Description

Definition at line 58 of file G4HadronXSDataTable.hh.

Constructor & Destructor Documentation

G4HadElementSelector()

G4HadElementSelector::G4HadElementSelector	(	G4DynamicParticle *	dp,
		G4CrossSectionDataStore *	xs,
		const G4Material *	mat,
		G4int	bins,
		G4double	emin,
		G4double	emax,
		G4bool	spline )

Definition at line 51 of file G4HadronXSDataTable.cc.

{
  std::size_t n = mat->GetNumberOfElements();
  nElmMinusOne = G4int(n - 1);
  theElementVector = mat->GetElementVector();
  if(nElmMinusOne > 0) {
    G4PhysicsVector* first = nullptr;
    xSections.resize(n, first);
    first = new G4PhysicsLogVector(emin,emax,bins,false);
    xSections[0] = first;
    for(std::size_t i=1; i<n; ++i) {
      xSections[i] = new G4PhysicsVector(*first);
    }
    std::vector<G4double> temp;
    temp.resize(n, 0.0);
    for(G4int j=0; j<=bins; ++j) {
      G4double cross = 0.0;
      G4double e = first->Energy(j);
      dp->SetKineticEnergy(e);
      for(std::size_t i=0; i<n; ++i) {
    cross += xs->GetCrossSection(dp, (*theElementVector)[i], mat);
        temp[i] = cross;
      }
      G4double fact = (cross > 0.0) ? 1.0/cross : 0.0;
      for(std::size_t i=0; i<n; ++i) {
    G4double y = (i<n-1) ? temp[i]*fact : 1.0; 
        xSections[i]->PutValue(j, y);
      }
    }
  }
}

~G4HadElementSelector()

G4HadElementSelector::~G4HadElementSelector

(

)

Definition at line 89 of file G4HadronXSDataTable.cc.

{
  if(nElmMinusOne > 0) {
    for(G4int i=0; i<=nElmMinusOne; ++i) { delete xSections[i]; }
  }
}

Member Function Documentation

Dump()

void G4HadElementSelector::Dump

(

)

Definition at line 98 of file G4HadronXSDataTable.cc.

{}

SelectRandomAtom()

const G4Element * G4HadElementSelector::SelectRandomAtom ( G4double e ) const

inline

Definition at line 70 of file G4HadronXSDataTable.hh.

  {
    const G4Element* element = (*theElementVector)[nElmMinusOne];
    if (nElmMinusOne > 0) {
      G4double x = G4UniformRand();
      for(G4int i=0; i<nElmMinusOne; ++i) {
    if (x <= xSections[i]->Value(e)) {
      element = (*theElementVector)[i];
      break;
    }
      }
    }
    return element;
  }

---

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

• G4HadronXSDataTable.hh
• G4HadronXSDataTable.cc