RooFoamGenerator.h

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitCore                                                       *
 *    File: $Id$
 * Authors:                                                                  *
 *   WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu       *
 *   DK, David Kirkby,    UC Irvine,         dkirkby@uci.edu                 *
 *                                                                           *
 * Copyright (c) 2000-2005, Regents of the University of California          *
 *                          and Stanford University. All rights reserved.    *
 *                                                                           *
 * 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 ROO_FOAM_GENERATOR
#define ROO_FOAM_GENERATOR
 
#include "RooAbsNumGenerator.h"
#include "RooPrintable.h"
#include "RooArgSet.h"
 
#include <TFoamIntegrand.h>
#include <TFoam.h>
 
class RooAbsReal;
class RooRealVar;
class RooDataSet;
 
class RooNumGenFactory ;
 
class RooFoamGenerator : public RooAbsNumGenerator {
public:
  RooFoamGenerator() = default;
  RooFoamGenerator(const RooAbsReal &func, const RooArgSet &genVars, const RooNumGenConfig& config, bool verbose=false, const RooAbsReal* maxFuncVal=nullptr);
  RooAbsNumGenerator* clone(const RooAbsReal& func, const RooArgSet& genVars, const RooArgSet& /*condVars*/,
             const RooNumGenConfig& config, bool verbose=false, const RooAbsReal* maxFuncVal=nullptr) const override {
    return new RooFoamGenerator(func,genVars,config,verbose,maxFuncVal) ;
  }
 
  const RooArgSet *generateEvent(UInt_t remaining, double& resampleRatio) override;
 
  TFoam& engine() { return *_tfoam; }
 
  bool canSampleConditional() const override { return false ; }
  bool canSampleCategories() const override { return false ; }
 
  std::string const& generatorName() const override;
 
protected:
 
  friend class RooNumGenFactory ;
  static void registerSampler(RooNumGenFactory& fact) ;
 
  std::unique_ptr<TFoamIntegrand> _binding; ///< Binding of RooAbsReal to TFoam function interface
  std::unique_ptr<TFoam> _tfoam;             ///< Instance of TFOAM generator
  std::vector<double> _xmin ;                ///< Lower bound of observables to be generated ;
  std::vector<double> _range;                ///< Range of observables to be generated ;
  std::vector<double> _vec;                  ///< Transfer array for FOAM output
};
 
#endif