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::auto_ptr< ChargedHadronVector >	return_type

Detailed Description

Definition at line 46 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Member Typedef Documentation

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

private

Definition at line 57 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 86 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 115 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 121 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 163 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References funct::abs(), reco::CompositeCandidate::addDaughter(), patPFMETCorrections_cff::algo, reco::tau::RecoTauVertexAssociator::associatedVertex(), reco::tau::atECALEntrance(), bField_, muons2muons_cfi::chargedHadron, reco::CompositePtrCandidate::daughterPtrVector(), DEFINE_EDM_PLUGIN, reco::deltaR(), 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_, reco::PFCandidate::muonRef(), reco::tau::RecoTauNamedPlugin::name(), convertSQLitetoXML_cfg::output, common_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::auto_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(chargedHadron);
  }

  return output.release();
}

Member Data Documentation

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::bField_

private

Definition at line 81 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtChargedHadron_

private

Definition at line 65 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtElectron_

private

Definition at line 67 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtNeutralHadron_

private

Definition at line 66 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtOther_

private

Definition at line 68 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtChargedHadron_

private

Definition at line 71 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtElectron_

private

Definition at line 73 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtNeutralHadron_

private

Definition at line 72 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtOther_

private

Definition at line 74 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

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

private

Definition at line 63 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsNeutralHadron_

private

Definition at line 70 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsPhoton_

private

Definition at line 76 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesNeutralHadron_

private

Definition at line 69 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesPhoton_

private

Definition at line 75 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeChargedHadronPt_

private

Definition at line 79 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeGammaEt_

private

Definition at line 78 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeNeutralHadronEt_

private

Definition at line 77 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

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

private

Definition at line 61 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::verbosity_

private

Definition at line 83 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

RecoTauVertexAssociator reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::vertexAssociator_

private

Definition at line 59 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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