PathToPlane2Order.h

Go to the documentation of this file.

#ifndef PathToPlane2Order_H
#define PathToPlane2Order_H

#include "DataFormats/GeometryVector/interface/Basic3DVector.h"
#include "TrackingTools/GeomPropagators/interface/HelixPlaneCrossing.h"
#include "TrackPropagation/RungeKutta/interface/CartesianState.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "DataFormats/GeometrySurface/interface/Plane.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"

class RKLocalFieldProvider;

class PathToPlane2Order {
public:

    typedef Plane::Scalar                                Scalar;
    typedef Basic3DVector<Scalar>                        Vector3D;
    typedef GloballyPositioned<Scalar>                   Frame;

    PathToPlane2Order( const RKLocalFieldProvider& fld, const Frame* fieldFrame) : 
      theField(fld), theFieldFrame(fieldFrame) {}

    std::pair<bool,double> operator()( const Plane& plane, 
                       const Vector3D& position,
                       const Vector3D& momentum,
                       double charge,
                       const PropagationDirection propDir = alongMomentum);

    std::pair<bool,double> operator()( const Plane& plane, 
                       const GlobalPoint& position,
                       const GlobalVector& momentum,
                       double charge,
                       const PropagationDirection propDir = alongMomentum) {
    return operator()( plane, theFieldFrame->toLocal(position).basicVector(), 
               theFieldFrame->toLocal(momentum).basicVector(), charge, propDir);
    }

private:

    const RKLocalFieldProvider& theField;
    const Frame*                theFieldFrame;
};

#endif