BeamSpotPdf Class Reference

#include <BeamSpotPdf.h>

Inheritance diagram for BeamSpotPdf:

Collaboration diagram for BeamSpotPdf:

Public Member Functions
	BeamSpotPdf ()
	BeamSpotPdf (const char *name, const char *title, RooAbsReal &x, RooAbsReal &y, RooAbsReal &z, RooAbsReal &vxx, RooAbsReal &vyy, RooAbsReal &vxy, RooAbsReal &mx, RooAbsReal &sx, RooAbsReal &ax, RooAbsReal &my, RooAbsReal &sy, RooAbsReal &ay, RooAbsReal &mz, RooAbsReal &sz, RooAbsReal &k, RooAbsReal &rho)
	BeamSpotPdf (const BeamSpotPdf &other, const char *name=0)
virtual TObject *	clone (const char *newname) const
virtual	~BeamSpotPdf ()
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const
	BeamSpotPdf ()
	BeamSpotPdf (const char *name, const char *title, RooAbsReal &_x, RooAbsReal &_y, RooAbsReal &_z, RooAbsReal &_vxx, RooAbsReal &_vyy, RooAbsReal &_vxy, RooAbsReal &_mx, RooAbsReal &_sx, RooAbsReal &_ax, RooAbsReal &_my, RooAbsReal &_sy, RooAbsReal &_ay, RooAbsReal &_mz, RooAbsReal &_sz, RooAbsReal &_k, RooAbsReal &_rho)
	BeamSpotPdf (const BeamSpotPdf &other, const char *name=0)
virtual TObject *	clone (const char *newname) const
virtual	~BeamSpotPdf ()
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const

Protected Member Functions
Double_t	evaluate () const
Double_t	evaluate () const

Protected Attributes
RooRealProxy	m_x
RooRealProxy	m_y
RooRealProxy	m_z
RooRealProxy	m_vxx
RooRealProxy	m_vyy
RooRealProxy	m_vxy
RooRealProxy	m_mx
RooRealProxy	m_sx
RooRealProxy	m_ax
RooRealProxy	m_my
RooRealProxy	m_sy
RooRealProxy	m_ay
RooRealProxy	m_mz
RooRealProxy	m_sz
RooRealProxy	m_k
RooRealProxy	m_rho
RooRealProxy	x
RooRealProxy	y
RooRealProxy	z
RooRealProxy	vxx
RooRealProxy	vyy
RooRealProxy	vxy
RooRealProxy	mx
RooRealProxy	sx
RooRealProxy	ax
RooRealProxy	my
RooRealProxy	sy
RooRealProxy	ay
RooRealProxy	mz
RooRealProxy	sz
RooRealProxy	k
RooRealProxy	rho

Detailed Description

Definition at line 21 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

Constructor & Destructor Documentation

BeamSpotPdf() [1/6]

BeamSpotPdf::BeamSpotPdf ( )

inline

Definition at line 28 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

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

BeamSpotPdf() [2/6]

BeamSpotPdf::BeamSpotPdf	(	const char *	name,
		const char *	title,
		RooAbsReal &	x,
		RooAbsReal &	y,
		RooAbsReal &	z,
		RooAbsReal &	vxx,
		RooAbsReal &	vyy,
		RooAbsReal &	vxy,
		RooAbsReal &	mx,
		RooAbsReal &	sx,
		RooAbsReal &	ax,
		RooAbsReal &	my,
		RooAbsReal &	sy,
		RooAbsReal &	ay,
		RooAbsReal &	mz,
		RooAbsReal &	sz,
		RooAbsReal &	k,
		RooAbsReal &	rho
	)

Definition at line 23 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx.

               :
   RooAbsPdf(name,title), 
   m_x("x","x",this,x),
   m_y("y","y",this,y),
   m_z("z","z",this,z),
   m_vxx("vxx","vxx",this,vxx),
   m_vyy("vyy","vyy",this,vyy),
   m_vxy("vxy","vxy",this,vxy),
   m_mx("mx","mx",this,mx),
   m_sx("sx","sx",this,sx),
   m_ax("ax","ax",this,ax),
   m_my("my","my",this,my),
   m_sy("sy","sy",this,sy),
   m_ay("ay","ay",this,ay),
   m_mz("mz","mz",this,mz),
   m_sz("sz","sz",this,sz),
   m_k("k","k",this,k),
   m_rho("rho","rho",this,rho)
 { 
 } 

