-
Notifications
You must be signed in to change notification settings - Fork 0
/
jec_2nd_pbp.h
247 lines (246 loc) · 29.7 KB
/
jec_2nd_pbp.h
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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
if (pseudorapidity >= -5.191 && pseudorapidity < -4.889) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((4.89709796e-01)+((5.51503399e-01)/((pow(log10(perp),2))+(2.89951353e+00))))+((4.01385463e-01)*exp(-((-2.24192156e-02)*pow(log10(perp)-(1.70978941e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((3.27038260e-01)+(3.21819561e-01)*(pow(log10(perp),-4.83096413e-01)));
}
if (pseudorapidity >= -4.889 && pseudorapidity < -4.716) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.20986174e-01)+((4.03848998e-01)/((pow(log10(perp),2))+(3.50422547e+00))))+((4.08218527e-01)*exp(-((-1.38254333e-02)*pow(log10(perp)-(1.49721617e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((3.11576203e-01)+(3.07153015e-01)*(pow(log10(perp),-2.63778004e-01)));
}
if (pseudorapidity >= -4.716 && pseudorapidity < -4.538) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.58601138e-01)+((1.82918712e-01)/((pow(log10(perp),2))+(4.36072077e+00))))+((4.17635486e-01)*exp(-((-3.93296677e-03)*pow(log10(perp)-(1.83756432e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.99447346e-01)+(2.93295764e-01)*(pow(log10(perp),-6.71913793e-02)));
}
if (pseudorapidity >= -4.538 && pseudorapidity < -4.363) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.58936406e-01)+((-8.50403286e-01)/((pow(log10(perp),2))+(7.69975566e+00))))+((4.39013262e-01)*exp(-((2.15401149e-02)*pow(log10(perp)-(5.19128479e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.84170599e-01)+(2.69638437e-01)*(pow(log10(perp),1.97911148e-01)));
}
if (pseudorapidity >= -4.363 && pseudorapidity < -4.191) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.03729069e-01)+((-3.70037272e+00)/((pow(log10(perp),2))+(1.51025495e+01))))+((4.12165497e-01)*exp(-((7.51979304e-02)*pow(log10(perp)-(1.29652486e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((1.97611962e-01)+(3.61410821e-01)*(pow(log10(perp),2.61089297e-01)));
}
if (pseudorapidity >= -4.191 && pseudorapidity < -4.013) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.10746913e+00)+((-7.66630584e+00)/((pow(log10(perp),2))+(1.88116668e+01))))+((2.68834583e-01)*exp(-((1.39632879e-01)*pow(log10(perp)-(1.02001153e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-3.27373122e-01)+(1.00284987e+00)*(pow(log10(perp),1.83238581e-01)));
}
if (pseudorapidity >= -4.013 && pseudorapidity < -3.839) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.15811655e+00)+((-7.92750886e+00)/((pow(log10(perp),2))+(2.13239745e+01))))+((1.85852972e-01)*exp(-((2.90955234e-01)*pow(log10(perp)-(1.36579487e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-6.14202657e-01)+(1.34747383e+00)*(pow(log10(perp),1.64252405e-01)));
}
if (pseudorapidity >= -3.839 && pseudorapidity < -3.664) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.92128704e-01)+((-3.29787901e+00)/((pow(log10(perp),2))+(1.41167149e+01))))+((3.29259712e-01)*exp(-((1.03849065e-01)*pow(log10(perp)-(1.09668115e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.13236237e-01)+(6.55524578e-01)*(pow(log10(perp),4.40377454e-01)));
}
if (pseudorapidity >= -3.664 && pseudorapidity < -3.489) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.25922016e-01)+((-4.33866708e+00)/((pow(log10(perp),2))+(2.55150895e+01))))+((3.42236602e-01)*exp(-((1.19973522e-01)*pow(log10(perp)-(1.56971902e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-4.65807848e-01)+(1.01563075e+00)*(pow(log10(perp),2.58761092e-01)));
}
if (pseudorapidity >= -3.489 && pseudorapidity < -3.314) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.11497916e-01)+((-1.74999562e+00)/((pow(log10(perp),2))+(7.45045241e+00))))+((3.92403703e-01)*exp(-((9.86075721e-02)*pow(log10(perp)-(9.91573173e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((8.89325633e-02)+(3.75398411e-01)*(pow(log10(perp),8.30562124e-01)));
}
if (pseudorapidity >= -3.314 && pseudorapidity < -3.139) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.03788662e+00)+((-6.17961818e+00)/((pow(log10(perp),2))+(1.72287940e+01))))+((2.88108253e-01)*exp(-((2.04915753e-01)*pow(log10(perp)-(1.37486064e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-5.46087822e-01)+(1.15871681e+00)*(pow(log10(perp),2.15052592e-01)));
}
if (pseudorapidity >= -3.139 && pseudorapidity < -2.964) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.92315895e-01)+((-4.28140540e+00)/((pow(log10(perp),2))+(2.07559286e+01))))+((2.43036420e-01)*exp(-((6.88985649e-02)*pow(log10(perp)-(1.41470957e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-3.35264796e-01)+(9.87217579e-01)*(pow(log10(perp),3.19488364e-01)));
}
if (pseudorapidity >= -2.964 && pseudorapidity < -2.853) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.00730531e+00)+((-4.10439351e+00)/((pow(log10(perp),2))+(1.93589751e+01))))+((1.98767118e-01)*exp(-((1.09902298e-01)*pow(log10(perp)-(9.51461609e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-4.02945800e-01)+(1.08779629e+00)*(pow(log10(perp),2.67163622e-01)));
}
if (pseudorapidity >= -2.853 && pseudorapidity < -2.65) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.16446811e+00)+((-4.80149644e+00)/((pow(log10(perp),2))+(2.04955740e+01))))+((6.59656620e-02)*exp(-((3.38602388e-01)*pow(log10(perp)-(5.60975485e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.77665560e-01)+(1.61220758e+00)*(pow(log10(perp),1.85746630e-01)));
}
if (pseudorapidity >= -2.65 && pseudorapidity < -2.5) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.70022127e-01)+((-3.35557019e+00)/((pow(log10(perp),2))+(1.73516666e+01))))+((2.06766255e-01)*exp(-((6.32975869e-02)*pow(log10(perp)-(9.16328517e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-3.29380414e-01)+(1.00159594e+00)*(pow(log10(perp),3.08412604e-01)));
}
if (pseudorapidity >= -2.5 && pseudorapidity < -2.322) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.91209637e-01)+((-4.98043688e-01)/((pow(log10(perp),2))+(7.09507154e+00))))+((3.71364177e-01)*exp(-((1.75891608e-02)*pow(log10(perp)-(7.39500534e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.26960044e-01)+(3.38376131e-01)*(pow(log10(perp),9.98562093e-01)));
}
if (pseudorapidity >= -2.322 && pseudorapidity < -2.172) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.85289241e-01)+((-2.58297265e+00)/((pow(log10(perp),2))+(1.38130833e+01))))+((2.74481453e-01)*exp(-((5.80935923e-02)*pow(log10(perp)-(1.26726953e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.47570500e-01)+(7.51518872e-01)*(pow(log10(perp),3.97050835e-01)));
}
if (pseudorapidity >= -2.172 && pseudorapidity < -2.043) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.21498983e+00)+((-4.82573411e+00)/((pow(log10(perp),2))+(1.58445360e+01))))+((6.32047718e-02)*exp(-((6.60601148e-01)*pow(log10(perp)-(9.66277429e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-8.20436979e-01)+(1.66413280e+00)*(pow(log10(perp),1.65352141e-01)));
}
if (pseudorapidity >= -2.043 && pseudorapidity < -1.93) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.17181487e+00)+((-5.04140536e+00)/((pow(log10(perp),2))+(2.04145329e+01))))+((6.56638118e-02)*exp(-((3.58382167e-01)*pow(log10(perp)-(8.24860697e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.86020398e-01)+(1.62087773e+00)*(pow(log10(perp),1.72390223e-01)));
}
if (pseudorapidity >= -1.93 && pseudorapidity < -1.83) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.42385703e-01)+((-2.40892873e+00)/((pow(log10(perp),2))+(2.59719739e+01))))+((1.70894498e-01)*exp(-((1.52235953e-02)*pow(log10(perp)-(-1.30398671e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-2.76854464e-01)+(9.90286199e-01)*(pow(log10(perp),3.25208678e-01)));
}
if (pseudorapidity >= -1.83 && pseudorapidity < -1.74) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.77148533e-01)+((-1.62111057e-01)/((pow(log10(perp),2))+(6.48871221e+00))))+((3.57047357e-01)*exp(-((1.19267084e-02)*pow(log10(perp)-(3.10504000e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((8.73666571e-02)+(3.78597004e-01)*(pow(log10(perp),1.01165186e+00)));
}
if (pseudorapidity >= -1.74 && pseudorapidity < -1.653) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.61767353e-01)+((5.14213181e-01)/((pow(log10(perp),2))+(3.74852751e+00))))+((3.59824863e-01)*exp(-((-1.50264603e-02)*pow(log10(perp)-(2.29656350e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.69992193e-01)+(2.35654414e-01)*(pow(log10(perp),1.70965022e+00)));
}
if (pseudorapidity >= -1.653 && pseudorapidity < -1.566) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.28242723e-01)+((7.88960667e-01)/((pow(log10(perp),2))+(4.46486882e+00))))+((3.56883975e-01)*exp(-((-1.84989614e-02)*pow(log10(perp)-(5.15084393e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((1.87044512e-01)+(2.18709283e-01)*(pow(log10(perp),1.75696501e+00)));
}
if (pseudorapidity >= -1.566 && pseudorapidity < -1.479) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.03761968e-01)+((-4.09161289e-01)/((pow(log10(perp),2))+(7.36497319e+00))))+((3.55970710e-01)*exp(-((2.21685210e-02)*pow(log10(perp)-(4.87618443e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((7.45064383e-02)+(4.00022750e-01)*(pow(log10(perp),9.30378542e-01)));
}
if (pseudorapidity >= -1.479 && pseudorapidity < -1.392) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.12561345e+00)+((-4.30842108e+00)/((pow(log10(perp),2))+(2.06402979e+01))))+((8.35624869e-02)*exp(-((2.05537337e-01)*pow(log10(perp)-(6.32818780e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.16035154e-01)+(1.51221460e+00)*(pow(log10(perp),2.14021763e-01)));
}
if (pseudorapidity >= -1.392 && pseudorapidity < -1.305) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.12342053e+00)+((-3.05663561e+00)/((pow(log10(perp),2))+(1.69178449e+01))))+((7.45935190e-02)*exp(-((2.84554091e-01)*pow(log10(perp)-(1.28883336e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.84269969e-01)+(1.61207594e+00)*(pow(log10(perp),1.93033876e-01)));
}
if (pseudorapidity >= -1.305 && pseudorapidity < -1.218) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.13134536e+00)+((-4.15321996e+00)/((pow(log10(perp),2))+(2.04740481e+01))))+((7.40522628e-02)*exp(-((2.50805459e-01)*pow(log10(perp)-(6.09734007e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.27228769e-01)+(1.53597720e+00)*(pow(log10(perp),2.06188501e-01)));
}
if (pseudorapidity >= -1.218 && pseudorapidity < -1.131) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.19874928e-01)+((-4.96480931e-01)/((pow(log10(perp),2))+(1.44819741e+01))))+((3.23144625e-01)*exp(-((4.48417078e-03)*pow(log10(perp)-(-5.19017521e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.43512772e-02)+(5.00979243e-01)*(pow(log10(perp),6.97873304e-01)));
}
if (pseudorapidity >= -1.131 && pseudorapidity < -1.044) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.29648299e-01)+((5.68318928e-01)/((pow(log10(perp),2))+(3.38234246e+00))))+((3.70469104e-01)*exp(-((-2.30212407e-02)*pow(log10(perp)-(1.88204940e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.21452123e-01)+(2.06109948e-01)*(pow(log10(perp),1.53919194e+00)));
}
if (pseudorapidity >= -1.044 && pseudorapidity < -0.957) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.74520673e-01)+((1.44694988e-01)/((pow(log10(perp),2))+(1.23983587e+01))))+((3.22904519e-01)*exp(-((1.00503286e-03)*pow(log10(perp)-(-1.09380474e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((3.21489198e-02)+(4.88191327e-01)*(pow(log10(perp),7.26144764e-01)));
}
if (pseudorapidity >= -0.957 && pseudorapidity < -0.879) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.02946385e-01)+((-1.77380037e+00)/((pow(log10(perp),2))+(2.76247127e+01))))+((1.73650505e-01)*exp(-((6.62543763e-03)*pow(log10(perp)-(-1.85818988e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-2.55245234e-01)+(9.55493386e-01)*(pow(log10(perp),3.50024100e-01)));
}
if (pseudorapidity >= -0.879 && pseudorapidity < -0.783) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.16520542e-01)+((-1.49897920e+00)/((pow(log10(perp),2))+(1.08749817e+01))))+((3.08338796e-01)*exp(-((3.86577489e-02)*pow(log10(perp)-(7.78467080e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-6.84623534e-02)+(6.13454781e-01)*(pow(log10(perp),5.77290701e-01)));
}
if (pseudorapidity >= -0.783 && pseudorapidity < -0.696) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.69986947e-01)+((3.77228570e-01)/((pow(log10(perp),2))+(2.53704383e+00))))+((3.54912016e-01)*exp(-((-2.32023724e-02)*pow(log10(perp)-(3.60666694e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.59825728e-01)+(2.79920846e-01)*(pow(log10(perp),1.32186660e+00)));
}
if (pseudorapidity >= -0.696 && pseudorapidity < -0.609) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.02984076e-01)+((6.84128083e-01)/((pow(log10(perp),2))+(2.90731153e+00))))+((3.59688267e-01)*exp(-((-3.54610883e-02)*pow(log10(perp)-(2.74778967e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.22312910e-01)+(1.85390190e-01)*(pow(log10(perp),1.89378045e+00)));
}
if (pseudorapidity >= -0.609 && pseudorapidity < -0.522) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.07039332e-01)+((-2.07793015e-01)/((pow(log10(perp),2))+(1.43650928e+01))))+((3.16140043e-01)*exp(-((2.24561623e-03)*pow(log10(perp)-(-6.94667629e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.61930596e-02)+(5.10500521e-01)*(pow(log10(perp),7.07745220e-01)));
}
if (pseudorapidity >= -0.522 && pseudorapidity < -0.435) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.13920913e+00)+((-4.10287852e+00)/((pow(log10(perp),2))+(2.01084075e+01))))+((6.32954040e-02)*exp(-((2.67012421e-01)*pow(log10(perp)-(6.06159378e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.49596863e-01)+(1.57059307e+00)*(pow(log10(perp),2.07601725e-01)));
}
if (pseudorapidity >= -0.435 && pseudorapidity < -0.348) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.42321368e-01)+((-1.17594682e+00)/((pow(log10(perp),2))+(2.65052929e+01))))+((2.13433603e-01)*exp(-((2.85400951e-03)*pow(log10(perp)-(-2.26534199e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.95306484e-01)+(8.51669856e-01)*(pow(log10(perp),3.90503666e-01)));
}
if (pseudorapidity >= -0.348 && pseudorapidity < -0.261) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.07647740e-01)+((6.62864064e-01)/((pow(log10(perp),2))+(2.95263456e+00))))+((3.60019410e-01)*exp(-((-3.34753603e-02)*pow(log10(perp)-(2.15670179e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.18090533e-01)+(2.09483169e-01)*(pow(log10(perp),1.56983401e+00)));
}
if (pseudorapidity >= -0.261 && pseudorapidity < -0.174) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.88822381e-01)+((2.69128931e-01)/((pow(log10(perp),2))+(1.80987056e+00))))+((3.45837101e-01)*exp(-((-2.28803199e-02)*pow(log10(perp)-(3.33338007e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.22936143e-01)+(2.96478350e-01)*(pow(log10(perp),1.44270469e+00)));
}
if (pseudorapidity >= -0.174 && pseudorapidity < -0.087) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.05050723e-01)+((3.14397066e-01)/((pow(log10(perp),2))+(2.68157098e+00))))+((3.36673206e-01)*exp(-((-2.29723150e-02)*pow(log10(perp)-(5.44133518e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((9.92793159e-02)+(2.84474607e-01)*(pow(log10(perp),1.77226806e+00)));
}
if (pseudorapidity >= -0.087 && pseudorapidity < 0.0) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.23321532e-01)+((6.63324730e-01)/((pow(log10(perp),2))+(2.95264036e+00))))+((3.48342625e-01)*exp(-((-3.02052019e-02)*pow(log10(perp)-(1.51951856e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.75013730e-01)+(1.26399506e-01)*(pow(log10(perp),3.30187951e+00)));
}
if (pseudorapidity >= 0.0 && pseudorapidity < 0.087) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.90192455e-01)+((7.28353567e-01)/((pow(log10(perp),2))+(5.03418488e+00))))+((3.06242501e-01)*exp(-((-3.11199751e-02)*pow(log10(perp)-(3.49782253e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.31931026e-01)+(3.30371943e-01)*(pow(log10(perp),1.34945959e+00)));
}
if (pseudorapidity >= 0.087 && pseudorapidity < 0.174) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.84278007e-01)+((-4.87804617e-01)/((pow(log10(perp),2))+(2.20043177e+01))))+((2.33921042e-01)*exp(-((-1.59039498e-03)*pow(log10(perp)-(-1.86925025e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.49677976e-01)+(7.69518793e-01)*(pow(log10(perp),4.32638862e-01)));
}
if (pseudorapidity >= 0.174 && pseudorapidity < 0.261) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.15707432e+00)+((-3.81866056e+00)/((pow(log10(perp),2))+(1.78504832e+01))))+((5.64615422e-02)*exp(-((4.07427939e-01)*pow(log10(perp)-(4.86490958e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.63685323e-01)+(1.61599743e+00)*(pow(log10(perp),1.78972517e-01)));
}
if (pseudorapidity >= 0.261 && pseudorapidity < 0.348) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.17442550e+00)+((-4.54208243e+00)/((pow(log10(perp),2))+(1.88716716e+01))))+((5.48298339e-02)*exp(-((4.19489090e-01)*pow(log10(perp)-(7.37668003e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.74578823e-01)+(1.63566661e+00)*(pow(log10(perp),1.61165760e-01)));
}
if (pseudorapidity >= 0.348 && pseudorapidity < 0.435) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.21034787e+00)+((-5.16446407e+00)/((pow(log10(perp),2))+(1.76299796e+01))))+((6.04915724e-02)*exp(-((6.39416542e-01)*pow(log10(perp)-(1.02197132e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.85606415e-01)+(1.63042382e+00)*(pow(log10(perp),1.65874329e-01)));
}
if (pseudorapidity >= 0.435 && pseudorapidity < 0.522) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.32878999e-01)+((-1.64670698e+00)/((pow(log10(perp),2))+(3.66754290e+01))))+((1.17175853e-01)*exp(-((1.18347453e-03)*pow(log10(perp)-(-3.69910470e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-3.68068421e-01)+(1.13201518e+00)*(pow(log10(perp),2.63136108e-01)));
}
if (pseudorapidity >= 0.522 && pseudorapidity < 0.609) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.11010587e-01)+((-1.39325940e-01)/((pow(log10(perp),2))+(1.65311336e+01))))+((2.98674964e-01)*exp(-((-6.00504070e-04)*pow(log10(perp)-(-1.44365439e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-2.44338946e-02)+(5.84430173e-01)*(pow(log10(perp),5.35962405e-01)));
}
if (pseudorapidity >= 0.609 && pseudorapidity < 0.696) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((8.99062129e-01)+((-1.64567810e+00)/((pow(log10(perp),2))+(3.05114001e+01))))+((1.58787821e-01)*exp(-((1.78997391e-03)*pow(log10(perp)-(-2.56874701e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-2.88296731e-01)+(1.00098629e+00)*(pow(log10(perp),3.26980879e-01)));
}
if (pseudorapidity >= 0.696 && pseudorapidity < 0.783) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.81456880e-01)+((-3.80367026e-01)/((pow(log10(perp),2))+(2.50045885e+01))))+((2.30720530e-01)*exp(-((-1.26704000e-03)*pow(log10(perp)-(-3.04070730e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.65009809e-01)+(7.95708565e-01)*(pow(log10(perp),3.79166917e-01)));
}
if (pseudorapidity >= 0.783 && pseudorapidity < 0.879) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((9.57687966e-01)+((-2.29730563e+00)/((pow(log10(perp),2))+(2.86357742e+01))))+((1.25928951e-01)*exp(-((5.52658179e-03)*pow(log10(perp)-(-1.32586744e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-3.74292263e-01)+(1.13661374e+00)*(pow(log10(perp),2.56348862e-01)));
}
if (pseudorapidity >= 0.879 && pseudorapidity < 0.957) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.16031803e+00)+((-3.77461711e+00)/((pow(log10(perp),2))+(1.62310826e+01))))+((7.51089616e-02)*exp(-((4.30005083e-01)*pow(log10(perp)-(5.54963886e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.88539708e-01)+(1.62616705e+00)*(pow(log10(perp),1.84399451e-01)));
}
if (pseudorapidity >= 0.957 && pseudorapidity < 1.044) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.10477869e+00)+((-2.36988296e+00)/((pow(log10(perp),2))+(1.42758041e+01))))+((9.04189235e-02)*exp(-((2.95678497e-01)*pow(log10(perp)-(4.59132461e-03),2))))+(0.5*(copysign(1,30-perp)+1))*((-8.15139489e-01)+(1.66086443e+00)*(pow(log10(perp),1.76542229e-01)));
}
if (pseudorapidity >= 1.044 && pseudorapidity < 1.131) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.09619599e+00)+((-2.07532670e+00)/((pow(log10(perp),2))+(1.45058157e+01))))+((9.51159359e-02)*exp(-((2.61341276e-01)*pow(log10(perp)-(-2.88024481e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.97079818e-01)+(1.64849286e+00)*(pow(log10(perp),1.82738233e-01)));
}
if (pseudorapidity >= 1.131 && pseudorapidity < 1.218) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.16701203e+00)+((-3.89888696e+00)/((pow(log10(perp),2))+(1.57775926e+01))))+((9.90196117e-02)*exp(-((4.86211838e-01)*pow(log10(perp)-(5.58344855e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.85658254e-01)+(1.60692095e+00)*(pow(log10(perp),1.87980542e-01)));
}
if (pseudorapidity >= 1.218 && pseudorapidity < 1.305) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.18177363e+00)+((-4.91110536e+00)/((pow(log10(perp),2))+(1.73253938e+01))))+((9.19108672e-02)*exp(-((4.46150064e-01)*pow(log10(perp)-(9.23598273e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.57162168e-01)+(1.56621055e+00)*(pow(log10(perp),1.70461870e-01)));
}
if (pseudorapidity >= 1.305 && pseudorapidity < 1.392) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.07880302e+00)+((-1.73011463e+00)/((pow(log10(perp),2))+(1.38605318e+01))))+((1.01489958e-01)*exp(-((2.79114352e-01)*pow(log10(perp)-(-3.43720409e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.80233159e-01)+(1.61135997e+00)*(pow(log10(perp),2.08077165e-01)));
}
if (pseudorapidity >= 1.392 && pseudorapidity < 1.479) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.79319062e-01)+((-3.98684224e-01)/((pow(log10(perp),2))+(2.70801719e+01))))+((2.39958801e-01)*exp(-((-1.42287117e-04)*pow(log10(perp)-(-3.23322897e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.87380002e-01)+(8.06755237e-01)*(pow(log10(perp),4.04325508e-01)));
}
if (pseudorapidity >= 1.479 && pseudorapidity < 1.566) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.33468459e-01)+((1.13516304e+00)/((pow(log10(perp),2))+(5.10589772e+00))))+((3.06756086e-01)*exp(-((-3.45641166e-02)*pow(log10(perp)-(1.94644358e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.37552498e-01)+(3.58207594e-01)*(pow(log10(perp),1.04680284e+00)));
}
if (pseudorapidity >= 1.566 && pseudorapidity < 1.653) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.04987219e-01)+((6.82360347e-01)/((pow(log10(perp),2))+(2.96863277e+00))))+((3.63541237e-01)*exp(-((-3.20831341e-02)*pow(log10(perp)-(2.75282974e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((2.27811664e-01)+(1.62172959e-01)*(pow(log10(perp),2.17018513e+00)));
}
if (pseudorapidity >= 1.653 && pseudorapidity < 1.74) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((4.44528473e-01)+((3.80496785e+00)/((pow(log10(perp),2))+(9.73574903e+00))))+((1.91246035e-01)*exp(-((-4.31619891e-02)*pow(log10(perp)-(-9.93909355e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((4.12394851e-02)+(5.65222532e-01)*(pow(log10(perp),5.79635250e-01)));
}
if (pseudorapidity >= 1.74 && pseudorapidity < 1.83) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.12289791e+00)+((-3.19562591e+00)/((pow(log10(perp),2))+(1.79310691e+01))))+((8.42046461e-02)*exp(-((2.80356024e-01)*pow(log10(perp)-(-1.56147388e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.60285376e-01)+(1.58443211e+00)*(pow(log10(perp),2.14429606e-01)));
}
if (pseudorapidity >= 1.83 && pseudorapidity < 1.93) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.12259593e+00)+((-3.33208296e+00)/((pow(log10(perp),2))+(1.81867116e+01))))+((8.72168476e-02)*exp(-((2.74274859e-01)*pow(log10(perp)-(4.65715094e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.49698317e-01)+(1.56039468e+00)*(pow(log10(perp),2.20494802e-01)));
}
if (pseudorapidity >= 1.93 && pseudorapidity < 2.043) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.55240762e-01)+((1.08846331e+00)/((pow(log10(perp),2))+(5.03911322e+00))))+((2.83678260e-01)*exp(-((-3.40218037e-02)*pow(log10(perp)-(-7.32535028e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((1.71540082e-01)+(3.25094307e-01)*(pow(log10(perp),1.21794334e+00)));
}
if (pseudorapidity >= 2.043 && pseudorapidity < 2.172) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.21307155e-01)+((6.84629808e-01)/((pow(log10(perp),2))+(2.88295665e+00))))+((3.40628118e-01)*exp(-((-3.70866275e-02)*pow(log10(perp)-(2.04496925e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.59378586e-01)+(1.51682927e-01)*(pow(log10(perp),2.95625185e+00)));
}
if (pseudorapidity >= 2.172 && pseudorapidity < 2.322) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((6.16283216e-01)+((1.91644282e-01)/((pow(log10(perp),2))+(1.74802003e+00))))+((3.40426384e-01)*exp(-((-1.42872362e-02)*pow(log10(perp)-(4.69626417e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((9.49630668e-02)+(3.02555489e-01)*(pow(log10(perp),1.61600779e+00)));
}
if (pseudorapidity >= 2.322 && pseudorapidity < 2.5) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.82437094e-01)+((4.23448081e-01)/((pow(log10(perp),2))+(4.18017534e+00))))+((3.63918384e-01)*exp(-((-6.67981621e-03)*pow(log10(perp)-(2.72583223e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((1.63610200e-01)+(2.47424610e-01)*(pow(log10(perp),1.61526975e+00)));
}
if (pseudorapidity >= 2.5 && pseudorapidity < 2.65) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.44208929e-01)+((-3.64322604e-01)/((pow(log10(perp),2))+(2.16740678e+01))))+((2.75281972e-01)*exp(-((1.96228548e-04)*pow(log10(perp)-(-2.26361582e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-8.41026345e-02)+(6.61538473e-01)*(pow(log10(perp),5.22261957e-01)));
}
if (pseudorapidity >= 2.65 && pseudorapidity < 2.853) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((4.67697198e-01)+((3.30183607e+00)/((pow(log10(perp),2))+(8.56003250e+00))))+((1.81612283e-01)*exp(-((-5.12744870e-02)*pow(log10(perp)-(-7.79443880e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-6.21516830e-02)+(7.09160912e-01)*(pow(log10(perp),3.77410250e-01)));
}
if (pseudorapidity >= 2.853 && pseudorapidity < 2.964) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.02176385e+00)+((-3.40133952e+00)/((pow(log10(perp),2))+(1.92662894e+01))))+((2.37977528e-01)*exp(-((7.59290990e-02)*pow(log10(perp)-(-5.87773362e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-4.35552430e-01)+(1.13360936e+00)*(pow(log10(perp),3.16563617e-01)));
}
if (pseudorapidity >= 2.964 && pseudorapidity < 3.139) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((3.43057094e-01)+((1.55507723e+00)/((pow(log10(perp),2))+(1.79395061e+01))))+((6.37527336e-01)*exp(-((7.82905680e-02)*pow(log10(perp)-(1.58055804e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.90375992e-01)+(4.57568532e-01)*(pow(log10(perp),6.03349382e-01)));
}
if (pseudorapidity >= 3.139 && pseudorapidity < 3.314) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.09433276e+00)+((-2.74804299e+00)/((pow(log10(perp),2))+(1.62393882e+01))))+((7.68711171e-02)*exp(-((1.95683221e-01)*pow(log10(perp)-(3.63157513e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.57988101e-01)+(1.56805250e+00)*(pow(log10(perp),1.94133792e-01)));
}
if (pseudorapidity >= 3.314 && pseudorapidity < 3.489) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.16426929e+00)+((-4.63055115e+00)/((pow(log10(perp),2))+(1.96684218e+01))))+((7.49404365e-02)*exp(-((3.56127659e-01)*pow(log10(perp)-(4.46490104e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.91001157e-01)+(1.63197001e+00)*(pow(log10(perp),1.80247167e-01)));
}
if (pseudorapidity >= 3.489 && pseudorapidity < 3.664) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.18603038e+00)+((-5.46060315e+00)/((pow(log10(perp),2))+(1.89991182e+01))))+((8.82724762e-02)*exp(-((3.45175520e-01)*pow(log10(perp)-(7.43460384e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.80600965e-01)+(1.60649942e+00)*(pow(log10(perp),1.71060793e-01)));
}
if (pseudorapidity >= 3.664 && pseudorapidity < 3.839) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.20468975e+00)+((-6.41816318e+00)/((pow(log10(perp),2))+(2.05285045e+01))))+((9.30837749e-02)*exp(-((3.37901288e-01)*pow(log10(perp)-(7.74741171e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.91436297e-01)+(1.62097467e+00)*(pow(log10(perp),1.62049943e-01)));
}
if (pseudorapidity >= 3.839 && pseudorapidity < 4.013) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.21287175e+00)+((-7.48373802e+00)/((pow(log10(perp),2))+(2.42925887e+01))))+((8.82739371e-02)*exp(-((2.78510585e-01)*pow(log10(perp)-(6.05859771e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.66883753e-01)+(1.61484000e+00)*(pow(log10(perp),1.54321806e-01)));
}
if (pseudorapidity >= 4.013 && pseudorapidity < 4.191) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.20727799e+00)+((-7.76820442e+00)/((pow(log10(perp),2))+(2.62590118e+01))))+((7.58526970e-02)*exp(-((2.99999633e-01)*pow(log10(perp)-(8.08901973e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-7.78892608e-01)+(1.63805965e+00)*(pow(log10(perp),1.30592495e-01)));
}
if (pseudorapidity >= 4.191 && pseudorapidity < 4.363) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((1.06091674e+00)+((-6.68215016e+00)/((pow(log10(perp),2))+(1.68725476e+01))))+((3.10168279e-01)*exp(-((8.56229723e-02)*pow(log10(perp)-(5.97088792e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((-1.73956417e-01)+(8.16741578e-01)*(pow(log10(perp),2.30962617e-01)));
}
if (pseudorapidity >= 4.363 && pseudorapidity < 4.538) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((7.99754030e-01)+((-3.36075778e+00)/((pow(log10(perp),2))+(1.48842586e+01))))+((4.08731757e-01)*exp(-((6.32560835e-02)*pow(log10(perp)-(1.04510022e+00),2))))+(0.5*(copysign(1,30-perp)+1))*((1.74883102e-01)+(3.86569145e-01)*(pow(log10(perp),2.87913819e-01)));
}
if (pseudorapidity >= 4.538 && pseudorapidity < 4.716) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.74112167e-01)+((1.59792412e-01)/((pow(log10(perp),2))+(4.63797346e+00))))+((4.23992964e-01)*exp(-((5.10664183e-03)*pow(log10(perp)-(7.81711782e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((2.95804409e-01)+(2.88270519e-01)*(pow(log10(perp),1.33122471e-01)));
}
if (pseudorapidity >= 4.716 && pseudorapidity < 4.889) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((5.18751009e-01)+((4.85746993e-01)/((pow(log10(perp),2))+(3.35306606e+00))))+((4.11300575e-01)*exp(-((-8.48878665e-03)*pow(log10(perp)-(4.88813676e-02),2))))+(0.5*(copysign(1,30-perp)+1))*((3.15929325e-01)+(3.12213113e-01)*(pow(log10(perp),-2.09277332e-01)));
}
if (pseudorapidity >= 4.889 && pseudorapidity < 5.191) {
perp_scale = (0.5*(copysign(1,perp-30)+1))*((4.83842155e-01)+((6.39372692e-01)/((pow(log10(perp),2))+(2.69172573e+00))))+((4.06283081e-01)*exp(-((-1.77982114e-02)*pow(log10(perp)-(1.49331245e-01),2))))+(0.5*(copysign(1,30-perp)+1))*((3.41594148e-01)+(3.36023875e-01)*(pow(log10(perp),-5.82776480e-01)));
}