#include <MELANCSpline_2D_fast.h>

Inheritance diagram for MELANCSpline_2D_fast:

Collaboration diagram for MELANCSpline_2D_fast:

Public Member Functions
	MELANCSpline_2D_fast ()

	MELANCSpline_2D_fast (const char name, const char title)

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

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

virtual TObject *	clone (const char *newname) const

virtual	~MELANCSpline_2D_fast ()

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

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 ()

unsigned int	npointsY () const

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) const

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

virtual std::vector< std::vector< T > >	getCoefficientsPerY (const std::vector< T > &kappaX, const TMatrix_t &xAinv, const Int_t &ybin, MELANCSplineCore::BoundaryCondition const &bcBegin, MELANCSplineCore::BoundaryCondition const &bcEnd, const Int_t xbin) 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
T	rangeYmin

T	rangeYmax

const BoundaryCondition	bcBeginX

const BoundaryCondition	bcEndX

const BoundaryCondition	bcBeginY

const BoundaryCondition	bcEndY

RooRealProxy	theYVar

std::vector< T >	YList

std::vector< std::vector< T > >	FcnList

std::vector< T >	kappaX

std::vector< T >	kappaY

std::vector< std::vector< 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 11 of file MELANCSpline_2D_fast.h.

Constructor & Destructor Documentation

◆ MELANCSpline_2D_fast() [1/4]

MELANCSpline_2D_fast::MELANCSpline_2D_fast ( )

◆ MELANCSpline_2D_fast() [2/4]

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

Definition at line 23 of file MELANCSpline_2D_fast.cc.

     :
   MELANCSplineCore(name, title),
   rangeYmin(1), rangeYmax(-1),
   bcBeginX(MELANCSplineCore::bcNaturalSpline), bcEndX(MELANCSplineCore::bcNaturalSpline),
   bcBeginY(MELANCSplineCore::bcNaturalSpline), bcEndY(MELANCSplineCore::bcNaturalSpline),
   theYVar("theYVar", "theYVar", this)
 {}

◆ MELANCSpline_2D_fast() [3/4]

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

