Integrator that wraps a G4ClassicalRK4 and below a minimum step size uses a drift. More...

#include <BDSIntegratorG4RK4MinStep.hh>

Inheritance diagram for BDSIntegratorG4RK4MinStep:

Collaboration diagram for BDSIntegratorG4RK4MinStep:

Public Member Functions
	BDSIntegratorG4RK4MinStep (G4Mag_EqRhs *eqOfMIn, G4double minimumStepSizeIn)

virtual void	Stepper (const G4double y[], const G4double dydx[], const G4double h, G4double yOut[], G4double yErr[])

Public Member Functions inherited from BDSIntegratorMag
	BDSIntegratorMag (G4Mag_EqRhs *eqOfMIn, G4int nVariablesIn)

virtual G4double	DistChord () const
	Estimate maximum distance of curved solution and chord. More...

virtual G4int	IntegratorOrder () const
	Geant4 requires that the integrator order must be supplied by the derived class. More...

BDSIntegratorMag &	operator= (const BDSIntegratorMag &)=delete
	Assignment and copy constructor not implemented nor used.

	BDSIntegratorMag (BDSIntegratorMag &)=delete
	Assignment and copy constructor not implemented nor used.

Public Member Functions inherited from BDSIntegratorDrift
void	AdvanceDriftMag (const G4double yIn[], const G4double h, G4double yOut[], G4double yErr[]) const
	Error array [6] all 0. More...

void	AdvanceDriftMag (const G4double yIn[], const G4ThreeVector &unitMomentum, const G4double h, G4double yOut[], G4double yErr[]) const

Public Member Functions inherited from BDSAuxiliaryNavigator
G4VPhysicalVolume *	LocateGlobalPointAndSetup (const G4ThreeVector &point, const G4ThreeVector *direction=nullptr, const G4bool pRelativeSearch=true, const G4bool ignoreDirection=true, G4bool useCurvilinear=true) const
	A wrapper for the underlying static navigator instance located within this class. More...

G4VPhysicalVolume *	LocateGlobalPointAndSetup (G4Step const *const step, G4bool useCurvilinear=true) const

BDSStep	ConvertToLocal (G4Step const *const step, G4bool useCurvilinear=true) const

BDSStep	ConvertToLocal (const G4ThreeVector &globalPosition, const G4ThreeVector &globalDirection, const G4double stepLength=0, const G4bool useCurvilinear=true, const G4double marginLength=1) const

