sim::FieldStepper Class Reference

#include <FieldStepper.h>

Inheritance diagram for sim::FieldStepper:

Public Member Functions
double	DistChord () const override
	FieldStepper (G4Mag_UsualEqRhs *eq, double del=0.0)
int	IntegratorOrder () const override
G4MagIntegratorStepper *	select (const std::string &s)
void	Stepper (const double y[], const double dydx[], double h, double yout[], double yerr[]) override
	~FieldStepper () override

Private Attributes
double	delta
G4Mag_UsualEqRhs *	theEquation
G4MagIntegratorStepper *	theStepper

Detailed Description

Definition at line 9 of file FieldStepper.h.

Constructor & Destructor Documentation

FieldStepper::FieldStepper	(	G4Mag_UsualEqRhs *	eq,
		double	del = 0.0
	)

Definition at line 21 of file FieldStepper.cc.

                                                            :
  G4MagIntegratorStepper(eq, 6), theEquation(eq), delta(del) {}

FieldStepper::~FieldStepper ( )

override

Definition at line 24 of file FieldStepper.cc.

{}

Member Function Documentation

double FieldStepper::DistChord ( ) const

override

Definition at line 30 of file FieldStepper.cc.

References theStepper.

{ return theStepper->DistChord(); }

int FieldStepper::IntegratorOrder ( ) const

override

Definition at line 32 of file FieldStepper.cc.

References theStepper.

{ return theStepper->IntegratorOrder(); }

G4MagIntegratorStepper * FieldStepper::select

(

const std::string &

s

)

Definition at line 35 of file FieldStepper.cc.

References delta, theEquation, and theStepper.

Referenced by sim::FieldBuilder::configureForVolume().

{
    if      (ss == "G4ClassicalRK4")       theStepper = new G4ClassicalRK4(theEquation);
    else if (ss == "G4NystromRK4")         theStepper = new G4NystromRK4(theEquation, delta);
    else if (ss == "G4SimpleRunge")        theStepper = new G4SimpleRunge(theEquation);
    else if (ss == "G4SimpleHeum")         theStepper = new G4SimpleHeum(theEquation);
    else if (ss == "G4CashKarpRKF45")      theStepper = new G4CashKarpRKF45(theEquation);
    else if (ss == "G4DormandPrince745")   theStepper = new G4DormandPrince745(theEquation);
    else if (ss == "G4BogackiShampine45")  theStepper = new G4BogackiShampine45(theEquation);
    else if (ss == "G4TsitourasRK45")      theStepper = new G4TsitourasRK45(theEquation);
    else if (ss == "G4ImplicitEuler")      theStepper = new G4ImplicitEuler(theEquation);
    else if (ss == "G4HelixExplicitEuler") theStepper = new G4HelixExplicitEuler(theEquation);
    else if (ss == "G4HelixImplicitEuler") theStepper = new G4HelixImplicitEuler(theEquation);
    else if (ss == "G4HelixSimpleRunge")   theStepper = new G4HelixSimpleRunge(theEquation);
    else if (ss == "G4HelixHeum")          theStepper = new G4HelixHeum(theEquation);
    else
    {
      edm::LogWarning("SimG4CoreMagneticField") 
        << " FieldStepper <" << ss << "> is not known, defaulting to G4ClassicalRK4 ";
      theStepper = new G4ClassicalRK4(theEquation);
    }
    edm::LogInfo("SimG4CoreMagneticField") 
      << "### FieldStepper: <" << ss << ">";
    return theStepper;
}

void FieldStepper::Stepper ( const double  y[], const double  dydx[], double  h, double  yout[], double  yerr[]  )

override

Definition at line 26 of file FieldStepper.cc.

References theStepper.

{ theStepper->Stepper(y,dydx,h,yout,yerr); }

Member Data Documentation

double sim::FieldStepper::delta

private

Definition at line 22 of file FieldStepper.h.

Referenced by select().

G4Mag_UsualEqRhs* sim::FieldStepper::theEquation

private

Definition at line 21 of file FieldStepper.h.

Referenced by select().

G4MagIntegratorStepper* sim::FieldStepper::theStepper

private

Definition at line 20 of file FieldStepper.h.

Referenced by DistChord(), IntegratorOrder(), select(), and Stepper().