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