Finds the origin of each heavy flavour hadron and associated jets to it. More...

#include <GenHFHadronMatcher.cc>

Inheritance diagram for GenHFHadronMatcher:

Public Member Functions
	GenHFHadronMatcher (const edm::ParameterSet &)
	constructor initialising producer products and config parameters More...

	~GenHFHadronMatcher ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
	description of the run-time parameters More...

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
int	analyzeMothers (const reco::Candidate thisParticle, int &topDaughterQId, int &topBarDaughterQId, std::vector< const reco::Candidate > &hadMothers, std::vector< std::vector< int > > &hadMothersIndices, std::set< const reco::Candidate * > *analyzedParticles, const int prevPartIndex)
	do a recursive search for the mother particles until the b-quark is found or the absolute mother is found More...

virtual void	beginJob ()

virtual void	beginLuminosityBlock (edm::LuminosityBlock &, edm::EventSetup const &)

virtual void	beginRun (edm::Run &, edm::EventSetup const &)

bool	checkForLoop (std::vector< const reco::Candidate * > &particleChain, const reco::Candidate *particle)

virtual void	endJob ()

virtual void	endLuminosityBlock (edm::LuminosityBlock &, edm::EventSetup const &)

virtual void	endRun (edm::Run &, edm::EventSetup const &)

std::vector< int >	findHadronJets (const reco::GenParticleCollection genParticles, const reco::JetFlavourInfoMatchingCollection jetFlavourInfos, std::vector< int > &hadIndex, std::vector< reco::GenParticle > &hadMothersGenPart, std::vector< std::vector< int > > &hadMothersIndices, std::vector< int > &hadLeptonIndex, std::vector< int > &hadLeptonHadIndex, std::vector< int > &hadLeptonViaTau, std::vector< int > &hadFlavour, std::vector< int > &hadFromTopWeakDecay, std::vector< int > &hadBHadronId)
	identify the jets that contain b-hadrons More...

int	findInMothers (int idx, std::vector< int > &mothChains, const std::vector< std::vector< int > > &hadMothersIndices, const std::vector< reco::GenParticle > &hadMothers, int status, int pdgId, bool pdgAbs, int stopId, int firstLast, bool verbose)
	helper function to find indices of particles with particular pdgId and status from the list of mothers of a given particle More...

bool	fixExtraSameFlavours (const unsigned int hadId, const std::vector< int > &hadIndices, const std::vector< reco::GenParticle > &hadMothers, const std::vector< std::vector< int > > &hadMothersIndices, const std::vector< int > &isFromTopWeakDecay, const std::vector< std::vector< int > > &LastQuarkIds, const std::vector< std::vector< int > > &LastQuarkMotherIds, std::vector< int > &lastQuarkIndices, std::vector< int > &hadronFlavour, std::set< int > &checkedHadronIds, const int lastQuarkIndex)

int	flavourSign (const int pdgId)
	Sign of the flavour (matter/antimatter) More...

std::string	getParticleName (int id) const

bool	hasHadronDaughter (const int flavour, const reco::Candidate *thisParticle)
	Check if the particle has bHadron among daughters. More...

int	idInList (std::vector< const reco::Candidate * > particleList, const reco::Candidate *particle)
	Check if the cpecified particle is already in the list of particles. More...

int	idInList (std::vector< int > list, const int value)

bool	isBaryonPdgId (const int flavour, const int pdgId)
	Check the pdgId if it represents a baryon of particular flavour. More...

bool	isHadron (const int flavour, const reco::Candidate *thisParticle)
	Check the pdgId of a given particle if it is a hadron. More...

bool	isHadronPdgId (const int flavour, const int pdgId)
	Check the pdgId if it represents a hadron of particular flavour. More...

bool	isMesonPdgId (const int flavour, const int pdgId)
	Check the pdgId if it represents a meson of particular flavour. More...

bool	isNeutralPdg (int pdgId)
	Check whether a given pdgId represents neutral particle. More...

virtual void	produce (edm::Event &, const edm::EventSetup &)

bool	putMotherIndex (std::vector< std::vector< int > > &hadMothersIndices, int partIndex, int mothIndex)
	puts mother index to the list of mothers of particle, if it isn't there already More...

Private Attributes
int	flavour_

std::string	flavourStr_

edm::EDGetTokenT < reco::GenParticleCollection >	genParticlesToken_

edm::EDGetTokenT < reco::JetFlavourInfoMatchingCollection >	jetFlavourInfosToken_

