TwoHitInBoxAlgorithms.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/numerics/Weight.h>
#include <cmath>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    template<class AHitInBoxAlgorithm, class AnotherHitInBoxAlgorithm>
    class TwoHitInBoxAlgorithm {
    public:
      using HoughBox = typename AHitInBoxAlgorithm::HoughBox;
 
      template<class AObject>
      Weight operator()(const AObject& object,
                        const HoughBox* box)
      {
        AHitInBoxAlgorithm hitInBoxAlgorithm;
        AnotherHitInBoxAlgorithm anotherHitInBoxAlgorithm;
 
        const Weight& firstResult = hitInBoxAlgorithm(object, box);
        const Weight& secondResult = anotherHitInBoxAlgorithm(object, box);
 
        if (std::isnan(firstResult) and std::isnan(secondResult)) {
          return NAN;
        } else if (std::isnan(firstResult)) {
          return secondResult;
        } else if (std::isnan(secondResult)) {
          return firstResult;
        } else {
          return firstResult + secondResult;
        }
      }
    };
  }
}