@@ -184,21 +184,31 @@ class GaisserH3a(CosmicRayFlux):
184184
185185 pdgids = PDGID_5COMP
186186 _funcs = (
187- lambda E : 0.7860 * E ** - 2.66 * exp (- E / (4e6 * 1 ))
188- + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 ))
189- + 0.00017 * E ** - 2.4 * exp (- E / (2e9 * 1 )),
190- lambda E : 0.3550 * E ** - 2.58 * exp (- E / (4e6 * 2 ))
191- + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 2 ))
192- + 0.00017 * E ** - 2.4 * exp (- E / (2e9 * 2 )),
193- lambda E : 0.2200 * E ** - 2.63 * exp (- E / (4e6 * 7 ))
194- + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 7 ))
195- + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 7 )),
196- lambda E : 0.1430 * E ** - 2.67 * exp (- E / (4e6 * 13 ))
197- + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 13 ))
198- + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 13 )),
199- lambda E : 0.2120 * E ** - 2.63 * exp (- E / (4e6 * 26 ))
200- + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 26 ))
201- + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 26 )),
187+ lambda E : (
188+ 0.7860 * E ** - 2.66 * exp (- E / (4e6 * 1 ))
189+ + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 ))
190+ + 0.00017 * E ** - 2.4 * exp (- E / (2e9 * 1 ))
191+ ),
192+ lambda E : (
193+ 0.3550 * E ** - 2.58 * exp (- E / (4e6 * 2 ))
194+ + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 2 ))
195+ + 0.00017 * E ** - 2.4 * exp (- E / (2e9 * 2 ))
196+ ),
197+ lambda E : (
198+ 0.2200 * E ** - 2.63 * exp (- E / (4e6 * 7 ))
199+ + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 7 ))
200+ + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 7 ))
201+ ),
202+ lambda E : (
203+ 0.1430 * E ** - 2.67 * exp (- E / (4e6 * 13 ))
204+ + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 13 ))
205+ + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 13 ))
206+ ),
207+ lambda E : (
208+ 0.2120 * E ** - 2.63 * exp (- E / (4e6 * 26 ))
209+ + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 26 ))
210+ + 0.000114 * E ** - 2.4 * exp (- E / (2e9 * 26 ))
211+ ),
202212 )
203213
204214
@@ -212,9 +222,9 @@ class GaisserH4a(CosmicRayFlux):
212222
213223 pdgids = PDGID_5COMP
214224 _funcs = (
215- lambda E : 0.7860 * E ** - 2.66 * exp ( - E / ( 4e6 * 1 ))
216- + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 ))
217- + 0.0200 * E ** - 2.6 * exp ( - E / 6e10 ),
225+ lambda E : (
226+ 0.7860 * E ** - 2.66 * exp ( - E / ( 4e6 * 1 )) + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 )) + 0.0200 * E ** - 2.6 * exp ( - E / 6e10 )
227+ ),
218228 lambda E : 0.3550 * E ** - 2.58 * exp (- E / (4e6 * 2 )) + 0.002 * E ** - 2.4 * exp (- E / (3e7 * 2 )),
219229 lambda E : 0.2200 * E ** - 2.63 * exp (- E / (4e6 * 7 )) + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 7 )),
220230 lambda E : 0.1430 * E ** - 2.67 * exp (- E / (4e6 * 13 )) + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 13 )),
@@ -233,14 +243,16 @@ class GaisserH4a_IT(CosmicRayFlux):
233243 # pylint: disable=invalid-name
234244 pdgids = PDGID_4COMP
235245 _funcs = (
236- lambda E : 0.7860 * E ** - 2.66 * exp ( - E / ( 4e6 * 1 ))
237- + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 ))
238- + 0.0200 * E ** - 2.6 * exp ( - E / 6e10 ),
246+ lambda E : (
247+ 0.7860 * E ** - 2.66 * exp ( - E / ( 4e6 * 1 )) + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 1 )) + 0.0200 * E ** - 2.6 * exp ( - E / 6e10 )
248+ ),
239249 lambda E : 0.3550 * E ** - 2.58 * exp (- E / (4e6 * 2 )) + 0.0020 * E ** - 2.4 * exp (- E / (3e7 * 2 )),
240- lambda E : 0.2200 * E ** - 2.63 * exp (- E / (4e6 * 7 ))
241- + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 7 ))
242- + 0.1430 * E ** - 2.67 * exp (- E / (4e6 * 13 ))
243- + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 13 )),
250+ lambda E : (
251+ 0.2200 * E ** - 2.63 * exp (- E / (4e6 * 7 ))
252+ + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 7 ))
253+ + 0.1430 * E ** - 2.67 * exp (- E / (4e6 * 13 ))
254+ + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 13 ))
255+ ),
244256 lambda E : 0.2120 * E ** - 2.63 * exp (- E / (4e6 * 26 )) + 0.00134 * E ** - 2.4 * exp (- E / (3e7 * 26 )),
245257 )
246258
@@ -304,9 +316,11 @@ class GlobalFitGST(CosmicRayFlux):
304316 lambda E : 0.32 * E ** - 2.58 * exp (- E / 1.2e5 / 2 ) + 0.0065 * E ** - 2.3 * exp (- E / 4e6 / 2 ),
305317 lambda E : 0.01 * E ** - 2.40 * exp (- E / 1.2e5 / 7 ) + 0.0006 * E ** - 2.3 * exp (- E / 4e6 / 7 ),
306318 lambda E : 0.013 * E ** - 2.40 * exp (- E / 1.2e5 / 13 ) + 0.0007 * E ** - 2.3 * exp (- E / 4e6 / 13 ),
307- lambda E : 0.006 * E ** - 2.30 * exp (- E / 1.2e5 / 26 )
308- + 0.00023 * E ** - 2.2 * exp (- E / 4e6 / 26 )
309- + 0.0000025 * E ** - 2.2 * exp (- E / 1.3e9 / 26 ),
319+ lambda E : (
320+ 0.006 * E ** - 2.30 * exp (- E / 1.2e5 / 26 )
321+ + 0.00023 * E ** - 2.2 * exp (- E / 4e6 / 26 )
322+ + 0.0000025 * E ** - 2.2 * exp (- E / 1.3e9 / 26 )
323+ ),
310324 )
311325
312326
@@ -320,13 +334,17 @@ class GlobalFitGST_IT(CosmicRayFlux): # pylint: disable=invalid-name
320334 _funcs = (
321335 lambda E : 0.7 * E ** - 2.66 * exp (- E / 1.2e5 ) + 0.015 * E ** - 2.4 * exp (- E / 4e6 ) + 0.0014 * E ** - 2.4 * exp (- E / 1.3e9 ),
322336 lambda E : 0.32 * E ** - 2.58 * exp (- E / 1.2e5 / 2 ) + 0.0065 * E ** - 2.3 * exp (- E / 4e6 / 2 ),
323- lambda E : 0.01 * E ** - 2.40 * exp (- E / 1.2e5 / 7 )
324- + 0.0006 * E ** - 2.3 * exp (- E / 4e6 / 7 )
325- + 0.013 * E ** - 2.40 * exp (- E / 1.2e5 / 13 )
326- + 0.0007 * E ** - 2.3 * exp (- E / 4e6 / 13 ),
327- lambda E : 0.006 * E ** - 2.30 * exp (- E / 1.2e5 / 26 )
328- + 0.00023 * E ** - 2.2 * exp (- E / 4e6 / 26 )
329- + 0.0000025 * E ** - 2.2 * exp (- E / 1.3e9 / 26 ),
337+ lambda E : (
338+ 0.01 * E ** - 2.40 * exp (- E / 1.2e5 / 7 )
339+ + 0.0006 * E ** - 2.3 * exp (- E / 4e6 / 7 )
340+ + 0.013 * E ** - 2.40 * exp (- E / 1.2e5 / 13 )
341+ + 0.0007 * E ** - 2.3 * exp (- E / 4e6 / 13 )
342+ ),
343+ lambda E : (
344+ 0.006 * E ** - 2.30 * exp (- E / 1.2e5 / 26 )
345+ + 0.00023 * E ** - 2.2 * exp (- E / 4e6 / 26 )
346+ + 0.0000025 * E ** - 2.2 * exp (- E / 1.3e9 / 26 )
347+ ),
330348 )
331349
332350
0 commit comments