Class derived from the TopGenEvent for ttbar events. More...

#include "AnalysisDataFormats/TopObjects/interface/TtGenEvent.h"

Inheritance diagram for TtGenEvent:

Public Member Functions
bool	fromGluonFusion () const
	check if the tops were produced from a pair of gluons More...

bool	fromQuarkAnnihilation () const
	check if the tops were produced from qqbar More...

std::pair< WDecay::LeptonType, WDecay::LeptonType >	fullLeptonicChannel () const

const reco::GenParticle *	hadronicDecayB (bool excludeTauLeptons=false) const
	get b of hadronic decay branch More...

const reco::GenParticle *	hadronicDecayQuark (bool invertFlavor=false) const
	get light quark of hadronic decay branch More...

const reco::GenParticle *	hadronicDecayQuarkBar () const
	get light anti-quark of hadronic decay branch More...

const reco::GenParticle *	hadronicDecayTop (bool excludeTauLeptons=false) const
	get top of hadronic decay branch More...

std::vector< const reco::GenParticle * >	hadronicDecayTopRadiation (bool excludeTauLeptons=false) const
	gluons as radiated from the hadronicly decaying top quark More...

const reco::GenParticle *	hadronicDecayW (bool excludeTauLeptons=false) const
	get W of hadronic decay branch More...

bool	isFullHadronic (bool excludeTauLeptons=false) const
	check if the event can be classified as full hadronic More...

bool	isFullLeptonic (bool excludeTauLeptons=false) const
	check if the event can be classified as full leptonic More...

bool	isFullLeptonic (WDecay::LeptonType typeA, WDecay::LeptonType typeB) const
	check if the event is full leptonic with the lepton being of typeA or typeB irrelevant of order; all leptons including taus are allowed More...

bool	isSemiLeptonic (bool excludeTauLeptons=false) const
	check if the event can be classified as semi-laptonic More...

bool	isSemiLeptonic (WDecay::LeptonType typeA) const
	check if the event is semi-leptonic with the lepton being of typeA; all leptons including taus are allowed More...

bool	isSemiLeptonic (WDecay::LeptonType typeA, WDecay::LeptonType typeB) const
	check if the event is semi-leptonic with the lepton being of typeA or typeB; all leptons including taus are allowed More...

bool	isTtBar () const
	check if the event can be classified as ttbar More...

const reco::GenParticle *	lepton (bool excludeTauLeptons=false) const
	get lepton for semi-leptonic or full leptonic decays More...

const reco::GenParticle *	leptonBar (bool excludeTauLeptons=false) const
	get anti-lepton for semi-leptonic or full leptonic decays More...

const reco::GenParticle *	leptonicDecayB (bool excludeTauLeptons=false) const
	get b of leptonic decay branch More...

const reco::GenParticle *	leptonicDecayTop (bool excludeTauLeptons=false) const
	get top of leptonic decay branch More...

std::vector< const reco::GenParticle * >	leptonicDecayTopRadiation (bool excludeTauLeptons=false) const
	gluons as radiated from the leptonicly decaying top quark More...

const reco::GenParticle *	leptonicDecayW (bool excludeTauLeptons=false) const
	get W of leptonic decay branch More...

const reco::GenParticle *	neutrino (bool excludeTauLeptons=false) const
	get neutrino for semi-leptonic or full leptonic decays More...

const reco::GenParticle *	neutrinoBar (bool excludeTauLeptons=false) const
	get anti-neutrino for semi-leptonic or full leptonic decays More...

WDecay::LeptonType	semiLeptonicChannel () const
	return decay channel; all leptons including taus are allowed More...

const reco::GenParticle *	singleLepton (bool excludeTauLeptons=false) const
	return single lepton if available; 0 else More...

const reco::GenParticle *	singleNeutrino (bool excludeTauLeptons=false) const
	return single neutrino if available; 0 else More...

const math::XYZTLorentzVector *	topPair () const
	return combined 4-vector of top and topBar More...

	TtGenEvent ()
	empty constructor More...

	TtGenEvent (reco::GenParticleRefProd &decaySubset, reco::GenParticleRefProd &initSubset)
	default constructor from decaySubset and initSubset More...

	~TtGenEvent () override
	default destructor More...

Public Member Functions inherited from TopGenEvent
const reco::GenParticle *	b () const
	return b quark if available; 0 else More...

