#include <RecoTauQualityCuts.h>

Public Types
typedef boost::function< bool(const PFCandidate &)>	CandQCutFunc

typedef std::vector< CandQCutFunc >	CandQCutFuncCollection

typedef std::map < PFCandidate::ParticleType, CandQCutFuncCollection >	CandQCutFuncMap

typedef boost::function< bool(const TrackBaseRef &)>	TrackQCutFunc

typedef std::vector < TrackQCutFunc >	TrackQCutFuncCollection

Public Member Functions
const CandQCutFunc &	candPredicate () const

bool	filterCand (const reco::PFCandidate &cand) const
	Filter a single PFCandidate. More...

template<typename PFCandRefType >
bool	filterCandRef (const PFCandRefType &cand) const
	Filter a PFCandidate held by a smart pointer or Ref. More...

template<typename Coll >
Coll	filterCandRefs (const Coll &refcoll, bool invert=false) const
	Filter a ref vector of PFCandidates. More...

bool	filterTrack (const reco::TrackBaseRef &track) const
	Filter a single Track. More...

bool	filterTrack (const reco::TrackRef &track) const

template<typename Coll >
Coll	filterTracks (const Coll &coll, bool invert=false) const
	Filter a collection of Tracks. More...

	RecoTauQualityCuts (const edm::ParameterSet &qcuts)

void	setLeadTrack (const reco::TrackRef &leadTrack) const
	Update the leading track. More...

void	setLeadTrack (const reco::PFCandidate &leadCand) const

void	setLeadTrack (const reco::PFCandidateRef &leadCand) const

void	setPV (const reco::VertexRef &vtx) const
	Update the primary vertex. More...

const TrackQCutFunc &	trackPredicate () const
	Get the predicate used to filter. More...

Private Attributes
CandQCutFunc	candPredicate_

CandQCutFuncMap	candQCuts_

reco::TrackBaseRef	leadTrack_

reco::VertexRef	pv_

TrackQCutFunc	trackPredicate_

TrackQCutFuncCollection	trackQCuts_

Detailed Description

Definition at line 34 of file RecoTauQualityCuts.h.

Member Typedef Documentation

typedef boost::function<bool (const PFCandidate&)> reco::tau::RecoTauQualityCuts::CandQCutFunc

Definition at line 40 of file RecoTauQualityCuts.h.

typedef std::vector<CandQCutFunc> reco::tau::RecoTauQualityCuts::CandQCutFuncCollection

Definition at line 41 of file RecoTauQualityCuts.h.

typedef std::map<PFCandidate::ParticleType, CandQCutFuncCollection> reco::tau::RecoTauQualityCuts::CandQCutFuncMap

Definition at line 42 of file RecoTauQualityCuts.h.

typedef boost::function<bool (const TrackBaseRef&)> reco::tau::RecoTauQualityCuts::TrackQCutFunc

Definition at line 38 of file RecoTauQualityCuts.h.

typedef std::vector<TrackQCutFunc> reco::tau::RecoTauQualityCuts::TrackQCutFuncCollection

Definition at line 39 of file RecoTauQualityCuts.h.

Constructor & Destructor Documentation

reco::tau::RecoTauQualityCuts::RecoTauQualityCuts ( const edm::ParameterSet & qcuts )

explicit

Definition at line 220 of file RecoTauQualityCuts.cc.

