G4hZiegler1985p Class Reference

#include <G4hZiegler1985p.hh>

Inheritance diagram for G4hZiegler1985p:

Public Member Functions
	G4hZiegler1985p ()
	~G4hZiegler1985p ()
G4bool	HasMaterial (const G4Material *material) override
G4double	StoppingPower (const G4Material *material, G4double kineticEnergy) override
G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const override
Public Member Functions inherited from G4VhElectronicStoppingPower
	G4VhElectronicStoppingPower ()
virtual	~G4VhElectronicStoppingPower ()
	G4VhElectronicStoppingPower (G4VhElectronicStoppingPower &)=delete
G4VhElectronicStoppingPower &	operator= (const G4VhElectronicStoppingPower &right)=delete

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 59 of file G4hZiegler1985p.hh.

Constructor & Destructor Documentation

G4hZiegler1985p()

G4hZiegler1985p::G4hZiegler1985p ( )

explicit

Definition at line 168 of file G4hZiegler1985p.cc.

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

~G4hZiegler1985p()

G4hZiegler1985p::~G4hZiegler1985p

(

)

Definition at line 174 of file G4hZiegler1985p.cc.

{;}

Member Function Documentation

ElectronicStoppingPower()

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

overridevirtual

Implements G4VhElectronicStoppingPower.

Definition at line 203 of file G4hZiegler1985p.cc.

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

Referenced by StoppingPower().

HasMaterial()

G4bool G4hZiegler1985p::HasMaterial ( const G4Material * material )

overridevirtual

Implements G4VhElectronicStoppingPower.

Definition at line 179 of file G4hZiegler1985p.cc.

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

StoppingPower()

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

overridevirtual

Implements G4VhElectronicStoppingPower.

Definition at line 187 of file G4hZiegler1985p.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;
}  

---

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

• G4hZiegler1985p.hh
• G4hZiegler1985p.cc