FieldStepper Class Reference

#include <FieldStepper.h>

Inheritance diagram for FieldStepper:

Public Member Functions
G4double	DistChord () const override
	FieldStepper (G4Mag_UsualEqRhs *eq, double del, const std::string &name)
G4int	IntegratorOrder () const override
void	Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[]) override
	~FieldStepper () override

Private Member Functions
void	selectStepper (const std::string &)

Private Attributes
double	theDelta
G4Mag_UsualEqRhs *	theEquation
G4MagIntegratorStepper *	theStepper

Detailed Description

Definition at line 8 of file FieldStepper.h.

Constructor & Destructor Documentation

FieldStepper::FieldStepper ( G4Mag_UsualEqRhs *  eq, double  del, const std::string &  name  )

explicit

Definition at line 19 of file FieldStepper.cc.

References selectStepper().

  : G4MagIntegratorStepper(eq, 6), theEquation(eq), theDelta(del) 
{
  selectStepper(nam);
}

FieldStepper::~FieldStepper ( )

override

Definition at line 25 of file FieldStepper.cc.

{}

Member Function Documentation

G4double FieldStepper::DistChord ( ) const

override

Definition at line 34 of file FieldStepper.cc.

References theStepper.

{ 
  return theStepper->DistChord(); 
}

G4int FieldStepper::IntegratorOrder ( ) const

override

Definition at line 39 of file FieldStepper.cc.

References theStepper.

{ 
  return theStepper->IntegratorOrder(); 
}

void FieldStepper::selectStepper ( const std::string &  ss )

private

Definition at line 44 of file FieldStepper.cc.

References theDelta, theEquation, and theStepper.

Referenced by FieldStepper().

{    
  if      (ss == "G4ClassicalRK4")       theStepper = new G4ClassicalRK4(theEquation);
  else if (ss == "G4NystromRK4")         theStepper = new G4NystromRK4(theEquation, theDelta);
  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::LogVerbatim("SimG4CoreMagneticField") 
    << "### FieldStepper: <" << ss << ">";
}

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

override

Definition at line 28 of file FieldStepper.cc.

References theStepper.

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

Member Data Documentation

double FieldStepper::theDelta

private

Definition at line 25 of file FieldStepper.h.

Referenced by selectStepper().

G4Mag_UsualEqRhs* FieldStepper::theEquation

private

Definition at line 24 of file FieldStepper.h.

Referenced by selectStepper().

G4MagIntegratorStepper* FieldStepper::theStepper

private

Definition at line 23 of file FieldStepper.h.

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