/data/refman/pasoursint/CMSSW_4_1_8_patch13/src/TrackPropagation/RungeKutta/interface/CartesianStateAdaptor.h

Go to the documentation of this file.

#ifndef CartesianStateAdaptor_H
#define CartesianStateAdaptor_H

class CartesianStateAdaptor {
public:

    typedef double                            Scalar;
    typedef Basic3DVector<Scalar>             Vector3D;

    CartesianStateAdaptor( const RKSmallVector<double,6>& rk) :
      pos_(rk[0], rk[1], rk[2]), mom_(rk[3], rk[4], rk[5]) {}

    const Vector3D& position() const { return pos_;}
    const Vector3D& momentum() const { return mom_;}

    static Vector3D position(const RKSmallVector<double,6>& rk) {
      return Vector3D(rk[0], rk[1], rk[2]);
    }

    static Vector3D momentum(const RKSmallVector<double,6>& rk) {
      return Vector3D(rk[3], rk[4], rk[5]);
    }

    static RKSmallVector<double,6> rkstate( const Vector3D& pos, const Vector3D& mom) {
        RKSmallVector<double,6> res;
        res[0] = pos.x();
        res[1] = pos.y();
        res[2] = pos.z();
        res[3] = mom.x();
        res[4] = mom.y();
        res[5] = mom.z();
        return res;
    }

private:

    Vector3D pos_;
    Vector3D mom_;

};

#endif