#include <TrackSelector.h>

Public Member Functions
bool	operator() (const reco::Track &track, const reco::btag::TrackIPData &ipData, const reco::Jet &jet, const GlobalPoint &pv) const

bool	operator() (const reco::CandidatePtr &track, const reco::btag::TrackIPData &ipData, const reco::Jet &jet, const GlobalPoint &pv) const

bool	operator() (const reco::TrackRef &track, const reco::btag::TrackIPData &ipData, const reco::Jet &jet, const GlobalPoint &pv) const

	TrackSelector (const edm::ParameterSet &params)

	~TrackSelector ()

Private Member Functions
bool	trackSelection (const reco::Track &track, const reco::btag::TrackIPData &ipData, const reco::Jet &jet, const GlobalPoint &pv) const

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

bool	useVariableJTA_

reco::btag::variableJTAParameters	varJTApars

Detailed Description

Definition at line 14 of file TrackSelector.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file TrackSelector.cc.

References edm::ParameterSet::getParameter(), quality, reco::TrackBase::qualityByName(), HLT_2018_cff::qualityClass, selectQuality, AlCaHLTBitMon_QueryRunRegistry::string, reco::TrackBase::undefQuality, useVariableJTA_, and varJTApars.

     : 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")),
       useVariableJTA_(params.existsAs<bool>("useVariableJTA") ? params.getParameter<bool>("useVariableJTA") : false) {
   std::string qualityClass = params.getParameter<std::string>("qualityClass");
   if (qualityClass == "any" || qualityClass == "Any" || qualityClass == "ANY" || qualityClass.empty()) {
     selectQuality = false;
     quality = reco::TrackBase::undefQuality;
   } else {
     selectQuality = true;
     quality = reco::TrackBase::qualityByName(qualityClass);
   }
   if (useVariableJTA_) {
     varJTApars = {params.getParameter<double>("a_dR"),
                   params.getParameter<double>("b_dR"),
                   params.getParameter<double>("a_pT"),
                   params.getParameter<double>("b_pT"),
                   params.getParameter<double>("min_pT"),
                   params.getParameter<double>("max_pT"),
                   params.getParameter<double>("min_pT_dRcut"),
                   params.getParameter<double>("max_pT_dRcut"),
                   params.getParameter<double>("max_pT_trackPTcut")};
   }
 }

reco::TrackSelector::~TrackSelector ( )

inline

Definition at line 17 of file TrackSelector.h.

References metsig::jet, operator()(), MetAnalyzer::pv(), and HLT_2018_cff::track.

17 {}

Member Function Documentation

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

Definition at line 49 of file TrackSelector.cc.

References reco::deltaR(), reco::deltaR2(), maxJetDeltaR, minPt, reco::LeafCandidate::momentum(), reco::TrackBase::momentum(), reco::IPTagInfo< Container, Base >::passVariableJTA(), reco::LeafCandidate::pt(), reco::TrackBase::pt(), trackSelection(), useVariableJTA_, and varJTApars.

Referenced by ~TrackSelector().

                                                             {
   bool jtaPassed = false;
   if (useVariableJTA_) {
     jtaPassed = TrackIPTagInfo::passVariableJTA(
         varJTApars, jet.pt(), track.pt(), reco::deltaR(jet.momentum(), track.momentum()));
   } else
     jtaPassed = true;
 
   return track.pt() >= minPt && reco::deltaR2(jet.momentum(), track.momentum()) < maxJetDeltaR * maxJetDeltaR &&
          jtaPassed && trackSelection(track, ipData, jet, pv);
 }

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

Definition at line 64 of file TrackSelector.cc.

References reco::deltaR(), reco::deltaR2(), metsig::jet, maxJetDeltaR, minPt, reco::Candidate::momentum(), reco::LeafCandidate::momentum(), reco::IPTagInfo< Container, Base >::passVariableJTA(), reco::Candidate::pt(), reco::LeafCandidate::pt(), MetAnalyzer::pv(), reco::btag::toTrack(), trackSelection(), useVariableJTA_, and varJTApars.

                                                             {
   bool jtaPassed = false;
   if (useVariableJTA_) {
     jtaPassed = TrackIPTagInfo::passVariableJTA(
         varJTApars, jet.pt(), track->pt(), reco::deltaR(jet.momentum(), track->momentum()));
   } else
     jtaPassed = true;
 
   return track->pt() >= minPt && reco::deltaR2(jet.momentum(), track->momentum()) < maxJetDeltaR * maxJetDeltaR &&
          jtaPassed && trackSelection(*reco::btag::toTrack(track), ipData, jet, pv);
 }

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

inline

Definition at line 29 of file TrackSelector.h.

References metsig::jet, MetAnalyzer::pv(), HLT_2018_cff::track, and trackSelection().

                                                         {
       return (*this)(*track, ipData, jet, pv);
     }

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

private

Definition at line 79 of file TrackSelector.cc.

References funct::abs(), reco::btag::TrackIPData::closestToJetAxis, reco::btag::TrackIPData::distanceToJetAxis, reco::TrackBase::hitPattern(), createfilelist::int, reco::btag::TrackIPData::ip2d, reco::btag::TrackIPData::ip3d, maxDecayLen, maxDistToAxis, maxNormChi2, minPixelHits, minTotalHits, reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfValidHits(), reco::HitPattern::numberOfValidPixelHits(), quality, reco::TrackBase::quality(), selectQuality, Measurement1D::significance(), sip2dSigMax, sip2dSigMin, sip2dValMax, sip2dValMin, sip3dSigMax, sip3dSigMin, sip3dValMax, sip3dValMin, and Measurement1D::value().

Referenced by operator()().

                                                                 {
   return (!selectQuality || track.quality(quality)) &&
          (minPixelHits <= 0 || track.hitPattern().numberOfValidPixelHits() >= (int)minPixelHits) &&
          (minTotalHits <= 0 || track.hitPattern().numberOfValidHits() >= (int)minTotalHits) &&
          track.normalizedChi2() < maxNormChi2 && 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;
 }