G4ElementSelector Class Reference

#include <G4ElementSelector.hh>

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

Detailed Description

Definition at line 61 of file G4ElementSelector.hh.

Constructor & Destructor Documentation

G4ElementSelector::G4ElementSelector

(

)

Definition at line 49 of file G4ElementSelector.cc.

{}

G4ElementSelector::~G4ElementSelector ( )

virtual

Definition at line 54 of file G4ElementSelector.cc.

{}

Member Function Documentation

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

virtual

Definition at line 60 of file G4ElementSelector.cc.

References G4lrint(), G4UniformRand, G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), G4Element::GetIsotope(), G4Element::GetIsotopeVector(), G4Track::GetMaterial(), G4Isotope::GetN(), G4Material::GetNumberOfElements(), G4Element::GetRelativeAbundanceVector(), G4Element::GetZ(), and G4Nucleus::SetParameters().

Referenced by G4HadronStoppingProcess::AtRestDoIt().

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

---

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

• G4ElementSelector.hh
• G4ElementSelector.cc