References reco::tau::qcuts::AND(), candPredicate_, candQCuts_, reco::PFCandidate::e, reco::tau::qcuts::etMin_cand(), edm::hlt::Exception, edm::ParameterSet::exists(), reco::PFCandidate::gamma, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), reco::PFCandidate::h, reco::PFCandidate::h0, leadTrack_, reco::tau::qcuts::mapAndCutByType(), reco::tau::qcuts::minTrackVertexWeight(), reco::PFCandidate::mu, reco::tau::qcuts::ptMin(), pv_, AlCaHLTBitMon_QueryRunRegistry::string, trackPredicate_, trackQCuts_, reco::tau::qcuts::trkChi2(), reco::tau::qcuts::trkLongitudinalImpactParameter(), reco::tau::qcuts::trkLongitudinalImpactParameterWrtTrack(), reco::tau::qcuts::trkPixelHits(), reco::tau::qcuts::trkTrackerHits(), and reco::tau::qcuts::trkTransverseImpactParameter().

 {
   // Setup all of our predicates
   CandQCutFuncCollection chargedHadronCuts;
   CandQCutFuncCollection gammaCuts;
   CandQCutFuncCollection neutralHadronCuts;
 
   // Make sure there are no extra passed options
   std::set<std::string> passedOptionSet;
   std::vector<std::string> passedOptions = qcuts.getParameterNames();
 
   BOOST_FOREACH(const std::string& option, passedOptions) {
     passedOptionSet.insert(option);
   }
 
   // Build all the QCuts for tracks
   if ( qcuts.exists("minTrackPt") ) {
     trackQCuts_.push_back(boost::bind(qcuts::ptMin, _1, qcuts.getParameter<double>("minTrackPt")));
     passedOptionSet.erase("minTrackPt");
   }
 
   if ( qcuts.exists("maxTrackChi2") ) {
     trackQCuts_.push_back(boost::bind(qcuts::trkChi2, _1, qcuts.getParameter<double>("maxTrackChi2")));
     passedOptionSet.erase("maxTrackChi2");
   }
 
   if ( qcuts.exists("minTrackPixelHits") ) {
     uint32_t minTrackPixelHits = qcuts.getParameter<uint32_t>("minTrackPixelHits");
     if ( minTrackPixelHits >= 1 ) {
       trackQCuts_.push_back(boost::bind(qcuts::trkPixelHits, _1, minTrackPixelHits));
     }
     passedOptionSet.erase("minTrackPixelHits");
   }
 
   if ( qcuts.exists("minTrackHits") ) {
     uint32_t minTrackHits = qcuts.getParameter<uint32_t>("minTrackHits");
     if ( minTrackHits >= 1 ) {
       trackQCuts_.push_back(boost::bind(qcuts::trkTrackerHits, _1, minTrackHits));
     }
     passedOptionSet.erase("minTrackHits");
   }
 
   // The impact parameter functions are bound to our member PV, since they
   // need it to compute the discriminant value.
   if ( qcuts.exists("maxTransverseImpactParameter") ) {
     trackQCuts_.push_back(boost::bind(qcuts::trkTransverseImpactParameter, _1, &pv_, qcuts.getParameter<double>("maxTransverseImpactParameter")));
     passedOptionSet.erase("maxTransverseImpactParameter");
   }
 
   if ( qcuts.exists("maxDeltaZ") ) {
     trackQCuts_.push_back(boost::bind(qcuts::trkLongitudinalImpactParameter, _1, &pv_, qcuts.getParameter<double>("maxDeltaZ")));
     passedOptionSet.erase("maxDeltaZ");
   }
 
   if ( qcuts.exists("maxDeltaZToLeadTrack") ) {
     trackQCuts_.push_back(boost::bind(qcuts::trkLongitudinalImpactParameterWrtTrack, _1, &leadTrack_, &pv_, qcuts.getParameter<double>("maxDeltaZToLeadTrack")));
     passedOptionSet.erase("maxDeltaZToLeadTrack");
   }
 
   // Require tracks to contribute a minimum weight to the associated vertex.
   if ( qcuts.exists("minTrackVertexWeight") ) {
     double minTrackVertexWeight = qcuts.getParameter<double>("minTrackVertexWeight");
     if ( minTrackVertexWeight > -1. ) {
       trackQCuts_.push_back(boost::bind(qcuts::minTrackVertexWeight, _1, &pv_, minTrackVertexWeight));
     }
     passedOptionSet.erase("minTrackVertexWeight");
   }
 
   // Build the QCuts for gammas
   if ( qcuts.exists("minGammaEt") ) {
     gammaCuts.push_back(boost::bind(qcuts::etMin_cand, _1, qcuts.getParameter<double>("minGammaEt")));
     passedOptionSet.erase("minGammaEt");
   }
 
   // Build QCuts for netural hadrons
   if ( qcuts.exists("minNeutralHadronEt") ) {
     neutralHadronCuts.push_back(boost::bind(qcuts::etMin_cand, _1, qcuts.getParameter<double>("minNeutralHadronEt")));
     passedOptionSet.erase("minNeutralHadronEt");
   }
 
   // Check if there are any remaining unparsed QCuts
   if ( passedOptionSet.size() ) {
     std::string unParsedOptions;
     bool thereIsABadParameter = false;
     BOOST_FOREACH( const std::string& option, passedOptionSet ) {
       // Workaround for HLT - TODO FIXME
       if ( option == "useTracksInsteadOfPFHadrons" ) {
         // Crash if true - no one should have this option enabled.
         if ( qcuts.getParameter<bool>("useTracksInsteadOfPFHadrons") ) {
           throw cms::Exception("DontUseTracksInQcuts")
             << "The obsolete exception useTracksInsteadOfPFHadrons "
             << "is set to true in the quality cut config." << std::endl;
         }
         continue;
       }
       
       // If we get to this point, there is a real unknown parameter
       thereIsABadParameter = true;
       
       unParsedOptions += option;
       unParsedOptions += "\n";
     }
     if ( thereIsABadParameter ) {
       throw cms::Exception("BadQualityCutConfig")
         << " The PSet passed to the RecoTauQualityCuts class had"
         << " the following unrecognized options: " << std::endl
         << unParsedOptions;
     }
   }
 
   // Make sure there are at least some quality cuts
   size_t nCuts = chargedHadronCuts.size() + gammaCuts.size() + neutralHadronCuts.size() + trackQCuts_.size();
   if ( !nCuts ) {
     throw cms::Exception("BadQualityCutConfig")
       << " No options were passed to the quality cut class!" << std::endl;
   }
 
   // Build final level predicate that works on Tracks
   trackPredicate_ = boost::bind(qcuts::AND, _1, boost::cref(trackQCuts_));
 
   // Map our QCut collections to the particle Ids they are associated to.
   candQCuts_[PFCandidate::h]     = chargedHadronCuts;
   candQCuts_[PFCandidate::gamma] = gammaCuts;
   candQCuts_[PFCandidate::h0]    = neutralHadronCuts;
   // We use the same qcuts for muons/electrons and charged hadrons.
   candQCuts_[PFCandidate::e]     = chargedHadronCuts;
   candQCuts_[PFCandidate::mu]    = chargedHadronCuts;
 
   // Build a final level predicate that works on any PFCand
   candPredicate_ = boost::bind(qcuts::mapAndCutByType, _1, boost::cref(candQCuts_));
 }

Member Function Documentation

const CandQCutFunc& reco::tau::RecoTauQualityCuts::candPredicate ( ) const

inline

Definition at line 59 of file RecoTauQualityCuts.h.

References candPredicate_.

59 { return candPredicate_; }

reco::tau::RecoTauQualityCuts::candPredicate_

CandQCutFunc candPredicate_

Definition: RecoTauQualityCuts.h:105

bool reco::tau::RecoTauQualityCuts::filterCand ( const reco::PFCandidate & cand ) const

Filter a single PFCandidate.

Definition at line 375 of file RecoTauQualityCuts.cc.

References candPredicate_, and trackPredicate_.

Referenced by filterCandRef().

 {
   auto track = getTrackRef(cand);
   if ( track.isNonnull() ) return (trackPredicate_(track) & candPredicate_(cand));
   return candPredicate_(cand);
 }

template<typename PFCandRefType >

bool reco::tau::RecoTauQualityCuts::filterCandRef ( const PFCandRefType & cand ) const

inline

Filter a PFCandidate held by a smart pointer or Ref.

Definition at line 81 of file RecoTauQualityCuts.h.

References filterCand().

Referenced by filterCandRefs().

81 { return filterCand(*cand); }

reco::tau::RecoTauQualityCuts::filterCand

bool filterCand(const reco::PFCandidate &cand) const

Filter a single PFCandidate.

Definition: RecoTauQualityCuts.cc:375

template<typename Coll >

Coll reco::tau::RecoTauQualityCuts::filterCandRefs	(	const Coll &	refcoll,
		bool	invert = `false`
	)		const

inline

Filter a ref vector of PFCandidates.

Definition at line 85 of file RecoTauQualityCuts.h.

References filterCandRef(), and convertSQLitetoXML_cfg::output.

Referenced by reco::tau::RecoTauVertexAssociator::getLeadTrack(), reco::tau::RecoTauBuilderCombinatoricPlugin::operator()(), reco::tau::RecoTauPiZeroTrivialPlugin::operator()(), reco::tau::RecoTauPiZeroCombinatoricPlugin::operator()(), reco::tau::RecoTauBuilderConePlugin::operator()(), reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::operator()(), reco::tau::RecoTauPiZeroStripPlugin::operator()(), and reco::tau::RecoTauPiZeroStripPlugin2::operator()().

   {
     Coll output;
     BOOST_FOREACH( const typename Coll::value_type cand, refcoll ) {
       if ( filterCandRef(cand)^invert ) output.push_back(cand);
     }
     return output;
   }

bool reco::tau::RecoTauQualityCuts::filterTrack ( const reco::TrackBaseRef & track ) const

Filter a single Track.

Definition at line 370 of file RecoTauQualityCuts.cc.

References trackPredicate_.

Referenced by filterTracks(), and reco::tau::PFRecoTauChargedHadronFromTrackPlugin::operator()().

 {
   return trackPredicate_(track);
 }

bool reco::tau::RecoTauQualityCuts::filterTrack ( const reco::TrackRef & track ) const

inline

Definition at line 63 of file RecoTauQualityCuts.h.

References filterTrack().

Referenced by filterTrack().

63 { return filterTrack(reco::TrackBaseRef(track)); }

edm::RefToBase< reco::Track >

reco::tau::RecoTauQualityCuts::filterTrack

bool filterTrack(const reco::TrackBaseRef &track) const

Filter a single Track.

Definition: RecoTauQualityCuts.cc:370

template<typename Coll >

Coll reco::tau::RecoTauQualityCuts::filterTracks	(	const Coll &	coll,
		bool	invert = `false`
	)		const

inline

Filter a collection of Tracks.

Definition at line 67 of file RecoTauQualityCuts.h.

References filterTrack(), and convertSQLitetoXML_cfg::output.

   {
     Coll output;
     BOOST_FOREACH( const typename Coll::value_type track, coll ) {
       if ( filterTrack(track)^invert ) output.push_back(track);
     }
     return output;
   }

void reco::tau::RecoTauQualityCuts::setLeadTrack ( const reco::TrackRef & leadTrack ) const

Update the leading track.

Definition at line 382 of file RecoTauQualityCuts.cc.

References leadTrack_.

 {
   leadTrack_ = reco::TrackBaseRef(leadTrack);
 }

void reco::tau::RecoTauQualityCuts::setLeadTrack ( const reco::PFCandidate & leadCand ) const

Definition at line 387 of file RecoTauQualityCuts.cc.

References leadTrack_.

 {
   leadTrack_ = getTrackRef(leadCand);
 }

void reco::tau::RecoTauQualityCuts::setLeadTrack ( const reco::PFCandidateRef & leadCand ) const

Update the leading track (using reference) If null, this will set the lead track ref null.

Definition at line 392 of file RecoTauQualityCuts.cc.

References edm::Ref< C, T, F >::isNonnull(), and leadTrack_.

 {
   if ( leadCand.isNonnull() ) {
     leadTrack_ = getTrackRef(*leadCand);
   } else {
     // Set null
     leadTrack_ = reco::TrackBaseRef();
   }
 }