BDSStep	ConvertToGlobalStep (const G4ThreeVector &localPosition, const G4ThreeVector &localDirection, const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertToLocal (const G4double globalPoint[3], const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertToLocal (const G4ThreeVector &globalPosition, const G4bool useCurvilinear=true) const
	Vector version - see notes above. More...

G4ThreeVector	ConvertToLocalNoSetup (const G4ThreeVector &globalPosition, const G4bool useCurvilinear=true) const
	Similar to above function but does NOT initialise the transforms. More...

G4ThreeVector	ConvertAxisToLocal (const G4ThreeVector &globalAxis, const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertAxisToLocal (const G4double globalPoint[3], const G4double globalAxis[3], const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertAxisToLocal (const G4ThreeVector &globalPoint, const G4ThreeVector &globalAxis, const G4bool useCurvilinear=true) const
	Vector version. More...

G4ThreeVector	ConvertAxisToGlobal (const G4ThreeVector &localAxis, const G4bool useCurvilinear=true) const

std::pair< G4ThreeVector, G4ThreeVector >	ConvertAxisToGlobal (const std::pair< G4ThreeVector, G4ThreeVector > &localAxis, const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertToGlobal (const G4ThreeVector &localPosition, const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertAxisToGlobal (const G4ThreeVector &globalPosition, const G4ThreeVector &localAxis, const G4bool useCurvilinear=true) const

G4ThreeVector	ConvertToGlobal (const G4ThreeVector &globalPosition, const G4ThreeVector &localPosition, const G4bool useCurvilinear=true) const

BDSStep	GlobalToCurvilinear (const G4double fieldArcLength, const G4ThreeVector &unitField, const G4double angle, const G4ThreeVector &position, const G4ThreeVector &unitMomentum, const G4double h, const G4bool useCurvilinearWorld, const G4double FCof, const G4double tilt=0)

BDSStep	GlobalToCurvilinear (const G4ThreeVector &position, const G4ThreeVector &unitMomentum, const G4double h, const G4bool useCurvilinearWorld)

BDSStep	CurvilinearToGlobal (const G4ThreeVector &localPosition, const G4ThreeVector &localMomentum, const G4bool useCurvilinearWorld)

BDSStep	CurvilinearToGlobal (const G4double fieldArcLength, const G4ThreeVector &unitField, const G4double angle, const G4ThreeVector &CLPosition, const G4ThreeVector &CLMomentum, const G4bool useCurvilinearWorld, const G4double FCof, const G4double tilt=0)

Private Member Functions
	BDSIntegratorG4RK4MinStep ()
	Private default constructor to enforce use of supplied constructor.

Private Attributes
G4MagIntegratorStepper *	g4integrator

G4double	minimumStepSize

Additional Inherited Members
Static Public Member Functions inherited from BDSAuxiliaryNavigator
static void	AttachWorldVolumeToNavigator (G4VPhysicalVolume *worldPVIn)
	Setup the navigator w.r.t. to a world volume - typically real world. More...

static void	AttachWorldVolumeToNavigatorCL (G4VPhysicalVolume *curvilinearWorldPVIn)

static void	RegisterCurvilinearBridgeWorld (G4VPhysicalVolume *curvilinearBridgeWorldPVIn)

static void	ResetNavigatorStates ()

Static Public Attributes inherited from BDSIntegratorMag
static G4double	thinElementLength = -1

static G4double	nominalMatrixRelativeMomCut = -1

static G4bool	currentTrackIsPrimary = false

Protected Member Functions inherited from BDSIntegratorMag
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 More...

void	SetDistChord (G4double distChordIn)
	Setter for distChord to private member. More...

Protected Attributes inherited from BDSIntegratorMag
G4Mag_EqRhs *	eqOfM

const G4int	nVariables
	Cache of the number of variables. More...

G4MagIntegratorStepper *	backupStepper

G4bool	zeroStrength

G4double	backupStepperMomLimit

Protected Attributes inherited from BDSAuxiliaryNavigator
G4AffineTransform	globalToLocal

G4AffineTransform	localToGlobal

G4AffineTransform	globalToLocalCL

G4AffineTransform	localToGlobalCL

G4bool	bridgeVolumeWasUsed

Static Protected Attributes inherited from BDSAuxiliaryNavigator
static G4Navigator *	auxNavigator = new G4Navigator()

static G4Navigator *	auxNavigatorCL = new G4Navigator()

static G4Navigator *	auxNavigatorCLB = new G4Navigator()

Detailed Description

Integrator that wraps a G4ClassicalRK4 and below a minimum step size uses a drift.

Author: Laurie Nevay

Definition at line 35 of file BDSIntegratorG4RK4MinStep.hh.

Constructor & Destructor Documentation

◆ BDSIntegratorG4RK4MinStep()

BDSIntegratorG4RK4MinStep::BDSIntegratorG4RK4MinStep	(	G4Mag_EqRhs *	eqOfMIn,
		G4double	minimumStepSizeIn
	)

Definition at line 26 of file BDSIntegratorG4RK4MinStep.cc.

◆ ~BDSIntegratorG4RK4MinStep()

BDSIntegratorG4RK4MinStep::~BDSIntegratorG4RK4MinStep ( )

virtual

Definition at line 33 of file BDSIntegratorG4RK4MinStep.cc.

Member Function Documentation

◆ Stepper()

void BDSIntegratorG4RK4MinStep::Stepper	(	const G4double	y[],
		const G4double	dydx[],
		const G4double	h,
		G4double	yOut[],
		G4double	yErr[]
	)

virtual

Check if the quadrupole has finite strength and use drift if not. If finite strength, convert to local curvilinear coordinates and check for paraxial approximation. If paraxial, use thick quadrupole matrix for transport, else use the G4ClassicalRK4 backup stepper.

Definition at line 38 of file BDSIntegratorG4RK4MinStep.cc.

References BDSIntegratorDrift::AdvanceDriftMag(), BDSIntegratorMag::eqOfM, BDS::IsFiniteStrength(), and BDSIntegratorMag::SetDistChord().

Here is the call graph for this function: