/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryParameters.h

Go to the documentation of this file.

#ifndef PerigeeTrajectoryParameters_H
#define PerigeeTrajectoryParameters_H

#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "DataFormats/Math/interface/Vector.h"
#include "DataFormats/TrackReco/interface/TrackBase.h"

class PerigeeTrajectoryParameters
{

public:

  PerigeeTrajectoryParameters() {}

  explicit PerigeeTrajectoryParameters(const AlgebraicVector &aVector, bool charged = true):
       theCurv(aVector[0]), theTheta(aVector[1]), thePhi(aVector[2]),
       theTip(aVector[3]), theLip(aVector[4]), theVector(asSVector<5>(aVector)), 
       vectorIsAvailable(true)
  {
    if ( charged )
      theCharge = theCurv>0 ? -1 : 1;
    else
      theCharge = 0;
  }

  explicit PerigeeTrajectoryParameters(const AlgebraicVector5 &aVector, bool charged = true):
       theCurv(aVector[0]), theTheta(aVector[1]), thePhi(aVector[2]),
       theTip(aVector[3]), theLip(aVector[4]), theVector(aVector), 
       vectorIsAvailable(true)
  {
    if ( charged )
      theCharge = theCurv>0 ? -1 : 1;
    else
      theCharge = 0;
  }

  PerigeeTrajectoryParameters(double aCurv, double aTheta, double aPhi,
                              double aTip, double aLip, bool charged = true):
    theCurv(aCurv), theTheta(aTheta), thePhi(aPhi), theTip(aTip), theLip(aLip),
    vectorIsAvailable(false)
  {
    if ( charged )
      theCharge = theCurv>0 ? -1 : 1;
    else
      theCharge = 0;
  }

  TrackCharge charge() const {return theCharge;}

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

  double theta() const {return theTheta;}

  double phi() const {return thePhi;}

  double transverseImpactParameter() const {return theTip;}

  double longitudinalImpactParameter() const {return theLip;}

  const AlgebraicVector  vector_old() const { return asHepVector(theVector); }

  const AlgebraicVector5 & vector() const
  {
    if (!vectorIsAvailable) {
      //theVector = AlgebraicVector5();
      theVector[0] = theCurv;
      theVector[1] = theTheta;
      theVector[2] = thePhi;
      theVector[3] = theTip;
      theVector[4] = theLip;
      vectorIsAvailable = true;
    }
    return theVector;
  }


private:
  double theCurv, theTheta, thePhi, theTip, theLip;
  TrackCharge theCharge;
  mutable AlgebraicVector5 theVector;
  mutable bool vectorIsAvailable;

};
#endif