|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
|
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "AnalysisDataFormats/TopObjects/interface/TtFullLeptonicEvent.h"
#include "AnalysisDataFormats/TopObjects/interface/TtFullLepEvtPartons.h"
#include "printParticle.h"
using ttevent::printParticle;
// print info via MessageLogger
void TtFullLeptonicEvent::print(const int verbosity) const {
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 << "+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++";
}
|