#include <MELANCSplineFactory_3D.h>

Collaboration diagram for MELANCSplineFactory_3D:

[legend]

Public Member Functions
	MELANCSplineFactory_3D (RooAbsReal &XVar_, RooAbsReal &YVar_, RooAbsReal &ZVar_, TString appendName_="", 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)

	~MELANCSplineFactory_3D ()

MELANCSpline_3D_fast *	getFunc ()

MELAFuncPdf *	getPDF ()

void	setEndConditions (MELANCSplineCore::BoundaryCondition const bcBegin, MELANCSplineCore::BoundaryCondition const bcEnd, const unsigned int direction)

void	setPoints (TTree *tree)

void	setPoints (const std::vector< splineQuadruplet_t > &pList)

template<typename inType >
void	setPoints (const std::vector< inType > &XList, const std::vector< inType > &YList, const std::vector< inType > &ZList, const std::vector< inType > &FcnList)

Protected Member Functions
const std::vector< splineQuadruplet_t >	getPoints (const std::vector< MELANCSplineCore::T > &XList, const std::vector< MELANCSplineCore::T > &YList, const std::vector< MELANCSplineCore::T > &ZList, const std::vector< MELANCSplineCore::T > &FcnList)

void	destroyPDF ()

void	initPDF (const std::vector< splineQuadruplet_t > &pList)

void	addUnique (std::vector< MELANCSplineCore::T > &list, MELANCSplineCore::T val)

Protected Attributes
TString	appendName

MELANCSplineCore::BoundaryCondition	bcBeginX

MELANCSplineCore::BoundaryCondition	bcEndX

MELANCSplineCore::BoundaryCondition	bcBeginY

MELANCSplineCore::BoundaryCondition	bcEndY

MELANCSplineCore::BoundaryCondition	bcBeginZ

MELANCSplineCore::BoundaryCondition	bcEndZ

RooAbsReal *	XVar

RooAbsReal *	YVar

RooAbsReal *	ZVar

MELANCSpline_3D_fast *	fcn

MELAFuncPdf *	PDF

Detailed Description

Definition at line 14 of file MELANCSplineFactory_3D.h.

Constructor & Destructor Documentation

◆ MELANCSplineFactory_3D()

