svhelper Namespace Reference

Functions
double	cov33 (const reco::Vertex &sv)
double	cov33 (const reco::VertexCompositePtrCandidate &sv)
void	updateVertexTime (TransientVertex &vtx)

Function Documentation

cov33() [1/2]

double svhelper::cov33 ( const reco::Vertex &  sv )

inline

Definition at line 48 of file TrackVertexArbitration.h.

References pfDeepBoostedJetPreprocessParams_cfi::sv.

Referenced by TrackVertexArbitration< VTX >::trackVertexArbitrator().

{ return sv.covariance(3, 3); }

cov33() [2/2]

double svhelper::cov33 ( const reco::VertexCompositePtrCandidate &  sv )

inline

Definition at line 49 of file TrackVertexArbitration.h.

References pfDeepBoostedJetPreprocessParams_cfi::sv.

{ return sv.vertexCovariance(3, 3); }

updateVertexTime()

void svhelper::updateVertexTime ( TransientVertex &  vtx )

inline

Definition at line 11 of file SVTimeHelpers.h.

References submitPVResolutionJobs::err, edm::isFinite(), hcalRecHitTable_cff::time, hltEgammaHLTExtra_cfi::trks, L1BJetProducer_cff::vtx, and w().

Referenced by TemplatedInclusiveVertexFinder< InputContainer, VTX >::produce(), and TrackVertexArbitration< VTX >::trackVertexArbitrator().

                                                     {
    const auto& trks = vtx.originalTracks();
    double meantime = 0., expv_x2 = 0., normw = 0., timecov = 0.;
    for (const auto& trk : trks) {
      if (edm::isFinite(trk.timeExt())) {
        const double time = trk.timeExt();
        const double inverr = 1.0 / trk.dtErrorExt();
        const double w = inverr * inverr;
        meantime += time * w;
        expv_x2 += time * time * w;
        normw += w;
      }
    }
    if (normw > 0.) {
      meantime = meantime / normw;
      expv_x2 = expv_x2 / normw;
      timecov = expv_x2 - meantime * meantime;
      auto err = vtx.positionError().matrix4D();
      err(3, 3) = timecov / (double)trks.size();
      vtx = TransientVertex(vtx.position(), meantime, err, vtx.originalTracks(), vtx.totalChiSquared());
    }
  }