Chi2MeasurementEstimatorBase.h

Go to the documentation of this file.

#ifndef CommonDet_Chi2MeasurementEstimatorBase_H
#define CommonDet_Chi2MeasurementEstimatorBase_H

#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
#include <limits>

class Chi2MeasurementEstimatorBase : public MeasurementEstimator {
public:
  explicit Chi2MeasurementEstimatorBase(double maxChi2,
                                        double nSigma = 3.,
                                        float maxDisp = std::numeric_limits<float>::max())
      : theMaxChi2(maxChi2), theNSigma(nSigma), theMaxDisplacement(maxDisp) {}

  template <typename... Args>
  Chi2MeasurementEstimatorBase(double maxChi2, double nSigma, float maxDisp, Args&&... args)
      : MeasurementEstimator(args...), theMaxChi2(maxChi2), theNSigma(nSigma), theMaxDisplacement(maxDisp) {}

  std::pair<bool, double> estimate(const TrajectoryStateOnSurface& ts, const TrackingRecHit&) const override = 0;

  bool estimate(const TrajectoryStateOnSurface& ts, const Plane& plane) const final;

  Local2DVector maximalLocalDisplacement(const TrajectoryStateOnSurface& ts, const Plane& plane) const final;

  double chiSquaredCut() const { return theMaxChi2; }
  double nSigmaCut() const { return theNSigma; }

protected:
  std::pair<bool, double> returnIt(double est) const {
    return est > chiSquaredCut() ? HitReturnType(false, est) : HitReturnType(true, est);
  }

private:
  const double theMaxChi2;
  const double theNSigma;
  const float theMaxDisplacement;
};

#endif