/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/RecoVertex/PrimaryVertexProducer/src/VertexHigherPtSquared.cc

Go to the documentation of this file.

#include "RecoVertex/PrimaryVertexProducer/interface/VertexHigherPtSquared.h"

using namespace reco;

bool 
VertexHigherPtSquared::operator() ( const TransientVertex & v1, 
                                    const TransientVertex & v2) const
{
  return (sumPtSquared(v1) > sumPtSquared(v2));
}


bool 
VertexHigherPtSquared::operator() ( const Vertex & v1, 
                                    const Vertex & v2) const
{
  return (sumPtSquared(v1) > sumPtSquared(v2));
}




double VertexHigherPtSquared::sumPtSquared(const Vertex & v) const 
{
  double sum = 0.;
  double pT;
  for (Vertex::trackRef_iterator it = v.tracks_begin(); it != v.tracks_end(); it++) {
    pT = (**it).pt();
    double epT=(**it).ptError(); pT=pT>epT ? pT-epT : 0;

    sum += pT*pT;
  }
  return sum;
}