G4hZiegler1977Nuclear Class Reference

#include <G4hZiegler1977Nuclear.hh>

Inheritance diagram for G4hZiegler1977Nuclear:

Public Member Functions
	G4hZiegler1977Nuclear ()
	~G4hZiegler1977Nuclear ()
G4double	NuclearStoppingPower (G4double kineticEnergy, G4double z1, G4double z2, G4double m1, G4double m2) const
Public Member Functions inherited from G4VhNuclearStoppingPower
	G4VhNuclearStoppingPower ()
virtual	~G4VhNuclearStoppingPower ()
void	SetNuclearStoppingFluctuationsOn ()
void	SetNuclearStoppingFluctuationsOff ()
virtual G4double	NuclearStoppingPower (G4double kineticEnergy, G4double z1, G4double z2, G4double m1, G4double m2) const =0

Additional Inherited Members
Protected Attributes inherited from G4VhNuclearStoppingPower
G4bool	lossFlucFlag

Detailed Description

Definition at line 57 of file G4hZiegler1977Nuclear.hh.

Constructor & Destructor Documentation

G4hZiegler1977Nuclear()

G4hZiegler1977Nuclear::G4hZiegler1977Nuclear

(

)

Definition at line 62 of file G4hZiegler1977Nuclear.cc.

                                            :G4VhNuclearStoppingPower()
{;}

~G4hZiegler1977Nuclear()

G4hZiegler1977Nuclear::~G4hZiegler1977Nuclear

(

)

Definition at line 67 of file G4hZiegler1977Nuclear.cc.

{;}

Member Function Documentation

NuclearStoppingPower()

G4double G4hZiegler1977Nuclear::NuclearStoppingPower ( G4double  kineticEnergy, G4double  z1, G4double  z2, G4double  m1, G4double  m2  ) const

virtual

Implements G4VhNuclearStoppingPower.

Definition at line 72 of file G4hZiegler1977Nuclear.cc.

{  
  G4double energy = kineticEnergy/keV ;  // energy in keV
  G4double ionloss ;
  
  G4double rm = (m1 + m2Local) * std::sqrt( std::pow(z1, 0.667) + std::pow(z2, 0.667) ) ;
  
  G4double er = 32.53 * m2Local * energy / ( z1 * z2 * rm ) ;  // reduced energy
  
  if ( er < 0.01 ) {
    ionloss = std::sqrt(er) * 1.593 ; 
    
  } else if ( er < 10.0 ) {
    ionloss = 1.7 * std::sqrt(er) * std::log(er + std::exp(1.0)) / 
      (1.0 + 6.8 * er + 3.4 * std::pow(er, 1.5)) ; 
    
  } else {
    ionloss = std::log(0.47 * er) * 0.5 / er  ;
  }
 
  // Stragling
  if(lossFlucFlag) {
    G4double sig = 4.0 * m1 * m2Local * std::sqrt( (std::pow(z1, 0.23) + std::pow(z2, 0.23)) / 
                                         (std::pow(z1, 0.667) + std::pow(z2, 0.667)) ) 
                 / ((m1 +m2Local)*(m1 + m2Local)*
                    (4.0 + 0.197*std::pow(er,-1.6991)+6.584*std::pow(er,-1.0494))) ;
 
    ionloss *= G4RandGauss::shoot(1.0,sig) ;
  }
  
  ionloss *= 8.462 * z1 * z2 * m1 / rm ; // Return to [ev/(10^15 atoms/cm^2]
 
  if ( ionloss < 0.0) ionloss = 0.0 ;
  
  return ionloss;
}

---

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

• G4hZiegler1977Nuclear.hh
• G4hZiegler1977Nuclear.cc