#include "RooFit.h"
#include "RooNumIntConfig.h"
#include "RooNumIntConfig.h"
#include "RooArgSet.h"
#include "RooAbsIntegrator.h"
#include "RooNumIntFactory.h"
#include "TClass.h"
ClassImp(RooNumIntConfig)
;
RooNumIntConfig& RooNumIntConfig::defaultConfig()
{
static RooNumIntConfig* _default = 0 ;
if (_default==0) {
_default = new RooNumIntConfig ;
RooNumIntFactory::instance() ;
}
return *_default ;
}
RooNumIntConfig::RooNumIntConfig() :
_epsAbs(1e-7),
_epsRel(1e-7),
_printEvalCounter(kFALSE),
_method1D("method1D","1D integration method"),
_method2D("method2D","2D integration method"),
_methodND("methodND","ND integration method"),
_method1DOpen("method1DOpen","1D integration method in open domain"),
_method2DOpen("method2DOpen","2D integration method in open domain"),
_methodNDOpen("methodNDOpen","ND integration method in open domain")
{
_method1D.defineType("N/A",0) ;
_method2D.defineType("N/A",0) ;
_methodND.defineType("N/A",0) ;
_method1DOpen.defineType("N/A",0) ;
_method2DOpen.defineType("N/A",0) ;
_methodNDOpen.defineType("N/A",0) ;
}
RooNumIntConfig::~RooNumIntConfig()
{
_configSets.Delete() ;
}
RooNumIntConfig::RooNumIntConfig(const RooNumIntConfig& other) :
TObject(other), RooPrintable(other),
_epsAbs(other._epsAbs),
_epsRel(other._epsRel),
_printEvalCounter(other._printEvalCounter),
_method1D(other._method1D),
_method2D(other._method2D),
_methodND(other._methodND),
_method1DOpen(other._method1DOpen),
_method2DOpen(other._method2DOpen),
_methodNDOpen(other._methodNDOpen)
{
TIterator* iter = other._configSets.MakeIterator() ;
RooArgSet* set ;
while((set=(RooArgSet*)iter->Next())) {
RooArgSet* setCopy = (RooArgSet*) set->snapshot() ;
setCopy->setName(set->GetName()) ;
_configSets.Add(setCopy);
}
delete iter ;
}
RooNumIntConfig& RooNumIntConfig::operator=(const RooNumIntConfig& other)
{
if (&other==this) {
return *this ;
}
_epsAbs = other._epsAbs ;
_epsRel = other._epsRel ;
_method1D.setIndex(other._method1D.getIndex()) ;
_method2D.setIndex(other._method2D.getIndex()) ;
_methodND.setIndex(other._methodND.getIndex()) ;
_method1DOpen.setIndex(other._method1DOpen.getIndex()) ;
_method2DOpen.setIndex(other._method2DOpen.getIndex()) ;
_methodNDOpen.setIndex(other._methodNDOpen.getIndex()) ;
_configSets.Delete() ;
TIterator* iter = other._configSets.MakeIterator() ;
RooArgSet* set ;
while((set=(RooArgSet*)iter->Next())) {
RooArgSet* setCopy = (RooArgSet*) set->snapshot() ;
setCopy->setName(set->GetName()) ;
_configSets.Add(setCopy);
}
delete iter ;
return *this ;
}
Bool_t RooNumIntConfig::addConfigSection(const RooAbsIntegrator* proto, const RooArgSet& defaultConfig)
{
TString name = proto->IsA()->GetName() ;
if (proto->canIntegrate1D()) {
_method1D.defineType(name) ;
if (proto->canIntegrateOpenEnded()) {
_method1DOpen.defineType(name) ;
}
}
if (proto->canIntegrate2D()) {
_method2D.defineType(name) ;
if (proto->canIntegrateOpenEnded()) {
_method2DOpen.defineType(name) ;
}
}
if (proto->canIntegrateND()) {
_methodND.defineType(name) ;
if (proto->canIntegrateOpenEnded()) {
_methodNDOpen.defineType(name) ;
}
}
RooArgSet* config = (RooArgSet*) defaultConfig.snapshot() ;
config->setName(name) ;
_configSets.Add(config) ;
return kFALSE ;
}
RooArgSet& RooNumIntConfig::getConfigSection(const char* name)
{
return const_cast<RooArgSet&>((const_cast<RooNumIntConfig*>(this)->getConfigSection(name))) ;
}
const RooArgSet& RooNumIntConfig::getConfigSection(const char* name) const
{
static RooArgSet dummy ;
RooArgSet* config = (RooArgSet*) _configSets.FindObject(name) ;
if (!config) {
cout << "RooNumIntConfig::getIntegrator: ERROR: no configuration stored for integrator '" << name << "'" << endl ;
return dummy ;
}
return *config ;
}
void RooNumIntConfig::setEpsAbs(Double_t newEpsAbs)
{
if (newEpsAbs<=0) {
cout << "RooNumIntConfig::setEpsAbs: ERROR: target absolute precision must be greater than zero" << endl ;
return ;
}
_epsAbs = newEpsAbs ;
}
void RooNumIntConfig::setEpsRel(Double_t newEpsRel)
{
if (newEpsRel<=0) {
cout << "RooNumIntConfig::setEpsRel: ERROR: target absolute precision must be greater than zero" << endl ;
return ;
}
_epsRel = newEpsRel ;
}
void RooNumIntConfig::printToStream(ostream &os, PrintOption opt, TString indent) const
{
switch(opt) {
case InLine:
os << indent << "RooNumIntConfig" ;
case OneLine:
os << indent << endl ;
break ;
case Shape:
break ;
case Standard:
os << indent << "Requested precision: " << _epsAbs << " absolute, " << _epsRel << " relative" << endl << endl ;
if (_printEvalCounter) {
os << indent << "Printing of function evaluation counter for each integration enabled" << endl << endl ;
}
os << indent << "1-D integration method: " << _method1D.getLabel() ;
if (_method1DOpen.getIndex()!=_method1D.getIndex()) {
os << " (" << _method1DOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
os << indent << "2-D integration method: " << _method2D.getLabel() ;
if (_method2DOpen.getIndex()!=_method2D.getIndex()) {
os << " (" << _method2DOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
os << indent << "N-D integration method: " << _methodND.getLabel() ;
if (_methodNDOpen.getIndex()!=_methodND.getIndex()) {
os << " (" << _methodNDOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
break ;
case Verbose:
os << indent << "Requested precision: " << _epsAbs << " absolute, " << _epsRel << " relative" << endl << endl ;;
if (_printEvalCounter) {
os << indent << "Printing of function evaluation counter for each integration enabled" << endl << endl ;
}
os << indent << "Selected integration methods:" << endl ;
os << indent << "1-D integration method: " << _method1D.getLabel() ;
if (_method1DOpen.getIndex()!=_method1D.getIndex()) {
os << " (" << _method1DOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
os << indent << "2-D integration method: " << _method2D.getLabel() ;
if (_method2DOpen.getIndex()!=_method2D.getIndex()) {
os << " (" << _method2DOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
os << indent << "N-D integration method: " << _methodND.getLabel() ;
if (_methodNDOpen.getIndex()!=_methodND.getIndex()) {
os << " (" << _methodNDOpen.getLabel() << " if open-ended)" << endl ;
} else {
os << endl ;
}
os << endl << "Available integration methods:" << endl << endl ;
TIterator* cIter = _configSets.MakeIterator() ;
RooArgSet* configSet ;
while ((configSet=(RooArgSet*)cIter->Next())) {
os << indent << "*** " << configSet->GetName() << " ***" << endl ;
os << indent << "Capabilities: " ;
const RooAbsIntegrator* proto = RooNumIntFactory::instance().getProtoIntegrator(configSet->GetName()) ;
if (proto->canIntegrate1D()) os << "[1-D] " ;
if (proto->canIntegrate2D()) os << "[2-D] " ;
if (proto->canIntegrateND()) os << "[N-D] " ;
if (proto->canIntegrateOpenEnded()) os << "[OpenEnded] " ;
os << endl ;
os << "Configuration: " << endl ;
configSet->writeToStream(os,kFALSE) ;
const char* depName = RooNumIntFactory::instance().getDepIntegratorName(configSet->GetName()) ;
if (strlen(depName)>0) {
os << indent << "(Depends on '" << depName << "')" << endl ;
}
os << endl ;
}
delete cIter ;
break ;
}
}
ROOT page - Class index - Class Hierarchy - Top of the page
This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.