TtFullLeptonicEvent Class Reference

Class derived from the TtEvent for the full leptonic decay channel. More...

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

Inheritance diagram for TtFullLeptonicEvent:

Public Member Functions
const reco::Candidate *	b (const std::string &key, const unsigned &cmb=0) const
	get b of the given hypothesis More...
const reco::Candidate *	b (const HypoClassKey &key, const unsigned &cmb=0) const
	get b of the given hypothesis More...
const reco::Candidate *	bBar (const std::string &key, const unsigned &cmb=0) const
	get anti-b of the given hypothesis More...
const reco::Candidate *	bBar (const HypoClassKey &key, const unsigned &cmb=0) const
	get anti-b of the given hypothesis More...
const reco::GenParticle *	genB () const
	get b of the TtGenEvent More...
const reco::GenParticle *	genBBar () const
	get anti-b of the TtGenEvent More...
const reco::GenParticle *	genLepton () const
	get lepton of the TtGenEvent More...
const reco::GenParticle *	genLeptonBar () const
	get anti-lepton of the TtGenEvent More...
const reco::GenParticle *	genNeutrino () const
	get neutrino of the TtGenEvent More...
const reco::GenParticle *	genNeutrinoBar () const
	get anti-neutrino of the TtGenEvent More...
const reco::GenParticle *	genTop () const
	get top of the TtGenEvent More...
const reco::GenParticle *	genTopBar () const
	get anti-top of the TtGenEvent More...
const reco::GenParticle *	genWMinus () const
	get Wminus of the TtGenEvent More...
const reco::GenParticle *	genWPlus () const
	get Wplus of the TtGenEvent More...
bool	isWrongCharge () const
	return if the kinematic solution of hypothesis 'cmb' is right or wrong charge if available; -1 else More...
const reco::Candidate *	lepton (const std::string &key, const unsigned &cmb=0) const
	get lepton of the given hypothesis More...
const reco::Candidate *	lepton (const HypoClassKey &key, const unsigned &cmb=0) const
	get lepton of the given hypothesis More...
const reco::Candidate *	leptonBar (const std::string &key, const unsigned &cmb=0) const
	get anti-lepton of the given hypothesis More...
const reco::Candidate *	leptonBar (const HypoClassKey &key, const unsigned &cmb=0) const
	get anti-lepton of the given hypothesis More...
const reco::Candidate *	neutrino (const std::string &key, const unsigned &cmb=0) const
	get neutrino of the given hypothesis More...
const reco::Candidate *	neutrino (const HypoClassKey &key, const unsigned &cmb=0) const
	get neutrino of the given hypothesis More...
const reco::Candidate *	neutrinoBar (const std::string &key, const unsigned &cmb=0) const
	get anti-neutrino of the given hypothesis More...
const reco::Candidate *	neutrinoBar (const HypoClassKey &key, const unsigned &cmb=0) const
	get anti-neutrino of the given hypothesis More...
void	print (const int verbosity=1) const
void	setSolWeight (const std::vector< double > &val)
	set weight of kKinSolution hypothesis More...
void	setWrongCharge (const bool &val)
	set right or wrong charge combination of kKinSolution hypothesis More...
double	solWeight (const unsigned &cmb=0) const
	return the weight of the kinematic solution of hypothesis 'cmb' if available; -1 else More...
const reco::Candidate *	top (const std::string &key, const unsigned &cmb=0) const
	get top of the given hypothesis More...
const reco::Candidate *	top (const HypoClassKey &key, const unsigned &cmb=0) const
	get top of the given hypothesis More...
const reco::Candidate *	topBar (const std::string &key, const unsigned &cmb=0) const
	get anti-top of the given hypothesis More...
const reco::Candidate *	topBar (const HypoClassKey &key, const unsigned &cmb=0) const
	get anti-top of the given hypothesis More...
	TtFullLeptonicEvent ()
	empty constructor More...
const reco::Candidate *	wMinus (const std::string &key, const unsigned &cmb=0) const
	get Wminus of the given hypothesis More...
const reco::Candidate *	wMinus (const HypoClassKey &key, const unsigned &cmb=0) const
	get Wminus of the given hypothesis More...
const reco::Candidate *	wPlus (const std::string &key, const unsigned &cmb=0) const
	get Wplus of the given hypothesis More...
