Inheritance diagram for reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin:

Public Member Functions
void	beginEvent () override
	Hook called at the beginning of the event. More...

return_type	operator() (const reco::Jet &) const override
	Build a collection of chargedHadrons from objects in the input jet. More...

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

	~PFRecoTauChargedHadronFromPFCandidatePlugin () override

Public Member Functions inherited from reco::tau::PFRecoTauChargedHadronBuilderPlugin
	PFRecoTauChargedHadronBuilderPlugin (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC)

	~PFRecoTauChargedHadronBuilderPlugin () 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 std::vector< reco::CandidatePtr >	CandPtrs

Private Attributes
double	bField_

const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bFieldToken_

double	dRmergeNeutralHadronWrtChargedHadron_

double	dRmergeNeutralHadronWrtElectron_

double	dRmergeNeutralHadronWrtNeutralHadron_

double	dRmergeNeutralHadronWrtOther_

double	dRmergePhotonWrtChargedHadron_

double	dRmergePhotonWrtElectron_

double	dRmergePhotonWrtNeutralHadron_

double	dRmergePhotonWrtOther_

std::vector< int >	inputParticleIds_

int	maxUnmatchedBlockElementsNeutralHadron_

int	maxUnmatchedBlockElementsPhoton_

int	minBlockElementMatchesNeutralHadron_

int	minBlockElementMatchesPhoton_

double	minMergeChargedHadronPt_

double	minMergeGammaEt_

double	minMergeNeutralHadronEt_

RecoTauQualityCuts *	qcuts_

int	verbosity_

RecoTauVertexAssociator	vertexAssociator_

Additional Inherited Members
Public Types inherited from reco::tau::PFRecoTauChargedHadronBuilderPlugin
typedef boost::ptr_vector< PFRecoTauChargedHadron >	ChargedHadronVector

typedef std::unique_ptr< ChargedHadronVector >	return_type

Detailed Description

Definition at line 43 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Member Typedef Documentation

◆ CandPtrs

typedef std::vector<reco::CandidatePtr> reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::CandPtrs

private

Definition at line 53 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Constructor & Destructor Documentation

◆ PFRecoTauChargedHadronFromPFCandidatePlugin()

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

explicit

Definition at line 83 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

         : PFRecoTauChargedHadronBuilderPlugin(pset, std::move(iC)),
           vertexAssociator_(pset.getParameter<edm::ParameterSet>("qualityCuts"), std::move(iC)),
           qcuts_(nullptr),
           bFieldToken_(iC.esConsumes()) {
       edm::ParameterSet qcuts_pset = pset.getParameterSet("qualityCuts").getParameterSet("signalQualityCuts");
       qcuts_ = new RecoTauQualityCuts(qcuts_pset);
  
       inputParticleIds_ = pset.getParameter<std::vector<int> >("chargedHadronCandidatesParticleIds");
  
       dRmergeNeutralHadronWrtChargedHadron_ = pset.getParameter<double>("dRmergeNeutralHadronWrtChargedHadron");
       dRmergeNeutralHadronWrtNeutralHadron_ = pset.getParameter<double>("dRmergeNeutralHadronWrtNeutralHadron");
       dRmergeNeutralHadronWrtElectron_ = pset.getParameter<double>("dRmergeNeutralHadronWrtElectron");
       dRmergeNeutralHadronWrtOther_ = pset.getParameter<double>("dRmergeNeutralHadronWrtOther");
       minBlockElementMatchesNeutralHadron_ = pset.getParameter<int>("minBlockElementMatchesNeutralHadron");
       maxUnmatchedBlockElementsNeutralHadron_ = pset.getParameter<int>("maxUnmatchedBlockElementsNeutralHadron");
       dRmergePhotonWrtChargedHadron_ = pset.getParameter<double>("dRmergePhotonWrtChargedHadron");
       dRmergePhotonWrtNeutralHadron_ = pset.getParameter<double>("dRmergePhotonWrtNeutralHadron");
       dRmergePhotonWrtElectron_ = pset.getParameter<double>("dRmergePhotonWrtElectron");
       dRmergePhotonWrtOther_ = pset.getParameter<double>("dRmergePhotonWrtOther");
       minBlockElementMatchesPhoton_ = pset.getParameter<int>("minBlockElementMatchesPhoton");
       maxUnmatchedBlockElementsPhoton_ = pset.getParameter<int>("maxUnmatchedBlockElementsPhoton");
       minMergeNeutralHadronEt_ = pset.getParameter<double>("minMergeNeutralHadronEt");
       minMergeGammaEt_ = pset.getParameter<double>("minMergeGammaEt");
       minMergeChargedHadronPt_ = pset.getParameter<double>("minMergeChargedHadronPt");
  
       verbosity_ = pset.getParameter<int>("verbosity");
     }

