BSpdfsFcn Class Reference

#include <BSpdfsFcn.h>

Inheritance diagram for BSpdfsFcn:

Public Member Functions
double	operator() (const std::vector< double > &) const override
void	SetData (const std::vector< BSTrkParameters > &a_BSvector)
void	SetPDFs (std::string usepdfs)
double	Up () const override

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)

---

Definition at line 22 of file BSpdfsFcn.h.

Member Function Documentation

operator()()

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

override

Definition at line 76 of file BSpdfsFcn.cc.

References f, fBSvector, fusepdfs, dqm-mbProfile::log, submitPVValidationJobs::params, PDFGauss_d(), PDFGauss_d_resolution(), and PDFGauss_z().

                                                                  {
  double f = 0.0;

  //std::cout << "fusepdfs=" << fusepdfs << " params.size="<<params.size() << std::endl;

  std::vector<BSTrkParameters>::const_iterator iparam = fBSvector.begin();

  double pdf = 0;

  for (iparam = fBSvector.begin(); iparam != fBSvector.end(); ++iparam) {
    if (fusepdfs == "PDFGauss_z") {
      pdf = PDFGauss_z(iparam->z0(), iparam->sigz0(), params);
    } else if (fusepdfs == "PDFGauss_d") {
      pdf = PDFGauss_d(iparam->z0(), iparam->d0(), iparam->sigd0(), iparam->phi0(), params);
    } else if (fusepdfs == "PDFGauss_d_resolution") {
      pdf = PDFGauss_d_resolution(iparam->z0(), iparam->d0(), iparam->phi0(), iparam->pt(), params);
    } else if (fusepdfs == "PDFGauss_d*PDFGauss_z") {
      //std::cout << "pdf= " << pdf << std::endl;
      pdf = PDFGauss_d(iparam->z0(), iparam->d0(), iparam->sigd0(), iparam->phi0(), params) *
            PDFGauss_z(iparam->z0(), iparam->sigz0(), params);
    } else if (fusepdfs == "PDFGauss_d_resolution*PDFGauss_z") {
      pdf = PDFGauss_d_resolution(iparam->z0(), iparam->d0(), iparam->phi0(), iparam->pt(), params) *
            PDFGauss_z(iparam->z0(), iparam->sigz0(), params);
    }

    f = log(pdf) + f;
  }

  f = -2.0 * f;
  return f;
}

PDFGauss_d()

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)

---

Definition at line 20 of file BSpdfsFcn.cc.

References funct::cos(), ztail::d, JetChargeProducer_cfi::exp, fPar_dxdz, fPar_dydz, fPar_SigmaBeam, fPar_X0, fPar_Y0, phi, Pi, mps_fire::result, funct::sin(), mathSSE::sqrt(), and z.

Referenced by operator()().

                                                                                                                {
  //---------------------------------------------------------------------------
  //  PDF for d0 distribution. This PDF is a simple gaussian in the
  //  beam reference frame.
  //---------------------------------------------------------------------------
  double fsqrt2pi = sqrt(2. * TMath::Pi());

  double sig = sqrt(parms[fPar_SigmaBeam] * parms[fPar_SigmaBeam] + sigmad * sigmad);

  double dprime =
      d - ((parms[fPar_X0] + z * parms[fPar_dxdz]) * sin(phi) - (parms[fPar_Y0] + z * parms[fPar_dydz]) * cos(phi));

  double result = (exp(-(dprime * dprime) / (2.0 * sig * sig))) / (sig * fsqrt2pi);

  return result;
}

PDFGauss_d_resolution()

double BSpdfsFcn::PDFGauss_d_resolution ( double  z, double  d, double  phi, double  pt, const std::vector< double > &  parms  ) const

private

Definition at line 38 of file BSpdfsFcn.cc.

References funct::cos(), ztail::d, JetChargeProducer_cfi::exp, fPar_c0, fPar_c1, fPar_dxdz, fPar_dydz, fPar_SigmaBeam, fPar_X0, fPar_Y0, phi, Pi, DiDispStaMuonMonitor_cfi::pt, mps_fire::result, funct::sin(), mathSSE::sqrt(), and z.

