G4hZiegler1985Nuclear Class Reference

#include <G4hZiegler1985Nuclear.hh>

Inheritance diagram for G4hZiegler1985Nuclear:

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

Protected Attributes
G4bool	lossFlucFlag

Detailed Description

Definition at line 58 of file G4hZiegler1985Nuclear.hh.

Constructor & Destructor Documentation

G4hZiegler1985Nuclear()

G4hZiegler1985Nuclear::G4hZiegler1985Nuclear ( )

explicit

Definition at line 63 of file G4hZiegler1985Nuclear.cc.

                                            :G4VhNuclearStoppingPower()
{;}

~G4hZiegler1985Nuclear()

G4hZiegler1985Nuclear::~G4hZiegler1985Nuclear

(

)

Definition at line 68 of file G4hZiegler1985Nuclear.cc.

{;}

Member Function Documentation

NuclearStoppingPower()

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

overridevirtual

Implements G4VhNuclearStoppingPower.

Definition at line 73 of file G4hZiegler1985Nuclear.cc.

{  
  G4double energy = kineticEnergy/keV ;  // energy in keV
  G4double ionloss ;
  
  G4double rm = (m1 + m2Local) * ( std::pow(z1, .23) + std::pow(z2, .23))  ;
 
  G4double er = 32.536 * m2Local * energy / ( z1 * z2 * rm ) ;  // reduced energy
 
  if ( er <= 30 ) {
    ionloss = 0.5*std::log(1+1.1383*er)/
             (er+0.01312*std::pow(er,0.21226)+0.19593*std::sqrt(er)) ;
 
  } else {
    ionloss = 0.5*std::log(er)/er ;
  }
  
  // Stragling
  if(lossFlucFlag) {
    G4double sig = 4.0 * m1 * m2Local / ((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;
}

References G4INCL::KinematicsUtils::energy(), keV, G4VhNuclearStoppingPower::lossFlucFlag, and G4INCL::DeJongSpin::shoot().

SetNuclearStoppingFluctuationsOff()

void G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOff ( )

inlineinherited

Definition at line 63 of file G4VhNuclearStoppingPower.hh.

{lossFlucFlag = false;}; 

References G4VhNuclearStoppingPower::lossFlucFlag.

Referenced by G4hNuclearStoppingModel::InitializeMe(), and G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOff().

SetNuclearStoppingFluctuationsOn()

void G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOn ( )

inlineinherited

Definition at line 62 of file G4VhNuclearStoppingPower.hh.

{lossFlucFlag = true;}; 

References G4VhNuclearStoppingPower::lossFlucFlag.

Referenced by G4hNuclearStoppingModel::SetNuclearStoppingFluctuationsOn().

Field Documentation

lossFlucFlag

G4bool G4VhNuclearStoppingPower::lossFlucFlag

protectedinherited

Definition at line 73 of file G4VhNuclearStoppingPower.hh.

Referenced by G4hICRU49Nuclear::NuclearStoppingPower(), NuclearStoppingPower(), G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOff(), and G4VhNuclearStoppingPower::SetNuclearStoppingFluctuationsOn().

---

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

• source/processes/electromagnetic/lowenergy/include/G4hZiegler1985Nuclear.hh
• source/processes/electromagnetic/lowenergy/src/G4hZiegler1985Nuclear.cc