Chi2ChargeMeasurementEstimator.h

Go to the documentation of this file.

#ifndef CommonDet_Chi2ChargeMeasurementEstimator_H
#define CommonDet_Chi2ChargeMeasurementEstimator_H

#include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
#include "DataFormats/SiStripDetId/interface/SiStripDetId.h"
#include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
class SiStripCluster;
class TrackingRecHit;
class Chi2ChargeMeasurementEstimator GCC11_FINAL : public Chi2MeasurementEstimator {
public:

  explicit Chi2ChargeMeasurementEstimator(double maxChi2, double nSigma,
    bool cutOnPixelCharge, bool cutOnStripCharge, double minGoodPixelCharge, double minGoodStripCharge,
    float pTChargeCutThreshold = 100000.) : 
    Chi2MeasurementEstimator( maxChi2, nSigma), cutOnPixelCharge_(cutOnPixelCharge),
    cutOnStripCharge_(cutOnStripCharge), minGoodPixelCharge_(minGoodPixelCharge),
    minGoodStripCharge_(minGoodStripCharge) {
      if (pTChargeCutThreshold>=0.) pTChargeCutThreshold_=pTChargeCutThreshold;
      else pTChargeCutThreshold_=100000;
    }

  using Chi2MeasurementEstimator::estimate;
  virtual std::pair<bool,double> estimate(const TrajectoryStateOnSurface&,
                     const TrackingRecHit&) const;


  virtual Chi2ChargeMeasurementEstimator* clone() const {
    return new Chi2ChargeMeasurementEstimator(*this);
  }
private:
  bool cutOnPixelCharge_;
  bool cutOnStripCharge_;
  double minGoodPixelCharge_; 
  double minGoodStripCharge_;
  float pTChargeCutThreshold_;
  inline double minGoodCharge(int subdet) const {return (subdet>2?minGoodStripCharge_:minGoodPixelCharge_);}

  bool thickSensors (const SiStripDetId& detid) const;
  
  float sensorThickness (const DetId& detid) const;
  bool checkClusterCharge(const OmniClusterRef::ClusterStripRef cluster, float chargeCut) const;
  bool checkCharge(const TrackingRecHit& aRecHit, int subdet, float chargeCut) const;

};

#endif