CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
BSpdfsFcn Class Reference

#include <BSpdfsFcn.h>

Inheritance diagram for BSpdfsFcn:

Public Member Functions

virtual double operator() (const std::vector< double > &) const
 
void SetData (std::vector< BSTrkParameters > a_BSvector)
 
void SetPDFs (std::string usepdfs)
 
virtual double Up () const
 

Private Member Functions

double PDFGauss_d (double z, double d, double sigmad, double phi, const std::vector< double > &parms) const
 
double PDFGauss_d_resolution (double z, double d, double phi, double pt, const std::vector< double > &parms) const
 
double PDFGauss_z (double z, double sigmaz, const std::vector< double > &parms) const
 

Private Attributes

std::vector< BSTrkParameters > fBSvector
 
std::string fusepdfs
 

Static Private Attributes

static const int fPar_c0 = 7
 
static const int fPar_c1 = 8
 
static const int fPar_dxdz = 4
 
static const int fPar_dydz = 5
 
static const int fPar_SigmaBeam = 6
 
static const int fPar_SigmaZ = 3
 
static const int fPar_X0 = 0
 
static const int fPar_Y0 = 1
 
static const int fPar_Z0 = 2
 

Detailed Description

class: BSpdfsFcn.h package: RecoVertex/BeamSpotProducer

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

version

Id:
BSpdfsFcn.h,v 1.1 2006/12/15 20:00:37 yumiceva Exp

Definition at line 23 of file BSpdfsFcn.h.

Member Function Documentation

double BSpdfsFcn::operator() ( const std::vector< double > &  params) const
virtual

Definition at line 88 of file BSpdfsFcn.cc.

References f, fBSvector, fusepdfs, funct::log(), PDFGauss_d(), PDFGauss_d_resolution(), and PDFGauss_z().

88  {
89 
90  double f = 0.0;
91 
92  //std::cout << "fusepdfs=" << fusepdfs << " params.size="<<params.size() << std::endl;
93 
94  std::vector<BSTrkParameters>::const_iterator iparam = fBSvector.begin();
95 
96  double pdf = 0;
97 
98  for( iparam = fBSvector.begin(); iparam != fBSvector.end(); ++iparam) {
99 
100 
101  if (fusepdfs == "PDFGauss_z") {
102  pdf = PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
103  }
104  else if (fusepdfs == "PDFGauss_d") {
105  pdf = PDFGauss_d( iparam->z0(), iparam->d0(),
106  iparam->sigd0(), iparam->phi0(),params);
107  }
108  else if (fusepdfs == "PDFGauss_d_resolution") {
109  pdf = PDFGauss_d_resolution( iparam->z0(), iparam->d0(),
110  iparam->phi0(), iparam->pt(),params);
111  }
112  else if (fusepdfs == "PDFGauss_d*PDFGauss_z") {
113  //std::cout << "pdf= " << pdf << std::endl;
114  pdf = PDFGauss_d( iparam->z0(), iparam->d0(),
115  iparam->sigd0(), iparam->phi0(),params)*
116  PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
117  }
118  else if (fusepdfs == "PDFGauss_d_resolution*PDFGauss_z") {
119  pdf = PDFGauss_d_resolution( iparam->z0(), iparam->d0(),
120  iparam->phi0(), iparam->pt(),params)*
121  PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
122  }
123 
124  f = log(pdf) + f;
125  }
126 
127  f= -2.0*f;
128  return f;
129 
130 }
BSpdfsFcn::fBSvector
std::vector< BSTrkParameters > fBSvector
Definition: BSpdfsFcn.h:50
BSpdfsFcn::PDFGauss_d
double PDFGauss_d(double z, double d, double sigmad, double phi, const std::vector< double > &parms) const
Definition: BSpdfsFcn.cc:21
BSpdfsFcn::PDFGauss_z
double PDFGauss_z(double z, double sigmaz, const std::vector< double > &parms) const
Definition: BSpdfsFcn.cc:67
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
BSpdfsFcn::PDFGauss_d_resolution
double PDFGauss_d_resolution(double z, double d, double phi, double pt, const std::vector< double > &parms) const
Definition: BSpdfsFcn.cc:43
BSpdfsFcn::fusepdfs
std::string fusepdfs
Definition: BSpdfsFcn.h:49
funct::log
Log< T >::type log(const T &t)
Definition: Log.h:22
double BSpdfsFcn::PDFGauss_d ( double  z,
double  d,
double  sigmad,
double  phi,
const std::vector< double > &  parms 
) const
private

