reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin Class Reference

Inheritance diagram for reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin:

Public Member Functions
virtual void	beginEvent ()
	Hook called at the beginning of the event. More...
return_type	operator() (const reco::PFJet &) const
	Build a collection of chargedHadrons from objects in the input jet. More...
	PFRecoTauChargedHadronFromPFCandidatePlugin (const edm::ParameterSet &, edm::ConsumesCollector &&iC)
virtual	~PFRecoTauChargedHadronFromPFCandidatePlugin ()
Public Member Functions inherited from reco::tau::PFRecoTauChargedHadronBuilderPlugin
	PFRecoTauChargedHadronBuilderPlugin (const edm::ParameterSet &pset, edm::ConsumesCollector &&iC)
virtual	~PFRecoTauChargedHadronBuilderPlugin ()
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 &)
virtual	~RecoTauEventHolderPlugin ()
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::PFCandidatePtr >	PFCandPtrs

Private Attributes
double	dRmergeNeutralHadronWrtChargedHadron_
double	dRmergeNeutralHadronWrtElectron_
double	dRmergeNeutralHadronWrtNeutralHadron_
double	dRmergeNeutralHadronWrtOther_
double	dRmergePhotonWrtChargedHadron_
double	dRmergePhotonWrtElectron_
double	dRmergePhotonWrtNeutralHadron_
double	dRmergePhotonWrtOther_
std::vector< int >	inputPdgIds_
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 41 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Member Typedef Documentation

typedef std::vector<reco::PFCandidatePtr> reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::PFCandPtrs

private

Definition at line 52 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 79 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References dRmergeNeutralHadronWrtChargedHadron_, dRmergeNeutralHadronWrtElectron_, dRmergeNeutralHadronWrtNeutralHadron_, dRmergeNeutralHadronWrtOther_, dRmergePhotonWrtChargedHadron_, dRmergePhotonWrtElectron_, dRmergePhotonWrtNeutralHadron_, dRmergePhotonWrtOther_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), inputPdgIds_, 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_(0)
{
  edm::ParameterSet qcuts_pset = pset.getParameterSet("qualityCuts").getParameterSet("signalQualityCuts");
  qcuts_ = new RecoTauQualityCuts(qcuts_pset);

  inputPdgIds_ = 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.exists("verbosity") ) ?
    pset.getParameter<int>("verbosity") : 0;
}

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

virtual

Definition at line 109 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References qcuts_.

{
  delete qcuts_;
}

Member Function Documentation

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

virtual

Hook called at the beginning of the event.

Reimplemented from reco::tau::PFRecoTauChargedHadronBuilderPlugin.

Definition at line 115 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References reco::tau::RecoTauEventHolderPlugin::evt(), reco::tau::RecoTauVertexAssociator::setEvent(), and vertexAssociator_.

{
  vertexAssociator_.setEvent(*evt());
}

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

virtual

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

Implements reco::tau::PFRecoTauChargedHadronBuilderPlugin.

Definition at line 166 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

References reco::tau::RecoTauVertexAssociator::associatedVertex(), gather_cfg::cout, reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRmergeNeutralHadronWrtChargedHadron_, dRmergeNeutralHadronWrtElectron_, dRmergeNeutralHadronWrtNeutralHadron_, dRmergeNeutralHadronWrtOther_, dRmergePhotonWrtChargedHadron_, dRmergePhotonWrtElectron_, dRmergePhotonWrtNeutralHadron_, dRmergePhotonWrtOther_, reco::PFCandidate::e, reco::tau::RecoTauQualityCuts::filterCandRefs(), reco::PFCandidate::gamma, reco::PFJet::getPFConstituents(), reco::PFCandidate::h, reco::PFCandidate::h0, inputPdgIds_, reco::PFRecoTauChargedHadron::kChargedPFCandidate, reco::PFRecoTauChargedHadron::kPFNeutralHadron, reco::PFRecoTauChargedHadron::kUndefined, maxUnmatchedBlockElementsNeutralHadron_, maxUnmatchedBlockElementsPhoton_, minBlockElementMatchesNeutralHadron_, minBlockElementMatchesPhoton_, minMergeChargedHadronPt_, minMergeGammaEt_, minMergeNeutralHadronEt_, reco::tau::RecoTauNamedPlugin::name(), convertSQLitetoXML_cfg::output, reco::tau::pfCandidates(), qcuts_, edm::refToPtr(), reco::tau::setChargedHadronP4(), reco::tau::RecoTauQualityCuts::setPV(), verbosity_, and vertexAssociator_.

