95 G4double velocityVal = initVelocity.mag();
111 const G4int nvar = 6 ;
113 G4double yTemp[7], yIn[7], yTemp2[7];
117 for(i=0;i<nvar;i++) yIn[i]=yInput[i];
132 for(i=0;i<nvar;i++) {
133 yErr[i] = yOut[i] - yTemp2[i] ;
157 distChord=
GetRadHelix()*(1-std::cos(0.5*Ang_curve));
184 G4cout<<
"In G4HelixMixedStepper: Choosing Stepper for small steps. Choice is ";
185 switch ( StepperNumber )
void AdvanceHelix(const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
CLHEP::Hep3Vector G4ThreeVector
G4double GetRadHelix() const
virtual void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])=0
void MagFieldEvaluate(const G4double y[], G4ThreeVector &Bfield)
G4MagIntegratorStepper * SetupStepper(G4Mag_EqRhs *EqRhs, G4int StepperName)
G4GLOB_DLL std::ostream G4cout
G4HelixMixedStepper(G4Mag_EqRhs *EqRhs, G4int fStepperNumber=-1, G4double Angle_threshold=-1.0)
G4double GetInverseCurve(const G4double Momentum, const G4double Bmag)
static const double twopi
void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
void SetAngCurve(const G4double Ang)
G4MagIntegratorStepper * fRK4Stepper
void SetCurve(const G4double Curve)
void DumbStepper(const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
G4double fAngle_threshold
void SetVerbose(G4int newvalue)
G4double GetAngCurve() const
G4double DistChord() const