Class derived from the TtEvent for the semileptonic decay channel. More...

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

Inheritance diagram for TtSemiLeptonicEvent:

Public Member Functions
const reco::Candidate *	hadronicDecayB (const std::string &key, const unsigned &cmb=0) const
	get hadronic b of the given hypothesis More...

const reco::Candidate *	hadronicDecayB (const HypoClassKey &key, const unsigned &cmb=0) const
	get hadronic b of the given hypothesis More...

const reco::GenParticle *	hadronicDecayB () const
	get hadronic b of the TtGenEvent More...

const reco::Candidate *	hadronicDecayQuark (const std::string &key, const unsigned &cmb=0) const
	get hadronic light quark of the given hypothesis More...

const reco::Candidate *	hadronicDecayQuark (const HypoClassKey &key, const unsigned &cmb=0) const
	get hadronic light quark of the given hypothesis More...

const reco::GenParticle *	hadronicDecayQuark () const
	get hadronic light quark of the TtGenEvent More...

const reco::Candidate *	hadronicDecayQuarkBar (const std::string &key, const unsigned &cmb=0) const
	get hadronic light quark of the given hypothesis More...

const reco::Candidate *	hadronicDecayQuarkBar (const HypoClassKey &key, const unsigned &cmb=0) const
	get hadronic light quark of the given hypothesis More...

const reco::GenParticle *	hadronicDecayQuarkBar () const
	get hadronic light quark of the TtGenEvent More...

const reco::Candidate *	hadronicDecayTop (const std::string &key, const unsigned &cmb=0) const
	get hadronic top of the given hypothesis More...

const reco::Candidate *	hadronicDecayTop (const HypoClassKey &key, const unsigned &cmb=0) const
	get hadronic top of the given hypothesis More...

const reco::GenParticle *	hadronicDecayTop () const
	get hadronic top of the TtGenEvent More...

const reco::Candidate *	hadronicDecayW (const std::string &key, const unsigned &cmb=0) const
	get hadronic W of the given hypothesis More...

const reco::Candidate *	hadronicDecayW (const HypoClassKey &key, const unsigned &cmb=0) const
	get hadronic W of the given hypothesis More...

const reco::GenParticle *	hadronicDecayW () const
	get hadronic W of the TtGenEvent More...

const reco::Candidate *	leptonicDecayB (const std::string &key, const unsigned &cmb=0) const
	get leptonic b of the given hypothesis More...

const reco::Candidate *	leptonicDecayB (const HypoClassKey &key, const unsigned &cmb=0) const
	get leptonic b of the given hypothesis More...

const reco::GenParticle *	leptonicDecayB () const
	get leptonic b of the TtGenEvent More...

const reco::Candidate *	leptonicDecayTop (const std::string &key, const unsigned &cmb=0) const
	get leptonic top of the given hypothesis More...

const reco::Candidate *	leptonicDecayTop (const HypoClassKey &key, const unsigned &cmb=0) const
	get leptonic top of the given hypothesis More...

const reco::GenParticle *	leptonicDecayTop () const
	get leptonic top of the TtGenEvent More...

const reco::Candidate *	leptonicDecayW (const std::string &key, const unsigned &cmb=0) const
	get leptonic W of the given hypothesis More...

const reco::Candidate *	leptonicDecayW (const HypoClassKey &key, const unsigned &cmb=0) const
	get leptonic W of the given hypothesis More...

const reco::GenParticle *	leptonicDecayW () const
	get leptonic W of the TtGenEvent More...

const int	numberOfRealNeutrinoSolutions (const HypoClassKey &key) const
	get number of real neutrino solutions for a given hypo class More...

const int	numberOfRealNeutrinoSolutions (const std::string &key) const
	get number of real neutrino solutions for a given hypo class More...

void	print (const int verbosity=1) const

void	setNumberOfRealNeutrinoSolutions (const HypoClassKey &key, const int &nr)
	set number of real neutrino solutions for a given hypo class More...

const reco::Candidate *	singleLepton (const std::string &key, const unsigned &cmb=0) const
	get leptonic light quark of the given hypothesis More...

const reco::Candidate *	singleLepton (const HypoClassKey &key, const unsigned &cmb=0) const
	get leptonic light quark of the given hypothesis More...

const reco::GenParticle *	singleLepton () const
	get lepton top of the TtGenEvent More...

const reco::Candidate *	singleNeutrino (const std::string &key, const unsigned &cmb=0) const
	get leptonic light quark of the given hypothesis More...

