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<limits>

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, float mpt) :
     m_maxSagitta(maxSag),
     m_minTolerance2(minToll*minToll),
     m_minPt2ForHitRecoveryInGluedDet(mpt*mpt)
     {}

  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;}
  float minPt2ForHitRecoveryInGluedDet() const { return m_minPt2ForHitRecoveryInGluedDet;}

private:
  /*
   *  why here? 
   * MeasurementEstimator is the only configurable item that percolates down to geometry event by event (actually hit by hit) and not at initialization time
   * It is therefore the natural candidate to collect all parameters that affect pattern-recongnition 
   * and require to be controlled with higher granularity than job level (such as iteration by iteration)
   */ 
  float m_maxSagitta=-1.; // maximal sagitta for linear approximation
  float m_minTolerance2=100.; // square of minimum tolerance ot be considered inside a detector
  float m_minPt2ForHitRecoveryInGluedDet=std::numeric_limits<float>::max();  // 0.81 to mitigate wrong preAmpl setting
};

#endif // Tracker_MeasurementEstimator_H