References dRmergeNeutralHadronWrtChargedHadron_, dRmergeNeutralHadronWrtElectron_, dRmergeNeutralHadronWrtNeutralHadron_, dRmergeNeutralHadronWrtOther_, dRmergePhotonWrtChargedHadron_, dRmergePhotonWrtElectron_, dRmergePhotonWrtNeutralHadron_, dRmergePhotonWrtOther_, inputParticleIds_, maxUnmatchedBlockElementsNeutralHadron_, maxUnmatchedBlockElementsPhoton_, minBlockElementMatchesNeutralHadron_, minBlockElementMatchesPhoton_, minMergeChargedHadronPt_, minMergeGammaEt_, minMergeNeutralHadronEt_, muonDTDigis_cfi::pset, qcuts_, and verbosity_.

◆ ~PFRecoTauChargedHadronFromPFCandidatePlugin()

reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::~PFRecoTauChargedHadronFromPFCandidatePlugin ( )

override

Definition at line 113 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

113 { delete qcuts_; }

References qcuts_.

Member Function Documentation

◆ beginEvent()

void reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::beginEvent ( )

overridevirtual

Hook called at the beginning of the event.

Reimplemented from reco::tau::PFRecoTauChargedHadronBuilderPlugin.

Definition at line 116 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

                                                                  {
       vertexAssociator_.setEvent(*evt());
  
       bField_ = evtSetup()->getData(bFieldToken_).inTesla(GlobalPoint(0, 0, 0)).z();
     }

References bField_, bFieldToken_, reco::tau::RecoTauEventHolderPlugin::evt(), reco::tau::RecoTauEventHolderPlugin::evtSetup(), edm::EventSetup::getData(), reco::tau::RecoTauVertexAssociator::setEvent(), and vertexAssociator_.

◆ operator()()

PFRecoTauChargedHadronFromPFCandidatePlugin::return_type reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::operator() ( const reco::Jet & ) const

overridevirtual

Build a collection of chargedHadrons from objects in the input jet.

Implements reco::tau::PFRecoTauChargedHadronBuilderPlugin.

