#include <TtFullLepHypGenMatch.h>

Inheritance diagram for TtFullLepHypGenMatch:

Public Member Functions
	TtFullLepHypGenMatch (const edm::ParameterSet &)

	~TtFullLepHypGenMatch () override

Public Member Functions inherited from TtFullLepHypothesis
	TtFullLepHypothesis (const edm::ParameterSet &)
	default constructor More...

	~TtFullLepHypothesis () 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

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::vector< ModuleDescription const * > &modules, 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

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 TtFullLepHypothesis
reco::ShallowClonePtrCandidate *	b_

reco::ShallowClonePtrCandidate *	bBar_

edm::EDGetTokenT< std::vector< pat::Electron > >	elecsToken_

bool	getMatch_

std::string	jetCorrectionLevel_

edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_

int	key_
	hypothesis key (to be set by the buildKey function) More...

reco::ShallowClonePtrCandidate *	lepton_

reco::ShallowClonePtrCandidate *	leptonBar_

edm::EDGetTokenT< std::vector< std::vector< int > > >	matchToken_
	input label for all necessary collections More...

edm::EDGetTokenT< std::vector< pat::MET > >	metsToken_

edm::EDGetTokenT< std::vector< pat::Muon > >	musToken_

reco::ShallowClonePtrCandidate *	neutrino_

reco::ShallowClonePtrCandidate *	neutrinoBar_

reco::LeafCandidate *	recNu
	candidates needed for the genmatch hypothesis More...

reco::LeafCandidate *	recNuBar

