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 9 of file getBestVertex.h.

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

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 30 of file getBestVertex.h.

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

                                                                                                                                         {

    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;

  }