const reco::GenParticle *	bBar () const
	return anti-b quark if available; 0 else More...

const reco::GenParticle *	candidate (int id, unsigned int parentId=0) const
	get candidate with given pdg id if available; 0 else More...

const reco::GenParticle *	daughterQuarkBarOfWMinus () const
	return anti-quark daughter of anti-W boson More...

const reco::GenParticle *	daughterQuarkBarOfWPlus () const
	return anti-quark daughter of W boson More...

const reco::GenParticle *	daughterQuarkOfTop (bool invertCharge=false) const
	return daughter quark of top quark (which can have flavor b, s or d) More...

const reco::GenParticle *	daughterQuarkOfTopBar () const
	return daughter quark of anti-top quark (which can have flavor b, s or d) More...

const reco::GenParticle *	daughterQuarkOfWMinus () const
	return quark daughter of anti-W boson More...

const reco::GenParticle *	daughterQuarkOfWPlus (bool invertQuarkCharge=false, bool invertBosonCharge=false) const
	return quark daughter quark of W boson More...

const reco::GenParticle *	eMinus () const
	return electron if available; 0 else More...

const reco::GenParticle *	ePlus () const
	return positron if available; 0 else More...

const reco::GenParticleCollection &	initialPartons () const
	return particles of initial partons More...

std::vector< const reco::GenParticle * >	lightQuarks (bool includingBQuarks=false) const
	return all light quarks or all quarks including b's More...

const reco::GenParticle *	muMinus () const
	return muon if available; 0 else More...

const reco::GenParticle *	muPlus () const
	return anti-muon if available; 0 else More...

int	numberOfBQuarks (bool fromTopQuark=true) const
	return number of b quarks in the decay chain More...

int	numberOfLeptons (bool fromWBoson=true) const
	return number of leptons in the decay chain More...

int	numberOfLeptons (WDecay::LeptonType type, bool fromWBoson=true) const
	return number of leptons in the decay chain More...

const reco::GenParticleCollection &	particles () const
	return particles of decay chain More...

void	print () const

std::vector< const reco::GenParticle * >	radiatedGluons (int pdgId) const
	return radiated gluons from particle with pdgId More...

const reco::GenParticle *	tauMinus () const
	return tau if available; 0 else More...

const reco::GenParticle *	tauPlus () const
	return anti-tau if available; 0 else More...

const reco::GenParticle *	top () const
	return top if available; 0 else More...

const reco::GenParticle *	topBar () const
	return anti-top if available; 0 else More...

	TopGenEvent ()
	empty constructor More...

	TopGenEvent (reco::GenParticleRefProd &decaySubset, reco::GenParticleRefProd &iniSubset)
	default constructor More...

std::vector< const reco::GenParticle * >	topSisters () const
	return number of top anti-top sisters More...

const reco::GenParticle *	wMinus () const
	return W minus if available; 0 else More...

const reco::GenParticle *	wPlus () const
	return W plus if available; 0 else More...

virtual	~TopGenEvent ()
	default destructor More...

Protected Attributes
math::XYZTLorentzVector	topPair_
	combined 4-vector of top and topBar More...

Protected Attributes inherited from TopGenEvent
reco::GenParticleRefProd	initPartons_
	reference to the list of initial partons (has to be kept in the event!) More...

reco::GenParticleRefProd	parts_
	reference to the top decay chain (has to be kept in the event!) More...

Private Member Functions
bool	isNumberOfLeptons (bool excludeTauLeptons, int nlep) const

Detailed Description

Class derived from the TopGenEvent for ttbar events.

The structure holds reference information to the generator particles of the decay chains for each top quark and of the initial partons and provides access and administration. The derived class contains a few additional getters with respect to its base class.

Definition at line 18 of file TtGenEvent.h.

Constructor & Destructor Documentation

◆ TtGenEvent() [1/2]

TtGenEvent::TtGenEvent ( )

inline

empty constructor

Definition at line 21 of file TtGenEvent.h.

21 {};

◆ TtGenEvent() [2/2]

TtGenEvent::TtGenEvent	(	reco::GenParticleRefProd &	decaySubset,
		reco::GenParticleRefProd &	initSubset
	)

default constructor from decaySubset and initSubset

