HitPixelLayersTPSelector Class Reference

#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_

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 19 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 31 of file HitPixelLayersTPSelector.h.

typedef TrackingParticleRefVector HitPixelLayersTPSelector::container

Definition at line 28 of file HitPixelLayersTPSelector.h.

Constructor & Destructor Documentation

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

inline

Definition at line 34 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"))
      {};

Member Function Documentation

const_iterator HitPixelLayersTPSelector::begin ( void  ) const

inline

Definition at line 116 of file HitPixelLayersTPSelector.h.

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

    {
      return selected_.begin();
    }

const_iterator HitPixelLayersTPSelector::end ( void  ) const

inline

Definition at line 122 of file HitPixelLayersTPSelector.h.

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

    {
      return selected_.end();
    }

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

inline

Definition at line 105 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 96 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 51 of file HitPixelLayersTPSelector.h.

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

    {
      selected_.clear();
      //Retrieve tracker topology from geometry
      edm::ESHandle<TrackerTopology> tTopoHand;
      iSetup.get<TrackerTopologyRcd>().get(tTopoHand);
      const TrackerTopology *tTopo=tTopoHand.product();


      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);
        }

    }
    }

size_t HitPixelLayersTPSelector::size ( void  ) const

inline

Definition at line 128 of file HitPixelLayersTPSelector.h.

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

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

    {
      return selected_.size();
    }

Member Data Documentation

bool HitPixelLayersTPSelector::chargedOnly_

Definition at line 144 of file HitPixelLayersTPSelector.h.

Referenced by select().

double HitPixelLayersTPSelector::lip_

Definition at line 141 of file HitPixelLayersTPSelector.h.

Referenced by select().

double HitPixelLayersTPSelector::maxRapidity_

Definition at line 139 of file HitPixelLayersTPSelector.h.

Referenced by select().

int HitPixelLayersTPSelector::minHit_

Definition at line 142 of file HitPixelLayersTPSelector.h.

Referenced by select().

double HitPixelLayersTPSelector::minRapidity_

Definition at line 138 of file HitPixelLayersTPSelector.h.

Referenced by select().

std::vector<int> HitPixelLayersTPSelector::pdgId_

Definition at line 147 of file HitPixelLayersTPSelector.h.

Referenced by select().

bool HitPixelLayersTPSelector::primaryOnly_

Definition at line 145 of file HitPixelLayersTPSelector.h.

Referenced by select().

double HitPixelLayersTPSelector::ptMin_

Definition at line 137 of file HitPixelLayersTPSelector.h.

Referenced by select().

container HitPixelLayersTPSelector::selected_

Definition at line 135 of file HitPixelLayersTPSelector.h.

Referenced by begin(), end(), select(), and size().

bool HitPixelLayersTPSelector::signalOnly_

Definition at line 143 of file HitPixelLayersTPSelector.h.

Referenced by select().

double HitPixelLayersTPSelector::tip_

Definition at line 140 of file HitPixelLayersTPSelector.h.

Referenced by select().

bool HitPixelLayersTPSelector::tpStatusBased_

Definition at line 146 of file HitPixelLayersTPSelector.h.

Referenced by select().

bool HitPixelLayersTPSelector::tripletSeedOnly_

Definition at line 136 of file HitPixelLayersTPSelector.h.

Referenced by select().