G4HadElementSelector Class Reference

#include <G4HadronXSDataTable.hh>

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

Private Member Functions
	G4HadElementSelector (G4HadElementSelector &)=delete
G4HadElementSelector &	operator= (const G4HadElementSelector &right)=delete

Private Attributes
G4int	nElmMinusOne
const G4ElementVector *	theElementVector
std::vector< G4PhysicsVector * >	xSections

Detailed Description

Definition at line 58 of file G4HadronXSDataTable.hh.

Constructor & Destructor Documentation

G4HadElementSelector() [1/2]

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.

{
  G4int n = mat->GetNumberOfElements();
  nElmMinusOne = 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(G4int 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(G4int 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(G4int i=0; i<n; ++i) {
    G4double y = (i<n-1) ? temp[i]*fact : 1.0; 
        xSections[i]->PutValue(j, y);
      }
    }
  }
}

References anonymous_namespace{G4HyperonSampler.cc}::bins, emax, G4PhysicsVector::Energy(), G4CrossSectionDataStore::GetCrossSection(), G4Material::GetElementVector(), G4Material::GetNumberOfElements(), CLHEP::detail::n, nElmMinusOne, G4DynamicParticle::SetKineticEnergy(), theElementVector, and xSections.

~G4HadElementSelector()

G4HadElementSelector::~G4HadElementSelector

(

)

Definition at line 89 of file G4HadronXSDataTable.cc.

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

References nElmMinusOne, and xSections.

G4HadElementSelector() [2/2]

G4HadElementSelector::G4HadElementSelector ( G4HadElementSelector &  )

privatedelete

Member Function Documentation

Dump()

void G4HadElementSelector::Dump

(

)

Definition at line 98 of file G4HadronXSDataTable.cc.

{}

operator=()

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

privatedelete

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

References G4UniformRand, nElmMinusOne, and xSections.

Field Documentation

nElmMinusOne

G4int G4HadElementSelector::nElmMinusOne

private

Definition at line 90 of file G4HadronXSDataTable.hh.

Referenced by G4HadElementSelector(), SelectRandomAtom(), and ~G4HadElementSelector().

theElementVector

const G4ElementVector* G4HadElementSelector::theElementVector

private

Definition at line 91 of file G4HadronXSDataTable.hh.

Referenced by G4HadElementSelector().

xSections

std::vector<G4PhysicsVector*> G4HadElementSelector::xSections

private

Definition at line 92 of file G4HadronXSDataTable.hh.

Referenced by G4HadElementSelector(), SelectRandomAtom(), and ~G4HadElementSelector().

---

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

• source/processes/hadronic/cross_sections/include/G4HadronXSDataTable.hh
• source/processes/hadronic/cross_sections/src/G4HadronXSDataTable.cc