MELANCSplineFactory_3D::MELANCSplineFactory_3D	(	RooAbsReal &	XVar_,
		RooAbsReal &	YVar_,
		RooAbsReal &	ZVar_,
		TString	appendName_ = `""`,
		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`
	)

Definition at line 6 of file MELANCSplineFactory_3D.cc.

   :
   appendName(appendName_),
   bcBeginX(bcBeginX_), bcEndX(bcEndX_),
   bcBeginY(bcBeginY_), bcEndY(bcEndY_),
   bcBeginZ(bcBeginZ_), bcEndZ(bcEndZ_),
   XVar(&XVar_), YVar(&YVar_), ZVar(&ZVar_),
   fcn(0),
   PDF(0)
 {}

◆ ~MELANCSplineFactory_3D()

MELANCSplineFactory_3D::~MELANCSplineFactory_3D ( )

Definition at line 23 of file MELANCSplineFactory_3D.cc.

                                                {
   destroyPDF();
 }

Member Function Documentation

◆ addUnique()

void MELANCSplineFactory_3D::addUnique	(	std::vector< MELANCSplineCore::T > &	list,
		MELANCSplineCore::T	val
	)

protected

Definition at line 27 of file MELANCSplineFactory_3D.cc.

                                                                                                {
   for (unsigned int ip=0; ip<list.size(); ip++){ if (list.at(ip)==val) return; }
   list.push_back(val);
 }

◆ destroyPDF()

void MELANCSplineFactory_3D::destroyPDF ( )

protected

Definition at line 61 of file MELANCSplineFactory_3D.cc.

61 { delete PDF; PDF=0; delete fcn; fcn=0; }

◆ getFunc()

MELANCSpline_3D_fast* MELANCSplineFactory_3D::getFunc ( )

inline

Definition at line 50 of file MELANCSplineFactory_3D.h.

50 { return fcn; }

◆ getPDF()

MELAFuncPdf* MELANCSplineFactory_3D::getPDF ( )

inline

Definition at line 51 of file MELANCSplineFactory_3D.h.

51 { return PDF; }

◆ getPoints()

const std::vector< splineQuadruplet_t > MELANCSplineFactory_3D::getPoints	(	const std::vector< MELANCSplineCore::T > &	XList,
		const std::vector< MELANCSplineCore::T > &	YList,
		const std::vector< MELANCSplineCore::T > &	ZList,
		const std::vector< MELANCSplineCore::T > &	FcnList
	)

protected

Definition at line 31 of file MELANCSplineFactory_3D.cc.

    {
   const unsigned int nX = XList.size();
   const unsigned int nY = YList.size();
   const unsigned int nZ = ZList.size();
   const unsigned int n = FcnList.size();
   if (nX*nY*nZ!=n){
     cerr << "MELANCSplineFactory_3D::getPoints: nX=" << nX << " x nY=" << nY << " x nZ=" << nZ << " != nFcn=" << n << endl;
     assert(0);
   }
  
   std::vector<splineQuadruplet_t> pList; pList.reserve(n);
   for (unsigned int ix=0; ix<nX; ix++){
     MELANCSplineCore::T xval = XList.at(ix);
     for (unsigned int iy=0; iy<nY; iy++){
       MELANCSplineCore::T yval = YList.at(iy);
       for (unsigned int iz=0; iz<nZ; iz++){
         MELANCSplineCore::T zval = ZList.at(iz);
         unsigned int ip = nZ*(nY*ix + iy) + iz;
         pList.push_back(splineQuadruplet_t(xval, yval, zval, FcnList.at(ip)));
       }
     }
   }
   return pList;
 }

◆ initPDF()

void MELANCSplineFactory_3D::initPDF ( const std::vector< splineQuadruplet_t > & pList )

protected

Definition at line 62 of file MELANCSplineFactory_3D.cc.

                                                                               {
   destroyPDF();
  
   const unsigned int n = pList.size();
   vector<MELANCSplineCore::T> XList;
   vector<MELANCSplineCore::T> YList;
   vector<MELANCSplineCore::T> ZList;
   vector<vector<vector<MELANCSplineCore::T>>> FcnList;
   for (unsigned int ip=0; ip<n; ip++){
     addUnique(XList, (pList.at(ip))[0]);
     addUnique(YList, (pList.at(ip))[1]);
     addUnique(ZList, (pList.at(ip))[2]);
   }
   FcnList.reserve(ZList.size());
   for (unsigned int iz=0; iz<ZList.size(); iz++){
     vector<vector<MELANCSplineCore::T>> dumz;
     dumz.reserve(YList.size());
     for (unsigned int iy=0; iy<YList.size(); iy++){
       vector<MELANCSplineCore::T> dumy;
       dumy.reserve(XList.size());
       for (unsigned int ix=0; ix<XList.size(); ix++){
         unsigned int ip = ZList.size()*(YList.size()*ix + iy) + iz;
         dumy.push_back((pList.at(ip))[3]); // Do not use unique here
       }
       dumz.push_back(dumy);
     }
     FcnList.push_back(dumz);
   }
  
   TString name = "Func";
   if (appendName!="") name = Form("%s_%s", name.Data(), appendName.Data());
   TString title=name;
   fcn = new MELANCSpline_3D_fast(
     name.Data(),
     title.Data(),
     *XVar, *YVar, *ZVar,
     XList, YList, ZList, FcnList,
     bcBeginX, bcEndX,
     bcBeginY, bcEndY,
     bcBeginZ, bcEndZ
   );
  
   name.Prepend("PDF_"); title=name;
   PDF = new MELAFuncPdf(
     name.Data(),
     title.Data(),
     *fcn
     );
 }

◆ setEndConditions()

void MELANCSplineFactory_3D::setEndConditions	(	MELANCSplineCore::BoundaryCondition const	bcBegin,
		MELANCSplineCore::BoundaryCondition const	bcEnd,
		const unsigned int	direction
	)

Definition at line 124 of file MELANCSplineFactory_3D.cc.

  {
   switch (direction){
   case 0:
     bcBeginX=bcBegin;
     bcEndX=bcEnd;
     break;
   case 1:
     bcBeginY=bcBegin;
     bcEndY=bcEnd;
     break;
   case 2:
     bcBeginZ=bcBegin;
     bcEndZ=bcEnd;
     break;
   default:
     // Do nothing
     break;
   }
 }

◆ setPoints() [1/3]

template<typename inType >

void MELANCSplineFactory_3D::setPoints	(	const std::vector< inType > &	XList,
		const std::vector< inType > &	YList,
		const std::vector< inType > &	ZList,
		const std::vector< inType > &	FcnList
	)

inline

Definition at line 61 of file MELANCSplineFactory_3D.h.

                                                                                                                                                                             {
     std::vector<MELANCSplineCore::T> transXList;
     std::vector<MELANCSplineCore::T> transYList;
     std::vector<MELANCSplineCore::T> transZList;
     std::vector<MELANCSplineCore::T> transFcnList;
     for (unsigned int ip=0; ip<XList.size(); ip++) transXList.push_back((MELANCSplineCore::T)XList.at(ip));
     for (unsigned int ip=0; ip<YList.size(); ip++) transYList.push_back((MELANCSplineCore::T)YList.at(ip));
     for (unsigned int ip=0; ip<ZList.size(); ip++) transZList.push_back((MELANCSplineCore::T)ZList.at(ip));
     for (unsigned int ip=0; ip<FcnList.size(); ip++) transFcnList.push_back((MELANCSplineCore::T)FcnList.at(ip));
     const std::vector<splineQuadruplet_t> pList = getPoints(transXList, transYList, transZList, transFcnList);
     setPoints(pList);
   }

◆ setPoints() [2/3]

void MELANCSplineFactory_3D::setPoints ( const std::vector< splineQuadruplet_t > & pList )

inline

Definition at line 60 of file MELANCSplineFactory_3D.h.

60 { initPDF(pList); }

◆ setPoints() [3/3]

void MELANCSplineFactory_3D::setPoints ( TTree * tree )

Definition at line 112 of file MELANCSplineFactory_3D.cc.

                                                  {
   vector<splineQuadruplet_t> pList;
   MELANCSplineCore::T x, y, z, fcn;
   tree->SetBranchAddress("X", &x);
   tree->SetBranchAddress("Y", &y);
   tree->SetBranchAddress("Z", &z);
   tree->SetBranchAddress("Fcn", &fcn);
   int n = tree->GetEntries();
   for (int ip=0; ip<n; ip++){ tree->GetEntry(ip); pList.push_back(splineQuadruplet_t(x, y, z, fcn)); }
   setPoints(pList);
 }

Member Data Documentation

◆ appendName

TString MELANCSplineFactory_3D::appendName

protected

Definition at line 16 of file MELANCSplineFactory_3D.h.

◆ bcBeginX

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcBeginX

protected

Definition at line 18 of file MELANCSplineFactory_3D.h.

◆ bcBeginY

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcBeginY

protected

Definition at line 20 of file MELANCSplineFactory_3D.h.

◆ bcBeginZ

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcBeginZ

protected

Definition at line 22 of file MELANCSplineFactory_3D.h.

◆ bcEndX

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcEndX

protected

Definition at line 19 of file MELANCSplineFactory_3D.h.

◆ bcEndY

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcEndY

protected

Definition at line 21 of file MELANCSplineFactory_3D.h.

◆ bcEndZ

MELANCSplineCore::BoundaryCondition MELANCSplineFactory_3D::bcEndZ

protected

Definition at line 23 of file MELANCSplineFactory_3D.h.

◆ fcn

MELANCSpline_3D_fast* MELANCSplineFactory_3D::fcn

protected

Definition at line 28 of file MELANCSplineFactory_3D.h.

◆ PDF

MELAFuncPdf* MELANCSplineFactory_3D::PDF

protected

Definition at line 29 of file MELANCSplineFactory_3D.h.

◆ XVar

RooAbsReal* MELANCSplineFactory_3D::XVar

protected

Definition at line 25 of file MELANCSplineFactory_3D.h.

◆ YVar

RooAbsReal* MELANCSplineFactory_3D::YVar

protected

Definition at line 26 of file MELANCSplineFactory_3D.h.

◆ ZVar

RooAbsReal* MELANCSplineFactory_3D::ZVar

protected

Definition at line 27 of file MELANCSplineFactory_3D.h.

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

MELA/interface/MELANCSplineFactory_3D.h
MELA/src/MELANCSplineFactory_3D.cc

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MELANCSplineFactory_3D()

◆ ~MELANCSplineFactory_3D()

Member Function Documentation

◆ addUnique()

◆ destroyPDF()

◆ getFunc()

◆ getPDF()

◆ getPoints()

◆ initPDF()

◆ setEndConditions()

◆ setPoints() [1/3]

◆ setPoints() [2/3]

◆ setPoints() [3/3]

Member Data Documentation

◆ appendName

◆ bcBeginX

◆ bcBeginY

◆ bcBeginZ

◆ bcEndX

◆ bcEndY

◆ bcEndZ

◆ fcn

◆ PDF

◆ XVar

◆ YVar

◆ ZVar