82 |
|
double |
83 |
|
ElectronSFWeight::at (const double &eta, const double &pt) |
84 |
|
{ |
85 |
< |
cout << eta << " " << pt << endl; |
86 |
< |
cout << "doing it\n"; |
85 |
> |
//cout << eta << " " << pt << endl; |
86 |
> |
//cout << "doing it\n"; |
87 |
|
if (cmsswRelease_ == "53X") |
88 |
|
{ |
89 |
+ |
if (id_ == "loose") |
90 |
+ |
{ |
91 |
+ |
if (fabs (eta) < 0.8) |
92 |
+ |
{ |
93 |
+ |
if (pt < 15.0 && pt >= 10.0) |
94 |
+ |
return 0.855; |
95 |
+ |
if (pt < 20.0 && pt >= 15.0) |
96 |
+ |
return 0.962; |
97 |
+ |
if (pt < 30.0 && pt >= 20.0) |
98 |
+ |
return 1.005; |
99 |
+ |
if (pt < 40.0 && pt >= 30.0) |
100 |
+ |
return 1.004; |
101 |
+ |
if (pt < 50.0 && pt >= 40.0) |
102 |
+ |
return 1.008; |
103 |
+ |
if (pt < 200.0 && pt >= 50.0) |
104 |
+ |
return 1.008; |
105 |
+ |
} |
106 |
+ |
if (fabs (eta) < 1.442 && fabs (eta) >= 0.8) |
107 |
+ |
{ |
108 |
+ |
if (pt < 15.0 && pt >= 10.0) |
109 |
+ |
return 0.858; |
110 |
+ |
if (pt < 20.0 && pt >= 15.0) |
111 |
+ |
return 0.962; |
112 |
+ |
if (pt < 30.0 && pt >= 20.0) |
113 |
+ |
return 0.981; |
114 |
+ |
if (pt < 40.0 && pt >= 30.0) |
115 |
+ |
return 0.991; |
116 |
+ |
if (pt < 50.0 && pt >= 40.0) |
117 |
+ |
return 0.994; |
118 |
+ |
if (pt < 200.0 && pt >= 50.0) |
119 |
+ |
return 0.999; |
120 |
+ |
} |
121 |
+ |
if (fabs (eta) < 1.556 && fabs (eta) >= 1.442) |
122 |
+ |
{ |
123 |
+ |
if (pt < 15.0 && pt >= 10.0) |
124 |
+ |
return 1.109; |
125 |
+ |
if (pt < 20.0 && pt >= 15.0) |
126 |
+ |
return 0.903; |
127 |
+ |
if (pt < 30.0 && pt >= 20.0) |
128 |
+ |
return 1.044; |
129 |
+ |
if (pt < 40.0 && pt >= 30.0) |
130 |
+ |
return 0.998; |
131 |
+ |
if (pt < 50.0 && pt >= 40.0) |
132 |
+ |
return 0.989; |
133 |
+ |
if (pt < 200.0 && pt >= 50.0) |
134 |
+ |
return 0.994; |
135 |
+ |
} |
136 |
+ |
if (fabs (eta) < 2.0 && fabs (eta) >= 1.556) |
137 |
+ |
{ |
138 |
+ |
if (pt < 15.0 && pt >= 10.0) |
139 |
+ |
return 0.838; |
140 |
+ |
if (pt < 20.0 && pt >= 15.0) |
141 |
+ |
return 0.939; |
142 |
+ |
if (pt < 30.0 && pt >= 20.0) |
143 |
+ |
return 0.980; |
144 |
+ |
if (pt < 40.0 && pt >= 30.0) |
145 |
+ |
return 0.992; |
146 |
+ |
if (pt < 50.0 && pt >= 40.0) |
147 |
+ |
return 1.004; |
148 |
+ |
if (pt < 200.0 && pt >= 50.0) |
149 |
+ |
return 1.006; |
150 |
+ |
} |
151 |
+ |
if (fabs (eta) < 2.5 && fabs (eta) >= 2.0) |
152 |
+ |
{ |
153 |
+ |
if (pt < 15.0 && pt >= 10.0) |
154 |
+ |
return 1.034; |
155 |
+ |
if (pt < 20.0 && pt >= 15.0) |
156 |
+ |
return 0.970; |
157 |
+ |
if (pt < 30.0 && pt >= 20.0) |
158 |
+ |
return 1.017; |
159 |
+ |
if (pt < 40.0 && pt >= 30.0) |
160 |
+ |
return 1.019; |
161 |
+ |
if (pt < 50.0 && pt >= 40.0) |
162 |
+ |
return 1.005; |
163 |
+ |
if (pt < 200.0 && pt >= 50.0) |
164 |
+ |
return 1.009; |
165 |
+ |
} |
166 |
+ |
} |
167 |
|
if (id_ == "tight") |
168 |
|
{ |
169 |
|
if (fabs (eta) < 0.8) |
244 |
|
} |
245 |
|
if (id_ == "mvaNonTrig") |
246 |
|
{ |
247 |
< |
cout << "mvaNonTrig\n"; |
247 |
> |
//cout << "mvaNonTrig\n"; |
248 |
|
if (fabs (eta) < 0.8) |
249 |
|
{ |
250 |
|
if (pt < 20.0) |