MELANCSplineCore Class Reference

#include <MELANCSplineCore.h>

Inheritance diagram for MELANCSplineCore:

Collaboration diagram for MELANCSplineCore:

Public Types
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

Public Member Functions
	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 TObject *	clone (const char *newname) const =0
virtual	~MELANCSplineCore ()
virtual void	setVerbosity (VerbosityLevel flag)
void	setEvalFloor (T val)
void	setIntFloor (T val)
void	doFloor (Bool_t flag)
virtual 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 =0
virtual Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const =0

Protected Member Functions
virtual void	emptyFcnList ()=0
void	getLeafDependents (RooRealProxy &proxy, RooArgSet &set)
void	addLeafDependents (RooArgSet &set)
unsigned int	npointsX () const
virtual Int_t	getWhichBin (const T &val, const Int_t whichDirection) const =0
virtual void	getKappas (std::vector< T > &kappas, const Int_t whichDirection)=0
virtual T	getTVar (const std::vector< T > &kappas, const T &val, const Int_t &bin, const Int_t whichDirection) const =0
virtual Bool_t	testRangeValidity (const T &val, const Int_t whichDirection) const =0
virtual void	cropValueForRange (T &val, const Int_t whichDirection) const =0
virtual T	interpolateFcn (Int_t code, const char *rangeName=0) const =0
virtual Double_t	evaluate () const =0
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
VerbosityLevel	verbosity
Bool_t	useFloor
T	floorEval
T	floorInt
T	rangeXmin
T	rangeXmax
RooRealProxy	theXVar
RooListProxy	leafDepsList
std::vector< T >	XList

Detailed Description

Definition at line 16 of file MELANCSplineCore.h.

Member Typedef Documentation

T

typedef Float_t MELANCSplineCore::T

Definition at line 18 of file MELANCSplineCore.h.

TMatrix_t

typedef TMatrixT<T> MELANCSplineCore::TMatrix_t

Definition at line 19 of file MELANCSplineCore.h.

TVector_t

typedef TVectorT<T> MELANCSplineCore::TVector_t

Definition at line 20 of file MELANCSplineCore.h.

Member Enumeration Documentation

BoundaryCondition

enum MELANCSplineCore::BoundaryCondition

Enumerator
bcApproximatedSlope
bcClamped
bcApproximatedSecondDerivative
bcNaturalSpline
bcQuadratic
bcQuadraticWithNullSlope
NBoundaryConditions

Definition at line 28 of file MELANCSplineCore.h.

                        {
    bcApproximatedSlope, // D1 is the same as deltaY/deltaX in the first/last segment
    bcClamped, // D1=0 at endpoint
    bcApproximatedSecondDerivative, // D2 is approximated
    bcNaturalSpline, // D2=0 at endpoint
    bcQuadratic, // Coefficient D=0
    bcQuadraticWithNullSlope, // Coefficient D=0 and D1=0 at endpoint
    NBoundaryConditions
  };

VerbosityLevel

enum MELANCSplineCore::VerbosityLevel

Enumerator
kSilent
kError
kVerbose

Definition at line 22 of file MELANCSplineCore.h.

                     {
    kSilent,
    kError,
    kVerbose
  };

Constructor & Destructor Documentation

MELANCSplineCore() [1/4]

MELANCSplineCore::MELANCSplineCore

(

)

MELANCSplineCore() [2/4]

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

Definition at line 24 of file MELANCSplineCore.cc.

    :
  RooAbsReal(name, title),
  verbosity(MELANCSplineCore::kSilent),
  useFloor(true), floorEval(0), floorInt(0),
  rangeXmin(1), rangeXmax(-1),
  theXVar("theXVar", "theXVar", this),
  leafDepsList("leafDepsList", "leafDepsList", this)
{}

MELANCSplineCore() [3/4]

MELANCSplineCore::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
	)

Definition at line 36 of file MELANCSplineCore.cc.

    :
  RooAbsReal(name, title),
  verbosity(MELANCSplineCore::kSilent),
  useFloor(inUseFloor), floorEval(inFloorEval), floorInt(inFloorInt),
  rangeXmin(1), rangeXmax(-1),
  theXVar("theXVar", "theXVar", this, inXVar),
  leafDepsList("leafDepsList", "leafDepsList", this),
  XList(inXList)
{}

MELANCSplineCore() [4/4]

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

Definition at line 54 of file MELANCSplineCore.cc.

    :
  RooAbsReal(other, name),
  verbosity(other.verbosity),
  useFloor(other.useFloor), floorEval(other.floorEval), floorInt(other.floorInt),
  rangeXmin(other.rangeXmin), rangeXmax(other.rangeXmax),
  theXVar("theXVar", this, other.theXVar),
  leafDepsList("leafDepsList", this, other.leafDepsList),
  XList(other.XList)
{}