const reco::Candidate *	wPlus (const HypoClassKey &key, const unsigned &cmb=0) const
	get Wplus of the given hypothesis More...
	~TtFullLeptonicEvent () 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	printParticle (edm::LogInfo &log, const char *name, const reco::Candidate *cand) const
	print pt, eta, phi, mass of a given candidate into an existing LogInfo 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::vector< double >	solWeight_
	result of kinematic solution More...
bool	wrongCharge_
	right/wrong charge booleans 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 full leptonic 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 TtFullLeptonicEvent.h.

Constructor & Destructor Documentation

TtFullLeptonicEvent::TtFullLeptonicEvent ( )

inline

empty constructor

Definition at line 28 of file TtFullLeptonicEvent.h.

{};

TtFullLeptonicEvent::~TtFullLeptonicEvent ( )

inlineoverride

default destructor

Definition at line 30 of file TtFullLeptonicEvent.h.

{};

Member Function Documentation

const reco::Candidate* TtFullLeptonicEvent::b ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get b of the given hypothesis

Definition at line 37 of file TtFullLeptonicEvent.h.

References b().

Referenced by b(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::b ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get b of the given hypothesis

Definition at line 39 of file TtFullLeptonicEvent.h.

References TtFullLepDaughter::B, and reco::Candidate::daughter().

{ return !isHypoValid(key,cmb) ? nullptr : top(key,cmb)->daughter(TtFullLepDaughter::B); };

const reco::Candidate* TtFullLeptonicEvent::bBar ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get anti-b of the given hypothesis

Definition at line 57 of file TtFullLeptonicEvent.h.

References bBar().

Referenced by bBar(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::bBar ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get anti-b of the given hypothesis

Definition at line 59 of file TtFullLeptonicEvent.h.

References TtFullLepDaughter::BBar, and reco::Candidate::daughter().

{ return !isHypoValid(key,cmb) ? nullptr : topBar(key,cmb)->daughter(TtFullLepDaughter::BBar  ); };

const reco::GenParticle* TtFullLeptonicEvent::genB ( ) const

inline

get b of the TtGenEvent

Definition at line 76 of file TtFullLeptonicEvent.h.

References b, and nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genBBar ( ) const

inline

get anti-b of the TtGenEvent

Definition at line 86 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genLepton ( ) const

inline

get lepton of the TtGenEvent

Definition at line 90 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genLeptonBar ( ) const

inline

get anti-lepton of the TtGenEvent

Definition at line 80 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genNeutrino ( ) const

inline

get neutrino of the TtGenEvent

Definition at line 82 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genNeutrinoBar ( ) const

inline

get anti-neutrino of the TtGenEvent

Definition at line 92 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genTop ( ) const

inline

get top of the TtGenEvent

Definition at line 74 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genTopBar ( ) const

inline

get anti-top of the TtGenEvent

Definition at line 84 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genWMinus ( ) const

inline

get Wminus of the TtGenEvent

Definition at line 88 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

const reco::GenParticle* TtFullLeptonicEvent::genWPlus ( ) const

inline

get Wplus of the TtGenEvent

Definition at line 78 of file TtFullLeptonicEvent.h.

References nano_cff::genEvent.

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

bool TtFullLeptonicEvent::isWrongCharge ( ) const

inline

return if the kinematic solution of hypothesis 'cmb' is right or wrong charge if available; -1 else

Definition at line 97 of file TtFullLeptonicEvent.h.

Referenced by print().

{ return wrongCharge_; }

const reco::Candidate* TtFullLeptonicEvent::lepton ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get lepton of the given hypothesis

Definition at line 65 of file TtFullLeptonicEvent.h.

References lepton().

Referenced by lepton(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::lepton ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get lepton of the given hypothesis

Definition at line 67 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::Lep.

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

const reco::Candidate* TtFullLeptonicEvent::leptonBar ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get anti-lepton of the given hypothesis

Definition at line 45 of file TtFullLeptonicEvent.h.

References leptonBar().

Referenced by leptonBar(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::leptonBar ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get anti-lepton of the given hypothesis

Definition at line 47 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::LepBar.

{ return !isHypoValid(key,cmb) ? nullptr : wPlus(key,cmb)->daughter(TtFullLepDaughter::LepBar); };

const reco::Candidate* TtFullLeptonicEvent::neutrino ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get neutrino of the given hypothesis

Definition at line 49 of file TtFullLeptonicEvent.h.

References neutrino().

Referenced by neutrino(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::neutrino ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get neutrino of the given hypothesis

Definition at line 51 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::Nu.

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

const reco::Candidate* TtFullLeptonicEvent::neutrinoBar ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get anti-neutrino of the given hypothesis

Definition at line 69 of file TtFullLeptonicEvent.h.

References neutrinoBar().

Referenced by neutrinoBar(), and print().

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

const reco::Candidate* TtFullLeptonicEvent::neutrinoBar ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get anti-neutrino of the given hypothesis

Definition at line 71 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::NuBar.

{ return !isHypoValid(key,cmb) ? nullptr : wMinus   (key,cmb)->daughter(TtFullLepDaughter::NuBar ); };

void TtFullLeptonicEvent::print

(

const int

verbosity = 1

)

const

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

Definition at line 7 of file TtFullLeptonicEvent.cc.

References b(), bBar(), TtEvent::evtHyp_, plotBeamSpotDB::first, TtEvent::genEvent(), TtEvent::genMatchSumDR(), TtEvent::genMatchSumPt(), TtEvent::isHypoValid(), isWrongCharge(), 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, lepton(), leptonBar(), cmsBatch::log, neutrino(), neutrinoBar(), TtEvent::numberOfAvailableHypoClasses(), TtEvent::numberOfAvailableHypos(), TtEvent::printParticle(), edm::second(), solWeight(), top(), topBar(), TtEvent::topPair(), wMinus(), and wPlus().

{
  if(verbosity%10 <= 0)
    return;

  edm::LogInfo log("TtFullLeptonicEvent");

  log << "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ \n";

  // get some information from the genEvent
  log << " TtGenEvent says: ";
  if( !this->genEvent()->isTtBar() )            log << "Not TtBar";
  else if( this->genEvent()->isFullHadronic() ) log << "Fully Hadronic TtBar";
  else if( this->genEvent()->isSemiLeptonic() ) log << "Semi-leptonic TtBar";
  else if( this->genEvent()->isFullLeptonic() ) {
    log << "Fully Leptonic TtBar, ";
    switch( this->genEvent()->fullLeptonicChannel().first ) {
    case WDecay::kElec : log << "Electron-"; break;
    case WDecay::kMuon : log << "Muon-"    ; break;
    case WDecay::kTau  : log << "Tau-"     ; break;
    default            : log << "Unknown-" ; break;
    }
    switch( this->genEvent()->fullLeptonicChannel().second ) {
    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: ";
  log << "  b    ";
  log << " bbar  ";
  log << " e1(+) ";
  log << " e2(-) ";
  log << " mu1(+)";
  log << " mu2(-)";
  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 not (yet) applicable to TtFullLeptonicEvent --> skipping"             ; continue;
    case kWMassMaxSumPt     : log << " WMassMaxSumPt not (yet) applicable to TtFullLeptonicEvent --> skipping"    ; continue;
    case kMaxSumPtWMass     : log << " MaxSumPtWMass not (yet) applicable to TtFullLeptonicEvent --> skipping"    ; continue;
    case kGenMatch          : log << " GenMatch"                                                                  ; break;
    case kMVADisc           : log << " MVADisc not (yet) applicable to TtFullLeptonicEvent --> skipping"          ; continue;
    case kKinFit            : log << " KinFit not (yet) applicable to TtFullLeptonicEvent --> skipping"           ; continue;
    case kKinSolution       : log << " KinSolution"                                                               ; break;
    case kWMassDeltaTopMass : log << " WMassDeltaTopMass not (yet) applicable to TtFullLeptonicEvent --> skipping"; continue;
    case kHitFit            : log << " HitFit not (yet) applicable to TtFullLeptonicEvent --> skipping"           ; continue;
    default                 : log << " Unknown TtEvent::HypoClassKey provided --> skipping"                       ; continue;
    }
    log << "-Hypothesis: \n";
    unsigned nOfHyp = this->numberOfAvailableHypos(hypKey);
    if(nOfHyp > 1) {
      log << " * Number of available 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 kKinSolution : log << " * Weight      : "  << this->solWeight(cmb)     << " \n"
                << " * isWrongCharge: " << this->isWrongCharge()    << " \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, "top         ", this->top  (hypKey, cmb));
      printParticle(log, "W plus      ", this->wPlus(hypKey, cmb));
      if(verbosity%10 >= 3) {
        printParticle(log, "b           ", this->b        (hypKey, cmb));
        printParticle(log, "leptonBar   ", this->leptonBar(hypKey, cmb));
        printParticle(log, "neutrino    ", this->neutrino (hypKey, cmb));
      }
      printParticle(log, "topBar      ", this->topBar(hypKey, cmb));
      printParticle(log, "W minus     ", this->wMinus(hypKey, cmb));
      if(verbosity%10 >= 3) {
        printParticle(log, "bBar        ", this->bBar       (hypKey, cmb));
        printParticle(log, "lepton      ", this->lepton     (hypKey, cmb));
        printParticle(log, "neutrinoBar ", this->neutrinoBar(hypKey, cmb));
      }
    }
      }
    }
  }

  log << "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++";
}

void TtFullLeptonicEvent::setSolWeight ( const std::vector< double > &  val )

inline

set weight of kKinSolution hypothesis

Definition at line 100 of file TtFullLeptonicEvent.h.

References heppy_batch::val.

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

{ solWeight_=val; }; 

void TtFullLeptonicEvent::setWrongCharge ( const bool &  val )

inline

set right or wrong charge combination of kKinSolution hypothesis

Definition at line 102 of file TtFullLeptonicEvent.h.

References reco::print(), heppy_batch::val, and HIPAlignmentAlgorithm_cfi::verbosity.

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

{ wrongCharge_=val; }; 

double TtFullLeptonicEvent::solWeight ( const unsigned &  cmb = 0 ) const

inline

return the weight of the kinematic solution of hypothesis 'cmb' if available; -1 else

Definition at line 95 of file TtFullLeptonicEvent.h.

Referenced by print().

{ return (cmb<solWeight_.size() ? solWeight_[cmb] : -1.); }    

const reco::Candidate* TtFullLeptonicEvent::top ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get top of the given hypothesis

Definition at line 33 of file TtFullLeptonicEvent.h.

References top().

Referenced by print(), and top().

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

const reco::Candidate* TtFullLeptonicEvent::top ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get top of the given hypothesis

Definition at line 35 of file TtFullLeptonicEvent.h.

References TtFullLepDaughter::Top.

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

const reco::Candidate* TtFullLeptonicEvent::topBar ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get anti-top of the given hypothesis

Definition at line 53 of file TtFullLeptonicEvent.h.

References topBar().

Referenced by print(), and topBar().

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

const reco::Candidate* TtFullLeptonicEvent::topBar ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get anti-top of the given hypothesis

Definition at line 55 of file TtFullLeptonicEvent.h.

References TtFullLepDaughter::TopBar.

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

const reco::Candidate* TtFullLeptonicEvent::wMinus ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get Wminus of the given hypothesis

Definition at line 61 of file TtFullLeptonicEvent.h.

References wMinus().

Referenced by print(), and wMinus().

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

const reco::Candidate* TtFullLeptonicEvent::wMinus ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get Wminus of the given hypothesis

Definition at line 63 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::WMinus.

{ return !isHypoValid(key,cmb) ? nullptr : topBar(key,cmb)->daughter(TtFullLepDaughter::WMinus); };

const reco::Candidate* TtFullLeptonicEvent::wPlus ( const std::string &  key, const unsigned &  cmb = 0  ) const

inline

get Wplus of the given hypothesis

Definition at line 41 of file TtFullLeptonicEvent.h.

References wPlus().

Referenced by print(), and wPlus().

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

const reco::Candidate* TtFullLeptonicEvent::wPlus ( const HypoClassKey &  key, const unsigned &  cmb = 0  ) const

inline

get Wplus of the given hypothesis

Definition at line 43 of file TtFullLeptonicEvent.h.

References reco::Candidate::daughter(), and TtFullLepDaughter::WPlus.

{ return !isHypoValid(key,cmb) ? nullptr : top(key,cmb)->daughter(TtFullLepDaughter::WPlus); };

Member Data Documentation

std::vector<double> TtFullLeptonicEvent::solWeight_

protected

result of kinematic solution

Definition at line 111 of file TtFullLeptonicEvent.h.

bool TtFullLeptonicEvent::wrongCharge_

protected

right/wrong charge booleans

Definition at line 113 of file TtFullLeptonicEvent.h.