EvalContext.h

Go to the documentation of this file.

/*
 * Project: RooFit
 * Authors:
 *   Jonas Rembser, CERN  12/2021
 *
 * Copyright (c) 2023, CERN
 *
 * Redistribution and use in source and binary forms,
 * with or without modification, are permitted according to the terms
 * listed in LICENSE (http://roofit.sourceforge.net/license.txt)
 */
 
#ifndef RooFit_Detail_EvalContext_h
#define RooFit_Detail_EvalContext_h
 
#include <RooAbsArg.h>
 
#include <ROOT/RSpan.hxx>
 
#include <TNamed.h>
#include <TObject.h>
 
#include <Math/Util.h>
 
#include <map>
#include <stdexcept>
#include <sstream>
 
template <class T>
class RooTemplateProxy;
 
namespace RooBatchCompute {
class Config;
}
 
/// \class RooFit::DataKey
/// To use as a key type for RooFit data maps and containers. A RooFit::DataKey
/// can be constructed with no runtime overhead from a RooAbsArg (or any
/// templated RooFit proxy for convenience). Compared to using the RooAbsArg
/// pointer directly, this has the advantage that one can easily change the way
/// the key is constructed from the object, just by changing the implementation
/// of the DataKey. For example, it is trivial to move from using the RooAbsArg
/// pointer to using the unique name pointer retrieved by RooAbsArg::namePtr().
 
namespace RooFit {
namespace Detail {
 
class DataKey {
public:
   inline DataKey(RooAbsArg const *arg) : _ptr{arg->namePtr()} {}
   inline DataKey(TNamed const *arg) : _ptr{arg} {}
   template <class T>
   inline DataKey(RooTemplateProxy<T> const &proxy) : DataKey{&*proxy}
   {
   }
 
   // Comparison operators that wrap the pointer comparisons.
   friend inline bool operator==(const DataKey &k1, const DataKey &k2) { return k1._ptr == k2._ptr; }
   friend inline bool operator!=(const DataKey &k1, const DataKey &k2) { return k1._ptr != k2._ptr; }
   friend inline bool operator<(const DataKey &k1, const DataKey &k2) { return k1._ptr < k2._ptr; }
 
   // Implementing pointer-style operators.
   inline TObject const &operator*() const { return *_ptr; }
   inline TObject const *operator->() const { return _ptr; }
 
private:
   TObject const *_ptr;
};
 
} // namespace Detail
} // namespace RooFit
 
namespace std {
 
template <>
struct hash<RooFit::Detail::DataKey> {
   std::size_t operator()(const RooFit::Detail::DataKey &k) const { return hash<TObject const *>{}(&*k); }
};
 
} // namespace std
 
namespace RooFit {
 
class EvalContext {
public:
   enum class OffsetMode { WithoutOffset, WithOffset, OnlyOffset };
 
   auto size() const { return _ctx.size(); }
   void resize(std::size_t n);
 
   inline void set(RooAbsArg const *arg, std::span<const double> const &span)
   {
      if (!arg->hasDataToken())
         return;
      std::size_t idx = arg->dataToken();
      _ctx[idx] = span;
   }
 
   void setConfig(RooAbsArg const *arg, RooBatchCompute::Config const &config);
 
   std::span<const double> at(RooAbsArg const *arg, RooAbsArg const *caller = nullptr);
 
   template <class T>
   inline std::span<const double> at(RooTemplateProxy<T> const &proxy)
   {
      return at(&proxy.arg(), proxy.owner());
   }
 
   RooBatchCompute::Config config(RooAbsArg const *arg) const;
   void enableVectorBuffers(bool enable) { _enableVectorBuffers = enable; }
   void resetVectorBuffers() { _bufferIdx = 0; }
   std::span<double> output() { return _currentOutput; }
 
   void setOutputWithOffset(RooAbsArg const *arg, ROOT::Math::KahanSum<double> val,
                            ROOT::Math::KahanSum<double> const &offset);
 
private:
   friend class Evaluator;
 
   OffsetMode _offsetMode = OffsetMode::WithoutOffset;
   std::span<double> _currentOutput;
   std::vector<std::span<const double>> _ctx;
   bool _enableVectorBuffers = false;
   std::vector<std::vector<double>> _buffers;
   std::size_t _bufferIdx = 0;
   std::vector<RooBatchCompute::Config> _cfgs;
};
 
} // namespace RooFit
 
#endif