GlobalTrajectoryParameters.h

Go to the documentation of this file.

#ifndef _TRACKER_GLOBALTRAJECTORYPARAMETERS_H_
#define _TRACKER_GLOBALTRAJECTORYPARAMETERS_H_

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include "DataFormats/TrajectoryState/interface/TrackCharge.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"

class MagneticField;

class GlobalTrajectoryParameters {
public:
// construct
  GlobalTrajectoryParameters() :
    theField(0), 
    theX(), theP(), 
    cachedCurvature_(0.0),
    theCharge(0),
    hasCurvature_(false)  {}  // we must initialize cache to non-NAN to avoid FPE

  GlobalTrajectoryParameters(const GlobalPoint& aX,
                             const GlobalVector& aP,
                             TrackCharge aCharge, 
                 const MagneticField* fieldProvider) :
    theField(fieldProvider),
    theX(aX), theP(aP),     
    cachedCurvature_(aCharge),
    theCharge(aCharge),
    hasCurvature_(false) {
  } // we must initialize cache to non-NAN to avoid FPE

  GlobalTrajectoryParameters(const GlobalPoint& aX,
                             const GlobalVector& direction,
                             float transverseCurvature, int, 
                 const MagneticField* fieldProvider);

  GlobalPoint position() const {
    return theX;
  }

  GlobalVector momentum() const {
    return theP;
  }

  GlobalVector direction() const {
    return theP.unit();
  }

  TrackCharge charge() const {
    return theCharge;
  }

  float signedInverseMomentum() const {
    return theCharge/theP.mag();
  }

  float signedInverseTransverseMomentum() const {
    return theCharge/theP.perp();
  }

  float transverseCurvature() const;

  AlgebraicVector6 vector() const {
    //AlgebraicVector6 v;
    //v[0] = theX.x();
    //v[1] = theX.y();
    //v[2] = theX.z();
    //v[3] = theP.x();
    //v[4] = theP.y();
    //v[5] = theP.z();
    return AlgebraicVector6(theX.x(),theX.y(),theX.z(),theP.x(),theP.y(),theP.z());
  }

 
  GlobalVector magneticFieldInInverseGeV( const GlobalPoint& x) const; 
  const MagneticField& magneticField() const {return *theField;}


private:
  const MagneticField* theField;
  GlobalPoint theX;
  GlobalVector theP;
  mutable float cachedCurvature_;
  signed char  theCharge;
  mutable bool hasCurvature_; 
  //mutable bool hasMagneticField_; mutable GlobalVector cachedMagneticField_; // 

};

#endif