G4hSRIM2000p Class Reference

#include <G4hSRIM2000p.hh>

Inheritance diagram for G4hSRIM2000p:

Collaboration diagram for G4hSRIM2000p:

Public Member Functions
	G4hSRIM2000p ()
	~G4hSRIM2000p ()
G4bool	HasMaterial (const G4Material *material)
G4double	StoppingPower (const G4Material *material, G4double kineticEnergy)
G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const
Public Member Functions inherited from G4VhElectronicStoppingPower
	G4VhElectronicStoppingPower ()
virtual	~G4VhElectronicStoppingPower ()

Additional Inherited Members
Protected Member Functions inherited from G4VhElectronicStoppingPower
G4double	HeEffChargeSquare (const G4double z, const G4double kineticEnergyHe) const
G4double	GetHeMassAMU () const

Detailed Description

Definition at line 58 of file G4hSRIM2000p.hh.

Constructor & Destructor Documentation

G4hSRIM2000p::G4hSRIM2000p

(

)

Definition at line 156 of file G4hSRIM2000p.cc.

                          :G4VhElectronicStoppingPower(), 
  protonMassAMU(1.007276)
{;}

G4hSRIM2000p::~G4hSRIM2000p

(

)

Definition at line 162 of file G4hSRIM2000p.cc.

{;}

Member Function Documentation

G4double G4hSRIM2000p::ElectronicStoppingPower ( G4double  z, G4double  kineticEnergy  ) const

virtual

Implements G4VhElectronicStoppingPower.

Definition at line 191 of file G4hSRIM2000p.cc.

{
  G4double ionloss ;
  G4int i = G4int(z) - 1 ;  // index of atom
  if(i < 0)  i = 0 ;
  if(i > 91) i = 91 ;
  
  // Proton kinetic energy for parametrisation in Ziegler's units (keV/amu)  
  
  G4double T = kineticEnergy/(keV*protonMassAMU) ; 
  
  
  
  G4double e = T ;
  if ( T < 25.0 ) e = 25.0 ;
  
  // universal approximation  
  G4double slow  = a[i][0] * std::pow(e, a[i][1]) + a[i][2] * std::pow(e, a[i][3])  ;
  G4double shigh = std::log( a[i][6]/e + a[i][7]*e ) * a[i][4] / std::pow(e, a[i][5]) ;
  ionloss = slow*shigh / (slow + shigh) ; 
    
  // low energy region
  if ( T < 25.0 ) {
    
    G4double  sLocal = 0.45 ;
    // light elements
    if(6.5 > z) sLocal = 0.25 ;
    // semiconductors
    if(5 == i || 13 == i || 31 == i) sLocal = 0.375 ;
    
    ionloss *= std::pow(T/25.0, sLocal) ;
  }
  
  if ( ionloss < 0.0) ionloss = 0.0 ;
  
  return ionloss;
}

Here is the caller graph for this function:

G4bool G4hSRIM2000p::HasMaterial ( const G4Material *  material )

virtual

Implements G4VhElectronicStoppingPower.

Definition at line 167 of file G4hSRIM2000p.cc.

{
  if(1 == (material->GetNumberOfElements())) return true;
  return false ;
}

Here is the call graph for this function:

G4double G4hSRIM2000p::StoppingPower ( const G4Material *  material, G4double  kineticEnergy  )

virtual

Implements G4VhElectronicStoppingPower.

Definition at line 175 of file G4hSRIM2000p.cc.

{
  G4double ionloss = 0.0 ;

  // pure material (normally not the case for this function)
  if(1 == (material->GetNumberOfElements())) {
    G4double z = material->GetZ() ;
    ionloss = ElectronicStoppingPower( z, kineticEnergy ) ;  
  }
  
  return ionloss;
}  

Here is the call graph for this function:

---

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

• source/geant4.10.03.p03/source/processes/electromagnetic/lowenergy/include/G4hSRIM2000p.hh
• source/geant4.10.03.p03/source/processes/electromagnetic/lowenergy/src/G4hSRIM2000p.cc