VarSet.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/varsets/BaseVarSet.h>
 
#include <tracking/trackFindingCDC/varsets/FixedSizeNamedFloatTuple.h>
 
#include <vector>
#include <string>
#include <limits>
#include <cmath>
 
namespace Belle2 {
  namespace TrackFindingCDC {
    template<class AVarNames>
    class VarSet : public BaseVarSet<typename AVarNames::Object> {
 
    private:
      using Super = BaseVarSet<typename AVarNames::Object>;
 
    public:
      using Object = typename Super::Object;
 
    private:
      static const size_t nVars = AVarNames::nVars;
 
    public:
      std::vector<Named<Float_t*>> getNamedVariables(const std::string& prefix) override
      {
        return m_variables.getNamedVariables(prefix);
      }
 
      MayBePtr<Float_t> find(const std::string& varName) override
      {
        return m_variables.find(varName);
      }
 
    protected:
      constexpr static int named(const char* name)
      {
        return index<nVars>(AVarNames::getName, name);
      }
 
      template<int I>
      Float_t get() const
      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables.get(I);
      }
 
      template<int I>
      Float_t& var()
      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return m_variables[I];
      }
 
      template<typename AFloat>
      struct AssignFinite {
        explicit AssignFinite(AFloat& value)
          : m_value(value)
        {
        }
 
        operator AFloat& ()
        {
          return m_value;
        }
 
        void operator=(const AFloat value)
        {
          m_value = value;
          if (not std::isfinite(value)) {
            m_value = std::copysign(std::numeric_limits<AFloat>::max(), value);
          }
        }
 
      private:
        AFloat& m_value;
      };
 
      template<int I>
      AssignFinite<Float_t> finitevar()
      {
        static_assert(I < nVars, "Requested variable index exceeds number of variables.");
        return AssignFinite<Float_t>(m_variables[I]);
      }
 
    private:
      FixedSizeNamedFloatTuple<AVarNames> m_variables;
    };
  }
}