MeasurementEstimator.h

Go to the documentation of this file.

#ifndef Tracker_MeasurementEstimator_H
#define Tracker_MeasurementEstimator_H

#include "DataFormats/GeometryVector/interface/Vector2DBase.h"
#include "DataFormats/GeometryVector/interface/LocalTag.h"
#include <utility>

class Plane;
class TrajectoryStateOnSurface;
class Surface;
class TrackingRecHit;

class MeasurementEstimator {
public:

  struct OpaquePayload { virtual ~OpaquePayload(){} int tag=0;};

  using Local2DVector = Vector2DBase< float, LocalTag>;


  MeasurementEstimator() {}
  MeasurementEstimator(float maxSag, float minToll) :
     m_maxSagitta(maxSag),
     m_minTolerance2(minToll*minToll){}

  virtual ~MeasurementEstimator() {}

  using HitReturnType     = std::pair<bool,double>;
  using SurfaceReturnType = bool;

  virtual HitReturnType estimate( const TrajectoryStateOnSurface& ts, 
                  const TrackingRecHit& hit) const = 0;

  /* verify the compatibility of the Hit with the Trajectory based
   * on hit properties other than those used in estimate 
   * (that usually computes the compatibility of the Trajectory with the Hit)
   * 
   */
  virtual bool preFilter(const TrajectoryStateOnSurface&, OpaquePayload const &) const { return true;}


  virtual SurfaceReturnType estimate( const TrajectoryStateOnSurface& ts, 
                      const Plane& plane) const = 0;

  virtual MeasurementEstimator* clone() const = 0;

  virtual Local2DVector 
  maximalLocalDisplacement( const TrajectoryStateOnSurface& ts,
                const Plane& plane) const=0;

  float maxSagitta() const { return m_maxSagitta;}
  float minTolerance2() const { return m_minTolerance2;}


private:
  float m_maxSagitta=-1.; // maximal sagitta for linear approximation
  float m_minTolerance2=100.; // square of minimum tolerance ot be considered inside a detector

};

#endif // Tracker_MeasurementEstimator_H