class: BSpdfsFcn.cc package: RecoVertex/BeamSpotProducer

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

version

Id:
BSpdfsFcn.cc,v 1.0 2006/09/19 17:13:31 yumiceva Exp

Definition at line 21 of file BSpdfsFcn.cc.

References funct::cos(), funct::exp(), fPar_dxdz, fPar_dydz, fPar_SigmaBeam, fPar_X0, fPar_Y0, Pi, query::result, funct::sin(), and mathSSE::sqrt().

Referenced by operator()().

22  {
23 
24  //---------------------------------------------------------------------------
25  // PDF for d0 distribution. This PDF is a simple gaussian in the
26  // beam reference frame.
27  //---------------------------------------------------------------------------
28  double fsqrt2pi = sqrt(2.* TMath::Pi());
29 
30  double sig = sqrt(parms[fPar_SigmaBeam]*parms[fPar_SigmaBeam] +
31  sigmad*sigmad);
32 
33  double dprime = d - ( ( parms[fPar_X0] + z*parms[fPar_dxdz] )*sin(phi)
34  - ( parms[fPar_Y0] + z*parms[fPar_dydz] )*cos(phi) );
35 
36  double result = (exp(-(dprime*dprime)/(2.0*sig*sig)))/(sig*fsqrt2pi);
37 
38  return result;
39 
40 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
BSpdfsFcn::fPar_dydz
static const int fPar_dydz
Definition: BSpdfsFcn.h:57
BSpdfsFcn::fPar_Y0
static const int fPar_Y0
Definition: BSpdfsFcn.h:53
BSpdfsFcn::fPar_X0
static const int fPar_X0
Definition: BSpdfsFcn.h:52
BSpdfsFcn::fPar_dxdz
static const int fPar_dxdz
Definition: BSpdfsFcn.h:56
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
BSpdfsFcn::fPar_SigmaBeam
static const int fPar_SigmaBeam
Definition: BSpdfsFcn.h:58
z
Definition: DDAxes.h:10
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
query.result
tuple result
Definition: query.py:137
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
phi
Definition: DDAxes.h:10
double BSpdfsFcn::PDFGauss_d_resolution ( double  z,
double  d,
double  phi,
double  pt,
const std::vector< double > &  parms 
) const
private

Definition at line 43 of file BSpdfsFcn.cc.

References funct::cos(), funct::exp(), fPar_c0, fPar_c1, fPar_dxdz, fPar_dydz, fPar_SigmaBeam, fPar_X0, fPar_Y0, Pi, ExpressReco_HICollisions_FallBack::pt, query::result, funct::sin(), and mathSSE::sqrt().

Referenced by operator()().

43  {
44 
45  //---------------------------------------------------------------------------
46  // PDF for d0 distribution. This PDF is a simple gaussian in the
47  // beam reference frame. The IP resolution is parametrize by a linear
48  // function as a function of 1/pt.
49  //---------------------------------------------------------------------------
50  double fsqrt2pi = sqrt(2.* TMath::Pi());
51 
52  double sigmad = parms[fPar_c0] + parms[fPar_c1]/pt;
53 
54  double sig = sqrt(parms[fPar_SigmaBeam]*parms[fPar_SigmaBeam] +
55  sigmad*sigmad);
56 
57  double dprime = d - ( ( parms[fPar_X0] + z*parms[fPar_dxdz] )*sin(phi)
58  - ( parms[fPar_Y0] + z*parms[fPar_dydz] )*cos(phi) );
59 
60  double result = (exp(-(dprime*dprime)/(2.0*sig*sig)))/(sig*fsqrt2pi);
61 
62  return result;
63 
64 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
BSpdfsFcn::fPar_dydz
static const int fPar_dydz
Definition: BSpdfsFcn.h:57
BSpdfsFcn::fPar_Y0
static const int fPar_Y0
Definition: BSpdfsFcn.h:53
BSpdfsFcn::fPar_X0
static const int fPar_X0
Definition: BSpdfsFcn.h:52
BSpdfsFcn::fPar_dxdz
static const int fPar_dxdz
Definition: BSpdfsFcn.h:56
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
BSpdfsFcn::fPar_c0
static const int fPar_c0
Definition: BSpdfsFcn.h:59
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
BSpdfsFcn::fPar_SigmaBeam
static const int fPar_SigmaBeam
Definition: BSpdfsFcn.h:58
z
Definition: DDAxes.h:10
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
query.result
tuple result
Definition: query.py:137
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
BSpdfsFcn::fPar_c1
static const int fPar_c1
Definition: BSpdfsFcn.h:60
ExpressReco_HICollisions_FallBack.pt
tuple pt
Definition: ExpressReco_HICollisions_FallBack.py:1070
phi
Definition: DDAxes.h:10
double BSpdfsFcn::PDFGauss_z ( double  z,
double  sigmaz,
const std::vector< double > &  parms 
) const
private

Definition at line 67 of file BSpdfsFcn.cc.

References funct::exp(), fPar_SigmaZ, fPar_Z0, Pi, query::result, and mathSSE::sqrt().

Referenced by operator()().

67  {
68 
69  //---------------------------------------------------------------------------
70  // PDF for z-vertex distribution. This distribution
71  // is parametrized by a simple normalized gaussian distribution.
72  //---------------------------------------------------------------------------
73  double fsqrt2pi = sqrt(2.* TMath::Pi());
74 
75  double sig = sqrt(sigmaz*sigmaz+parms[fPar_SigmaZ]*parms[fPar_SigmaZ]);
76  //double sig = sqrt(sigmaz*sigmaz+parms[1]*parms[1]);
77  double result = (exp(-((z-parms[fPar_Z0])*(z-parms[fPar_Z0]))/(2.0*sig*sig)))/(sig*fsqrt2pi);
78  //double result = (exp(-((z-parms[0])*(z-parms[0]))/(2.0*sig*sig)))/(sig*fsqrt2pi);
79 
80  return result;
81 
82 }
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
BSpdfsFcn::fPar_SigmaZ
static const int fPar_SigmaZ
Definition: BSpdfsFcn.h:55
z
Definition: DDAxes.h:10
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
query.result
tuple result
Definition: query.py:137
BSpdfsFcn::fPar_Z0
static const int fPar_Z0
Definition: BSpdfsFcn.h:54
void BSpdfsFcn::SetData ( std::vector< BSTrkParameters >  a_BSvector)
inline

Definition at line 27 of file BSpdfsFcn.h.

References fBSvector.

27  {
28 
29  fBSvector = a_BSvector;
30 
31  };
BSpdfsFcn::fBSvector
std::vector< BSTrkParameters > fBSvector
Definition: BSpdfsFcn.h:50
void BSpdfsFcn::SetPDFs ( std::string  usepdfs)
inline

Definition at line 33 of file BSpdfsFcn.h.

References fusepdfs.

33  {
34  fusepdfs = usepdfs;
35  }
BSpdfsFcn::fusepdfs
std::string fusepdfs
Definition: BSpdfsFcn.h:49
virtual double BSpdfsFcn::Up ( ) const
inlinevirtual

Definition at line 38 of file BSpdfsFcn.h.

38 {return 1.;}

Member Data Documentation

std::vector< BSTrkParameters > BSpdfsFcn::fBSvector
private

Definition at line 50 of file BSpdfsFcn.h.

Referenced by operator()(), and SetData().

const int BSpdfsFcn::fPar_c0 = 7
staticprivate

Definition at line 59 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_c1 = 8
staticprivate

Definition at line 60 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_dxdz = 4
staticprivate

Definition at line 56 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_dydz = 5
staticprivate

Definition at line 57 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_SigmaBeam = 6
staticprivate

Definition at line 58 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_SigmaZ = 3
staticprivate

Definition at line 55 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

const int BSpdfsFcn::fPar_X0 = 0
staticprivate

Definition at line 52 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_Y0 = 1
staticprivate

Definition at line 53 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_Z0 = 2
staticprivate

Definition at line 54 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

std::string BSpdfsFcn::fusepdfs
private

Definition at line 49 of file BSpdfsFcn.h.

Referenced by operator()(), and SetPDFs().

Generated for CMSSW Reference Manual by doxygen 1.8.5