#include <Chi2Strip1DEstimator.h>
Public Member Functions | |
Chi2Strip1DEstimator (double maxChi2, double nSigma=3.) | |
virtual Chi2Strip1DEstimator * | clone () const |
virtual std::pair< bool, double > | estimate (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const |
A Chi2 MeasurementEstimator that only uses the X coordinate in the measurement frame (the one perpendicular to the strip). Ported from ORCA.
Definition at line 16 of file Chi2Strip1DEstimator.h.
Chi2Strip1DEstimator::Chi2Strip1DEstimator | ( | double | maxChi2, |
double | nSigma = 3. |
||
) | [inline, explicit] |
Definition at line 20 of file Chi2Strip1DEstimator.h.
Referenced by clone().
: Chi2MeasurementEstimatorBase( maxChi2, nSigma) {}
virtual Chi2Strip1DEstimator* Chi2Strip1DEstimator::clone | ( | void | ) | const [inline, virtual] |
Implements MeasurementEstimator.
Definition at line 26 of file Chi2Strip1DEstimator.h.
References Chi2Strip1DEstimator().
{ return new Chi2Strip1DEstimator(*this); }
pair< bool, double > Chi2Strip1DEstimator::estimate | ( | const TrajectoryStateOnSurface & | ts, |
const TransientTrackingRecHit & | hit | ||
) | const [virtual] |
Returns pair( true, value) if the TrajectoryStateOnSurface is compatible with the RecHit, and pair( false, value) if it is not compatible. The TrajectoryStateOnSurface must be on the same Surface as the RecHit. For an estimator where there is no value computed, e.g. fixed window estimator, only the first(bool) part is of interest.
Implements Chi2MeasurementEstimatorBase.
Definition at line 9 of file Chi2Strip1DEstimator.cc.
References funct::C, TransientTrackingRecHit::detUnit(), GeomDetType::isTrackerPixel(), TrajectoryStateOnSurface::localError(), TrackingRecHit::localPosition(), TrajectoryStateOnSurface::localPosition(), TrackingRecHit::localPositionError(), m, Topology::measurementError(), Topology::measurementPosition(), LocalTrajectoryError::positionError(), dttmaxenums::R, alignCSCRings::r, GeomDetUnit::topology(), GeomDetUnit::type(), and x.
{ if(//hit.isMatched() || hit.detUnit()->type().isTrackerPixel() ) return HitReturnType(false,0.); const StripTopology* topology = dynamic_cast<const StripTopology*>(&(hit.detUnit()->topology())); double m = topology->measurementPosition(hit.localPosition()).x(); double x = topology->measurementPosition(state.localPosition()).x(); double V = topology->measurementError(hit.localPosition(), hit.localPositionError()).uu(); double C = topology->measurementError(state.localPosition(), state.localError().positionError()).uu(); double r = m - x; double R = V + C; double est = r*r/R; return returnIt(est); }