#include <MELANCSpline_1D_fast.h>

Inheritance diagram for MELANCSpline_1D_fast:

Collaboration diagram for MELANCSpline_1D_fast:

Public Member Functions
	MELANCSpline_1D_fast ()

	MELANCSpline_1D_fast (const char name, const char title)

	MELANCSpline_1D_fast (const char name, const char title, RooAbsReal &inXVar, const std::vector< T > &inXList, const std::vector< T > &inFcnList, MELANCSplineCore::BoundaryCondition const bcBeginX_=MELANCSplineCore::bcNaturalSpline, MELANCSplineCore::BoundaryCondition const bcEndX_=MELANCSplineCore::bcNaturalSpline, Bool_t inUseFloor=true, T inFloorEval=0, T inFloorInt=0)

	MELANCSpline_1D_fast (const MELANCSpline_1D_fast &other, const char *name=0)

virtual TObject *	clone (const char *newname) const

virtual	~MELANCSpline_1D_fast ()

void	setRangeValidity (const T valmin, const T valmax, const Int_t whichDirection=0)

virtual Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const

virtual Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const

Public Member Functions inherited from MELANCSplineCore
	MELANCSplineCore ()

	MELANCSplineCore (const char name, const char title)

	MELANCSplineCore (const char name, const char title, RooAbsReal &inXVar, const std::vector< T > &inXList, Bool_t inUseFloor=true, T inFloorEval=1e-15, T inFloorInt=1e-10)

	MELANCSplineCore (const MELANCSplineCore &other, const char *name=0)

virtual	~MELANCSplineCore ()

virtual void	setVerbosity (VerbosityLevel flag)

void	setEvalFloor (T val)

void	setIntFloor (T val)

void	doFloor (Bool_t flag)

Protected Member Functions
virtual void	emptyFcnList ()

virtual Int_t	getWhichBin (const T &val, const Int_t whichDirection) const

virtual T	getTVar (const std::vector< T > &kappas, const T &val, const Int_t &bin, const Int_t whichDirection) const

virtual void	getKappas (std::vector< T > &kappas, const Int_t whichDirection)

Bool_t	testRangeValidity (const T &val, const Int_t whichDirection=0) const

void	cropValueForRange (T &val, const Int_t whichDirection=0) const

virtual T	interpolateFcn (Int_t code, const char *rangeName=0) const

virtual Double_t	evaluate () const

Protected Member Functions inherited from MELANCSplineCore
void	getLeafDependents (RooRealProxy &proxy, RooArgSet &set)

void	addLeafDependents (RooArgSet &set)

unsigned int	npointsX () const

virtual void	getBArray (const std::vector< T > &kappas, const std::vector< T > &fcnList, std::vector< T > &BArray, BoundaryCondition const &bcBegin, BoundaryCondition const &bcEnd) const

virtual void	getAArray (const std::vector< T > &kappas, std::vector< std::vector< T >> &AArray, BoundaryCondition const &bcBegin, BoundaryCondition const &bcEnd) const

virtual std::vector< std::vector< T > >	getCoefficientsAlongDirection (const std::vector< T > &kappas, const TMatrix_t &Ainv, const std::vector< T > &fcnList, BoundaryCondition const &bcBegin, BoundaryCondition const &bcEnd, const Int_t pickBin) const

virtual std::vector< T >	getCoefficients (const TVector_t &S, const std::vector< T > &kappas, const std::vector< T > &fcnList, const Int_t &bin) const

virtual T	evalSplineSegment (const std::vector< T > &coefs, const T &kappa, const T &tup, const T &tdn, Bool_t doIntegrate=false) const

Protected Attributes
const BoundaryCondition	bcBeginX

const BoundaryCondition	bcEndX

std::vector< T >	FcnList

std::vector< T >	kappaX

std::vector< std::vector< T > >	coefficients

Protected Attributes inherited from MELANCSplineCore
VerbosityLevel	verbosity

Bool_t	useFloor

T	floorEval

T	floorInt

T	rangeXmin

T	rangeXmax

RooRealProxy	theXVar

RooListProxy	leafDepsList

std::vector< T >	XList

Additional Inherited Members
Public Types inherited from MELANCSplineCore
enum	VerbosityLevel { kSilent, kError, kVerbose }

enum	BoundaryCondition { bcApproximatedSlope, bcClamped, bcApproximatedSecondDerivative, bcNaturalSpline, bcQuadratic, bcQuadraticWithNullSlope, NBoundaryConditions }

typedef Float_t	T

typedef TMatrixT< T >	TMatrix_t

typedef TVectorT< T >	TVector_t

Detailed Description

Definition at line 12 of file MELANCSpline_1D_fast.h.

Constructor & Destructor Documentation

◆ MELANCSpline_1D_fast() [1/4]

MELANCSpline_1D_fast::MELANCSpline_1D_fast ( )

