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

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

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)

virtual	~ProducerBase () noexcept(false)

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

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)

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, 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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 7 of file TtSemiLepHypGenMatch.h.

Constructor & Destructor Documentation

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

override

Definition at line 11 of file TtSemiLepHypGenMatch.cc.

11 { }

Member Function Documentation

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 14 of file TtSemiLepHypGenMatch.cc.

Referenced by buildKey().

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

void TtSemiLepHypGenMatch::buildKey ( )

inlineoverrideprivatevirtual

build the event hypothesis key

Implements TtSemiLepHypothesis.

Definition at line 17 of file TtSemiLepHypGenMatch.h.

References buildHypo(), findMatchingLepton(), TtGenEvtProducer_cfi::genEvt, TtSemiLepHypothesis::key_, and TtEvent::kGenMatch.

17 { key_= TtSemiLeptonicEvent::kGenMatch; };

TtEvent::kGenMatch

Definition: TtEvent.h:31

TtSemiLepHypothesis::key_

int key_

hypothesis key (to be set by the buildKey function)

Definition: TtSemiLepHypothesis.h:97

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 75 of file TtSemiLepHypGenMatch.cc.

References deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::LeafCandidate::eta(), mps_fire::i, TtGenEvent::isSemiLeptonic(), TtGenEvent::isTtBar(), TtSemiLepHitFitProducer_Electrons_cfi::leps, TtSemiLepHypothesis::leptonType(), electronTrackIsolations_cfi::minDR, reco::LeafCandidate::phi(), and TtGenEvent::singleLepton().

Referenced by buildHypo(), and buildKey().

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

Member Data Documentation

edm::EDGetTokenT<TtGenEvent> TtSemiLepHypGenMatch::genEvtToken_

protected

Definition at line 30 of file TtSemiLepHypGenMatch.h.

Referenced by buildHypo().

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation