class RooGaussKronrodIntegrator1D: public RooAbsIntegrator

RooGaussKronrodIntegrator1D implements the Gauss-Kronrod integration algorithm. An Gaussian quadrature method for numerical integration in which error is estimation based on evaluation at special points known as "Kronrod points." By suitably picking these points, abscissas from previous iterations can be reused as part of the new set of points, whereas usual Gaussian quadrature would require recomputation of all abscissas at each iteration. This class automatically handles (-inf,+inf) integrals by dividing the integration in three regions (-inf,-1), (-1,1), (1,inf) and calculating the 1st and 3rd term using a x -> 1/x coordinate transformation This class embeds the Gauss-Kronrod integrator from the GNU Scientific Library version 1.5 and applies the 10-, 21-, 43- and 87-point rule in succession until the required target precision is reached

Function Members (Methods)

public:

	RooGaussKronrodIntegrator1D()
	RooGaussKronrodIntegrator1D(const RooGaussKronrodIntegrator1D&)
	RooGaussKronrodIntegrator1D(const RooAbsFunc& function, const RooNumIntConfig& config)
	RooGaussKronrodIntegrator1D(const RooAbsFunc& function, Double_t xmin, Double_t xmax, const RooNumIntConfig& config)
virtual	~RooGaussKronrodIntegrator1D()
void	TObject::AbstractMethod(const char* method) const
virtual void	TObject::AppendPad(Option_t* option = "")
virtual void	TObject::Browse(TBrowser* b)
Double_t	RooAbsIntegrator::calculate(const Double_t* yvec = 0)
virtual Bool_t	canIntegrate1D() const
virtual Bool_t	canIntegrate2D() const
virtual Bool_t	canIntegrateND() const
virtual Bool_t	canIntegrateOpenEnded() const
virtual Bool_t	checkLimits() const
static TClass*	Class()
virtual const char*	TObject::ClassName() const
virtual void	TObject::Clear(Option_t* = "")
virtual TObject*	TObject::Clone(const char* newname = "") const
virtual RooAbsIntegrator*	clone(const RooAbsFunc& function, const RooNumIntConfig& config) const
virtual Int_t	TObject::Compare(const TObject* obj) const
virtual void	TObject::Copy(TObject& object) const
virtual void	TObject::Delete(Option_t* option = "")MENU
virtual Int_t	TObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual void	TObject::Draw(Option_t* option = "")
virtual void	TObject::DrawClass() constMENU
virtual TObject*	TObject::DrawClone(Option_t* option = "") constMENU
virtual void	TObject::Dump() constMENU
virtual void	TObject::Error(const char* method, const char* msgfmt) const
virtual void	TObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual void	TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual void	TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual void	TObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*	TObject::FindObject(const char* name) const
virtual TObject*	TObject::FindObject(const TObject* obj) const
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
virtual const char*	TObject::GetIconName() const
virtual const char*	TObject::GetName() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
virtual const char*	TObject::GetTitle() const
virtual UInt_t	TObject::GetUniqueID() const
virtual Bool_t	TObject::HandleTimer(TTimer* timer)
virtual ULong_t	TObject::Hash() const
virtual void	TObject::Info(const char* method, const char* msgfmt) const
virtual Bool_t	TObject::InheritsFrom(const char* classname) const
virtual Bool_t	TObject::InheritsFrom(const TClass* cl) const
virtual void	TObject::Inspect() constMENU
virtual Double_t	integral(const Double_t* yvec = 0)
const RooAbsFunc*	RooAbsIntegrator::integrand() const
Double_t	RooAbsIntegrator::integrand(const Double_t* x) const
void	TObject::InvertBit(UInt_t f)
virtual TClass*	IsA() const
virtual Bool_t	TObject::IsEqual(const TObject* obj) const
virtual Bool_t	TObject::IsFolder() const
Bool_t	TObject::IsOnHeap() const
virtual Bool_t	TObject::IsSortable() const
Bool_t	RooAbsIntegrator::isValid() const
Bool_t	TObject::IsZombie() const
virtual void	TObject::ls(Option_t* option = "") const
void	TObject::MayNotUse(const char* method) const
virtual Bool_t	TObject::Notify()
void	TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static void	TObject::operator delete(void* ptr)
static void	TObject::operator delete(void* ptr, void* vp)
static void	TObject::operator delete[](void* ptr)
static void	TObject::operator delete[](void* ptr, void* vp)
void*	TObject::operator new(size_t sz)
void*	TObject::operator new(size_t sz, void* vp)
void*	TObject::operator new[](size_t sz)
void*	TObject::operator new[](size_t sz, void* vp)
TObject&	TObject::operator=(const TObject& rhs)
virtual void	TObject::Paint(Option_t* option = "")
virtual void	TObject::Pop()
virtual void	TObject::Print(Option_t* option = "") const
Bool_t	RooAbsIntegrator::printEvalCounter() const
virtual Int_t	TObject::Read(const char* name)
virtual void	TObject::RecursiveRemove(TObject* obj)
void	TObject::ResetBit(UInt_t f)
virtual void	TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual void	TObject::SavePrimitive(ostream& out, Option_t* option = "")
void	TObject::SetBit(UInt_t f)
void	TObject::SetBit(UInt_t f, Bool_t set)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
virtual Bool_t	setLimits(Double_t* xmin, Double_t* xmax)
static void	TObject::SetObjectStat(Bool_t stat)
void	RooAbsIntegrator::setPrintEvalCounter(Bool_t value)
virtual void	TObject::SetUniqueID(UInt_t uid)
virtual Bool_t	setUseIntegrandLimits(Bool_t flag)
virtual void	ShowMembers(TMemberInspector& insp)
virtual void	Streamer(TBuffer& b)
void	StreamerNVirtual(TBuffer& b)
virtual void	TObject::SysError(const char* method, const char* msgfmt) const
Bool_t	TObject::TestBit(UInt_t f) const
Int_t	TObject::TestBits(UInt_t f) const
virtual void	TObject::UseCurrentStyle()
virtual void	TObject::Warning(const char* method, const char* msgfmt) const
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const