const reco::Candidate *	singleNeutrino (const HypoClassKey &key, const unsigned &cmb=0) const
	get leptonic light quark of the given hypothesis More...

const reco::GenParticle *	singleNeutrino () const
	get neutrino of the TtGenEvent More...

	TtSemiLeptonicEvent ()
	empty constructor More...

	~TtSemiLeptonicEvent () override
	default destructor More...

Public Member Functions inherited from TtEvent
void	addEventHypo (const HypoClassKey &key, const HypoCombPair &hyp)
	add new hypotheses More...

int	correspondingHypo (const std::string &key1, const unsigned &hyp1, const std::string &key2) const
	return the hypothesis in hypothesis class 'key2', which corresponds to hypothesis 'hyp1' in hypothesis class 'key1' More...

int	correspondingHypo (const HypoClassKey &key1, const unsigned &hyp1, const HypoClassKey &key2) const
	return the hypothesis in hypothesis class 'key2', which corresponds to hypothesis 'hyp1' in hypothesis class 'key1' More...

const reco::CompositeCandidate &	eventHypo (const HypoClassKey &key, const unsigned &cmb=0) const

double	fitChi2 (const unsigned &cmb=0) const
	return the chi2 of the kinematic fit of hypothesis 'cmb' if available; -1 else More...

double	fitProb (const unsigned &cmb=0) const
	return the fit probability of hypothesis 'cmb' if available; -1 else More...

const edm::RefProd< TtGenEvent > &	genEvent () const
	get TtGenEvent More...

double	genMatchSumDR (const unsigned &cmb=0) const
	return the sum dr of the generator match if available; -1 else More...

double	genMatchSumPt (const unsigned &cmb=0) const
	return the sum pt of the generator match if available; -1 else More...

double	hitFitChi2 (const unsigned &cmb=0) const
	return the hitfit chi2 of hypothesis 'cmb' if available; -1 else More...

double	hitFitMT (const unsigned &cmb=0) const
	return the hitfit top mass of hypothesis 'cmb' if available; -1 else More...

double	hitFitProb (const unsigned &cmb=0) const
	return the hitfit probability of hypothesis 'cmb' if available; -1 else More...

double	hitFitSigMT (const unsigned &cmb=0) const
	return the hitfit top mass uncertainty of hypothesis 'cmb' if available; -1 else More...

bool	isHypoAvailable (const std::string &key, const unsigned &cmb=0) const

bool	isHypoAvailable (const HypoClassKey &key, const unsigned &cmb=0) const
	check if hypothesis 'cmb' is available within the hypothesis class More...

bool	isHypoClassAvailable (const std::string &key) const
	check if hypothesis class 'key' was added to the event structure More...

bool	isHypoClassAvailable (const HypoClassKey &key) const
	check if hypothesis class 'key' was added to the event structure More...

bool	isHypoValid (const std::string &key, const unsigned &cmb=0) const
	check if hypothesis 'cmb' within the hypothesis class was valid; if not it lead to an empty CompositeCandidate More...

bool	isHypoValid (const HypoClassKey &key, const unsigned &cmb=0) const
	check if hypothesis 'cmb' within the hypothesis class was valid; if not it lead to an empty CompositeCandidate More...

std::vector< int >	jetLeptonCombination (const std::string &key, const unsigned &cmb=0) const
	return the vector of jet lepton combinatorics for a given hypothesis and class More...

std::vector< int >	jetLeptonCombination (const HypoClassKey &key, const unsigned &cmb=0) const
	return the vector of jet lepton combinatorics for a given hypothesis and class More...

std::pair< WDecay::LeptonType, WDecay::LeptonType >	lepDecays () const
	get leptonic decay channels More...

double	mvaDisc (const unsigned &cmb=0) const
	return the mva discriminant value of hypothesis 'cmb' if available; -1 else More...

std::string	mvaMethod () const
	return the label of the mva method in use for the jet parton association (if kMVADisc is not available the string is empty) More...

unsigned int	numberOfAvailableHypoClasses () const
	return number of available hypothesis classes More...

unsigned int	numberOfAvailableHypos (const std::string &key) const
	return number of available hypotheses within a given hypothesis class More...

unsigned int	numberOfAvailableHypos (const HypoClassKey &key) const
	return number of available hypotheses within a given hypothesis class More...

int	numberOfConsideredJets (const std::string &key) const
	return number of jets that were considered when building a given hypothesis More...

int	numberOfConsideredJets (const HypoClassKey &key) const
	return number of jets that were considered when building a given hypothesis More...

void	setFitChi2 (const std::vector< double > &val)
	set chi2 of kKinFit hypothesis More...