Definition at line 160 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

                                   {
       if (verbosity_) {
         edm::LogPrint("TauChHadronFromPF") << "<PFRecoTauChargedHadronFromPFCandidatePlugin::operator()>:";
         edm::LogPrint("TauChHadronFromPF") << " pluginName = " << name();
       }
  
       ChargedHadronVector output;
  
       // Get the candidates passing our quality cuts
       qcuts_->setPV(vertexAssociator_.associatedVertex(jet));
       CandPtrs candsVector = qcuts_->filterCandRefs(pfCandidates(jet, inputParticleIds_));
  
       for (CandPtrs::iterator cand = candsVector.begin(); cand != candsVector.end(); ++cand) {
         if (verbosity_) {
           edm::LogPrint("TauChHadronFromPF")
               << "processing PFCandidate: Pt = " << (*cand)->pt() << ", eta = " << (*cand)->eta()
               << ", phi = " << (*cand)->phi() << " (pdgId = " << (*cand)->pdgId() << ", charge = " << (*cand)->charge()
               << ")";
         }
  
         PFRecoTauChargedHadron::PFRecoTauChargedHadronAlgorithm algo = PFRecoTauChargedHadron::kUndefined;
         if (std::abs((*cand)->charge()) > 0.5)
           algo = PFRecoTauChargedHadron::kChargedPFCandidate;
         else
           algo = PFRecoTauChargedHadron::kPFNeutralHadron;
         std::unique_ptr<PFRecoTauChargedHadron> chargedHadron(new PFRecoTauChargedHadron(**cand, algo));
  
         const reco::PFCandidate* pfCand = dynamic_cast<const reco::PFCandidate*>(&**cand);
         if (pfCand) {
           if (pfCand->trackRef().isNonnull())
             chargedHadron->track_ = edm::refToPtr(pfCand->trackRef());
           else if (pfCand->muonRef().isNonnull() && pfCand->muonRef()->innerTrack().isNonnull())
             chargedHadron->track_ = edm::refToPtr(pfCand->muonRef()->innerTrack());
           else if (pfCand->muonRef().isNonnull() && pfCand->muonRef()->globalTrack().isNonnull())
             chargedHadron->track_ = edm::refToPtr(pfCand->muonRef()->globalTrack());
           else if (pfCand->muonRef().isNonnull() && pfCand->muonRef()->outerTrack().isNonnull())
             chargedHadron->track_ = edm::refToPtr(pfCand->muonRef()->outerTrack());
           else if (pfCand->gsfTrackRef().isNonnull())
             chargedHadron->track_ = edm::refToPtr(pfCand->gsfTrackRef());
         }  // TauReco@MiniAOD: Tracks only available dynamically, so no possiblity to save ref here; checked by code downstream
  
         chargedHadron->positionAtECALEntrance_ = atECALEntrance(&**cand, bField_);
         chargedHadron->chargedPFCandidate_ = (*cand);
         chargedHadron->addDaughter(*cand);
  
         int pdgId = std::abs(chargedHadron->chargedPFCandidate_->pdgId());
  
         if (chargedHadron->pt() > minMergeChargedHadronPt_) {
           for (const auto& jetConstituent : jet.daughterPtrVector()) {
             // CV: take care of not double-counting energy in case "charged" PFCandidate is in fact a PFNeutralHadron
             if (jetConstituent == chargedHadron->chargedPFCandidate_)
               continue;
  
             int jetConstituentPdgId = std::abs(jetConstituent->pdgId());
             if (!(jetConstituentPdgId == 130 || jetConstituentPdgId == 22))
               continue;
  
             double dR = deltaR(atECALEntrance(jetConstituent.get(), bField_),
                                atECALEntrance(chargedHadron->chargedPFCandidate_.get(), bField_));
             double dRmerge = -1.;
             int minBlockElementMatches = 1000;
             int maxUnmatchedBlockElements = 0;
             double minMergeEt = 1.e+6;
             if (jetConstituentPdgId == 130) {
               if (pdgId == 211)
                 dRmerge = dRmergeNeutralHadronWrtChargedHadron_;
               else if (pdgId == 130)
                 dRmerge = dRmergeNeutralHadronWrtNeutralHadron_;
               else if (pdgId == 11)
                 dRmerge = dRmergeNeutralHadronWrtElectron_;
               else
                 dRmerge = dRmergeNeutralHadronWrtOther_;
               minBlockElementMatches = minBlockElementMatchesNeutralHadron_;
               maxUnmatchedBlockElements = maxUnmatchedBlockElementsNeutralHadron_;
               minMergeEt = minMergeNeutralHadronEt_;
             } else if (jetConstituentPdgId == 22) {
               if (pdgId == 211)
                 dRmerge = dRmergePhotonWrtChargedHadron_;
               else if (pdgId == 130)
                 dRmerge = dRmergePhotonWrtNeutralHadron_;
               else if (pdgId == 11)
                 dRmerge = dRmergePhotonWrtElectron_;
               else
                 dRmerge = dRmergePhotonWrtOther_;
               minBlockElementMatches = minBlockElementMatchesPhoton_;
               maxUnmatchedBlockElements = maxUnmatchedBlockElementsPhoton_;
               minMergeEt = minMergeGammaEt_;
             }
  
             if (jetConstituent->et() > minMergeEt) {
               if (dR < dRmerge) {
                 chargedHadron->neutralPFCandidates_.push_back(jetConstituent);
                 chargedHadron->addDaughter(jetConstituent);
               } else {
                 // TauReco@MiniAOD: No access to PF blocks at MiniAOD level, but the code below seems to have very minor impact
                 const reco::PFCandidate* pfJetConstituent =
                     dynamic_cast<const reco::PFCandidate*>(jetConstituent.get());
                 if (pfCand != nullptr && pfJetConstituent != nullptr) {
                   if (isMatchedByBlockElement(*pfJetConstituent,
                                               *pfCand,
                                               minBlockElementMatches,
                                               minBlockElementMatches,
                                               maxUnmatchedBlockElements)) {
                     chargedHadron->neutralPFCandidates_.push_back(jetConstituent);
                     chargedHadron->addDaughter(jetConstituent);
                   }
                 }
               }
             }
           }
         }
  
         setChargedHadronP4(*chargedHadron);
  
         if (verbosity_) {
           edm::LogPrint("TauChHadronFromPF") << *chargedHadron;
         }
         // Update the vertex
         if (chargedHadron->daughterPtr(0).isNonnull())
           chargedHadron->setVertex(chargedHadron->daughterPtr(0)->vertex());
         output.push_back(std::move(chargedHadron));
       }
  
       return output.release();
     }