/usr/include/root/RooAdaptiveGaussKronrodIntegrator1D.h is in libroot-roofit-dev 5.34.00-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 | /*****************************************************************************
* Project: RooFit *
* Package: RooFitCore *
* File: $Id: RooAdaptiveGaussKronrodIntegrator1D.h,v 1.5 2007/05/11 09:11:30 verkerke Exp $
* 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) *
*****************************************************************************/
#ifndef ROO_ADAPTIVE_GAUSS_KRONROD_INTEGRATOR_1D
#define ROO_ADAPTIVE_GAUSS_KRONROD_INTEGRATOR_1D
#include "RooAbsIntegrator.h"
#include "RooNumIntConfig.h"
double RooAdaptiveGaussKronrodIntegrator1D_GSL_GlueFunction(double x, void *data) ;
class RooAdaptiveGaussKronrodIntegrator1D : public RooAbsIntegrator {
public:
// Constructors, assignment etc
RooAdaptiveGaussKronrodIntegrator1D() ;
RooAdaptiveGaussKronrodIntegrator1D(const RooAbsFunc& function, const RooNumIntConfig& config) ;
RooAdaptiveGaussKronrodIntegrator1D(const RooAbsFunc& function, Double_t xmin, Double_t xmax,
const RooNumIntConfig& config) ;
virtual RooAbsIntegrator* clone(const RooAbsFunc& function, const RooNumIntConfig& config) const ;
virtual ~RooAdaptiveGaussKronrodIntegrator1D();
virtual Bool_t checkLimits() const;
virtual Double_t integral(const Double_t *yvec=0) ;
using RooAbsIntegrator::setLimits ;
Bool_t setLimits(Double_t* xmin, Double_t* xmax);
virtual Bool_t setUseIntegrandLimits(Bool_t flag) {
// If flag is true, intergration limits are taken from definition in input function binding
_useIntegrandLimits = flag ; return kTRUE ;
}
virtual Bool_t canIntegrate1D() const {
// We can integrate 1-dimensional functions
return kTRUE ;
}
virtual Bool_t canIntegrate2D() const {
// We can not integrate 2-dimensional functions
return kFALSE ;
}
virtual Bool_t canIntegrateND() const {
// We can not integrate >2-dimensional functions
return kFALSE ;
}
virtual Bool_t canIntegrateOpenEnded() const {
// We can integrate over open-ended domains
return kTRUE ;
}
protected:
friend class RooNumIntFactory ;
static void registerIntegrator(RooNumIntFactory& fact) ;
enum DomainType { Closed, OpenLo, OpenHi, Open } ;
mutable DomainType _domainType ;
friend double RooAdaptiveGaussKronrodIntegrator1D_GSL_GlueFunction(double x, void *data) ;
Bool_t initialize();
Bool_t _useIntegrandLimits;
Double_t* xvec(Double_t& xx) {
// Return contents of xx in internal array pointer
_x[0] = xx ; return _x ;
}
Double_t *_x ; //! Current coordinate
Double_t _epsAbs ; // Absolute precision
Double_t _epsRel ; // Relative precision
Int_t _methodKey ; // GSL method key
Int_t _maxSeg ; // Maximum number of segments
void* _workspace ; // GSL workspace
mutable Double_t _xmin; //! Lower integration bound
mutable Double_t _xmax; //! Upper integration bound
ClassDef(RooAdaptiveGaussKronrodIntegrator1D,0) // 1-dimensional adaptive Gauss-Kronrod numerical integration engine
};
#endif
|