GenGauss3D Class Reference

#include <GenGauss3D.h>

Inheritance diagram for GenGauss3D:

Collaboration diagram for GenGauss3D:

Public Member Functions
	GenGauss3D ()
	GenGauss3D (const char *name, const char *title, RooArgList &_xvec, RooArgList &_muvex, RooAbsReal &_vxx, RooAbsReal &_vxy, RooAbsReal &_vxz, RooAbsReal &_vyy, RooAbsReal &_vyz, RooAbsReal &_vzz, TMatrixDSym &dummyMatrix)
	GenGauss3D (const GenGauss3D &other, const char *name=0)
virtual TObject *	clone (const char *newname) const
virtual	~GenGauss3D ()
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const
void	generateEvent (Int_t code)
void	updateCov ()

Protected Member Functions
Double_t	evaluate () const

Protected Attributes
RooRealProxy	xvec
RooRealProxy	muvec
RooRealProxy	vxx
RooRealProxy	vxy
RooRealProxy	vxz
RooRealProxy	vyy
RooRealProxy	vyz
RooRealProxy	vzz

Detailed Description

Definition at line 22 of file GenGauss3D.h.

Constructor & Destructor Documentation

GenGauss3D() [1/3]

GenGauss3D::GenGauss3D ( )

inline

Definition at line 24 of file GenGauss3D.h.

{} ; 

GenGauss3D() [2/3]

GenGauss3D::GenGauss3D	(	const char *	name,
		const char *	title,
		RooArgList &	_xvec,
		RooArgList &	_muvex,
		RooAbsReal &	_vxx,
		RooAbsReal &	_vxy,
		RooAbsReal &	_vxz,
		RooAbsReal &	_vyy,
		RooAbsReal &	_vyz,
		RooAbsReal &	_vzz,
		TMatrixDSym &	dummyMatrix )

GenGauss3D() [3/3]

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

Definition at line 50 of file GenGauss3D.cxx.

                                                                 :  
   RooMultiVarGaussian(other,name), 
   vxx("vxx",this,other.vxx),
   vxy("vxy",this,other.vxy),
   vxz("vxz",this,other.vxz),
   vyy("vyy",this,other.vyy),
   vyz("vyz",this,other.vyz),
   vzz("vzz",this,other.vzz)
   //nUpdates(0),
   //nChanged(0),
   //oldVxx(-1),
   //oldVxy(-1),
   //oldVyy(-1),
   //oldVzz(-1)
 { 
 } 

~GenGauss3D()

virtual GenGauss3D::~GenGauss3D ( )

inlinevirtual

Definition at line 37 of file GenGauss3D.h.

{ }

Member Function Documentation

analyticalIntegral()

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

Definition at line 139 of file GenGauss3D.cxx.

 { 
   // RETURN ANALYTICAL INTEGRAL DEFINED BY RETURN CODE ASSIGNED BY getAnalyticalIntegral
   // THE MEMBER FUNCTION x.min(rangeName) AND x.max(rangeName) WILL RETURN THE INTEGRATION
   // BOUNDARIES FOR EACH OBSERVABLE x
 
   // assert(code==1) ; 
   // return (x.max(rangeName)-x.min(rangeName)) ; 
 
    // Since we need to update _cov here, must get rid of const
   GenGauss3D* self =  const_cast<GenGauss3D*>(this);
   self->updateCov();
 
   // Prevent caching of partial integrals - these may depend on the
   // particular covariance matrix!
   _anaIntCache.clear();
 
   return RooMultiVarGaussian::analyticalIntegral(code,rangeName) ; 
} 

clone()

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

inlinevirtual

Definition at line 36 of file GenGauss3D.h.

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

evaluate()

Double_t GenGauss3D::evaluate ( ) const

protected

