PerigeeTrajectoryParameters.h

Go to the documentation of this file.

#ifndef PerigeeTrajectoryParameters_H
#define PerigeeTrajectoryParameters_H

#include "DataFormats/TrajectoryState/interface/TrackCharge.h"
#include "DataFormats/Math/interface/Vector.h"
#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"

class PerigeeTrajectoryParameters {
public:
  PerigeeTrajectoryParameters() {}

  explicit PerigeeTrajectoryParameters(const AlgebraicVector5 &aVector, bool charged = true) : theVector(aVector) {
    if (charged)
      theCharge = theVector[0] > 0 ? -1 : 1;
    else
      theCharge = 0;
  }

  PerigeeTrajectoryParameters(double aCurv, double aTheta, double aPhi, double aTip, double aLip, bool charged = true) {
    theVector[0] = aCurv;
    theVector[1] = aTheta;
    theVector[2] = aPhi;
    theVector[3] = aTip;
    theVector[4] = aLip;

    if (charged)
      theCharge = aCurv > 0 ? -1 : 1;
    else
      theCharge = 0;
  }

  TrackCharge charge() const { return theCharge; }

  double transverseCurvature() const { return ((charge() != 0) ? theVector[0] : 0.); }

  double theta() const { return theVector[1]; }

  double phi() const { return theVector[2]; }

  double transverseImpactParameter() const { return theVector[3]; }

  double longitudinalImpactParameter() const { return theVector[4]; }

  const AlgebraicVector5 &vector() const { return theVector; }

private:
  AlgebraicVector5 theVector;
  TrackCharge theCharge;
};
#endif