G4StopElementSelector Class Reference

#include <G4StopElementSelector.hh>

Public Member Functions
	G4StopElementSelector ()
	~G4StopElementSelector ()
G4Element *	GetElement (const G4Material *aMaterial)
G4double	GetMuonCaptureRate (G4double Z, G4double A)
G4double	GetMuonDecayRate (G4double Z, G4double A)

Detailed Description

Definition at line 54 of file G4StopElementSelector.hh.

Constructor & Destructor Documentation

G4StopElementSelector::G4StopElementSelector

(

)

Definition at line 53 of file G4StopElementSelector.cc.

{ }

G4StopElementSelector::~G4StopElementSelector

(

)

Definition at line 59 of file G4StopElementSelector.cc.

{ }

Member Function Documentation

G4Element * G4StopElementSelector::GetElement

(

const G4Material *

aMaterial

)

Definition at line 65 of file G4StopElementSelector.cc.

References G4UniformRand, G4Material::GetAtomicNumDensityVector(), G4Material::GetElementVector(), and G4Material::GetNumberOfElements().

Referenced by G4MuonMinusCaptureAtRest::AtRestDoIt().

{
  // Fermi-Teller Z-low of mu- capture and exceptions 
  // for halogens and oxigen.
  // N.C.Mukhopadhyay Phys. Rep. 30 (1977) 1.
  G4int i;
  G4double Z;
  const G4int numberOfElements = aMaterial->GetNumberOfElements();
  const G4ElementVector* theElementVector = aMaterial->GetElementVector();

  if(1 == numberOfElements) return (*theElementVector)[0];
    
  const G4double* theAtomicNumberDensity = aMaterial->GetAtomicNumDensityVector();

  G4double sum = 0.0;
  for ( i=0; i < numberOfElements; i++ ) {

    Z = (*theElementVector)[i]->GetZ();

      // Halogens
    if( (9.0 == Z) || (17.0 == Z) || (35.0 == Z) || (53.0 == Z) || (85.0 == Z) ) {
      sum += 0.66 * Z * theAtomicNumberDensity[i] ; 

      // Oxigen
    } else if( 8.0 == Z ) {
      sum += 0.56 * Z * theAtomicNumberDensity[i] ; 

      // Others
    } else {
      sum +=        Z * theAtomicNumberDensity[i] ; 
    }
  }

  G4double random = G4UniformRand() * sum;
  sum = 0.0 ;
  i   = -1;

  // Selection of element
  do {
    i++;
    Z = (*theElementVector)[i]->GetZ();

      // Galogens
    if( (9.0 == Z) || (17.0 == Z) || (35.0 == Z) || (53.0 == Z) || (85.0 == Z) ) {
      sum += 0.66 * Z * theAtomicNumberDensity[i] ; 

      // Oxigen
    } else if( 8.0 == Z ) {
      sum += 0.56 * Z * theAtomicNumberDensity[i] ; 

      // Others
    } else {
      sum +=        Z * theAtomicNumberDensity[i] ; 
    }
  } while ( (sum < random) && (i < numberOfElements - 1) );

  return (*theElementVector)[i];
}

G4double G4StopElementSelector::GetMuonCaptureRate	(	G4double	Z,
		G4double	A
	)

Definition at line 126 of file G4StopElementSelector.cc.

References G4MuonMinusBoundDecay::GetMuonCaptureRate().

Referenced by G4MuonMinusCaptureAtRest::AtRestDoIt().

{
  return G4MuonMinusBoundDecay::GetMuonCaptureRate(G4int(Z),G4int(A));
}

G4double G4StopElementSelector::GetMuonDecayRate	(	G4double	Z,
		G4double	A
	)

Definition at line 133 of file G4StopElementSelector.cc.

References G4MuonMinusBoundDecay::GetMuonDecayRate().

Referenced by G4MuonMinusCaptureAtRest::AtRestDoIt().

{
  return G4MuonMinusBoundDecay::GetMuonDecayRate(G4int(Z));
}

---

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

• G4StopElementSelector.hh
• G4StopElementSelector.cc