bool	noBBbarResonances_

bool	onlyJetClusteredHadrons_

edm::ESHandle< ParticleDataTable >	pdt_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Finds the origin of each heavy flavour hadron and associated jets to it.

Starting from each consituent of each jet, tracks back in chain to find heavy flavour hadrons. From each hadron traces back until finds the b quark and its mother. For each hadron identifies the jet to which it was injected as a ghost hadron.

The description of the run-time parameters can be found at fillDescriptions()

The description of the products can be found at GenHFHadronMatcher()

Definition at line 58 of file GenHFHadronMatcher.cc.

Constructor & Destructor Documentation

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

explicit

constructor initialising producer products and config parameters

Warning: Definition of b-hadron and anti-b-hadron: The term b-hadron and anti-b-hadron is in reference to the quark content and not to distinguish particles from anti-particles. Here a b-hadron contains a b-quark and an anti-b-hadron contains an anti-b-quark. For mesons this means an inversion with respect to the PDG definition, as mesons actually contain anti-b-quarks and anti-mesons contain b-quarks.

Definition at line 139 of file GenHFHadronMatcher.cc.

GenHFHadronMatcher::~GenHFHadronMatcher ( )

Definition at line 169 of file GenHFHadronMatcher.cc.

Member Function Documentation

int GenHFHadronMatcher::analyzeMothers	(	const reco::Candidate *	thisParticle,
		int &	topDaughterQId,
		int &	topBarDaughterQId,
		std::vector< const reco::Candidate * > &	hadMothers,
		std::vector< std::vector< int > > &	hadMothersIndices,
		std::set< const reco::Candidate * > *	analyzedParticles,
		const int	prevPartIndex
	)

private

do a recursive search for the mother particles until the b-quark is found or the absolute mother is found

the treatment of b-bar resonances depends on the global parameter noBBbarResonances_

Parameters

[in]	thisParticle	current particle from which starts the search of the hadron and all its mothers up to proton
[out]	hadron	the last hadron in the decay chain [that decays weekly]
[out]	lepton	lepton found in the current decay chain
[out]	topDaughterQId	Id of the top quark daughter b(c) quark
[out]	topBarDaughterQId	Id of the antitop quark daughter b(c) quark
[out]	hadMothers	list of all processed particles ending with proton
[out]	hadMothersIndices	list of i-vectors containing j-indices representing particles that are mothers of each i-particle from hadMothers
[out]	analyzedParticles	list of particles analysed in the chain [used for infinite loop detection]
[out]	prevPartIndex	index of the previous particle in the current chain [used for infinite loop detection]

Returns: index of hadron in the hadMothers list [-1 if no hadron found]

