TrackL1NormResidualLessThan.cpp

Go to the documentation of this file.

 #include "NSL/Tracker/TrackL1NormResidualLessThan.h"

NSL::Selections::Track::L1NormResidualLessThan::L1NormResidualLessThan(float max, NAIA::TrTrack::Fit fit) {
  m_matcher = std::make_shared<boolMatcher>([=](Event &event) {
    const auto &testContainer = event.trTrackBase->TrTrackHitPos;

    if (!event.trTrackBase->FitIDExists(fit, NAIA::TrTrack::Span::InnerL1))
      return false;

    int ndof_inner = -3;
    for (int il = 1; il < 8; il++) {
      if (NAIA::ContainsKeys(testContainer, il, NAIA::TrTrack::Y)) {
        ndof_inner++;
      }
    }

    float residual =
        event.trTrackBase->TrChiSq[fit][NAIA::TrTrack::Span::InnerL1][NAIA::TrTrack::Side::Y] * (ndof_inner + 1) -
        event.trTrackBase->TrChiSq[fit][NAIA::TrTrack::Span::InnerOnly][NAIA::TrTrack::Side::Y] * ndof_inner;

    return residual < max;
  });
}