HelixForwardPlaneCrossing.h

Go to the documentation of this file.

#ifndef HelixForwardPlaneCrossing_H_
#define HelixForwardPlaneCrossing_H_

#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "TrackingTools/GeomPropagators/interface/HelixPlaneCrossing.h"
#include "DataFormats/GeometrySurface/interface/Plane.h"
#include "FWCore/Utilities/interface/Likely.h"
#include <limits>

class HelixForwardPlaneCrossing final : public HelixPlaneCrossing {
public:
  HelixForwardPlaneCrossing(const PositionType& point,
                const DirectionType& direction,
                const float curvature,
                const PropagationDirection propDir = alongMomentum);
  // destructor
  ~HelixForwardPlaneCrossing() override {}

  std::pair<bool,double> pathLength(const Plane& plane) override {
    //
    // Protect against p_z=0 and calculate path length
    //
    if UNLIKELY( std::abs(theCosTheta)<std::numeric_limits<float>::min()  )  return std::pair<bool,double>(false,0);
    
    double dS = (plane.position().z()-theZ0) / theCosTheta;
    
    // negative logic to avoid checking for anyDirection...
    return std::make_pair( ! ( 
                  ( (thePropDir==alongMomentum) & (dS<0.) ) |
                  ( (thePropDir==oppositeToMomentum) & (dS>0.) )  
                   ) , dS);

  }

  PositionType position(double s) const override;

  DirectionType direction(double s) const override;

private:
  //
  // double precision vectors for internal use
  //
  typedef Basic3DVector<double>  PositionTypeDouble;
  typedef Basic3DVector<double>  DirectionTypeDouble;

  const double theX0,theY0,theZ0;
  double theCosPhi0,theSinPhi0;
  double theCosTheta,theSinTheta;
  const double theRho;

  const PropagationDirection thePropDir;

  mutable double theCachedS;
  mutable double theCachedDPhi;
  mutable double theCachedSDPhi;
  mutable double theCachedCDPhi;

  static const float theNumericalPrecision;
};

#endif