G4hZiegler1985p Class Reference

#include <G4hZiegler1985p.hh>

Inheritance diagram for G4hZiegler1985p:

Public Member Functions
	G4hZiegler1985p ()
	~G4hZiegler1985p ()
G4bool	HasMaterial (const G4Material *material)
G4double	StoppingPower (const G4Material *material, G4double kineticEnergy)
G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const

---

Detailed Description

Definition at line 59 of file G4hZiegler1985p.hh.

---

Constructor & Destructor Documentation

G4hZiegler1985p::G4hZiegler1985p

(

)

Definition at line 168 of file G4hZiegler1985p.cc.

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

G4hZiegler1985p::~G4hZiegler1985p

(

)

Definition at line 174 of file G4hZiegler1985p.cc.

{;}

---

Member Function Documentation

G4double G4hZiegler1985p::ElectronicStoppingPower	(	G4double	z,
		G4double	kineticEnergy
	)			const [virtual]

Implements G4VhElectronicStoppingPower.

Definition at line 203 of file G4hZiegler1985p.cc.

Referenced by StoppingPower().

{
  G4double ionloss ;
  G4int i = G4int(z) - 1 ;  // index of atom
  if(i < 0)  i = 0 ;
  if(i > 91) i = 91 ;
  
  // The data and the fit from: 
  // J.F.Ziegler, J.P.Biersack, U.Littmark The Stoping and
  // Range of Ions in Solids, Vol.1, Pergamon Press, 1985
  // 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;
}

G4bool G4hZiegler1985p::HasMaterial

(

const G4Material *

material

)

[virtual]

Implements G4VhElectronicStoppingPower.

Definition at line 179 of file G4hZiegler1985p.cc.

References G4Material::GetNumberOfElements().

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

G4double G4hZiegler1985p::StoppingPower	(	const G4Material *	material,
		G4double	kineticEnergy
	)			[virtual]

Implements G4VhElectronicStoppingPower.

Definition at line 187 of file G4hZiegler1985p.cc.

References ElectronicStoppingPower(), G4Material::GetNumberOfElements(), and G4Material::GetZ().

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

---

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

• G4hZiegler1985p.hh
• G4hZiegler1985p.cc

---