RangeInBox.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/numerics/Weight.h>
#include <cmath>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    template<class AHitInBoxAlgorithm>
    class RangeInBox {
    public:
      using HoughBox = typename AHitInBoxAlgorithm::HoughBox;
 
      template<class ARangeObject>
      Weight operator()(const ARangeObject& rangeObject,
                        const HoughBox* box)
      {
        AHitInBoxAlgorithm hitInBoxAlgorithm;
        double sumOfWeights = 0;
        double numberOfPassedItems = 0;
        for (const auto& item : rangeObject) {
          const double hitWeight = hitInBoxAlgorithm(item, box);
          if (not std::isnan(hitWeight)) {
            numberOfPassedItems++;
            sumOfWeights += hitWeight;
          }
        }
 
        if (numberOfPassedItems > minimalRatio * static_cast<double>(rangeObject.size())) {
          return sumOfWeights;
        } else {
          return NAN;
        }
      }
 
      const double minimalRatio = 0.6;
    };
  }
}