~MELANCSplineCore()

virtual MELANCSplineCore::~MELANCSplineCore ( )

inlinevirtual

Definition at line 54 of file MELANCSplineCore.h.

{}

Member Function Documentation

addLeafDependents()

void MELANCSplineCore::addLeafDependents ( RooArgSet &  set )

protected

Definition at line 280 of file MELANCSplineCore.cc.

                                                      {
  TIterator* iter = set.createIterator();
  RooAbsArg* absarg;
  while ((absarg = (RooAbsArg*)iter->Next())){ if (dynamic_cast<RooRealVar*>(absarg)) leafDepsList.add(*absarg); }
  delete iter;
}

analyticalIntegral()

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

pure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

clone()

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

pure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

cropValueForRange()

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

protectedpure virtual

Implemented in MELANCSpline_1D_fast, MELANCSpline_3D_fast, and MELANCSpline_2D_fast.

doFloor()

void MELANCSplineCore::doFloor

(

Bool_t

flag

)

Definition at line 70 of file MELANCSplineCore.cc.

{ useFloor = flag; }

emptyFcnList()

virtual void MELANCSplineCore::emptyFcnList ( )

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

evalSplineSegment()

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

protectedvirtual

Definition at line 265 of file MELANCSplineCore.cc.

                                                                                                                                                                                                                   {
  DefaultAccumulator<MELANCSplineCore::T> res;
  for (unsigned int ic=0; ic<coefs.size(); ic++){
    if (doIntegrate) res += coefs.at(ic)*(pow(tup, (int)(ic+1))-pow(tdn, (int)(ic+1)))/((MELANCSplineCore::T)(ic+1));
    else res += coefs.at(ic)*pow(tup, (int)ic);
  }
  if (doIntegrate) res /= kappa;
  return res;
}

evaluate()

virtual Double_t MELANCSplineCore::evaluate ( ) const

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

getAArray()

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

protectedvirtual

Definition at line 159 of file MELANCSplineCore.cc.

                                                                                                                                                                                         {
  AArray.clear();
  Int_t npoints = kappas.size();
  for (int i=0; i<npoints; i++){
    vector<MELANCSplineCore::T> Ai(npoints, 0.);
    if (npoints==1) Ai[0]=1;
    else if (i==0){
      switch(bcBegin){
      case bcNaturalSpline:
      case bcApproximatedSecondDerivative:
        Ai[0]=2; Ai[1]=kappas.at(1)/kappas.at(0);
        break;
      case bcClamped:
      case bcApproximatedSlope:
      case bcQuadraticWithNullSlope:
        Ai[0]=1;
        break;
      case bcQuadratic:
        Ai[0]=1; Ai[1]=kappas.at(1)/kappas.at(0);
        break;
      default:
        cerr << "MELANCSplineCore::getAArray: bcBegin " << bcBegin << " is not implemented!" << endl;
        assert(0);
      }
    }
    else if (i==npoints-1){
      switch(bcEnd){
      case bcNaturalSpline:
      case bcApproximatedSecondDerivative:
        Ai[npoints-2]=1; Ai[npoints-1]=2.*kappas.at(npoints-1)/kappas.at(npoints-2);
        break;
      case bcClamped:
      case bcApproximatedSlope:
      case bcQuadraticWithNullSlope:
        Ai[npoints-1]=1;
        break;
      case bcQuadratic:
        Ai[npoints-2]=1; Ai[npoints-1]=kappas.at(npoints-1)/kappas.at(npoints-2);
        break;
      default:
        cerr << "MELANCSplineCore::getAArray: bcEnd " << bcEnd << " is not implemented!" << endl;
        assert(0);
      }
    }
    else{
      if ((i==1 && bcBegin==bcQuadraticWithNullSlope) || (i==npoints-2 && bcEnd==bcQuadraticWithNullSlope)) Ai[i]=1;
      else{
        MELANCSplineCore::T kappa_j = kappas.at(i);
        MELANCSplineCore::T kappa_jmo = kappas.at(i-1);
        MELANCSplineCore::T kappa_jpo = kappas.at(i+1);
 
        Ai[i-1]=1;
        Ai[i]=2.*kappa_j/kappa_jmo*(1.+kappa_j/kappa_jmo);
        Ai[i+1]=kappa_j*kappa_jpo/pow(kappa_jmo, 2);
      }
    }
    AArray.push_back(Ai);
  }
}

getAnalyticalIntegral()

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

pure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

getBArray()

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

protectedvirtual

