PerigeeConversions.h

Go to the documentation of this file.

#ifndef PerigeeConversions_H
#define PerigeeConversions_H

#include "DataFormats/Math/interface/AlgebraicROOTObjects.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryParameters.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryError.h"

class TrajectoryStateClosestToPoint;

namespace PerigeeConversions {
  typedef FreeTrajectoryState       FTS;
  PerigeeTrajectoryParameters ftsToPerigeeParameters(const FTS& originalFTS,
    const GlobalPoint& referencePoint, double& pt);

  PerigeeTrajectoryError ftsToPerigeeError(const FTS& originalFTS);


  GlobalPoint positionFromPerigee(const PerigeeTrajectoryParameters& parameters,
    const GlobalPoint& referencePoint);

   GlobalVector momentumFromPerigee(const AlgebraicVector3& momentum, 
    const TrackCharge& charge, const GlobalPoint& referencePoint,
    const MagneticField* field);

  GlobalVector momentumFromPerigee (const PerigeeTrajectoryParameters& parameters,
                    double pt,
                    const GlobalPoint& referencePoint);


  CurvilinearTrajectoryError curvilinearError(const PerigeeTrajectoryError& perigeeError,
    const GlobalTrajectoryParameters& gtp);


  TrajectoryStateClosestToPoint trajectoryStateClosestToPoint
    (const AlgebraicVector3& momentum, const GlobalPoint& referencePoint,
     const TrackCharge& charge, const AlgebraicSymMatrix66& theCovarianceMatrix,
     const MagneticField* field);


  AlgebraicMatrix66  jacobianParameters2Cartesian
    (const AlgebraicVector3& momentum, const GlobalPoint& position,
     const TrackCharge& charge, const MagneticField* field);


  AlgebraicMatrix55 jacobianCurvilinear2Perigee(const FreeTrajectoryState& fts);

  AlgebraicMatrix55 jacobianPerigee2Curvilinear(const GlobalTrajectoryParameters& gtp);


}

#endif