Private Member Functions
void	buildHypo (edm::Event &evt, const edm::Handle< std::vector< pat::Electron > > &elecs, const edm::Handle< std::vector< pat::Muon > > &mus, const edm::Handle< std::vector< pat::Jet > > &jets, const edm::Handle< std::vector< pat::MET > > &mets, std::vector< int > &match, 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...

void	buildMatchingNeutrinos (edm::Event &, const edm::Handle< std::vector< pat::MET > > &)

template<typename O >
int	findMatchingLepton (const reco::GenParticle *, const edm::Handle< std::vector< O > > &)

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	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from TtFullLepHypothesis
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...

int	key () const
	return key 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...

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>
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<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 7 of file TtFullLepHypGenMatch.h.

Constructor & Destructor Documentation

◆ TtFullLepHypGenMatch()

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

explicit

Definition at line 5 of file TtFullLepHypGenMatch.cc.

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

◆ ~TtFullLepHypGenMatch()

TtFullLepHypGenMatch::~TtFullLepHypGenMatch ( )

override

Definition at line 8 of file TtFullLepHypGenMatch.cc.

8 {}

Member Function Documentation

◆ buildHypo()

void TtFullLepHypGenMatch::buildHypo	(	edm::Event &	evt,
		const edm::Handle< std::vector< pat::Electron > > &	elecs,
		const edm::Handle< std::vector< pat::Muon > > &	mus,
		const edm::Handle< std::vector< pat::Jet > > &	jets,
		const edm::Handle< std::vector< pat::MET > > &	mets,
		std::vector< int > &	match,
		const unsigned int	iComb
	)

overrideprivatevirtual

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

Implements TtFullLepHypothesis.

Definition at line 10 of file TtFullLepHypGenMatch.cc.

                                                                {
   // -----------------------------------------------------
   // add jets
   // -----------------------------------------------------
   for (unsigned idx = 0; idx < match.size(); ++idx) {
     if (isValid(match[idx], jets)) {
       switch (idx) {
         case TtFullLepEvtPartons::B:
           setCandidate(jets, match[idx], b_, jetCorrectionLevel_);
           break;
         case TtFullLepEvtPartons::BBar:
           setCandidate(jets, match[idx], bBar_, jetCorrectionLevel_);
           break;
       }
     }
   }
  
   // -----------------------------------------------------
   // add leptons
   // -----------------------------------------------------
   // get genEvent
   edm::Handle<TtGenEvent> genEvt;
   evt.getByToken(genEvtToken_, genEvt);
  
   // push back fake indices if no leptons in genevent
   if (!genEvt->isFullLeptonic() || !genEvt->lepton() || !genEvt->leptonBar()) {
     match.push_back(-1);
     match.push_back(-1);
     match.push_back(-1);
     match.push_back(-1);
   } else if (genEvt->isFullLeptonic(WDecay::kElec, WDecay::kElec) && elecs->size() >= 2) {
     //search indices for electrons
     int iLepBar = findMatchingLepton(genEvt->leptonBar(), elecs);
     setCandidate(elecs, iLepBar, leptonBar_);
     match.push_back(iLepBar);
     int iLep = findMatchingLepton(genEvt->lepton(), elecs);
     setCandidate(elecs, iLep, lepton_);
     match.push_back(iLep);
  
     // fake indices for muons
     match.push_back(-1);
     match.push_back(-1);
   } else if (genEvt->isFullLeptonic(WDecay::kElec, WDecay::kMuon) && !elecs->empty() && !mus->empty()) {
     if (genEvt->leptonBar()->isElectron()) {
       // push back index for e+
       int iLepBar = findMatchingLepton(genEvt->leptonBar(), elecs);
       setCandidate(elecs, iLepBar, leptonBar_);
       match.push_back(iLepBar);
       // push back fake indices for e- and mu+
       match.push_back(-1);
       match.push_back(-1);
       // push back index for mu-
       int iLep = findMatchingLepton(genEvt->lepton(), mus);
       setCandidate(mus, iLep, lepton_);
       match.push_back(iLep);
     } else {
       // push back fake index for e+
       match.push_back(-1);
       // push back index for e-
       int iLepBar = findMatchingLepton(genEvt->leptonBar(), mus);
       setCandidate(mus, iLepBar, leptonBar_);
       match.push_back(iLepBar);
       // push back index for mu+
       int iLep = findMatchingLepton(genEvt->lepton(), elecs);
       setCandidate(elecs, iLep, lepton_);
       match.push_back(iLep);
       // push back fake index for mu-
       match.push_back(-1);
     }
   } else if (genEvt->isFullLeptonic(WDecay::kMuon, WDecay::kMuon) && mus->size() >= 2) {
     // fake indices for electrons
     match.push_back(-1);
     match.push_back(-1);
  
     //search indices for electrons
     int iLepBar = findMatchingLepton(genEvt->leptonBar(), mus);
     setCandidate(mus, iLepBar, leptonBar_);
     match.push_back(iLepBar);
     int iLep = findMatchingLepton(genEvt->lepton(), mus);
     setCandidate(mus, iLep, lepton_);
     match.push_back(iLep);
   } else {  //this 'else' should happen if at least one genlepton is a tau
     match.push_back(-1);
     match.push_back(-1);
     match.push_back(-1);
     match.push_back(-1);
   }
  
   // -----------------------------------------------------
   // add met and neutrinos
   // -----------------------------------------------------
   if (!mets->empty()) {
     //setCandidate(mets, 0, met_);
     buildMatchingNeutrinos(evt, mets);
   }
 }

References TtFullLepEvtPartons::B, TtFullLepHypothesis::b_, TtFullLepEvtPartons::BBar, TtFullLepHypothesis::bBar_, buildMatchingNeutrinos(), singleTopDQM_cfi::elecs, findMatchingLepton(), TtGenEvtProducer_cfi::genEvt, genEvtToken_, edm::Event::getByToken(), training_settings::idx, TtFullLepHypothesis::isValid(), TtFullLepHypothesis::jetCorrectionLevel_, singleTopDQM_cfi::jets, WDecay::kElec, WDecay::kMuon, TtFullLepHypothesis::lepton_, TtFullLepHypothesis::leptonBar_, match(), singleTopDQM_cfi::mets, and TtFullLepHypothesis::setCandidate().

◆ buildKey()

void TtFullLepHypGenMatch::buildKey ( )

inlineoverrideprivatevirtual

build the event hypothesis key

Implements TtFullLepHypothesis.

Definition at line 14 of file TtFullLepHypGenMatch.h.

14 { key_ = TtFullLeptonicEvent::kGenMatch; };

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

◆ buildMatchingNeutrinos()

void TtFullLepHypGenMatch::buildMatchingNeutrinos	(	edm::Event &	evt,
		const edm::Handle< std::vector< pat::MET > > &	mets
	)

private

Definition at line 128 of file TtFullLepHypGenMatch.cc.

                                                                                                           {
   // get genEvent
   edm::Handle<TtGenEvent> genEvt;
   evt.getByToken(genEvtToken_, genEvt);
  
   if (genEvt->isTtBar() && genEvt->isFullLeptonic() && genEvt->neutrino() && genEvt->neutrinoBar()) {
     double momXNu = genEvt->neutrino()->px();
     double momYNu = genEvt->neutrino()->py();
     double momXNuBar = genEvt->neutrinoBar()->px();
     double momYNuBar = genEvt->neutrinoBar()->py();
  
     double momXMet = mets->at(0).px();
     double momYMet = mets->at(0).py();
  
     double momXNeutrino = 0.5 * (momXNu - momXNuBar + momXMet);
     double momYNeutrino = 0.5 * (momYNu - momYNuBar + momYMet);
     double momXNeutrinoBar = momXMet - momXNeutrino;
     double momYNeutrinoBar = momYMet - momYNeutrino;
  
     math::XYZTLorentzVector recNuFM(
         momXNeutrino, momYNeutrino, 0, sqrt(momXNeutrino * momXNeutrino + momYNeutrino * momYNeutrino));
     recNu = new reco::LeafCandidate(0, recNuFM);
  
     math::XYZTLorentzVector recNuBarFM(momXNeutrinoBar,
                                        momYNeutrinoBar,
                                        0,
                                        sqrt(momXNeutrinoBar * momXNeutrinoBar + momYNeutrinoBar * momYNeutrinoBar));
     recNuBar = new reco::LeafCandidate(0, recNuBarFM);
   }
 }

References TtGenEvtProducer_cfi::genEvt, genEvtToken_, edm::Event::getByToken(), singleTopDQM_cfi::mets, TtFullLepHypothesis::recNu, TtFullLepHypothesis::recNuBar, and mathSSE::sqrt().

Referenced by buildHypo().

◆ findMatchingLepton()

template<typename O >

int TtFullLepHypGenMatch::findMatchingLepton	(	const reco::GenParticle *	genLep,
		const edm::Handle< std::vector< O > > &	leps
	)

private

Definition at line 114 of file TtFullLepHypGenMatch.cc.

                                                                                    {
   int idx = -1;
   double minDR = -1;
   for (unsigned i = 0; i < leps->size(); ++i) {
     double dR = deltaR(genLep->eta(), genLep->phi(), (*leps)[i].eta(), (*leps)[i].phi());
     if (minDR < 0 || dR < minDR) {
       minDR = dR;
       idx = i;
     }
   }
   return idx;
 }

References PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, reco::LeafCandidate::eta(), mps_fire::i, training_settings::idx, TtSemiLepHitFitProducer_Electrons_cfi::leps, PDWG_EXOHSCP_cff::minDR, and reco::LeafCandidate::phi().

Referenced by buildHypo().

Member Data Documentation

◆ genEvtToken_

edm::EDGetTokenT<TtGenEvent> TtFullLepHypGenMatch::genEvtToken_

protected

Definition at line 29 of file TtFullLepHypGenMatch.h.

Referenced by buildHypo(), and buildMatchingNeutrinos().

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members