Definition at line 72 of file MELANCSplineCore.cc.

                                                                                                                                                                                                                                 {
  BArray.clear();
  int npoints=kappas.size();
  if (npoints!=(int)fcnList.size()){
    coutE(InputArguments) << "MELANCSplineCore::getBArray: Dim(kappas)=" << npoints << " != Dim(fcnList)=" << fcnList.size() << endl;
    assert(0);
  }
  if (
    bcBegin==bcQuadraticWithNullSlope || bcEnd==bcQuadraticWithNullSlope
    ||
    bcBegin==bcQuadratic || bcEnd==bcQuadratic
    ){
    int nthr=1+(bcBegin==bcQuadraticWithNullSlope || bcBegin==bcQuadratic ? 1 : 0) + (bcEnd==bcQuadraticWithNullSlope || bcEnd==bcQuadratic ? 1 : 0);
    if (npoints<=nthr){
      cerr << "Npoints " << npoints << " <= " << nthr << endl;
      assert(0);
    }
  }
  if (npoints>1){
    MELANCSplineCore::T bcval=0, ecval=0;
    // First point constraint
    switch (bcBegin){
    case bcApproximatedSecondDerivative:
      bcval=(npoints<3 ? 0 : ((fcnList.at(2)-fcnList.at(1))*kappas.at(1)-(fcnList.at(1)-fcnList.at(0))*kappas.at(0))/(0.5/kappas.at(0)+0.5/kappas.at(1)));
    case bcNaturalSpline:
      BArray.push_back(3.*(fcnList.at(1)-fcnList.at(0)) - bcval/2./pow(kappas.at(0), 2));
      break;
    case bcApproximatedSlope:
      bcval=(fcnList.at(1)-fcnList.at(0))*kappas.at(0);
    case bcClamped:
      BArray.push_back(bcval/kappas.at(0));
      break;
    case bcQuadratic:
      bcval=2.*(fcnList.at(1)-fcnList.at(0));
      BArray.push_back(bcval);
      break;
    case bcQuadraticWithNullSlope:
      bcval=2.*(fcnList.at(1)-fcnList.at(0));
      BArray.push_back(0);
      BArray.push_back(bcval*kappas.at(0)/kappas.at(1));
      break;
    default:
      cerr << "MELANCSplineCore::getBArray: bcBegin " << bcBegin << " is not implemented!" << endl;
      assert(0);
    }
    // Intermediate point constraint is always D0, D1 and D2 continuous
    for (int j=1; j<npoints-1; j++){
      if (j==1 && bcBegin==bcQuadraticWithNullSlope) continue;
      if (j==npoints-2 && bcEnd==bcQuadraticWithNullSlope) continue;
      MELANCSplineCore::T val_j = fcnList.at(j);
      MELANCSplineCore::T val_jpo = fcnList.at(j+1);
      MELANCSplineCore::T val_jmo = fcnList.at(j-1);
      MELANCSplineCore::T kappa_j = kappas.at(j);
      MELANCSplineCore::T kappa_jmo = kappas.at(j-1);
      MELANCSplineCore::T rsq = pow(kappa_j/kappa_jmo, 2);
      MELANCSplineCore::T Bval = val_jpo*rsq + val_j*(1.-rsq) - val_jmo;
      Bval *= 3.;
      BArray.push_back(Bval);
    }
    // Last point constraint
    switch (bcEnd){
    case bcApproximatedSecondDerivative:
      ecval=(npoints<3 ? 0 : ((fcnList.at(npoints-1)-fcnList.at(npoints-2))*kappas.at(npoints-2)-(fcnList.at(npoints-2)-fcnList.at(npoints-3))*kappas.at(npoints-3))/(0.5/kappas.at(npoints-3)+0.5/kappas.at(npoints-2)));
    case bcNaturalSpline:
      BArray.push_back(3.*(fcnList.at(npoints-1)-fcnList.at(npoints-2)) + ecval/2./pow(kappas.at(npoints-1), 2));
      break;
    case bcApproximatedSlope:
      ecval=(fcnList.at(npoints-1)-fcnList.at(npoints-2))*kappas.at(npoints-2);
    case bcClamped:
      BArray.push_back(ecval/kappas.at(npoints-1));
      break;
    case bcQuadratic:
      ecval=2.*(fcnList.at(npoints-1)-fcnList.at(npoints-2));
      BArray.push_back(ecval);
      break;
    case bcQuadraticWithNullSlope:
      ecval=2.*(fcnList.at(npoints-1)-fcnList.at(npoints-2));
      BArray.push_back(ecval);
      BArray.push_back(0);
      break;
    default:
      cerr << "MELANCSplineCore::getAArray: bcEnd " << bcEnd << " is not implemented!" << endl;
      assert(0);
    }
  }
  else if (npoints==1) BArray.push_back(0);
}

getCoefficients()

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

protectedvirtual

