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)

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::InputTag	vertices_

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

Constructor & Destructor Documentation

pf2pat::IPCutPFCandidateSelectorDefinition::IPCutPFCandidateSelectorDefinition ( const edm::ParameterSet & cfg )

inline

Definition at line 23 of file IPCutPFCandidateSelectorDefinition.h.

                                                                        :
       vertices_( 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 30 of file IPCutPFCandidateSelectorDefinition.h.

References d0Cut_, d0SigCut_, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), dzCut_, reco::TrackBase::dzError(), dzSigCut_, edm::Event::getByLabel(), combine::key, reco::Vertex::position(), pf2pat::PFCandidateSelectorDefinition::selected_, vertices_, reco::Vertex::xError(), reco::Vertex::yError(), and reco::Vertex::zError().

Referenced by Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Views.LineDecayView.LineDecayContainer::mousePressEvent(), Vispa.Gui.PortConnection.PointToPointConnection::mousePressEvent(), Vispa.Gui.VispaWidget.VispaWidget::mousePressEvent(), and Vispa.Views.AbstractView.AbstractView::restoreSelection().

                                  {
       selected_.clear();
     
       edm::Handle<reco::VertexCollection> vertices;
       e.getByLabel(vertices_, 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 );
       selected_.back().setSourceCandidatePtr( ptrToMother );
     }
       }
     }