Referenced by operator()().

                                                                                     {
  //---------------------------------------------------------------------------
  //  PDF for d0 distribution. This PDF is a simple gaussian in the
  //  beam reference frame. The IP resolution is parametrize by a linear
  //  function as a function of 1/pt.
  //---------------------------------------------------------------------------
  double fsqrt2pi = sqrt(2. * TMath::Pi());

  double sigmad = parms[fPar_c0] + parms[fPar_c1] / pt;

  double sig = sqrt(parms[fPar_SigmaBeam] * parms[fPar_SigmaBeam] + sigmad * sigmad);

  double dprime =
      d - ((parms[fPar_X0] + z * parms[fPar_dxdz]) * sin(phi) - (parms[fPar_Y0] + z * parms[fPar_dydz]) * cos(phi));

  double result = (exp(-(dprime * dprime) / (2.0 * sig * sig))) / (sig * fsqrt2pi);

  return result;
}

PDFGauss_z()

double BSpdfsFcn::PDFGauss_z ( double  z, double  sigmaz, const std::vector< double > &  parms  ) const

private

Definition at line 60 of file BSpdfsFcn.cc.

References JetChargeProducer_cfi::exp, fPar_SigmaZ, fPar_Z0, Pi, mps_fire::result, mathSSE::sqrt(), and z.

Referenced by operator()().

                                                                                          {
  //---------------------------------------------------------------------------
  //  PDF for z-vertex distribution. This distribution
  // is parametrized by a simple normalized gaussian distribution.
  //---------------------------------------------------------------------------
  double fsqrt2pi = sqrt(2. * TMath::Pi());

  double sig = sqrt(sigmaz * sigmaz + parms[fPar_SigmaZ] * parms[fPar_SigmaZ]);
  //double sig = sqrt(sigmaz*sigmaz+parms[1]*parms[1]);
  double result = (exp(-((z - parms[fPar_Z0]) * (z - parms[fPar_Z0])) / (2.0 * sig * sig))) / (sig * fsqrt2pi);
  //double result = (exp(-((z-parms[0])*(z-parms[0]))/(2.0*sig*sig)))/(sig*fsqrt2pi);

  return result;
}

SetData()

void BSpdfsFcn::SetData ( const std::vector< BSTrkParameters > &  a_BSvector )

inline

Definition at line 25 of file BSpdfsFcn.h.

References fBSvector.

{ fBSvector = a_BSvector; };

SetPDFs()

void BSpdfsFcn::SetPDFs ( std::string  usepdfs )

inline

Definition at line 27 of file BSpdfsFcn.h.

References fusepdfs.

{ fusepdfs = usepdfs; }

Up()

double BSpdfsFcn::Up ( ) const

inlineoverride

Definition at line 30 of file BSpdfsFcn.h.

{ return 1.; }

Member Data Documentation

fBSvector

std::vector<BSTrkParameters> BSpdfsFcn::fBSvector

private

Definition at line 39 of file BSpdfsFcn.h.

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

fPar_c0

const int BSpdfsFcn::fPar_c0 = 7

staticprivate

Definition at line 48 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

fPar_c1

const int BSpdfsFcn::fPar_c1 = 8

staticprivate

Definition at line 49 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

fPar_dxdz

const int BSpdfsFcn::fPar_dxdz = 4

staticprivate

Definition at line 45 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

fPar_dydz

const int BSpdfsFcn::fPar_dydz = 5

staticprivate

Definition at line 46 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

fPar_SigmaBeam

const int BSpdfsFcn::fPar_SigmaBeam = 6

staticprivate

Definition at line 47 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

fPar_SigmaZ

const int BSpdfsFcn::fPar_SigmaZ = 3

staticprivate

Definition at line 44 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

fPar_X0

const int BSpdfsFcn::fPar_X0 = 0

staticprivate

Definition at line 41 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

fPar_Y0

const int BSpdfsFcn::fPar_Y0 = 1

staticprivate

Definition at line 42 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

fPar_Z0

const int BSpdfsFcn::fPar_Z0 = 2

staticprivate

Definition at line 43 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

fusepdfs

std::string BSpdfsFcn::fusepdfs

private

Definition at line 38 of file BSpdfsFcn.h.

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