Log of /trunk/roofit/roofitcore/inc/RooIntegrator1D.h
Parent Directory
Revision
24278 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Sun Jun 15 15:21:16 2008 UTC (6 years, 7 months ago) by
wouter
File length: 4215 byte(s)
Diff to
previous 23134
o RooMPSentinel WHERE RooGenContext RooLinearMorph RooHistFunc RooInvTransform RooFracRemainder
RooIntegrator2D RooGaussKronrodIntegrator1D RooFitResult RooFFTConvPdf RooHist RooFormula RooIntegrator1D
RooMCStudy RooGenericPdf RooIntegratorBinding RooMPSentinel RooImproperIntegrator1D RooExtendPdf RooInt
RooLinTransBinning RooLinkedList RooLinkedListElem RooGenCategory RooList RooGrid RooLinkedListIter
RooMCIntegrator RooHistPdf RooListProxy RooFormulaVar RooHistError RooLinearVar RooGenProdProj RooHashTable
- Upgrade inline documention: make sure that (almost) every member function has some descriptive comments, every
data member has a descriptive comment, make sure that each class has a description.
- Fix class description header so that it is correctly identified by THtml
Revision
20161 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Oct 1 16:28:25 2007 UTC (7 years, 3 months ago) by
wouter
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3864 byte(s)
Diff to
previous 18735
o RooAbsCache
- New abstract base class for cache objects associated with RooAbsArgs.
Caches register themselves with the owning RooAbsArg and propagate
all necessary state changes and meta information to objects in the cache
o RooCacheManager
- Concrete template implementation of RooAbsCache with an interface
and functionality styled after RooNormListManager. The template class
provide no special handling for cache contents that includes RooAbsArg
objects and does not forward hook calls. It is thus mostly intended
to cache basic types indexed on a normalization/integration configuration
o RooObjCacheManager
- Class derived from RooCacheManager<RooAbsCacheElement> that implements
all proper call forwarding and handling of state and other meta information
to RooAbsArg objects embedded in the cache elements. To facilitate this
functionality the RooAbsCacheElement defines a pure virtual method
containedArgs() that must be implemented to expose its contents to
the cache manager. The default behavior on a serverRedirect hook is
to sterilize the cache. OperMode changes are not forwarded by default
to the cache contents as cached objects do not necessarily have the
same operation mode (e.g. if the cache contains an integral of its owner).
Any calls to optimizeCacheMode(), printCompactTreeHook(), and findConstantNodes()
are forwarded. Moreover, optimizeCacheMode() calls are remembered and
any objects inserted in the cache after such a call was received will be
optimized upon insertion in the cache
o RooAbsCacheElement
- Abstract base class for objects to be managed by RooObjCacheManager. Each
object must in addition to its payload implement the containedArgs()
member function that exposes the RooAbsArg object content to the manager
for processing. The argument passed to containedArgs() indicates the use
case for the returned list of objects, allowing the cache payload to
return different lists depending upon use. This is for most applications
not required.
o RooAbsCachedPdf
- Abstract base class derived from RooAbsPdf for p.d.f.s that implement
value caching in a grid over the entire domain of the observables.
The base class implements the overall caching logistics and manages
the RooDataHists and RooHistPdfs that respresent the cached contents.
Derived classes must inputBaseName() defined a unique name prefix for
all cache objects, actualObservales() and actualParameters() to
determined the actual observables and parameters given a set of
available observables and fillCacheObject() which should fill the
RooDataHist cache for a given obs/params configuration with the
p.d.fs values.
o RooAbsSelfCachedPdf
- Abstract base class derived from RooAbsCachedPdf that implements
all of the required logistics in RooCachedPdf and only require
derived classes from RooAbsSelfCachedPdf to implement evaluate().
Any class derived from RooAbsPdf can thus be trivially switched
to derived from RooAbsSelfCachedPdf and gain the caching functionality.
o RooCachedPdf
- Concrete class derived from RooAbsCachedPdf that implements
the caching of an externally provided p.d.f in the constructor.
o RooFFTConvPdf
- Concrete class derived from RooAbsCachedPdf that implement generic
numeric convolution of two input p.d.f.s using FFT technique
and convolution theorem. Actual FFTs are calculated by FFTW
through TVirtualFFT interface in ROOT. Present implementation only supports
convolution in one dimension, but this will be extended later
as there are no fundamental limitations in RooAbsCachedPdf nor FFTW.
o RooHistPdf
- Moved from RooFit to RooFitCore because it is used by caching classes
o RooProdPdf
- Use new RooObjCacheManager to store partial objects
- Fix memory leak (already patched in 15.0.14h)
o RooDataHist
- Use new RooCacheManager template to cache partial bin volumes
- Add new set(),setAllWeights() method to update current weight and all weights
respectively, for use in cached p.d.f.s
- Add sliceIterator() that iterates over all bins in a lower dimensional
slice of the dataset
o RooDataHistSliceIter
- Class derived from TIterator that implements the iteration
over a slice in a reduced number of dimensions in a RooDataHist
o RooMsgService
- Speed up processing of debug streams in cases where no debug streams are active
- add macros series oocoutX and oocxcoutX that allow to associate message with
a different object than *this
o RooRealBinding
- Interface new setValue method that allows value to be set outside default
domain, but inside given alternate named range. For use in numeric integrator
bindings
o RooAbsRealLValue
- Add new inRange() method with range name specifier to check against chosen
range instead of default range
o RooCmdConfig
- Fix dangling return pointer problem in decodeStringOnTheFly() method
o RooNameSet
- Fix bug in refoll()
o RooChi2Var, RooNLLVar
- Add processing of AddCoefRange() named argument
o RooConvCoefVar
- Hold p.d.f in proxy without value/shape state propagation to obtain
correct dependency of RooConvCoefVar through separate set proxy
with parameters that do have state propagation
- In getAnalyticalIntegral() ask for integration coefficient specific
to coefficient (following interface change in RooAbsAnaConv)
o RooEfficiency
- Follow interface change in RooCategoryProxy
o RooAbsAnaConvPdf
- Use new RooObjCacheManager to hold normalization and coefficient representations
- Extend getCoefAnalyticalIntegral() interface so that each coefficient
can have a separate integration configuration
- Modify makeCoefVars() such that it constructs the RooConvCoefVars with
actual dependents for each coefficients, rather than inserting dependencies
of any of the coefficients in each coefficient. This logic requires support
from the implementation class that is interface through a new member function
coefVars()
o RooAbsReal
- Add Boolean data member treeVar that indicates is value buffer is also
branch buffer of a tree. This information is already present in the RooAbsReal
but is relatively expensive to retrieve and is always needed in the inner
loop of data set operations. This extra member significantly speeds up operations
- Remove superfluous _norm[] token in name returned by integralNameSuffix() if
no normalization set is supplied.
- Add a fillDataHist() method that fills a given RooDataHist with a sampling
of the given function. Added for support of caching p.d.f.s
- Downgrading warning on non-dependent projection variables in plotOn() to
a debug level message
- Remove function optimizeDirty(),doConstOpt(),undoConstOpt(),findCacheableBranches()
findUnusedDatasetVariables() and findRedundantCacheServers() as part of reorganization
of optimization functionality.
o RooAbsCollection
- Use iterator semantics instead of vector semantics to process RooLinkedList
contents in assignment operator
- Add utility function contentsString() returning stl string with names
of contents
o RooAddPdf
- Rewrite cache management using new RooObjCacheManager that holds projection
integrals, range integrals and supplemental normalization lists
o RooAbsOptGoodnessOfFit
- Improve handling of normalization ranges and reference ranges for test statistics
containing RooAddPdfs. New named argument in fitTo() and RooNllVar::ctor() pass
reference range for fractions in RooAddPdf components that do not have an explicit
fixed definition already. Name all normalization ranges explicitly after the
subrange to which the're applied to avoid name clashes and conflicts if multiple
normalization ranges are involved.
- Rewrite informational messages through RooMsgService
- Implement new constOptimizeTestStatistic() function that signals all test statistic
components of a RooAbsArg expression tree to perform their constant term optimization
- Implement new optimizeCaching() driver function that changes default lazy-evaluation-with-caching
scheme to test statistic scheme where only parameters objects propagate state changes and
where all components that depend directly or indirectly on observables have lazy evaluation
disabled. Changes to p.d.f objects are delegated to RooAbsPdf::optimizeCacheMode(). Changes
to dataset are delegated to RooAbsData::optimizeReadingWithCaching()
- Implement new optimizeConstantTerms() function that is back end for constOptimizeTestStatistic
implementation.
o RooAICRegistry, RooNormSetCache
- Remove tracking of allocation cycles of RooArgSets as this scheme is superseded by
memory pool mechanism for RooArgSet allocation
o RooMappedCategory
- Follow interface change in RooCategoryProxy.
o RooRealMPFE
- Follow reorganization of optimization functionality: implement new constOptimizeTestStatistic
function for forwarding to slave processes.
o RooRealVar
- Add new setVal() function with range name to facilite setting of values that are outside
default range, but inside specified alternate range. Needed for numeric integration and
certain RooAddPdf coefficient projection.
o RooAbsArg
- Make state management functions public [is/set][Value/ShapeDirty](), setOperMode90
- Add static _flipAClean boolean to support setACleanAdirty() static method.
- Implement mechanism for registration of cache objects with base class RooAbsCache.
Forward all calls to operModeHook(), redirectServersHook(), optimizeCacheMode()
and printCompactTreeHook() to all caches for further processing
- Define constOptimize() function and provide default implementation that forwards
to all servers
- Implement optimizeCacheMode() function that puts expression tree in optimized
caching for test-statistics operations (only cache/lazy evaluate nodes that
do not (indirectly) depend on observables
- Implement findConstantNodes() function that identifies all nodes that depend
exclusively on constant observables and >0 observables that are eligible
for precalculation and caching in test statistics operations by
RooAbsOptTestStatistic::optimizeConstantTerms()
- Define constOptimizeTestStatistics that signal an expression tree that any
embedded test statistic objects should perform their constant term optimization.
Default implementation in RooAbsArg forwards call to all servers.
o RooAbsPdf
- Use new RooObjCacheManager to manage normalization integrals cache.
- Process new SumCoefRange() option in fitTo() that sets the coeffient
interpretation range for all embedded RooAddPdf objects of which
the interpretation range was not explicitly set
- Extend interface of fixAddCoefNormalizationSet(): add extra argument
that controls if AddCoef normalization set are forcibly overruled
for objects that already have an explicit range setting.
o RooRealIntegral
- Add some informational messaging through RooMsgService, rewrite existing ones
- Speed up numeric summations by precreating the involved RooSuperCategory
in the ctor rather than in evaluate().
- Fix bug in dependency insertion: explicitly require servers to be value
server to avoid pulling in fake dependents.
- Remove obsolete data members
- Implement new setACleanADirty() function that atomically flips state of
all AClean objects to ADirty for the duration of the integration step.
Replaces the much more expensive prepareACleanFunc() that did this
only for the objects involved.
- Process factorizing terms in integral only if the list is not empty
- Process Jacobian terms in integral only if the list is not empty
o RooAddGenContext
- Adapt to cache management chanches in RooAddPdf
o RooNumIntFactory
- Use STL classes for internal memory management
o RooSetProxyy
- Implement operator new() and operator() delete as required because
RooArgSet implements them
o RooAbsGoodnessOfFit
- Store and forward AddCoefRange information for embedded RooAddPdf objects
- Implement constOptimizeTestStatistic() and forward to components
o RooErrorVar
- Follow inRange() interface change in RooAbsRealLValue
o RooConvGenContext
- Adapt code to make it work with both RooAbsAnaConv objects and
the newly introduced RooFFTConvPdf
o RooMinuit
- Follow reorganization of optimization functionality
o RooProjectedPdf
- Store projection using new RooObjCacheManager
o RooGlobalFunc
- Add RooFit::SumCoefRange() named argument for RooAbsPdf::fitTo()
and RooNLLVar::ctor() and RooChi2Var::ctor()
o RooClassFactory
- Extend existing makePdf(),makeFunction and makeClass() functions
to also accept RooAbsCategory arguments, to accept an optional string
inserted as the return value expression inserted in the
evaluate() method of the class created, as well as as option string
specifying one or more analytical integrations of the given function
expression that will be advertised in the generated getAnalyticalIntegral()
code and implemented in the generated analyticalIntegral() code
- Add makeAndCompileFunction() and makeAndCompilePdf methods that
generate the code and immediate compile it through ROOTs ACliC interface.
These functions have a simpler interface because the list of arguments
of the function to be created can be specified as a const RooArgSet&
rather than a string with a list of names, thus making it easier
for on-the-fly definition of compiled objects.
- Add defineFunction() and definePdf() methods define and compile
a function/pdf using the above functions and immediately instantiate
an object of the given type passing the specified RooArgSet arguments
to the constructor of the object. These functions make a compiled
equivalent of RooFormulaVar based on the TFormula interpreter, i.e.
RooAbsReal* f = new RooFormulaVar("f","D(1-2w)",RooArgSet(D,w)) ;
becomes
RooAbsReal* f = RooClassFactory::defineFunction("f","D(1-2w)",RooArgSet(D,w)) ;
o RooAbsCategory
- Add Boolean data member treeVar that indicates is value buffer is also
branch buffer of a tree. This information is already present in the RooAbsReal
but is relatively expensive to retrieve and is always needed in the inner
loop of data set operations. This extra member significantly speeds up operations
- Construct explit clean name for TBranch name if node is cached, equivalent
to what is done in RooAbsReal. Needed because new constOptimizer can also
cache RooAbsCategory objects
o RooArgSet
- Implement operator new() and operator delete() to implement memory pool scheme
for RooArgSets. The goal of this implementation is to insure that two RooArgSet
pointers are _never_ the same. This allows to RooNormSet cache manager to safely
and quickly index integration and normalization configurations by evaluating
the contents of a RooArgSet only the first time it sees it.
o RooSimultaneous
- Implement caching of partial integrals using new RooObjCacheManager
- Following interface change in RooCategoryProxy ;
o RooAbsData
- Defined pure virtual methods optimizeReadingWithCache() and allClientsCached()
o RooTreeData
- Implement optimizeReadingWithCache() which determines given a p.d.f to be
evaluated and a list of the nodes that are cached for it, which observables
in the dataset are no longer required because they are either not an observable
of the p.d.f, or they exclusively feed nodes that are cached. Reading of these
observables is disabled through SetBranchStatus().
- Move allClientsCached() utility function here from RooAbsReal.
o RooCategoryProxy
- Remove const char*() conversion operator in favor of a const char* label() method
to remove possibly ambiguities in C++ math expressions involving category proxies.
o LinkDef3.h
- Add newly defined classes
Revision
15668 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Jul 3 15:37:11 2006 UTC (8 years, 6 months ago) by
wverkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3891 byte(s)
Diff to
previous 11181
o RooAbsCollection
- Fix bug in printLatex() : mangle both upper and lowercase options for sibling lists instead
of just uppercase options
o RooAbsGoodnessOfFit
- Forgot to put certain print statements in if (_verbose) conditionals (C.Roat)
- Order of arguments in calls to create() sometimes wrong (verbose/ncpu swapped) (C.Roat)
o RooAbsOptGoodnessOfFit
- Fix more mistakes with if (_verbose) conditionals (C.Roat)
o RooAbsRealLValue
- Remove forgotten debug print statement
o RooAdaptiveGaussKronrodIntegrator1D,RooGaussKronrodIntegator1D,RooIntegrator1D,RooIntegrator2D,
RooSegmentedIntegrator1D,RooSegmentedIntegrator2D,RooImproperIntegrator1D,RooMCIntegator
- Change initialization process of RooNumIntFactory from dynamic process on loading to
static process in RooNumIntFactor ctor as current model has problems
o RooAddPdf
- Add forgotton initializer for _refCoefRangeName in one of the constructors
o RooArgSet
- Declare writToFile() and writeToStream() as 'const'
o RooBanner
- Increment version tag to v2.08
o RooChi2Var
- Fix copy/paste error in print statements in ctor
o RooFitResult
- Put ROOT version dependent #ifdef around TMatrix forward declaration
to make it work with both ROOT4 and ROOT5
o RooHist,RooPlot
- New implementation of pullHist/residualHist from G.Raven
o RooNumIntConfig
- Add forgotten initializers of RooCategory datamembers to ctor (fixes crashes
observed by J. Rademackers)
o RooNumIntFactory
- Perform static initialization of all numeric integrator engines in the factory constructor
(fixes oberved crashes of all numeric integration and configuration business)
o RooRealVar
- Back out change that setting asymmetric error clear symmetric error and vice versa.
It causes just too much confusion...
o RooSharedPropertiesList
- Explicitly check for and handle input null pointers
Revision
11139 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Feb 17 14:32:38 2005 UTC (9 years, 11 months ago) by
wverkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3796 byte(s)
Diff to
previous 10648
o All numeric integrators
- Change way integrators are registered: Old method involves ROOT function call that assumes gROOT exist
leading to static initialization order fiasco. New method registers initializer functions with
RooNumIntFactory that will perform actual initialization when factory is used for the 1st time
o RooNumIntConfig
- Add methods registerInitializerFunction() to allow registration of initializer functions for
each numeric integrator engine that exists
- Add processInitializers() method to process registration of numeric integrators. Call this
function in get/set/createProtoIntegrator().
- Move static pointer to singleton instance inside instance() function so that initialization
to zero is guaranteed when instance() code is executed
o RooNumIntConfig
- Move static pointer to default instance inside default() function so that initialization
to zero of pointer is guaranteed when default() code is executed
o RooIntegrator1D
- Add configuration option 'extrapolation' that allows to switch off wynn-epsilon extrapolation
and calculate error the simple way (difference in result between last two iterations).
Default remains extrapolation on.
o RooProdPdf
- Change named argument constructor to be less confusing. Constructor now takes name,title and
const RooArgSet& with regular PDFs followed by zero or more Conditional() arguments that
specificy conditional PDFS
o RooGlobalFunc
- Remove global functions Partial() and Full(), add Conditional()
o RooProdGenContext
- Follow changes in RooProdPdf
Revision
10648 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Nov 29 20:24:44 2004 UTC (10 years, 1 month ago) by
wverkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3775 byte(s)
Diff to
previous 8538
o All Classs
- Remove '#include BaBar/BaBar.hh' that was erroneously
introduced by the BaBar IOstream migration utilities
o Several class
- Initializion order changes, removal of unused variables,
addition of missing return statements to make all code
compile cleanly with '-Wall -Wno-parentheses'
o RooAbsFunc
- Make resetNumCall() const
o RooAbsIntegrator
- Add calculate() method that wrape integral() between
profiling code the track the number of integrand
function calls
o RooIntegrator1D,RooImproperIntegrator1D,RooSegmentedIntegrator1D,
RooSegmentedIntegrator2D,RooMCIntegrator
- Change interface to take RooNumIntConfig class instead
of now deprecated RooIntegratorConfig class
- Add member function that advertise capability of the integrator:
canIntegrate1D,2D,ND,OpenEnded()
o RooAbsReal
- Change normalization integral management classes
(createIntegral(), (set)defaultIntegratorConfig(),
(set)specialIntegratorConfig()) to work with new
RooNumIntConfig class.
o RooArgSet
- Add new convenience functions get/setRealValue(),
get/setCatLabel/Index(), get/setStringValue() to
manage common content of RooArgSet.
o RooConvCoefVar, RooConvGenContext
- Support class name change RooConvolutedPdf -> RooAbsAnaConfPdf
o RooCustomizer
- Add function setCloneBrancSet() that allows to specify
an external set that takes ownership of the cloned branch
nodes so that customized objects can live beyond the lifetime
of the RooCustomizer instance
o RooImproperIntegrator1D
- Add code to handle limit scenario changes (closed to open ended
and vice versa) on the fly
o RooIntegrator1D
- Add fixSteps option to terminate integral calculation after
fixed number of steps w/o extrapolation. This is more robust
for certain types of integrals
o RooRealIntegral
- Change interface to take RooNumIntConfig class instead of
RooIntegratorConfig class
- Delegate instantiation of numeric integration engines to
new class RooNumIntFactory
o RooNumIntConfig
- New class to store configuration of numeric integator engines.
This class has a brief common section with common properties
like absolute and relative precision convergence criteria,
and a section with implementation specific details for each
engine in the form of a RooArgSet. The default settings for
each existing engine are automatically added to the defaultConfig()
instance during run-time initialization
o RooNumIntFactory
- Factory class for numeric integration engines. Each existing
integrator engine automatically registers with the factory
stating its capabilities (1D,2D,ND / Open-ended). The
createIntegrator() function instantiates an integrator
appropriate for a given problem. In case more than one
integration engine can solve the given problem the (user
changable) information in the RooNumIntConfig passed along
with the call is used.
o RooGaussKronrodIntegrator1D
- New integrator engine based on Gaussian quadrature techniques.
This implementation, taken from GSL 1.5 implements the Gauss-Kronrod
techniques where the error on the integral is take from the difference
between the Gauss and Kronrod techniques. These techniques share common
sampling points allowing to calculate the error without additional cost.
This integrator applies the 10,21,43, and 87 point rules in succession
until the desired accuracy is reached. Open-ended integrals are
automatically handled through a x -> 1/x coordinate transformation
o RooAdaptiveGaussKronrodIntegrator1D
- New integrator engine based on Gaussian quadrature techniques.
This implementation, taken from GSL 1.5 implements the Gauss-Kronrod
techniques where the error on the integral is take from the difference
between the Gauss and Kronrod techniques. These techniques share common
sampling points allowing to calculate the error without additional cost.
This integrator applies the 10,21,31,41,51 or 61 point rule and
recursively subdivides the integration domain until the desired accuracy
is reached. For integrands with integrable singulaties the Wynn epsilon
rule can be selected to speed up the converges of these integrals. Open
ended integrals are automatically handles through a x -> 1/x
coordinate transformation
o RooNumConvPdf
- New operator PDF class that implements numeric convolution of
any two PDFs. The calculation of the convolution integrand
is implemented in the auxiliary class RooConvIntegrandBinding,
which in turn can be integrated with any 1-D numeric integration
engine. The default integration method for convolution integrals
is the adaptive Gauss-Kronrod method, which appears to give the
most stable results. By default the integration is performed
from -inf to +inf, but can optionally be restricted to a sliding
window using the setConvolutionWindow class, useful for relatively
narrow resolution PDFs.
o RooAbsAnaConvPdf
- New name of RooConvolutedPdf to clearly distinguish its functionality
from RooNumConvPdf
Revision
6593 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri May 9 20:48:23 2003 UTC (11 years, 8 months ago) by
wverkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3279 byte(s)
Diff to
previous 6572
Implementation of segmented 1,2D integrators, some streamlining
of intergrators class structure
o RooAbsFunc
- Add call counter as profiling hook
o RooAbsIntegrator
- Change interface of integral() to allow passing of
coordinates of non-integrated dimensions
o RooImproperIntegrator1D
- Follow interface change in RooAbsIntegrator
o RooIntegator1D
- Follow interface change in RooAbsIntegrator,
explicitly allow partial integration
- Add optional fast convergence criterium for 'zero' integrals
o RooIntegrator2D
- Use new universal integrator binding
o RooIntegratorConfig
- Add #segments field for 1D integrals
- Add min #steps for 'zero' integrals
o RooMCIntegrator
- Follow interface change in RooAbsIntegrator
o RooSegmentedIntegrator1D
- New piecewise 1D integrator. Uses array of regular
1D integrator to integrate pieces
o RooSegmentedIntegrator2D
- New piecewise 2D integrator. Uses array of
segmented 1D integrators
o RooIntegratorBinding
- New universal integrator binding, replaces RooIntegrator1DBinding
Revision
6572 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed May 7 21:06:25 2003 UTC (11 years, 8 months ago) by
wverkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 3306 byte(s)
Diff to
previous 5212
Implementation of alternate 2D integration strategy
o RooAbsPdf
- Make interface for normalization integratator configuration
a bit easier to use
o RooAbsReal
- Add new forms of createIntegral() that can take a RooIntegratorConfig
as argument
o RooIntegrator1D
- Change code such that it can work with multi-dimensional functions.
The integrated dimension is always the first. The coordinates
in the remaining dimensions are optionally passed to integral()
o RooIntegrator1DBinding
- New RooAbsFunc binding for RooIntegrator1D
o RooIntegrator2D
- New 2D integrator derived from RooIntegrator1D. Instead of
integrating RooAbsFunc, integrate a 1-D integral of RooAbsFunc,
which is effectively becomes the 2D integral of the RooAbsFunc
o RooIntegratorConfig
- Add new flag useMCFor2D() to select integration technique
for 2D numeric integrals. MC integration is default to
retain backward compatibility.
o RooRealIntegral
- Implement optional use of new RooIntegrator2D class
o RooTreeData
- Small fix in fillHistogram() for weighted datasets.
Revision
4227 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri Mar 22 22:43:58 2002 UTC (12 years, 10 months ago) by
verkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2549 byte(s)
Diff to
previous 2795
o RooAbsArg
- Keep track of cloning with attributes. add isCloneOf()
method to test if objects are clones
o RooAbsCollection
- Implement check of 'unique-name rule' in snapshot() method.
Disabled for the moment as we violate this rule in benign ways
in several places in the core code. This needs to be fixed first
o RooAbsPdf
- Add methods to specify numerical integrator configuration to be
used in normalization
There is one global default configuration (RooIntegratorConfig object)
and each PDF can optionally have a private configuration that
overrules the default one.
o RooAddPdf
- Add on-the-fly transformations for coefficients between
normalization sets. Users can now specify a
'reference normalization' via fixCoefNormalization() in
which the fractions will always be interpreted.
If no such reference set is specified the default behaviour
will prevail: the interpretation depends on the current normalization.
o RooAbsReal
- Change plotOn() such that when a PDF containing RooAddPdfs
is projected with data, an appropriate reference normalization
is set for each, if no user-supplied reference normalization
was set. This ensures intuitive behaviour of plotOn() for RooAddPdfs
- Add easy-to-use wrapper function for createProjection()
o RooErrorVar
- Add readFromStream/writToStream()
- Declare fundamental
o RooImproperIntegrator1D,RooIntegrator1D
- Add constructors from RooIntegratorConfig
o RooMCIntegrator
- Add constructor from RooIntegratorConfig
- Make #iterations and #call per iteration for grid refinement
and integration configurable parameters
o RooProdPdf
- Complete rewrite. Now _all_ permutations of integration/normalization
are correctly performed. New implementation uses a product
of partial integrals created with pdf->createIntegral()
o RooRealIntegral
- Add optional RooIntegratorConfig to ctor and pass to numerical
integrator ctors in various places
o RooSetProxy
- Fix server management for owning set proxies
o Many classes
- Anticipate potential failure of snapshot() due to unique-name
rule enforcement and handle such situations appropriately
Revision
2795 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Sat Sep 15 00:26:04 2001 UTC (13 years, 4 months ago) by
david
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2320 byte(s)
Diff to
previous 2747
Finish the optimization of RooRealIntegral's numerical engine usage,
so that whenever possible, an existing engine is reused instead of
creating a new engine (which has a high memory allocation overhead).
RooIntegrator1D:
- add setLimits() method to allow integration range to be changed
after construction and independently of the integrand's range
RooImproperIntegrator1D:
- clean up handling of different integration strategies by
introducing LimitsCase enum and limitsCase() method
- implement logic to re-use existing RooIntegrator1D objects whenever
possible in checkLimits()
RooMCIntegrator:
- implement logic to re-use ourself as long as all limits are still
finite in checkLimits()
RooIntegratorTest:
- use GenericPdf instead of FormulaVar since a FormulaVar is no
longer numerically integrable (by construction)
Revision
2747 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri Aug 24 23:55:15 2001 UTC (13 years, 5 months ago) by
david
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2266 byte(s)
Diff to
previous 2654
RooRealIntegral:
- add logic to reuse an existing numerical integration engine when possible
RooArgSet:
- replace Delete() with _list.Delete()
- create addOwned() member fn that allows objects to be added to a list
with ownership transfered to the list (like addClone())
RooFitResult:
- change _globalCorr and the elements of _corrMatrix to use "copied"
ArgSet's and let the ArgSet dtor handle the cleanup (ie, remove Delete()'s
from RooFitResult)
Revision
2654 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed Aug 8 23:11:25 2001 UTC (13 years, 5 months ago) by
david
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2228 byte(s)
Diff to
previous 2629
RooAbsIntegrator:
- no longer inherits from TObject
- uses RooAbsFunc interface instead of RooAbsReal
- remove (AbsReal,mode) binding functionality (now performed by
class RooRealAnalytic)
- add counter for number of integrand evaluations
RooIntegrator1D:
- use RooAbsFunc interface instead of RooAbsIntegrator::eval()
- fix array size bug (already fixed in RFT after this routine was copied)
- add a constructor option to select the summation and refinement rule
to use and change the maxSteps, eps defaults to zero (defaults depend on
the rule being used and are now applied in the initialize() method)
- implement a new addMidpoint() summation rule that can be used to
estimate certain improper integrals
- add check for finite integration limits
RooRealIntegrator:
- add isValid() member fn to flag ctor problems
- add check that all integration vars are lvalues
- use new RooIntegrator1D interface
- use new RooRealBinding,RooRealAnalytic classes
- remove maxSteps, eps ctor parameters (they were not being used)
- use new RooInvTransform and RooImproperIntegrator classes for 1-dim
numerical integrals with one or both limits infinite
- change logic for calling initNumIntegrator() so that it correctly
propagates shape/value changes
- make initNumIntegrator() const so it can be called from evaluate()
and make the data members it modifies mutable
- support 1D numerical integration with either limit infinite via
RooImproperIntegrator1D
RooAbsArg,RooAbsRealLValue,RooAbsCategoryLValue:
- add isLValue() member fns
RooInvTransform:
- new RooAbsFunc implementation that performs the change-of-vars
transformation x -> 1/x on an input RooAbsFunc.
RooImproperIntegrator1D:
- new RooAbsIntegrator implementation that performs an integration
where one or both limits are infinite by splitting up the range
into up to 3 subranges that can each be handled by RooIntegrator1D.
RooNumber:
- new class to encapsulate platform-specific numerics (eg, infinity)
RooAbsRealLValue,RooIntegrator1D,RooImproperIntegrator1D,RooLinearVar,
RooRealIntegrator,RooRealVar:
- use RooNumber statics instead of (non-portable) preprocessor INFINITY
and remove unecessary include dependencies on RooAbsRealLValue.hh
RooIntegratorTest.cc:
- test macro that is not compiled into the library
Revision
2629 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Aug 2 23:54:24 2001 UTC (13 years, 5 months ago) by
david
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2148 byte(s)
Diff to
previous 2215
RooRealIntegral:
- rename Unity to PassThrough
- change function arg type from RooAbsPdf to RooAbsReal
- change (RooAbsReal&) to const_cast (why does the RooRealProxy ctor
need a non-const ref?)
RooAbsIntegrator,RooIntegrator1D:
- change signatures from AbsPdf to AbsReal
Revision
2215 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed May 16 07:41:08 2001 UTC (13 years, 8 months ago) by
verkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2147 byte(s)
Diff to
previous 2200
o RooAbsReal/RooAbsPdf
- Move analytical integration support from AbsReal to AbsPdf
(getAnalyticalIntegral/analyticalIntegral/tryIntegral)
- Add more tryIntegral functions that take 2,3,4 RooArgProxies
o RooAbsRealLValue
- Rename isSafeForIntegration() to isJacobianOK()
- Add jacobian() member that return actual jacobian value
o RooRealIntegral
- Change analytical/numerical integration logic to handle
RooAbsRealLValue objects with a non-unit but constant
jacobian correctly
o RooLinearVar
- Implements most general case of RooAbsRealLValue with
constant but non-unit jacobian: X*slope + offset.
X must be a RooRealVar and slope and offset (RooAbsReal)
must not depend on X
isJacobianOK() only returns true for a given dependent
list of slope (the jacobian) has no dependents, except
for RooAbsCategoryLValues
o Some miscellaneous small bug fixes
Revision
2158 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu May 10 00:16:09 2001 UTC (13 years, 8 months ago) by
verkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2142 byte(s)
Diff to
previous 2033
o RooAbsArg
- Small change in 'Standard' printing logic
- isXXXDirty() will always return kFALSE if isDerived is kFALSE
o RooAbsCategory
- Absorb caching/evaluation code of RooDerivedCategory
- Add defineType(const char*) (without index code)
- Add clearTypes() ;
o RooDerivedCategory
- Retired
o RooAbsCategoryLValue
- New abstract base class for LValue categories
(e.g. RooCategory, RooSuperCategory)
o RooCategory:
- Add methods that publicize the defineType/clearType
members of RooAbsCategory
- Implement getIndex()/Label() that bypasses caching mechanism
o RooMappedCategory
- Complete rewrite using wildcard regular expression mapping
of labels. Hard-coded number mapping deprecated
o RooMapCatEntry
- New auxiliary class for RooMappedCategory
o RooSuperCategory
- New, functionality similar to version in RooFitTools,
but rewritten from scratch
o RooAbsReal
- Absorb caching/evaluation code of RooDerivedReal
- Absorb analytical integration support code of RooDerivedReal
- Absorb range checking code from RooRealVar
o RooAbsRealLValue
- New abstract base class for LValue reals
(.e.g RooRealVar, RooUnblind*)
o RooDerivedReal
- Retired
o RooRealVar
- Implement concrete getVal() that bypasses caching mechanism
o RooProdPdf
- New, functionality similar to version in RooFitTools,
but rewritten from scratch
Revision
2004 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Added
Thu Apr 19 01:42:27 2001 UTC (13 years, 9 months ago) by
verkerke
Original Path:
trunk/roofitcore/inc/RooIntegrator1D.h
File length: 2079 byte(s)
o RooDerivedReal:
- Extend analytical integral interface so that function can,
given a list of reals and categories, indicate which subset it
can analytically integrate / efficiency sum.
o RooRealIntegral
- Can now perform combined analytical/numerical integration.
The RooDerivedReal is asked for the largest part it can
analytically integrate/sum, the rest is taken care of numerically.
The actual numerical integration algorithm is migrated into
a dedicated class.
RRI now also handles higher dimensional numerical integrals,
provided there is a integrator class available for that
o RooAbsIntegrator (new)
Abstract base class for numerical integrators
o RooIntegrator1D (new)
Implementation of 1-dimensional numerical integrator that
used to be implemented in RooRealIntegral
This form allows you to request diffs between any two revisions of this file.
For each of the two "sides" of the diff,
enter a numeric revision.
