RooSplineProduct Class Reference

#include <RooSplineProduct.h>

Inheritance diagram for RooSplineProduct:

Public Member Functions
	RooSplineProduct ()
	RooSplineProduct (const char *name, const char *title, RooRealVar &x, const RooCubicSplineFun &sp1, const RooCubicSplineFun &sp2)
	~RooSplineProduct ()
	RooSplineProduct (const RooSplineProduct &other, const char *name=0)
TObject *	clone (const char *newname) const
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName) const
std::complex< double >	productAnalyticalIntegral (Double_t umin, Double_t umax, Double_t scale, Double_t offset, const std::complex< double > &z) const
Public Member Functions inherited from RooAbsGaussModelEfficiency
	RooAbsGaussModelEfficiency ()
	RooAbsGaussModelEfficiency (const char *name, const char *title, const char *unit="")
	RooAbsGaussModelEfficiency (const RooAbsGaussModelEfficiency &other, const char *name=0)
	~RooAbsGaussModelEfficiency ()

Private Member Functions
void	init ()
Double_t	evaluate () const
std::complex< double >	gaussIntegralE (bool left, const RooGaussModelAcceptance::M_n< 7U > &dM, const RooGaussModelAcceptance::K_n &K, double offset, double *sc) const

Private Attributes
RooRealProxy	_x
RooCubicSplineFun	_sp1
RooCubicSplineFun	_sp2
RooListProxy	_coefList1
RooListProxy	_coefList2

Detailed Description

Definition at line 25 of file RooSplineProduct.h.

Constructor & Destructor Documentation

RooSplineProduct() [1/3]

RooSplineProduct::RooSplineProduct

(

)

Definition at line 41 of file RooSplineProduct.cpp.

{
}

RooSplineProduct() [2/3]

RooSplineProduct::RooSplineProduct	(	const char *	name,
		const char *	title,
		RooRealVar &	x,
		const RooCubicSplineFun &	sp1,
		const RooCubicSplineFun &	sp2
	)

Definition at line 46 of file RooSplineProduct.cpp.

                                   :
   RooAbsGaussModelEfficiency(name, title),
   _x("x", "Dependent", this, x),
   _sp1(sp1,"spline1"),
   _sp2(sp2,"spline2"),
   _coefList1("coefficients1","List of coefficients for first spline",this),
   _coefList2("coefficients2","List of coefficients for second spline",this)
{
  _coefList1.add(sp1.coefficients());
  _coefList2.add(sp2.coefficients());
  init();
}

~RooSplineProduct()

RooSplineProduct::~RooSplineProduct

(

)

Definition at line 75 of file RooSplineProduct.cpp.

{
}

RooSplineProduct() [3/3]

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

Definition at line 64 of file RooSplineProduct.cpp.

                                                                                  :
  RooAbsGaussModelEfficiency(other, name),
  _x("x", this, other._x),
  _sp1(other._sp1,"spline1"),
  _sp2(other._sp2,"spline2"),
  _coefList1("coefList1", this, other._coefList1),
  _coefList2("coefList2", this, other._coefList2)
{
}

Member Function Documentation

analyticalIntegral()

Double_t RooSplineProduct::analyticalIntegral	(	Int_t	code,
		const char *	rangeName
	)		const

Definition at line 91 of file RooSplineProduct.cpp.

{
  return 0;
}

clone()

TObject* RooSplineProduct::clone ( const char *  newname ) const

inline

Definition at line 33 of file RooSplineProduct.h.

{ return new RooSplineProduct(*this, newname); }

evaluate()

Double_t RooSplineProduct::evaluate ( ) const

private

Definition at line 80 of file RooSplineProduct.cpp.

{
  return _sp1._aux.evaluate(_x, _coefList1) * _sp2._aux.evaluate(_x, _coefList2);
}

gaussIntegralE()

complex< double > RooSplineProduct::gaussIntegralE ( bool  left, const RooGaussModelAcceptance::M_n< 7U > &  dM, const RooGaussModelAcceptance::K_n &  K, double  offset, double *  sc  ) const

