HistFactoryImpl.h

Go to the documentation of this file.

/*
 * Project: RooFit
 * Authors:
 *   Jonas Rembser, CERN 2023
 *
 * 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 HistFactoryImplHelpers_h
#define HistFactoryImplHelpers_h
 
#include <RooGlobalFunc.h>
#include <RooWorkspace.h>
 
#include <ROOT/RSpan.hxx>
 
namespace RooStats::HistFactory {
 
namespace Constraint {
 
enum Type {
   Gaussian,
   Poisson
};
std::string Name(Type type);
Type GetType(const std::string &Name);
 
} // namespace Constraint
 
namespace Detail {
 
namespace MagicConstants {
 
constexpr double defaultGammaMin = 0;
constexpr double defaultShapeFactorGammaMax = 1000;
constexpr double defaultShapeSysGammaMax = 10;
constexpr double defaultStatErrorGammaMax = 10;
constexpr double minShapeUncertainty = 0.0;
 
} // namespace MagicConstants
 
template <class Arg_t, class... Params_t>
Arg_t &getOrCreate(RooWorkspace &ws, std::string const &name, Params_t &&...params)
{
   Arg_t *arg = static_cast<Arg_t *>(ws.obj(name));
   if (arg)
      return *arg;
   Arg_t newArg(name.c_str(), name.c_str(), std::forward<Params_t>(params)...);
   ws.import(newArg, RooFit::RecycleConflictNodes(true), RooFit::Silence(true));
   return *static_cast<Arg_t *>(ws.obj(name));
}
 
void configureConstrainedGammas(RooArgList const &gammas, std::span<const double> relSigmas, double minSigma);
 
struct CreateGammaConstraintsOutput {
   std::vector<std::unique_ptr<RooAbsPdf>> constraints;
   std::vector<RooRealVar *> globalObservables;
};
 
CreateGammaConstraintsOutput createGammaConstraints(RooArgList const &paramList, std::span<const double> relSigmas,
                                                    double minSigma, Constraint::Type type);
 
} // namespace Detail
} // namespace RooStats::HistFactory
 
#endif