TtSemiLepHypGenMatch Class Reference

#include <TtSemiLepHypGenMatch.h>

Inheritance diagram for TtSemiLepHypGenMatch:

Public Member Functions
	TtSemiLepHypGenMatch (const edm::ParameterSet &)
	~TtSemiLepHypGenMatch () override
Public Member Functions inherited from TtSemiLepHypothesis
	TtSemiLepHypothesis (const edm::ParameterSet &)
	default constructor More...
	~TtSemiLepHypothesis () override
	default destructor More...
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDProducer () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
std::vector< bool > const &	recordProvenanceList () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Protected Attributes
edm::EDGetTokenT< TtGenEvent >	genEvtToken_
Protected Attributes inherited from TtSemiLepHypothesis
bool	getMatch_
reco::ShallowClonePtrCandidate *	hadronicB_
std::string	jetCorrectionLevel_
edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_
	input label for all necessary collections More...
int	key_
	hypothesis key (to be set by the buildKey function) More...
edm::EDGetTokenT< edm::View< reco::RecoCandidate > >	lepsToken_
reco::ShallowClonePtrCandidate *	lepton_
reco::ShallowClonePtrCandidate *	leptonicB_
reco::ShallowClonePtrCandidate *	lightQ_
reco::ShallowClonePtrCandidate *	lightQBar_
edm::EDGetTokenT< std::vector< std::vector< int > > >	matchToken_
edm::EDGetTokenT< std::vector< pat::MET > >	metsToken_
reco::ShallowClonePtrCandidate *	neutrino_
int	neutrinoSolutionType_
	algorithm used to calculate neutrino solutions (see cfi for further details) More...
edm::EDGetTokenT< int >	nJetsConsideredToken_
int	numberOfRealNeutrinoSolutions_

Private Member Functions
void	buildHypo (edm::Event &, const edm::Handle< edm::View< reco::RecoCandidate > > &, const edm::Handle< std::vector< pat::MET > > &, const edm::Handle< std::vector< pat::Jet > > &, std::vector< int > &, const unsigned int iComb) override
	build event hypothesis from the reco objects of a semi-leptonic event More...
void	buildKey () override
	build the event hypothesis key More...
int	findMatchingLepton (const edm::Handle< TtGenEvent > &genEvt, const edm::Handle< edm::View< reco::RecoCandidate > > &)
	find index of the candidate nearest to the singleLepton of the generator event in the collection; return -1 if this fails More...

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex > >
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from TtSemiLepHypothesis
void	buildHypo (const edm::Handle< edm::View< reco::RecoCandidate > > &leps, const edm::Handle< std::vector< pat::MET > > &mets, const edm::Handle< std::vector< pat::Jet > > &jets, std::vector< int > &jetPartonAssociation)
	minimalistic build function for simple hypotheses More...
reco::CompositeCandidate	hypo ()
	return event hypothesis More...
bool	isValid (const int &idx, const edm::Handle< std::vector< pat::Jet > > &jets)
	check if index is in valid range of selected jets More...
std::string	jetCorrectionLevel (const std::string &quarkType)
	helper function to construct the proper correction level string for corresponding quarkType More...
int	key () const
WDecay::LeptonType	leptonType (const reco::RecoCandidate *cand)
	determine lepton type of reco candidate and return a corresponding WDecay::LeptonType; the type is kNone if it is whether a muon nor an electron More...
void	produce (edm::Event &, const edm::EventSetup &) override
	produce the event hypothesis as CompositeCandidate and Key More...
void	resetCandidates ()
	reset candidate pointers before hypo build process More...
template<typename C >
void	setCandidate (const edm::Handle< C > &handle, const int &idx, reco::ShallowClonePtrCandidate *&clone)
	use one object in a collection to set a ShallowClonePtrCandidate More...
void	setCandidate (const edm::Handle< std::vector< pat::Jet > > &handle, const int &idx, reco::ShallowClonePtrCandidate *&clone, const std::string &correctionLevel)
	use one object in a jet collection to set a ShallowClonePtrCandidate with proper jet corrections More...
void	setNeutrino (const edm::Handle< std::vector< pat::MET > > &met, const edm::Handle< edm::View< reco::RecoCandidate > > &leps, const int &idx, const int &type)
	set neutrino, using mW = 80.4 to calculate the neutrino pz More...
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 6 of file TtSemiLepHypGenMatch.h.

Constructor & Destructor Documentation

TtSemiLepHypGenMatch()

TtSemiLepHypGenMatch::TtSemiLepHypGenMatch ( const edm::ParameterSet &  cfg )

explicit

Definition at line 6 of file TtSemiLepHypGenMatch.cc.

    : TtSemiLepHypothesis(cfg), genEvtToken_(consumes<TtGenEvent>(edm::InputTag("genEvt"))) {}

~TtSemiLepHypGenMatch()

TtSemiLepHypGenMatch::~TtSemiLepHypGenMatch ( )

override

Definition at line 9 of file TtSemiLepHypGenMatch.cc.

{}

Member Function Documentation

buildHypo()

void TtSemiLepHypGenMatch::buildHypo ( edm::Event &  evt, const edm::Handle< edm::View< reco::RecoCandidate > > &  leps, const edm::Handle< std::vector< pat::MET > > &  mets, const edm::Handle< std::vector< pat::Jet > > &  jets, std::vector< int > &  match, const unsigned int  iComb  )

overrideprivatevirtual

