BSpdfsFcn Class Reference

#include <BSpdfsFcn.h>

Inheritance diagram for BSpdfsFcn:

Public Member Functions
virtual double	operator() (const std::vector< double > &) const
void	SetData (const 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)

---

Definition at line 22 of file BSpdfsFcn.h.

Member Function Documentation

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

virtual

Definition at line 87 of file BSpdfsFcn.cc.

References f, fBSvector, fusepdfs, dqm-mbProfile::log, 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;

}

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(), create_public_lumi_plots::exp, fPar_dxdz, fPar_dydz, fPar_SigmaBeam, fPar_X0, fPar_Y0, Pi, mps_fire::result, funct::sin(), and mathSSE::sqrt().

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;

}

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

private

Definition at line 42 of file BSpdfsFcn.cc.

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

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;

}

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

private

Definition at line 66 of file BSpdfsFcn.cc.

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

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;

}

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

inline

Definition at line 26 of file BSpdfsFcn.h.

References fBSvector.

                                                                 {
        
        fBSvector = a_BSvector;
        
    };

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

inline

Definition at line 32 of file BSpdfsFcn.h.

References fusepdfs.

                                    {
        fusepdfs = usepdfs;
    }

virtual double BSpdfsFcn::Up ( ) const

inlinevirtual

Definition at line 37 of file BSpdfsFcn.h.

{return 1.;}

Member Data Documentation

std::vector< BSTrkParameters > BSpdfsFcn::fBSvector

private

Definition at line 49 of file BSpdfsFcn.h.

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

const int BSpdfsFcn::fPar_c0 = 7

staticprivate

Definition at line 58 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_c1 = 8

staticprivate

Definition at line 59 of file BSpdfsFcn.h.

Referenced by PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_dxdz = 4

staticprivate

Definition at line 55 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_dydz = 5

staticprivate

Definition at line 56 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_SigmaBeam = 6

staticprivate

Definition at line 57 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_SigmaZ = 3

staticprivate

Definition at line 54 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

const int BSpdfsFcn::fPar_X0 = 0

staticprivate

Definition at line 51 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_Y0 = 1

staticprivate

Definition at line 52 of file BSpdfsFcn.h.

Referenced by PDFGauss_d(), and PDFGauss_d_resolution().

const int BSpdfsFcn::fPar_Z0 = 2

staticprivate

Definition at line 53 of file BSpdfsFcn.h.

Referenced by PDFGauss_z().

std::string BSpdfsFcn::fusepdfs

private

Definition at line 48 of file BSpdfsFcn.h.

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