void	setFitProb (const std::vector< double > &val)
	set fit probability of kKinFit hypothesis More...

void	setGenEvent (const edm::Handle< TtGenEvent > &evt)
	set TtGenEvent More...

void	setGenMatchSumDR (const std::vector< double > &val)
	set sum dr of kGenMatch hypothesis More...

void	setGenMatchSumPt (const std::vector< double > &val)
	set sum pt of kGenMatch hypothesis More...

void	setHitFitChi2 (const std::vector< double > &val)
	set chi2 of kHitFit hypothesis More...

void	setHitFitMT (const std::vector< double > &val)
	set fitted top mass of kHitFit hypothesis More...

void	setHitFitProb (const std::vector< double > &val)
	set fit probability of kHitFit hypothesis More...

void	setHitFitSigMT (const std::vector< double > &val)
	set fitted top mass uncertainty of kHitFit hypothesis More...

void	setLepDecays (const WDecay::LeptonType &lepDecTop1, const WDecay::LeptonType &lepDecTop2)
	set leptonic decay channels More...

void	setMvaDiscriminators (const std::vector< double > &val)
	set mva discriminant values of kMVADisc hypothesis More...

void	setMvaMethod (const std::string &name)
	set label of mva method for kMVADisc hypothesis More...

void	setNumberOfConsideredJets (const HypoClassKey &key, const unsigned int nJets)
	set number of jets considered when building a given hypothesis More...

const reco::Candidate *	topPair (const std::string &key, const unsigned &cmb=0) const
	get combined 4-vector of top and topBar of the given hypothesis More...

const reco::Candidate *	topPair (const HypoClassKey &key, const unsigned &cmb=0) const
	get combined 4-vector of top and topBar of the given hypothesis More...

const math::XYZTLorentzVector *	topPair () const
	get combined 4-vector of top and topBar from the TtGenEvent More...

	TtEvent ()
	empty constructor More...

virtual	~TtEvent ()
	default destructor More...

Protected Attributes
std::map< HypoClassKey, int >	numberOfRealNeutrinoSolutions_
	number of real neutrino solutions for all hypo classes More...

Protected Attributes inherited from TtEvent
std::map< HypoClassKey, std::vector< HypoCombPair > >	evtHyp_

std::vector< double >	fitChi2_
	result of kinematic fit More...

std::vector< double >	fitProb_
	result of kinematic fit More...

edm::RefProd< TtGenEvent >	genEvt_
	reference to TtGenEvent (has to be kept in the event!) More...

std::vector< double >	genMatchSumDR_
	result of gen match More...

std::vector< double >	genMatchSumPt_
	result of gen match More...

std::vector< double >	hitFitChi2_

std::vector< double >	hitFitMT_
	result of hitfit More...

std::vector< double >	hitFitProb_

std::vector< double >	hitFitSigMT_

std::pair< WDecay::LeptonType, WDecay::LeptonType >	lepDecays_
	leptonic decay channels More...

std::vector< double >	mvaDisc_
	MVA discriminants. More...

std::string	mvaMethod_
	label of the MVA method More...

std::map< HypoClassKey, int >	nJetsConsidered_
	number of jets considered when building the hypotheses More...

Additional Inherited Members
Public Types inherited from TtEvent
enum	HypoClassKey { kGeom, kWMassMaxSumPt, kMaxSumPtWMass, kGenMatch, kMVADisc, kKinFit, kKinSolution, kWMassDeltaTopMass, kHitFit }
	supported classes of event hypotheses More...

typedef std::pair < reco::CompositeCandidate, std::vector< int > >	HypoCombPair
	pair of hypothesis and lepton jet combinatorics for a given hypothesis More...

Protected Member Functions inherited from TtEvent
HypoClassKey	hypoClassKeyFromString (const std::string &label) const
	return the corresponding enum value from a string More...

Detailed Description

Class derived from the TtEvent for the semileptonic decay channel.

The structure holds information on the leptonic decay channels, all event hypotheses of different classes (user defined during production) and a reference to the TtGenEvent (if available). It provides access and administration.

Definition at line 24 of file TtSemiLeptonicEvent.h.

Constructor & Destructor Documentation

TtSemiLeptonicEvent::TtSemiLeptonicEvent ( )

inline

empty constructor

Definition at line 27 of file TtSemiLeptonicEvent.h.

27 {};

TtSemiLeptonicEvent::~TtSemiLeptonicEvent ( )

inlineoverride

default destructor

Definition at line 29 of file TtSemiLeptonicEvent.h.

29 {};