Definition at line 644 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 237 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::beginLuminosityBlock	(	edm::LuminosityBlock &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 258 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::beginRun	(	edm::Run &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 247 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::checkForLoop	(	std::vector< const reco::Candidate * > &	particleChain,
		const reco::Candidate *	particle
	)

private

void GenHFHadronMatcher::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 242 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::endLuminosityBlock	(	edm::LuminosityBlock &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 263 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::endRun	(	edm::Run &	,
		edm::EventSetup const &
	)

privatevirtual

Definition at line 253 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

description of the run-time parameters

Definition at line 179 of file GenHFHadronMatcher.cc.

std::vector< int > GenHFHadronMatcher::findHadronJets	(	const reco::GenParticleCollection *	genParticles,
		const reco::JetFlavourInfoMatchingCollection *	jetFlavourInfos,
		std::vector< int > &	hadIndex,
		std::vector< reco::GenParticle > &	hadMothers,
		std::vector< std::vector< int > > &	hadMothersIndices,
		std::vector< int > &	hadLeptonIndex,
		std::vector< int > &	hadLeptonHadIndex,
		std::vector< int > &	hadLeptonViaTau,
		std::vector< int > &	hadFlavour,
		std::vector< int > &	hadFromTopWeakDecay,
		std::vector< int > &	hadBHadronId
	)

private

identify the jets that contain b-hadrons

For each hadron all mothers from all levels and chains are analysed to find the quark or gluon from which the hadron has originated. This is done by searching through the generator particle decay tree starting from the hadron, performing checks for flavour and kinematic consistency. The hadrons that are not last in the decay chain (don't decay weakly) are skipped. Additionally for each hadron it is checked whether it comes from the top weak decay branch and whether it comes from some other b-hadron decay

b-bbar (c-cbar) resonances can either be considered as hadrons or not, depending on the configuration.

Parameters

[out]	hadIndex	vector of indices of found hadrons in hadMothers
[out]	hadMothers	vector of all mothers at all levels of each found hadron
[out]	hadMothersIndices	connection between each particle from hadMothers and its mothers
[out]	hadLeptonIndex	index of lepton among the hadMothers
[out]	hadLeptonHadIndex	index of hadron associated to each lepton
[out]	hadLeptonViaTau	whether lepton is from direct hadron->lepton or hadron->tau->lepton
[out]	hadFlavour	flavour of each found hadron
[out]	hadFromTopWeakDecay	whether each hadron contains the top quark in its decay chain [works only for B-Hadrons]
[out]	hadBHadronId	for each hadron - index of the ancestor b-hadron [-1 if hadron doesn't come from another b-hdaron]

Returns: vector of jet indices that were matched to each hadron [by the jet clustering algorithm]

Definition at line 291 of file GenHFHadronMatcher.cc.

int GenHFHadronMatcher::findInMothers	(	int	idx,
		std::vector< int > &	mothChains,
		const std::vector< std::vector< int > > &	hadMothersIndices,
		const std::vector< reco::GenParticle > &	hadMothers,
		int	status,
		int	pdgId,
		bool	pdgAbs = `false`,
		int	stopId = `-1`,
		int	firstLast = `0`,
		bool	verbose = `false`
	)

private

helper function to find indices of particles with particular pdgId and status from the list of mothers of a given particle

Parameters

[in]	idx	index of particle, mothers of which should be searched
[in]	mothChains	vector of indices where the found mothers are stored
[in]	hadMothersIndices	list of indices pointing to mothers of each particle from list of mothers
[in]	hadMothers	vector of all hadron mother particles of all levels
[in]	status	status of mother that is being looked for
[in]	pdgId	pdgId of mother that is being looked for [flavour*111111 used to identify hadrons of particular flavour]
[in]	pdgAbs	whether the sign of pdgId should be taken into account
[in]	stopId	id of the particle in the hadMothers array after which the checking should stop
[in]	firstLast	should all(0), first(1) or last(2) occurances of the searched particle be stored
[in]	verbose	option to print every particle that is processed during the search

Returns: index of the found particle in the hadMothers array [-1 if the specified particle not found]

Definition at line 791 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::fixExtraSameFlavours	(	const unsigned int	hadId,
		const std::vector< int > &	hadIndices,
		const std::vector< reco::GenParticle > &	hadMothers,
		const std::vector< std::vector< int > > &	hadMothersIndices,
		const std::vector< int > &	isFromTopWeakDecay,
		const std::vector< std::vector< int > > &	LastQuarkIds,
		const std::vector< std::vector< int > > &	LastQuarkMotherIds,
		std::vector< int > &	lastQuarkIndices,
		std::vector< int > &	hadronFlavour,
		std::set< int > &	checkedHadronIds,
		const int	lastQuarkIndex
	)

private

Finds hadrons that have the same flavour and origin and resolve this ambiguity fixExtraSameFlavours

Parameters

hadId	Index of the hadron being checked
hadIndex	Vector of indices of each hadron pointing to the hadMothers
hadMothers	Vector of gen particles (contain hadrons and all its ancestors)
hadMothersIndices	Vector of indices for each particle from hadMothers
isFromTopWeakDecay	Vector of values showing whether particle comes from the top weak decay
LastQuarkIds	Vector of indices of last quarks for each hadron
LastQuarkMotherIds	Vector of indices of mothers for each last quark from LastQuarkIds
lastQuakIndices	Vector of indices pointing to particular index from the LastQuarkIds and LastQuarkMotherIds to be used for each hadron
lastQuarkIndex	Index from the LastQuarkIds and LastQuarkMotherIds for this particular hadron with index hadId

Returns: Whether other mother with unique association has been found for the hadrons

Definition at line 936 of file GenHFHadronMatcher.cc.

int GenHFHadronMatcher::flavourSign ( const int pdgId )

private

Sign of the flavour (matter/antimatter)

Parameters

[in] pdgId pdgId to be checked

Returns: +1/-1/0 matter/antimatter/undefined

Definition at line 591 of file GenHFHadronMatcher.cc.

std::string GenHFHadronMatcher::getParticleName ( int id ) const

private

bool GenHFHadronMatcher::hasHadronDaughter	(	const int	flavour,
		const reco::Candidate *	thisParticle
	)

private

Check if the particle has bHadron among daughters.

Parameters

[in]	flavour	flavour of a hadron that is being searched (5-B, 4-C)
[in]	thisParticle	a particle that is to be analysed

Returns: whether the particle has a hadron among its daughters

Definition at line 611 of file GenHFHadronMatcher.cc.

int GenHFHadronMatcher::idInList	(	std::vector< const reco::Candidate * >	particleList,
		const reco::Candidate *	particle
	)

private

Check if the cpecified particle is already in the list of particles.

Parameters

[in]	particleList	list of particles to be checked
[in]	particle	particle that should be checked

Returns: the index of the particle in the list [-1 if particle not found]

Definition at line 493 of file GenHFHadronMatcher.cc.

int GenHFHadronMatcher::idInList	(	std::vector< int >	list,
		const int	value
	)

private

Definition at line 501 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::isBaryonPdgId	(	const int	flavour,
		const int	pdgId
	)

private

Check the pdgId if it represents a baryon of particular flavour.

Parameters

[in]	flavour	flavour of a hadron that is being searched (5-B, 4-C)
[in]	pdgId	pdgId to be checked

Returns: true if the pdgId represents a baryon of specified flavour

Definition at line 572 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::isHadron	(	const int	flavour,
		const reco::Candidate *	thisParticle
	)

private

Check the pdgId of a given particle if it is a hadron.

Parameters

[in]	flavour	flavour of a hadron that is being searched (5-B, 4-C)
[in]	thisParticle	a particle that is to be analysed

Returns: whether the particle is a hadron of specified flavour

Definition at line 518 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::isHadronPdgId	(	const int	flavour,
		const int	pdgId
	)

private

Check the pdgId if it represents a hadron of particular flavour.

Parameters

[in]	flavour	flavour of a hadron that is being searched (5-B, 4-C)
[in]	pdgId	pdgId to be checked

Returns: true if the pdgId represents a hadron of specified flavour

Definition at line 532 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::isMesonPdgId	(	const int	flavour,
		const int	pdgId
	)

private

Check the pdgId if it represents a meson of particular flavour.

Parameters

[in]	flavour	flavour of a hadron that is being searched (5-B, 4-C)
[in]	pdgId	pdgId to be checked

Returns: true if the pdgId represents a meson of specified flavour

Definition at line 548 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::isNeutralPdg ( int pdgId )

private

Check whether a given pdgId represents neutral particle.

Parameters

[in]	pdgId
[in]	thisParticle	- a particle that is to be analysed

Returns: if the particle has a hadron among its daughters

Definition at line 912 of file GenHFHadronMatcher.cc.

void GenHFHadronMatcher::produce	(	edm::Event &	evt,
		const edm::EventSetup &	setup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 197 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::putMotherIndex	(	std::vector< std::vector< int > > &	hadMothersIndices,
		int	partIndex,
		int	mothIndex
	)

private

puts mother index to the list of mothers of particle, if it isn't there already

Parameters

[in]	hadMothersIndices	vector of indices of mothers for each particle
[in]	partIndex	index of the particle for which the mother index should be stored
[in]	mothIndex	index of mother that should be stored for current particle

Returns: whether the particle index was alreade in the list

Definition at line 738 of file GenHFHadronMatcher.cc.

Member Data Documentation

int GenHFHadronMatcher::flavour_

private

Definition at line 110 of file GenHFHadronMatcher.cc.

std::string GenHFHadronMatcher::flavourStr_

private

Definition at line 114 of file GenHFHadronMatcher.cc.

edm::EDGetTokenT<reco::GenParticleCollection> GenHFHadronMatcher::genParticlesToken_

private

Definition at line 108 of file GenHFHadronMatcher.cc.

edm::EDGetTokenT<reco::JetFlavourInfoMatchingCollection> GenHFHadronMatcher::jetFlavourInfosToken_

private

Definition at line 109 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::noBBbarResonances_

private

Definition at line 111 of file GenHFHadronMatcher.cc.

bool GenHFHadronMatcher::onlyJetClusteredHadrons_

private

Definition at line 112 of file GenHFHadronMatcher.cc.

edm::ESHandle<ParticleDataTable> GenHFHadronMatcher::pdt_

private

Definition at line 116 of file GenHFHadronMatcher.cc.

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation