ResolutionFunctions.hh

// -*- C++ -*-
#ifndef RIVET_ResolutionFunctions_HH
#define RIVET_ResolutionFunctions_HH
 
#include "Rivet/Jet.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Math/Vector4.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Tools/RivetYODA.hh"
 
 
namespace Rivet {
 
  // TODO: define functions as well as functors? Not needed (yet)
 
  // (largely) virtual base class for resolution functions that require refdata,
  // which should be loaded at construction (in init) only.
  // TODO: restrict to Jet/Particle (3/4-mom?)
  template <typename T>
  class ResolutionFunctor {
  public:

    virtual ~ResolutionFunctor() = default;
    virtual double resolution(const T& input) const = 0;
 

 
    void _cacheRefData() const {
      // For thread safety, as in Analysis.cc
      static std::mutex cache_ref_data_mutex_guard{};
      std::lock_guard<mutex> lock(cache_ref_data_mutex_guard);
 
      if (_refdata.empty()) {
        MSG_TRACE("Getting refdata cache for file " << _refFileName);
        _refdata = getRefData(_refFileName);
      }
    }

    // Based on Analysis.cc
    template <typename T2 = YODA::Estimate1D>
    const T2& refData(const string& hname) const {
      _cacheRefData();
      if (!_refdata[hname]) {
        MSG_ERROR("Can't find reference histogram " << hname);
        throw Exception("Reference data " + hname + " not found.");
      }
      try {
        return dynamic_cast<T2&>(*_refdata[hname]);
      }
      catch (std::bad_cast& e) {
        throw Exception("Requested object type "s + typeid(T2).name() + ", but reference data \""s + hname
                        + "\" actually has type " + _refdata[hname]->type() + ".\n");
      }
    }

 
    Log& getLog() const {
      return Log::getLog("Rivet.ResolutionFunctor");
    }
 
  protected:
 
    mutable std::map<std::string, YODA::AnalysisObjectPtr> _refdata;
    string _refFileName;
  };
 
  template <typename T>
  using ResolutionFunctorPtr = unique_ptr<ResolutionFunctor<T>>;
 
  // TODO: Is there an useful argument for a ConstResolutionFunctor?
 
  // TODO: We could auto-generate smearing functions from pairs of
  // (pT, phi) res functions -- add machinery?
}
 
#endif