reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin Class Reference

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
const edm::Event *	evt () const
edm::Event *	evt ()
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_
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

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

private

Definition at line 53 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 82 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References dRmergeNeutralHadronWrtChargedHadron_, dRmergeNeutralHadronWrtElectron_, dRmergeNeutralHadronWrtNeutralHadron_, dRmergeNeutralHadronWrtOther_, dRmergePhotonWrtChargedHadron_, dRmergePhotonWrtElectron_, dRmergePhotonWrtNeutralHadron_, dRmergePhotonWrtOther_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), inputParticleIds_, maxUnmatchedBlockElementsNeutralHadron_, maxUnmatchedBlockElementsPhoton_, minBlockElementMatchesNeutralHadron_, minBlockElementMatchesPhoton_, minMergeChargedHadronPt_, minMergeGammaEt_, minMergeNeutralHadronEt_, qcuts_, and verbosity_.

        : PFRecoTauChargedHadronBuilderPlugin(pset, std::move(iC)),
          vertexAssociator_(pset.getParameter<edm::ParameterSet>("qualityCuts"), std::move(iC)),
          qcuts_(nullptr) {
      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");
    }

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

override

Definition at line 111 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References qcuts_.

{ delete qcuts_; }

Member Function Documentation

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

overridevirtual

Hook called at the beginning of the event.

Reimplemented from reco::tau::PFRecoTauChargedHadronBuilderPlugin.

Definition at line 114 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References bField_, reco::PFCandidate::elementsInBlocks(), reco::tau::RecoTauEventHolderPlugin::evt(), reco::tau::RecoTauEventHolderPlugin::evtSetup(), edm::EventSetup::get(), MagneticField::inTesla(), trackerHitRTTI::isMatched(), reco::tau::RecoTauVertexAssociator::setEvent(), vertexAssociator_, and PV3DBase< T, PVType, FrameType >::z().

                                                                 {
      vertexAssociator_.setEvent(*evt());

      edm::ESHandle<MagneticField> pSetup;
      evtSetup()->get<IdealMagneticFieldRecord>().get(pSetup);
      bField_ = pSetup->inTesla(GlobalPoint(0, 0, 0)).z();
    }

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.

References funct::abs(), reco::CompositeCandidate::addDaughter(), ecalcalib_dqm_sourceclient-live_cfg::algo, reco::tau::RecoTauVertexAssociator::associatedVertex(), reco::tau::atECALEntrance(), bField_, muons2muons_cfi::chargedHadron, reco::CompositePtrCandidate::daughterPtrVector(), DEFINE_EDM_PLUGIN, reco::deltaR(), HGC3DClusterGenMatchSelector_cfi::dR, dRmergeNeutralHadronWrtChargedHadron_, dRmergeNeutralHadronWrtElectron_, dRmergeNeutralHadronWrtNeutralHadron_, dRmergeNeutralHadronWrtOther_, dRmergePhotonWrtChargedHadron_, dRmergePhotonWrtElectron_, dRmergePhotonWrtNeutralHadron_, dRmergePhotonWrtOther_, reco::tau::RecoTauQualityCuts::filterCandRefs(), reco::PFCandidate::gsfTrackRef(), inputParticleIds_, edm::Ref< C, T, F >::isNonnull(), reco::PFRecoTauChargedHadron::kChargedPFCandidate, reco::PFRecoTauChargedHadron::kPFNeutralHadron, reco::PFRecoTauChargedHadron::kUndefined, maxUnmatchedBlockElementsNeutralHadron_, maxUnmatchedBlockElementsPhoton_, minBlockElementMatchesNeutralHadron_, minBlockElementMatchesPhoton_, minMergeChargedHadronPt_, minMergeGammaEt_, minMergeNeutralHadronEt_, eostools::move(), reco::PFCandidate::muonRef(), reco::tau::RecoTauNamedPlugin::name(), convertSQLitetoXML_cfg::output, EgammaValidation_cff::pdgId, reco::tau::pfCandidates(), qcuts_, edm::refToPtr(), reco::tau::setChargedHadronP4(), reco::tau::RecoTauQualityCuts::setPV(), reco::PFCandidate::trackRef(), verbosity_, and vertexAssociator_.

                                  {
      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();
    }

Member Data Documentation

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::bField_

private

Definition at line 77 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by beginEvent(), and operator()().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtChargedHadron_

private

Definition at line 61 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtElectron_

private

Definition at line 63 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtNeutralHadron_

private

Definition at line 62 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtOther_

private

Definition at line 64 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtChargedHadron_

private

Definition at line 67 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtElectron_

private

Definition at line 69 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtNeutralHadron_

private

Definition at line 68 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtOther_

private

Definition at line 70 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

std::vector<int> reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::inputParticleIds_

private

Definition at line 59 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsNeutralHadron_

private

Definition at line 66 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsPhoton_

private

Definition at line 72 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesNeutralHadron_

private

Definition at line 65 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesPhoton_

private

Definition at line 71 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeChargedHadronPt_

private

Definition at line 75 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeGammaEt_

private

Definition at line 74 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeNeutralHadronEt_

private

Definition at line 73 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

RecoTauQualityCuts* reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::qcuts_

private

Definition at line 57 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), PFRecoTauChargedHadronFromPFCandidatePlugin(), and ~PFRecoTauChargedHadronFromPFCandidatePlugin().

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::verbosity_

private

Definition at line 79 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by operator()(), and PFRecoTauChargedHadronFromPFCandidatePlugin().

RecoTauVertexAssociator reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::vertexAssociator_

private

Definition at line 55 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Referenced by beginEvent(), and operator()().