#include <TrackFilterForPVFinding.h>

Inheritance diagram for TrackFilterForPVFinding:

Public Member Functions
bool	operator() (const reco::TransientTrack &tracks) const

std::vector< reco::TransientTrack >	select (const std::vector< reco::TransientTrack > &tracks) const

	TrackFilterForPVFinding (const edm::ParameterSet &conf)

Public Member Functions inherited from TrackFilterForPVFindingBase
	TrackFilterForPVFindingBase ()

	TrackFilterForPVFindingBase (const edm::ParameterSet &conf)

virtual	~TrackFilterForPVFindingBase ()

Private Attributes
float	maxD0Sig_

float	maxNormChi2_

float	minPt_

int	minPxLayers_

int	minSiLayers_

reco::TrackBase::TrackQuality	quality_

Detailed Description

Description: track selection for PV finding

Definition at line 14 of file TrackFilterForPVFinding.h.

Constructor & Destructor Documentation

TrackFilterForPVFinding::TrackFilterForPVFinding ( const edm::ParameterSet & conf )

Definition at line 4 of file TrackFilterForPVFinding.cc.

References edm::ParameterSet::getParameter(), maxD0Sig_, maxNormChi2_, minPt_, minPxLayers_, minSiLayers_, quality_, reco::TrackBase::qualityByName(), AlCaHLTBitMon_QueryRunRegistry::string, and reco::TrackBase::undefQuality.

 {
   maxD0Sig_    = conf.getParameter<double>("maxD0Significance");
   minPt_       = conf.getParameter<double>("minPt");
   maxNormChi2_ = conf.getParameter<double>("maxNormalizedChi2");
   minSiLayers_    = conf.getParameter<int>("minSiliconLayersWithHits");  
   minPxLayers_    = conf.getParameter<int>("minPixelLayersWithHits");  
 
   // the next few lines are taken from RecoBTag/SecondaryVertex/interface/TrackSelector.h"
   std::string qualityClass =
     conf.getParameter<std::string>("trackQuality");
   if (qualityClass == "any" || qualityClass == "Any" ||
       qualityClass == "ANY" || qualityClass == "") {
     quality_ = reco::TrackBase::undefQuality;
   } else {
     quality_ = reco::TrackBase::qualityByName(qualityClass);
   }
 
 }

Member Function Documentation

bool TrackFilterForPVFinding::operator() ( const reco::TransientTrack & tracks ) const

Definition at line 26 of file TrackFilterForPVFinding.cc.

References TrajectoryStateOnSurface::globalMomentum(), reco::TransientTrack::hitPattern(), reco::TransientTrack::impactPointState(), TrajectoryStateClosestToBeamLine::isValid(), maxD0Sig_, maxNormChi2_, minPt_, minPxLayers_, minSiLayers_, reco::TransientTrack::normalizedChi2(), reco::HitPattern::pixelLayersWithMeasurement(), reco::TrackBase::quality(), quality_, Measurement1D::significance(), reco::TransientTrack::stateAtBeamLine(), reco::TransientTrack::track(), reco::HitPattern::trackerLayersWithMeasurement(), PV3DBase< T, PVType, FrameType >::transverse(), TrajectoryStateClosestToBeamLine::transverseImpactParameter(), and reco::TrackBase::undefQuality.

 {
     if (!tk.stateAtBeamLine().isValid()) return false;
     bool IPSigCut = tk.stateAtBeamLine().transverseImpactParameter().significance()<maxD0Sig_;
     bool pTCut    = tk.impactPointState().globalMomentum().transverse() > minPt_;
     bool normChi2Cut  = tk.normalizedChi2() < maxNormChi2_;
     bool nPxLayCut = tk.hitPattern().pixelLayersWithMeasurement() >= minPxLayers_;
     bool nSiLayCut =  tk.hitPattern().trackerLayersWithMeasurement() >= minSiLayers_;
     bool trackQualityCut = (quality_==reco::TrackBase::undefQuality)|| tk.track().quality(quality_);
 
     return IPSigCut && pTCut && normChi2Cut && nPxLayCut && nSiLayCut && trackQualityCut;
 }

std::vector< reco::TransientTrack > TrackFilterForPVFinding::select ( const std::vector< reco::TransientTrack > & tracks ) const

virtual

Implements TrackFilterForPVFindingBase.

Definition at line 42 of file TrackFilterForPVFinding.cc.

Referenced by HITrackFilterForPVFinding::select().

 {
   std::vector <reco::TransientTrack> seltks;
   for (std::vector<reco::TransientTrack>::const_iterator itk = tracks.begin();
        itk != tracks.end(); itk++) {
     if ( operator()(*itk) ) seltks.push_back(*itk);  //  calls the filter function for single tracks
   }
   return seltks;
 }