Inheritance diagram for reco::tau::RecoTauBuilderConePlugin:

Public Member Functions
return_type	operator() (const reco::JetBaseRef &jet, const std::vector< reco::PFRecoTauChargedHadron > &chargedHadrons, const std::vector< RecoTauPiZero > &piZeros, const std::vector< CandidatePtr > &regionalExtras) const override

	RecoTauBuilderConePlugin (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC)

	~RecoTauBuilderConePlugin () override

Public Member Functions inherited from reco::tau::RecoTauBuilderPlugin
void	beginEvent () override

const edm::Handle< edm::View< reco::Candidate > > &	getPFCands () const
	Hack to be able to convert Ptrs to Refs. More...

reco::VertexRef	primaryVertex (const reco::JetBaseRef &jet) const
	Get primary vertex associated to this jet. More...

reco::VertexRef	primaryVertex (const reco::PFTau &tau, bool useJet=false) const
	Get primary vertex associated to this tau. More...

	RecoTauBuilderPlugin (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC)

	~RecoTauBuilderPlugin () override

Public Member Functions inherited from reco::tau::RecoTauEventHolderPlugin
edm::Event *	evt ()

const edm::Event *	evt () const

const edm::EventSetup *	evtSetup () const

	RecoTauEventHolderPlugin (const edm::ParameterSet &pset)

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

	~RecoTauEventHolderPlugin () override

Public Member Functions inherited from reco::tau::RecoTauNamedPlugin
const std::string &	name () const

	RecoTauNamedPlugin (const edm::ParameterSet &pset)

virtual	~RecoTauNamedPlugin ()

Private Types
typedef StringObjectFunction< reco::Jet >	JetFunc

Private Member Functions
void	setTauQuantities (reco::PFTau &aTau, double minAbsPhotonSumPt_insideSignalCone=2.5, double minRelPhotonSumPt_insideSignalCone=0., double minAbsPhotonSumPt_outsideSignalCone=1.e+9, double minRelPhotonSumPt_outsideSignalCone=1.e+9) const

Private Attributes
JetFunc	isoConeChargedHadrons_

JetFunc	isoConeNeutralHadrons_

JetFunc	isoConePiZeros_

double	leadObjecPtThreshold_

JetFunc	matchingCone_

int	maxSignalConeChargedHadrons_

double	minAbsPhotonSumPt_insideSignalCone_

double	minRelPhotonSumPt_insideSignalCone_

std::unique_ptr< RecoTauQualityCuts >	qcuts_

JetFunc	signalConeChargedHadrons_

JetFunc	signalConeNeutralHadrons_

JetFunc	signalConePiZeros_

StringObjectFunction< reco::PFTau >	signalConeSizeToStore_

bool	usePFLeptonsAsChargedHadrons_

Additional Inherited Members
Public Types inherited from reco::tau::RecoTauBuilderPlugin
typedef boost::ptr_vector< reco::PFTau >	output_type

typedef std::unique_ptr< output_type >	return_type

Detailed Description

Definition at line 35 of file RecoTauBuilderConePlugin.cc.

Member Typedef Documentation

◆ JetFunc

typedef StringObjectFunction<reco::Jet> reco::tau::RecoTauBuilderConePlugin::JetFunc

private

Definition at line 53 of file RecoTauBuilderConePlugin.cc.

Constructor & Destructor Documentation

◆ RecoTauBuilderConePlugin()

reco::tau::RecoTauBuilderConePlugin::RecoTauBuilderConePlugin	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &&	iC
	)

explicit

Definition at line 77 of file RecoTauBuilderConePlugin.cc.

         : RecoTauBuilderPlugin(pset, std::move(iC)),
           qcuts_(new RecoTauQualityCuts(pset.getParameterSet("qualityCuts").getParameterSet("signalQualityCuts"))),
           usePFLeptonsAsChargedHadrons_(pset.getParameter<bool>("usePFLeptons")),
           leadObjecPtThreshold_(pset.getParameter<double>("leadObjectPt")),
           matchingCone_(pset.getParameter<std::string>("matchingCone")),
           signalConeChargedHadrons_(pset.getParameter<std::string>("signalConeChargedHadrons")),
           isoConeChargedHadrons_(pset.getParameter<std::string>("isoConeChargedHadrons")),
           signalConePiZeros_(pset.getParameter<std::string>("signalConePiZeros")),
           signalConeSizeToStore_(  //MB: stored inside PFTau and used to compte photon-pt-out-of-signal-cone and DM hypothsis
               pset.getParameter<std::string>("signalConePiZeros")),
           isoConePiZeros_(pset.getParameter<std::string>("isoConePiZeros")),
           signalConeNeutralHadrons_(pset.getParameter<std::string>("signalConeNeutralHadrons")),
           isoConeNeutralHadrons_(pset.getParameter<std::string>("isoConeNeutralHadrons")),
           maxSignalConeChargedHadrons_(pset.getParameter<int>("maxSignalConeChargedHadrons")),
           minAbsPhotonSumPt_insideSignalCone_(pset.getParameter<double>("minAbsPhotonSumPt_insideSignalCone")),
           minRelPhotonSumPt_insideSignalCone_(pset.getParameter<double>("minRelPhotonSumPt_insideSignalCone"))
  
     {}

◆ ~RecoTauBuilderConePlugin()

reco::tau::RecoTauBuilderConePlugin::~RecoTauBuilderConePlugin ( )

inlineoverride

Definition at line 38 of file RecoTauBuilderConePlugin.cc.

38 {}

Member Function Documentation

◆ operator()()

RecoTauBuilderConePlugin::return_type reco::tau::RecoTauBuilderConePlugin::operator()	(	const reco::JetBaseRef &	,
		const std::vector< reco::PFRecoTauChargedHadron > &	,
		const std::vector< RecoTauPiZero > &	,
		const std::vector< CandidatePtr > &
	)		const

overridevirtual

Construct one or more PFTaus from the a PFJet and its asscociated reconstructed PiZeros and regional extras i.e. objects in a 0.8 cone about the jet

Implements reco::tau::RecoTauBuilderPlugin.

