Class derived from the TopGenEvent for single-top events. More...

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

Inheritance diagram for StGenEvent:

Public Member Functions
const reco::GenParticle *	associatedB () const
	return associated b More...

const reco::GenParticle *	decayB () const
	return decay b More...

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

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

const reco::GenParticle *	singleTop () const
	return single Top More...

const reco::GenParticle *	singleW () const
	return single W More...

	StGenEvent ()
	empty constructor More...

	StGenEvent (reco::GenParticleRefProd &, reco::GenParticleRefProd &)
	default constructor More...

	~StGenEvent () 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...

Additional Inherited Members
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...

Detailed Description

Class derived from the TopGenEvent for single-top 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 17 of file StGenEvent.h.

Constructor & Destructor Documentation

◆ StGenEvent() [1/2]

StGenEvent::StGenEvent ( )

empty constructor

Definition at line 9 of file StGenEvent.cc.

9 {}

◆ StGenEvent() [2/2]

StGenEvent::StGenEvent	(	reco::GenParticleRefProd &	parts,
		reco::GenParticleRefProd &	inits
	)

default constructor

Definition at line 11 of file StGenEvent.cc.

References TopGenEvent::initPartons_, and TopGenEvent::parts_.

                                                                                    {
   parts_ = parts;
   initPartons_ = inits;
 }

◆ ~StGenEvent()

StGenEvent::~StGenEvent ( )

override

default destructor

Definition at line 16 of file StGenEvent.cc.

16 {}

Member Function Documentation

◆ associatedB()

const reco::GenParticle * StGenEvent::associatedB ( ) const

return associated b

Definition at line 34 of file StGenEvent.cc.

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

                                                      {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton()) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton();
     for (unsigned int i = 0; i < parts_->size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::bID &&
           reco::flavour(singleLep) == reco::flavour(partsColl[i])) {
         // ... but it should be the opposite!
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

◆ decayB()

const reco::GenParticle * StGenEvent::decayB ( ) const

return decay b

Definition at line 18 of file StGenEvent.cc.

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

Referenced by StEvtSolution::getGenBottom().

                                                 {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton()) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton();
     for (unsigned int i = 0; i < parts_->size(); ++i) {
       if (std::abs(partsColl[i].pdgId()) == TopDecayID::bID &&
           reco::flavour(singleLep) == -reco::flavour(partsColl[i])) {
         // ... but it should be the opposite!
         cand = &partsColl[i];
       }
     }
   }
   return cand;
 }

◆ singleLepton()

const reco::GenParticle * StGenEvent::singleLepton ( ) const

return single lepton if available; 0 else

Definition at line 50 of file StGenEvent.cc.

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

Referenced by associatedB(), decayB(), StEvtSolution::getGenLepton(), singleTop(), and singleW().

                                                       {
   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) {
       cand = &partsColl[i];
     }
   }
   return cand;
 }

◆ singleNeutrino()

const reco::GenParticle * StGenEvent::singleNeutrino ( ) const

return single neutrino if available; 0 else

Definition at line 62 of file StGenEvent.cc.

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

Referenced by StEvtSolution::getGenNeutrino().

                                                         {
   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) {
       cand = &partsColl[i];
     }
   }
   return cand;
 }

◆ singleTop()

const reco::GenParticle * StGenEvent::singleTop ( ) const

return single Top

Definition at line 90 of file StGenEvent.cc.

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

Referenced by StEvtSolution::getGenLept().

                                                    {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton()) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton();
     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;
 }

◆ singleW()

const reco::GenParticle * StGenEvent::singleW ( ) const

return single W

Definition at line 74 of file StGenEvent.cc.

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

Referenced by StEvtSolution::getGenLepW().

                                                  {
   const reco::GenParticle* cand = nullptr;
   if (singleLepton()) {
     const reco::GenParticleCollection& partsColl = *parts_;
     const reco::GenParticle& singleLep = *singleLepton();
     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;
 }

Public Member Functions

Additional Inherited Members