RooAbsIntegrator.cxx

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitCore                                                       *
 * @(#)root/roofitcore:$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)             *
 *****************************************************************************/
 
/**
\file RooAbsIntegrator.cxx
\class RooAbsIntegrator
\ingroup Roofitcore
 
Abstract interface for integrators of real-valued
functions that implement the RooAbsFunc interface.
**/
 
#include "Riostream.h"
 
#include "RooAbsIntegrator.h"
#include "RooMsgService.h"
#include "TClass.h"
 
using std::endl;
 
RooAbsIntegrator::RooAbsIntegrator() = default;
 
////////////////////////////////////////////////////////////////////////////////
/// Copy constructor
 
RooAbsIntegrator::RooAbsIntegrator(const RooAbsFunc& function, bool doPrintEvalCounter) :
  _function(&function), _valid(function.isValid()), _printEvalCounter(doPrintEvalCounter)
{
}
 
 
 
////////////////////////////////////////////////////////////////////////////////
/// Calculate integral value with given array of parameter values
 
double RooAbsIntegrator::calculate(const double *yvec)
{
  integrand()->resetNumCall() ;
 
  integrand()->saveXVec() ;
  double ret = integral(yvec) ;
  integrand()->restoreXVec() ;
 
  oocxcoutD(static_cast<TObject*>(nullptr), NumericIntegration) << "RooAbsIntegrator::calculate(" << _function->getName() << ") number of function calls = " << integrand()->numCall()<<", result  = "<<ret << std::endl ;
  return ret ;
}
 
 
 
////////////////////////////////////////////////////////////////////////////////
/// Interface to set limits on integration
 
bool RooAbsIntegrator::setLimits(double xmin, double xmax)
{
  return setLimits(&xmin,&xmax) ;
}
 
 
 
////////////////////////////////////////////////////////////////////////////////
/// Interface function that allows to defer limit definition to integrand definition
 
bool RooAbsIntegrator::setUseIntegrandLimits(bool)
{
  return false ;
}