private

Definition at line 97 of file RooSplineProduct.cpp.

{
  RooCubicSplineKnot::S2_edge s2_jk( _sp1._aux.S2_jk_edge( left, _coefList1, _coefList2 ), offset );
  std::complex<double> sum(0,0);
  for (int j=0;j<3;++j) for (int k=0;k<3-j;++k) sum += dM(j)*s2_jk(j,k)*K(k)*sc[j+k];
  return sum;
}

getAnalyticalIntegral()

Int_t RooSplineProduct::getAnalyticalIntegral	(	RooArgSet &	allVars,
		RooArgSet &	analVars,
		const char *	rangeName
	)		const

Definition at line 86 of file RooSplineProduct.cpp.

{
  return 0;
}

init()

void RooSplineProduct::init ( )

private

Definition at line 33 of file RooSplineProduct.cpp.

                            {
    // TODO: current implementation of product of two splines
    // only for splines with the same knot vector
    assert(_sp1.knots() == _sp2.knots());
    assert(_sp1.knotSize() == _sp2.knotSize());
}

productAnalyticalIntegral()

std::complex< double > RooSplineProduct::productAnalyticalIntegral ( Double_t  umin, Double_t  umax, Double_t  scale, Double_t  offset, const std::complex< double > &  z  ) const

virtual

Implements RooAbsGaussModelEfficiency.

Definition at line 109 of file RooSplineProduct.cpp.

{
    RooGaussModelAcceptance::K_n K(z);
    assert(_sp1.knotSize()>1);
    double lo = scale*umin+offset;
    double hi = scale*umax+offset+1.e-7;
    typedef RooGaussModelAcceptance::M_n<7U> M_n;
    std::vector<M_n> M; M.reserve( _sp1.knotSize() );
    for (unsigned int i=0;i<_sp1.knotSize();++i) {
        double x = (_sp1.u(i)-offset)/scale ;
        if (lo>=_sp1.u(i)) x = umin ; // take M[i] if lo<=u(i) else M_n(lo)
        if (_sp1.u(i)>=hi) x = umax ; // take M[i+1] if u(i+1)<=hi else M_n(hi)
        M.push_back( M_n( x, z ) );
    }
    double sc[7]; for (int i=0;i<7;++i) sc[i] = pow(scale,i);
    std::complex<double> sum(0,0);
    if (lo<_sp1.u(0)) sum += gaussIntegralE(true,  M.front()-M_n( umin,z), K, offset, sc);
    for (unsigned i=0;i<_sp1.knotSize()-1 && _sp1.u(i)<hi ;++i) {
        M_n dM = M[i+1]-M[i];
        RooCubicSplineKnot::S2_jk s2_jk( _sp1._aux.S2_jk_sum(i, _coefList1, _coefList2), offset );
        for (int j=0;j<7;++j) for (int k=0;k<7-j;++k) sum += dM(j)*s2_jk(j,k)*K(k)*sc[j+k];
    }
    if (hi>_sp1.u(_sp1.knotSize()-1)) sum += gaussIntegralE(false, M_n(umax,z)-M.back(), K, offset, sc);
    return sum;
}

Member Data Documentation

_coefList1

RooListProxy RooSplineProduct::_coefList1

private

Definition at line 48 of file RooSplineProduct.h.

_coefList2

RooListProxy RooSplineProduct::_coefList2

private

Definition at line 49 of file RooSplineProduct.h.

_sp1

RooCubicSplineFun RooSplineProduct::_sp1

private

Definition at line 46 of file RooSplineProduct.h.

_sp2

RooCubicSplineFun RooSplineProduct::_sp2

private

Definition at line 47 of file RooSplineProduct.h.

_x

RooRealProxy RooSplineProduct::_x

private

Definition at line 45 of file RooSplineProduct.h.

---

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

• Mint/RooSplineProduct.h
• src/Mojito/TimeDependent/RooSplineProduct.cpp