#include <TtFullLepHypGenMatch.h>

Inheritance diagram for TtFullLepHypGenMatch:

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

	~TtFullLepHypGenMatch ()

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

	~TtFullLepHypothesis ()
	default destructor More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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
virtual 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)
	build event hypothesis from the reco objects of a semi-leptonic event More...

virtual void	buildKey ()
	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
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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 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...

virtual void	produce (edm::Event &, const edm::EventSetup &)
	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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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

Constructor & Destructor Documentation

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

explicit

Definition at line 5 of file TtFullLepHypGenMatch.cc.

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

TtFullLepHypGenMatch::~TtFullLepHypGenMatch ( )

Definition at line 11 of file TtFullLepHypGenMatch.cc.

11 { }

Member Function Documentation

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
	)

privatevirtual

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

Implements TtFullLepHypothesis.

Definition at line 14 of file TtFullLepHypGenMatch.cc.

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

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

virtual void TtFullLepHypGenMatch::buildKey ( )

inlineprivatevirtual

build the event hypothesis key

Implements TtFullLepHypothesis.

Definition at line 17 of file TtFullLepHypGenMatch.h.

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

17 { key_= TtFullLeptonicEvent::kGenMatch; };

TtFullLepHypothesis::key_

int key_

hypothesis key (to be set by the buildKey function)

Definition: TtFullLepHypothesis.h:89

TtEvent::kGenMatch

Definition: TtEvent.h:31

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

private

Definition at line 138 of file TtFullLepHypGenMatch.cc.

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

Referenced by buildHypo().

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

template<typename O >

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

private

Definition at line 123 of file TtFullLepHypGenMatch.cc.

References deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, reco::LeafCandidate::eta(), i, customizeTrackingMonitorSeedNumber::idx, HLT_25ns14e33_v1_cff::minDR, and reco::LeafCandidate::phi().

Referenced by buildHypo().

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

Member Data Documentation

edm::EDGetTokenT<TtGenEvent> TtFullLepHypGenMatch::genEvtToken_

protected

Definition at line 35 of file TtFullLepHypGenMatch.h.

Referenced by buildHypo(), and buildMatchingNeutrinos().

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation