/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/TrackingTools/TrajectoryParametrization/interface/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