Inheritance diagram for PFRecoTauDiscriminationByIsolation:

Public Member Functions
void	beginEvent (const edm::Event &evt, const edm::EventSetup &evtSetup) override

double	discriminate (const PFTauRef &pfTau) override

	PFRecoTauDiscriminationByIsolation (const edm::ParameterSet &pset)

	~PFRecoTauDiscriminationByIsolation ()

Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual double	discriminate (const TauRef &tau)=0

virtual void	endEvent (edm::Event &evt)

void	produce (edm::Event &, const edm::EventSetup &)

	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)

	TauDiscriminationProducerBase ()

virtual	~TauDiscriminationProducerBase ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
bool	applyDeltaBeta_

bool	applyOccupancyCut_

bool	applyRelativeSumPtCut_

bool	applyRhoCorrection_

bool	applySumPtCut_

std::vector< reco::PFCandidatePtr >	chargedPFCandidatesInEvent_

double	customIsoCone_

double	deltaBetaCollectionCone_

double	deltaBetaFactorThisEvent_

std::auto_ptr< TFormula >	deltaBetaFormula_

bool	includeGammas_

bool	includeTracks_

std::vector< PFCandidatePtr >	isoCharged_

std::vector< PFCandidatePtr >	isoNeutral_

std::vector< PFCandidatePtr >	isoPU_

uint32_t	maximumOccupancy_

double	maximumRelativeSumPt_

double	maximumSumPt_

std::string	moduleLabel_

edm::EDGetTokenT < reco::PFCandidateCollection >	pfCand_token

edm::InputTag	pfCandSrc_

std::auto_ptr < tau::RecoTauQualityCuts >	pileupQcutsGeneralQCuts_

std::auto_ptr < tau::RecoTauQualityCuts >	pileupQcutsPUTrackSelection_

std::auto_ptr < tau::RecoTauQualityCuts >	qcuts_

edm::ParameterSet	qualityCutsPSet_

edm::EDGetTokenT< double >	rho_token

double	rhoConeSize_

double	rhoCorrectionThisEvent_

edm::InputTag	rhoProducer_

double	rhoThisEvent_

double	rhoUEOffsetCorrection_

bool	storeRawOccupancy_

bool	storeRawPUsumPt_

bool	storeRawSumPt_

int	verbosity_

edm::EDGetTokenT < reco::VertexCollection >	vertex_token

std::auto_ptr < tau::RecoTauVertexAssociator >	vertexAssociator_

edm::InputTag	vertexSrc_

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
typedef std::vector< TauType >	TauCollection

typedef edm::Ref< TauCollection >	TauRef

typedef edm::RefProd < TauCollection >	TauRefProd

Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
std::string	moduleLabel_

double	prediscriminantFailValue_

edm::EDGetTokenT< TauCollection >	Tau_token

edm::InputTag	TauProducer_

Detailed Description

Definition at line 25 of file PFRecoTauDiscriminationByIsolation.cc.

Constructor & Destructor Documentation

PFRecoTauDiscriminationByIsolation::PFRecoTauDiscriminationByIsolation ( const edm::ParameterSet & pset )

inlineexplicit

Definition at line 28 of file PFRecoTauDiscriminationByIsolation.cc.

References edm::hlt::Exception, edm::ParameterSet::exists(), reco::tau::factorizePUQCuts(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), and AlCaHLTBitMon_QueryRunRegistry::string.

     : PFTauDiscriminationProducerBase(pset),
       moduleLabel_(pset.getParameter<std::string>("@module_label")),
       qualityCutsPSet_(pset.getParameter<edm::ParameterSet>("qualityCuts")) 
   {
     includeTracks_ = pset.getParameter<bool>(
       "ApplyDiscriminationByTrackerIsolation");
     includeGammas_ = pset.getParameter<bool>(
       "ApplyDiscriminationByECALIsolation");
 
     applyOccupancyCut_ = pset.getParameter<bool>("applyOccupancyCut");
     maximumOccupancy_ = pset.getParameter<uint32_t>("maximumOccupancy");
 
     applySumPtCut_ = pset.getParameter<bool>("applySumPtCut");
     maximumSumPt_ = pset.getParameter<double>("maximumSumPtCut");
 
     applyRelativeSumPtCut_ = pset.getParameter<bool>(
       "applyRelativeSumPtCut");
     maximumRelativeSumPt_ = pset.getParameter<double>(
       "relativeSumPtCut");
 
     storeRawOccupancy_ = pset.exists("storeRawOccupancy") ?
       pset.getParameter<bool>("storeRawOccupancy") : false;
     storeRawSumPt_ = pset.exists("storeRawSumPt") ?
       pset.getParameter<bool>("storeRawSumPt") : false;
     storeRawPUsumPt_ = pset.exists("storeRawPUsumPt") ?
       pset.getParameter<bool>("storeRawPUsumPt") : false;
 
     // Sanity check on requested options.  We can't apply cuts and store the
     // raw output at the same time
     if ( applySumPtCut_ || applyOccupancyCut_ || applyRelativeSumPtCut_ ) {
       if ( storeRawSumPt_ || storeRawOccupancy_ || storeRawPUsumPt_ ) {
     throw cms::Exception("BadIsoConfig") 
       << "A 'store raw' and a 'apply cut' option have been set to true "
       << "simultaneously.  These options are mutually exclusive.";
       }
     }
     
     // Can only store one type
     int numStoreOptions = 0;
     if ( storeRawSumPt_     ) ++numStoreOptions;
     if ( storeRawOccupancy_ ) ++numStoreOptions;
     if ( storeRawPUsumPt_   ) ++numStoreOptions;
     if ( numStoreOptions > 1 ) {
       throw cms::Exception("BadIsoConfig") 
     << "Both 'store sum pt' and 'store occupancy' options are set."
     << " These options are mutually exclusive.";
     }
 
     if ( pset.exists("customOuterCone") ) {
       customIsoCone_ = pset.getParameter<double>("customOuterCone");
     } else {
       customIsoCone_ = -1;
     }
 
     // Get the quality cuts specific to the isolation region
     edm::ParameterSet isolationQCuts = qualityCutsPSet_.getParameterSet(
       "isolationQualityCuts");
 
     qcuts_.reset(new tau::RecoTauQualityCuts(isolationQCuts));
 
     vertexAssociator_.reset(
                 new tau::RecoTauVertexAssociator(qualityCutsPSet_,consumesCollector()));
 
     applyDeltaBeta_ = pset.exists("applyDeltaBetaCorrection") ?
       pset.getParameter<bool>("applyDeltaBetaCorrection") : false;
 
     if ( applyDeltaBeta_ ) {
       // Factorize the isolation QCuts into those that are used to
       // select PU and those that are not.
       std::pair<edm::ParameterSet, edm::ParameterSet> puFactorizedIsoQCuts =
     reco::tau::factorizePUQCuts(isolationQCuts);
 
       // Determine the pt threshold for the PU tracks
       // First check if the user specifies explicitly the cut.
       if ( pset.exists("deltaBetaPUTrackPtCutOverride") ) {
     puFactorizedIsoQCuts.second.addParameter<double>(
       "minTrackPt",
       pset.getParameter<double>("deltaBetaPUTrackPtCutOverride"));
       } else {
     // Secondly take it from the minGammaEt
     puFactorizedIsoQCuts.second.addParameter<double>(
           "minTrackPt",
       isolationQCuts.getParameter<double>("minGammaEt"));
       }
 
       pileupQcutsPUTrackSelection_.reset(new tau::RecoTauQualityCuts(
         puFactorizedIsoQCuts.first));
 
       pileupQcutsGeneralQCuts_.reset(new tau::RecoTauQualityCuts(
         puFactorizedIsoQCuts.second));
 
       pfCandSrc_ = pset.getParameter<edm::InputTag>("particleFlowSrc");
       pfCand_token=consumes<reco::PFCandidateCollection>(pfCandSrc_);
       vertexSrc_ = pset.getParameter<edm::InputTag>("vertexSrc");
       vertex_token=consumes<reco::VertexCollection>(vertexSrc_);
       deltaBetaCollectionCone_ = pset.getParameter<double>(
         "isoConeSizeForDeltaBeta");
       std::string deltaBetaFactorFormula =
     pset.getParameter<string>("deltaBetaFactor");
       deltaBetaFormula_.reset(
         new TFormula("DB_corr", deltaBetaFactorFormula.c_str()));
     }
 
     applyRhoCorrection_ = pset.exists("applyRhoCorrection") ?
       pset.getParameter<bool>("applyRhoCorrection") : false;
     if ( applyRhoCorrection_ ) {
       rhoProducer_ = pset.getParameter<edm::InputTag>("rhoProducer");
       rho_token=consumes<double>(rhoProducer_);
       rhoConeSize_ = pset.getParameter<double>("rhoConeSize");
       rhoUEOffsetCorrection_ =
     pset.getParameter<double>("rhoUEOffsetCorrection");
     }
 
     verbosity_ = ( pset.exists("verbosity") ) ?
       pset.getParameter<int>("verbosity") : 0;
   }

