IPTools.h

Go to the documentation of this file.

#ifndef TrackingTools_IPTools_h
#define TrackingTools_IPTools_h
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/Measurement1D.h"
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"
#include <utility>
#include "DataFormats/CLHEP/interface/Migration.h"
#include "TrackingTools/PatternTools/interface/TransverseImpactPointExtrapolator.h"
#include "RecoVertex/VertexTools/interface/VertexDistance.h"
#include "RecoVertex/VertexTools/interface/VertexDistance3D.h"
#include "RecoVertex/VertexTools/interface/VertexDistanceXY.h"
#include "RecoVertex/VertexPrimitives/interface/VertexState.h"


 
namespace IPTools
{
        std::pair<bool,Measurement1D> absoluteImpactParameter3D(const reco::TransientTrack & transientTrack, const  reco::Vertex & vertex);

        std::pair<bool,Measurement1D> absoluteTransverseImpactParameter(const reco::TransientTrack & transientTrack, const  reco::Vertex & vertex);

    std::pair<bool,Measurement1D> signedTransverseImpactParameter(const reco::TransientTrack & track,
             const GlobalVector & direction, const  reco::Vertex & vertex);

    std::pair<bool,Measurement1D> signedImpactParameter3D(const reco::TransientTrack & track,
             const GlobalVector & direction, const  reco::Vertex & vertex);

        std::pair<bool,Measurement1D> absoluteImpactParameter(const TrajectoryStateOnSurface & tsos  , const  reco::Vertex & vertex, VertexDistance & distanceComputer) ;


    inline  TrajectoryStateOnSurface transverseExtrapolate(const TrajectoryStateOnSurface & track,const GlobalPoint & vertexPosition, const MagneticField * field)
    {
         TransverseImpactPointExtrapolator extrapolator(field);
         return extrapolator.extrapolate(track, vertexPosition);
    }



    TrajectoryStateOnSurface closestApproachToJet(const TrajectoryStateOnSurface & state, const reco::Vertex & vertex,
         const GlobalVector& aJetDirection,const MagneticField * field);

    GlobalVector linearImpactParameter(const TrajectoryStateOnSurface & aTSOS, const GlobalPoint & point);

        
        std::pair<bool,Measurement1D> linearizedSignedImpactParameter3D(const TrajectoryStateOnSurface & state,
                 const GlobalVector & direction, const  reco::Vertex & vertex);
        
    inline std::pair<bool,Measurement1D> linearizedSignedImpactParameter3D(const reco::TransientTrack & transientTrack,
                 const GlobalVector & direction, const  reco::Vertex & vertex)
        {
      // extrapolate to the point of closest approach to the jet axis
      TrajectoryStateOnSurface closestToJetState = closestApproachToJet(transientTrack.impactPointState(), vertex, direction,transientTrack.field());
     return linearizedSignedImpactParameter3D(closestToJetState,direction,vertex);
    }
    
       std::pair<bool,Measurement1D> signedDecayLength3D(const TrajectoryStateOnSurface & state,
                 const GlobalVector & direction, const  reco::Vertex & vertex);

       inline std::pair<bool,Measurement1D> signedDecayLength3D(const reco::TransientTrack & transientTrack,
                 const GlobalVector & direction, const  reco::Vertex & vertex)
        {
          // extrapolate to the point of closest approach to the jet axis
          TrajectoryStateOnSurface closestToJetState = closestApproachToJet(transientTrack.impactPointState(), vertex, direction,transientTrack.field());
         return signedDecayLength3D(closestToJetState,direction,vertex);
        }
 

    std::pair<double,Measurement1D> jetTrackDistance(const reco::TransientTrack & track, const GlobalVector & direction,
                 const reco::Vertex & vertex);
}

#endif