BeamSpotPdf() [3/6]

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

Definition at line 62 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx.

                                                                    :  
   RooAbsPdf(other,name), 
   m_x("x",this,other.m_x),
   m_y("y",this,other.m_y),
   m_z("z",this,other.m_z),
   m_vxx("vxx",this,other.m_vxx),
   m_vyy("vyy",this,other.m_vyy),
   m_vxy("vxy",this,other.m_vxy),
   m_mx("mx",this,other.m_mx),
   m_sx("sx",this,other.m_sx),
   m_ax("ax",this,other.m_ax),
   m_my("my",this,other.m_my),
   m_sy("sy",this,other.m_sy),
   m_ay("ay",this,other.m_ay),
   m_mz("mz",this,other.m_mz),
   m_sz("sz",this,other.m_sz),
   m_k("k",this,other.m_k),
   m_rho("rho",this,other.m_rho)
 { 
 } 

~BeamSpotPdf() [1/2]

virtual BeamSpotPdf::~BeamSpotPdf ( )

inlinevirtual

Definition at line 49 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

:
 

BeamSpotPdf() [4/6]

BeamSpotPdf::BeamSpotPdf ( )

inline

Definition at line 26 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

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

BeamSpotPdf() [5/6]

BeamSpotPdf::BeamSpotPdf	(	const char *	name,
		const char *	title,
		RooAbsReal &	_x,
		RooAbsReal &	_y,
		RooAbsReal &	_z,
		RooAbsReal &	_vxx,
		RooAbsReal &	_vyy,
		RooAbsReal &	_vxy,
		RooAbsReal &	_mx,
		RooAbsReal &	_sx,
		RooAbsReal &	_ax,
		RooAbsReal &	_my,
		RooAbsReal &	_sy,
		RooAbsReal &	_ay,
		RooAbsReal &	_mz,
		RooAbsReal &	_sz,
		RooAbsReal &	_k,
		RooAbsReal &	_rho
	)

BeamSpotPdf() [6/6]

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

~BeamSpotPdf() [2/2]

virtual BeamSpotPdf::~BeamSpotPdf ( )

inlinevirtual

Definition at line 46 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

:
 

Member Function Documentation

analyticalIntegral() [1/2]

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

analyticalIntegral() [2/2]

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