Definition at line 10 of file TtGenEvent.cc.

                                                                                               {
   parts_ = decaySubset;
   initPartons_ = initSubset;
   if (top() && topBar())
     topPair_ = math::XYZTLorentzVector(top()->p4() + topBar()->p4());
 }

References TopDecaySubset_cfi::decaySubset, TopGenEvent::initPartons_, TopInitSubset_cfi::initSubset, p4, TopGenEvent::parts_, TopGenEvent::top(), TopGenEvent::topBar(), and topPair_.

◆ ~TtGenEvent()

TtGenEvent::~TtGenEvent ( )

inlineoverride

default destructor

Definition at line 25 of file TtGenEvent.h.

25 {};

Member Function Documentation

◆ fromGluonFusion()

bool TtGenEvent::fromGluonFusion ( ) const

check if the tops were produced from a pair of gluons

Definition at line 17 of file TtGenEvent.cc.

                                        {
   const reco::GenParticleCollection& initPartsColl = *initPartons_;
   if (initPartsColl.size() == 2)
     if (initPartsColl[0].pdgId() == 21 && initPartsColl[1].pdgId() == 21)
       return true;
   return false;
 }

References TopGenEvent::initPartons_, and EgammaValidation_cff::pdgId.

◆ fromQuarkAnnihilation()

bool TtGenEvent::fromQuarkAnnihilation ( ) const

check if the tops were produced from qqbar

Definition at line 25 of file TtGenEvent.cc.

                                              {
   const reco::GenParticleCollection& initPartsColl = *initPartons_;
   if (initPartsColl.size() == 2)
     if (std::abs(initPartsColl[0].pdgId()) < TopDecayID::tID && initPartsColl[0].pdgId() == -initPartsColl[1].pdgId())
       return true;
   return false;
 }

References funct::abs(), TopGenEvent::initPartons_, EgammaValidation_cff::pdgId, and TopDecayID::tID.

◆ fullLeptonicChannel()

std::pair< WDecay::LeptonType, WDecay::LeptonType > TtGenEvent::fullLeptonicChannel ( ) const

Definition at line 46 of file TtGenEvent.cc.

                                                                                   {
   WDecay::LeptonType typeA = WDecay::kNone, typeB = WDecay::kNone;
   if (isFullLeptonic()) {
     if (lepton()) {
       if (std::abs(lepton()->pdgId()) == TopDecayID::elecID)
         typeA = WDecay::kElec;
       if (std::abs(lepton()->pdgId()) == TopDecayID::muonID)
         typeA = WDecay::kMuon;
       if (std::abs(lepton()->pdgId()) == TopDecayID::tauID)
         typeA = WDecay::kTau;
     }
     if (leptonBar()) {
       if (std::abs(leptonBar()->pdgId()) == TopDecayID::elecID)
         typeB = WDecay::kElec;
       if (std::abs(leptonBar()->pdgId()) == TopDecayID::muonID)
         typeB = WDecay::kMuon;
       if (std::abs(leptonBar()->pdgId()) == TopDecayID::tauID)
         typeB = WDecay::kTau;
     }
   }
   return (std::pair<WDecay::LeptonType, WDecay::LeptonType>(typeA, typeB));
 }

References funct::abs(), TopDecayID::elecID, isFullLeptonic(), WDecay::kElec, WDecay::kMuon, WDecay::kNone, WDecay::kTau, lepton(), leptonBar(), TopDecayID::muonID, EgammaValidation_cff::pdgId, and TopDecayID::tauID.

Referenced by isFullLeptonic().

◆ hadronicDecayB()

const reco::GenParticle * TtGenEvent::hadronicDecayB ( bool excludeTauLeptons = false ) const

get b of hadronic decay branch

Definition at line 181 of file TtGenEvent.cc.

                                                                               {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::bID &&
           reco::flavour(singleLep) == reco::flavour(partsColl[i])) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), TopDecayID::bID, reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, and singleLepton().

Referenced by TtSemiEvtSolution::getGenHadb().

◆ hadronicDecayQuark()

const reco::GenParticle * TtGenEvent::hadronicDecayQuark ( bool invertFlavor = false ) const

get light quark of hadronic decay branch

