CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Types | Public Member Functions
res::HelperJet Class Reference

#include <Jet.h>

Public Types

enum  Flavor { kUds, kB }
 

Public Member Functions

double a (double pt, double eta, Flavor flav)
 
double b (double pt, double eta, Flavor flav)
 
double c (double pt, double eta, Flavor flav)
 
double d (double pt, double eta, Flavor flav)
 
double et (double pt, double eta, Flavor flav)
 
double eta (double pt, double eta, Flavor flav)
 
 HelperJet ()
 
double phi (double pt, double eta, Flavor flav)
 
double pt (double pt, double eta, Flavor flav)
 
double theta (double pt, double eta, Flavor flav)
 
 ~HelperJet ()
 

Detailed Description

Definition at line 7 of file Jet.h.

Member Enumeration Documentation

enum res::HelperJet::Flavor
Enumerator
kUds 
kB 

Definition at line 9 of file Jet.h.

9 { kUds, kB };

Constructor & Destructor Documentation

res::HelperJet::HelperJet ( )
inline

Definition at line 11 of file Jet.h.

11 {};
res::HelperJet::~HelperJet ( )
inline

Definition at line 12 of file Jet.h.

12 {};

Member Function Documentation

double res::HelperJet::a ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 42 of file Jet.h.

References funct::exp(), and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix().

