58 G4cout <<
" G4RToEConvForPositron::G4RToEConvForPositron() ";
80 const G4double cbr1=0.02, cbr2=-5.7e-5, cbr3=1., cbr4=0.072;
84 if( std::fabs(AtomicNumber-
Z)>0.1 ) {
101 -(6.*taul+1.5*tsq-taul*(1.-tsq/3.)/t2-tsq*(0.5-tsq/12.)/
103 dEdx = (std::log(2.*taul+4.)-2.*
ionpotlog+f)/beta2;
104 dEdx = twopi_mc2_rcl2*
Z*dEdx;
105 G4double clow = dEdx*std::sqrt(taul);
106 dEdx = clow/std::sqrt(KineticEnergy/
Mass);
114 - (6.*tau+1.5*tsq-tau*(1.-tsq/3.)/t2-tsq*(0.5-tsq/12.)/
116 dEdx = (std::log(2.*tau+4.)-2.*
ionpotlog+f)/beta2;
117 dEdx = twopi_mc2_rcl2*
Z*dEdx;
121 *(cbr3+cbr4*std::log(KineticEnergy/Thigh));
122 cbrem =
Z*(
Z+1.)*cbrem*tau/beta2;
124 dEdx += twopi_mc2_rcl2*cbrem;
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
G4int GetVerboseLevel() const
virtual ~G4RToEConvForPositron()
G4GLOB_DLL std::ostream G4cout
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
static constexpr double GeV
static constexpr double MeV
virtual G4double ComputeLoss(G4double AtomicNumber, G4double KineticEnergy)
static constexpr double keV
const G4ParticleDefinition * theParticle