Range.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/trackFindingVXD/filterMap/filterFramework/TBranchLeafType.h>
#include <TBranch.h>
#include <TTree.h>
 
#include <framework/logging/Logger.h>
 
namespace Belle2 {
  template<typename InfType, typename SupType>
  class Range {
    const char* c_infSuffix = "_inf";
    const char* c_supSuffix = "_sup";
 
    InfType m_inf;
    SupType m_sup;
  public:
 
    Range(InfType inf, SupType sup): m_inf(inf), m_sup(sup) {};
 
    Range(): m_inf(0.), m_sup(0.) {};
    template<class VariableType>
    inline bool contains(const VariableType& x) const { return m_inf < x && x < m_sup ;};
 
    void persist(TTree* t, const std::string& branchName, const std::string& variableName)
    {
 
      std::string leafList;
      leafList += variableName;
      leafList += c_infSuffix;
      leafList += "/";
      leafList += TBranchLeafType(m_inf);
      leafList += ":";
      leafList += variableName;
      leafList += c_supSuffix;
      leafList += "/";
      leafList += TBranchLeafType(m_sup);
      TBranch* branch = new TBranch(t, branchName.c_str(), & m_inf, leafList.c_str());
      t->GetListOfBranches()->Add(branch);
    }
 
    void setBranchAddress(TTree* t, const std::string& branchName,
                          const std::string& /*variableName*/)
    {
      if (t->SetBranchAddress(branchName.c_str(), & m_inf) < 0) B2FATAL("Range: branch address not valid");
    }
 
    InfType getInf(void) const { return m_inf; } ;
 
    SupType getSup(void) const { return m_sup; } ;
 
 
    std::string getNameAndReference(std::vector<std::pair<char, void*>>* pointers = nullptr, const std::string& varname = "X")
    {
      std::string minVal = std::to_string(m_inf);
      std::string maxVal = std::to_string(m_sup);
      // if pointer to vector is provided fill it
      if (pointers != nullptr) {
        // use the position in the vector as unique identifier
        minVal = "#" + std::to_string(pointers->size());
        (*pointers).push_back({TBranchLeafType(m_inf), &m_inf});
        maxVal = "#" + std::to_string(pointers->size());
        (*pointers).push_back({TBranchLeafType(m_sup), &m_sup});
      }
      return ("(" + minVal + " < " + varname + " < " + maxVal + ")");
    }
  };
}