reco::tau::qcuts Namespace Reference

Functions
bool	AND (const TrackBaseRef &track, const RecoTauQualityCuts::TrackQCutFuncCollection &cuts)
bool	AND_cand (const Candidate &cand, const RecoTauQualityCuts::CandQCutFuncCollection &cuts)
bool	etMin_cand (const Candidate &cand, double cut)
bool	mapAndCutByType (const Candidate &cand, const RecoTauQualityCuts::CandQCutFuncMap &funcMap)
bool	minPackedCandVertexWeight (const pat::PackedCandidate &pCand, const reco::VertexRef *pv, double cut)
bool	minTrackVertexWeight (const TrackBaseRef &track, const reco::VertexRef *pv, double cut)
bool	minTrackVertexWeight (const TrackRef &track, const reco::VertexRef *pv, double cut)
bool	minTrackVertexWeight_cand (const Candidate &cand, const reco::VertexRef *pv, double cut)
bool	ptMin (const TrackBaseRef &track, double cut)
bool	ptMin_cand (const Candidate &cand, double cut)
bool	trkChi2 (const Track *track, double cut)
bool	trkChi2_cand (const Candidate &cand, double cut)
bool	trkLongitudinalImpactParameter (const TrackBase *track, const reco::VertexRef *pv, double cut)
bool	trkLongitudinalImpactParameter_cand (const Candidate &cand, const reco::VertexRef *pv, double cut)
bool	trkLongitudinalImpactParameterWrtTrack (const Track *track, const Track *leadTrack, const reco::VertexRef *pv, double cut)
	DZ cut, with respect to the current lead rack. More...
bool	trkLongitudinalImpactParameterWrtTrack_cand (const Candidate &cand, const reco::Track *leadTrack, const reco::VertexRef *pv, double cut)
bool	trkPixelHits (const Track *track, int cut)
bool	trkPixelHits_cand (const Candidate &cand, int cut)
bool	trkTrackerHits (const Track *track, int cut)
bool	trkTrackerHits_cand (const Candidate &cand, int cut)
bool	trkTransverseImpactParameter (const Track *track, const reco::VertexRef *pv, double cut)
bool	trkTransverseImpactParameter_cand (const Candidate &cand, const reco::VertexRef *pv, double cut)

Function Documentation

bool reco::tau::qcuts::AND	(	const TrackBaseRef &	track,
		const RecoTauQualityCuts::TrackQCutFuncCollection &	cuts
	)

Definition at line 265 of file RecoTauQualityCuts.cc.

References patCandidates_cff::func.

Referenced by BinaryOP::evaluate(), egHLT::EgObjTrigCut< T >::pass(), and egHLT::EgEvtTrigCut< T >::pass().

{
  for(auto const& func : cuts ) {
    if ( !func(track) ) return false;
  }
  return true;
}

bool reco::tau::qcuts::AND_cand	(	const Candidate &	cand,
		const RecoTauQualityCuts::CandQCutFuncCollection &	cuts
	)

Definition at line 273 of file RecoTauQualityCuts.cc.

References patCandidates_cff::func.

Referenced by mapAndCutByType().

{
  for(auto const& func : cuts ) {
    if ( !func(cand) ) return false;
  }
  return true;
}

bool reco::tau::qcuts::etMin_cand	(	const Candidate &	cand,
		double	cut
	)

Definition at line 71 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::Candidate::et(), and LogDebug.

{
  LogDebug("TauQCuts") << "<etMin_cand>: Et = " << cand.et() << ", cut = " << cut ;
  return (cand.et() > cut);
}

bool reco::tau::qcuts::mapAndCutByType	(	const Candidate &	cand,
		const RecoTauQualityCuts::CandQCutFuncMap &	funcMap
	)

Definition at line 282 of file RecoTauQualityCuts.cc.

References funct::abs(), AND_cand(), egammaForCoreTracking_cff::cuts, and reco::Candidate::pdgId().

{
  // Find the cuts that for this particle type
  RecoTauQualityCuts::CandQCutFuncMap::const_iterator cuts = funcMap.find(std::abs(cand.pdgId()));
  // Return false if we dont' know how to deal with this particle type
  if ( cuts == funcMap.end() ) return false; 
  return AND_cand(cand, cuts->second); // Otherwise AND all the cuts
}

