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