Geant4: G4hZiegler1985Nuclear Class Reference

G4hZiegler1985Nuclear::G4hZiegler1985Nuclear ( )

Definition at line 63 of file G4hZiegler1985Nuclear.cc.

00063                                             :G4VhNuclearStoppingPower()
00064 {;}

G4hZiegler1985Nuclear::~G4hZiegler1985Nuclear ( )

Definition at line 68 of file G4hZiegler1985Nuclear.cc.

00069 {;}

G4double G4hZiegler1985Nuclear::NuclearStoppingPower	(	G4double	kineticEnergy,
		G4double	z1,
		G4double	z2,
		G4double	m1,
		G4double	m2
	)			const `[virtual]`

Implements G4VhNuclearStoppingPower.

Definition at line 73 of file G4hZiegler1985Nuclear.cc.

References G4VhNuclearStoppingPower::lossFlucFlag.

00076 {  
00077   G4double energy = kineticEnergy/keV ;  // energy in keV
00078   G4double ionloss ;
00079   
00080   G4double rm = (m1 + m2Local) * ( std::pow(z1, .23) + std::pow(z2, .23))  ;
00081 
00082   G4double er = 32.536 * m2Local * energy / ( z1 * z2 * rm ) ;  // reduced energy
00083 
00084   if ( er <= 30 ) {
00085     ionloss = 0.5*std::log(1+1.1383*er)/
00086              (er+0.01312*std::pow(er,0.21226)+0.19593*std::sqrt(er)) ;
00087 
00088   } else {
00089     ionloss = 0.5*std::log(er)/er ;
00090   }
00091   
00092   // Stragling
00093   if(lossFlucFlag) {
00094     G4double sig = 4.0 * m1 * m2Local / ((m1 + m2Local)*(m1 + m2Local)*
00095                   (4.0 + 0.197*std::pow(er,-1.6991)+6.584*std::pow(er,-1.0494))) ;
00096 
00097 
00098     ionloss *= G4RandGauss::shoot(1.0,sig) ;
00099   }
00100 
00101   ionloss *= 8.462 * z1 * z2 * m1 / rm ; // Return to [ev/(10^15 atoms/cm^2]
00102 
00103   if ( ionloss < 0.0) ionloss = 0.0 ;
00104   
00105   return ionloss;
00106 }


Public Member Functions
	G4hZiegler1985Nuclear ()
	~G4hZiegler1985Nuclear ()
G4double	NuclearStoppingPower (G4double kineticEnergy, G4double z1, G4double z2, G4double m1, G4double m2) const

G4hZiegler1985Nuclear Class Reference

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation