#include <TrackSelector.h>

Public Member Functions
bool	operator() (const reco::Track &track, const reco::TrackIPTagInfo::TrackIPData &ipData, const reco::Jet &jet, const GlobalPoint &pv) const
	TrackSelector (const edm::ParameterSet &params)
	~TrackSelector ()
Private Attributes
double	maxDecayLen
double	maxDistToAxis
double	maxJetDeltaR
double	maxNormChi2
unsigned int	minPixelHits
double	minPt
unsigned int	minTotalHits
reco::TrackBase::TrackQuality	quality
bool	selectQuality
double	sip2dSigMax
double	sip2dSigMin
double	sip2dValMax
double	sip2dValMin
double	sip3dSigMax
double	sip3dSigMin
double	sip3dValMax
double	sip3dValMin

Detailed Description

Definition at line 13 of file TrackSelector.h.

Constructor & Destructor Documentation

TrackSelector::TrackSelector ( const edm::ParameterSet & params )

Definition at line 16 of file TrackSelector.cc.

References edm::ParameterSet::getParameter(), quality, reco::TrackBase::qualityByName(), selectQuality, and reco::TrackBase::undefQuality.

                                                          :
        minPixelHits(params.getParameter<unsigned int>("pixelHitsMin")),
        minTotalHits(params.getParameter<unsigned int>("totalHitsMin")),
        minPt(params.getParameter<double>("ptMin")),
        maxNormChi2(params.getParameter<double>("normChi2Max")),
        maxJetDeltaR(params.getParameter<double>("jetDeltaRMax")),
        maxDistToAxis(params.getParameter<double>("maxDistToAxis")),
        maxDecayLen(params.getParameter<double>("maxDecayLen")),
        sip2dValMin(params.getParameter<double>("sip2dValMin")),
        sip2dValMax(params.getParameter<double>("sip2dValMax")),
        sip2dSigMin(params.getParameter<double>("sip2dSigMin")),
        sip2dSigMax(params.getParameter<double>("sip2dSigMax")),
        sip3dValMin(params.getParameter<double>("sip3dValMin")),
        sip3dValMax(params.getParameter<double>("sip3dValMax")),
        sip3dSigMin(params.getParameter<double>("sip3dSigMin")),
        sip3dSigMax(params.getParameter<double>("sip3dSigMax"))
{
        std::string qualityClass =
                        params.getParameter<std::string>("qualityClass");
        if (qualityClass == "any" || qualityClass == "Any" ||
            qualityClass == "ANY" || qualityClass == "") {
                selectQuality = false;
                quality = reco::TrackBase::undefQuality;
        } else {
                selectQuality = true;
                quality = reco::TrackBase::qualityByName(qualityClass);
        }
}

reco::TrackSelector::~TrackSelector ( ) [inline]

Definition at line 16 of file TrackSelector.h.

{}

Member Function Documentation

bool TrackSelector::operator()	(	const reco::Track &	track,
		const reco::TrackIPTagInfo::TrackIPData &	ipData,
		const reco::Jet &	jet,
		const GlobalPoint &	pv
	)		const

Definition at line 46 of file TrackSelector.cc.

References abs, reco::TrackIPTagInfo::TrackIPData::closestToJetAxis, reco::TrackIPTagInfo::TrackIPData::distanceToJetAxis, reco::TrackBase::hitPattern(), reco::TrackIPTagInfo::TrackIPData::ip2d, reco::TrackIPTagInfo::TrackIPData::ip3d, maxDecayLen, maxDistToAxis, maxJetDeltaR, maxNormChi2, minPixelHits, minPt, minTotalHits, reco::LeafCandidate::momentum(), reco::TrackBase::momentum(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidHits(), reco::HitPattern::numberOfValidPixelHits(), reco::TrackBase::pt(), quality, reco::TrackBase::quality(), selectQuality, Measurement1D::significance(), sip2dSigMax, sip2dSigMin, sip2dValMax, sip2dValMin, sip3dSigMax, sip3dSigMin, sip3dValMax, sip3dValMin, and Measurement1D::value().

{
        return (!selectQuality || track.quality(quality)) &&
               (minPixelHits <= 0 ||
                track.hitPattern().numberOfValidPixelHits() >= (int)minPixelHits) &&
               (minTotalHits <= 0 ||
                track.hitPattern().numberOfValidHits() >= (int)minTotalHits) &&
               track.pt() >= minPt &&
               track.normalizedChi2() < maxNormChi2 &&
               VectorUtil::DeltaR(jet.momentum(),
                                  track.momentum()) < maxJetDeltaR &&
               std::abs(ipData.distanceToJetAxis.value()) <= maxDistToAxis &&
               (ipData.closestToJetAxis - pv).mag() <= maxDecayLen &&
               ipData.ip2d.value()        >= sip2dValMin &&
               ipData.ip2d.value()        <= sip2dValMax &&
               ipData.ip2d.significance() >= sip2dSigMin &&
               ipData.ip2d.significance() <= sip2dSigMax &&
               ipData.ip3d.value()        >= sip3dValMin &&
               ipData.ip3d.value()        <= sip3dValMax &&
               ipData.ip3d.significance() >= sip3dSigMin &&
               ipData.ip3d.significance() <= sip3dSigMax;
}