G4ElementSelector Class Reference

#include <G4ElementSelector.hh>

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

Detailed Description

Definition at line 60 of file G4ElementSelector.hh.

Constructor & Destructor Documentation

G4ElementSelector()

G4ElementSelector::G4ElementSelector

(

)

Definition at line 48 of file G4ElementSelector.cc.

{}

~G4ElementSelector()

G4ElementSelector::~G4ElementSelector ( )

virtual

Definition at line 53 of file G4ElementSelector.cc.

{}

Member Function Documentation

SelectZandA()

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 = G4lrint((*theElementVector)[i]->GetZ()); 
 
      // 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; }
    }
  }
  
  G4Element* elm = (*theElementVector)[i];
  G4int Z = G4lrint(elm->GetZ());
 
  // 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;
}

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

---

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

• G4ElementSelector.hh
• G4ElementSelector.cc