RooAcceptReject.h

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitCore                                                       *
 *    File: $Id: RooAcceptReject.h,v 1.16 2007/05/11 09:11:30 verkerke Exp $
 * 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_ACCEPT_REJECT
#define ROO_ACCEPT_REJECT
 
#include "RooAbsNumGenerator.h"
#include "RooPrintable.h"
#include "RooArgSet.h"
 
class RooAbsReal;
class RooRealVar;
class RooDataSet;
class RooRealBinding;
class RooNumGenFactory ;
 
class RooAcceptReject : public RooAbsNumGenerator {
public:
  RooAcceptReject() {
    // coverity[UNINIT_CTOR]
  } ;
  RooAcceptReject(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 RooAcceptReject(func,genVars,config,verbose,maxFuncVal) ;
  }
 
  const RooArgSet *generateEvent(UInt_t remaining, double& resampleRatio) override;
  double getFuncMax() override ;
 
 
  // Advertisement of capabilities
  bool canSampleConditional() const override { return true ; }
  bool canSampleCategories() const override { return true ; }
 
  std::string const& generatorName() const override;
 
protected:
 
  friend class RooNumGenFactory ;
  static void registerSampler(RooNumGenFactory& fact) ;
 
  void addEventToCache();
  const RooArgSet *nextAcceptedEvent();
 
  double _maxFuncVal, _funcSum;       ///< Maximum function value found, and sum of all samples made
  UInt_t _realSampleDim,_catSampleMult; ///< Number of real and discrete dimensions to be sampled
  UInt_t _minTrials;                    ///< Minimum number of max.finding trials, total number of samples
  UInt_t _totalEvents;                  ///< Total number of function samples
  UInt_t _eventsUsed;                   ///< Accepted number of function samples
 
  UInt_t _minTrialsArray[4];            ///< Minimum number of trials samples for 1,2,3 dimensional problems
};
 
#endif