Definition at line 179 of file RecoTauBuilderConePlugin.cc.

                                                              {
       //std::cout << "<RecoTauBuilderConePlugin::operator()>:" << std::endl;
       //std::cout << " jet: Pt = " << jet->pt() << ", eta = " << jet->eta() << ", phi = " << jet->phi() << std::endl;
  
       // Get access to our cone tools
       using namespace cone;
       // Define output.  We only produce one tau per jet for the cone algo.
       output_type output;
  
       // Our tau builder - the true indicates to automatically copy gamma candidates
       // from the pizeros.
       RecoTauConstructor tau(jet,
                              getPFCands(),
                              true,
                              &signalConeSizeToStore_,
                              minAbsPhotonSumPt_insideSignalCone_,
                              minRelPhotonSumPt_insideSignalCone_);
       // Setup our quality cuts to use the current vertex (supplied by base class)
       qcuts_->setPV(primaryVertex(jet));
  
       typedef std::vector<CandidatePtr> CandPtrs;
  
       // Get the PF Charged hadrons + quality cuts
       CandPtrs pfchs;
       if (!usePFLeptonsAsChargedHadrons_) {
         pfchs = qcuts_->filterCandRefs(pfCandidatesByPdgId(*jet, 211));
       } else {
         // Check if we want to include electrons in muons in "charged hadron"
         // collection.  This is the preferred behavior, as the PF lepton selections
         // are very loose.
         pfchs = qcuts_->filterCandRefs(pfChargedCands(*jet));
       }
  
       // CV: sort collection of PF Charged hadrons by descending Pt
       std::sort(pfchs.begin(), pfchs.end(), SortPFCandsDescendingPt());
  
       // Get the PF gammas
       CandPtrs pfGammas = qcuts_->filterCandRefs(pfCandidatesByPdgId(*jet, 22));
       // Neutral hadrons
       CandPtrs pfnhs = qcuts_->filterCandRefs(pfCandidatesByPdgId(*jet, 130));
  
       // All the extra junk
       CandPtrs regionalJunk = qcuts_->filterCandRefs(regionalExtras);
  
       /***********************************************
    ******     Lead Candidate Finding    **********
    ***********************************************/
  
       // Define our matching cone and filters
       double matchingCone = matchingCone_(*jet);
       CandPtrDRFilter matchingConeFilter(jet->p4(), 0, matchingCone);
  
       // Find the maximum PFCharged hadron in the matching cone.  The call to
       // PFCandidates always a sorted list, so we can just take the first if it
       // if it exists.
       CandidatePtr leadPFCH;
       CandPtrs::iterator leadPFCH_iter = std::find_if(pfchs.begin(), pfchs.end(), matchingConeFilter);
  
       if (leadPFCH_iter != pfchs.end()) {
         leadPFCH = *leadPFCH_iter;
         // Set leading candidate
         tau.setleadChargedHadrCand(leadPFCH);
       } else {
         // If there is no leading charged candidate at all, return nothing - the
         // producer class that owns the plugin will build a null tau if desired.
         return output.release();
       }
  
       // Find the leading neutral candidate
       CandidatePtr leadPFGamma;
       CandPtrs::iterator leadPFGamma_iter = std::find_if(pfGammas.begin(), pfGammas.end(), matchingConeFilter);
  
       if (leadPFGamma_iter != pfGammas.end()) {
         leadPFGamma = *leadPFGamma_iter;
         // Set leading neutral candidate
         tau.setleadNeutralCand(leadPFGamma);
       }
  
       CandidatePtr leadPFCand;
       // Always use the leadPFCH if it is above our threshold
       if (leadPFCH.isNonnull() && leadPFCH->pt() > leadObjecPtThreshold_) {
         leadPFCand = leadPFCH;
       } else if (leadPFGamma.isNonnull() && leadPFGamma->pt() > leadObjecPtThreshold_) {
         // Otherwise use the lead Gamma if it is above threshold
         leadPFCand = leadPFGamma;
       } else {
         // If both are too low PT, just take the charged one
         leadPFCand = leadPFCH;
       }
  
       tau.setleadCand(leadPFCand);
  
       // Our cone axis is defined about the lead charged hadron
       reco::Candidate::LorentzVector coneAxis = leadPFCH->p4();
  
       /***********************************************
    ******     Cone Construction         **********
    ***********************************************/
  
       // Define the signal and isolation cone sizes for this jet and build filters
       // to select elements in the given DeltaR regions
  
       CandPtrDRFilter signalConePFCHFilter(coneAxis, -0.1, signalConeChargedHadrons_(*jet));
       CandPtrDRFilter signalConePFNHFilter(coneAxis, -0.1, signalConeNeutralHadrons_(*jet));
       PiZeroDRFilter signalConePiZeroFilter(coneAxis, -0.1, signalConePiZeros_(*jet));
  
       CandPtrDRFilter isoConePFCHFilter(coneAxis, signalConeChargedHadrons_(*jet), isoConeChargedHadrons_(*jet));
       CandPtrDRFilter isoConePFGammaFilter(coneAxis, signalConePiZeros_(*jet), isoConePiZeros_(*jet));
       CandPtrDRFilter isoConePFNHFilter(coneAxis, signalConeNeutralHadrons_(*jet), isoConeNeutralHadrons_(*jet));
       PiZeroDRFilter isoConePiZeroFilter(coneAxis, signalConePiZeros_(*jet), isoConePiZeros_(*jet));
  
       // Additionally make predicates to select the different PF object types
       // of the regional junk objects to add to the iso cone.
       typedef xclean::PredicateAND<xclean::FilterCandByAbsPdgId, CandPtrDRFilter> RegionalJunkConeAndIdFilter;
  
       xclean::FilterCandByAbsPdgId pfchCandSelector(211);
       xclean::FilterCandByAbsPdgId pfgammaCandSelector(22);
       xclean::FilterCandByAbsPdgId pfnhCandSelector(130);
  
       // Predicate to select the regional junk in the iso cone by PF id
       RegionalJunkConeAndIdFilter pfChargedJunk(pfchCandSelector,  // select charged stuff from junk
                                                 isoConePFCHFilter  // only take those in iso cone
       );
  
       RegionalJunkConeAndIdFilter pfGammaJunk(pfgammaCandSelector,  // select gammas from junk
                                               isoConePFGammaFilter  // only take those in iso cone
       );
  
       RegionalJunkConeAndIdFilter pfNeutralJunk(pfnhCandSelector,  // select neutral stuff from junk
                                                 isoConePFNHFilter  // select stuff in iso cone
       );
  
       // Build filter iterators select the signal charged stuff.
       CandPtrDRFilterIter signalPFCHCands_begin(signalConePFCHFilter, pfchs.begin(), pfchs.end());
       CandPtrDRFilterIter signalPFCHCands_end(signalConePFCHFilter, pfchs.end(), pfchs.end());
       CandPtrs signalPFCHs;
       int numSignalPFCHs = 0;
       CandPtrs isolationPFCHs;
       int numIsolationPFCHs = 0;
       for (CandPtrDRFilterIter iter = signalPFCHCands_begin; iter != signalPFCHCands_end; ++iter) {
         if (numSignalPFCHs < maxSignalConeChargedHadrons_ || maxSignalConeChargedHadrons_ == -1) {
           //std::cout << "adding signalPFCH #" << numSignalPFCHs << ": Pt = " << (*iter)->pt() << ", eta = " << (*iter)->eta() << ", phi = " << (*iter)->phi() << std::endl;
           signalPFCHs.push_back(*iter);
           ++numSignalPFCHs;
         } else {
           //std::cout << "maxSignalConeChargedHadrons reached"
           //      << " --> adding isolationPFCH #" << numIsolationPFCHs << ": Pt = " << (*iter)->pt() << ", eta = " << (*iter)->eta() << ", phi = " << (*iter)->phi() << std::endl;
           isolationPFCHs.push_back(*iter);
           ++numIsolationPFCHs;
         }
       }
       CandPtrs::const_iterator signalPFCHs_begin = signalPFCHs.begin();
       CandPtrs::const_iterator signalPFCHs_end = signalPFCHs.end();
  
       // Cross clean pi zero content using signal cone charged hadron constituents.
       xclean::CrossCleanPiZeros<CandPtrDRFilterIter> piZeroXCleaner(signalPFCHCands_begin, signalPFCHCands_end);
       std::vector<reco::RecoTauPiZero> cleanPiZeros = piZeroXCleaner(piZeros);
  
       // For the rest of the constituents, we need to filter any constituents that
       // are already contained in the pizeros (i.e. electrons)
       xclean::CrossCleanPtrs<PiZeroList::const_iterator> pfCandXCleaner(cleanPiZeros.begin(), cleanPiZeros.end());
  
       auto isolationPFCHCands_begin(boost::make_filter_iterator(
           xclean::makePredicateAND(isoConePFCHFilter, pfCandXCleaner), pfchs.begin(), pfchs.end()));
       auto isolationPFCHCands_end(boost::make_filter_iterator(
           xclean::makePredicateAND(isoConePFCHFilter, pfCandXCleaner), pfchs.end(), pfchs.end()));
       for (auto iter = isolationPFCHCands_begin; iter != isolationPFCHCands_end; ++iter) {
         //std::cout << "adding isolationPFCH #" << numIsolationPFCHs << ": Pt = " << (*iter)->pt() << ", eta = " << (*iter)->eta() << ", phi = " << (*iter)->phi() << std::endl;
         isolationPFCHs.push_back(*iter);
         ++numIsolationPFCHs;
       }
       CandPtrs::const_iterator isolationPFCHs_begin = isolationPFCHs.begin();
       CandPtrs::const_iterator isolationPFCHs_end = isolationPFCHs.end();
  
       // Build signal charged hadrons
       tau.addPFCands(
           RecoTauConstructor::kSignal, RecoTauConstructor::kChargedHadron, signalPFCHs_begin, signalPFCHs_end);
  
       tau.addPFCands(RecoTauConstructor::kSignal,
                      RecoTauConstructor::kNeutralHadron,
                      boost::make_filter_iterator(
                          xclean::makePredicateAND(signalConePFNHFilter, pfCandXCleaner), pfnhs.begin(), pfnhs.end()),
                      boost::make_filter_iterator(
                          xclean::makePredicateAND(signalConePFNHFilter, pfCandXCleaner), pfnhs.end(), pfnhs.end()));
  
       // Build signal PiZeros
       tau.addPiZeros(RecoTauConstructor::kSignal,
                      PiZeroDRFilterIter(signalConePiZeroFilter, cleanPiZeros.begin(), cleanPiZeros.end()),
                      PiZeroDRFilterIter(signalConePiZeroFilter, cleanPiZeros.end(), cleanPiZeros.end()));
  
       // Build isolation charged hadrons
       tau.addPFCands(
           RecoTauConstructor::kIsolation, RecoTauConstructor::kChargedHadron, isolationPFCHs_begin, isolationPFCHs_end);
  
       // Add all the stuff in the isolation cone that wasn't in the jet constituents
       tau.addPFCands(RecoTauConstructor::kIsolation,
                      RecoTauConstructor::kChargedHadron,
                      boost::make_filter_iterator(pfChargedJunk, regionalJunk.begin(), regionalJunk.end()),
                      boost::make_filter_iterator(pfChargedJunk, regionalJunk.end(), regionalJunk.end()));
  
       // Build isolation neutral hadrons
       tau.addPFCands(RecoTauConstructor::kIsolation,
                      RecoTauConstructor::kNeutralHadron,
                      boost::make_filter_iterator(
                          xclean::makePredicateAND(isoConePFNHFilter, pfCandXCleaner), pfnhs.begin(), pfnhs.end()),
                      boost::make_filter_iterator(
                          xclean::makePredicateAND(isoConePFNHFilter, pfCandXCleaner), pfnhs.end(), pfnhs.end()));
  
       // Add regional stuff not in jet
       tau.addPFCands(RecoTauConstructor::kIsolation,
                      RecoTauConstructor::kNeutralHadron,
                      boost::make_filter_iterator(pfNeutralJunk, regionalJunk.begin(), regionalJunk.end()),
                      boost::make_filter_iterator(pfNeutralJunk, regionalJunk.end(), regionalJunk.end()));
  
       // Build isolation PiZeros
       tau.addPiZeros(RecoTauConstructor::kIsolation,
                      PiZeroDRFilterIter(isoConePiZeroFilter, cleanPiZeros.begin(), cleanPiZeros.end()),
                      PiZeroDRFilterIter(isoConePiZeroFilter, cleanPiZeros.end(), cleanPiZeros.end()));
  
       // Add regional stuff not in jet
       tau.addPFCands(RecoTauConstructor::kIsolation,
                      RecoTauConstructor::kGamma,
                      boost::make_filter_iterator(pfGammaJunk, regionalJunk.begin(), regionalJunk.end()),
                      boost::make_filter_iterator(pfGammaJunk, regionalJunk.end(), regionalJunk.end()));
  
       // Put our built tau in the output - 'false' indicates don't build the
       // leading candidates, we already did that explicitly above.
  
       std::unique_ptr<reco::PFTau> tauPtr = tau.get(false);
  
       // Set event vertex position for tau
       reco::VertexRef primaryVertexRef = primaryVertex(*tauPtr);
       if (primaryVertexRef.isNonnull())
         tauPtr->setVertex(primaryVertexRef->position());
  
       // Set missing tau quantities
       setTauQuantities(*tauPtr, minAbsPhotonSumPt_insideSignalCone_, minRelPhotonSumPt_insideSignalCone_);
  
       output.push_back(std::move(tauPtr));

References reco::tau::RecoTauBuilderPlugin::getPFCands(), edm::Ptr< T >::isNonnull(), edm::Ref< C, T, F >::isNonnull(), isoConeChargedHadrons_, isoConeNeutralHadrons_, isoConePiZeros_, metsig::jet, reco::tau::RecoTauConstructor::kChargedHadron, reco::tau::RecoTauConstructor::kGamma, reco::tau::RecoTauConstructor::kIsolation, reco::tau::RecoTauConstructor::kNeutralHadron, reco::tau::RecoTauConstructor::kSignal, leadObjecPtThreshold_, reco::tau::xclean::makePredicateAND(), HLT_FULL_cff::matchingCone, matchingCone_, maxSignalConeChargedHadrons_, minAbsPhotonSumPt_insideSignalCone_, minRelPhotonSumPt_insideSignalCone_, eostools::move(), convertSQLitetoXML_cfg::output, reco::tau::pfCandidatesByPdgId(), reco::tau::pfChargedCands(), reco::tau::pfGammas(), reco::tau::RecoTauBuilderPlugin::primaryVertex(), qcuts_, setTauQuantities(), signalConeChargedHadrons_, signalConeNeutralHadrons_, signalConePiZeros_, signalConeSizeToStore_, jetUpdater_cfi::sort, metsig::tau, and usePFLeptonsAsChargedHadrons_.

◆ setTauQuantities()

void reco::tau::RecoTauBuilderConePlugin::setTauQuantities	(	reco::PFTau &	aTau,
		double	minAbsPhotonSumPt_insideSignalCone = `2.5`,
		double	minRelPhotonSumPt_insideSignalCone = `0.`,
		double	minAbsPhotonSumPt_outsideSignalCone = `1.e+9`,
		double	minRelPhotonSumPt_outsideSignalCone = `1.e+9`
	)		const

private

Definition at line 117 of file RecoTauBuilderConePlugin.cc.

                                                                                                       {
       //MB: Set tau quantities which are computed by RecoTauConstructor::get()
       //    method with PFRecoTauChargedHadrons not used here
  
       // Set charge
       int charge = 0;
       int leadCharge = aTau.leadChargedHadrCand().isNonnull() ? aTau.leadChargedHadrCand()->charge() : 0;
       const std::vector<reco::CandidatePtr>& pfChs = aTau.signalChargedHadrCands();
       for (const reco::CandidatePtr& pfCh : pfChs) {
         charge += pfCh->charge();
       }
       charge = charge == 0 ? leadCharge : charge;
       aTau.setCharge(charge);
  
       // Set PDG id
       aTau.setPdgId(aTau.charge() < 0 ? 15 : -15);
  
       // Set Decay Mode
       PFTau::hadronicDecayMode dm = PFTau::kNull;
       unsigned int nPiZeros = 0;
       const std::vector<RecoTauPiZero>& piZeros = aTau.signalPiZeroCandidates();
       for (const RecoTauPiZero& piZero : piZeros) {
         double photonSumPt_insideSignalCone = 0.;
         double photonSumPt_outsideSignalCone = 0.;
         int numPhotons = piZero.numberOfDaughters();
         for (int idxPhoton = 0; idxPhoton < numPhotons; ++idxPhoton) {
           const reco::Candidate* photon = piZero.daughter(idxPhoton);
           double dR = deltaR(photon->p4(), aTau.p4());
           if (dR < aTau.signalConeSize()) {
             photonSumPt_insideSignalCone += photon->pt();
           } else {
             photonSumPt_outsideSignalCone += photon->pt();
           }
         }
         if (photonSumPt_insideSignalCone > minAbsPhotonSumPt_insideSignalCone ||
             photonSumPt_insideSignalCone > (minRelPhotonSumPt_insideSignalCone * aTau.pt()) ||
             photonSumPt_outsideSignalCone > minAbsPhotonSumPt_outsideSignalCone ||
             photonSumPt_outsideSignalCone > (minRelPhotonSumPt_outsideSignalCone * aTau.pt()))
           ++nPiZeros;
       }
       // Find the maximum number of PiZeros our parameterization can hold
       const unsigned int maxPiZeros = PFTau::kOneProngNPiZero;
       // Determine our track index
       unsigned int nCharged = pfChs.size();
       if (nCharged > 0) {
         unsigned int trackIndex = (nCharged - 1) * (maxPiZeros + 1);
         // Check if we handle the given number of tracks
         if (trackIndex >= PFTau::kRareDecayMode)
           dm = PFTau::kRareDecayMode;
         else
           dm = static_cast<PFTau::hadronicDecayMode>(trackIndex + std::min(nPiZeros, maxPiZeros));
       } else {
         dm = PFTau::kNull;
       }
       aTau.setDecayMode(dm);
       return;
     }