Definition at line 34 of file MELANCSpline_2D_fast.cc.

     :
   MELANCSplineCore(name, title, inXVar, inXList, inUseFloor, inFloorEval, inFloorInt),
   rangeYmin(1), rangeYmax(-1),
   bcBeginX(bcBeginX_), bcEndX(bcEndX_),
   bcBeginY(bcBeginY_), bcEndY(bcEndY_),
   theYVar("theYVar", "theYVar", this, inYVar),
   YList(inYList),
   FcnList(inFcnList)
 {
   if (npointsX()>1 && npointsY()>1){
     // Prepare A and kappa arrays for x and y coordinates
     int npoints;
     Double_t det;
  
     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); } }
     det=0;
     TMatrix_t xAinv = xAtrans.Invert(&det);
     if (det==0.){
       coutE(InputArguments) << "MELANCSpline_2D_fast::interpolateFcn: Matrix xA could not be inverted. Something is wrong with the x coordinates of points!" << endl;
       assert(0);
     }
  
     vector<vector<MELANCSplineCore::T>> yA; getKappas(kappaY, 1); getAArray(kappaY, yA, bcBeginY, bcEndY);
     npoints=kappaY.size();
     TMatrix_t yAtrans(npoints, npoints);
     for (int i=0; i<npoints; i++){ for (int j=0; j<npoints; j++){ yAtrans[i][j]=yA.at(i).at(j); } }
     det=0;
     TMatrix_t yAinv = yAtrans.Invert(&det);
     if (det==0.){
       coutE(InputArguments) << "MELANCSpline_2D_fast::interpolateFcn: Matrix yA could not be inverted. Something is wrong with the y coordinates of points!" << endl;
       assert(0);
     }
  
     // Get the grid of coefficients
     vector<vector<vector<MELANCSplineCore::T>>> coefsAlongY; // [Ax(y),Bx(y),Cx(y),Dx(y)][xbin][ybin]
     int npoldim=0;
     int nxbins=0;
     for (unsigned int j=0; j<npointsY(); j++){
       vector<vector<MELANCSplineCore::T>> xcoefsAtYj = getCoefficientsPerY(kappaX, xAinv, j, bcBeginX, bcEndX, -1); // [ix][Ax,Bx,Cx,Dx] at each y_j
       if (j==0){
         nxbins=xcoefsAtYj.size();
         npoldim=xcoefsAtYj.at(0).size();
         for (int ipow=0; ipow<npoldim; ipow++){
           vector<vector<MELANCSplineCore::T>> dum_xycoefarray;
           for (int ix=0; ix<nxbins; ix++){
             vector<MELANCSplineCore::T> dum_ycoefarray;
             dum_xycoefarray.push_back(dum_ycoefarray);
           }
           coefsAlongY.push_back(dum_xycoefarray);
         }
       }
       if (nxbins!=(int)xcoefsAtYj.size() || npoldim!=(int)xcoefsAtYj.at(0).size()){
         coutE(InputArguments) << "MELANCSpline_2D_fast::interpolateFcn: nxbins!=(int)xcoefsAtYj.size() || npoldim!=(int)xcoefsAtYj.at(0).size()!" << endl;
         assert(0);
       }
       for (int ix=0; ix<nxbins; ix++){
         for (int ipow=0; ipow<npoldim; ipow++) coefsAlongY.at(ipow).at(ix).push_back(xcoefsAtYj.at(ix).at(ipow));
       }
     }
  
     for (int ix=0; ix<nxbins; ix++){
       // Get the x coefficients interpolated across y
       vector<vector<vector<MELANCSplineCore::T>>> xCoefs;
       for (int ic=0; ic<npoldim; ic++){
         vector<vector<MELANCSplineCore::T>> yCoefs = getCoefficientsAlongDirection(kappaY, yAinv, coefsAlongY.at(ic).at(ix), bcBeginY, bcEndY, -1); // [iy][A,B,C,D]
         xCoefs.push_back(yCoefs);
       }
       coefficients.push_back(xCoefs);
     }
   }
   else assert(0);
  
   RooArgSet leafset;
   getLeafDependents(theXVar, leafset);
   getLeafDependents(theYVar, leafset);
   addLeafDependents(leafset);
  
   emptyFcnList();
 }

◆ MELANCSpline_2D_fast() [4/4]

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

Definition at line 132 of file MELANCSpline_2D_fast.cc.

     :
   MELANCSplineCore(other, name),
   rangeYmin(other.rangeYmin), rangeYmax(other.rangeYmax),
   bcBeginX(other.bcBeginX), bcEndX(other.bcEndX),
   bcBeginY(other.bcBeginY), bcEndY(other.bcEndY),
   theYVar("theYVar", this, other.theYVar),
   YList(other.YList),
   FcnList(other.FcnList),
   kappaX(other.kappaX),
   kappaY(other.kappaY),
   coefficients(other.coefficients)
 {}

◆ ~MELANCSpline_2D_fast()

virtual MELANCSpline_2D_fast::~MELANCSpline_2D_fast ( )

inlinevirtual

Definition at line 54 of file MELANCSpline_2D_fast.h.

54 {}

Member Function Documentation

◆ analyticalIntegral()

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

virtual

Implements MELANCSplineCore.

Definition at line 347 of file MELANCSpline_2D_fast.cc.

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

◆ clone()

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

inlinevirtual

Implements MELANCSplineCore.

Definition at line 53 of file MELANCSpline_2D_fast.h.

53 { return new MELANCSpline_2D_fast(*this, newname); }

◆ cropValueForRange()