Definition at line 144 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx.

 { 
   // RETURN ANALYTICAL INTEGRAL DEFINED BY RETURN CODE ASSIGNED BY getAnalyticalIntegral
   // THE MEMBER FUNCTION m_x.min(rangeName) AND m_x.max(rangeName) WILL RETURN THE INTEGRATION
   // BOUNDARIES FOR EACH OBSERVABLE m_x
 
   double pi = TMath::Pi();
   double k2 = m_k*m_k;
   double wxx = m_sx*m_sx + k2*m_vxx;
   double wyy = m_sy*m_sy + k2*m_vyy;
   double wxy = m_rho*m_sx*m_sy + k2*m_vxy;
   double wzz = m_sz*m_sz;
   double detw = wxx*wyy - wxy*wxy;
  
   if (code==1) {
     //integrating over x,y,z
    
  //   return sqrt((0.5*pi)*(0.5*pi)*(0.5*pi))*sqrt(wxx)*sqrt(wyy)*m_sz*
  //     (TMath::Erf((m_mx-m_x.min(rangeName))/sqrt(2*wxx))-TMath::Erf((m_mx-m_x.max(rangeName))/sqrt(2*wxx)))
  //     *(TMath::Erf((m_my-m_y.min(rangeName))/sqrt(2*wyy))-TMath::Erf((m_my-m_y.max(rangeName))/sqrt(2*wyy)))
  //     *(TMath::Erf((m_mz-m_z.min(rangeName))/sqrt(2*wzz))-TMath::Erf((m_mz-m_z.max(rangeName))/sqrt(2*wzz)));
  //
     //The next line is used if we want to approximate the integral with that over -inf to +inf     
     return 2*pi*sqrt(detw)*sqrt(2*pi)*m_sz;
   } 
   
   
   if (code==2) {
     //integrating over x,y only
     double dz = m_z - m_mz;
     //erf approximation
    //
    //   return 0.5*pi*sqrt(wxx)*sqrt(wyy)*
    //   (TMath::Erf((m_mx-m_x.min(rangeName))/sqrt(2*wxx))-TMath::Erf((m_mx-m_x.max(rangeName))/sqrt(2*wxx)))
    //   *(TMath::Erf((m_my-m_y.min(rangeName))/sqrt(2*wyy))-TMath::Erf((m_my-m_y.max(rangeName))/sqrt(2*wyy)))
    //   *exp(-0.5*dz*dz/(m_sz*m_sz));
    //
       
     //inf approximation
     return 2*pi*sqrt(detw)*exp(-0.5*dz*dz/(m_sz*m_sz));
   }
   if (code==3){
     //integrating over x,z
     double dy = m_y-m_my;
     return 0.5*pi*sqrt(wxx)*sqrt(wzz)*
       (TMath::Erf((m_mx-m_x.min(rangeName))/sqrt(2*wxx))-TMath::Erf((m_mx-m_x.max(rangeName))/sqrt(2*wxx)))
       *(TMath::Erf((m_mz-m_z.min(rangeName))/sqrt(2*wzz))-TMath::Erf((m_mz-m_z.max(rangeName))/sqrt(2*wzz)))
       *exp(-0.5*dy*dy/(wyy));
   }
   if (code==4){
     //integrating over y,z
     double dx = m_x-m_mx;
     return 0.5*pi*sqrt(wyy)*sqrt(wzz)*
       (TMath::Erf((m_my-m_y.min(rangeName))/sqrt(2*wyy))-TMath::Erf((m_my-m_y.max(rangeName))/sqrt(2*wyy)))
       *(TMath::Erf((m_mz-m_z.min(rangeName))/sqrt(2*wzz))-TMath::Erf((m_mz-m_z.max(rangeName))/sqrt(2*wzz)))
       *exp(-0.5*dx*dx/(wxx));
   }
   
   return 0 ; 
 } 

clone() [1/2]

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

inlinevirtual

Definition at line 45 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

:
 

clone() [2/2]

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

inlinevirtual

Definition at line 48 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

:
 

evaluate() [1/2]

Double_t BeamSpotPdf::evaluate ( ) const

protected

evaluate() [2/2]

Double_t BeamSpotPdf::evaluate ( ) const

protected

Definition at line 85 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx.

 { 
   // ENTER EXPRESSION IN TERMS OF VARIABLE ARGUMENTS HERE 
   double dx = m_x - (m_mx + m_ax*(m_z-m_mz));
   double dy = m_y - (m_my + m_ay*(m_z-m_mz));
   double dz = m_z - m_mz;
   double k2 = m_k*m_k;
   double wxx = m_sx*m_sx + k2*m_vxx;
   double wyy = m_sy*m_sy + k2*m_vyy;
   double wxy = m_rho*m_sx*m_sy + k2*m_vxy;
   double detw = wxx*wyy - wxy*wxy;
   double iwxx = wyy/detw;
   double iwyy = wxx/detw;
   double iwxy = -wxy/detw;
   double arg = dx*iwxx*dx + 2*dx*iwxy*dy + dy*iwyy*dy;
   arg += dz*dz/(m_sz*m_sz);
   return exp(-0.5*arg);
 }

getAnalyticalIntegral() [1/2]

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