Definition at line 219 of file MELANCSplineCore.cc.

                                                                                                                                                                                       {
  DefaultAccumulator<MELANCSplineCore::T> A, B, C, D;
  vector<MELANCSplineCore::T> res;
 
  const int fcnsize = fcnList.size();
  if (fcnsize>bin){
    A=fcnList.at(bin);
    B=S[bin];
    if (fcnsize>(bin+1)){
      DefaultAccumulator<MELANCSplineCore::T> dFcn = fcnList.at(bin+1);
      dFcn -= A;
 
      C += MELANCSplineCore::T(3.)*dFcn;
      C -= 2.*B;
      C -= S[bin+1]*kappas.at(bin+1)/kappas.at(bin);
 
      D += MELANCSplineCore::T(-2.)*dFcn;
      D += B;
      D += S[bin+1]*kappas.at(bin+1)/kappas.at(bin);
    }
  }
 
  res.push_back(A);
  res.push_back(B);
  res.push_back(C);
  res.push_back(D);
  return res;
}

getCoefficientsAlongDirection()

vector< vector< MELANCSplineCore::T > > MELANCSplineCore::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

protectedvirtual

Definition at line 247 of file MELANCSplineCore.cc.

                                                                                                                                                                                                                                                                                        {
  vector<MELANCSplineCore::T> BArray;
  getBArray(kappas, fcnList, BArray, bcBegin, bcEnd);
 
  Int_t npoints = BArray.size();
  TVector_t Btrans(npoints);
  for (int i=0; i<npoints; i++) Btrans[i]=BArray.at(i);
  TVector_t Strans = Ainv*Btrans;
 
  vector<vector<MELANCSplineCore::T>> coefs;
  for (Int_t bin=0; bin<(npoints>1 ? npoints-1 : 1); bin++){
    if (pickBin>=0 && bin!=pickBin) continue;
    vector<MELANCSplineCore::T> coef = getCoefficients(Strans, kappas, fcnList, bin);
    coefs.push_back(coef);
  }
  return coefs;
}

getKappas()

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

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

getLeafDependents()

void MELANCSplineCore::getLeafDependents ( RooRealProxy &  proxy, RooArgSet &  set  )

protected

Definition at line 275 of file MELANCSplineCore.cc.

                                                                           {
  RooArgSet deps;
  proxy.absArg()->leafNodeServerList(&deps, 0, true);
  set.add(deps);
}

getTVar()

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

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

getWhichBin()

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

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

interpolateFcn()

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

protectedpure virtual

Implemented in MELANCSpline_3D_fast, MELANCSpline_2D_fast, and MELANCSpline_1D_fast.

npointsX()

unsigned int MELANCSplineCore::npointsX ( ) const

inlineprotected

Definition at line 84 of file MELANCSplineCore.h.

{ return XList.size(); }

setEvalFloor()

void MELANCSplineCore::setEvalFloor

(

MELANCSplineCore::T

val

)

Definition at line 68 of file MELANCSplineCore.cc.

{ floorEval = val; }

setIntFloor()

void MELANCSplineCore::setIntFloor

(

MELANCSplineCore::T

val

)

Definition at line 69 of file MELANCSplineCore.cc.

{ floorInt = val; }

setRangeValidity()

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

pure virtual

Implemented in MELANCSpline_1D_fast, MELANCSpline_3D_fast, and MELANCSpline_2D_fast.

setVerbosity()

void MELANCSplineCore::setVerbosity ( VerbosityLevel  flag )

virtual

Definition at line 67 of file MELANCSplineCore.cc.

{ verbosity = flag; }

testRangeValidity()

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

protectedpure virtual

Implemented in MELANCSpline_1D_fast, MELANCSpline_3D_fast, and MELANCSpline_2D_fast.

Member Data Documentation

floorEval

T MELANCSplineCore::floorEval

protected

Definition at line 69 of file MELANCSplineCore.h.

floorInt

T MELANCSplineCore::floorInt

protected

Definition at line 70 of file MELANCSplineCore.h.

leafDepsList

RooListProxy MELANCSplineCore::leafDepsList

protected

Definition at line 76 of file MELANCSplineCore.h.

rangeXmax

T MELANCSplineCore::rangeXmax

protected

Definition at line 73 of file MELANCSplineCore.h.

rangeXmin

T MELANCSplineCore::rangeXmin

protected

Definition at line 72 of file MELANCSplineCore.h.

theXVar

RooRealProxy MELANCSplineCore::theXVar

protected

Definition at line 75 of file MELANCSplineCore.h.

useFloor

Bool_t MELANCSplineCore::useFloor

protected

Definition at line 68 of file MELANCSplineCore.h.

verbosity

VerbosityLevel MELANCSplineCore::verbosity

protected

Definition at line 67 of file MELANCSplineCore.h.

XList

std::vector<T> MELANCSplineCore::XList

protected

Definition at line 77 of file MELANCSplineCore.h.

---

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

• MELA/interface/MELANCSplineCore.h
• MELA/src/MELANCSplineCore.cc