protected:

virtual void	TObject::DoError(int level, const char* location, const char* fmt, va_list va) const
Bool_t	initialize()
void	TObject::MakeZombie()
static void	registerIntegrator(RooNumIntFactory& fact)
Double_t*	xvec(Double_t& xx)

Data Members

public:

enum TObject::EStatusBits {	kCanDelete
	kMustCleanup
	kObjInCanvas
	kIsReferenced
	kHasUUID
	kCannotPick
	kNoContextMenu
	kInvalidObject
};
enum TObject::[unnamed] {	kIsOnHeap
	kNotDeleted
	kZombie
	kBitMask
	kSingleKey
	kOverwrite
	kWriteDelete
};

protected:

Double_t	_epsAbs	Absolute precision
Double_t	_epsRel	Relative precision
const RooAbsFunc*	RooAbsIntegrator::_function	Pointer to function binding of integrand
Bool_t	RooAbsIntegrator::_printEvalCounter	If true print number of function evaluation required for integration
Bool_t	_useIntegrandLimits	Use limits in function binding?
Bool_t	RooAbsIntegrator::_valid	Is integrator in valid state?
Double_t*	_x	! do not persist
Double_t	_xmax	! Upper integration bound
Double_t	_xmin	! Lower integration bound

Class Charts

Inheritance Inherited Members Includes Libraries

Function documentation

void registerIntegrator(RooNumIntFactory& fact)

 Register RooGaussKronrodIntegrator1D, its parameters and capabilities with RooNumIntConfig

RooGaussKronrodIntegrator1D()

 coverity[UNINIT_CTOR]
 Default constructor

RooGaussKronrodIntegrator1D(const RooAbsFunc& function, const RooNumIntConfig& config)

 Construct integral on 'function' using given configuration object. The integration
 range is taken from the definition in the function binding

RooGaussKronrodIntegrator1D(const RooAbsFunc& function, Double_t xmin, Double_t xmax, const RooNumIntConfig& config)

 Construct integral on 'function' using given configuration object in the given range

RooAbsIntegrator* clone(const RooAbsFunc& function, const RooNumIntConfig& config) const

 Clone integrator with given function and configuration. Needed for RooNumIntFactory

Bool_t initialize()

 Perform one-time initialization of integrator

~RooGaussKronrodIntegrator1D()

 Destructor

Bool_t setLimits(Double_t* xmin, Double_t* xmax)

 Change our integration limits. Return kTRUE if the new limits are
 ok, or otherwise kFALSE. Always returns kFALSE and does nothing
 if this object was constructed to always use our integrand's limits.

Bool_t checkLimits() const

 Check that our integration range is finite and otherwise return kFALSE.
 Update the limits from the integrand if requested.

Double_t integral(const Double_t* yvec = 0)

 Calculate and return integral

RooGaussKronrodIntegrator1D()

 Constructors, assignment etc

Bool_t setUseIntegrandLimits(Bool_t flag)

{_useIntegrandLimits = flag ; return kTRUE ; }

Bool_t canIntegrate1D() const

{ return kTRUE ; }

Bool_t canIntegrate2D() const

{ return kFALSE ; }

Bool_t canIntegrateND() const

{ return kFALSE ; }

Bool_t canIntegrateOpenEnded() const

{ return kTRUE ; }

Double_t* xvec(Double_t& xx)

{ _x[0] = xx ; return _x ; }