CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Jet.h
Go to the documentation of this file.
1 #ifndef TopObjetcResolutionsJet_h
2 #define TopObjetcResolutionsJet_h
3 
4 namespace res{
5  class HelperJet {
6 
7  public:
8  enum Flavor {kUds, kB};
9 
10  HelperJet(){};
11  ~HelperJet(){};
12 
13  inline double pt (double pt, double eta, Flavor flav);
14  inline double eta(double pt, double eta, Flavor flav);
15  inline double phi(double pt, double eta, Flavor flav);
16  inline double a(double pt, double eta, Flavor flav);
17  inline double b(double pt, double eta, Flavor flav);
18  inline double c(double pt, double eta, Flavor flav);
19  inline double d(double pt, double eta, Flavor flav);
20  inline double et(double pt, double eta, Flavor flav);
21  inline double theta(double pt, double eta, Flavor flav);
22  };
23 }
24 
25 inline double res::HelperJet::pt(double pt, double eta, Flavor flav)
26 {
27  double res = 0.29*sqrt(pt);
28  if(fabs(eta)<1.4) res+= 6.68;
29  else res+=-3.14*fabs(eta)+11.89;
30  if(flav==kB){
31  res=0.33*sqrt(pt);
32  if(fabs(eta)<1.4) res+= 6.57;
33  else res+=-1.09*fabs(eta)+8.50;
34  }
35  return res;
36 }
37 
38 inline double res::HelperJet::a(double pt, double eta, Flavor flav)
39 {
40  double res = 0.29*sqrt(pt);
41  if(flav==kB){
42  if(fabs(eta)<0.17) res = 0.142188+0.288157*exp(-(0.0275628*pt));
43  else if(fabs(eta)<0.35) res = 0.141215+0.277502*exp(-(0.0282327*pt));
44  else if(fabs(eta)<0.5) res = 0.140493+0.26136*exp(-(0.0254784*pt));
45  else if(fabs(eta)<0.7) res = 0.142608+0.371256*exp(-(0.0299132*pt));
46  else if(fabs(eta)<0.9) res = 0.138817+0.237167*exp(-(0.0268413*pt));
47  else if(fabs(eta)<1.15) res = 0.137322+0.293834*exp(-(0.0264771*pt));
48  else if(fabs(eta)<1.4) res = 0.147184+0.40834*exp(-(0.0366769*pt));
49  else if(fabs(eta)<1.7) res = 0.141529+0.195563*exp(-(0.023933*pt));
50  else if(fabs(eta)<2.1) res = 0.129132+0.261012*exp(-(0.0323658*pt));
51  else if(fabs(eta)<2.5) res = 0.120065+0.239351*exp(-(0.0356176*pt));
52  }else{
53  if(fabs(eta)<0.17) res = 0.142145+0.274541*exp(-(0.0286487*pt));
54  else if(fabs(eta)<0.35) res = 0.146112+0.284777*exp(-(0.031039*pt));
55  else if(fabs(eta)<0.5) res = 0.146703+0.343982*exp(-(0.0333917*pt));
56  else if(fabs(eta)<0.7) res = 0.149067+0.241197*exp(-(0.0303415*pt));
57  else if(fabs(eta)<0.9) res = 0.141233+0.304052*exp(-(0.0314644*pt));
58  else if(fabs(eta)<1.15) res = 0.144531+0.266871*exp(-(0.0283385*pt));
59  else if(fabs(eta)<1.4) res = 0.156593+0.41251*exp(-(0.0411896*pt));
60  else if(fabs(eta)<1.7) res = 0.13452+0.269*exp(-(0.0230171*pt));
61  else if(fabs(eta)<2.1) res = 0.0969416+0.211636*exp(-(0.0166117*pt));
62  else if(fabs(eta)<2.5) res = 0.0727398+0.205132*exp(-(0.01661*pt));
63  }
64  return res;
65 }
66 
67 inline double res::HelperJet::b(double pt, double eta, Flavor flav)
68 {
69  double res = 0.29*sqrt(pt);
70  if(flav==kB){
71  if(fabs(eta)<0.17) res = 2.02521+0.174301*exp(-(-0.0136646*pt));
72  else if(fabs(eta)<0.35) res = 0.690872+1.32174*exp(-(-0.00459989*pt));
73  else if(fabs(eta)<0.5) res = -956.522+958.486*exp(-(-7.57551e-06*pt));
74  else if(fabs(eta)<0.7) res = 2.01017+0.0670861*exp(-(-0.0194879*pt));
75  else if(fabs(eta)<0.9) res = 1.92701+0.124384*exp(-(-0.0144246*pt));
76  else if(fabs(eta)<1.15) res = 1.94854+0.034138*exp(-(-0.0227995*pt));
77  else if(fabs(eta)<1.4) res = 2.07599+0.0208406*exp(-(-0.0214439*pt));
78  else if(fabs(eta)<1.7) res = 1.99002+0.0077067*exp(-(-0.0278379*pt));
79  else if(fabs(eta)<2.1) res = 1.93181+0.00196985*exp(-(-0.0390155*pt));
80  else if(fabs(eta)<2.5) res = 1.75611+0.0287099*exp(-(-0.0202162*pt));
81  }else{
82  if(fabs(eta)<0.17) res = -644.18+645.342*exp(-(-1.7462e-05*pt));
83  else if(fabs(eta)<0.35) res = -1.34986+2.48937*exp(-(-0.00396448*pt));
84  else if(fabs(eta)<0.5) res = -4.8078+5.98195*exp(-(-0.0015*pt));
85  else if(fabs(eta)<0.7) res = -9.55949+10.7419*exp(-(-0.000891802*pt));
86  else if(fabs(eta)<0.9) res = 0.374264+0.884391*exp(-(-0.00642501*pt));
87  else if(fabs(eta)<1.15) res = -9.51561+10.7429*exp(-(-0.000639173*pt));
88  else if(fabs(eta)<1.4) res = -0.158886+1.43841*exp(-(-0.00347282*pt));
89  else if(fabs(eta)<1.7) res = 0.28876+0.992659*exp(-(-0.00454673*pt));
90  else if(fabs(eta)<2.1) res = 0.292201+0.922609*exp(-(-0.00536289*pt));
91  else if(fabs(eta)<2.5) res = -657.479+658.533*exp(-(-1.40481e-05*pt));
92  }
93  return res;
94 }
95 
96 inline double res::HelperJet::c(double pt, double eta, Flavor flav)
97 {
98  double res = 0.29*sqrt(pt);
99  if(flav==kB){
100  if(fabs(eta)<0.17) res = -2048.11+2051.27*exp(-(1.81337e-06*pt));
101  else if(fabs(eta)<0.35) res = 2.86314+0.598484*exp(-(0.0322582*pt));
102  else if(fabs(eta)<0.5) res = 2.54785+0.77807*exp(-(0.0100979*pt));
103  else if(fabs(eta)<0.7) res = 2.78865+6.82558*exp(-(0.0677953*pt));
104  else if(fabs(eta)<0.9) res = 2.63261+2.66562*exp(-(0.044099*pt));
105  else if(fabs(eta)<1.15) res = 2.51691+2.75456*exp(-(0.0425201*pt));
106  else if(fabs(eta)<1.4) res = 2.60458+0.623776*exp(-(0.0148415*pt));
107  else if(fabs(eta)<1.7) res = 2.39005+1.60119*exp(-(0.0154536*pt));
108  else if(fabs(eta)<2.1) res = 2.49783+0.408896*exp(-(0.04005*pt));
109  else if(fabs(eta)<2.5) res = 2.02143+0.000354016*exp(-(-0.0476218*pt));
110  }else{
111  if(fabs(eta)<0.17) res = 1.72736+0.212526*exp(-(-0.00672611*pt));
112  else if(fabs(eta)<0.35) res = -3480.32+3482.17*exp(-(-7.08939e-07*pt));
113  else if(fabs(eta)<0.5) res = 1.89173+0.0381423*exp(-(-0.0151248*pt));
114  else if(fabs(eta)<0.7) res = 2+0.00575323*exp(-(-0.0233459*pt));
115  else if(fabs(eta)<0.9) res = 1.92618+0.00541047*exp(-(-0.0269683*pt));
116  else if(fabs(eta)<1.15) res = 1.95535+0.0194328*exp(-(-0.0175058*pt));
117  else if(fabs(eta)<1.4) res = 1.99855+0.00403398*exp(-(-0.027822*pt));
118  else if(fabs(eta)<1.7) res = 2.02356+1.33745e-15*exp(-(-0.201028*pt));
119  else if(fabs(eta)<2.1) res = 1.70305+0.0428122*exp(-(-0.0206051*pt));
120  else if(fabs(eta)<2.5) res = 0.964128+0.370407*exp(-(-0.0119056*pt));
121  }
122  return res;
123 }
124 
125 inline double res::HelperJet::d(double pt, double eta, Flavor flav)
126 {
127  double res = 0.29*sqrt(pt);
128  if(flav==kB){
129  if(fabs(eta)<0.17) res = 0.144545+0.374307*exp(-(0.0314642*pt));
130  else if(fabs(eta)<0.35) res = 0.141366+0.241457*exp(-(0.0263273*pt));
131  else if(fabs(eta)<0.5) res = 0.142509+0.324478*exp(-(0.0288208*pt));
132  else if(fabs(eta)<0.7) res = 0.139809+0.292392*exp(-(0.0257797*pt));
133  else if(fabs(eta)<0.9) res = 0.137762+0.258213*exp(-(0.0265974*pt));
134  else if(fabs(eta)<1.15) res = 0.135831+0.338346*exp(-(0.0277574*pt));
135  else if(fabs(eta)<1.4) res = 0.148253+0.53999*exp(-(0.0415757*pt));
136  else if(fabs(eta)<1.7) res = 0.145785+0.297049*exp(-(0.0319928*pt));
137  else if(fabs(eta)<2.1) res = 0.126895+0.291297*exp(-(0.0328976*pt));
138  else if(fabs(eta)<2.5) res = 0.119259+0.228403*exp(-(0.0341018*pt));
139  }else{
140  if(fabs(eta)<0.17) res = 0.139738+0.310387*exp(-(0.0296992*pt));
141  else if(fabs(eta)<0.35) res = 0.142229+0.274556*exp(-(0.0286575*pt));
142  else if(fabs(eta)<0.5) res = 0.14344+0.313425*exp(-(0.0294567*pt));
143  else if(fabs(eta)<0.7) res = 0.147284+0.252005*exp(-(0.0303533*pt));
144  else if(fabs(eta)<0.9) res = 0.13541+0.261376*exp(-(0.0260026*pt));
145  else if(fabs(eta)<1.15) res = 0.144023+0.27422*exp(-(0.0282282*pt));
146  else if(fabs(eta)<1.4) res = 0.158694+0.352571*exp(-(0.0385597*pt));
147  else if(fabs(eta)<1.7) res = 0.135936+0.24773*exp(-(0.0226282*pt));
148  else if(fabs(eta)<2.1) res = 0.0957325+0.213108*exp(-(0.016424*pt));
149  else if(fabs(eta)<2.5) res = 0.0622281+0.185045*exp(-(0.0129321*pt));
150  }
151  return res;
152 }
153 
154 inline double res::HelperJet::theta(double pt, double eta, Flavor flav)
155 {
156  double res = 0.29*sqrt(pt);
157  if(flav==kB){
158  if(fabs(eta)<0.17) res = 0.0348657+0.0966067*exp(-(0.0388496*pt));
159  else if(fabs(eta)<0.35) res = 0.0327521+0.0901303*exp(-(0.0354149*pt));
160  else if(fabs(eta)<0.5) res = 0.030283+0.0970732*exp(-(0.036806*pt));
161  else if(fabs(eta)<0.7) res = 0.0266898+0.0994089*exp(-(0.0379773*pt));
162  else if(fabs(eta)<0.9) res = 0.022614+0.091974*exp(-(0.0384986*pt));
163  else if(fabs(eta)<1.15) res = 0.0173219+0.0923928*exp(-(0.0395886*pt));
164  else if(fabs(eta)<1.4) res = 0.0129851+0.077473*exp(-(0.0363236*pt));
165  else if(fabs(eta)<1.7) res = 0.00967747+0.0667463*exp(-(0.0383101*pt));
166  else if(fabs(eta)<2.1) res = 0.00744593+0.050027*exp(-(0.0432323*pt));
167  else if(fabs(eta)<2.5) res = 0.00589309+0.0346904*exp(-(0.0484705*pt));
168  }else{
169  if(fabs(eta)<0.17) res = 0.0384999+0.105382*exp(-(0.0440247*pt));
170  else if(fabs(eta)<0.35) res = 0.0384758+0.086093*exp(-(0.0436464*pt));
171  else if(fabs(eta)<0.5) res = 0.0325808+0.0860212*exp(-(0.0371074*pt));
172  else if(fabs(eta)<0.7) res = 0.0304022+0.094116*exp(-(0.0402575*pt));
173  else if(fabs(eta)<0.9) res = 0.0246649+0.0883167*exp(-(0.0386519*pt));
174  else if(fabs(eta)<1.15) res = 0.0197668+0.080217*exp(-(0.0368397*pt));
175  else if(fabs(eta)<1.4) res = 0.0148531+0.0692569*exp(-(0.0354383*pt));
176  else if(fabs(eta)<1.7) res = 0.0113007+0.0601861*exp(-(0.0373138*pt));
177  else if(fabs(eta)<2.1) res = 0.00860103+0.0522148*exp(-(0.0450405*pt));
178  else if(fabs(eta)<2.5) res = 0.00692992+0.031671*exp(-(0.0490789*pt));
179  }
180  return res;
181 }
182 
183 inline double res::HelperJet::phi(double pt, double eta, Flavor flav)
184 {
185  double res = 0.29*sqrt(pt);
186  if(flav==kB){
187  if(fabs(eta)<0.17) res = 0.0217379+0.157228*exp(-(0.0237049*pt));
188  else if(fabs(eta)<0.35) res = 0.0246735+0.165694*exp(-(0.0269728*pt));
189  else if(fabs(eta)<0.5) res = 0.0226496+0.162875*exp(-(0.0250225*pt));
190  else if(fabs(eta)<0.7) res = 0.0232004+0.167887*exp(-(0.0259427*pt));
191  else if(fabs(eta)<0.9) res = 0.0220083+0.183889*exp(-(0.0273017*pt));
192  else if(fabs(eta)<1.15) res = 0.0216601+0.18482*exp(-(0.027148*pt));
193  else if(fabs(eta)<1.4) res = 0.0225191+0.19558*exp(-(0.0271747*pt));
194  else if(fabs(eta)<1.7) res = 0.0208365+0.19241*exp(-(0.0262627*pt));
195  else if(fabs(eta)<2.1) res = 0.0233749+0.180684*exp(-(0.032052*pt));
196  else if(fabs(eta)<2.5) res = 0.0231891+0.165484*exp(-(0.0389916*pt));
197  }else{
198  if(fabs(eta)<0.17) res = 0.025945+0.16529*exp(-(0.027385*pt));
199  else if(fabs(eta)<0.35) res = 0.0257294+0.155366*exp(-(0.0268639*pt));
200  else if(fabs(eta)<0.5) res = 0.0231924+0.166508*exp(-(0.0264153*pt));
201  else if(fabs(eta)<0.7) res = 0.0252736+0.176977*exp(-(0.0281594*pt));
202  else if(fabs(eta)<0.9) res = 0.0269299+0.174505*exp(-(0.0295757*pt));
203  else if(fabs(eta)<1.15) res = 0.0250191+0.178321*exp(-(0.0284006*pt));
204  else if(fabs(eta)<1.4) res = 0.0242345+0.174682*exp(-(0.027202*pt));
205  else if(fabs(eta)<1.7) res = 0.0242234+0.199735*exp(-(0.0291262*pt));
206  else if(fabs(eta)<2.1) res = 0.0300017+0.175617*exp(-(0.0357093*pt));
207  else if(fabs(eta)<2.5) res = 0.0282978+0.141994*exp(-(0.0383663*pt));
208  }
209  return res;
210 }
211 
212 inline double res::HelperJet::et(double pt, double eta, Flavor flav)
213 {
214  double res = 0.29*sqrt(pt);
215  if(flav==kB){
216  if(fabs(eta)<0.17) res = 10.1445+0.0445196*pt;
217  else if(fabs(eta)<0.35) res = 8.944+0.0576391*pt;
218  else if(fabs(eta)<0.5) res = 8.8462+0.0630703*pt;
219  else if(fabs(eta)<0.7) res = 8.90812+0.057605*pt;
220  else if(fabs(eta)<0.9) res = 8.42307+0.0623024*pt;
221  else if(fabs(eta)<1.15) res = 8.64735+0.0619689*pt;
222  else if(fabs(eta)<1.4) res = 8.90574+0.0595526*pt;
223  else if(fabs(eta)<1.7) res = 8.74415+0.0612412*pt;
224  else if(fabs(eta)<2.1) res = 7.80525+0.0511843*pt;
225  else if(fabs(eta)<2.5) res = 6.07877+0.0569574*pt;
226  }else{
227  if(fabs(eta)<0.17) res = 8.04345+0.0705188*pt;
228  else if(fabs(eta)<0.35) res = 7.77565+0.0722247*pt;
229  else if(fabs(eta)<0.5) res = 8.53558+0.0647725*pt;
230  else if(fabs(eta)<0.7) res = 7.99374+0.0681428*pt;
231  else if(fabs(eta)<0.9) res = 7.51695+0.0755128*pt;
232  else if(fabs(eta)<1.15) res = 7.70804+0.0748981*pt;
233  else if(fabs(eta)<1.4) res = 7.2996+0.0852397*pt;
234  else if(fabs(eta)<1.7) res = 7.95923+0.0715722*pt;
235  else if(fabs(eta)<2.1) res = 7.68012+0.0543317*pt;
236  else if(fabs(eta)<2.5) res = 6.59141+0.0452384*pt;
237  }
238  return res;
239 }
240 
241 inline double res::HelperJet::eta(double pt, double eta, Flavor flav)
242 {
243  double res=-1.53e-4*pt+0.05;
244  if(flav==kB){
245  if(fabs(eta)<0.17) res = 0.0354411+0.130211*exp(-(0.0423321*pt));
246  else if(fabs(eta)<0.35) res = 0.0341742+0.115784*exp(-(0.0375787*pt));
247  else if(fabs(eta)<0.5) res = 0.0336537+0.143109*exp(-(0.0416702*pt));
248  else if(fabs(eta)<0.7) res = 0.031894+0.140572*exp(-(0.0405564*pt));
249  else if(fabs(eta)<0.9) res = 0.0301507+0.126005*exp(-(0.0384231*pt));
250  else if(fabs(eta)<1.15) res = 0.0268018+0.140956*exp(-(0.03872*pt));
251  else if(fabs(eta)<1.4) res = 0.0242941+0.131887*exp(-(0.0337869*pt));
252  else if(fabs(eta)<1.7) res = 0.0223985+0.142949*exp(-(0.0351412*pt));
253  else if(fabs(eta)<2.1) res = 0.0229694+0.145993*exp(-(0.0391759*pt));
254  else if(fabs(eta)<2.5) res = 0.0243171+0.145419*exp(-(0.0430824*pt));
255  }else{
256  if(fabs(eta)<0.17) res = 0.038575+0.129026*exp(-(0.0427157*pt));
257  else if(fabs(eta)<0.35) res = 0.0400221+0.127222*exp(-(0.0458212*pt));
258  else if(fabs(eta)<0.5) res = 0.0358472+0.119962*exp(-(0.0395586*pt));
259  else if(fabs(eta)<0.7) res = 0.0359612+0.136094*exp(-(0.04212*pt));
260  else if(fabs(eta)<0.9) res = 0.0327581+0.127339*exp(-(0.0393077*pt));
261  else if(fabs(eta)<1.15) res = 0.0305095+0.116238*exp(-(0.0347537*pt));
262  else if(fabs(eta)<1.4) res = 0.027529+0.117893*exp(-(0.032858*pt));
263  else if(fabs(eta)<1.7) res = 0.0253065+0.112751*exp(-(0.0307028*pt));
264  else if(fabs(eta)<2.1) res = 0.0266957+0.131998*exp(-(0.0380667*pt));
265  else if(fabs(eta)<2.5) res = 0.0298413+0.13651*exp(-(0.0458475*pt));
266  }
267  return res;
268 }
269 
270 #endif
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:25
create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1087
res::HelperJet::phi
double phi(double pt, double eta, Flavor flav)
Definition: Jet.h:183
res::HelperJet::theta
double theta(double pt, double eta, Flavor flav)
Definition: Jet.h:154
res::HelperJet::et
double et(double pt, double eta, Flavor flav)
Definition: Jet.h:212
res::HelperJet::~HelperJet
~HelperJet()
Definition: Jet.h:11
res::HelperJet::b
double b(double pt, double eta, Flavor flav)
Definition: Jet.h:67
res::HelperJet::d
double d(double pt, double eta, Flavor flav)
Definition: Jet.h:125
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:241
res::HelperJet::HelperJet
HelperJet()
Definition: Jet.h:10
res::HelperJet::c
double c(double pt, double eta, Flavor flav)
Definition: Jet.h:96
res::HelperJet::a
double a(double pt, double eta, Flavor flav)
Definition: Jet.h:38
Generated for CMSSW Reference Manual by doxygen 1.8.5