PerigeeConversions.h

Go to the documentation of this file.

#ifndef PerigeeConversions_H
#define PerigeeConversions_H

#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryParameters.h"
#include "TrackingTools/TrajectoryParametrization/interface/PerigeeTrajectoryError.h"

class TrajectoryStateClosestToPoint;

class PerigeeConversions
{
  typedef FreeTrajectoryState       FTS;

public:

  PerigeeTrajectoryParameters ftsToPerigeeParameters(const FTS& originalFTS,
    const GlobalPoint& referencePoint, double& pt) const;

  PerigeeTrajectoryError ftsToPerigeeError(const FTS& originalFTS) const;

//   PerigeeTrajectoryParameters helixToPerigeeParameters
//     (const reco::helix::Parameters & helixPar, const GlobalPoint& referencePoint) const;
// 
//   PerigeeTrajectoryError helixToPerigeeError(const reco::helix::Parameters & helixPar, 
//  const reco::helix::Covariance & helixCov) const;



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

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

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

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

  TrackCharge chargeFromPerigee(const PerigeeTrajectoryParameters& perigee) const;

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


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

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


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




  AlgebraicMatrix jacobianCurvilinear2Perigee_old(const FreeTrajectoryState& fts) const;

  AlgebraicMatrix jacobianPerigee2Curvilinear_old(const GlobalTrajectoryParameters& gtp) const;

  AlgebraicMatrix55 jacobianCurvilinear2Perigee(const FreeTrajectoryState& fts) const;

  AlgebraicMatrix55 jacobianPerigee2Curvilinear(const GlobalTrajectoryParameters& gtp) const;


//   AlgebraicMatrix jacobianHelix2Perigee(const reco::helix::Parameters & helixPar, 
//  const reco::helix::Covariance & helixCov) const;


};
#endif