ActsTrk::FitterHelperFunctions::ATLASOutlierFinder Struct Reference

Outlier finder using a Chi2 cut. More...

#include <FitterHelperFunctions.h>

Collaboration diagram for ActsTrk::FitterHelperFunctions::ATLASOutlierFinder:

Public Member Functions
template<typename trajectory_t >
bool	operator() (typename Acts::MultiTrajectory< trajectory_t >::ConstTrackStateProxy state) const
	Classify a measurement as a valid one or an outlier. More...

Public Attributes
double	StateChiSquaredPerNumberDoFCut = std::numeric_limits<double>::max()

Detailed Description

Outlier finder using a Chi2 cut.

Definition at line 39 of file FitterHelperFunctions.h.

Member Function Documentation

operator()()

template<typename trajectory_t >

bool ActsTrk::FitterHelperFunctions::ATLASOutlierFinder::operator() ( typename Acts::MultiTrajectory< trajectory_t >::ConstTrackStateProxy  state ) const

inline

Classify a measurement as a valid one or an outlier.

Template Parameters

track_state_t

Type of the track state

Parameters

state

The track state to classify

Return values

False	if the measurement is not an outlier
True	if the measurement is an outlier

Definition at line 48 of file FitterHelperFunctions.h.

                                                                                                {
      // can't determine an outlier w/o a measurement or predicted parameters
      if (not state.hasCalibrated() or not state.hasPredicted()) {
        return false;
      }
      return Acts::visit_measurement(
          state.calibratedSize(),
    [&] (auto N) -> bool {
      constexpr size_t kMeasurementSize = decltype(N)::value;
 
      typename Acts::TrackStateTraits<kMeasurementSize, true>::Measurement calibrated{
        state.template calibrated<Acts::MultiTrajectoryTraits::MeasurementSizeMax>().data()};
      
      typename Acts::TrackStateTraits<kMeasurementSize, true>::MeasurementCovariance
        calibratedCovariance{state.template calibratedCovariance<Acts::MultiTrajectoryTraits::MeasurementSizeMax>().data()};
      
      // Take the projector (measurement mapping function)
            const auto H =
                state.projector()
                    .template topLeftCorner<kMeasurementSize, Acts::BoundIndices::eBoundSize>()
                    .eval();
      
      const auto residual = calibrated - H * state.predicted();
      double chi2 = (residual.transpose() * (calibratedCovariance + H * state.predictedCovariance() * H.transpose()).inverse() * residual).value();     
            return (chi2 > StateChiSquaredPerNumberDoFCut * kMeasurementSize);
          });
    }

Member Data Documentation

StateChiSquaredPerNumberDoFCut

double ActsTrk::FitterHelperFunctions::ATLASOutlierFinder::StateChiSquaredPerNumberDoFCut = std::numeric_limits<double>::max()

Definition at line 40 of file FitterHelperFunctions.h.

---

The documentation for this struct was generated from the following file:

• FitterHelperFunctions.h