#include <MELANCSpline_3D_fast.h>

Inheritance diagram for MELANCSpline_3D_fast:

Collaboration diagram for MELANCSpline_3D_fast:

Public Member Functions
	MELANCSpline_3D_fast ()

	MELANCSpline_3D_fast (const char name, const char title)

	MELANCSpline_3D_fast (const char name, const char title, RooAbsReal &inXVar, RooAbsReal &inYVar, RooAbsReal &inZVar, const std::vector< T > &inXList, const std::vector< T > &inYList, const std::vector< T > &inZList, const std::vector< 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, MELANCSplineCore::BoundaryCondition const bcBeginZ_=MELANCSplineCore::bcNaturalSpline, MELANCSplineCore::BoundaryCondition const bcEndZ_=MELANCSplineCore::bcNaturalSpline, Bool_t inUseFloor=true, T inFloorEval=0, T inFloorInt=0)

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

virtual TObject *	clone (const char *newname) const

virtual	~MELANCSpline_3D_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

unsigned int	npointsZ () 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 > >	getCoefficientsPerYPerZ (const std::vector< T > &kappaX, const TMatrix_t &xAinv, const Int_t &ybin, const Int_t &zbin, 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

T	rangeZmin

T	rangeZmax

const BoundaryCondition	bcBeginX

const BoundaryCondition	bcEndX

const BoundaryCondition	bcBeginY

const BoundaryCondition	bcEndY

const BoundaryCondition	bcBeginZ

const BoundaryCondition	bcEndZ

RooRealProxy	theYVar

RooRealProxy	theZVar

std::vector< T >	YList

std::vector< T >	ZList

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

std::vector< T >	kappaX

std::vector< T >	kappaY

std::vector< T >	kappaZ

std::vector< std::vector< 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_3D_fast.h.

Constructor & Destructor Documentation

◆ MELANCSpline_3D_fast() [1/4]

MELANCSpline_3D_fast::MELANCSpline_3D_fast ( )

◆ MELANCSpline_3D_fast() [2/4]

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

Definition at line 26 of file MELANCSpline_3D_fast.cc.

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

◆ MELANCSpline_3D_fast() [3/4]

MELANCSpline_3D_fast::MELANCSpline_3D_fast	(	const char *	name,
		const char *	title,
		RooAbsReal &	inXVar,
		RooAbsReal &	inYVar,
		RooAbsReal &	inZVar,
		const std::vector< T > &	inXList,
		const std::vector< T > &	inYList,
		const std::vector< T > &	inZList,
		const std::vector< 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`,
		MELANCSplineCore::BoundaryCondition const	bcBeginZ_ = `MELANCSplineCore::bcNaturalSpline`,
		MELANCSplineCore::BoundaryCondition const	bcEndZ_ = `MELANCSplineCore::bcNaturalSpline`,
		Bool_t	inUseFloor = `true`,
		T	inFloorEval = `0`,
		T	inFloorInt = `0`
	)

Definition at line 40 of file MELANCSpline_3D_fast.cc.

     :
   MELANCSplineCore(name, title, inXVar, inXList, inUseFloor, inFloorEval, inFloorInt),
   rangeYmin(1), rangeYmax(-1),
   rangeZmin(1), rangeZmax(-1),
   bcBeginX(bcBeginX_), bcEndX(bcEndX_),
   bcBeginY(bcBeginY_), bcEndY(bcEndY_),
   bcBeginZ(bcBeginZ_), bcEndZ(bcEndZ_),
   theYVar("theYVar", "theYVar", this, inYVar),
   theZVar("theZVar", "theZVar", this, inZVar),
   YList(inYList),
   ZList(inZList),
   FcnList(inFcnList)
 {
   if (npointsX()>1 && npointsY()>1 && npointsZ()>1){
     // Prepare A and kappa arrays for x, y and z 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_3D_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_3D_fast::interpolateFcn: Matrix yA could not be inverted. Something is wrong with the y coordinates of points!" << endl;
       assert(0);
     }
  
     vector<vector<MELANCSplineCore::T>> zA; getKappas(kappaZ, 2); getAArray(kappaZ, zA, bcBeginZ, bcEndZ);
     npoints=kappaZ.size();
     TMatrix_t zAtrans(npoints, npoints);
     for (int i=0; i<npoints; i++){ for (int j=0; j<npoints; j++){ zAtrans[i][j]=zA.at(i).at(j); } }
     det=0;
     TMatrix_t zAinv = zAtrans.Invert(&det);
     if (det==0.){
       coutE(InputArguments) << "MELANCSpline_3D_fast::interpolateFcn: Matrix zA could not be inverted. Something is wrong with the z coordinates of points!" << endl;
       assert(0);
     }
  
     //cout << "MELANCSpline_3D_fast::MELANCSpline_3D_fast: Initial kappa arrays are set up" << endl;
  
     // Get the grid of coefficients
     int nxpoldim=0;
     int nxbins=0;
     int nybins=0;
     int nypoldim=0;
     std::vector<
       std::vector<std::vector<
       std::vector<std::vector<MELANCSplineCore::T>>
       >>
       > coefsAlongZ; // [ix][A_x,B_x,C_x,D_x][iy][A_x_y,B_x_y,C_x_y,D_x_y][z_k] at each z_k
     for (unsigned int k=0; k<npointsZ(); k++){
       //cout << "Finding coefficients along z line " << k << endl;
  
       std::vector<std::vector<
         std::vector<std::vector<MELANCSplineCore::T>>
         >> coefficients_perZ; // [ix][A_x,B_x,C_x,D_x][iy][A_x_y,B_x_y,C_x_y,D_x_y] at each z_k
  
       vector<vector<vector<MELANCSplineCore::T>>> coefsAlongY; // [xbin][Ax(y),Bx(y),Cx(y),Dx(y)][ybin] in each z
       for (unsigned int j=0; j<npointsY(); j++){
         vector<vector<MELANCSplineCore::T>> xcoefsAtYjZk = getCoefficientsPerYPerZ(kappaX, xAinv, j, k, bcBeginX, bcEndX, -1); // [ix][Ax,Bx,Cx,Dx] at each y_j z_k
         //cout << "\tCoefficients in y line " << j << " are found" << endl;
         if (j==0){
           if (k==0){
             nxbins=xcoefsAtYjZk.size();
             nxpoldim=xcoefsAtYjZk.at(0).size();
           }
           for (int ix=0; ix<nxbins; ix++){
             vector<vector<MELANCSplineCore::T>> dum_xycoefarray;
             for (int icx=0; icx<nxpoldim; icx++){
               vector<MELANCSplineCore::T> dum_ycoefarray;
               dum_xycoefarray.push_back(dum_ycoefarray);
             }
             coefsAlongY.push_back(dum_xycoefarray);
           }
         }
         if (nxbins!=(int)xcoefsAtYjZk.size() || nxpoldim!=(int)xcoefsAtYjZk.at(0).size()){
           coutE(InputArguments) << "MELANCSpline_3D_fast::interpolateFcn: nxbins!=(int)xcoefsAtYjZk.size() || nxpoldim!=(int)xcoefsAtYjZk.at(0).size()!" << endl;
           assert(0);
         }
         for (int ix=0; ix<nxbins; ix++){
           for (int icx=0; icx<nxpoldim; icx++) coefsAlongY.at(ix).at(icx).push_back(xcoefsAtYjZk.at(ix).at(icx));
         }
       }
  
       //cout << "\tCoefficients in each y line are found" << endl;
  
       for (int ix=0; ix<nxbins; ix++){
         // Get the x coefficients interpolated across y
         vector<vector<vector<MELANCSplineCore::T>>> xCoefs;
         for (int icx=0; icx<nxpoldim; icx++){
           vector<vector<MELANCSplineCore::T>> yCoefs = getCoefficientsAlongDirection(kappaY, yAinv, coefsAlongY.at(ix).at(icx), bcBeginY, bcEndY, -1); // [iy][A,B,C,D]
           xCoefs.push_back(yCoefs);
         }
         coefficients_perZ.push_back(xCoefs);
       }
  
       if (k==0){
         nybins = coefficients_perZ.at(0).at(0).size();
         nypoldim = coefficients_perZ.at(0).at(0).at(0).size();
         for (int ix=0; ix<nxbins; ix++){
           vector<vector<vector<vector<MELANCSplineCore::T>>>> xCoefs;
           for (int icx=0; icx<nxpoldim; icx++){
             vector<vector<vector<MELANCSplineCore::T>>> xCoefsAlongY;
             for (int iy=0; iy<nybins; iy++){
               vector<vector<MELANCSplineCore::T>> yCoefs;
               for (int icy=0; icy<nypoldim; icy++){
                 vector<MELANCSplineCore::T> yCoefAtZj;
                 yCoefs.push_back(yCoefAtZj);
               }
               xCoefsAlongY.push_back(yCoefs);
             }
             xCoefs.push_back(xCoefsAlongY);
           }
           coefsAlongZ.push_back(xCoefs);
         }
       }
       for (int ix=0; ix<nxbins; ix++){
         for (int icx=0; icx<nxpoldim; icx++){
           for (int iy=0; iy<nybins; iy++){
             for (int icy=0; icy<nypoldim; icy++){
               coefsAlongZ.at(ix).at(icx).at(iy).at(icy).push_back(coefficients_perZ.at(ix).at(icx).at(iy).at(icy));
             }
           }
         }
       }
     }
  
     //cout << "Finding the final coefficients" << endl;
     for (int ix=0; ix<nxbins; ix++){
       vector<vector<vector<vector<vector<MELANCSplineCore::T>>>>> xCoefs;
       for (int icx=0; icx<nxpoldim; icx++){
         vector<vector<vector<vector<MELANCSplineCore::T>>>> xCoefsAlongY;
         for (int iy=0; iy<nybins; iy++){
           vector<vector<vector<MELANCSplineCore::T>>> yCoefs;
           for (int icy=0; icy<nypoldim; icy++){
             vector<vector<MELANCSplineCore::T>> yCoefsAlongZ = getCoefficientsAlongDirection(kappaZ, zAinv, coefsAlongZ.at(ix).at(icx).at(iy).at(icy), bcBeginZ, bcEndZ, -1); // [iz][A,B,C,D]
             yCoefs.push_back(yCoefsAlongZ);
           }
           xCoefsAlongY.push_back(yCoefs);
         }
         xCoefs.push_back(xCoefsAlongY);
       }
       coefficients.push_back(xCoefs);
     }
  
   }
   else assert(0);
  
   RooArgSet leafset;
   getLeafDependents(theXVar, leafset);
   getLeafDependents(theYVar, leafset);
   getLeafDependents(theZVar, leafset);
   addLeafDependents(leafset);
  
   emptyFcnList();
 }

◆ MELANCSpline_3D_fast() [4/4]

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

Definition at line 229 of file MELANCSpline_3D_fast.cc.

     :
   MELANCSplineCore(other, name),
   rangeYmin(other.rangeYmin), rangeYmax(other.rangeYmax),
   rangeZmin(other.rangeZmin), rangeZmax(other.rangeZmax),
   bcBeginX(other.bcBeginX), bcEndX(other.bcEndX),
   bcBeginY(other.bcBeginY), bcEndY(other.bcEndY),
   bcBeginZ(other.bcBeginZ), bcEndZ(other.bcEndZ),
   theYVar("theYVar", this, other.theYVar),
   theZVar("theZVar", this, other.theZVar),
   YList(other.YList),
   ZList(other.ZList),
   FcnList(other.FcnList),
   kappaX(other.kappaX),
   kappaY(other.kappaY),
   kappaZ(other.kappaZ),
   coefficients(other.coefficients)
 {}

◆ ~MELANCSpline_3D_fast()

virtual MELANCSpline_3D_fast::~MELANCSpline_3D_fast ( )

inlinevirtual

Definition at line 69 of file MELANCSpline_3D_fast.h.

69 {}

Member Function Documentation

◆ analyticalIntegral()

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

virtual

Implements MELANCSplineCore.

Definition at line 491 of file MELANCSpline_3D_fast.cc.

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

◆ clone()

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

inlinevirtual

Implements MELANCSplineCore.

Definition at line 68 of file MELANCSpline_3D_fast.h.

68 { return new MELANCSpline_3D_fast(*this, newname); }

◆ cropValueForRange()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 534 of file MELANCSpline_3D_fast.cc.

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

◆ emptyFcnList()

virtual void MELANCSpline_3D_fast::emptyFcnList ( )

inlineprotectedvirtual

Implements MELANCSplineCore.

Definition at line 77 of file MELANCSpline_3D_fast.h.

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

◆ evaluate()

Double_t MELANCSpline_3D_fast::evaluate ( ) const

protectedvirtual

Implements MELANCSplineCore.

Definition at line 438 of file MELANCSpline_3D_fast.cc.

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

◆ getAnalyticalIntegral()

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

virtual

Implements MELANCSplineCore.

Definition at line 447 of file MELANCSpline_3D_fast.cc.

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

◆ getCoefficientsPerYPerZ()

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

protectedvirtual

Definition at line 426 of file MELANCSpline_3D_fast.cc.

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

◆ getKappas()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 356 of file MELANCSpline_3D_fast.cc.

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

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

◆ getWhichBin()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 385 of file MELANCSpline_3D_fast.cc.

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 250 of file MELANCSpline_3D_fast.cc.

                                                                                             {
   DefaultAccumulator<MELANCSplineCore::T> res;
  
   if (verbosity==MELANCSplineCore::kVerbose) cout << "MELANCSpline_3D_fast(" << GetName() << ")::interpolateFcn begin with code: " << code << endl;
  
   // Get bins
   vector<Int_t> varprime; varprime.push_back(2); varprime.push_back(3); varprime.push_back(5);
   const Int_t ndims=(const Int_t)varprime.size();
   vector<Int_t> varbin;
   vector<Int_t> varbinmin;
   vector<Int_t> varbinmax;
   vector<MELANCSplineCore::T> tvar;
   vector<MELANCSplineCore::T> tvarmin;
   vector<MELANCSplineCore::T> tvarmax;
   vector<const vector<MELANCSplineCore::T>*> varkappa; varkappa.push_back(&kappaX); varkappa.push_back(&kappaY); varkappa.push_back(&kappaZ);
   vector<const RooRealProxy*> varcoord; varcoord.push_back(&theXVar); varcoord.push_back(&theYVar); varcoord.push_back(&theZVar);
   for (Int_t idim=0; idim<ndims; idim++){
     Int_t binval=-1;
     Int_t binmin=-1;
     Int_t binmax=-1;
     MELANCSplineCore::T tval=0;
     MELANCSplineCore::T tmin=0;
     MELANCSplineCore::T tmax=0;
     if (code==0 || code%varprime.at(idim)!=0){
       if (!testRangeValidity(*(varcoord.at(idim)), idim)) return 0;
       binval = getWhichBin(*(varcoord.at(idim)), idim);
       tval = getTVar(*(varkappa.at(idim)), *(varcoord.at(idim)), binval, idim);
     }
     else{
       MELANCSplineCore::T coordmin = varcoord.at(idim)->min(rangeName); cropValueForRange(coordmin, idim);
       MELANCSplineCore::T coordmax = varcoord.at(idim)->max(rangeName); cropValueForRange(coordmax, idim);
       const std::vector<MELANCSplineCore::T>& kappadim = *(varkappa.at(idim));
       binmin = getWhichBin(coordmin, idim);
       tmin = getTVar(kappadim, coordmin, binmin, idim);
       binmax = getWhichBin(coordmax, idim);
       tmax = getTVar(kappadim, coordmax, binmax, idim);
     }
     varbin.push_back(binval);
     varbinmin.push_back(binmin);
     varbinmax.push_back(binmax);
     tvar.push_back(tval);
     tvarmin.push_back(tmin);
     tvarmax.push_back(tmax);
   }
  
   for (int ix=0; ix<(int)coefficients.size(); ix++){
     if (
       (varbin[0]>=0 && ix!=varbin[0])
       ||
       (varbinmin[0]>=0 && varbinmax[0]>=varbinmin[0] && !(varbinmin[0]<=ix && ix<=varbinmax[0]))
       ) continue;
     MELANCSplineCore::T txlow=0, txhigh=1;
     if (code>0 && code%varprime[0]==0){
       if (ix==varbinmin[0]) txlow=tvarmin[0];
       if (ix==varbinmax[0]) txhigh=tvarmax[0];
     }
     else txhigh=tvar[0];
     // Get the x coefficients interpolated across y
     vector<MELANCSplineCore::T> xCoefs;
     for (int icx=0; icx<(int)coefficients.at(ix).size(); icx++){
       DefaultAccumulator<MELANCSplineCore::T> theXCoef;
       for (int iy=0; iy<(int)coefficients.at(ix).at(icx).size(); iy++){
         if (
           (varbin[1]>=0 && iy!=varbin[1])
           ||
           (varbinmin[1]>=0 && varbinmax[1]>=varbinmin[1] && !(varbinmin[1]<=iy && iy<=varbinmax[1]))
           ) continue;
         MELANCSplineCore::T tylow=0, tyhigh=1;
         if (code>0 && code%varprime[1]==0){
           if (iy==varbinmin[1]) tylow=tvarmin[1];
           if (iy==varbinmax[1]) tyhigh=tvarmax[1];
         }
         else tyhigh=tvar[1];
         // Get the y coefficients interpolated across z
         vector<MELANCSplineCore::T> yCoefs;
         for (int icy=0; icy<(int)coefficients.at(ix).at(icx).at(iy).size(); icy++){
           DefaultAccumulator<MELANCSplineCore::T> theYCoef;
           for (int iz=0; iz<(int)coefficients.at(ix).at(icx).at(iy).at(icy).size(); iz++){
             if (
               (varbin[2]>=0 && iz!=varbin[2])
               ||
               (varbinmin[2]>=0 && varbinmax[2]>=varbinmin[2] && !(varbinmin[2]<=iz && iz<=varbinmax[2]))
               ) continue;
  
             MELANCSplineCore::T tzlow=0, tzhigh=1;
             if (code>0 && code%varprime[2]==0){
               if (iz==varbinmin[2]) tzlow=tvarmin[2];
               if (iz==varbinmax[2]) tzhigh=tvarmax[2];
             }
             else tzhigh=tvar[2];
  
             theYCoef += evalSplineSegment(coefficients.at(ix).at(icx).at(iy).at(icy).at(iz), varkappa[2]->at(iz), tzhigh, tzlow, (code>0 && code%varprime[2]==0));
           }
           yCoefs.push_back(theYCoef);
         }
         theXCoef += evalSplineSegment(yCoefs, varkappa[1]->at(iy), tyhigh, tylow, (code>0 && code%varprime[1]==0));
       }
       xCoefs.push_back(theXCoef);
     }
     // Evaluate value of spline at x with coefficients evaluated at y
     res += evalSplineSegment(xCoefs, varkappa[0]->at(ix), txhigh, txlow, (code>0 && code%varprime[0]==0));
   }
  
   return res;
 }

◆ npointsY()

unsigned int MELANCSpline_3D_fast::npointsY ( ) const

inlineprotected

Definition at line 79 of file MELANCSpline_3D_fast.h.

79 { return YList.size(); }

◆ npointsZ()

unsigned int MELANCSpline_3D_fast::npointsZ ( ) const

inlineprotected

Definition at line 80 of file MELANCSpline_3D_fast.h.

80 { return ZList.size(); }

◆ setRangeValidity()

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

virtual

Implements MELANCSplineCore.

Definition at line 517 of file MELANCSpline_3D_fast.cc.

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

◆ testRangeValidity()

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

protectedvirtual

Implements MELANCSplineCore.

Definition at line 501 of file MELANCSpline_3D_fast.cc.

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

Member Data Documentation

◆ bcBeginX

const BoundaryCondition MELANCSpline_3D_fast::bcBeginX

protected

Definition at line 18 of file MELANCSpline_3D_fast.h.

◆ bcBeginY

const BoundaryCondition MELANCSpline_3D_fast::bcBeginY

protected

Definition at line 20 of file MELANCSpline_3D_fast.h.

◆ bcBeginZ

const BoundaryCondition MELANCSpline_3D_fast::bcBeginZ

protected

Definition at line 22 of file MELANCSpline_3D_fast.h.

◆ bcEndX

const BoundaryCondition MELANCSpline_3D_fast::bcEndX

protected

Definition at line 19 of file MELANCSpline_3D_fast.h.

◆ bcEndY

const BoundaryCondition MELANCSpline_3D_fast::bcEndY

protected

Definition at line 21 of file MELANCSpline_3D_fast.h.

◆ bcEndZ

const BoundaryCondition MELANCSpline_3D_fast::bcEndZ

protected

Definition at line 23 of file MELANCSpline_3D_fast.h.

◆ coefficients

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

protected

Definition at line 39 of file MELANCSpline_3D_fast.h.

◆ FcnList

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

protected

Definition at line 30 of file MELANCSpline_3D_fast.h.

◆ kappaX

std::vector<T> MELANCSpline_3D_fast::kappaX

protected

Definition at line 32 of file MELANCSpline_3D_fast.h.

◆ kappaY

std::vector<T> MELANCSpline_3D_fast::kappaY

protected

Definition at line 33 of file MELANCSpline_3D_fast.h.

◆ kappaZ

std::vector<T> MELANCSpline_3D_fast::kappaZ

protected

Definition at line 34 of file MELANCSpline_3D_fast.h.

◆ rangeYmax

T MELANCSpline_3D_fast::rangeYmax

protected

Definition at line 14 of file MELANCSpline_3D_fast.h.

◆ rangeYmin

T MELANCSpline_3D_fast::rangeYmin

protected

Definition at line 13 of file MELANCSpline_3D_fast.h.

◆ rangeZmax

T MELANCSpline_3D_fast::rangeZmax

protected

Definition at line 16 of file MELANCSpline_3D_fast.h.

◆ rangeZmin

T MELANCSpline_3D_fast::rangeZmin

protected

Definition at line 15 of file MELANCSpline_3D_fast.h.

◆ theYVar

RooRealProxy MELANCSpline_3D_fast::theYVar

protected

Definition at line 25 of file MELANCSpline_3D_fast.h.

◆ theZVar

RooRealProxy MELANCSpline_3D_fast::theZVar

protected

Definition at line 26 of file MELANCSpline_3D_fast.h.

◆ YList

std::vector<T> MELANCSpline_3D_fast::YList

protected

Definition at line 27 of file MELANCSpline_3D_fast.h.

◆ ZList

std::vector<T> MELANCSpline_3D_fast::ZList

protected

Definition at line 28 of file MELANCSpline_3D_fast.h.

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

MELA/interface/MELANCSpline_3D_fast.h
MELA/src/MELANCSpline_3D_fast.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MELANCSpline_3D_fast() [1/4]

◆ MELANCSpline_3D_fast() [2/4]

◆ MELANCSpline_3D_fast() [3/4]

◆ MELANCSpline_3D_fast() [4/4]

◆ ~MELANCSpline_3D_fast()

Member Function Documentation

◆ analyticalIntegral()

◆ clone()

◆ cropValueForRange()

◆ emptyFcnList()

◆ evaluate()

◆ getAnalyticalIntegral()

◆ getCoefficientsPerYPerZ()

◆ getKappas()

◆ getTVar()

◆ getWhichBin()

◆ interpolateFcn()

◆ npointsY()

◆ npointsZ()

◆ setRangeValidity()

◆ testRangeValidity()

Member Data Documentation

◆ bcBeginX

◆ bcBeginY

◆ bcBeginZ

◆ bcEndX

◆ bcEndY

◆ bcEndZ

◆ coefficients

◆ FcnList

◆ kappaX

◆ kappaY

◆ kappaZ

◆ rangeYmax

◆ rangeYmin

◆ rangeZmax

◆ rangeZmin

◆ theYVar

◆ theZVar

◆ YList

◆ ZList