CartesianLorentzForce Class Reference

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	~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

typedef RKDerivative<double, 6> CartesianLorentzForce::Base

Definition at line 12 of file CartesianLorentzForce.h.

typedef Base::Scalar CartesianLorentzForce::Scalar

Definition at line 13 of file CartesianLorentzForce.h.

typedef Base::Vector CartesianLorentzForce::Vector

Definition at line 14 of file CartesianLorentzForce.h.

Constructor & Destructor Documentation

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

inline

Definition at line 16 of file CartesianLorentzForce.h.

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

Member Function Documentation

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

inlineoverridevirtual

Derivative d(pos)/ds is simply normalized momentum

Lorentz force in absence of electric field

Implements RKDerivative< double, 6 >.

Definition at line 26 of file CartesianLorentzForce.h.

References Basic3DVector< T >::cross(), RKLocalFieldProvider::inTesla(), isotrackApplyRegressor::k, CartesianStateAdaptor::momentum(), CartesianStateAdaptor::position(), CartesianStateAdaptor::rkstate(), command_line::start, theCharge, theField, and Basic3DVector< T >::unit().

                                                                                                        {
  // 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);
}

Member Data Documentation

float CartesianLorentzForce::theCharge

private

Definition at line 22 of file CartesianLorentzForce.h.

Referenced by operator()().

const RKLocalFieldProvider& CartesianLorentzForce::theField

private

Definition at line 21 of file CartesianLorentzForce.h.

Referenced by operator()().