42  {
43  double res = 0.29 * sqrt(pt);
44  if (flav == kB) {
45  if (fabs(eta) < 0.17)
46  res = 0.142188 + 0.288157 * exp(-(0.0275628 * pt));
47  else if (fabs(eta) < 0.35)
48  res = 0.141215 + 0.277502 * exp(-(0.0282327 * pt));
49  else if (fabs(eta) < 0.5)
50  res = 0.140493 + 0.26136 * exp(-(0.0254784 * pt));
51  else if (fabs(eta) < 0.7)
52  res = 0.142608 + 0.371256 * exp(-(0.0299132 * pt));
53  else if (fabs(eta) < 0.9)
54  res = 0.138817 + 0.237167 * exp(-(0.0268413 * pt));
55  else if (fabs(eta) < 1.15)
56  res = 0.137322 + 0.293834 * exp(-(0.0264771 * pt));
57  else if (fabs(eta) < 1.4)
58  res = 0.147184 + 0.40834 * exp(-(0.0366769 * pt));
59  else if (fabs(eta) < 1.7)
60  res = 0.141529 + 0.195563 * exp(-(0.023933 * pt));
61  else if (fabs(eta) < 2.1)
62  res = 0.129132 + 0.261012 * exp(-(0.0323658 * pt));
63  else if (fabs(eta) < 2.5)
64  res = 0.120065 + 0.239351 * exp(-(0.0356176 * pt));
65  } else {
66  if (fabs(eta) < 0.17)
67  res = 0.142145 + 0.274541 * exp(-(0.0286487 * pt));
68  else if (fabs(eta) < 0.35)
69  res = 0.146112 + 0.284777 * exp(-(0.031039 * pt));
70  else if (fabs(eta) < 0.5)
71  res = 0.146703 + 0.343982 * exp(-(0.0333917 * pt));
72  else if (fabs(eta) < 0.7)
73  res = 0.149067 + 0.241197 * exp(-(0.0303415 * pt));
74  else if (fabs(eta) < 0.9)
75  res = 0.141233 + 0.304052 * exp(-(0.0314644 * pt));
76  else if (fabs(eta) < 1.15)
77  res = 0.144531 + 0.266871 * exp(-(0.0283385 * pt));
78  else if (fabs(eta) < 1.4)
79  res = 0.156593 + 0.41251 * exp(-(0.0411896 * pt));
80  else if (fabs(eta) < 1.7)
81  res = 0.13452 + 0.269 * exp(-(0.0230171 * pt));
82  else if (fabs(eta) < 2.1)
83  res = 0.0969416 + 0.211636 * exp(-(0.0166117 * pt));
84  else if (fabs(eta) < 2.5)
85  res = 0.0727398 + 0.205132 * exp(-(0.01661 * pt));
86  }
87  return res;
88 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::b ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 90 of file Jet.h.

References alignCSCRings::e, funct::exp(), and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix().

90  {
91  double res = 0.29 * sqrt(pt);
92  if (flav == kB) {
93  if (fabs(eta) < 0.17)
94  res = 2.02521 + 0.174301 * exp(-(-0.0136646 * pt));
95  else if (fabs(eta) < 0.35)
96  res = 0.690872 + 1.32174 * exp(-(-0.00459989 * pt));
97  else if (fabs(eta) < 0.5)
98  res = -956.522 + 958.486 * exp(-(-7.57551e-06 * pt));
99  else if (fabs(eta) < 0.7)
100  res = 2.01017 + 0.0670861 * exp(-(-0.0194879 * pt));
101  else if (fabs(eta) < 0.9)
102  res = 1.92701 + 0.124384 * exp(-(-0.0144246 * pt));
103  else if (fabs(eta) < 1.15)
104  res = 1.94854 + 0.034138 * exp(-(-0.0227995 * pt));
105  else if (fabs(eta) < 1.4)
106  res = 2.07599 + 0.0208406 * exp(-(-0.0214439 * pt));
107  else if (fabs(eta) < 1.7)
108  res = 1.99002 + 0.0077067 * exp(-(-0.0278379 * pt));
109  else if (fabs(eta) < 2.1)
110  res = 1.93181 + 0.00196985 * exp(-(-0.0390155 * pt));
111  else if (fabs(eta) < 2.5)
112  res = 1.75611 + 0.0287099 * exp(-(-0.0202162 * pt));
113  } else {
114  if (fabs(eta) < 0.17)
115  res = -644.18 + 645.342 * exp(-(-1.7462e-05 * pt));
116  else if (fabs(eta) < 0.35)
117  res = -1.34986 + 2.48937 * exp(-(-0.00396448 * pt));
118  else if (fabs(eta) < 0.5)
119  res = -4.8078 + 5.98195 * exp(-(-0.0015 * pt));
120  else if (fabs(eta) < 0.7)
121  res = -9.55949 + 10.7419 * exp(-(-0.000891802 * pt));
122  else if (fabs(eta) < 0.9)
123  res = 0.374264 + 0.884391 * exp(-(-0.00642501 * pt));
124  else if (fabs(eta) < 1.15)
125  res = -9.51561 + 10.7429 * exp(-(-0.000639173 * pt));
126  else if (fabs(eta) < 1.4)
127  res = -0.158886 + 1.43841 * exp(-(-0.00347282 * pt));
128  else if (fabs(eta) < 1.7)
129  res = 0.28876 + 0.992659 * exp(-(-0.00454673 * pt));
130  else if (fabs(eta) < 2.1)
131  res = 0.292201 + 0.922609 * exp(-(-0.00536289 * pt));
132  else if (fabs(eta) < 2.5)
133  res = -657.479 + 658.533 * exp(-(-1.40481e-05 * pt));
134  }
135  return res;
136 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::c ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 138 of file Jet.h.

References alignCSCRings::e, funct::exp(), and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix().

138  {
139  double res = 0.29 * sqrt(pt);
140  if (flav == kB) {
141  if (fabs(eta) < 0.17)
142  res = -2048.11 + 2051.27 * exp(-(1.81337e-06 * pt));
143  else if (fabs(eta) < 0.35)
144  res = 2.86314 + 0.598484 * exp(-(0.0322582 * pt));
145  else if (fabs(eta) < 0.5)
146  res = 2.54785 + 0.77807 * exp(-(0.0100979 * pt));
147  else if (fabs(eta) < 0.7)
148  res = 2.78865 + 6.82558 * exp(-(0.0677953 * pt));
149  else if (fabs(eta) < 0.9)
150  res = 2.63261 + 2.66562 * exp(-(0.044099 * pt));
151  else if (fabs(eta) < 1.15)
152  res = 2.51691 + 2.75456 * exp(-(0.0425201 * pt));
153  else if (fabs(eta) < 1.4)
154  res = 2.60458 + 0.623776 * exp(-(0.0148415 * pt));
155  else if (fabs(eta) < 1.7)
156  res = 2.39005 + 1.60119 * exp(-(0.0154536 * pt));
157  else if (fabs(eta) < 2.1)
158  res = 2.49783 + 0.408896 * exp(-(0.04005 * pt));
159  else if (fabs(eta) < 2.5)
160  res = 2.02143 + 0.000354016 * exp(-(-0.0476218 * pt));
161  } else {
162  if (fabs(eta) < 0.17)
163  res = 1.72736 + 0.212526 * exp(-(-0.00672611 * pt));
164  else if (fabs(eta) < 0.35)
165  res = -3480.32 + 3482.17 * exp(-(-7.08939e-07 * pt));
166  else if (fabs(eta) < 0.5)
167  res = 1.89173 + 0.0381423 * exp(-(-0.0151248 * pt));
168  else if (fabs(eta) < 0.7)
169  res = 2 + 0.00575323 * exp(-(-0.0233459 * pt));
170  else if (fabs(eta) < 0.9)
171  res = 1.92618 + 0.00541047 * exp(-(-0.0269683 * pt));
172  else if (fabs(eta) < 1.15)
173  res = 1.95535 + 0.0194328 * exp(-(-0.0175058 * pt));
174  else if (fabs(eta) < 1.4)
175  res = 1.99855 + 0.00403398 * exp(-(-0.027822 * pt));
176  else if (fabs(eta) < 1.7)
177  res = 2.02356 + 1.33745e-15 * exp(-(-0.201028 * pt));
178  else if (fabs(eta) < 2.1)
179  res = 1.70305 + 0.0428122 * exp(-(-0.0206051 * pt));
180  else if (fabs(eta) < 2.5)
181  res = 0.964128 + 0.370407 * exp(-(-0.0119056 * pt));
182  }
183  return res;
184 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::d ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 186 of file Jet.h.

References funct::exp(), and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix().

186  {
187  double res = 0.29 * sqrt(pt);
188  if (flav == kB) {
189  if (fabs(eta) < 0.17)
190  res = 0.144545 + 0.374307 * exp(-(0.0314642 * pt));
191  else if (fabs(eta) < 0.35)
192  res = 0.141366 + 0.241457 * exp(-(0.0263273 * pt));
193  else if (fabs(eta) < 0.5)
194  res = 0.142509 + 0.324478 * exp(-(0.0288208 * pt));
195  else if (fabs(eta) < 0.7)
196  res = 0.139809 + 0.292392 * exp(-(0.0257797 * pt));
197  else if (fabs(eta) < 0.9)
198  res = 0.137762 + 0.258213 * exp(-(0.0265974 * pt));
199  else if (fabs(eta) < 1.15)
200  res = 0.135831 + 0.338346 * exp(-(0.0277574 * pt));
201  else if (fabs(eta) < 1.4)
202  res = 0.148253 + 0.53999 * exp(-(0.0415757 * pt));
203  else if (fabs(eta) < 1.7)
204  res = 0.145785 + 0.297049 * exp(-(0.0319928 * pt));
205  else if (fabs(eta) < 2.1)
206  res = 0.126895 + 0.291297 * exp(-(0.0328976 * pt));
207  else if (fabs(eta) < 2.5)
208  res = 0.119259 + 0.228403 * exp(-(0.0341018 * pt));
209  } else {
210  if (fabs(eta) < 0.17)
211  res = 0.139738 + 0.310387 * exp(-(0.0296992 * pt));
212  else if (fabs(eta) < 0.35)
213  res = 0.142229 + 0.274556 * exp(-(0.0286575 * pt));
214  else if (fabs(eta) < 0.5)
215  res = 0.14344 + 0.313425 * exp(-(0.0294567 * pt));
216  else if (fabs(eta) < 0.7)
217  res = 0.147284 + 0.252005 * exp(-(0.0303533 * pt));
218  else if (fabs(eta) < 0.9)
219  res = 0.13541 + 0.261376 * exp(-(0.0260026 * pt));
220  else if (fabs(eta) < 1.15)
221  res = 0.144023 + 0.27422 * exp(-(0.0282282 * pt));
222  else if (fabs(eta) < 1.4)
223  res = 0.158694 + 0.352571 * exp(-(0.0385597 * pt));
224  else if (fabs(eta) < 1.7)
225  res = 0.135936 + 0.24773 * exp(-(0.0226282 * pt));
226  else if (fabs(eta) < 2.1)
227  res = 0.0957325 + 0.213108 * exp(-(0.016424 * pt));
228  else if (fabs(eta) < 2.5)
229  res = 0.0622281 + 0.185045 * exp(-(0.0129321 * pt));
230  }
231  return res;
232 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::et ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 330 of file Jet.h.

References DiDispStaMuonMonitor_cfi::pt, and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix().

330  {
331  double res = 0.29 * sqrt(pt);
332  if (flav == kB) {
333  if (fabs(eta) < 0.17)
334  res = 10.1445 + 0.0445196 * pt;
335  else if (fabs(eta) < 0.35)
336  res = 8.944 + 0.0576391 * pt;
337  else if (fabs(eta) < 0.5)
338  res = 8.8462 + 0.0630703 * pt;
339  else if (fabs(eta) < 0.7)
340  res = 8.90812 + 0.057605 * pt;
341  else if (fabs(eta) < 0.9)
342  res = 8.42307 + 0.0623024 * pt;
343  else if (fabs(eta) < 1.15)
344  res = 8.64735 + 0.0619689 * pt;
345  else if (fabs(eta) < 1.4)
346  res = 8.90574 + 0.0595526 * pt;
347  else if (fabs(eta) < 1.7)
348  res = 8.74415 + 0.0612412 * pt;
349  else if (fabs(eta) < 2.1)
350  res = 7.80525 + 0.0511843 * pt;
351  else if (fabs(eta) < 2.5)
352  res = 6.07877 + 0.0569574 * pt;
353  } else {
354  if (fabs(eta) < 0.17)
355  res = 8.04345 + 0.0705188 * pt;
356  else if (fabs(eta) < 0.35)
357  res = 7.77565 + 0.0722247 * pt;
358  else if (fabs(eta) < 0.5)
359  res = 8.53558 + 0.0647725 * pt;
360  else if (fabs(eta) < 0.7)
361  res = 7.99374 + 0.0681428 * pt;
362  else if (fabs(eta) < 0.9)
363  res = 7.51695 + 0.0755128 * pt;
364  else if (fabs(eta) < 1.15)
365  res = 7.70804 + 0.0748981 * pt;
366  else if (fabs(eta) < 1.4)
367  res = 7.2996 + 0.0852397 * pt;
368  else if (fabs(eta) < 1.7)
369  res = 7.95923 + 0.0715722 * pt;
370  else if (fabs(eta) < 2.1)
371  res = 7.68012 + 0.0543317 * pt;
372  else if (fabs(eta) < 2.5)
373  res = 6.59141 + 0.0452384 * pt;
374  }
375  return res;
376 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::eta ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 378 of file Jet.h.

References funct::exp().

Referenced by Particle.Particle::__str__(), StKinFitter::addKinFitInfo(), Jet.Jet::jetID(), Jet.Jet::puJetId(), and CovarianceMatrix::setupMatrix().

378  {
379  double res = -1.53e-4 * pt + 0.05;
380  if (flav == kB) {
381  if (fabs(eta) < 0.17)
382  res = 0.0354411 + 0.130211 * exp(-(0.0423321 * pt));
383  else if (fabs(eta) < 0.35)
384  res = 0.0341742 + 0.115784 * exp(-(0.0375787 * pt));
385  else if (fabs(eta) < 0.5)
386  res = 0.0336537 + 0.143109 * exp(-(0.0416702 * pt));
387  else if (fabs(eta) < 0.7)
388  res = 0.031894 + 0.140572 * exp(-(0.0405564 * pt));
389  else if (fabs(eta) < 0.9)
390  res = 0.0301507 + 0.126005 * exp(-(0.0384231 * pt));
391  else if (fabs(eta) < 1.15)
392  res = 0.0268018 + 0.140956 * exp(-(0.03872 * pt));
393  else if (fabs(eta) < 1.4)
394  res = 0.0242941 + 0.131887 * exp(-(0.0337869 * pt));
395  else if (fabs(eta) < 1.7)
396  res = 0.0223985 + 0.142949 * exp(-(0.0351412 * pt));
397  else if (fabs(eta) < 2.1)
398  res = 0.0229694 + 0.145993 * exp(-(0.0391759 * pt));
399  else if (fabs(eta) < 2.5)
400  res = 0.0243171 + 0.145419 * exp(-(0.0430824 * pt));
401  } else {
402  if (fabs(eta) < 0.17)
403  res = 0.038575 + 0.129026 * exp(-(0.0427157 * pt));
404  else if (fabs(eta) < 0.35)
405  res = 0.0400221 + 0.127222 * exp(-(0.0458212 * pt));
406  else if (fabs(eta) < 0.5)
407  res = 0.0358472 + 0.119962 * exp(-(0.0395586 * pt));
408  else if (fabs(eta) < 0.7)
409  res = 0.0359612 + 0.136094 * exp(-(0.04212 * pt));
410  else if (fabs(eta) < 0.9)
411  res = 0.0327581 + 0.127339 * exp(-(0.0393077 * pt));
412  else if (fabs(eta) < 1.15)
413  res = 0.0305095 + 0.116238 * exp(-(0.0347537 * pt));
414  else if (fabs(eta) < 1.4)
415  res = 0.027529 + 0.117893 * exp(-(0.032858 * pt));
416  else if (fabs(eta) < 1.7)
417  res = 0.0253065 + 0.112751 * exp(-(0.0307028 * pt));
418  else if (fabs(eta) < 2.1)
419  res = 0.0266957 + 0.131998 * exp(-(0.0380667 * pt));
420  else if (fabs(eta) < 2.5)
421  res = 0.0298413 + 0.13651 * exp(-(0.0458475 * pt));
422  }
423  return res;
424 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::phi ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 282 of file Jet.h.

References funct::exp(), and mathSSE::sqrt().

Referenced by Particle.Particle::__str__(), StKinFitter::addKinFitInfo(), ntupleDataFormat.Track::phiPull(), and CovarianceMatrix::setupMatrix().

282  {
283  double res = 0.29 * sqrt(pt);
284  if (flav == kB) {
285  if (fabs(eta) < 0.17)
286  res = 0.0217379 + 0.157228 * exp(-(0.0237049 * pt));
287  else if (fabs(eta) < 0.35)
288  res = 0.0246735 + 0.165694 * exp(-(0.0269728 * pt));
289  else if (fabs(eta) < 0.5)
290  res = 0.0226496 + 0.162875 * exp(-(0.0250225 * pt));
291  else if (fabs(eta) < 0.7)
292  res = 0.0232004 + 0.167887 * exp(-(0.0259427 * pt));
293  else if (fabs(eta) < 0.9)
294  res = 0.0220083 + 0.183889 * exp(-(0.0273017 * pt));
295  else if (fabs(eta) < 1.15)
296  res = 0.0216601 + 0.18482 * exp(-(0.027148 * pt));
297  else if (fabs(eta) < 1.4)
298  res = 0.0225191 + 0.19558 * exp(-(0.0271747 * pt));
299  else if (fabs(eta) < 1.7)
300  res = 0.0208365 + 0.19241 * exp(-(0.0262627 * pt));
301  else if (fabs(eta) < 2.1)
302  res = 0.0233749 + 0.180684 * exp(-(0.032052 * pt));
303  else if (fabs(eta) < 2.5)
304  res = 0.0231891 + 0.165484 * exp(-(0.0389916 * pt));
305  } else {
306  if (fabs(eta) < 0.17)
307  res = 0.025945 + 0.16529 * exp(-(0.027385 * pt));
308  else if (fabs(eta) < 0.35)
309  res = 0.0257294 + 0.155366 * exp(-(0.0268639 * pt));
310  else if (fabs(eta) < 0.5)
311  res = 0.0231924 + 0.166508 * exp(-(0.0264153 * pt));
312  else if (fabs(eta) < 0.7)
313  res = 0.0252736 + 0.176977 * exp(-(0.0281594 * pt));
314  else if (fabs(eta) < 0.9)
315  res = 0.0269299 + 0.174505 * exp(-(0.0295757 * pt));
316  else if (fabs(eta) < 1.15)
317  res = 0.0250191 + 0.178321 * exp(-(0.0284006 * pt));
318  else if (fabs(eta) < 1.4)
319  res = 0.0242345 + 0.174682 * exp(-(0.027202 * pt));
320  else if (fabs(eta) < 1.7)
321  res = 0.0242234 + 0.199735 * exp(-(0.0291262 * pt));
322  else if (fabs(eta) < 2.1)
323  res = 0.0300017 + 0.175617 * exp(-(0.0357093 * pt));
324  else if (fabs(eta) < 2.5)
325  res = 0.0282978 + 0.141994 * exp(-(0.0383663 * pt));
326  }
327  return res;
328 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::pt ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 26 of file Jet.h.

References kB, and mathSSE::sqrt().

Referenced by Particle.Particle::__str__(), DiObject.DiMuon::__str__(), StKinFitter::addKinFitInfo(), Photon.Photon::calScaledIsoValueExp(), Photon.Photon::calScaledIsoValueLin(), Photon.Photon::calScaledIsoValueQuadr(), Electron.Electron::mvaIDLoose(), Electron.Electron::mvaIDRun2(), Electron.Electron::mvaIDTight(), Electron.Electron::ptErr(), ntupleDataFormat.Track::ptPull(), Lepton.Lepton::relIso(), Lepton.Lepton::relIsoFromEA(), Lepton.Lepton::relIsoR(), and Jet.Jet::setCorrP4().

26  {
27  double res = 0.29 * sqrt(pt);
28  if (fabs(eta) < 1.4)
29  res += 6.68;
30  else
31  res += -3.14 * fabs(eta) + 11.89;
32  if (flav == kB) {
33  res = 0.33 * sqrt(pt);
34  if (fabs(eta) < 1.4)
35  res += 6.57;
36  else
37  res += -1.09 * fabs(eta) + 8.50;
38  }
39  return res;
40 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
double res::HelperJet::theta ( double  pt,
double  eta,
Flavor  flav 
)
inline

Definition at line 234 of file Jet.h.

References funct::exp(), and mathSSE::sqrt().

Referenced by CovarianceMatrix::setupMatrix(), and Tau.Tau::zImpact().

234  {
235  double res = 0.29 * sqrt(pt);
236  if (flav == kB) {
237  if (fabs(eta) < 0.17)
238  res = 0.0348657 + 0.0966067 * exp(-(0.0388496 * pt));
239  else if (fabs(eta) < 0.35)
240  res = 0.0327521 + 0.0901303 * exp(-(0.0354149 * pt));
241  else if (fabs(eta) < 0.5)
242  res = 0.030283 + 0.0970732 * exp(-(0.036806 * pt));
243  else if (fabs(eta) < 0.7)
244  res = 0.0266898 + 0.0994089 * exp(-(0.0379773 * pt));
245  else if (fabs(eta) < 0.9)
246  res = 0.022614 + 0.091974 * exp(-(0.0384986 * pt));
247  else if (fabs(eta) < 1.15)
248  res = 0.0173219 + 0.0923928 * exp(-(0.0395886 * pt));
249  else if (fabs(eta) < 1.4)
250  res = 0.0129851 + 0.077473 * exp(-(0.0363236 * pt));
251  else if (fabs(eta) < 1.7)
252  res = 0.00967747 + 0.0667463 * exp(-(0.0383101 * pt));
253  else if (fabs(eta) < 2.1)
254  res = 0.00744593 + 0.050027 * exp(-(0.0432323 * pt));
255  else if (fabs(eta) < 2.5)
256  res = 0.00589309 + 0.0346904 * exp(-(0.0484705 * pt));
257  } else {
258  if (fabs(eta) < 0.17)
259  res = 0.0384999 + 0.105382 * exp(-(0.0440247 * pt));
260  else if (fabs(eta) < 0.35)
261  res = 0.0384758 + 0.086093 * exp(-(0.0436464 * pt));
262  else if (fabs(eta) < 0.5)
263  res = 0.0325808 + 0.0860212 * exp(-(0.0371074 * pt));
264  else if (fabs(eta) < 0.7)
265  res = 0.0304022 + 0.094116 * exp(-(0.0402575 * pt));
266  else if (fabs(eta) < 0.9)
267  res = 0.0246649 + 0.0883167 * exp(-(0.0386519 * pt));
268  else if (fabs(eta) < 1.15)
269  res = 0.0197668 + 0.080217 * exp(-(0.0368397 * pt));
270  else if (fabs(eta) < 1.4)
271  res = 0.0148531 + 0.0692569 * exp(-(0.0354383 * pt));
272  else if (fabs(eta) < 1.7)
273  res = 0.0113007 + 0.0601861 * exp(-(0.0373138 * pt));
274  else if (fabs(eta) < 2.1)
275  res = 0.00860103 + 0.0522148 * exp(-(0.0450405 * pt));
276  else if (fabs(eta) < 2.5)
277  res = 0.00692992 + 0.031671 * exp(-(0.0490789 * pt));
278  }
279  return res;
280 }
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
funct::exp
Exp< T >::type exp(const T &t)
Definition: Exp.h:22
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
Generated for CMSSW Reference Manual by doxygen 1.8.5