Logo ROOT  
Reference Guide
 
Loading...
Searching...
No Matches
RooNumIntFactory.cxx
Go to the documentation of this file.
1/*****************************************************************************
2 * Project: RooFit *
3 * Package: RooFitCore *
4 * @(#)root/roofitcore:$Id$
5 * Authors: *
6 * WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu *
7 * DK, David Kirkby, UC Irvine, dkirkby@uci.edu *
8 * *
9 * Copyright (c) 2000-2005, Regents of the University of California *
10 * and Stanford University. All rights reserved. *
11 * *
12 * Redistribution and use in source and binary forms, *
13 * with or without modification, are permitted according to the terms *
14 * listed in LICENSE (http://roofit.sourceforge.net/license.txt) *
15 *****************************************************************************/
16
17/**
18\file RooNumIntFactory.cxx
19\class RooNumIntFactory
20\ingroup Roofitcore
21
22RooNumIntFactory is a factory to instantiate numeric integrators
23from a given function binding and a given configuration. The factory
24searches for a numeric integrator registered with the factory that
25has the ability to perform the numeric integration. The choice of
26method may depend on the number of dimensions integrated,
27the nature of the integration limits (closed or open ended) and
28the preference of the caller as encoded in the configuration object.
29**/
30
31#include "TClass.h"
32#include "TSystem.h"
33#include "Riostream.h"
34
35#include "RooFit.h"
36
37#include "RooNumIntFactory.h"
38#include "RooArgSet.h"
39#include "RooAbsFunc.h"
40#include "RooNumIntConfig.h"
41#include "RooNumber.h"
42
43#include "RooIntegrator1D.h"
44#include "RooBinIntegrator.h"
45#include "RooIntegrator2D.h"
49#include "RooMCIntegrator.h"
50//#include "RooGaussKronrodIntegrator1D.h"
51//#include "RooAdaptiveGaussKronrodIntegrator1D.h"
53
54#include "RooMsgService.h"
55
56using namespace std ;
57
59
60
61
62////////////////////////////////////////////////////////////////////////////////
63/// Register all known integrators by calling
64/// their static registration functions
65void RooNumIntFactory::init() {
73 // GSL integrator is now in RooFitMore and it register itself
74 //RooAdaptiveGaussKronrodIntegrator1D::registerIntegrator(*this) ;
75 //RooGaussKronrodIntegrator1D::registerIntegrator(*this) ;
77
78 RooNumIntConfig::defaultConfig().method1D().setLabel("RooIntegrator1D") ;
79 RooNumIntConfig::defaultConfig().method1DOpen().setLabel("RooImproperIntegrator1D") ;
80 RooNumIntConfig::defaultConfig().method2D().setLabel("RooAdaptiveIntegratorND") ;
81 RooNumIntConfig::defaultConfig().methodND().setLabel("RooAdaptiveIntegratorND") ;
82
83 //if GSL is available load (and register GSL integrator)
84#ifdef R__HAS_MATHMORE
85 int iret = gSystem->Load("libRooFitMore");
86 if (iret < 0) {
87 oocoutE((TObject*)nullptr, Integration) << " RooNumIntFactory::Init : libRooFitMore cannot be loaded. GSL integrators will not beavailable ! " << std::endl;
88 }
89#endif
90}
91
92
93////////////////////////////////////////////////////////////////////////////////
94/// Static method returning reference to singleton instance of factory
95
97{
98 static unique_ptr<RooNumIntFactory> instance;
99
100 if (!instance) {
101 // This is needed to break a deadlock. During init(),
102 // other functions may call back to this one here. So we need to construct first,
103 // and ensure that we can return an instance while we are waiting for init()
104 // to finish.
105 instance.reset(new RooNumIntFactory);
106 instance->init();
107 }
108
109 return *instance;
110}
111
112
113
114////////////////////////////////////////////////////////////////////////////////
115/// Method accepting registration of a prototype numeric integrator along with a RooArgSet of its
116/// default configuration options and an optional list of names of other numeric integrators
117/// on which this integrator depends. Returns true if integrator was previously registered
118
120{
121 TString name = proto->IsA()->GetName() ;
122
124 //cout << "RooNumIntFactory::storeIntegrator() ERROR: integrator '" << name << "' already registered" << endl ;
125 return kTRUE ;
126 }
127
128 // Add to factory
129 _map[name.Data()] = std::make_pair(unique_ptr<RooAbsIntegrator>(proto), std::string(depName));
130
131 // Add default config to master config
133
134 return kFALSE ;
135}
136
137
138
139////////////////////////////////////////////////////////////////////////////////
140/// Return prototype integrator with given (class) name
141
143{
144 auto item = _map.find(name);
145
146 return item == _map.end() ? nullptr : item->second.first.get();
147}
148
149
150
151////////////////////////////////////////////////////////////////////////////////
152/// Get list of class names of integrators needed by integrator named 'name'
153
154const char* RooNumIntFactory::getDepIntegratorName(const char* name) const
155{
156 auto item = _map.find(name);
157
158 return item == _map.end() ? nullptr : item->second.second.c_str();
159}
160
161
162
163////////////////////////////////////////////////////////////////////////////////
164/// Construct a numeric integrator instance that operates on function 'func' and is configured
165/// with 'config'. If ndimPreset is greater than zero that number is taken as the dimensionality
166/// of the integration, otherwise it is queried from 'func'. This function iterators over list
167/// of available prototype integrators and returns an clone attached to the given function of
168/// the first class that matches the specifications of the requested integration considering
169/// the number of dimensions, the nature of the limits (open ended vs closed) and the user
170/// preference stated in 'config'
171
173{
174 // First determine dimensionality and domain of integrand
175 Int_t ndim = ndimPreset>0 ? ndimPreset : ((Int_t)func.getDimension()) ;
176
177 Bool_t openEnded = kFALSE ;
178 Int_t i ;
179 for (i=0 ; i<ndim ; i++) {
182 openEnded = kTRUE ;
183 }
184 }
185
186 // Find method defined configuration
187 TString method ;
188 switch(ndim) {
189 case 1:
190 method = openEnded ? config.method1DOpen().getCurrentLabel() : config.method1D().getCurrentLabel() ;
191 break ;
192
193 case 2:
194 method = openEnded ? config.method2DOpen().getCurrentLabel() : config.method2D().getCurrentLabel() ;
195 break ;
196
197 default:
198 method = openEnded ? config.methodNDOpen().getCurrentLabel() : config.methodND().getCurrentLabel() ;
199 break ;
200 }
201
202 // If distribution is binned and not open-ended override with bin integrator
203 if (isBinned & !openEnded) {
204 method = "RooBinIntegrator" ;
205 }
206
207 // Check that a method was defined for this case
208 if (!method.CompareTo("N/A")) {
209 oocoutE((TObject*)0,Integration) << "RooNumIntFactory::createIntegrator: No integration method has been defined for "
210 << (openEnded?"an open ended ":"a ") << ndim << "-dimensional integral" << endl ;
211 return 0 ;
212 }
213
214 // Retrieve proto integrator and return clone configured for the requested integration task
215 const RooAbsIntegrator* proto = getProtoIntegrator(method) ;
216 RooAbsIntegrator* engine = proto->clone(func,config) ;
217 if (config.printEvalCounter()) {
218 engine->setPrintEvalCounter(kTRUE) ;
219 }
220 return engine ;
221}
#define oocoutE(o, a)
int Int_t
Definition RtypesCore.h:45
const Bool_t kFALSE
Definition RtypesCore.h:101
const Bool_t kTRUE
Definition RtypesCore.h:100
#define ClassImp(name)
Definition Rtypes.h:364
char name[80]
Definition TGX11.cxx:110
R__EXTERN TSystem * gSystem
Definition TSystem.h:559
const char * proto
Definition civetweb.c:16613
virtual const char * getCurrentLabel() const
Return label string of current state.
Abstract interface for evaluating a real-valued function of one real variable and performing numerica...
Definition RooAbsFunc.h:27
virtual Double_t getMinLimit(UInt_t dimension) const =0
virtual Double_t getMaxLimit(UInt_t dimension) const =0
UInt_t getDimension() const
Definition RooAbsFunc.h:33
RooAbsIntegrator is the abstract interface for integrators of real-valued functions that implement th...
virtual RooAbsIntegrator * clone(const RooAbsFunc &function, const RooNumIntConfig &config) const =0
void setPrintEvalCounter(Bool_t value)
static void registerIntegrator(RooNumIntFactory &fact)
Register RooAdaptiveIntegratorND, its parameters, dependencies and capabilities with RooNumIntFactory...
RooArgSet is a container object that can hold multiple RooAbsArg objects.
Definition RooArgSet.h:35
static void registerIntegrator(RooNumIntFactory &fact)
Register RooBinIntegrator, is parameters and capabilities with RooNumIntFactory.
virtual Bool_t setLabel(const char *label, bool printError=true) override
Set value by specifying the name of the desired state.
static void registerIntegrator(RooNumIntFactory &fact)
Register RooImproperIntegrator1D, its parameters and capabilities with RooNumIntFactory.
static void registerIntegrator(RooNumIntFactory &fact)
Register RooIntegrator1D, is parameters and capabilities with RooNumIntFactory.
static void registerIntegrator(RooNumIntFactory &fact)
Register RooIntegrator2D, is parameters and capabilities with RooNumIntFactory.
static void registerIntegrator(RooNumIntFactory &fact)
This function registers class RooMCIntegrator, its configuration options and its capabilities with Ro...
RooNumIntConfig holds the configuration parameters of the various numeric integrators used by RooReal...
RooCategory & method2DOpen()
RooCategory & method2D()
RooCategory & methodND()
RooCategory & methodNDOpen()
Bool_t printEvalCounter() const
RooCategory & method1D()
RooCategory & method1DOpen()
Bool_t addConfigSection(const RooAbsIntegrator *proto, const RooArgSet &defaultConfig)
Add a configuration section for a particular integrator.
static RooNumIntConfig & defaultConfig()
Return reference to instance of default numeric integrator configuration object.
RooNumIntFactory is a factory to instantiate numeric integrators from a given function binding and a ...
std::map< std::string, std::pair< std::unique_ptr< RooAbsIntegrator >, std::string > > _map
RooAbsIntegrator * createIntegrator(RooAbsFunc &func, const RooNumIntConfig &config, Int_t ndim=0, Bool_t isBinned=kFALSE) const
Construct a numeric integrator instance that operates on function 'func' and is configured with 'conf...
void init()
Register all known integrators by calling their static registration functions.
const RooAbsIntegrator * getProtoIntegrator(const char *name) const
Return prototype integrator with given (class) name.
const char * getDepIntegratorName(const char *name) const
Get list of class names of integrators needed by integrator named 'name'.
Bool_t storeProtoIntegrator(RooAbsIntegrator *proto, const RooArgSet &defConfig, const char *depName="")
Method accepting registration of a prototype numeric integrator along with a RooArgSet of its default...
static RooNumIntFactory & instance()
Static method returning reference to singleton instance of factory.
static Int_t isInfinite(Double_t x)
Return true if x is infinite by RooNumBer internal specification.
Definition RooNumber.cxx:58
static void registerIntegrator(RooNumIntFactory &fact)
Register RooSegmentedIntegrator1D, its parameters, dependencies and capabilities with RooNumIntFactor...
static void registerIntegrator(RooNumIntFactory &fact)
Register RooSegmentedIntegrator2D, its parameters, dependencies and capabilities with RooNumIntFactor...
Mother of all ROOT objects.
Definition TObject.h:41
Basic string class.
Definition TString.h:136
int CompareTo(const char *cs, ECaseCompare cmp=kExact) const
Compare a string to char *cs2.
Definition TString.cxx:442
virtual int Load(const char *module, const char *entry="", Bool_t system=kFALSE)
Load a shared library.
Definition TSystem.cxx:1855