utility.hh

Go to the documentation of this file.

 
#pragma once
 
#include <memory>
 
#include <boost/core/demangle.hpp>
 
#include "fastjet/PseudoJet.hh"
 
namespace HEJ{
 
 
  template<class T, class... Args>
  std::unique_ptr<T> make_unique(Args&&... a){
    return std::unique_ptr<T>{new T{std::forward<Args>(a)...}};
  }
 
  template<typename T, typename... U>
  constexpr
  std::array<T, 1 + sizeof...(U)> make_array(T t, U&&... rest){
    return {{t, std::forward<U>(rest)...}};
  }
 
 
  inline
  std::string join(
      std::string const & /* delim */
  ){
    return "";
  }
 
  inline
  std::string join(
      std::string const & /* delim */, std::string const & str
  ){
    return str;
  }
 
 
  template<typename... Strings>
  std::string join(
      std::string const & delim,
      std::string const & first, std::string const & second,
      Strings&&... rest
  ){
    return join(delim, first + delim + second, std::forward<Strings>(rest)...);
  }
 
  template<typename T>
  std::string type_string(T&&){
    return boost::core::demangle(typeid(T).name());
  }
 
  template<typename... T>
  constexpr void ignore(T&&...) {}
 
  inline
  bool nearby_ep(double a, double b, double ep){
    assert(ep > 0);
    return std::abs(a-b) < ep;
  }
 
  inline
  bool nearby_ep(
      fastjet::PseudoJet const & pa, fastjet::PseudoJet const & pb,
      double ep
  ){
    assert(ep > 0);
    for(size_t i = 0; i < 4; ++i){
      if(!nearby_ep(pa[i], pb[i], ep)) return false;
    }
    return true;
  }
 
  inline
  bool nearby(
      fastjet::PseudoJet const & pa, fastjet::PseudoJet const & pb, double const norm = 1.
  ){
    return nearby_ep(pa, pb, 1e-7*norm);
  }
 
}