20#include "BDSIntegratorSextupole.hh"
21#include "BDSMagnetStrength.hh"
23#include "BDSUtilities.hh"
26#include "G4Mag_EqRhs.hh"
27#include "G4ThreeVector.hh"
29#include "CLHEP/Units/SystemOfUnits.h"
35 G4Mag_EqRhs* eqOfMIn):
39 bDoublePrime = brho * (*strength)[
"k2"] / CLHEP::m3;
42 G4cout << __METHOD_NAME__ <<
"B'' = " << bDoublePrime << G4endl;
51 G4ThreeVector mom = G4ThreeVector(yIn[3], yIn[4], yIn[5]);
52 G4double momUnit = mom.mag();
55 if (std::abs(kappa) < 1e-12)
62 G4ThreeVector pos = G4ThreeVector( yIn[0], yIn[1], yIn[2]);
66 G4ThreeVector localMomUnit = localMom.unit();
68 G4double x0 = localPos.x();
69 G4double y0 = localPos.y();
71 G4double xp = localMomUnit.x();
72 G4double yp = localMomUnit.y();
73 G4double zp = localMomUnit.z();
76 const G4double halfH = 0.5*h;
77 x0 = x0 + localMomUnit.x()*halfH;
78 y0 = y0 + localMomUnit.y()*halfH;
80 G4double x02My02 = x0*x0 - y0*y0;
84 localA.setX(zp*x02My02);
85 localA.setY(-2*zp*x0*y0);
86 localA.setZ(xp*x02My02-2*yp*x0*y0);
BDSStep ConvertToLocal(G4Step const *const step, G4bool useCurvilinear=true) const
void AdvanceDriftMag(const G4double yIn[], const G4double h, G4double yOut[], G4double yErr[]) const
Error array [6] all 0.
Common functionality for Euler integrators.
void AdvanceChord(const G4double h, G4ThreeVector &localPos, G4ThreeVector &localMom, const G4ThreeVector &localA)
void SetDistChord(G4double distChordIn)
Setter for distChord to private member.
void ConvertToGlobal(const G4ThreeVector &localPos, const G4ThreeVector &localMom, G4double yOut[], G4double yErr[], const G4double momScaling=1.0)
scaling of momentum in case localMom is a unit vector
virtual void AdvanceHelix(const G4double yIn[], G4double h, G4double yOut[], G4double yErr[])
Calculate the new particle coordinates.
G4double bDoublePrime
2nd derivative of the field
Efficient storage of magnet strengths.
A simple class to represent the positions of a step.
G4ThreeVector PostStepPoint() const
Accessor.
G4ThreeVector PreStepPoint() const
Accessor.
G4bool IsFinite(G4double value, G4double tolerance=std::numeric_limits< double >::epsilon())