7 if (verbosity % 10 <= 0)
12 log <<
"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ \n";
15 log <<
" TtGenEvent says: ";
18 else if (this->
genEvent()->isFullHadronic())
19 log <<
"Fully Hadronic TtBar";
20 else if (this->
genEvent()->isSemiLeptonic())
21 log <<
"Semi-leptonic TtBar";
22 else if (this->
genEvent()->isFullLeptonic()) {
23 log <<
"Fully Leptonic TtBar, ";
60 log <<
" - JetLepComb: ";
70 typedef std::map<HypoClassKey, std::vector<HypoCombPair> >::const_iterator EventHypo;
71 for (EventHypo hyp =
evtHyp_.begin(); hyp !=
evtHyp_.end(); ++hyp) {
74 log <<
"------------------------------------------------------------ \n";
77 log <<
" Geom not (yet) applicable to TtFullLeptonicEvent --> skipping";
80 log <<
" WMassMaxSumPt not (yet) applicable to TtFullLeptonicEvent --> skipping";
83 log <<
" MaxSumPtWMass not (yet) applicable to TtFullLeptonicEvent --> skipping";
89 log <<
" MVADisc not (yet) applicable to TtFullLeptonicEvent --> skipping";
92 log <<
" KinFit not (yet) applicable to TtFullLeptonicEvent --> skipping";
95 log <<
" KinSolution";
98 log <<
" WMassDeltaTopMass not (yet) applicable to TtFullLeptonicEvent --> skipping";
101 log <<
" HitFit not (yet) applicable to TtFullLeptonicEvent --> skipping";
104 log <<
" Unknown TtEvent::HypoClassKey provided --> skipping";
107 log <<
"-Hypothesis: \n";
110 log <<
" * Number of available jet combinations: " << nOfHyp <<
"\n";
112 log <<
" The following was found to be the best one:\n";
117 for (
unsigned cmb = 0; cmb < nOfHyp; cmb++) {
120 log <<
" * Not valid! \n";
124 log <<
" * JetLepComb:";
126 for (
unsigned int iJet = 0; iJet < jets.size(); iJet++) {
127 log <<
" " << jets[iJet] <<
" ";
133 log <<
" * Sum(DeltaR) : " << this->
genMatchSumDR(cmb) <<
" \n" 137 log <<
" * Weight : " << this->
solWeight(cmb) <<
" \n" 144 if (verbosity % 10 >= 2) {
145 log <<
" * Candidates (pt; eta; phi; mass):\n";
146 if (verbosity % 10 >= 3)
150 if (verbosity % 10 >= 3) {
157 if (verbosity % 10 >= 3) {
167 log <<
"+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++";
HypoClassKey
supported classes of event hypotheses
const reco::Candidate * bBar(const std::string &key, const unsigned &cmb=0) const
get anti-b of the given hypothesis
const edm::RefProd< TtGenEvent > & genEvent() const
get TtGenEvent
const math::XYZTLorentzVector * topPair() const
get combined 4-vector of top and topBar from the TtGenEvent
unsigned int numberOfAvailableHypos(const std::string &key) const
return number of available hypotheses within a given hypothesis class
const reco::Candidate * top(const std::string &key, const unsigned &cmb=0) const
get top of the given hypothesis
unsigned int numberOfAvailableHypoClasses() const
return number of available hypothesis classes
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
const reco::Candidate * lepton(const std::string &key, const unsigned &cmb=0) const
get lepton of the given hypothesis
U second(std::pair< T, U > const &p)
double genMatchSumDR(const unsigned &cmb=0) const
return the sum dr of the generator match if available; -1 else
void print(const int verbosity=1) const
const reco::Candidate * neutrino(const std::string &key, const unsigned &cmb=0) const
get neutrino of the given hypothesis
bool isWrongCharge() const
return if the kinematic solution of hypothesis 'cmb' is right or wrong charge if available; -1 else ...
const reco::Candidate * neutrinoBar(const std::string &key, const unsigned &cmb=0) const
get anti-neutrino of the given hypothesis
const reco::Candidate * b(const std::string &key, const unsigned &cmb=0) const
get b of the given hypothesis
const reco::Candidate * topBar(const std::string &key, const unsigned &cmb=0) const
get anti-top of the given hypothesis
double solWeight(const unsigned &cmb=0) const
return the weight of the kinematic solution of hypothesis 'cmb' if available; -1 else ...
const reco::Candidate * wPlus(const std::string &key, const unsigned &cmb=0) const
get Wplus of the given hypothesis
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 Composite...
double genMatchSumPt(const unsigned &cmb=0) const
return the sum pt of the generator match if available; -1 else
const reco::Candidate * wMinus(const std::string &key, const unsigned &cmb=0) const
get Wminus of the given hypothesis
std::map< HypoClassKey, std::vector< HypoCombPair > > evtHyp_
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
const reco::Candidate * leptonBar(const std::string &key, const unsigned &cmb=0) const
get anti-lepton of the given hypothesis