◆ MELANCSpline_1D_fast() [2/4]

MELANCSpline_1D_fast::MELANCSpline_1D_fast	(	const char *	name,
		const char *	title
	)

Definition at line 20 of file MELANCSpline_1D_fast.cc.

     :
   MELANCSplineCore(name, title),
   bcBeginX(MELANCSplineCore::bcNaturalSpline), bcEndX(MELANCSplineCore::bcNaturalSpline)
 {}

◆ MELANCSpline_1D_fast() [3/4]

MELANCSpline_1D_fast::MELANCSpline_1D_fast	(	const char *	name,
		const char *	title,
		RooAbsReal &	inXVar,
		const std::vector< T > &	inXList,
		const std::vector< T > &	inFcnList,
		MELANCSplineCore::BoundaryCondition const	bcBeginX_ = `MELANCSplineCore::bcNaturalSpline`,
		MELANCSplineCore::BoundaryCondition const	bcEndX_ = `MELANCSplineCore::bcNaturalSpline`,
		Bool_t	inUseFloor = `true`,
		T	inFloorEval = `0`,
		T	inFloorInt = `0`
	)

Definition at line 28 of file MELANCSpline_1D_fast.cc.

     :
   MELANCSplineCore(name, title, inXVar, inXList, inUseFloor, inFloorEval, inFloorInt),
   bcBeginX(bcBeginX_), bcEndX(bcEndX_),
   FcnList(inFcnList)
 {
   if (npointsX()>1){
     int npoints;
  
     vector<vector<MELANCSplineCore::T>> xA; getKappas(kappaX, 0); getAArray(kappaX, xA, bcBeginX, bcEndX);
     npoints=kappaX.size();
     TMatrix_t xAtrans(npoints, npoints);
     for (int i=0; i<npoints; i++){ for (int j=0; j<npoints; j++){ xAtrans[i][j]=xA.at(i).at(j); } }
     Double_t det=0;
     TMatrix_t xAinv = xAtrans.Invert(&det);
     if (det==0.){
       coutE(InputArguments) << "MELANCSpline_1D::interpolateFcn: Matrix xA could not be inverted. Something is wrong with the x coordinates of points!" << endl;
       assert(0);
     }
  
     coefficients = getCoefficientsAlongDirection(kappaX, xAinv, FcnList, bcBeginX, bcEndX, -1);
   }
   else assert(0);
  
   RooArgSet leafset;
   getLeafDependents(theXVar, leafset);
   addLeafDependents(leafset);
  
   emptyFcnList();
 }

◆ MELANCSpline_1D_fast() [4/4]

MELANCSpline_1D_fast::MELANCSpline_1D_fast	(	const MELANCSpline_1D_fast &	other,
		const char *	name = `0`
	)

Definition at line 69 of file MELANCSpline_1D_fast.cc.

     :
   MELANCSplineCore(other, name),
   bcBeginX(other.bcBeginX), bcEndX(other.bcEndX),
   FcnList(other.FcnList),
   kappaX(other.kappaX),
   coefficients(other.coefficients)
 {}

◆ ~MELANCSpline_1D_fast()

virtual MELANCSpline_1D_fast::~MELANCSpline_1D_fast ( )

inlinevirtual

Definition at line 42 of file MELANCSpline_1D_fast.h.

42 {}

Member Function Documentation

◆ analyticalIntegral()

Double_t MELANCSpline_1D_fast::analyticalIntegral	(	Int_t	code,
		const char *	rangeName = `0`
	)		const

virtual

Implements MELANCSplineCore.

Definition at line 197 of file MELANCSpline_1D_fast.cc.

                                                                                         {
   Double_t value = interpolateFcn(code, rangeName);
   if (useFloor && value<floorInt){
     if (verbosity>=MELANCSplineCore::kError) coutE(Integration) << "MELANCSpline_1D_fast ERROR::MELANCSpline_1D_fast(" << GetName() << ") integration returned " << value << " for code = " << code << endl;
     value = floorInt;
   }
   if (verbosity==MELANCSplineCore::kVerbose){ cout << "MELANCSpline_1D_fast(" << GetName() << ")::analyticalIntegral = " << value << " for code = " << code << endl; }
   return value;
 }

◆ clone()

virtual TObject* MELANCSpline_1D_fast::clone ( const char * newname ) const

inlinevirtual

Implements MELANCSplineCore.

Definition at line 41 of file MELANCSpline_1D_fast.h.

41 { return new MELANCSpline_1D_fast(*this, newname); }

◆ cropValueForRange()