Definition at line 107 of file GenGauss3D.cxx.

 { 
   // ENTER EXPRESSION IN TERMS OF VARIABLE ARGUMENTS HERE 
 
   // Since we need to update _cov here, must get rid of const
   GenGauss3D* self =  const_cast<GenGauss3D*>(this);
   self->updateCov();
 
   return RooMultiVarGaussian::evaluate();
 } 

generateEvent()

void GenGauss3D::generateEvent

(

Int_t

code

)

Definition at line 160 of file GenGauss3D.cxx.

                                          {
   _genCache.clear() ;
   updateCov();
   _genCache.clear() ;
   return RooMultiVarGaussian::generateEvent(code);
 }

getAnalyticalIntegral()

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

Definition at line 120 of file GenGauss3D.cxx.

 { 
   // LIST HERE OVER WHICH VARIABLES ANALYTICAL INTEGRATION IS SUPPORTED, 
   // ASSIGN A NUMERIC CODE FOR EACH SUPPORTED (SET OF) PARAMETERS 
   // THE EXAMPLE BELOW ASSIGNS CODE 1 TO INTEGRATION OVER VARIABLE X
   // YOU CAN ALSO IMPLEMENT MORE THAN ONE ANALYTICAL INTEGRAL BY REPEATING THE matchArgs 
   // EXPRESSION MULTIPLE TIMES
 
   // if (matchArgs(allVars,analVars,x)) return 1 ; 
 
   // Since we need to update _cov here, must get rid of const
   GenGauss3D* self =  const_cast<GenGauss3D*>(this);
   self->updateCov();
 
   return RooMultiVarGaussian::getAnalyticalIntegral(allVars,analVars,rangeName) ; 
 } 

updateCov()

void GenGauss3D::updateCov

(

)

Definition at line 68 of file GenGauss3D.cxx.

                            {
   //nUpdates++;
   //if (vxx==oldVxx &&
   //    vyy==oldVyy &&
   //    vxy==oldVxy &&
   //    vzz==oldVzz)
   //  return;
   //oldVxx = vxx;
   //oldVyy = vyy;
   //oldVxy = vxy;
   //oldVzz = vzz;
 
   //nChanged++;
   //if (nUpdates % 10000 == 0) {
   //  cout << nUpdates << " calls to updateCov" << endl;
   //  cout << nChanged << " calls required recalculation of of covariance matrix" << endl << endl;
   //}
 
   _cov(0,0) = vxx;
   _cov(0,1) = vxy;
   _cov(0,2) = vxz;
   _cov(1,0) = vxy;
   _cov(1,1) = vyy;
   _cov(1,2) = vyz;
   _cov(2,0) = vxz;
   _cov(2,1) = vyz;
   _cov(2,2) = vzz;
   _det = _cov.Determinant();
 
   int n = _cov.GetNrows();
   for (int i=0;i<n;i++) {
     for  (int j=0;j<n;j++) {
       _covI(i,j)=_cov(i,j);
     }
   }
   _covI.Invert();
 }

Member Data Documentation

muvec

RooRealProxy GenGauss3D::muvec

protected

Definition at line 49 of file GenGauss3D.h.

vxx

RooRealProxy GenGauss3D::vxx

protected

Definition at line 50 of file GenGauss3D.h.

vxy

RooRealProxy GenGauss3D::vxy

protected

Definition at line 51 of file GenGauss3D.h.

vxz

RooRealProxy GenGauss3D::vxz

protected

Definition at line 52 of file GenGauss3D.h.

vyy

RooRealProxy GenGauss3D::vyy

protected

Definition at line 53 of file GenGauss3D.h.

vyz

RooRealProxy GenGauss3D::vyz

protected

Definition at line 54 of file GenGauss3D.h.

vzz

RooRealProxy GenGauss3D::vzz

protected

Definition at line 55 of file GenGauss3D.h.

xvec

RooRealProxy GenGauss3D::xvec

protected

Definition at line 48 of file GenGauss3D.h.

---

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

• GenGauss3D.h
• GenGauss3D.cxx