Definition at line 153 of file TtGenEvent.cc.

                                                                              {
   const reco::GenParticle* cand = nullptr;
   // catch W boson and check its daughters for a quark;
   // make sure the decay is semi-leptonic first; this
   // only makes sense if taus are not excluded from the
   // decision
   if (singleLepton(false)) {
     for (reco::GenParticleCollection::const_iterator w = parts_->begin(); w != parts_->end(); ++w) {
       if (std::abs(w->pdgId()) == TopDecayID::WID) {
         // make sure that the particle is a W daughter
         for (reco::GenParticle::const_iterator wd = w->begin(); wd != w->end(); ++wd) {
           // make sure that the parton is a quark
           if (std::abs(wd->pdgId()) < TopDecayID::tID) {
             if (invertFlavor ? reco::flavour(*wd) < 0 : reco::flavour(*wd) > 0) {
               cand = dynamic_cast<const reco::GenParticle*>(&(*wd));
               if (cand == nullptr) {
                 throw edm::Exception(edm::errors::InvalidReference, "Not a GenParticle");
               }
               break;
             }
           }
         }
       }
     }
   }
   return cand;
 }

References funct::abs(), Exception, reco::flavour(), edm::errors::InvalidReference, TopGenEvent::parts_, singleLepton(), TopDecayID::tID, w, and TopDecayID::WID.

Referenced by TtSemiEvtSolution::getGenHadp(), and hadronicDecayQuarkBar().

◆ hadronicDecayQuarkBar()

const reco::GenParticle* TtGenEvent::hadronicDecayQuarkBar ( ) const

inline

get light anti-quark of hadronic decay branch

Definition at line 78 of file TtGenEvent.h.

78 { return hadronicDecayQuark(true); };

References hadronicDecayQuark().

Referenced by TtSemiEvtSolution::getGenHadq().

◆ hadronicDecayTop()

const reco::GenParticle * TtGenEvent::hadronicDecayTop ( bool excludeTauLeptons = false ) const

get top of hadronic decay branch

Definition at line 212 of file TtGenEvent.cc.

                                                                                 {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::tID &&
           reco::flavour(singleLep) == reco::flavour(partsColl[i])) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, singleLepton(), and TopDecayID::tID.

Referenced by TtSemiEvtSolution::getGenHadt(), and hadronicDecayTopRadiation().

◆ hadronicDecayTopRadiation()

std::vector< const reco::GenParticle * > TtGenEvent::hadronicDecayTopRadiation ( bool excludeTauLeptons = false ) const

gluons as radiated from the hadronicly decaying top quark

Definition at line 282 of file TtGenEvent.cc.

                                                                                                     {
   if (hadronicDecayTop(excludeTauLeptons)) {
     return (hadronicDecayTop(excludeTauLeptons)->pdgId() > 0 ? radiatedGluons(TopDecayID::tID)
                                                              : radiatedGluons(-TopDecayID::tID));
   }
   std::vector<const reco::GenParticle*> rad;
   return (rad);
 }

References hadronicDecayTop(), EgammaValidation_cff::pdgId, TopGenEvent::radiatedGluons(), and TopDecayID::tID.

◆ hadronicDecayW()

const reco::GenParticle * TtGenEvent::hadronicDecayW ( bool excludeTauLeptons = false ) const

get W of hadronic decay branch