PFRecoTauDiscriminationByIsolation::~PFRecoTauDiscriminationByIsolation ( )

inline

Definition at line 146 of file PFRecoTauDiscriminationByIsolation.cc.

146 {}

Member Function Documentation

void PFRecoTauDiscriminationByIsolation::beginEvent	(	const edm::Event &	evt,
		const edm::EventSetup &	evtSetup
	)

overridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 211 of file PFRecoTauDiscriminationByIsolation.cc.

References i, and reco::tau::pfCandidates().

 {
   // NB: The use of the PV in this context is necessitated by its use in
   // applying quality cuts to the different objects in the isolation cone
   // The vertex associator contains the logic to select the appropriate vertex
   // We need to pass it the event so it can load the vertices.
   vertexAssociator_->setEvent(event);
 
   // If we are applying the delta beta correction, we need to get the PF
   // candidates from the event so we can find the PU tracks.
   if ( applyDeltaBeta_ ) {
     // Collect all the PF pile up tracks
     edm::Handle<reco::PFCandidateCollection> pfCandidates;
     event.getByToken(pfCand_token, pfCandidates);
     chargedPFCandidatesInEvent_.clear();
     chargedPFCandidatesInEvent_.reserve(pfCandidates->size());
     size_t numPFCandidates = pfCandidates->size();
     for ( size_t i = 0; i < numPFCandidates; ++i ) {
       reco::PFCandidatePtr pfCandidate(pfCandidates, i);
       if ( pfCandidate->charge() != 0 ) {
         chargedPFCandidatesInEvent_.push_back(pfCandidate);
       }
     }
     // Count all the vertices in the event, to parameterize the DB
     // correction factor
     edm::Handle<reco::VertexCollection> vertices;
     event.getByToken(vertex_token, vertices);
     size_t nVtxThisEvent = vertices->size();
     deltaBetaFactorThisEvent_ = deltaBetaFormula_->Eval(nVtxThisEvent);
   }
 
   if ( applyRhoCorrection_ ) {
     edm::Handle<double> rhoHandle_;
     event.getByToken(rho_token, rhoHandle_);
     rhoThisEvent_ = (*rhoHandle_ - rhoUEOffsetCorrection_)*
       (3.14159)*rhoConeSize_*rhoConeSize_;
   }
 }

double PFRecoTauDiscriminationByIsolation::discriminate ( const PFTauRef & pfTau )

override

Definition at line 251 of file PFRecoTauDiscriminationByIsolation.cc.

