G4ElementSelector Class Reference

#include <G4ElementSelector.hh>

Public Member Functions
	G4ElementSelector ()
virtual const G4Element *	SelectZandA (const G4Track &track, G4Nucleus *)
virtual	~G4ElementSelector ()

Private Member Functions
	G4ElementSelector (const G4ElementSelector &)
G4ElementSelector &	operator= (const G4ElementSelector &right)

Private Attributes
std::vector< G4double >	prob

Detailed Description

Definition at line 60 of file G4ElementSelector.hh.

Constructor & Destructor Documentation

G4ElementSelector() [1/2]

G4ElementSelector::G4ElementSelector

(

)

Definition at line 48 of file G4ElementSelector.cc.

{}

~G4ElementSelector()

G4ElementSelector::~G4ElementSelector ( )

virtual

Definition at line 53 of file G4ElementSelector.cc.

{}

G4ElementSelector() [2/2]

G4ElementSelector::G4ElementSelector ( const G4ElementSelector &  )

private

Member Function Documentation

operator=()

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

private

SelectZandA()

const G4Element * G4ElementSelector::SelectZandA ( const G4Track &  track, G4Nucleus *  target  )

virtual

Definition at line 59 of file G4ElementSelector.cc.

{
  // Fermi-Teller Z-low of mu- capture and exceptions 
  // for halogens and oxigen.
  // N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1.
 
  size_t i = 0;
  const G4Material* mat = track.GetMaterial();
  size_t numberOfElements = mat->GetNumberOfElements();
  const G4ElementVector* theElementVector = mat->GetElementVector();
 
  if(1 < numberOfElements) {
    if(numberOfElements > prob.size()) { prob.resize(numberOfElements, 0.0); }
    
    const G4double* theAtomNumDensity = mat->GetAtomicNumDensityVector();
 
    G4double sum = 0.0;
    for (i=0; i < numberOfElements; ++i) {
 
      G4int Z = (*theElementVector)[i]->GetZasInt(); 
 
      // Halogens
      if( (9 == Z) || (17 == Z) || (35 == Z) || (53 == Z) || (85 == Z) ) {
    sum += 0.66 * Z * theAtomNumDensity[i]; 
 
    // Oxigen
      } else if( 8 == Z ) {
    sum += 0.56 * Z * theAtomNumDensity[i]; 
 
    // Others
      } else {
    sum += Z * theAtomNumDensity[i]; 
      }
      prob[i] = sum;
    }
 
    sum *= G4UniformRand();
    for (i=0; i < numberOfElements; ++i) {
      if(sum <= prob[i]) { break; }
    }
  }
  
  const G4Element* elm = (*theElementVector)[i];
  G4int Z = elm->GetZasInt();
 
  // select isotope
  const G4IsotopeVector* isv = elm->GetIsotopeVector();
  size_t ni = isv->size();
  i = 0;
 
  if(1 < ni) {
 
    const G4double* ab = elm->GetRelativeAbundanceVector();
    G4double y = G4UniformRand();
    for(i=0; i<ni; ++i) {
      y -= ab[i];
      if(y <= 0.0) { break; }
    }
  }
 
  G4int A = elm->GetIsotope(i)->GetN();
  target->SetParameters(A, Z);
 
  return elm;
}

References A, ab, G4UniformRand, G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), G4Element::GetIsotope(), G4Element::GetIsotopeVector(), G4Track::GetMaterial(), G4Isotope::GetN(), G4Material::GetNumberOfElements(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZasInt(), prob, G4Nucleus::SetParameters(), and Z.

Referenced by G4HadronStoppingProcess::AtRestDoIt(), and G4MuonMinusAtomicCapture::AtRestDoIt().

Field Documentation

prob

std::vector<G4double> G4ElementSelector::prob

private

Definition at line 76 of file G4ElementSelector.hh.

Referenced by SelectZandA().

---

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

• source/processes/hadronic/stopping/include/G4ElementSelector.hh
• source/processes/hadronic/stopping/src/G4ElementSelector.cc