getBestVertex.h File Reference

#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include <limits>

Go to the source code of this file.

Typedefs
using	Point = math::XYZPoint

Functions
Point	getBestVertex (reco::Track const &trk, reco::VertexCollection const &vertices, const size_t minNtracks=2)
Point	getBestVertex_withError (reco::Track const &trk, reco::VertexCollection const &vertices, Point &error, const size_t minNtracks=2)

Typedef Documentation

using Point = math::XYZPoint

Definition at line 6 of file getBestVertex.h.

Function Documentation

Point getBestVertex ( reco::Track const &  trk, reco::VertexCollection const &  vertices, const size_t  minNtracks = 2  )

inline

Definition at line 8 of file getBestVertex.h.

References funct::abs(), PVValHelper::dz, reco::TrackBase::dz(), SiStripPI::max, PDWG_EXOHSCP_cff::tracks, and bphysicsOniaDQM_cfi::vertex.

Referenced by TrackingNtuple::fillTracks(), and MultiTrackSelector::produce().

                                                        {
  Point p_dz(0, 0, -99999);
  float dzmin = std::numeric_limits<float>::max();
  for (auto const& vertex : vertices) {
    size_t tracks = vertex.tracksSize();
    if (tracks < minNtracks) {
      continue;
    }
    float dz = std::abs(trk.dz(vertex.position()));
    if (dz < dzmin) {
      p_dz = vertex.position();
      dzmin = dz;
    }
  }

  return p_dz;
}

Point getBestVertex_withError ( reco::Track const &  trk, reco::VertexCollection const &  vertices, Point &  error, const size_t  minNtracks = 2  )

inline

Definition at line 28 of file getBestVertex.h.

References funct::abs(), PVValHelper::dz, reco::TrackBase::dz(), SiStripPI::max, PDWG_EXOHSCP_cff::tracks, and bphysicsOniaDQM_cfi::vertex.

                                                                  {
  Point p_dz(0, 0, -99999);
  float dzmin = std::numeric_limits<float>::max();
  for (auto const& vertex : vertices) {
    size_t tracks = vertex.tracksSize();
    if (tracks < minNtracks) {
      continue;
    }
    float dz = std::abs(trk.dz(vertex.position()));
    if (dz < dzmin) {
      p_dz = vertex.position();
      error.SetXYZ(vertex.xError(), vertex.yError(), vertex.zError());
      dzmin = dz;
    }
  }

  return p_dz;
}