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	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 22 of file IPCutPFCandidateSelectorDefinition.h.

Constructor & Destructor Documentation

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

inline

Definition at line 24 of file IPCutPFCandidateSelectorDefinition.h.

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

Member Function Documentation

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

inline

Definition at line 31 of file IPCutPFCandidateSelectorDefinition.h.

References d0Cut_, d0SigCut_, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), dzCut_, reco::TrackBase::dzError(), dzSigCut_, edm::Event::getByToken(), relval_steps::key, reco::Vertex::position(), pf2pat::PFCandidateSelectorDefinition::selected_, HLT_25ns10e33_v2_cff::vertices, verticesToken_, reco::Vertex::xError(), reco::Vertex::yError(), and reco::Vertex::zError().

Referenced by Vispa.Views.AbstractView.AbstractView::restoreSelection().

                                  {
       selected_.clear();
 
       edm::Handle<reco::VertexCollection> vertices;
       e.getByToken(verticesToken_, vertices);
       if (vertices->empty()) return;
       const reco::Vertex &vtx = (*vertices)[0];
 
       unsigned key=0;
       for( collection::const_iterator pfc = hc->begin();
        pfc != hc->end(); ++pfc, ++key) {
 
     bool passing = true;
     const reco::Track *tk = 0;
     if (pfc->gsfTrackRef().isNonnull())   tk = pfc->gsfTrackRef().get();
     else if (pfc->trackRef().isNonnull()) tk = pfc->trackRef().get();
 
     if (tk != 0) {
       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;
     }
 
     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 );
       }
     }
       }
     }