Definition at line 104 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx.

 { 
   //Previously, we checked to see if the integration range was larger than 5 sigma from the mean, and used 
   //the integral from -inf to +inf to approximate it. Uncomment the next block of code to implement that procedure:
   
   
 //  int mult = 5;
// 
  // if (matchArgs(allVars,analVars,m_x,m_y,m_z)){
    // if ((m_x.min(rangeName) < m_mx - mult*m_sx) && (m_x.max(rangeName) > m_mx + mult*m_sx)){
      // if ((m_y.min(rangeName) < m_my - mult*m_sy) && (m_y.max(rangeName) > m_my + mult*m_sy)){
    // if ((m_z.min(rangeName) < m_mz - mult*m_sz) && (m_z.max(rangeName) > m_mz + mult*m_sz)){
//     return 1 ;
//   }
//       }
//     }
//     
//   }
//   if (matchArgs(allVars,analVars,m_x,m_y)){
//     if ((m_x.min(rangeName) < m_mx - mult*m_sx ) && (m_x.max(rangeName) > m_mx + mult*m_sx)){
//       if ((m_y.min(rangeName) < m_my - mult*m_sy) && (m_y.max(rangeName) > m_my + mult*m_sy)){
//   return 2 ; 
//       }
//     }
//     }     
   
 
   //Since we don't care about the range of integration, we aren't using this rangeName variable.
   //The following line does nothing but silence the compiler warning.
   //(void) rangeName;
 
   //Now we can use the Erf approximation over the entire integration range:
     if (matchArgs(allVars,analVars,m_x,m_y,m_z)) return 1;
     if (matchArgs(allVars,analVars,m_x,m_y) ) return 2;
     if (matchArgs(allVars,analVars,m_x,m_z) ) return 3;
     if (matchArgs(allVars,analVars,m_y,m_z) ) return 4;
   
   return 0 ; 
 } 

getAnalyticalIntegral() [2/2]

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

Member Data Documentation

ax

RooRealProxy BeamSpotPdf::ax

protected

Definition at line 61 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

ay

RooRealProxy BeamSpotPdf::ay

protected

Definition at line 64 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

k

RooRealProxy BeamSpotPdf::k

protected

Definition at line 67 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

m_ax

RooRealProxy BeamSpotPdf::m_ax

protected

Definition at line 64 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_ay

RooRealProxy BeamSpotPdf::m_ay

protected

Definition at line 67 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_k

RooRealProxy BeamSpotPdf::m_k

protected

Definition at line 70 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_mx

RooRealProxy BeamSpotPdf::m_mx

protected

Definition at line 62 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_my

RooRealProxy BeamSpotPdf::m_my

protected

Definition at line 65 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_mz

RooRealProxy BeamSpotPdf::m_mz

protected

Definition at line 68 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_rho

RooRealProxy BeamSpotPdf::m_rho

protected

Definition at line 71 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_sx

RooRealProxy BeamSpotPdf::m_sx

protected

Definition at line 63 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_sy

RooRealProxy BeamSpotPdf::m_sy

protected

Definition at line 66 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_sz

RooRealProxy BeamSpotPdf::m_sz

protected

Definition at line 69 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_vxx

RooRealProxy BeamSpotPdf::m_vxx

protected

Definition at line 59 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_vxy

RooRealProxy BeamSpotPdf::m_vxy

protected

Definition at line 61 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_vyy

RooRealProxy BeamSpotPdf::m_vyy

protected

Definition at line 60 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_x

RooRealProxy BeamSpotPdf::m_x

protected

Definition at line 56 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_y

RooRealProxy BeamSpotPdf::m_y

protected

Definition at line 57 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

m_z

RooRealProxy BeamSpotPdf::m_z

protected

Definition at line 58 of file InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h.

mx

RooRealProxy BeamSpotPdf::mx

protected

Definition at line 59 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

my

RooRealProxy BeamSpotPdf::my

protected

Definition at line 62 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

mz

RooRealProxy BeamSpotPdf::mz

protected

Definition at line 65 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

rho

RooRealProxy BeamSpotPdf::rho

protected

Definition at line 68 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

sx

RooRealProxy BeamSpotPdf::sx

protected

Definition at line 60 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

sy

RooRealProxy BeamSpotPdf::sy

protected

Definition at line 63 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

sz

RooRealProxy BeamSpotPdf::sz

protected

Definition at line 66 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

vxx

RooRealProxy BeamSpotPdf::vxx

protected

Definition at line 56 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

vxy

RooRealProxy BeamSpotPdf::vxy

protected

Definition at line 58 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

vyy

RooRealProxy BeamSpotPdf::vyy

protected

Definition at line 57 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

x

RooRealProxy BeamSpotPdf::x

protected

Definition at line 53 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

y

RooRealProxy BeamSpotPdf::y

protected

Definition at line 54 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

z

RooRealProxy BeamSpotPdf::z

protected

Definition at line 55 of file InDetExample/InDetBeamSpotExample/roofit/BeamSpotPdf.h.

---

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

• InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.h
• InDetCalibAlgs/InDetBeamSpotFinder/src/BeamSpotPdf.cxx