Derivative calculation for the 6D cartesian case. More...

#include <CartesianLorentzForce.h>

Inheritance diagram for CartesianLorentzForce:

Public Types
typedef RKDerivative< double, 6 >	Base

typedef Base::Scalar	Scalar

typedef Base::Vector	Vector

Public Types inherited from RKDerivative< double, 6 >
typedef double	Scalar

typedef RKSmallVector< double, N >	Vector

Public Member Functions
	CartesianLorentzForce (const RKLocalFieldProvider &field, float ch)

Vector	operator() (Scalar z, const Vector &state) const override

Public Member Functions inherited from RKDerivative< double, 6 >
virtual Vector	operator() (Scalar startPar, const Vector &startState) const=0

virtual	~RKDerivative ()

Private Attributes
float	theCharge

const RKLocalFieldProvider &	theField

Detailed Description

Derivative calculation for the 6D cartesian case.

Definition at line 10 of file CartesianLorentzForce.h.

Member Typedef Documentation

◆ Base

typedef RKDerivative<double, 6> CartesianLorentzForce::Base

Definition at line 12 of file CartesianLorentzForce.h.

◆ Scalar

typedef Base::Scalar CartesianLorentzForce::Scalar

Definition at line 13 of file CartesianLorentzForce.h.

◆ Vector

typedef Base::Vector CartesianLorentzForce::Vector

Definition at line 14 of file CartesianLorentzForce.h.

Constructor & Destructor Documentation

◆ CartesianLorentzForce()

CartesianLorentzForce::CartesianLorentzForce	(	const RKLocalFieldProvider &	field,
		float	ch
	)

inline

Definition at line 16 of file CartesianLorentzForce.h.

16 : theField(field), theCharge(ch) {}

Member Function Documentation

◆ operator()()

CartesianLorentzForce::Vector CartesianLorentzForce::operator()	(	Scalar	z,
		const Vector &	state
	)		const

inlineoverride

Derivative d(pos)/ds is simply normalized momentum

Lorentz force in absence of electric field

Definition at line 26 of file CartesianLorentzForce.h.

                                                                                                         {
   // derivatives in case S is the free parameter
   CartesianStateAdaptor start(state);
   auto bfield = theField.inTesla(RKLocalFieldProvider::LocalPoint(start.position()));
   constexpr float k = 2.99792458e-3;  // conversion to [cm]
  
   auto dpos = start.momentum().unit();
  
   auto dmom = (k * theCharge) * dpos.cross(bfield);
  
   return CartesianStateAdaptor::rkstate(dpos, dmom);
 }

References RKLocalFieldProvider::inTesla(), dqmdumpme::k, CartesianStateAdaptor::rkstate(), command_line::start, theCharge, and theField.

Member Data Documentation

◆ theCharge

float CartesianLorentzForce::theCharge

private

Definition at line 22 of file CartesianLorentzForce.h.

Referenced by operator()().

◆ theField

const RKLocalFieldProvider& CartesianLorentzForce::theField