3 |
|
|
4 |
|
#include "VHbbAnalysis/HbbAnalyzer/interface/VHbbCandidate.h" |
5 |
|
|
6 |
+ |
#include <iostream> |
7 |
+ |
|
8 |
|
class VHbbCandidateTools { |
9 |
|
public: |
10 |
+ |
|
11 |
+ |
VHbbCandidateTools(bool verbose = false): verbose_(verbose){} |
12 |
+ |
|
13 |
|
float deltaPhi(float in2, float in1){ |
14 |
|
float dphi = in2-in1; |
15 |
|
if ( dphi > M_PI ) { |
21 |
|
} |
22 |
|
|
23 |
|
VHbbCandidate getHZmumuCandidate(const VHbbCandidate & in, bool & ok){ |
24 |
+ |
if (verbose_){ |
25 |
+ |
std::cout <<" getHZmumuCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<std::endl; |
26 |
+ |
} |
27 |
|
ok = false; |
28 |
|
VHbbCandidate temp=in; |
29 |
|
if (temp.V.muons.size()!=2) return in ; |
44 |
|
return temp; |
45 |
|
} |
46 |
|
VHbbCandidate getHZeeCandidate(const VHbbCandidate & in, bool & ok){ |
47 |
+ |
if (verbose_){ |
48 |
+ |
std::cout <<" getHZeeCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<std::endl; |
49 |
+ |
} |
50 |
+ |
|
51 |
|
ok = false; |
52 |
|
VHbbCandidate temp=in; |
53 |
|
if (temp.V.electrons.size()!=2) return in ; |
73 |
|
return temp; |
74 |
|
} |
75 |
|
VHbbCandidate getHZnnCandidate(const VHbbCandidate & in, bool & ok){ |
76 |
+ |
if (verbose_){ |
77 |
+ |
std::cout <<" getHZnnCandidate input mu "<<in.V.muons.size()<<" e "<<in.V.electrons.size()<<" met "<<in.V.mets.size()<<std::endl; |
78 |
+ |
} |
79 |
+ |
|
80 |
|
ok = false; |
81 |
|
VHbbCandidate temp=in; |
82 |
|
if (temp.V.mets.size()!=1) return in; |
84 |
|
if (temp.V.electrons.size()!=0) return in ; |
85 |
|
|
86 |
|
temp.V.fourMomentum = temp.V.mets[0].fourMomentum; |
87 |
< |
|
87 |
> |
if (verbose_) { |
88 |
> |
std::cout <<" debug met "<< temp.V.mets[0].metSig << " " << temp.V.mets[0].sumEt<< std::endl; |
89 |
> |
} |
90 |
|
if (temp.V.mets[0].metSig<5) return in; |
91 |
< |
if (temp.V.mets[0].sumEt<150) return in; |
91 |
> |
if (temp.V.mets[0].sumEt<50) return in; |
92 |
|
// if (temp.H.fourMomentum.Pt()<150)return in; |
93 |
|
// if (temp.H.firstJet().csv< 0.9) return in; |
94 |
|
// if (temp.H.secondJet().csv<0.5) return in; |
123 |
|
return in; |
124 |
|
} |
125 |
|
|
126 |
+ |
public: |
127 |
+ |
bool verbose_; |
128 |
+ |
|
129 |
|
}; |
130 |
|
|
131 |
|
|