|
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
|
#include "DQMOffline/Trigger/interface/EgHLTDQMCut.h"
using namespace egHLT;
bool EgTrigTagProbeCut::pass(const OffEle& theEle, const OffEvt& evt) const {
//first we check if our electron passes our id
if (((theEle.*cutCodeFunc_)() & cutCode_) != 0x0)
return false;
//new we check that there is another tag in the event (this electron may be a tag, we are not going to test this, all we care about is that another electron in the event is a tag)
int nrTags = 0;
const OffEle* tagEle = nullptr;
const std::vector<OffEle>& eles = evt.eles();
//we are looking for an *additional* tag
for (auto const& ele : eles) {
if (((ele.*cutCodeFunc_)() & cutCode_) == 0x0 && (bitsToPass_ & ele.trigBits()) == bitsToPass_) {
//now a check that the tag is not the same as the probe
if (reco::deltaR2(theEle.eta(), theEle.phi(), ele.eta(), ele.phi()) >
0.1 * 0.1) { //not in a cone of 0.1 of probe electron
nrTags++;
tagEle = &ele;
}
}
}
if (nrTags ==
1) { //we are requiring one and only one additional tag (the theEle is automatically excluded from the tag list)
float mass = (theEle.p4() + tagEle->p4()).mag();
if (mass > minMass_ && mass < maxMass_)
return true; //mass requirements
}
return false;
}
bool EgTrigTagProbeCut_New::pass(const OffEle& theEle, const OffEvt& evt) const {
//looking at only Et>20, since probe is required to pass SC17
//if(theEle.et()>20/* && theEle.et()<20*/){
//first we check if our probe electron passes WP80 and the second leg of our T&P trigger
if (((theEle.*cutCodeFunc_)() & cutCode_) != 0x0 || (bit2ToPass_ & theEle.trigBits()) != bit2ToPass_)
return false;
//now we check that there is a WP80 tag electron that passes the first leg of the trigger(this electron may be a tag, we are not going to test this, all we care about is that another electron in the event is a tag)
int nrTags = 0;
const OffEle* tagEle = nullptr;
const std::vector<OffEle>& eles = evt.eles();
//we are looking for an *additional* tag
for (auto const& ele : eles) {
if (((ele.*cutCodeFunc_)() & cutCode_) == 0x0 && (bit1ToPass_ & ele.trigBits()) == bit1ToPass_) {
//now a check that the tag is not the same as the probe
if (reco::deltaR2(theEle.eta(), theEle.phi(), ele.eta(), ele.phi()) >
0.1 * 0.1) { //not in a cone of 0.1 of probe electron
nrTags++;
tagEle = &ele;
}
}
}
if (nrTags ==
1) { //we are requiring one and only one additional tag (the theEle is automatically excluded from the tag list)
float mass = (theEle.p4() + tagEle->p4()).mag();
if (mass > minMass_ && mass < maxMass_)
return true; //mass requirements
}
//}//if 10<pt<20
return false;
}
//same for photons
bool EgTrigTagProbeCut_NewPho::pass(const OffPho& thePho, const OffEvt& evt) const {
//looking at only Et>20, since probe is required to pass SC17
//if(theEle.et()>20/* && theEle.et()<20*/){
//first we check if our probe electron passes WP80 and the second leg of our T&P trigger
if (((thePho.*cutCodeFunc_)() & cutCode_) != 0x0 || (bit2ToPass_ & thePho.trigBits()) != bit2ToPass_)
return false;
//now we check that there is a WP80 tag electron that passes the first leg of the trigger(this electron may be a tag, we are not going to test this, all we care about is that another electron in the event is a tag)
int nrTags = 0;
const OffPho* tagPho = nullptr;
const std::vector<OffPho>& phos = evt.phos();
//we are looking for an *additional* tag
for (auto const& pho : phos) {
if (((pho.*cutCodeFunc_)() & cutCode_) == 0x0 && (bit1ToPass_ & pho.trigBits()) == bit1ToPass_) {
//now a check that the tag is not the same as the probe
if (reco::deltaR2(thePho.eta(), thePho.phi(), pho.eta(), pho.phi()) >
0.1 * 0.1) { //not in a cone of 0.1 of probe "photon"
nrTags++;
tagPho = &pho;
}
}
}
if (nrTags ==
1) { //we are requiring one and only one additional tag (the thePho is automatically excluded from the tag list)
float mass = (thePho.p4() + tagPho->p4()).mag();
if (mass > minMass_ && mass < maxMass_)
return true; //mass requirements
}
//}//if 10<pt<20
return false;
}
bool EgDiEleCut::pass(const OffEle& obj, const OffEvt& evt) const {
const std::vector<OffEle>& eles = evt.eles();
for (auto const& ele : eles) {
if (&ele != &obj) { //different electrons
int diEleCutCode = (obj.*cutCodeFunc_)() | (ele.*cutCodeFunc_)();
if ((diEleCutCode & cutCode_) == 0x0)
return true;
}
}
return false;
}
bool EgDiPhoCut::pass(const OffPho& obj, const OffEvt& evt) const {
const std::vector<OffPho>& phos = evt.phos();
for (auto const& pho : phos) {
if (&pho != &obj) { //different phos
int diPhoCutCode = (obj.*cutCodeFunc_)() | (pho.*cutCodeFunc_)();
if ((diPhoCutCode & cutCode_) == 0x0)
return true;
}
}
return false;
}
|