build event hypothesis from the reco objects of a semi-leptonic event

Implements TtSemiLepHypothesis.

Definition at line 11 of file TtSemiLepHypGenMatch.cc.

                                                               {
  // -----------------------------------------------------
  // get genEvent (to distinguish between uds and c quarks
  // and for the lepton matching)
  // -----------------------------------------------------
  edm::Handle<TtGenEvent> genEvt;
  evt.getByToken(genEvtToken_, genEvt);
 
  // -----------------------------------------------------
  // add jets
  // -----------------------------------------------------
  for (unsigned idx = 0; idx < match.size(); ++idx) {
    if (isValid(match[idx], jets)) {
      switch (idx) {
        case TtSemiLepEvtPartons::LightQ:
          if (std::abs(genEvt->hadronicDecayQuark()->pdgId()) == 4)
            setCandidate(jets, match[idx], lightQ_, jetCorrectionLevel("cQuark"));
          else
            setCandidate(jets, match[idx], lightQ_, jetCorrectionLevel("udsQuark"));
          break;
        case TtSemiLepEvtPartons::LightQBar:
          if (std::abs(genEvt->hadronicDecayQuarkBar()->pdgId()) == 4)
            setCandidate(jets, match[idx], lightQBar_, jetCorrectionLevel("cQuark"));
          else
            setCandidate(jets, match[idx], lightQBar_, jetCorrectionLevel("udsQuark"));
          break;
        case TtSemiLepEvtPartons::HadB:
          setCandidate(jets, match[idx], hadronicB_, jetCorrectionLevel("bQuark"));
          break;
        case TtSemiLepEvtPartons::LepB:
          setCandidate(jets, match[idx], leptonicB_, jetCorrectionLevel("bQuark"));
          break;
      }
    }
  }
 
  // -----------------------------------------------------
  // add lepton
  // -----------------------------------------------------
  int iLepton = findMatchingLepton(genEvt, leps);
  if (iLepton < 0)
    return;
  setCandidate(leps, iLepton, lepton_);
  match.push_back(iLepton);
 
  // -----------------------------------------------------
  // add neutrino
  // -----------------------------------------------------
  if (mets->empty())
    return;
  if (neutrinoSolutionType_ == -1)
    setCandidate(mets, 0, neutrino_);
  else
    setNeutrino(mets, leps, iLepton, neutrinoSolutionType_);
}

References funct::abs(), findMatchingLepton(), TtGenEvtProducer_cfi::genEvt, genEvtToken_, edm::Event::getByToken(), TtSemiLepEvtPartons::HadB, TtSemiLepHypothesis::hadronicB_, heavyIonCSV_trainingSettings::idx, TtSemiLepHypothesis::isValid(), TtSemiLepHypothesis::jetCorrectionLevel(), singleTopDQM_cfi::jets, TtSemiLepEvtPartons::LepB, TtSemiLepHitFitProducer_Electrons_cfi::leps, TtSemiLepHypothesis::lepton_, TtSemiLepHypothesis::leptonicB_, TtSemiLepEvtPartons::LightQ, TtSemiLepHypothesis::lightQ_, TtSemiLepEvtPartons::LightQBar, TtSemiLepHypothesis::lightQBar_, match(), singleTopDQM_cfi::mets, TtSemiLepHypothesis::neutrino_, TtSemiLepHypothesis::neutrinoSolutionType_, TtSemiLepHypothesis::setCandidate(), and TtSemiLepHypothesis::setNeutrino().

buildKey()

void TtSemiLepHypGenMatch::buildKey ( )

inlineoverrideprivatevirtual

build the event hypothesis key

Implements TtSemiLepHypothesis.

Definition at line 13 of file TtSemiLepHypGenMatch.h.

{ key_ = TtSemiLeptonicEvent::kGenMatch; };

References TtSemiLepHypothesis::key_, and TtEvent::kGenMatch.

findMatchingLepton()

int TtSemiLepHypGenMatch::findMatchingLepton ( const edm::Handle< TtGenEvent > &  genEvt, const edm::Handle< edm::View< reco::RecoCandidate > > &  leps  )

private

find index of the candidate nearest to the singleLepton of the generator event in the collection; return -1 if this fails

Definition at line 73 of file TtSemiLepHypGenMatch.cc.

                                                                                                 {
  int genIdx = -1;
 
  // jump out with -1 when the collection is empty
  if (leps->empty())
    return genIdx;
 
  if (genEvt->isTtBar() && genEvt->isSemiLeptonic(leptonType(&(leps->front()))) && genEvt->singleLepton()) {
    double minDR = -1;
    for (unsigned i = 0; i < leps->size(); ++i) {
      double dR =
          deltaR(genEvt->singleLepton()->eta(), genEvt->singleLepton()->phi(), (*leps)[i].eta(), (*leps)[i].phi());
      if (minDR < 0 || dR < minDR) {
        minDR = dR;
        genIdx = i;
      }
    }
  }
  return genIdx;
}

References PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, TtGenEvtProducer_cfi::genEvt, mps_fire::i, TtSemiLepHitFitProducer_Electrons_cfi::leps, TtSemiLepHypothesis::leptonType(), and HLT_FULL_cff::minDR.

Referenced by buildHypo().

Member Data Documentation

genEvtToken_

edm::EDGetTokenT<TtGenEvent> TtSemiLepHypGenMatch::genEvtToken_

protected

Definition at line 25 of file TtSemiLepHypGenMatch.h.

Referenced by buildHypo().