#include <G4hZiegler1985p.hh>

Inheritance diagram for G4hZiegler1985p:

Public Member Functions
G4double	ElectronicStoppingPower (G4double z, G4double kineticEnergy) const override

	G4hZiegler1985p ()

G4bool	HasMaterial (const G4Material *material) override

G4double	StoppingPower (const G4Material *material, G4double kineticEnergy) override

	~G4hZiegler1985p ()

Protected Member Functions
G4double	GetHeMassAMU () const

G4double	HeEffChargeSquare (const G4double z, const G4double kineticEnergyHe) const

Private Attributes
const G4double	protonMassAMU

const G4double	theHeMassAMU

Static Private Attributes
static const G4double	a [92][8]

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.

175{;}

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

References a, keV, and protonMassAMU.

Referenced by StoppingPower().

◆ GetHeMassAMU()

G4double G4VhElectronicStoppingPower::GetHeMassAMU ( ) const

inlineprotectedinherited

Definition at line 86 of file G4VhElectronicStoppingPower.hh.

86{return theHeMassAMU;};

G4VhElectronicStoppingPower::theHeMassAMU

const G4double theHeMassAMU

Definition: G4VhElectronicStoppingPower.hh:86

References G4VhElectronicStoppingPower::theHeMassAMU.

Referenced by G4VhElectronicStoppingPower::HeEffChargeSquare().

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

References eplot::material.

◆ HeEffChargeSquare()

G4double G4VhElectronicStoppingPower::HeEffChargeSquare	(	const G4double	z,
		const G4double	kineticEnergyHe
	)		const

protectedinherited

This method returns He effective charge square parametrised according to J.F.Ziegler, J.P. Biersack, U. Littmark The Stopping and Range of Ions in Matter, Vol.1, Pergamon Press, 1985

Definition at line 69 of file G4VhElectronicStoppingPower.cc.

{
  // The aproximation of He effective charge from:
  // J.F.Ziegler, J.P. Biersack, U. Littmark
  // The Stopping and Range of Ions in Matter,
  // Vol.1, Pergamon Press, 1985
 
  static const G4double c[6] = {0.2865,  0.1266, -0.001429,
                0.02402,-0.01135, 0.001475} ;
 
  G4double e = std::log( std::max( 1.0, kineticEnergyHe/(keV*GetHeMassAMU()))) ;
  G4double x = c[0] ;
  G4double y = 1.0 ;
  for (G4int i=1; i<6; i++) {
    y *= e ;
    x += y * c[i] ;
  }
 
  G4double w = 7.6 -  e ;
  w = 1.0 + (0.007 + 0.00005*z) * G4Exp( -w*w ) ;
  w = 4.0 * (1.0 - G4Exp(-x)) * w * w ;
  return w;
}

References G4Exp(), G4VhElectronicStoppingPower::GetHeMassAMU(), keV, and G4INCL::Math::max().

Referenced by G4hICRU49He::ElectronicStoppingPower(), and G4hICRU49He::StoppingPower().

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