FilterSelector.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Nils Braun                                               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/findlets/base/Findlet.h>
 
#include <tracking/trackFindingCDC/numerics/Weight.h>
 
#include <tracking/trackFindingCDC/utilities/WeightedRelation.h>
#include <tracking/trackFindingCDC/utilities/Algorithms.h>
 
#include <algorithm>
#include <vector>
#include <string>
 
namespace Belle2 {
  class ModuleParamList;
  namespace TrackFindingCDC {
    template <class ACollectorItem, class ACollectionItem, class AFilter>
    class FilterSelector :
      public Findlet<WeightedRelation<ACollectorItem, const ACollectionItem>> {
    public:
      using WeightedRelationItem = WeightedRelation<ACollectorItem, const ACollectionItem>;
 
      using Super = Findlet<WeightedRelation<ACollectorItem, const ACollectionItem>>;
 
      FilterSelector() : Super()
      {
        Super::addProcessingSignalListener(&m_filter);
      }
 
      void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override
      {
        Super::exposeParameters(moduleParamList, prefix);
        m_filter.exposeParameters(moduleParamList, prefix);
      }
 
      void apply(std::vector<WeightedRelationItem>& weightedRelations) override
      {
        for (WeightedRelationItem& weightedRelation : weightedRelations) {
          const Weight weight = m_filter(weightedRelation);
          weightedRelation.setWeight(weight);
        }
 
        const auto& weightIsNan = [](const WeightedRelationItem & item) {
          return std::isnan(item.getWeight());
        };
 
        // Erase all items with a weight of NAN
        erase_remove_if(weightedRelations, weightIsNan);
 
        std::sort(weightedRelations.begin(), weightedRelations.end());
      }
 
    private:
      AFilter m_filter;
    };
  }
}