3 |
|
|
4 |
|
namespace HbbAnalysis { |
5 |
|
|
6 |
+ |
void HbbEvent::retrieveLeadJets(const std::vector<HbbAnalysis::Jet> & aVec, |
7 |
+ |
const std::vector<HbbAnalysis::Jet> & aNonbVec, |
8 |
+ |
HbbAnalysis::Jet & aLead, |
9 |
+ |
HbbAnalysis::Jet & aNtoLead) |
10 |
+ |
{ |
11 |
+ |
if (aVec.size() == 0) return; |
12 |
+ |
if (aNonbVec.size() == 0) { |
13 |
+ |
aLead = aVec.at(0); |
14 |
+ |
if (aVec.size() > 1) aNtoLead = aVec.at(1); |
15 |
+ |
} |
16 |
+ |
else { |
17 |
+ |
if (aVec.at(0).recoVars().pT>aNonbVec.at(0).recoVars().pT){ |
18 |
+ |
aLead = aVec.at(0); |
19 |
+ |
if (aVec.size() > 1){ |
20 |
+ |
if (aVec.at(1).recoVars().pT>aNonbVec.at(0).recoVars().pT) aNtoLead = aVec.at(1); |
21 |
+ |
else aNtoLead = aNonbVec.at(0); |
22 |
+ |
} |
23 |
+ |
} |
24 |
+ |
else { |
25 |
+ |
aLead = aNonbVec.at(0); |
26 |
+ |
if (aNonbVec.size() > 1){ |
27 |
+ |
if (aNonbVec.at(1).recoVars().pT > aVec.at(0).recoVars().pT) aNtoLead = aNonbVec.at(1); |
28 |
+ |
else aNtoLead = aVec.at(0); |
29 |
+ |
} |
30 |
+ |
} |
31 |
+ |
} |
32 |
+ |
|
33 |
+ |
} |
34 |
+ |
|
35 |
|
|
36 |
|
} |