void MELANCSpline_1D_fast::cropValueForRange	(	T &	val,
		const Int_t	whichDirection = `0`
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 220 of file MELANCSpline_1D_fast.cc.

                                                                       {
   if (testRangeValidity(val)) return;
   const T* range[2];
   range[0] = &rangeXmin;
   range[1] = &rangeXmax;
   if (val<*(range[0])) val = *(range[0]);
   if (val>*(range[1])) val = *(range[1]);
 }

◆ emptyFcnList()

virtual void MELANCSpline_1D_fast::emptyFcnList ( )

inlineprotectedvirtual

Implements MELANCSplineCore.

Definition at line 50 of file MELANCSpline_1D_fast.h.

50 { std::vector<T> tmp; FcnList.swap(tmp); }

◆ evaluate()

Double_t MELANCSpline_1D_fast::evaluate ( ) const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 170 of file MELANCSpline_1D_fast.cc.

                                              {
   Double_t value = interpolateFcn(0);
   if (useFloor && value<floorEval){
     if (verbosity>=MELANCSplineCore::kError) coutE(Eval) << "MELANCSpline_1D_fast ERROR::MELANCSpline_1D_fast(" << GetName() << ") evaluation returned " << value << " at x = " << theXVar << endl;
     value = floorEval;
   }
   if (verbosity==MELANCSplineCore::kVerbose){
     cout << "MELANCSpline_1D_fast(" << GetName() << ")::evaluate = " << value << " at x = " << theXVar << endl;
     RooArgSet Xdeps; theXVar.absArg()->leafNodeServerList(&Xdeps, 0, true);
     TIterator* iter = Xdeps.createIterator();
     RooAbsArg* var;
     while ((var = (RooAbsArg*)iter->Next())){
       cout << var->GetName() << " value = " << dynamic_cast<RooAbsReal*>(var)->getVal() << endl;
     }
     delete iter;
     cout << endl;
   }
   return value;
 }

◆ getAnalyticalIntegral()

Int_t MELANCSpline_1D_fast::getAnalyticalIntegral	(	RooArgSet &	allVars,
		RooArgSet &	analVars,
		const char *	rangeName = `0`
	)		const

virtual

Implements MELANCSplineCore.

Definition at line 189 of file MELANCSpline_1D_fast.cc.

                                                                                                              {
   if (_forceNumInt) return 0;
   Int_t code=0;
   if (dynamic_cast<RooRealVar*>(theXVar.absArg())!=0){
     if (matchArgs(allVars, analVars, theXVar)) code=2;
   }
   return code;
 }

◆ getKappas()

void MELANCSpline_1D_fast::getKappas	(	std::vector< T > &	kappas,
		const Int_t	whichDirection
	)

protectedvirtual

Implements MELANCSplineCore.

Definition at line 124 of file MELANCSpline_1D_fast.cc.

                                                                                     {
   kappas.clear();
   MELANCSplineCore::T kappa=1;
  
   Int_t npoints;
   vector<MELANCSplineCore::T> const* coord;
   npoints=npointsX();
   coord=&XList;
  
   for (Int_t j=0; j<npoints-1; j++){
     MELANCSplineCore::T val_j = coord->at(j);
     MELANCSplineCore::T val_jpo = coord->at(j+1);
     MELANCSplineCore::T val_diff = (val_jpo-val_j);
     if (fabs(val_diff)>MELANCSplineCore::T(0)) kappa = 1./val_diff;
     else kappa = 0;
     kappas.push_back(kappa);
   }
   kappas.push_back(kappa); // Push the same kappa_(N-1)=kappa_(N-2) at the end point
 }

◆ getTVar()

MELANCSplineCore::T MELANCSpline_1D_fast::getTVar	(	const std::vector< T > &	kappas,
		const T &	val,
		const Int_t &	bin,
		const Int_t	whichDirection
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 165 of file MELANCSpline_1D_fast.cc.

                                                                                                                                                           {
   const MELANCSplineCore::T& K=kappas.at(bin);
   return (val-XList.at(bin))*K;
 }

◆ getWhichBin()

Int_t MELANCSpline_1D_fast::getWhichBin	(	const T &	val,
		const Int_t	whichDirection
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 143 of file MELANCSpline_1D_fast.cc.

                                                                                        {
   Int_t bin=-1;
   MELANCSplineCore::T valj, valjpo;
   Int_t npoints;
   vector<MELANCSplineCore::T> const* coord;
   coord=&XList;
   npoints=npointsX();
  
   if (npoints<=1) bin=0;
   else{
     valjpo = coord->at(0);
     for (Int_t j=0; j<npoints-1; j++){
       valj = coord->at(j);
       valjpo = coord->at(j+1);
       if (val<valjpo && val>=valj){ bin=j; break; }
     }
     if (bin==-1 && val>=valjpo) bin=npoints-2;
     else if (bin==-1) bin=0;
   }
  
   return bin;
 }

◆ interpolateFcn()

MELANCSplineCore::T MELANCSpline_1D_fast::interpolateFcn	(	Int_t	code,
		const char *	rangeName = `0`
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 80 of file MELANCSpline_1D_fast.cc.

                                                                                             {
   DefaultAccumulator<MELANCSplineCore::T> res;
  
   if (verbosity==MELANCSplineCore::kVerbose) cout << "MELANCSpline_1D_fast(" << GetName() << ")::interpolateFcn begin with code: " << code << endl;
  
   // Get bins
   Int_t xbin=-1, xbinmin=-1, xbinmax=-1;
   MELANCSplineCore::T tx=0, txmin=0, txmax=0;
   if (code==0 || code%2!=0){ // Case to just compute the value at x
     if (!testRangeValidity(theXVar)) return 0;
     xbin = getWhichBin(theXVar, 0);
     tx = getTVar(kappaX, theXVar, xbin, 0);
   }
   else{ // Case to integrate along x
     MELANCSplineCore::T coordmin = theXVar.min(rangeName); cropValueForRange(coordmin);
     MELANCSplineCore::T coordmax = theXVar.max(rangeName); cropValueForRange(coordmax);
     xbinmin = getWhichBin(coordmin, 0);
     txmin = getTVar(kappaX, coordmin, xbinmin, 0);
     xbinmax = getWhichBin(coordmax, 0);
     txmax = getTVar(kappaX, coordmax, xbinmax, 0);
   }
  
   int nxbins = (int)coefficients.size();
   for (int ix=0; ix<nxbins; ix++){
     if (
       (xbin>=0 && ix!=xbin)
       ||
       (xbinmin>=0 && xbinmax>=xbinmin && !(xbinmin<=ix && ix<=xbinmax))
       ) continue;
  
     MELANCSplineCore::T txlow=0, txhigh=1;
     if (code>0 && code%2==0){
       if (ix==xbinmin) txlow=txmin;
       if (ix==xbinmax) txhigh=txmax;
     }
     else txhigh=tx;
  
     // Get the x coefficients at bin ix and evaluate value of spline at x
     res += evalSplineSegment(coefficients.at(ix), kappaX.at(ix), txhigh, txlow, (code>0 && code%2==0));
   }
  
   return res;
 }

◆ setRangeValidity()

void MELANCSpline_1D_fast::setRangeValidity	(	const T	valmin,
		const T	valmax,
		const Int_t	whichDirection = `0`
	)

virtual

Implements MELANCSplineCore.

Definition at line 213 of file MELANCSpline_1D_fast.cc.

                                                                                         {
   T* range[2];
   range[0] = &rangeXmin;
   range[1] = &rangeXmax;
   *(range[0])=valmin;
   *(range[1])=valmax;
 }

◆ testRangeValidity()

Bool_t MELANCSpline_1D_fast::testRangeValidity	(	const T &	val,
		const Int_t	whichDirection = `0`
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 207 of file MELANCSpline_1D_fast.cc.

                                                                                {
   const T* range[2];
   range[0] = &rangeXmin;
   range[1] = &rangeXmax;
   return (*(range[0])>*(range[1]) || (val>=*(range[0]) && val<=*(range[1])));
 }

Member Data Documentation

◆ bcBeginX

const BoundaryCondition MELANCSpline_1D_fast::bcBeginX

protected

Definition at line 14 of file MELANCSpline_1D_fast.h.

◆ bcEndX

const BoundaryCondition MELANCSpline_1D_fast::bcEndX

protected

Definition at line 15 of file MELANCSpline_1D_fast.h.

◆ coefficients

std::vector<std::vector<T> > MELANCSpline_1D_fast::coefficients

protected

Definition at line 20 of file MELANCSpline_1D_fast.h.

◆ FcnList

std::vector<T> MELANCSpline_1D_fast::FcnList

protected

Definition at line 17 of file MELANCSpline_1D_fast.h.

◆ kappaX

std::vector<T> MELANCSpline_1D_fast::kappaX

protected

Definition at line 19 of file MELANCSpline_1D_fast.h.

The documentation for this class was generated from the following files:

MELA/interface/MELANCSpline_1D_fast.h
MELA/src/MELANCSpline_1D_fast.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MELANCSpline_1D_fast() [1/4]

◆ MELANCSpline_1D_fast() [2/4]

◆ MELANCSpline_1D_fast() [3/4]

◆ MELANCSpline_1D_fast() [4/4]

◆ ~MELANCSpline_1D_fast()

Member Function Documentation

◆ analyticalIntegral()

◆ clone()

◆ cropValueForRange()

◆ emptyFcnList()

◆ evaluate()

◆ getAnalyticalIntegral()

◆ getKappas()

◆ getTVar()

◆ getWhichBin()

◆ interpolateFcn()

◆ setRangeValidity()

◆ testRangeValidity()

Member Data Documentation

◆ bcBeginX

◆ bcEndX

◆ coefficients

◆ FcnList

◆ kappaX