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);

}  // namespace PerigeeConversions

#endif