#include <PFTauExtractor.h>

Inheritance diagram for PFTauExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &evt, const edm::EventSetup &es, const reco::Candidate &candidate) const override

reco::IsoDeposit	deposit (const edm::Event &evt, const edm::EventSetup &es, const reco::Track &track) const override

void	fillVetos (const edm::Event &, const edm::EventSetup &, const reco::TrackCollection &) override
	definition of pure virtual functions inherited from IsoDepositExtractor base-class More...

	PFTauExtractor (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

	~PFTauExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const

virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Member Functions
template<typename T >
reco::IsoDeposit	depositFromObject (const edm::Event &, const edm::EventSetup &, const T &) const

Private Attributes
edm::EDGetTokenT< edm::View< reco::Candidate > >	candidateSourceToken_

double	dRIsoCone_

double	dRmatchPFTau_

double	dRVetoCone_

double	dRvetoPFTauSignalConeConstituents_

double	maxDxyTrack_

double	maxDzTrack_

edm::EDGetTokenT< reco::PFTauCollection >	tauSourceToken_
	configuration parameters More...

Detailed Description

Definition at line 20 of file PFTauExtractor.h.

Constructor & Destructor Documentation

◆ PFTauExtractor()

PFTauExtractor::PFTauExtractor	(	const edm::ParameterSet &	cfg,
		edm::ConsumesCollector &&	iC
	)

explicit

Definition at line 12 of file PFTauExtractor.cc.

                                                                                     {
   tauSourceToken_ = iC.consumes<reco::PFTauCollection>(cfg.getParameter<edm::InputTag>("tauSource"));
   candidateSourceToken_ =
       iC.mayConsume<edm::View<reco::Candidate> >(cfg.getParameter<edm::InputTag>("candidateSource"));
   maxDxyTrack_ = cfg.getParameter<double>("Diff_r");
   maxDzTrack_ = cfg.getParameter<double>("Diff_z");
   dRmatchPFTau_ = cfg.getParameter<double>("dRmatchPFTau");
   dRVetoCone_ = cfg.getParameter<double>("DR_Veto");
   dRIsoCone_ = cfg.getParameter<double>("DR_Max");
   dRvetoPFTauSignalConeConstituents_ = cfg.getParameter<double>("dRvetoPFTauSignalConeConstituents");
 }

References candidateSourceToken_, looper::cfg, dRIsoCone_, dRmatchPFTau_, dRVetoCone_, dRvetoPFTauSignalConeConstituents_, maxDxyTrack_, maxDzTrack_, and tauSourceToken_.

◆ ~PFTauExtractor()

PFTauExtractor::~PFTauExtractor ( )

inlineoverride

Definition at line 23 of file PFTauExtractor.h.

23 {}

Member Function Documentation

◆ deposit() [1/2]

reco::IsoDeposit PFTauExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Candidate &	track
	)		const

inlineoverridevirtual

make single IsoDeposit based on a candidate as input purely virtual: have to implement in concrete implementations

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 30 of file PFTauExtractor.h.

                                                                           {
     return depositFromObject(evt, es, candidate);
   }

References depositFromObject().

◆ deposit() [2/2]

reco::IsoDeposit PFTauExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

inlineoverridevirtual

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 27 of file PFTauExtractor.h.

                                                                                                               {
     return depositFromObject(evt, es, track);
   }

References depositFromObject(), and HLT_FULL_cff::track.

◆ depositFromObject()

template<typename T >

reco::IsoDeposit PFTauExtractor::depositFromObject	(	const edm::Event &	evt,
		const edm::EventSetup &	es,
		const T &	tauCandidate
	)		const

private

private member function for computing the IsoDeposits in case of reco::Track as well as in case of reco::Canididate input

Definition at line 25 of file PFTauExtractor.cc.

                                                                                 {
   //--- create IsoDeposit to be returned;
   //    set "direction" of IsoDeposit cone to jet-axis of tauCandidate
   reco::isodeposit::Direction tauCandidateDirection(tauCandidate.eta(), tauCandidate.phi());
   reco::IsoDeposit isoDeposit(tauCandidateDirection);
   isoDeposit.addCandEnergy(tauCandidate.pt());
  
   //--- find PFTau closest to tauDirection
   edm::Handle<reco::PFTauCollection> pfTaus;
   evt.getByToken(tauSourceToken_, pfTaus);
  
   double dR_min = -1.;
   const reco::PFTau* pfTau_matched = nullptr;
   for (reco::PFTauCollection::const_iterator pfTau = pfTaus->begin(); pfTau != pfTaus->end(); ++pfTau) {
     double dR = deltaR(pfTau->eta(), pfTau->phi(), tauCandidate.eta(), tauCandidate.phi());
     if (pfTau_matched == nullptr || dR < dR_min) {
       dR_min = dR;
       pfTau_matched = &(*pfTau);
     }
   }
  
   //--- compute IsoDeposit for matched PFTau
   if (pfTau_matched != nullptr && dR_min < dRmatchPFTau_) {
     edm::Handle<edm::View<reco::Candidate> > candidates;
     evt.getByToken(candidateSourceToken_, candidates);
  
     const reco::Particle::Point& tauVertex = pfTau_matched->vertex();
     double dRsignalCone_max = 0.;
     for (edm::View<reco::Candidate>::const_iterator candidate = candidates->begin(); candidate != candidates->end();
          ++candidate) {
       double dR = deltaR(candidate->momentum(), pfTau_matched->momentum());
  
       //--- check that candidate is inbetween veto and isolation cone,
       //    and is compatible with originating from the same primary event vertex as the PFTau
       if (dR > dRVetoCone_ && dR < dRIsoCone_ && (candidate->vertex() - tauVertex).Rho() < maxDxyTrack_ &&
           fabs(candidate->vertex().z() - tauVertex.z()) < maxDzTrack_) {
         //--- check that the candidate is not associated to one of the tau decay products
         //    within the signal cone of the PFTau
         bool isSignalCone = false;
         for (std::vector<reco::CandidatePtr>::const_iterator tauSignalConeConstituent =
                  pfTau_matched->signalCands().begin();
              tauSignalConeConstituent != pfTau_matched->signalCands().end();
              ++tauSignalConeConstituent) {
           double dR = deltaR(candidate->momentum(), (*tauSignalConeConstituent)->momentum());
           if (dR <= dRvetoPFTauSignalConeConstituents_)
             isSignalCone = true;
         }
  
         if (!isSignalCone) {
           reco::isodeposit::Direction candidateDirection(candidate->eta(), candidate->phi());
           isoDeposit.addDeposit(candidateDirection, candidate->pt());
         }
       }
     }
  
     //--- set size of veto cone of IsoDeposit to largest distance
     //    of any tau decay product within the signal cone of the PFTau
     //    (add a small positive number in order to avoid issues
     //     with rounding errors and "<" versus "<=" comparisson)
     reco::IsoDeposit::Veto isoDepositVeto;
     isoDepositVeto.vetoDir = tauCandidateDirection;
     isoDepositVeto.dR = dRsignalCone_max + 1.e-3;
     isoDeposit.setVeto(isoDepositVeto);
   } else {
     edm::LogWarning("PFTauExtractor::depositFromObject")
         << " Failed to match PFTau to tauCandidate direction given by"
         << " eta = " << tauCandidate.eta() << ", phi = " << tauCandidate.phi()
         << " --> skipping computation of IsoDeposit !!";
   }
  
   return isoDeposit;
 }