Selects PFCandidates basing on their compatibility with vertex. More...

#include "CommonTools/ParticleFlow/interface/IPCutPFCandidateSelectorDefinition.h"

Inheritance diagram for pf2pat::IPCutPFCandidateSelectorDefinition:

Public Member Functions
	IPCutPFCandidateSelectorDefinition (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)

void	select (const HandleToCollection &hc, const edm::Event &e, const edm::EventSetup &s)

Public Member Functions inherited from pf2pat::PFCandidateSelectorDefinition
const_iterator	begin () const

const_iterator	end () const

	PFCandidateSelectorDefinition ()

const container &	selected () const

size_t	size () const

Private Attributes
double	d0Cut_

double	d0SigCut_

double	dtCut_

double	dtSigCut_

double	dzCut_

double	dzSigCut_

edm::EDGetTokenT < reco::VertexCollection >	verticesToken_

Additional Inherited Members
Public Types inherited from pf2pat::PFCandidateSelectorDefinition
typedef reco::PFCandidateCollection	collection

typedef boost::transform_iterator < Pointer, container::const_iterator >	const_iterator

typedef std::vector < reco::PFCandidate >	container

typedef edm::Handle< collection >	HandleToCollection

Protected Attributes inherited from pf2pat::PFCandidateSelectorDefinition
container	selected_

Detailed Description

Selects PFCandidates basing on their compatibility with vertex.

Author: Giovanni Petrucciani

Version

Id:: IPCutPFCandidateSelectorDefinition.h,v 1.2 2011/04/06 12:12:38 rwolf Exp

Definition at line 24 of file IPCutPFCandidateSelectorDefinition.h.

Constructor & Destructor Documentation

pf2pat::IPCutPFCandidateSelectorDefinition::IPCutPFCandidateSelectorDefinition	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 25 of file IPCutPFCandidateSelectorDefinition.h.

         : verticesToken_(iC.consumes<reco::VertexCollection>(cfg.getParameter<edm::InputTag>("vertices"))),
           d0Cut_(cfg.getParameter<double>("d0Cut")),
           dzCut_(cfg.getParameter<double>("dzCut")),
           dtCut_(cfg.getParameter<double>("dtCut")),
           d0SigCut_(cfg.getParameter<double>("d0SigCut")),
           dzSigCut_(cfg.getParameter<double>("dzSigCut")),
           dtSigCut_(cfg.getParameter<double>("dtSigCut")) {}

Member Function Documentation

void pf2pat::IPCutPFCandidateSelectorDefinition::select	(	const HandleToCollection &	hc,
		const edm::Event &	e,
		const edm::EventSetup &	s
	)

inline

Definition at line 34 of file IPCutPFCandidateSelectorDefinition.h.

References d0, d0Cut_, d0SigCut_, dt, dtCut_, dtSigCut_, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), PVValHelper::dz, reco::TrackBase::dz(), dzCut_, reco::TrackBase::dzError(), dzSigCut_, edm::Event::getByToken(), submitPVResolutionJobs::key, compareTotals::passing, reco::Vertex::position(), pf2pat::PFCandidateSelectorDefinition::selected_, mathSSE::sqrt(), reco::Vertex::t(), reco::Vertex::tError(), beam_dqm_sourceclient-live_cfg::vertices, verticesToken_, reco::Vertex::xError(), reco::Vertex::yError(), and reco::Vertex::zError().

                                                                                          {
       selected_.clear();
 
       edm::Handle<reco::VertexCollection> vertices;
       e.getByToken(verticesToken_, vertices);
       if (vertices->empty())
         return;
       const reco::Vertex &vtx = (*vertices)[0];
       double vt = vtx.t();
       double vte = vtx.tError();
 
       unsigned key = 0;
       for (collection::const_iterator pfc = hc->begin(); pfc != hc->end(); ++pfc, ++key) {
         bool passing = true;
         const reco::Track *tk = nullptr;
         if (pfc->gsfTrackRef().isNonnull())
           tk = pfc->gsfTrackRef().get();
         else if (pfc->trackRef().isNonnull())
           tk = pfc->trackRef().get();
 
         if (tk != nullptr) {
           double d0 = fabs(tk->dxy(vtx.position()));
           double dz = fabs(tk->dz(vtx.position()));
           double d0e = hypot(tk->dxyError(), hypot(vtx.xError(), vtx.yError()));
           double dze = hypot(tk->dzError(), vtx.zError());
           if (d0Cut_ > 0 && d0 > d0Cut_)
             passing = false;
           if (dzCut_ > 0 && dz > dzCut_)
             passing = false;
           if (d0SigCut_ > 0 && d0e > 0 && d0 / d0e > d0SigCut_)
             passing = false;
           if (dzSigCut_ > 0 && dze > 0 && dz / dze > dzSigCut_)
             passing = false;
         }
         double pfct = pfc->time();
         double pfcte = pfc->timeError();
         double dt = fabs(pfct - vt);
         double dte = std::sqrt(pfcte * pfcte + vte * vte);
         if (dtCut_ > 0 && pfcte > 0 && vte > 0 && dt > dtCut_)
           passing = false;
         if (dtSigCut_ > 0 && pfcte > 0 && vte > 0 && dt / dte > dtSigCut_)
           passing = false;
 
         if (passing) {
           selected_.push_back(reco::PFCandidate(*pfc));
           reco::PFCandidatePtr ptrToMother(hc, key);
 
           if (pfc->numberOfSourceCandidatePtrs() > 0) {
             selected_.back().setSourceCandidatePtr(edm::Ptr<reco::PFCandidate>(pfc->sourceCandidatePtr(0)));
           } else {
             selected_.back().setSourceCandidatePtr(ptrToMother);
           }
         }
       }
     }