void MELANCSpline_2D_fast::cropValueForRange	(	T &	val,
		const Int_t	whichDirection
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 382 of file MELANCSpline_2D_fast.cc.

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

◆ emptyFcnList()

virtual void MELANCSpline_2D_fast::emptyFcnList ( )

inlineprotectedvirtual

Implements MELANCSplineCore.

Definition at line 62 of file MELANCSpline_2D_fast.h.

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

◆ evaluate()

Double_t MELANCSpline_2D_fast::evaluate ( ) const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 313 of file MELANCSpline_2D_fast.cc.

                                              {
   Double_t value = interpolateFcn(0);
   if (useFloor && value<floorEval){
     if (verbosity>=MELANCSplineCore::kError) coutE(Eval) << "MELANCSpline_2D_fast ERROR::MELANCSpline_2D_fast(" << GetName() << ") evaluation returned " << value << " at (x, y) = (" << theXVar << ", " << theYVar << ")" << endl;
     value = floorEval;
   }
   if (verbosity==MELANCSplineCore::kVerbose){ cout << "MELANCSpline_2D_fast(" << GetName() << ")::evaluate = " << value << " at (x, y) = (" << theXVar << ", " << theYVar << ")" << endl; }
   return value;
 }

◆ getAnalyticalIntegral()

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

virtual

Implements MELANCSplineCore.

Definition at line 322 of file MELANCSpline_2D_fast.cc.

                                                                                                              {
   if (_forceNumInt) return 0;
  
   Int_t code=1;
  
   RooArgSet Xdeps, Ydeps;
   RooRealVar* rrv_x = dynamic_cast<RooRealVar*>(theXVar.absArg());
   RooRealVar* rrv_y = dynamic_cast<RooRealVar*>(theYVar.absArg());
   if (rrv_x==0) theXVar.absArg()->leafNodeServerList(&Xdeps, 0, true);
   if (rrv_y==0) theYVar.absArg()->leafNodeServerList(&Ydeps, 0, true);
  
   if (rrv_x!=0){
     if (Ydeps.find(*rrv_x)==0 || rrv_y!=0){
       if (matchArgs(allVars, analVars, theXVar)) code*=2;
     }
   }
   if (rrv_y!=0){
     if (Xdeps.find(*rrv_y)==0 || rrv_x!=0){
       if (matchArgs(allVars, analVars, theYVar)) code*=3;
     }
   }
  
   if (code==1) code=0;
   return code;
 }

◆ getCoefficientsPerY()

vector< vector< MELANCSplineCore::T > > MELANCSpline_2D_fast::getCoefficientsPerY	(	const std::vector< T > &	kappaX,
		const TMatrix_t &	xAinv,
		const Int_t &	ybin,
		MELANCSplineCore::BoundaryCondition const &	bcBegin,
		MELANCSplineCore::BoundaryCondition const &	bcEnd,
		const Int_t	xbin
	)		const

protectedvirtual

Definition at line 306 of file MELANCSpline_2D_fast.cc.

                                                                                                                                                                                                                                                                                          {
   vector<MELANCSplineCore::T> fcnList;
   for (unsigned int bin=0; bin<npointsX(); bin++){ fcnList.push_back(FcnList.at(ybin).at(bin)); }
   vector<vector<MELANCSplineCore::T>> coefs = getCoefficientsAlongDirection(kappaX, xAinv, fcnList, bcBegin, bcEnd, xbin);
   return coefs;
 }

◆ getKappas()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 245 of file MELANCSpline_2D_fast.cc.

                                                                                                  {
   kappas.clear();
   MELANCSplineCore::T kappa=1;
  
   Int_t npoints;
   vector<MELANCSplineCore::T> const* coord;
   if (whichDirection==0){
     npoints=npointsX();
     coord=&XList;
   }
   else{
     npoints=npointsY();
     coord=&YList;
   }
  
   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_2D_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 298 of file MELANCSpline_2D_fast.cc.

                                                                                                                                                                        {
   const MELANCSplineCore::T& K=kappas.at(bin);
   vector<MELANCSplineCore::T> const* coord;
   if (whichDirection==0) coord=&XList;
   else coord=&YList;
   return (val-coord->at(bin))*K;
 }

◆ getWhichBin()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 270 of file MELANCSpline_2D_fast.cc.

                                                                                                     {
   Int_t bin=-1;
   MELANCSplineCore::T valj, valjpo;
   Int_t npoints;
   vector<MELANCSplineCore::T> const* coord;
   if (whichDirection==0){
     coord=&XList;
     npoints=npointsX();
   }
   else{
     coord=&YList;
     npoints=npointsY();
   }
  
   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_2D_fast::interpolateFcn	(	Int_t	code,
		const char *	rangeName = `0`
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 149 of file MELANCSpline_2D_fast.cc.

                                                                                             {
   DefaultAccumulator<MELANCSplineCore::T> res;
  
   if (verbosity==MELANCSplineCore::kVerbose){ cout << "MELANCSpline_2D_fast(" << GetName() << ")::interpolateFcn begin with code: " << code << endl; }
  
   // Get bins
   Int_t xbin=-1, xbinmin=-1, xbinmax=-1, ybin=-1, ybinmin=-1, ybinmax=-1;
   MELANCSplineCore::T tx=0, txmin=0, txmax=0, ty=0, tymin=0, tymax=0;
   if (code==0 || code%2!=0){ // Case to just compute the value at x
     if (!testRangeValidity(theXVar, 0)) 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, 0);
     MELANCSplineCore::T coordmax = theXVar.max(rangeName); cropValueForRange(coordmax, 0);
     xbinmin = getWhichBin(coordmin, 0);
     txmin = getTVar(kappaX, coordmin, xbinmin, 0);
     xbinmax = getWhichBin(coordmax, 0);
     txmax = getTVar(kappaX, coordmax, xbinmax, 0);
   }
   if (code==0 || code%3!=0){ // Case to just compute the value at y
     if (!testRangeValidity(theYVar, 1)) return 0;
     ybin = getWhichBin(theYVar, 1);
     ty = getTVar(kappaY, theYVar, ybin, 1);
   }
   else{ // Case to integrate along y
     MELANCSplineCore::T coordmin = theYVar.min(rangeName); cropValueForRange(coordmin, 1);
     MELANCSplineCore::T coordmax = theYVar.max(rangeName); cropValueForRange(coordmax, 1);
     ybinmin = getWhichBin(coordmin, 1);
     tymin = getTVar(kappaY, coordmin, ybinmin, 1);
     ybinmax = getWhichBin(coordmax, 1);
     tymax = getTVar(kappaY, coordmax, ybinmax, 1);
   }
  
   for (int ix=0; ix<(int)coefficients.size(); 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;
  
     if (verbosity==MELANCSplineCore::kVerbose){
       if (code==0 || code%2!=0) cout << "Evaluating tx=" << txhigh << " in bin " << ix << endl;
       else cout << "Evaluating tx[" << txlow << ", " << txhigh << "] in bin " << ix << endl;
     }
  
     // Get the x coefficients interpolated across y
     vector<MELANCSplineCore::T> xCoefs;
     for (int ic=0; ic<(int)coefficients.at(ix).size(); ic++){
       const vector<vector<MELANCSplineCore::T>>& yCoefs = coefficients.at(ix).at(ic);
  
       if (verbosity==MELANCSplineCore::kVerbose) cout << "\tCoefficient " << ic << ":\n";
  
       DefaultAccumulator<MELANCSplineCore::T> theCoef;
       for (int iy=0; iy<(int)yCoefs.size(); iy++){
         if (
           (ybin>=0 && iy!=ybin)
           ||
           (ybinmin>=0 && ybinmax>=ybinmin && !(ybinmin<=iy && iy<=ybinmax))
           ) continue;
  
         MELANCSplineCore::T tylow=0, tyhigh=1;
         if (code>0 && code%3==0){
           if (iy==ybinmin) tylow=tymin;
           if (iy==ybinmax) tyhigh=tymax;
         }
         else tyhigh=ty;
  
         if (verbosity==MELANCSplineCore::kVerbose){
           if (code==0 || code%3!=0) cout << "\tEvaluating ty=" << tyhigh << " in bin " << iy << endl;
           else cout << "\tEvaluating ty[" << tylow << ", " << tyhigh << "] in bin " << iy << endl;
         }
  
         theCoef += evalSplineSegment(yCoefs.at(iy), kappaY.at(iy), tyhigh, tylow, (code>0 && code%3==0));
       }
  
       //if (code==0) cout << "\tCoefficient is " << theCoef << endl;
  
       xCoefs.push_back(theCoef);
     }
  
     // Evaluate value of spline at x with coefficients evaluated at y
     res += evalSplineSegment(xCoefs, kappaX.at(ix), txhigh, txlow, (code>0 && code%2==0));
   }
  
   return res;
 }

◆ npointsY()

unsigned int MELANCSpline_2D_fast::npointsY ( ) const

inlineprotected

Definition at line 64 of file MELANCSpline_2D_fast.h.

64 { return YList.size(); }

◆ setRangeValidity()

void MELANCSpline_2D_fast::setRangeValidity	(	const T	valmin,
		const T	valmax,
		const Int_t	whichDirection
	)

virtual

Implements MELANCSplineCore.

Definition at line 369 of file MELANCSpline_2D_fast.cc.

                                                                                                      {
   T* range[2];
   if (whichDirection==0){
     range[0] = &rangeXmin;
     range[1] = &rangeXmax;
   }
   else{
     range[0] = &rangeYmin;
     range[1] = &rangeYmax;
   }
   *(range[0])=valmin;
   *(range[1])=valmax;
 }

◆ testRangeValidity()

Bool_t MELANCSpline_2D_fast::testRangeValidity	(	const T &	val,
		const Int_t	whichDirection
	)		const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 357 of file MELANCSpline_2D_fast.cc.

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

Member Data Documentation

◆ bcBeginX

const BoundaryCondition MELANCSpline_2D_fast::bcBeginX

protected

Definition at line 16 of file MELANCSpline_2D_fast.h.

◆ bcBeginY

const BoundaryCondition MELANCSpline_2D_fast::bcBeginY

protected

Definition at line 18 of file MELANCSpline_2D_fast.h.

◆ bcEndX

const BoundaryCondition MELANCSpline_2D_fast::bcEndX

protected

Definition at line 17 of file MELANCSpline_2D_fast.h.

◆ bcEndY

const BoundaryCondition MELANCSpline_2D_fast::bcEndY

protected

Definition at line 19 of file MELANCSpline_2D_fast.h.

◆ coefficients

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

protected

Definition at line 28 of file MELANCSpline_2D_fast.h.

◆ FcnList

std::vector<std::vector<T> > MELANCSpline_2D_fast::FcnList

protected

Definition at line 24 of file MELANCSpline_2D_fast.h.

◆ kappaX

std::vector<T> MELANCSpline_2D_fast::kappaX

protected

Definition at line 26 of file MELANCSpline_2D_fast.h.

◆ kappaY

std::vector<T> MELANCSpline_2D_fast::kappaY

protected

Definition at line 27 of file MELANCSpline_2D_fast.h.

◆ rangeYmax

T MELANCSpline_2D_fast::rangeYmax

protected

Definition at line 14 of file MELANCSpline_2D_fast.h.

◆ rangeYmin

T MELANCSpline_2D_fast::rangeYmin

protected

Definition at line 13 of file MELANCSpline_2D_fast.h.

◆ theYVar

RooRealProxy MELANCSpline_2D_fast::theYVar

protected

Definition at line 21 of file MELANCSpline_2D_fast.h.

◆ YList

std::vector<T> MELANCSpline_2D_fast::YList

protected

Definition at line 22 of file MELANCSpline_2D_fast.h.

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

MELA/interface/MELANCSpline_2D_fast.h
MELA/src/MELANCSpline_2D_fast.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MELANCSpline_2D_fast() [1/4]

◆ MELANCSpline_2D_fast() [2/4]

◆ MELANCSpline_2D_fast() [3/4]

◆ MELANCSpline_2D_fast() [4/4]

◆ ~MELANCSpline_2D_fast()

Member Function Documentation

◆ analyticalIntegral()

◆ clone()

◆ cropValueForRange()

◆ emptyFcnList()

◆ evaluate()

◆ getAnalyticalIntegral()

◆ getCoefficientsPerY()

◆ getKappas()

◆ getTVar()

◆ getWhichBin()

◆ interpolateFcn()

◆ npointsY()

◆ setRangeValidity()

◆ testRangeValidity()

Member Data Documentation

◆ bcBeginX

◆ bcBeginY

◆ bcEndX

◆ bcEndY

◆ coefficients

◆ FcnList

◆ kappaX

◆ kappaY

◆ rangeYmax

◆ rangeYmin

◆ theYVar

◆ YList