22 double phi,
const std::vector<double>& parms)
const {
36 double result = (
exp(-(dprime*dprime)/(2.0*sig*sig)))/(sig*fsqrt2pi);
60 double result = (
exp(-(dprime*dprime)/(2.0*sig*sig)))/(sig*fsqrt2pi);
77 double result = (
exp(-((z-parms[
fPar_Z0])*(z-parms[fPar_Z0]))/(2.0*sig*sig)))/(sig*fsqrt2pi);
94 std::vector<BSTrkParameters>::const_iterator iparam =
fBSvector.begin();
102 pdf =
PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
104 else if (
fusepdfs ==
"PDFGauss_d") {
106 iparam->sigd0(), iparam->phi0(),params);
108 else if (
fusepdfs ==
"PDFGauss_d_resolution") {
110 iparam->phi0(), iparam->pt(),params);
112 else if (
fusepdfs ==
"PDFGauss_d*PDFGauss_z") {
115 iparam->sigd0(), iparam->phi0(),params)*
116 PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
118 else if (
fusepdfs ==
"PDFGauss_d_resolution*PDFGauss_z") {
120 iparam->phi0(), iparam->pt(),params)*
121 PDFGauss_z( iparam->z0(), iparam->sigz0(),params);
std::vector< BSTrkParameters > fBSvector
static const int fPar_dydz
static const int fPar_dxdz
Sin< T >::type sin(const T &t)
double PDFGauss_d(double z, double d, double sigmad, double phi, const std::vector< double > &parms) const
Exp< T >::type exp(const T &t)
double PDFGauss_z(double z, double sigmaz, const std::vector< double > &parms) const
static const int fPar_SigmaZ
static const int fPar_SigmaBeam
virtual double operator()(const std::vector< double > &) const
Cos< T >::type cos(const T &t)
double PDFGauss_d_resolution(double z, double d, double phi, double pt, const std::vector< double > &parms) const
Log< T >::type log(const T &t)