51 G4cout <<
" G4RToEConvForElectron::G4RToEConvForElectron() ";
70 static G4double taul, ionpot, ionpotlog;
71 const G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
72 const G4double Tlow=10.*keV, Thigh=1.*GeV;
78 if( std::fabs(AtomicNumber-Z)>0.1 ) {
81 ionpot = 1.6e-5*MeV*std::exp(0.9*std::log(Z))/Mass;
82 ionpotlog = std::log(ionpot);
94 G4double f = 1.-beta2+std::log(tsq/2.)
95 +(0.5+0.25*tsq+(1.+2.*taul)*std::log(0.5))/(t1*t1);
96 dEdx = (std::log(2.*taul+4.)-2.*ionpotlog+f)/beta2;
97 dEdx = twopi_mc2_rcl2*Z*dEdx;
98 G4double clow = dEdx*std::sqrt(taul);
99 dEdx = clow/std::sqrt(KineticEnergy/Mass);
106 G4double f = 1.-beta2+std::log(tsq/2.)
107 +(0.5+0.25*tsq+(1.+2.*tau)*std::log(0.5))/(t1*t1);
108 dEdx = (std::log(2.*tau+4.)-2.*ionpotlog+f)/beta2;
109 dEdx = twopi_mc2_rcl2*Z*dEdx;
113 *(cbr3+cbr4*std::log(KineticEnergy/Thigh));
114 cbrem = Z*(Z+1.)*cbrem*tau/beta2;
116 cbrem *= bremfactor ;
118 dEdx += twopi_mc2_rcl2*cbrem;
G4DLLIMPORT std::ostream G4cout
G4double GetPDGMass() const
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
static G4ParticleTable * GetParticleTable()
virtual ~G4RToEConvForElectron()
virtual G4double ComputeLoss(G4double AtomicNumber, G4double KineticEnergy) const
G4int GetVerboseLevel() const
const G4ParticleDefinition * theParticle