Member Function Documentation

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayB	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic b of the given hypothesis

Definition at line 40 of file TtSemiLeptonicEvent.h.

References hadronicDecayB(), and TtEvent::hypoClassKeyFromString().

                                                                                            {
     return hadronicDecayB(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayB	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic b of the given hypothesis

Definition at line 44 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtSemiLepDaughter::HadB, hadronicDecayTop(), and TtEvent::isHypoValid().

                                                                                               {
     return !isHypoValid(key, cmb) ? nullptr : hadronicDecayTop(key, cmb)->daughter(TtSemiLepDaughter::HadB);
   };

const reco::GenParticle* TtSemiLeptonicEvent::hadronicDecayB ( ) const

inline

get hadronic b of the TtGenEvent

Definition at line 117 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), TtEvent::genEvt_, and hadronicDecayB().

Referenced by hadronicDecayB(), and print().

117 { return (!genEvt_ ? nullptr : this->genEvent()->hadronicDecayB()); };

TtEvent::genEvent

const edm::RefProd< TtGenEvent > & genEvent() const

get TtGenEvent

Definition: TtEvent.h:63

TtSemiLeptonicEvent::hadronicDecayB

const reco::GenParticle * hadronicDecayB() const

get hadronic b of the TtGenEvent

Definition: TtSemiLeptonicEvent.h:117

TtEvent::genEvt_

edm::RefProd< TtGenEvent > genEvt_

reference to TtGenEvent (has to be kept in the event!)

Definition: TtEvent.h:188

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayQuark	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic light quark of the given hypothesis

Definition at line 56 of file TtSemiLeptonicEvent.h.

References hadronicDecayQuark(), and TtEvent::hypoClassKeyFromString().

                                                                                                {
     return hadronicDecayQuark(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayQuark	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic light quark of the given hypothesis

Definition at line 60 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtSemiLepDaughter::HadP, hadronicDecayW(), and TtEvent::isHypoValid().

                                                                                                   {
     return !isHypoValid(key, cmb) ? nullptr : hadronicDecayW(key, cmb)->daughter(TtSemiLepDaughter::HadP);
   };

const reco::GenParticle* TtSemiLeptonicEvent::hadronicDecayQuark ( ) const

inline

get hadronic light quark of the TtGenEvent

Definition at line 121 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), and TtEvent::genEvt_.

Referenced by hadronicDecayQuark(), and print().

                                                     {
     return (!genEvt_ ? nullptr : this->genEvent()->hadronicDecayQuark());
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayQuarkBar	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic light quark of the given hypothesis

Definition at line 64 of file TtSemiLeptonicEvent.h.

References hadronicDecayQuarkBar(), and TtEvent::hypoClassKeyFromString().

                                                                                                   {
     return hadronicDecayQuarkBar(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayQuarkBar	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic light quark of the given hypothesis

Definition at line 68 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtSemiLepDaughter::HadQ, hadronicDecayW(), and TtEvent::isHypoValid().

                                                                                                      {
     return !isHypoValid(key, cmb) ? nullptr : hadronicDecayW(key, cmb)->daughter(TtSemiLepDaughter::HadQ);
   };

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

inline

get hadronic light quark of the TtGenEvent

Definition at line 125 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), and TtEvent::genEvt_.

Referenced by hadronicDecayQuarkBar(), and print().

                                                        {
     return (!genEvt_ ? nullptr : this->genEvent()->hadronicDecayQuarkBar());
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayTop	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic top of the given hypothesis

Definition at line 32 of file TtSemiLeptonicEvent.h.

References hadronicDecayTop(), and TtEvent::hypoClassKeyFromString().

                                                                                              {
     return hadronicDecayTop(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayTop	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic top of the given hypothesis

Definition at line 36 of file TtSemiLeptonicEvent.h.

References reco::CompositeCandidate::daughter(), TtEvent::eventHypo(), TtSemiLepDaughter::HadTop, and TtEvent::isHypoValid().

                                                                                                 {
     return !isHypoValid(key, cmb) ? nullptr : eventHypo(key, cmb).daughter(TtSemiLepDaughter::HadTop);
   };

const reco::GenParticle* TtSemiLeptonicEvent::hadronicDecayTop ( ) const

inline

get hadronic top of the TtGenEvent

Definition at line 113 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), and TtEvent::genEvt_.

Referenced by hadronicDecayB(), hadronicDecayTop(), hadronicDecayW(), and print().

                                                   {
     return (!genEvt_ ? nullptr : this->genEvent()->hadronicDecayTop());
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayW	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic W of the given hypothesis

Definition at line 48 of file TtSemiLeptonicEvent.h.

References hadronicDecayW(), and TtEvent::hypoClassKeyFromString().

                                                                                            {
     return hadronicDecayW(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::hadronicDecayW	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get hadronic W of the given hypothesis

Definition at line 52 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), hadronicDecayTop(), TtSemiLepDaughter::HadW, and TtEvent::isHypoValid().

                                                                                               {
     return !isHypoValid(key, cmb) ? nullptr : hadronicDecayTop(key, cmb)->daughter(TtSemiLepDaughter::HadW);
   };

const reco::GenParticle* TtSemiLeptonicEvent::hadronicDecayW ( ) const

inline

get hadronic W of the TtGenEvent

Definition at line 119 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), TtEvent::genEvt_, and hadronicDecayW().

Referenced by hadronicDecayQuark(), hadronicDecayQuarkBar(), hadronicDecayW(), and print().

119 { return (!genEvt_ ? nullptr : this->genEvent()->hadronicDecayW()); };

TtEvent::genEvent

const edm::RefProd< TtGenEvent > & genEvent() const

get TtGenEvent

Definition: TtEvent.h:63

TtSemiLeptonicEvent::hadronicDecayW

const reco::GenParticle * hadronicDecayW() const

get hadronic W of the TtGenEvent

Definition: TtSemiLeptonicEvent.h:119

TtEvent::genEvt_

edm::RefProd< TtGenEvent > genEvt_

reference to TtGenEvent (has to be kept in the event!)

Definition: TtEvent.h:188

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayB	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic b of the given hypothesis

Definition at line 80 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and leptonicDecayB().

                                                                                            {
     return leptonicDecayB(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayB	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic b of the given hypothesis

Definition at line 84 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtEvent::isHypoValid(), TtSemiLepDaughter::LepB, and leptonicDecayTop().

                                                                                               {
     return !isHypoValid(key, cmb) ? nullptr : leptonicDecayTop(key, cmb)->daughter(TtSemiLepDaughter::LepB);
   };

const reco::GenParticle* TtSemiLeptonicEvent::leptonicDecayB ( ) const

inline

get leptonic b of the TtGenEvent

Definition at line 133 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), TtEvent::genEvt_, and leptonicDecayB().

Referenced by leptonicDecayB(), and print().

133 { return (!genEvt_ ? nullptr : this->genEvent()->leptonicDecayB()); };

TtEvent::genEvent

const edm::RefProd< TtGenEvent > & genEvent() const

get TtGenEvent

Definition: TtEvent.h:63

TtSemiLeptonicEvent::leptonicDecayB

const reco::GenParticle * leptonicDecayB() const

get leptonic b of the TtGenEvent

Definition: TtSemiLeptonicEvent.h:133

TtEvent::genEvt_

edm::RefProd< TtGenEvent > genEvt_

reference to TtGenEvent (has to be kept in the event!)

Definition: TtEvent.h:188

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayTop	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic top of the given hypothesis

Definition at line 72 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and leptonicDecayTop().

                                                                                              {
     return leptonicDecayTop(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayTop	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic top of the given hypothesis

Definition at line 76 of file TtSemiLeptonicEvent.h.

References reco::CompositeCandidate::daughter(), TtEvent::eventHypo(), TtEvent::isHypoValid(), and TtSemiLepDaughter::LepTop.

                                                                                                 {
     return !isHypoValid(key, cmb) ? nullptr : eventHypo(key, cmb).daughter(TtSemiLepDaughter::LepTop);
   };

const reco::GenParticle* TtSemiLeptonicEvent::leptonicDecayTop ( ) const

inline

get leptonic top of the TtGenEvent

Definition at line 129 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), and TtEvent::genEvt_.

Referenced by leptonicDecayB(), leptonicDecayTop(), leptonicDecayW(), and print().

                                                   {
     return (!genEvt_ ? nullptr : this->genEvent()->leptonicDecayTop());
   };

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayW	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic W of the given hypothesis

Definition at line 88 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and leptonicDecayW().

                                                                                            {
     return leptonicDecayW(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::leptonicDecayW	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic W of the given hypothesis

Definition at line 92 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtEvent::isHypoValid(), leptonicDecayTop(), and TtSemiLepDaughter::LepW.

                                                                                               {
     return !isHypoValid(key, cmb) ? nullptr : leptonicDecayTop(key, cmb)->daughter(TtSemiLepDaughter::LepW);
   };

const reco::GenParticle* TtSemiLeptonicEvent::leptonicDecayW ( ) const

inline

get leptonic W of the TtGenEvent

Definition at line 135 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), TtEvent::genEvt_, and leptonicDecayW().

Referenced by leptonicDecayW(), print(), singleLepton(), and singleNeutrino().

135 { return (!genEvt_ ? nullptr : this->genEvent()->leptonicDecayW()); };

TtEvent::genEvent

const edm::RefProd< TtGenEvent > & genEvent() const

get TtGenEvent

Definition: TtEvent.h:63

TtSemiLeptonicEvent::leptonicDecayW

const reco::GenParticle * leptonicDecayW() const

get leptonic W of the TtGenEvent

Definition: TtSemiLeptonicEvent.h:135

TtEvent::genEvt_

edm::RefProd< TtGenEvent > genEvt_

reference to TtGenEvent (has to be kept in the event!)

Definition: TtEvent.h:188

const int TtSemiLeptonicEvent::numberOfRealNeutrinoSolutions ( const HypoClassKey & key ) const

inline

get number of real neutrino solutions for a given hypo class

Definition at line 146 of file TtSemiLeptonicEvent.h.

References numberOfRealNeutrinoSolutions_.

Referenced by numberOfRealNeutrinoSolutions(), and print().

                                                                          {
     return (numberOfRealNeutrinoSolutions_.find(key) == numberOfRealNeutrinoSolutions_.end()
                 ? -999
                 : numberOfRealNeutrinoSolutions_.find(key)->second);
   };

const int TtSemiLeptonicEvent::numberOfRealNeutrinoSolutions ( const std::string & key ) const

inline

get number of real neutrino solutions for a given hypo class

Definition at line 152 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and numberOfRealNeutrinoSolutions().

                                                                       {
     return numberOfRealNeutrinoSolutions(hypoClassKeyFromString(key));
   };

void TtSemiLeptonicEvent::print ( const int verbosity = 1 ) const

print full content of the structure as formated LogInfo to the MessageLogger output for debugging

Definition at line 9 of file TtSemiLeptonicEvent.cc.

References TtEvent::evtHyp_, TtEvent::fitChi2(), TtEvent::fitProb(), TtEvent::genEvent(), TtEvent::genEvt_, TtEvent::genMatchSumDR(), TtEvent::genMatchSumPt(), TtSemiLepEvtPartons::HadB, hadronicDecayB(), hadronicDecayQuark(), hadronicDecayQuarkBar(), hadronicDecayTop(), hadronicDecayW(), TtEvent::hitFitChi2(), TtEvent::hitFitMT(), TtEvent::hitFitProb(), TtEvent::hitFitSigMT(), TtEvent::isHypoValid(), TtEvent::jetLeptonCombination(), fwrapper::jets, WDecay::kElec, TtEvent::kGenMatch, TtEvent::kGeom, TtEvent::kHitFit, TtEvent::kKinFit, TtEvent::kKinSolution, TtEvent::kMaxSumPtWMass, WDecay::kMuon, TtEvent::kMVADisc, WDecay::kTau, TtEvent::kWMassDeltaTopMass, TtEvent::kWMassMaxSumPt, TtSemiLepEvtPartons::LepB, TtSemiLepEvtPartons::Lepton, leptonicDecayB(), leptonicDecayTop(), leptonicDecayW(), TtSemiLepEvtPartons::LightQ, TtSemiLepEvtPartons::LightQBar, log, TtEvent::mvaDisc(), TtEvent::mvaMethod(), TtEvent::numberOfAvailableHypoClasses(), TtEvent::numberOfAvailableHypos(), TtEvent::numberOfConsideredJets(), numberOfRealNeutrinoSolutions(), ttevent::printParticle(), singleLepton(), singleNeutrino(), and TtEvent::topPair().

                                                          {
   if (verbosity % 10 <= 0)
     return;
 
   edm::LogInfo log("TtSemiLeptonicEvent");
 
   log << "++++++++++++++++++++++++++++++++++++++++++++++++++ \n";
 
   // get some information from the genEvent (if available)
   if (!genEvt_)
     log << " TtGenEvent not available! \n";
   else {
     log << " TtGenEvent says: ";
     if (!this->genEvent()->isTtBar())
       log << "Not TtBar";
     else if (this->genEvent()->isFullHadronic())
       log << "Fully Hadronic TtBar";
     else if (this->genEvent()->isFullLeptonic())
       log << "Fully Leptonic TtBar";
     else if (this->genEvent()->isSemiLeptonic()) {
       log << "Semi-leptonic TtBar, ";
       switch (this->genEvent()->semiLeptonicChannel()) {
         case WDecay::kElec:
           log << "Electron";
           break;
         case WDecay::kMuon:
           log << "Muon";
           break;
         case WDecay::kTau:
           log << "Tau";
           break;
         default:
           log << "Unknown";
           break;
       }
       log << " Channel";
     }
     log << "\n";
   }
 
   // get number of available hypothesis classes
   log << " Number of available event hypothesis classes: " << this->numberOfAvailableHypoClasses() << " \n";
 
   // create a legend for the jetLepComb
   log << " - JetLepComb: ";
   for (unsigned idx = 0; idx < 5; idx++) {
     switch (idx) {
       case TtSemiLepEvtPartons::LightQ:
         log << "LightP ";
         break;
       case TtSemiLepEvtPartons::LightQBar:
         log << "LightQ ";
         break;
       case TtSemiLepEvtPartons::HadB:
         log << " HadB  ";
         break;
       case TtSemiLepEvtPartons::LepB:
         log << " LepB  ";
         break;
       case TtSemiLepEvtPartons::Lepton:
         log << "Lepton ";
         break;
     }
   }
   log << "\n";
 
   // get details from the hypotheses
   typedef std::map<HypoClassKey, std::vector<HypoCombPair> >::const_iterator EventHypo;
   for (EventHypo hyp = evtHyp_.begin(); hyp != evtHyp_.end(); ++hyp) {
     HypoClassKey hypKey = (*hyp).first;
     // header for each hypothesis
     log << "-------------------------------------------------- \n";
     switch (hypKey) {
       case kGeom:
         log << " Geom";
         break;
       case kWMassMaxSumPt:
         log << " WMassMaxSumPt";
         break;
       case kMaxSumPtWMass:
         log << " MaxSumPtWMass";
         break;
       case kGenMatch:
         log << " GenMatch";
         break;
       case kMVADisc:
         log << " MVADisc";
         break;
       case kKinFit:
         log << " KinFit";
         break;
       case kKinSolution:
         log << " KinSolution not (yet) applicable to TtSemiLeptonicEvent --> skipping";
         continue;
       case kWMassDeltaTopMass:
         log << " WMassDeltaTopMass";
         break;
       case kHitFit:
         log << " HitFit";
         break;
       default:
         log << " Unknown TtEvent::HypoClassKey provided --> skipping";
         continue;
     }
     log << "-Hypothesis: \n";
     log << " * Number of real neutrino solutions: " << this->numberOfRealNeutrinoSolutions(hypKey) << "\n";
     log << " * Number of considered jets        : " << this->numberOfConsideredJets(hypKey) << "\n";
     unsigned nOfHyp = this->numberOfAvailableHypos(hypKey);
     if (nOfHyp > 1) {
       log << " * Number of stored jet combinations: " << nOfHyp << "\n";
       if (verbosity < 10)
         log << " The following was found to be the best one:\n";
     }
     // if verbosity level is smaller than 10, never show more than the best jet combination
     if (verbosity < 10)
       nOfHyp = 1;
     for (unsigned cmb = 0; cmb < nOfHyp; cmb++) {
       // check if hypothesis is valid
       if (!this->isHypoValid(hypKey, cmb))
         log << " * Not valid! \n";
       // get meta information for valid hypothesis
       else {
         // jetLepComb
         log << " * JetLepComb:";
         std::vector<int> jets = this->jetLeptonCombination(hypKey, cmb);
         for (unsigned int iJet = 0; iJet < jets.size(); iJet++) {
           log << "   " << jets[iJet] << "   ";
         }
         log << "\n";
         // specialties for some hypotheses
         switch (hypKey) {
           case kGenMatch:
             log << " * Sum(DeltaR) : " << this->genMatchSumDR(cmb) << " \n"
                 << " * Sum(DeltaPt): " << this->genMatchSumPt(cmb) << " \n";
             break;
           case kMVADisc:
             log << " * Method  : " << this->mvaMethod() << " \n"
                 << " * Discrim.: " << this->mvaDisc(cmb) << " \n";
             break;
           case kKinFit:
             log << " * Chi^2      : " << this->fitChi2(cmb) << " \n"
                 << " * Prob(Chi^2): " << this->fitProb(cmb) << " \n";
             break;
           case kHitFit:
             log << " * Chi^2      : " << this->hitFitChi2(cmb) << " \n"
                 << " * Prob(Chi^2): " << this->hitFitProb(cmb) << " \n"
                 << " * Top mass   : " << this->hitFitMT(cmb) << " +/- " << this->hitFitSigMT(cmb) << " \n";
             break;
           default:
             break;
         }
         // kinematic quantities of particles (if last digit of verbosity level > 1)
         if (verbosity % 10 >= 2) {
           log << " * Candidates (pt; eta; phi; mass):\n";
           if (verbosity % 10 >= 3)
             printParticle(log, "top pair", this->topPair(hypKey, cmb));
           printParticle(log, "hadronic top", this->hadronicDecayTop(hypKey, cmb));
           printParticle(log, "hadronic W  ", this->hadronicDecayW(hypKey, cmb));
           if (verbosity % 10 >= 3) {
             printParticle(log, "hadronic b  ", this->hadronicDecayB(hypKey, cmb));
             printParticle(log, "hadronic p  ", this->hadronicDecayQuark(hypKey, cmb));
             printParticle(log, "hadronic q  ", this->hadronicDecayQuarkBar(hypKey, cmb));
           }
           printParticle(log, "leptonic top", this->leptonicDecayTop(hypKey, cmb));
           printParticle(log, "leptonic W  ", this->leptonicDecayW(hypKey, cmb));
           if (verbosity % 10 >= 3) {
             printParticle(log, "leptonic b  ", this->leptonicDecayB(hypKey, cmb));
             printParticle(log, "lepton      ", this->singleLepton(hypKey, cmb));
             printParticle(log, "neutrino    ", this->singleNeutrino(hypKey, cmb));
           }
         }
       }
     }
   }
 
   log << "++++++++++++++++++++++++++++++++++++++++++++++++++";
 }

void TtSemiLeptonicEvent::setNumberOfRealNeutrinoSolutions	(	const HypoClassKey &	key,
		const int &	nr
	)

inline

set number of real neutrino solutions for a given hypo class

Definition at line 157 of file TtSemiLeptonicEvent.h.

References submitPVResolutionJobs::key, and numberOfRealNeutrinoSolutions_.

Referenced by TtEvtBuilder< C >::fillSpecific().

                                                                                 {
     numberOfRealNeutrinoSolutions_[key] = nr;
   };

const reco::Candidate* TtSemiLeptonicEvent::singleLepton	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic light quark of the given hypothesis

Definition at line 104 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and singleLepton().

                                                                                          {
     return singleLepton(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::singleLepton	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic light quark of the given hypothesis

Definition at line 108 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtEvent::isHypoValid(), TtSemiLepDaughter::Lep, and leptonicDecayW().

                                                                                             {
     return !isHypoValid(key, cmb) ? nullptr : leptonicDecayW(key, cmb)->daughter(TtSemiLepDaughter::Lep);
   };

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

inline

get lepton top of the TtGenEvent

Definition at line 137 of file TtSemiLeptonicEvent.h.

References TtEvent::genEvent(), TtEvent::genEvt_, and singleLepton().

Referenced by print(), and singleLepton().

137 { return (!genEvt_ ? nullptr : this->genEvent()->singleLepton()); };

TtEvent::genEvent

const edm::RefProd< TtGenEvent > & genEvent() const

get TtGenEvent

Definition: TtEvent.h:63

TtSemiLeptonicEvent::singleLepton

const reco::GenParticle * singleLepton() const

get lepton top of the TtGenEvent

Definition: TtSemiLeptonicEvent.h:137

TtEvent::genEvt_

edm::RefProd< TtGenEvent > genEvt_

reference to TtGenEvent (has to be kept in the event!)

Definition: TtEvent.h:188

const reco::Candidate* TtSemiLeptonicEvent::singleNeutrino	(	const std::string &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic light quark of the given hypothesis

Definition at line 96 of file TtSemiLeptonicEvent.h.

References TtEvent::hypoClassKeyFromString(), and singleNeutrino().

                                                                                            {
     return singleNeutrino(hypoClassKeyFromString(key), cmb);
   };

const reco::Candidate* TtSemiLeptonicEvent::singleNeutrino	(	const HypoClassKey &	key,
		const unsigned &	cmb = `0`
	)		const

inline

get leptonic light quark of the given hypothesis

Definition at line 100 of file TtSemiLeptonicEvent.h.

References reco::CompositeRefCandidateT< D >::daughter(), TtEvent::isHypoValid(), leptonicDecayW(), and TtSemiLepDaughter::Nu.

                                                                                               {
     return !isHypoValid(key, cmb) ? nullptr : leptonicDecayW(key, cmb)->daughter(TtSemiLepDaughter::Nu);
   };