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, double ch)
virtual Vector	operator() (Scalar z, const Vector &state) const
Public Member Functions inherited from RKDerivative< double, 6 >
virtual	~RKDerivative ()

Private Attributes
double	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 13 of file CartesianLorentzForce.h.

typedef Base::Scalar CartesianLorentzForce::Scalar

Definition at line 14 of file CartesianLorentzForce.h.

typedef Base::Vector CartesianLorentzForce::Vector

Definition at line 15 of file CartesianLorentzForce.h.

Constructor & Destructor Documentation

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

inline

Definition at line 17 of file CartesianLorentzForce.h.

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

Member Function Documentation

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

inlinevirtual

Derivative d(pos)/ds is simply normalized momentum

Lorentz force in absence of electric field

Implements RKDerivative< double, 6 >.

Definition at line 33 of file CartesianLorentzForce.h.

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

{
    // derivatives in case S is the free parameter
    CartesianStateAdaptor start(state);
    RKLocalFieldProvider::Vector bfield = theField.inTesla( RKLocalFieldProvider::LocalPoint(start.position()));
    double k = 2.99792458e-3; // conversion to [cm]

    CartesianStateAdaptor::Vector3D dpos = start.momentum().unit();

    CartesianStateAdaptor::Vector3D dmom = (k*theCharge) * dpos.cross( bfield);

    return CartesianStateAdaptor::rkstate( dpos, dmom);
}

Member Data Documentation

double CartesianLorentzForce::theCharge

private

Definition at line 25 of file CartesianLorentzForce.h.

Referenced by operator()().

const RKLocalFieldProvider& CartesianLorentzForce::theField

private

Definition at line 24 of file CartesianLorentzForce.h.

Referenced by operator()().