70 if( angleThreshold < 0.0 )
127 constexpr G4int nvar = 6 ;
128 constexpr G4int nvarMax = 8 ;
129 G4double yTemp[nvarMax], yIn[nvarMax], yTemp2[nvarMax];
138 for(
G4int i=0; i<nvar; ++i)
161 for(
G4int i=0; i<nvar; ++i)
163 yErr[i] = yOut[i] - yTemp2[i];
189 distChord=
GetRadHelix()*(1-std::cos(0.5*Ang_curve));
209 G4cout <<
"In HelixMixedStepper::Number of calls to smallStepStepper = "
220 switch ( StepperNumber )
319 G4cout <<
" chosen as stepper for small steps in G4HelixMixedStepper."
static constexpr double twopi
static constexpr double pi
CLHEP::Hep3Vector G4ThreeVector
G4GLOB_DLL std::ostream G4cout
G4HelixMixedStepper(G4Mag_EqRhs *EqRhs, G4int StepperNumber=-1, G4double Angle_threshold=-1.0)
void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
G4double DistChord() const
void DumbStepper(const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
G4double fAngle_threshold
G4MagIntegratorStepper * fRK4Stepper
G4MagIntegratorStepper * SetupStepper(G4Mag_EqRhs *EqRhs, G4int StepperName)
void SetCurve(const G4double Curve)
void AdvanceHelix(const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
G4double GetRadHelix() const
void MagFieldEvaluate(const G4double y[], G4ThreeVector &Bfield)
G4double GetInverseCurve(const G4double Momentum, const G4double Bmag)
void SetAngCurve(const G4double Ang)
G4double GetAngCurve() const
virtual void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])=0