#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.

{
  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);
    first->SetSpline(spline);
    xSections[0] = first;
    for(G4int i=1; i<n; ++i) {
      xSections[i] = new G4PhysicsVector(*first);
    }
    std::vector<double> 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);
      }
    }
  }
}

◆ ~G4HadElementSelector()

G4HadElementSelector::~G4HadElementSelector ( )

Definition at line 90 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 99 of file G4HadronXSDataTable.cc.

100{}

◆ 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: