CurvRep.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2015 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/hough/axes/StandardAxes.h>
#include <tracking/trackFindingCDC/utilities/Functional.h>
 
#include <type_traits>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    class CurvBinsSpec {
    public:
      CurvBinsSpec(double lowerBound, double upperBound, long nBins, int nOverlap, int nWidth);
 
      DiscreteCurv::Array constructArray() const
      { return constructLinearArray(); }
 
      DiscreteCurv::Array constructInvLinearArray() const;
 
      DiscreteCurv::Array constructLinearArray() const;
 
      DiscreteCurvWithArcLength2DCache::Array constructCacheArray() const;
 
      long getNPositions() const;
 
      double getBinWidth() const;
 
      double getOverlap() const;
 
      int getNOverlap() const
      {
        return m_nOverlap;
      }
 
    private:
      double m_lowerBound;
 
      double m_upperBound;
 
      long m_nBins;
 
      int m_nOverlap = 1;
 
      int m_nWidth = 3;
    };
 
    struct GetLowerCurv {
      template <class AHoughBox, class SFINAE = std::enable_if_t<AHoughBox::template HasType<DiscreteCurv>::value > >
      float operator()(const AHoughBox& houghBox) const
      {
        return static_cast<float>(houghBox.template getLowerBound<DiscreteCurv>());
      }
 
      template <class AHoughBox, class SFINAE = std::enable_if_t<AHoughBox::template HasType<ContinuousCurv>::value > >
      double operator()(const AHoughBox& houghBox) const
      {
        return static_cast<float>(houghBox.template getLowerBound<ContinuousCurv>());
      }
    };
 
    struct GetUpperCurv {
      template <class AHoughBox, class SFINAE = std::enable_if_t<AHoughBox::template HasType<DiscreteCurv>::value> >
      float operator()(const AHoughBox& houghBox) const
      {
        return static_cast<float>(houghBox.template getUpperBound<DiscreteCurv>());
      }
 
      template <class AHoughBox, class SFINAE = std::enable_if_t<AHoughBox::template HasType<ContinuousCurv>::value> >
      double operator()(const AHoughBox& houghBox) const
      {
        return static_cast<float>(houghBox.template getUpperBound<ContinuousCurv>());
      }
    };
 
    template <class AHoughBox>
    float getLowerCurv(const AHoughBox& houghBox)
    {
      return getIfApplicable<float>(GetLowerCurv(), houghBox, 0.0);
    }
 
    template <class AHoughBox>
    float getUpperCurv(const AHoughBox& houghBox)
    {
      return getIfApplicable<float>(GetUpperCurv(), houghBox, 0.0);
    }
  }
}