50 fElectroMagCof =
eplus*pcharge;
53 fMassCof = particleMass*particleMass ;
68 G4double pSquared = y[3]*y[3] + y[4]*y[4] + y[5]*y[5] ;
70 G4double Energy = std::sqrt( pSquared + fMassCof );
73 G4double pModuleInverse = 1.0/std::sqrt(pSquared) ;
78 G4double cof1 = fElectroMagCof*pModuleInverse ;
82 dydx[0] = y[3]*pModuleInverse ;
83 dydx[1] = y[4]*pModuleInverse ;
84 dydx[2] = y[5]*pModuleInverse ;
86 dydx[3] = cof1*(cof2*Field[0] - (y[4]*Field[5] - y[5]*Field[4])) ;
88 dydx[4] = cof1*(cof2*Field[1] - (y[5]*Field[3] - y[3]*Field[5])) ;
90 dydx[5] = cof1*(cof2*Field[2] - (y[3]*Field[4] - y[4]*Field[3])) ;
95 dydx[7] = inverse_velocity;
G4double GetCharge() const
void EvaluateRhsGivenB(const G4double y[], const G4double Field[], G4double dydx[]) const
void SetChargeMomentumMass(G4ChargeState particleCharge, G4double MomentumXc, G4double mass)