svhelper Namespace Reference

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

Function Documentation

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

inline

Definition at line 51 of file TrackVertexArbitration.h.

References reco::Vertex::covariance().

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

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

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

inline

Definition at line 52 of file TrackVertexArbitration.h.

References reco::VertexCompositePtrCandidate::vertexCovariance().

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

void svhelper::updateVertexTime ( TransientVertex &  vtx )

inline

Definition at line 11 of file SVTimeHelpers.h.

References edm::isFinite(), GlobalErrorBase< T, ErrorWeightType >::matrix4D(), TransientVertex::originalTracks(), TransientVertex::position(), TransientVertex::positionError(), ntuplemaker::time, TransientVertex::totalChiSquared(), 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());
    }
  }