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 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, float  ch  )

inline

Definition at line 17 of file CartesianLorentzForce.h.

References RKDerivative< T, N >::operator()(), and z.

                                                                        : 
    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 33 of file CartesianLorentzForce.h.

References constexpr, Basic3DVector< T >::cross(), gen::k, CartesianStateAdaptor::momentum(), CartesianStateAdaptor::position(), CartesianStateAdaptor::rkstate(), command_line::start, AlignmentTrackSelector_cfi::theCharge, 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 25 of file CartesianLorentzForce.h.

const RKLocalFieldProvider& CartesianLorentzForce::theField

private

Definition at line 24 of file CartesianLorentzForce.h.