92 "Should not be called. Stepper must do all the work." );
110 distChord=
GetRadHelix()*(1-std::cos(0.5*Ang_curve));
112 else if(Ang_curve<
twopi)
G4double GetAngCurve() const
void AdvanceHelix(const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
void DumbStepper(const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
G4double GetRadHelix() const
virtual G4int IntegratorOrder() const
static constexpr double twopi
G4ExactHelixStepper(G4Mag_EqRhs *EqRhs)
void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
G4ThreeVector fBfieldValue
G4double DistChord() const
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
static constexpr double pi
void MagFieldEvaluate(const G4double y[], G4ThreeVector &Bfield)