#include <HitPixelLayersTPSelector.h>

Public Types
typedef TrackingParticleCollection	collection

typedef container::const_iterator	const_iterator

typedef TrackingParticleRefVector	container

Public Member Functions
const_iterator	begin () const

const_iterator	end () const

bool	goodHitPattern (const std::vector< bool > &hitpattern)

	HitPixelLayersTPSelector (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC)

std::vector< bool >	pixelHitPattern (const TrackingParticleRef &simTrack, const TrackerTopology *tTopo)

void	select (const edm::Handle< collection > &TPCH, const edm::Event &iEvent, const edm::EventSetup &iSetup)

size_t	size () const

Public Attributes
bool	chargedOnly_

double	lip_

double	maxRapidity_

int	minHit_

double	minRapidity_

std::vector< int >	pdgId_

bool	primaryOnly_

double	ptMin_

container	selected_

bool	signalOnly_

double	tip_

bool	tpStatusBased_

bool	tripletSeedOnly_

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	tTopoToken_

Detailed Description

Selector to select only tracking particles that leave hits in three pixel layers Additional selection done on pt, rapidity, impact parameter, min hits, pdg id, etc.

Inspired by CommonTools.RecoAlgos.TrackingParticleSelector.h

Definition at line 21 of file HitPixelLayersTPSelector.h.

Member Typedef Documentation

typedef TrackingParticleCollection HitPixelLayersTPSelector::collection

Definition at line 24 of file HitPixelLayersTPSelector.h.

typedef container::const_iterator HitPixelLayersTPSelector::const_iterator

Definition at line 30 of file HitPixelLayersTPSelector.h.

typedef TrackingParticleRefVector HitPixelLayersTPSelector::container

Definition at line 27 of file HitPixelLayersTPSelector.h.

Constructor & Destructor Documentation

HitPixelLayersTPSelector::HitPixelLayersTPSelector	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 33 of file HitPixelLayersTPSelector.h.

       : tripletSeedOnly_(iConfig.getParameter<bool>("tripletSeedOnly")),
         ptMin_(iConfig.getParameter<double>("ptMin")),
         minRapidity_(iConfig.getParameter<double>("minRapidity")),
         maxRapidity_(iConfig.getParameter<double>("maxRapidity")),
         tip_(iConfig.getParameter<double>("tip")),
         lip_(iConfig.getParameter<double>("lip")),
         minHit_(iConfig.getParameter<int>("minHit")),
         signalOnly_(iConfig.getParameter<bool>("signalOnly")),
         chargedOnly_(iConfig.getParameter<bool>("chargedOnly")),
         primaryOnly_(iConfig.getParameter<bool>("primaryOnly")),
         tpStatusBased_(iConfig.getParameter<bool>("tpStatusBased")),
         pdgId_(iConfig.getParameter<std::vector<int> >("pdgId")),
         tTopoToken_(iC.esConsumes()){};

Member Function Documentation

const_iterator HitPixelLayersTPSelector::begin ( void ) const

inline

Definition at line 110 of file HitPixelLayersTPSelector.h.

References edm::RefVector< C, T, F >::begin(), and selected_.

110 { return selected_.begin(); }

HitPixelLayersTPSelector::selected_

container selected_

Definition: HitPixelLayersTPSelector.h:120

edm::RefVector::begin

const_iterator begin() const

Initialize an iterator over the RefVector.

Definition: RefVector.h:223

const_iterator HitPixelLayersTPSelector::end ( void ) const

inline

Definition at line 113 of file HitPixelLayersTPSelector.h.

References edm::RefVector< C, T, F >::end(), and selected_.

113 { return selected_.end(); }

HitPixelLayersTPSelector::selected_

container selected_

Definition: HitPixelLayersTPSelector.h:120

edm::RefVector::end

const_iterator end() const

Termination of iteration.

Definition: RefVector.h:228

bool HitPixelLayersTPSelector::goodHitPattern ( const std::vector< bool > & hitpattern )

inline

Definition at line 101 of file HitPixelLayersTPSelector.h.

Referenced by select().

                                                          {
     if ((hitpattern[0] && hitpattern[1] && hitpattern[2]) || (hitpattern[0] && hitpattern[1] && hitpattern[3]) ||
         (hitpattern[0] && hitpattern[3] && hitpattern[4]))
       return true;
     else
       return false;
   }

std::vector<bool> HitPixelLayersTPSelector::pixelHitPattern	(	const TrackingParticleRef &	simTrack,
		const TrackerTopology *	tTopo
	)

inline

Definition at line 93 of file HitPixelLayersTPSelector.h.

Referenced by select().

                                                                                                      {
     std::vector<bool> hitpattern(5, false);  // PXB 0,1,2  PXF 0,1
     // This currently will always return false, since we can no loger use the sim hits to check for triplets.  This would need to be fixed if we want to enable this feature, but it's not being used at the moment, since tripletSeedOnly is always set to False  - Matt Nguyen, 24/7/2013
 
     return hitpattern;
   }

void HitPixelLayersTPSelector::select	(	const edm::Handle< collection > &	TPCH,
		const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

inline

Definition at line 50 of file HitPixelLayersTPSelector.h.

References chargedOnly_, edm::RefVector< C, T, F >::clear(), edm::EventSetup::getData(), goodHitPattern(), mps_fire::i, lip_, maxRapidity_, minHit_, minRapidity_, pdgId_, pixelHitPattern(), primaryOnly_, edm::Handle< T >::product(), ptMin_, edm::RefVector< C, T, F >::push_back(), selected_, signalOnly_, mathSSE::sqrt(), tip_, tpStatusBased_, tripletSeedOnly_, and tTopoToken_.

                                                                                                       {
     selected_.clear();
     //Retrieve tracker topology from geometry
     const TrackerTopology* tTopo = &iSetup.getData(tTopoToken_);
 
     const collection& tpc = *(TPCH.product());
 
     for (TrackingParticleCollection::size_type i = 0; i < tpc.size(); i++) {
       TrackingParticleRef tpr(TPCH, i);
 
       // quickly reject if it is from pile-up
       if (signalOnly_ && !(tpr->eventId().bunchCrossing() == 0 && tpr->eventId().event() == 0))
         continue;
       if (chargedOnly_ && tpr->charge() == 0)
         continue;  //select only if charge!=0
       if (tpStatusBased_ && primaryOnly_ && tpr->status() != 1)
         continue;  // TP status based sel primary
       if ((!tpStatusBased_) && primaryOnly_ && tpr->parentVertex()->nSourceTracks() != 0)
         continue;  // vertex based sel for primary
 
       // loop over specified PID values
       bool testId = false;
       unsigned int idSize = pdgId_.size();
       if (idSize == 0)
         testId = true;
       else
         for (unsigned int it = 0; it != idSize; ++it) {
           if (tpr->pdgId() == pdgId_[it])
             testId = true;
         }
 
       // selection criteria
       if (tpr->numberOfTrackerLayers() >= minHit_ && sqrt(tpr->momentum().perp2()) >= ptMin_ &&
           tpr->momentum().eta() >= minRapidity_ && tpr->momentum().eta() <= maxRapidity_ &&
           sqrt(tpr->vertex().perp2()) <= tip_ && fabs(tpr->vertex().z()) <= lip_ && testId) {
         if (tripletSeedOnly_ && !goodHitPattern(pixelHitPattern(tpr, tTopo)))
           continue;  //findable triplet seed
         selected_.push_back(tpr);
       }
     }
   }