Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4HelixImplicitEuler Class Reference

#include <G4HelixImplicitEuler.hh>

+ Inheritance diagram for G4HelixImplicitEuler:

Public Member Functions

 G4HelixImplicitEuler (G4Mag_EqRhs *EqRhs)
 
 ~G4HelixImplicitEuler ()
 
void DumbStepper (const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
 
G4int IntegratorOrder () const
 
- Public Member Functions inherited from G4MagHelicalStepper
 G4MagHelicalStepper (G4Mag_EqRhs *EqRhs)
 
virtual ~G4MagHelicalStepper ()
 
 G4MagHelicalStepper (const G4MagHelicalStepper &)=delete
 
G4MagHelicalStepperoperator= (const G4MagHelicalStepper &)=delete
 
virtual void Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
 
virtual void DumbStepper (const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])=0
 
G4double DistChord () const
 
- Public Member Functions inherited from G4MagIntegratorStepper
 G4MagIntegratorStepper (G4EquationOfMotion *Equation, G4int numIntegrationVariables, G4int numStateVariables=12, G4bool isFSAL=false)
 
virtual ~G4MagIntegratorStepper ()=default
 
 G4MagIntegratorStepper (const G4MagIntegratorStepper &)=delete
 
G4MagIntegratorStepperoperator= (const G4MagIntegratorStepper &)=delete
 
virtual void Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])=0
 
virtual G4double DistChord () const =0
 
void NormaliseTangentVector (G4double vec[6])
 
void NormalisePolarizationVector (G4double vec[12])
 
void RightHandSide (const G4double y[], G4double dydx[]) const
 
void RightHandSide (const G4double y[], G4double dydx[], G4double field[]) const
 
G4int GetNumberOfVariables () const
 
G4int GetNumberOfStateVariables () const
 
virtual G4int IntegratorOrder () const =0
 
G4int IntegrationOrder ()
 
G4EquationOfMotionGetEquationOfMotion ()
 
const G4EquationOfMotionGetEquationOfMotion () const
 
void SetEquationOfMotion (G4EquationOfMotion *newEquation)
 
unsigned long GetfNoRHSCalls ()
 
void ResetfNORHSCalls ()
 
G4bool IsFSAL () const
 

Additional Inherited Members

- Protected Member Functions inherited from G4MagHelicalStepper
void LinearStep (const G4double yIn[], G4double h, G4double yHelix[]) const
 
void AdvanceHelix (const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
 
void MagFieldEvaluate (const G4double y[], G4ThreeVector &Bfield)
 
G4double GetInverseCurve (const G4double Momentum, const G4double Bmag)
 
void SetAngCurve (const G4double Ang)
 
G4double GetAngCurve () const
 
void SetCurve (const G4double Curve)
 
G4double GetCurve () const
 
void SetRadHelix (const G4double Rad)
 
G4double GetRadHelix () const
 
- Protected Member Functions inherited from G4MagIntegratorStepper
void SetIntegrationOrder (G4int order)
 
void SetFSAL (G4bool flag=true)
 

Detailed Description

Definition at line 45 of file G4HelixImplicitEuler.hh.

Constructor & Destructor Documentation

◆ G4HelixImplicitEuler()

G4HelixImplicitEuler::G4HelixImplicitEuler ( G4Mag_EqRhs EqRhs)

Definition at line 42 of file G4HelixImplicitEuler.cc.

◆ ~G4HelixImplicitEuler()

G4HelixImplicitEuler::~G4HelixImplicitEuler ( )

Definition at line 47 of file G4HelixImplicitEuler.cc.

48{
49}

Member Function Documentation

◆ DumbStepper()

void G4HelixImplicitEuler::DumbStepper ( const G4double  y[],
G4ThreeVector  Bfld,
G4double  h,
G4double  yout[] 
)
virtual

Implements G4MagHelicalStepper.

Definition at line 52 of file G4HelixImplicitEuler.cc.

56{
57 const G4int nvar = 6 ;
58 G4double yTemp[6], yTemp2[6];
59 G4ThreeVector Bfld_endpoint;
60
61 // Step forward like in the explicit euler case
62 //
63 AdvanceHelix( yIn, Bfld, h, yTemp);
64
65 // now obtain the new field value at the new point
66 //
67 MagFieldEvaluate(yTemp, Bfld_endpoint);
68
69 // and also advance along a helix for this field value
70 //
71 AdvanceHelix( yIn, Bfld_endpoint, h, yTemp2);
72
73 // we take the average
74 //
75 for( G4int i = 0; i < nvar; ++i )
76 {
77 yOut[i] = 0.5 * ( yTemp[i] + yTemp2[i] );
78 }
79
80 // NormaliseTangentVector( yOut );
81}
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
void AdvanceHelix(const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
void MagFieldEvaluate(const G4double y[], G4ThreeVector &Bfield)

◆ IntegratorOrder()

G4int G4HelixImplicitEuler::IntegratorOrder ( ) const
inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 60 of file G4HelixImplicitEuler.hh.

60{ return 2; }

The documentation for this class was generated from the following files: