CurvilinearTrajectoryParameters.h

Go to the documentation of this file.

#ifndef _TRACKER_CURVILINEARTRAJECTORYPARAMETERS_H_
#define _TRACKER_CURVILINEARTRAJECTORYPARAMETERS_H_
 
#include "DataFormats/TrajectoryState/interface/TrackCharge.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
 
#include <cmath>
 
class CurvilinearTrajectoryParameters {
public:
 
  CurvilinearTrajectoryParameters() {}
 
  CurvilinearTrajectoryParameters(const AlgebraicVector5& v, bool charged = true)
      : theQbp(charged ? v[0] : 0), thelambda(v[1]), thephi(v[2]), thexT(v[3]), theyT(v[4]) {}
 
  CurvilinearTrajectoryParameters(double aQbp, double alambda, double aphi, double axT, double ayT, bool charged = true)
      : theQbp(charged ? aQbp : 0), thelambda(alambda), thephi(aphi), thexT(axT), theyT(ayT) {}
 
  CurvilinearTrajectoryParameters(const GlobalPoint& aX, const GlobalVector& aP, TrackCharge aCharge);
 
  TrackCharge charge() const { return (0 == Qbp()) ? 0 : (Qbp() > 0 ? 1 : -1); }
 
  double signedInverseMomentum() const { return Qbp(); }
 
  AlgebraicVector5 vector() const { return AlgebraicVector5(theQbp, thelambda, thephi, thexT, theyT); }
 
  double Qbp() const { return theQbp; }
  double lambda() const { return thelambda; }
  double phi() const { return thephi; }
  double xT() const { return thexT; }
  double yT() const { return theyT; }
 
  bool updateP(double dP);
 
private:
  double theQbp;
  double thelambda;
  double thephi;
  double thexT;
  double theyT;
};
 
#endif