Reference Manual

Go to the documentation of this file.

00001 
00002 #include "TrackPropagation/RungeKutta/interface/CartesianLorentzForce.h"
00003 #include "TrackPropagation/RungeKutta/interface/CartesianStateAdaptor.h"
00004 
00005 CartesianLorentzForce::Vector
00006 CartesianLorentzForce::operator()( Scalar z, const Vector& state) const
00007 {
00008     // derivatives in case S is the free parameter
00009     CartesianStateAdaptor start(state);
00010     RKLocalFieldProvider::Vector bfield = theField.inTesla( RKLocalFieldProvider::LocalPoint(start.position()));
00011     double k = 2.99792458e-3; // conversion to [cm]
00012 
00014     CartesianStateAdaptor::Vector3D dpos = start.momentum().unit();
00015 
00017     CartesianStateAdaptor::Vector3D dmom = k*theCharge * dpos.cross( bfield);
00018 
00019     return CartesianStateAdaptor::rkstate( dpos, dmom);
00020 }

/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/TrackPropagation/RungeKutta/src/CartesianLorentzForce.cc