bool reco::tau::qcuts::minPackedCandVertexWeight	(	const pat::PackedCandidate &	pCand,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 207 of file RecoTauQualityCuts.cc.

References pat::PackedCandidate::eta(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), edm::Ref< C, T, F >::key(), LogDebug, pat::PackedCandidate::phi(), pat::PackedCandidate::pt(), pat::PackedCandidate::pvAssociationQuality(), jets_cff::quality, pat::PackedCandidate::UsedInFitLoose, pat::PackedCandidate::UsedInFitTight, and pat::PackedCandidate::vertexRef().

Referenced by reco::tau::RecoTauQualityCuts::filterChargedCand(), and minTrackVertexWeight_cand().

                                                                                                     {

  if ( pv->isNull() ) {
    edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in " <<
      "RecoTauQualityCuts is invalid. - minPackedCandVertexWeight";
    return false;
  }
  //there is some low granular information on track weight in the vertex available with packed cands
  double weight = -9.9;
  if( pCand.vertexRef().isNonnull() && pCand.vertexRef().key() == pv->key() ){
    int quality = pCand.pvAssociationQuality();
    if( quality == pat::PackedCandidate::UsedInFitTight ) weight = 0.6;//0.6 as proxy for weight above 0.5
    else if( quality == pat::PackedCandidate::UsedInFitLoose ) weight = 0.1;//0.6 as proxy for weight below 0.5
  }
  LogDebug("TauQCuts") << " packedCand: Pt = " << pCand.pt() << ", eta = " << pCand.eta() << ", phi = " << pCand.phi() ;
  LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y() << ", z = " << (*pv)->position().z() ;
  LogDebug("TauQCuts") << "--> trackWeight from packedCand = " << weight << " (cut = " << cut << ")" ;
  return (weight >= cut);
}

bool reco::tau::qcuts::minTrackVertexWeight	(	const TrackBaseRef &	track,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 181 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::TrackBase::eta(), edm::Ref< C, T, F >::isNull(), LogDebug, reco::TrackBase::phi(), and reco::TrackBase::pt().

Referenced by minTrackVertexWeight_cand().

{
  if ( pv->isNull() ) {
    edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in " <<
      "RecoTauQualityCuts is invalid. - minTrackVertexWeight";
    return false;
  }
  LogDebug("TauQCuts") << " track: Pt = " << track->pt() << ", eta = " << track->eta() << ", phi = " << track->phi() ;
  LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y() << ", z = " << (*pv)->position().z() ;
  LogDebug("TauQCuts") << "--> trackWeight = " << (*pv)->trackWeight(track) << " (cut = " << cut << ")" ;
  return ((*pv)->trackWeight(track) >= cut);
}

bool reco::tau::qcuts::minTrackVertexWeight	(	const TrackRef &	track,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 194 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, edm::Ref< C, T, F >::isNull(), and LogDebug.

{
  if ( pv->isNull() ) {
    edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in " <<
      "RecoTauQualityCuts is invalid. - minTrackVertexWeight";
    return false;
  }
  LogDebug("TauQCuts") << " track: Pt = " << track->pt() << ", eta = " << track->eta() << ", phi = " << track->phi() ;
  LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y() << ", z = " << (*pv)->position().z() ;
  LogDebug("TauQCuts") << "--> trackWeight = " << (*pv)->trackWeight(track) << " (cut = " << cut << ")" ;
  return ((*pv)->trackWeight(track) >= cut);
}

bool reco::tau::qcuts::minTrackVertexWeight_cand	(	const Candidate &	cand,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 227 of file RecoTauQualityCuts.cc.

References reco::Candidate::charge(), TkAlMuonSelectors_cfi::cut, LogDebug, minPackedCandVertexWeight(), minTrackVertexWeight(), and HiIsolationCommonParameters_cff::track.

{
  auto track = getTrackRef(cand);
  if ( track.isNonnull() ) {
    return minTrackVertexWeight(track, pv, cut);
  }
  auto gsfTrack = getGsfTrackRef(cand);
  if ( gsfTrack.isNonnull() ) {
    return minTrackVertexWeight(gsfTrack, pv, cut);
  }

  const pat::PackedCandidate* pCand = dynamic_cast<const pat::PackedCandidate*>(&cand);
  if( pCand != nullptr && cand.charge() != 0) {
    return minPackedCandVertexWeight(*pCand, pv, cut);
  }
  LogDebug("TauQCuts") << "<minTrackVertexWeight_cand>: weight = N/A, cut = " << cut ;
  return false;
}

bool reco::tau::qcuts::ptMin	(	const TrackBaseRef &	track,
		double	cut
	)

Definition at line 59 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, LogDebug, and reco::TrackBase::pt().

{
  LogDebug("TauQCuts") << "<ptMin>: Pt = " << track->pt() << ", cut = " << cut ;
  return (track->pt() > cut);
}

bool reco::tau::qcuts::ptMin_cand	(	const Candidate &	cand,
		double	cut
	)

Definition at line 65 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, LogDebug, and reco::Candidate::pt().

{
  LogDebug("TauQCuts") << "<ptMin_cand>: Pt = " << cand.pt() << ", cut = " << cut ;
  return (cand.pt() > cut);
}

bool reco::tau::qcuts::trkChi2	(	const Track *	track,
		double	cut
	)

Definition at line 246 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, LogDebug, and reco::TrackBase::normalizedChi2().

Referenced by trkChi2_cand().

{
  LogDebug("TauQCuts") << "<trkChi2>: chi^2 = " << track->normalizedChi2() << ", cut = " << cut ;
  return (track->normalizedChi2() <= cut);
}

bool reco::tau::qcuts::trkChi2_cand	(	const Candidate &	cand,
		double	cut
	)

Definition at line 252 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, getTrack(), LogDebug, HiIsolationCommonParameters_cff::track, and trkChi2().

{
  auto track = getTrack(cand);
  if ( track ) {
    LogDebug("TauQCuts") << "<trkChi2_cand>: chi^2 = " << track->normalizedChi2() << ", cut = " << cut ;
    return trkChi2(track, cut);
  } else {
    LogDebug("TauQCuts") << "<trkChi2_cand>: chi^2 = N/A, cut = " << cut ;
    return false;
  }
}

bool reco::tau::qcuts::trkLongitudinalImpactParameter	(	const TrackBase *	track,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 139 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::TrackBase::dz(), reco::TrackBase::eta(), edm::Ref< C, T, F >::isNull(), LogDebug, reco::TrackBase::phi(), and reco::TrackBase::pt().

Referenced by trkLongitudinalImpactParameter_cand().

{
  if ( pv->isNull() ) {
    edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in " <<
      "RecoTauQualityCuts is invalid. - trkLongitudinalImpactParameter";
    return false;
  }
  LogDebug("TauQCuts") << " track: Pt = " << track->pt() << ", eta = " << track->eta() << ", phi = " << track->phi() ;
  LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y() << ", z = " << (*pv)->position().z() ;
  LogDebug("TauQCuts") << "--> dz = " << std::fabs(track->dz((*pv)->position())) << " (cut = " << cut << ")" ;
  return (std::fabs(track->dz((*pv)->position())) <= cut);
}

bool reco::tau::qcuts::trkLongitudinalImpactParameter_cand	(	const Candidate &	cand,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 152 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, getTrack(), LogDebug, HiIsolationCommonParameters_cff::track, and trkLongitudinalImpactParameter().

{
  auto track = getTrack(cand);
  if ( track ) {
    return trkLongitudinalImpactParameter(track, pv, cut);
  } else {
    LogDebug("TauQCuts") << "<trkLongitudinalImpactParameter_cand>: dZ = N/A, cut = " << cut ;
    return false;
  }
}

bool reco::tau::qcuts::trkLongitudinalImpactParameterWrtTrack	(	const Track *	track,
		const Track *	leadTrack,
		const reco::VertexRef *	pv,
		double	cut
	)

DZ cut, with respect to the current lead rack.

Definition at line 164 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, and reco::TrackBase::dz().

Referenced by trkLongitudinalImpactParameterWrtTrack_cand().

{
  if (!leadTrack) {
    edm::LogError("QCutsNoValidLeadTrack") << "Lead track Ref in " <<
      "RecoTauQualityCuts is invalid. - trkLongitudinalImpactParameterWrtTrack";
    return false;
  }
  return (std::fabs(track->dz((*pv)->position()) - leadTrack->dz((*pv)->position())) <= cut);
}

bool reco::tau::qcuts::trkLongitudinalImpactParameterWrtTrack_cand	(	const Candidate &	cand,
		const reco::Track *	leadTrack,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 174 of file RecoTauQualityCuts.cc.

References getTrack(), HiIsolationCommonParameters_cff::track, and trkLongitudinalImpactParameterWrtTrack().

{
  auto track = getTrack(cand);
  if ( track ) return trkLongitudinalImpactParameterWrtTrack(track, leadTrack, pv, cut);
  else return false;
}

bool reco::tau::qcuts::trkPixelHits	(	const Track *	track,
		int	cut
	)

Definition at line 77 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::TrackBase::hitPattern(), LogDebug, and reco::HitPattern::numberOfValidPixelHits().

Referenced by trkPixelHits_cand().

{
  // For some reason, the number of hits is signed
  LogDebug("TauQCuts") << "<trkPixelHits>: #Pxl hits = " << track->hitPattern().numberOfValidPixelHits() << ", cut = " << cut ;
  return (track->hitPattern().numberOfValidPixelHits() >= cut);
}

bool reco::tau::qcuts::trkPixelHits_cand	(	const Candidate &	cand,
		int	cut
	)

Definition at line 84 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, getTrack(), LogDebug, HiIsolationCommonParameters_cff::track, and trkPixelHits().

{
  // For some reason, the number of hits is signed
  auto track = getTrack(cand);
  if ( track ) {
    LogDebug("TauQCuts") << "<trkPixelHits_cand>: #Pxl hits = " << trkPixelHits(track, cut) << ", cut = " << cut ;
    return trkPixelHits(track, cut);
  } else {
    LogDebug("TauQCuts") << "<trkPixelHits_cand>: #Pxl hits = N/A, cut = " << cut ;
    return false;
  }
}

bool reco::tau::qcuts::trkTrackerHits	(	const Track *	track,
		int	cut
	)

Definition at line 97 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::TrackBase::hitPattern(), LogDebug, and reco::HitPattern::numberOfValidHits().

Referenced by trkTrackerHits_cand().

{
  LogDebug("TauQCuts") << "<trkTrackerHits>: #Trk hits = " << track->hitPattern().numberOfValidHits() << ", cut = " << cut ;
  return (track->hitPattern().numberOfValidHits() >= cut);
}

bool reco::tau::qcuts::trkTrackerHits_cand	(	const Candidate &	cand,
		int	cut
	)

Definition at line 103 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, getTrack(), LogDebug, HiIsolationCommonParameters_cff::track, and trkTrackerHits().

{
  auto track = getTrack(cand);
  if ( track ) {
    LogDebug("TauQCuts") << "<trkTrackerHits>: #Trk hits = " << track->hitPattern().numberOfValidHits() << ", cut = " << cut ;
    return trkTrackerHits(track, cut);
  } else {
    LogDebug("TauQCuts") << "<trkTrackerHits>: #Trk hits = N/A, cut = " << cut ;
    return false;
  }
}

bool reco::tau::qcuts::trkTransverseImpactParameter	(	const Track *	track,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 115 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, reco::TrackBase::dxy(), reco::TrackBase::eta(), edm::Ref< C, T, F >::isNull(), LogDebug, reco::TrackBase::phi(), and reco::TrackBase::pt().

Referenced by trkTransverseImpactParameter_cand().

{
  if ( pv->isNull() ) {
    edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in " <<
      "RecoTauQualityCuts is invalid. - trkTransverseImpactParameter";
    return false;
  }
  LogDebug("TauQCuts") << " track: Pt = " << track->pt() << ", eta = " << track->eta() << ", phi = " << track->phi() ;
  LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y() << ", z = " << (*pv)->position().z() ;
  LogDebug("TauQCuts") << "--> dxy = " << std::fabs(track->dxy((*pv)->position())) << " (cut = " << cut << ")" ;
  return (std::fabs(track->dxy((*pv)->position())) <= cut);
}

bool reco::tau::qcuts::trkTransverseImpactParameter_cand	(	const Candidate &	cand,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 128 of file RecoTauQualityCuts.cc.

References TkAlMuonSelectors_cfi::cut, getTrack(), LogDebug, HiIsolationCommonParameters_cff::track, and trkTransverseImpactParameter().

{
  auto track = getTrack(cand);
  if ( track ) {
    return trkTransverseImpactParameter(track, pv, cut);
  } else {
    LogDebug("TauQCuts") << "<trkTransverseImpactParameter_cand>: dXY = N/A, cut = " << cut ;
    return false;
  }
}