{
  if ( verbosity_ ) {
    std::cout << "<PFRecoTauChargedHadronFromPFCandidatePlugin::operator()>:" << std::endl;
    std::cout << " pluginName = " << name() << std::endl;
  }

  ChargedHadronVector output;

  // Get the candidates passing our quality cuts
  qcuts_->setPV(vertexAssociator_.associatedVertex(jet));
  PFCandPtrs candsVector = qcuts_->filterCandRefs(pfCandidates(jet, inputPdgIds_));

  for ( PFCandPtrs::iterator cand = candsVector.begin();
    cand != candsVector.end(); ++cand ) {
    if ( verbosity_ ) {
      std::cout << "processing PFCandidate: Pt = " << (*cand)->pt() << ", eta = " << (*cand)->eta() << ", phi = " << (*cand)->phi() 
        << " (type = " << getPFCandidateType((*cand)->particleId()) << ", charge = " << (*cand)->charge() << ")" << std::endl;
    }
    
    PFRecoTauChargedHadron::PFRecoTauChargedHadronAlgorithm algo = PFRecoTauChargedHadron::kUndefined;
    if ( fabs((*cand)->charge()) > 0.5 ) algo = PFRecoTauChargedHadron::kChargedPFCandidate;
    else algo = PFRecoTauChargedHadron::kPFNeutralHadron;
    std::auto_ptr<PFRecoTauChargedHadron> chargedHadron(new PFRecoTauChargedHadron(**cand, algo));
    if ( (*cand)->trackRef().isNonnull() ) chargedHadron->track_ = edm::refToPtr((*cand)->trackRef());
    else if ( (*cand)->muonRef().isNonnull() && (*cand)->muonRef()->innerTrack().isNonnull()  ) chargedHadron->track_ = edm::refToPtr((*cand)->muonRef()->innerTrack());
    else if ( (*cand)->muonRef().isNonnull() && (*cand)->muonRef()->globalTrack().isNonnull() ) chargedHadron->track_ = edm::refToPtr((*cand)->muonRef()->globalTrack());
    else if ( (*cand)->muonRef().isNonnull() && (*cand)->muonRef()->outerTrack().isNonnull()  ) chargedHadron->track_ = edm::refToPtr((*cand)->muonRef()->outerTrack());
    else if ( (*cand)->gsfTrackRef().isNonnull() ) chargedHadron->track_ = edm::refToPtr((*cand)->gsfTrackRef());
    chargedHadron->chargedPFCandidate_ = (*cand);
    chargedHadron->addDaughter(*cand);

    chargedHadron->positionAtECALEntrance_ = (*cand)->positionAtECALEntrance();

    reco::PFCandidate::ParticleType chargedPFCandidateType = chargedHadron->chargedPFCandidate_->particleId();

    if ( chargedHadron->pt() > minMergeChargedHadronPt_ ) {
      std::vector<reco::PFCandidatePtr> jetConstituents = jet.getPFConstituents();
      for ( std::vector<reco::PFCandidatePtr>::const_iterator jetConstituent = jetConstituents.begin();
        jetConstituent != jetConstituents.end(); ++jetConstituent ) {
    // CV: take care of not double-counting energy in case "charged" PFCandidate is in fact a PFNeutralHadron
    if ( (*jetConstituent) == chargedHadron->chargedPFCandidate_ ) continue;
    
    reco::PFCandidate::ParticleType jetConstituentType = (*jetConstituent)->particleId();
    if ( !(jetConstituentType == reco::PFCandidate::h0 || jetConstituentType == reco::PFCandidate::gamma) ) continue;

    double dR = deltaR((*jetConstituent)->positionAtECALEntrance(), chargedHadron->positionAtECALEntrance_);
    double dRmerge = -1.;      
    int minBlockElementMatches = 1000;
    int maxUnmatchedBlockElements = 0;
    double minMergeEt = 1.e+6;
    if ( jetConstituentType == reco::PFCandidate::h0 ) {
      if      ( chargedPFCandidateType == reco::PFCandidate::h  ) dRmerge = dRmergeNeutralHadronWrtChargedHadron_;
      else if ( chargedPFCandidateType == reco::PFCandidate::h0 ) dRmerge = dRmergeNeutralHadronWrtNeutralHadron_;
      else if ( chargedPFCandidateType == reco::PFCandidate::e  ) dRmerge = dRmergeNeutralHadronWrtElectron_;
      else                                                        dRmerge = dRmergeNeutralHadronWrtOther_;
      minBlockElementMatches = minBlockElementMatchesNeutralHadron_;
      maxUnmatchedBlockElements = maxUnmatchedBlockElementsNeutralHadron_;
      minMergeEt = minMergeNeutralHadronEt_;
    } else if ( jetConstituentType == reco::PFCandidate::gamma ) {
      if      ( chargedPFCandidateType == reco::PFCandidate::h  ) dRmerge = dRmergePhotonWrtChargedHadron_;
      else if ( chargedPFCandidateType == reco::PFCandidate::h0 ) dRmerge = dRmergePhotonWrtNeutralHadron_;
      else if ( chargedPFCandidateType == reco::PFCandidate::e  ) dRmerge = dRmergePhotonWrtElectron_;
      else                                                        dRmerge = dRmergePhotonWrtOther_;
      minBlockElementMatches = minBlockElementMatchesPhoton_;
      maxUnmatchedBlockElements = maxUnmatchedBlockElementsPhoton_;
      minMergeEt = minMergeGammaEt_;
    }
    if ( (*jetConstituent)->et() > minMergeEt && 
         (dR < dRmerge || isMatchedByBlockElement(**jetConstituent, *chargedHadron->chargedPFCandidate_, minBlockElementMatches, minBlockElementMatches, maxUnmatchedBlockElements)) ) {
      chargedHadron->neutralPFCandidates_.push_back(*jetConstituent);
      chargedHadron->addDaughter(*jetConstituent);
    }
      }
    }

    setChargedHadronP4(*chargedHadron);

    if ( verbosity_ ) {
      chargedHadron->print(std::cout);
    }
    // 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::dRmergeNeutralHadronWrtChargedHadron_

private

Definition at line 60 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtElectron_

private

Definition at line 62 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtNeutralHadron_

private

Definition at line 61 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergeNeutralHadronWrtOther_

private

Definition at line 63 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtChargedHadron_

private

Definition at line 66 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtElectron_

private

Definition at line 68 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtNeutralHadron_

private

Definition at line 67 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::dRmergePhotonWrtOther_

private

Definition at line 69 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

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

private

Definition at line 58 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsNeutralHadron_

private

Definition at line 65 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::maxUnmatchedBlockElementsPhoton_

private

Definition at line 71 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesNeutralHadron_

private

Definition at line 64 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minBlockElementMatchesPhoton_

private

Definition at line 70 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeChargedHadronPt_

private

Definition at line 74 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeGammaEt_

private

Definition at line 73 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

double reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::minMergeNeutralHadronEt_

private

Definition at line 72 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

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

private

Definition at line 56 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

int reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::verbosity_

private

Definition at line 76 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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

RecoTauVertexAssociator reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::vertexAssociator_

private

Definition at line 54 of file PFRecoTauChargedHadronFromPFCandidatePlugin.cc.

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