Definition at line 196 of file TtGenEvent.cc.

                                                                               {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::WID &&
           reco::flavour(singleLep) != -reco::flavour(partsColl[i])) {
         // PDG Id:13=mu- 24=W+ (+24)->(-13) (-24)->(+13) opposite sign
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, singleLepton(), and TopDecayID::WID.

Referenced by TtSemiEvtSolution::getGenHadW().

◆ isFullHadronic()

bool TtGenEvent::isFullHadronic ( bool excludeTauLeptons = false ) const

inline

check if the event can be classified as full hadronic

Definition at line 34 of file TtGenEvent.h.

                                                             {
     return isTtBar() ? isNumberOfLeptons(excludeTauLeptons, 0) : false;
   }

References isNumberOfLeptons(), and isTtBar().

Referenced by TtHadEvtSolution::setGenEvt().

◆ isFullLeptonic() [1/2]

bool TtGenEvent::isFullLeptonic ( bool excludeTauLeptons = false ) const

inline

check if the event can be classified as full leptonic

Definition at line 42 of file TtGenEvent.h.

                                                             {
     return isTtBar() ? isNumberOfLeptons(excludeTauLeptons, 2) : false;
   }

References isNumberOfLeptons(), and isTtBar().

Referenced by fullLeptonicChannel(), and TtDilepEvtSolution::setGenEvt().

◆ isFullLeptonic() [2/2]

bool TtGenEvent::isFullLeptonic	(	WDecay::LeptonType	typeA,
		WDecay::LeptonType	typeB
	)		const

inline

check if the event is full leptonic with the lepton being of typeA or typeB irrelevant of order; all leptons including taus are allowed

Definition at line 107 of file TtGenEvent.h.

                                                                                            {
   return ((fullLeptonicChannel().first == typeA && fullLeptonicChannel().second == typeB) ||
           (fullLeptonicChannel().first == typeB && fullLeptonicChannel().second == typeA));
 }

References first, fullLeptonicChannel(), and edm::second().

◆ isNumberOfLeptons()

bool TtGenEvent::isNumberOfLeptons	(	bool	excludeTauLeptons,
		int	nlep
	)		const

inlineprivate

check whether the number of leptons among the daughters of the W boson is nlep or not; there is an option to exclude taus from the list of leptons to consider

Definition at line 102 of file TtGenEvent.h.

                                                                  {
     return excludeTauLeptons ? (numberOfLeptons() - numberOfLeptons(WDecay::kTau)) == nlep : numberOfLeptons() == nlep;
   }

References WDecay::kTau, and TopGenEvent::numberOfLeptons().

Referenced by isFullHadronic(), isFullLeptonic(), and isSemiLeptonic().

◆ isSemiLeptonic() [1/3]

bool TtGenEvent::isSemiLeptonic ( bool excludeTauLeptons = false ) const

inline

check if the event can be classified as semi-laptonic

Definition at line 38 of file TtGenEvent.h.

                                                             {
     return isTtBar() ? isNumberOfLeptons(excludeTauLeptons, 1) : false;
   }

References isNumberOfLeptons(), and isTtBar().

Referenced by semiLeptonicChannel(), TtSemiEvtSolution::setGenEvt(), singleLepton(), and singleNeutrino().

◆ isSemiLeptonic() [2/3]

bool TtGenEvent::isSemiLeptonic ( WDecay::LeptonType typeA ) const

inline

check if the event is semi-leptonic with the lepton being of typeA; all leptons including taus are allowed

Definition at line 49 of file TtGenEvent.h.

49 { return semiLeptonicChannel() == typeA ? true : false; };

References semiLeptonicChannel(), and funct::true.

◆ isSemiLeptonic() [3/3]

bool TtGenEvent::isSemiLeptonic	(	WDecay::LeptonType	typeA,
		WDecay::LeptonType	typeB
	)		const

inline

check if the event is semi-leptonic with the lepton being of typeA or typeB; all leptons including taus are allowed

Definition at line 51 of file TtGenEvent.h.

                                                                             {
     return (semiLeptonicChannel() == typeA || semiLeptonicChannel() == typeB) ? true : false;
   };

References semiLeptonicChannel(), and funct::true.

◆ isTtBar()

bool TtGenEvent::isTtBar ( ) const

inline

check if the event can be classified as ttbar

Definition at line 28 of file TtGenEvent.h.

28 { return (top() && topBar()); }

References TopGenEvent::top(), and TopGenEvent::topBar().

Referenced by isFullHadronic(), isFullLeptonic(), and isSemiLeptonic().

◆ lepton()

const reco::GenParticle * TtGenEvent::lepton ( bool excludeTauLeptons = false ) const

get lepton for semi-leptonic or full leptonic decays

Definition at line 69 of file TtGenEvent.cc.

                                                                       {
   const reco::GenParticle* cand = nullptr;
   const reco::GenParticleCollection& partsColl = *parts_;
   for (unsigned int i = 0; i < partsColl.size(); ++i) {
     if (reco::isLepton(partsColl[i]) && partsColl[i].mother() &&
         std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
       if (reco::flavour(partsColl[i]) > 0) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, reco::isLepton(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by fullLeptonicChannel(), and TtDilepEvtSolution::getGenLepm().

◆ leptonBar()

const reco::GenParticle * TtGenEvent::leptonBar ( bool excludeTauLeptons = false ) const

get anti-lepton for semi-leptonic or full leptonic decays

Definition at line 83 of file TtGenEvent.cc.

                                                                          {
   const reco::GenParticle* cand = nullptr;
   const reco::GenParticleCollection& partsColl = *parts_;
   for (unsigned int i = 0; i < partsColl.size(); ++i) {
     if (reco::isLepton(partsColl[i]) && partsColl[i].mother() &&
         std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
       if (reco::flavour(partsColl[i]) < 0) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, reco::isLepton(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by fullLeptonicChannel(), and TtDilepEvtSolution::getGenLepp().

◆ leptonicDecayB()

const reco::GenParticle * TtGenEvent::leptonicDecayB ( bool excludeTauLeptons = false ) const

get b of leptonic decay branch

Definition at line 227 of file TtGenEvent.cc.

                                                                               {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::bID &&
           reco::flavour(singleLep) != reco::flavour(partsColl[i])) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), TopDecayID::bID, reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, and singleLepton().

Referenced by TtSemiEvtSolution::getGenLepb().

◆ leptonicDecayTop()

const reco::GenParticle * TtGenEvent::leptonicDecayTop ( bool excludeTauLeptons = false ) const

get top of leptonic decay branch

Definition at line 258 of file TtGenEvent.cc.

                                                                                 {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::tID &&
           reco::flavour(singleLep) != reco::flavour(partsColl[i])) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, singleLepton(), and TopDecayID::tID.

Referenced by TtSemiEvtSolution::getGenLept(), and leptonicDecayTopRadiation().

◆ leptonicDecayTopRadiation()

std::vector< const reco::GenParticle * > TtGenEvent::leptonicDecayTopRadiation ( bool excludeTauLeptons = false ) const

gluons as radiated from the leptonicly decaying top quark

Definition at line 273 of file TtGenEvent.cc.

                                                                                                     {
   if (leptonicDecayTop(excludeTauLeptons)) {
     return (leptonicDecayTop(excludeTauLeptons)->pdgId() > 0 ? radiatedGluons(TopDecayID::tID)
                                                              : radiatedGluons(-TopDecayID::tID));
   }
   std::vector<const reco::GenParticle*> rad;
   return (rad);
 }

References leptonicDecayTop(), EgammaValidation_cff::pdgId, TopGenEvent::radiatedGluons(), and TopDecayID::tID.

◆ leptonicDecayW()

const reco::GenParticle * TtGenEvent::leptonicDecayW ( bool excludeTauLeptons = false ) const

get W of leptonic decay branch

Definition at line 242 of file TtGenEvent.cc.

                                                                               {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton(excludeTauLeptons);
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::WID &&
           reco::flavour(singleLep) == -reco::flavour(partsColl[i])) {
         // PDG Id:13=mu- 24=W+ (+24)->(-13) (-24)->(+13) opposite sign
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, TopGenEvent::parts_, EgammaValidation_cff::pdgId, singleLepton(), and TopDecayID::WID.

Referenced by TtSemiEvtSolution::getGenLepW().

◆ neutrino()

const reco::GenParticle * TtGenEvent::neutrino ( bool excludeTauLeptons = false ) const

get neutrino for semi-leptonic or full leptonic decays

Definition at line 111 of file TtGenEvent.cc.

                                                                         {
   const reco::GenParticle* cand = nullptr;
   const reco::GenParticleCollection& partsColl = *parts_;
   for (unsigned int i = 0; i < partsColl.size(); ++i) {
     if (reco::isNeutrino(partsColl[i]) && partsColl[i].mother() &&
         std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
       if (reco::flavour(partsColl[i]) > 0) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, reco::isNeutrino(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by TtDilepEvtSolution::getGenN().

◆ neutrinoBar()

const reco::GenParticle * TtGenEvent::neutrinoBar ( bool excludeTauLeptons = false ) const

get anti-neutrino for semi-leptonic or full leptonic decays

Definition at line 125 of file TtGenEvent.cc.

                                                                            {
   const reco::GenParticle* cand = nullptr;
   const reco::GenParticleCollection& partsColl = *parts_;
   for (unsigned int i = 0; i < partsColl.size(); ++i) {
     if (reco::isNeutrino(partsColl[i]) && partsColl[i].mother() &&
         std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
       if (reco::flavour(partsColl[i]) < 0) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), reco::flavour(), mps_fire::i, reco::isNeutrino(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by TtDilepEvtSolution::getGenNbar().

◆ semiLeptonicChannel()

WDecay::LeptonType TtGenEvent::semiLeptonicChannel ( ) const

return decay channel; all leptons including taus are allowed

Definition at line 33 of file TtGenEvent.cc.

                                                        {
   WDecay::LeptonType type = WDecay::kNone;
   if (isSemiLeptonic() && singleLepton()) {
     if (std::abs(singleLepton()->pdgId()) == TopDecayID::elecID)
       type = WDecay::kElec;
     if (std::abs(singleLepton()->pdgId()) == TopDecayID::muonID)
       type = WDecay::kMuon;
     if (std::abs(singleLepton()->pdgId()) == TopDecayID::tauID)
       type = WDecay::kTau;
   }
   return type;
 }

References funct::abs(), TopDecayID::elecID, isSemiLeptonic(), WDecay::kElec, WDecay::kMuon, WDecay::kNone, WDecay::kTau, TopDecayID::muonID, EgammaValidation_cff::pdgId, singleLepton(), and TopDecayID::tauID.

Referenced by isSemiLeptonic().

◆ singleLepton()

const reco::GenParticle * TtGenEvent::singleLepton ( bool excludeTauLeptons = false ) const

return single lepton if available; 0 else

Definition at line 97 of file TtGenEvent.cc.

                                                                             {
   const reco::GenParticle* cand = nullptr;
   if (isSemiLeptonic(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (reco::isLepton(partsColl[i]) && partsColl[i].mother() &&
           std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), mps_fire::i, reco::isLepton(), isSemiLeptonic(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by TtSemiEvtSolution::getGenLepl(), hadronicDecayB(), hadronicDecayQuark(), hadronicDecayTop(), hadronicDecayW(), leptonicDecayB(), leptonicDecayTop(), leptonicDecayW(), and semiLeptonicChannel().

◆ singleNeutrino()

const reco::GenParticle * TtGenEvent::singleNeutrino ( bool excludeTauLeptons = false ) const

return single neutrino if available; 0 else

Definition at line 139 of file TtGenEvent.cc.

                                                                               {
   const reco::GenParticle* cand = nullptr;
   if (isSemiLeptonic(excludeTauLeptons)) {
     const reco::GenParticleCollection& partsColl = *parts_;
     for (unsigned int i = 0; i < partsColl.size(); ++i) {
       if (reco::isNeutrino(partsColl[i]) && partsColl[i].mother() &&
           std::abs(partsColl[i].mother()->pdgId()) == TopDecayID::WID) {
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

References funct::abs(), mps_fire::i, reco::isNeutrino(), isSemiLeptonic(), TopGenEvent::parts_, EgammaValidation_cff::pdgId, and TopDecayID::WID.

Referenced by TtSemiEvtSolution::getGenLepn().

◆ topPair()

const math::XYZTLorentzVector* TtGenEvent::topPair ( ) const

inline

return combined 4-vector of top and topBar

Definition at line 93 of file TtGenEvent.h.

93 { return isTtBar() ? &topPair_ : nullptr; };

Member Data Documentation

◆ topPair_

math::XYZTLorentzVector TtGenEvent::topPair_

protected

combined 4-vector of top and topBar

Definition at line 93 of file TtGenEvent.h.

Referenced by TtGenEvent().

Public Member Functions

Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ TtGenEvent() [1/2]

◆ TtGenEvent() [2/2]

◆ ~TtGenEvent()

Member Function Documentation

◆ fromGluonFusion()

◆ fromQuarkAnnihilation()

◆ fullLeptonicChannel()

◆ hadronicDecayB()

◆ hadronicDecayQuark()

◆ hadronicDecayQuarkBar()

◆ hadronicDecayTop()

◆ hadronicDecayTopRadiation()

◆ hadronicDecayW()

◆ isFullHadronic()

◆ isFullLeptonic() [1/2]

◆ isFullLeptonic() [2/2]

◆ isNumberOfLeptons()

◆ isSemiLeptonic() [1/3]

◆ isSemiLeptonic() [2/3]

◆ isSemiLeptonic() [3/3]

◆ isTtBar()

◆ lepton()

◆ leptonBar()

◆ leptonicDecayB()

◆ leptonicDecayTop()

◆ leptonicDecayTopRadiation()

◆ leptonicDecayW()

◆ neutrino()

◆ neutrinoBar()

◆ semiLeptonicChannel()

◆ singleLepton()

◆ singleNeutrino()

◆ topPair()

Member Data Documentation

◆ topPair_