References gather_cfg::cout, alcazmumu_cfi::filter, edm::Ref< C, T, F >::isNonnull(), and moduleLabel_().

 {
   if ( verbosity_ ) {
     std::cout << "<PFRecoTauDiscriminationByIsolation::discriminate (moduleLabel = " << moduleLabel_ <<")>:" << std::endl;
     std::cout << " tau: Pt = " << pfTau->pt() << ", eta = " << pfTau->eta() << ", phi = " << pfTau->phi() << std::endl;
   }
 
   // collect the objects we are working with (ie tracks, tracks+gammas, etc)
   isoCharged_.clear();
   isoCharged_.reserve(pfTau->isolationPFChargedHadrCands().size());
   isoNeutral_.clear();
   isoNeutral_.reserve(pfTau->isolationPFGammaCands().size());
   isoPU_.clear();
   isoPU_.reserve(chargedPFCandidatesInEvent_.size());
 
   // Get the primary vertex associated to this tau
   reco::VertexRef pv = vertexAssociator_->associatedVertex(*pfTau);
   // Let the quality cuts know which the vertex to use when applying selections
   // on dz, etc.
   if ( verbosity_ ) {
     if ( pv.isNonnull() ) {
       std::cout << "pv: x = " << pv->position().x() << ", y = " << pv->position().y() << ", z = " << pv->position().z() << std::endl;
     } else {
       std::cout << "pv: N/A" << std::endl;
     }
     if ( pfTau->leadPFChargedHadrCand().isNonnull() ) {
       std::cout << "leadPFChargedHadron:" 
         << " Pt = " << pfTau->leadPFChargedHadrCand()->pt() << "," 
         << " eta = " << pfTau->leadPFChargedHadrCand()->eta() << "," 
         << " phi = " << pfTau->leadPFChargedHadrCand()->phi() << std::endl;
     } else {
       std::cout << "leadPFChargedHadron: N/A" << std::endl; 
     }
   }
 
   // CV: isolation is not well defined in case primary vertex or leading charged hadron do not exist
   if ( !(pv.isNonnull() && pfTau->leadPFChargedHadrCand().isNonnull()) ) return 0.;
 
   qcuts_->setPV(pv);
   qcuts_->setLeadTrack(pfTau->leadPFChargedHadrCand());
 
   if ( applyDeltaBeta_ ) {
     pileupQcutsGeneralQCuts_->setPV(pv);
     pileupQcutsGeneralQCuts_->setLeadTrack(pfTau->leadPFChargedHadrCand());
     pileupQcutsPUTrackSelection_->setPV(pv);
     pileupQcutsPUTrackSelection_->setLeadTrack(pfTau->leadPFChargedHadrCand());
   }
 
   // Load the tracks if they are being used.
   if ( includeTracks_ ) {
     BOOST_FOREACH( const reco::PFCandidatePtr& cand, pfTau->isolationPFChargedHadrCands() ) {
       if ( qcuts_->filterCandRef(cand) ) {
         isoCharged_.push_back(cand);
       }
     }
   }
   if ( includeGammas_ ) {
     BOOST_FOREACH( const reco::PFCandidatePtr& cand, pfTau->isolationPFGammaCands() ) {
       if ( qcuts_->filterCandRef(cand) ) {
         isoNeutral_.push_back(cand);
       }
     }
   }
 
   typedef reco::tau::cone::DeltaRPtrFilter<PFCandidatePtr> DRFilter;
 
   // If desired, get PU tracks.
   if ( applyDeltaBeta_ ) {
     // First select by inverted the DZ/track weight cuts. True = invert
     //if ( verbosity_ ) {
     //  std::cout << "Initial PFCands: " << chargedPFCandidatesInEvent_.size() << std::endl;
     //}
 
     std::vector<PFCandidatePtr> allPU =
       pileupQcutsPUTrackSelection_->filterCandRefs(
           chargedPFCandidatesInEvent_, true);
     //if ( verbosity_ ) {
     //  std::cout << "After track cuts: " << allPU.size() << std::endl;
     //}
 
     // Now apply the rest of the cuts, like pt, and TIP, tracker hits, etc
     std::vector<PFCandidatePtr> cleanPU =
       pileupQcutsGeneralQCuts_->filterCandRefs(allPU);
     //if ( verbosity_ ) {
     //  std::cout << "After cleaning cuts: " << cleanPU.size() << std::endl;
     //}
 
     // Only select PU tracks inside the isolation cone.
     DRFilter deltaBetaFilter(pfTau->p4(), 0, deltaBetaCollectionCone_);
     BOOST_FOREACH(const reco::PFCandidatePtr& cand, cleanPU) {
       if ( deltaBetaFilter(cand) ) {
         isoPU_.push_back(cand);
       }
     }
     //if ( verbosity_ ) {
     //  std::cout << "After cone cuts: " << isoPU.size() << std::endl;
     //}
   }
 
   // Check if we want a custom iso cone
   if ( customIsoCone_ >= 0. ) {
     //std::cout << "<PFRecoTauDiscriminationByIsolation::discriminate (moduleLabel = " << moduleLabel_ <<")>:" << std::endl;
     //std::cout << " customIsoCone = " << customIsoCone_ << std::endl;
     DRFilter filter(pfTau->p4(), 0, customIsoCone_);
     std::vector<PFCandidatePtr> isoCharged_filter;
     std::vector<PFCandidatePtr> isoNeutral_filter;
     // Remove all the objects not in our iso cone
     BOOST_FOREACH( const PFCandidatePtr& isoObject, isoCharged_ ) {
       if ( filter(isoObject) ) isoCharged_filter.push_back(isoObject);
     }
     BOOST_FOREACH( const PFCandidatePtr& isoObject, isoNeutral_ ) {
       if ( filter(isoObject) ) isoNeutral_filter.push_back(isoObject);
     }
     isoCharged_ = isoCharged_filter;
     isoNeutral_ = isoNeutral_filter;
   }
 
   bool failsOccupancyCut     = false;
   bool failsSumPtCut         = false;
   bool failsRelativeSumPtCut = false;
 
 //--- nObjects requirement
   int neutrals = isoNeutral_.size();
 
   if ( applyDeltaBeta_ ) {
     neutrals -= TMath::Nint(deltaBetaFactorThisEvent_*isoPU_.size());
   }
   if ( neutrals < 0 ) {
     neutrals = 0;
   }
 
   size_t nOccupants = isoCharged_.size() + neutrals;
 
   failsOccupancyCut = ( nOccupants > maximumOccupancy_ );
 
   double totalPt = 0.0;
   double puPt = 0.0;
 //--- Sum PT requirement
   if ( applySumPtCut_ || applyRelativeSumPtCut_ || storeRawSumPt_ || storeRawPUsumPt_ ) {
     double chargedPt = 0.0;
     double neutralPt = 0.0;
     BOOST_FOREACH ( const PFCandidatePtr& isoObject, isoCharged_ ) {
       chargedPt += isoObject->pt();
     }
     BOOST_FOREACH ( const PFCandidatePtr& isoObject, isoNeutral_ ) {
       neutralPt += isoObject->pt();
     }
     BOOST_FOREACH ( const PFCandidatePtr& isoObject, isoPU_ ) {
       puPt += isoObject->pt();
     }
     if ( verbosity_ ) {
       std::cout << "chargedPt = " << chargedPt << std::endl;
       std::cout << "neutralPt = " << neutralPt << std::endl;
       std::cout << "puPt = " << puPt << " (delta-beta corr. = " << (deltaBetaFactorThisEvent_*puPt) << ")" << std::endl;
     }
 
     if ( applyDeltaBeta_ ) {
       neutralPt -= deltaBetaFactorThisEvent_*puPt;
     }
 
     if ( applyRhoCorrection_ ) {
       neutralPt -= rhoThisEvent_;
     }
 
     if ( neutralPt < 0.0 ) {
       neutralPt = 0.0;
     }
 
     totalPt = chargedPt + neutralPt;
     if ( verbosity_ ) {
       std::cout << "totalPt = " << totalPt << " (cut = " << maximumSumPt_ << ")" << std::endl;
     }
 
     failsSumPtCut = (totalPt > maximumSumPt_);
 
 //--- Relative Sum PT requirement
     failsRelativeSumPtCut = (totalPt > (pfTau->pt()*maximumRelativeSumPt_));
   }
 
   bool fails = (applyOccupancyCut_ && failsOccupancyCut) ||
     (applySumPtCut_ && failsSumPtCut) ||
     (applyRelativeSumPtCut_ && failsRelativeSumPtCut);
 
   // We did error checking in the constructor, so this is safe.
   if ( storeRawSumPt_ ) {
     return totalPt;
   } else if ( storeRawPUsumPt_ ) {
     if ( applyDeltaBeta_ ) return puPt;
     else if ( applyRhoCorrection_ ) return rhoThisEvent_;
     else return 0.;
   } else if ( storeRawOccupancy_ ) {
     return nOccupants;
   } else {
     return (fails ? 0. : 1.);
   }
 }