From feaea6075a06ece119a8d2829118d988c1feb06b Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Tue, 23 Jan 2024 14:20:07 -0700 Subject: [PATCH 01/13] fixing 4-comp fluxes for GST and GSF --- src/simweights/_fluxes.py | 25 +++++++++++++++++++++++++ 1 file changed, 25 insertions(+) diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index a7fb35a..0969c4b 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -310,6 +310,18 @@ class GlobalFitGST(CosmicRayFlux): ) +class GlobalFitGST_IT(CosmicRayFlux): + r"""GlobalFitGST for four components [p, He, O, Fe]. + The Oxygen group is the sum of Nitrogen and Aluminum groups of GlobalFitGST.""" + + pdgids = PDGID_4COMP + _funcs = ( + 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), + lambda E: 0.32 * E**-2.58 * exp(-E / 1.2e5 / 2) + 0.0065 * E**-2.3 * exp(-E / 4e6 / 2), + lambda E: 0.01 * E**-2.40 * exp(-E / 1.2e5 / 7) + 0.0006 * E**-2.3 * exp(-E / 4e6 / 7) +0.013 * E**-2.40 * exp(-E / 1.2e5 / 13) + 0.0007 * E**-2.3 * exp(-E / 4e6 / 13), + lambda E: 0.006 * E**-2.30 * exp(-E / 1.2e5 / 26)+ 0.00023 * E**-2.2 * exp(-E / 4e6 / 26) + 0.0000025 * E**-2.2 * exp(-E / 1.3e9 / 26), + ) + class GlobalSplineFitBase(CosmicRayFlux): r"""Data-driven spline fit of the cosmic ray spectrum by Dembinski et. al. \ [#GSFDembinski]. @@ -352,6 +364,19 @@ def __init__(self: GlobalSplineFit5Comp) -> None: super().__init__() +class GlobalSplineFit_IT(GlobalSplineFitBase): + r"""Sum of the flux of the GSF model for the standard 4 components injected by IceCube. + [(H), (He), (Li, Be, B, C, N, O, F, Ne), (Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni)] + GSF is a Data-driven spline fit of the cosmic ray spectrum by Dembinski et. al. \ [#GSFDembinski]. + """ + + pdgids = PDGID_4COMP + groups = ((1, 1), (2, 2), (3, 10), (11, 28)) + + def __init__(self: GlobalSplineFit4Comp_IT) -> None: + super().__init__() + + class FixedFractionFlux(CosmicRayFlux): """Total energy per particle flux flux split among mass groups with a constant fraction. From 754260865d3152ba0cf41d35c85c193d4411eebe Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Tue, 23 Jan 2024 14:23:50 -0700 Subject: [PATCH 02/13] fixing 4-comp fluxes for GST and GSF --- src/simweights/_fluxes.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index 0969c4b..cf22ad9 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -373,7 +373,7 @@ class GlobalSplineFit_IT(GlobalSplineFitBase): pdgids = PDGID_4COMP groups = ((1, 1), (2, 2), (3, 10), (11, 28)) - def __init__(self: GlobalSplineFit4Comp_IT) -> None: + def __init__(self: GlobalSplineFit_IT) -> None: super().__init__() From fbf081b5c1c9ae9a5808ae635af4844b40807848 Mon Sep 17 00:00:00 2001 From: "pre-commit-ci[bot]" <66853113+pre-commit-ci[bot]@users.noreply.github.com> Date: Tue, 23 Jan 2024 21:35:43 +0000 Subject: [PATCH 03/13] [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --- src/simweights/_fluxes.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index cf22ad9..231ddc1 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -318,10 +318,16 @@ class GlobalFitGST_IT(CosmicRayFlux): _funcs = ( 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), lambda E: 0.32 * E**-2.58 * exp(-E / 1.2e5 / 2) + 0.0065 * E**-2.3 * exp(-E / 4e6 / 2), - lambda E: 0.01 * E**-2.40 * exp(-E / 1.2e5 / 7) + 0.0006 * E**-2.3 * exp(-E / 4e6 / 7) +0.013 * E**-2.40 * exp(-E / 1.2e5 / 13) + 0.0007 * E**-2.3 * exp(-E / 4e6 / 13), - lambda E: 0.006 * E**-2.30 * exp(-E / 1.2e5 / 26)+ 0.00023 * E**-2.2 * exp(-E / 4e6 / 26) + 0.0000025 * E**-2.2 * exp(-E / 1.3e9 / 26), + lambda E: 0.01 * E**-2.40 * exp(-E / 1.2e5 / 7) + + 0.0006 * E**-2.3 * exp(-E / 4e6 / 7) + + 0.013 * E**-2.40 * exp(-E / 1.2e5 / 13) + + 0.0007 * E**-2.3 * exp(-E / 4e6 / 13), + lambda E: 0.006 * E**-2.30 * exp(-E / 1.2e5 / 26) + + 0.00023 * E**-2.2 * exp(-E / 4e6 / 26) + + 0.0000025 * E**-2.2 * exp(-E / 1.3e9 / 26), ) + class GlobalSplineFitBase(CosmicRayFlux): r"""Data-driven spline fit of the cosmic ray spectrum by Dembinski et. al. \ [#GSFDembinski]. From 5cc2fa7db14b05477b6eda61d73e734d734a743c Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 14:39:54 -0700 Subject: [PATCH 04/13] progress on #28 --- tests/flux_values.json | 44 ++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 44 insertions(+) diff --git a/tests/flux_values.json b/tests/flux_values.json index 81874c2..a80514f 100644 --- a/tests/flux_values.json +++ b/tests/flux_values.json @@ -129,6 +129,28 @@ 2.2713129512457712e-18, 3.871036212755032e-22, 1.8597326553441085e-24 ] }, + "GlobalFitGST_IT": { + "2212": [ + 0.03607536338982165, 8.303604540195788e-5, 1.8863983153209037e-7, + 3.1276596162531363e-10, 5.209459034660495e-13, 4.1532872002068553e-16, + 8.179387207124814e-19, 1.6295037238885218e-21, 6.3885346081344595e-27 + ], + "1000020040": [ + 0.023762074349477715, 6.617029782920994e-5, 1.878720144259654e-7, + 4.685738192425947e-10, 9.255593536028597e-13, 1.4792625219449634e-15, + 9.643447774808006e-23, 6.700436854786504e-74, 0.0 + ], + "1000080160": [ + 0.003971477965627884, 1.6135814814332014e-05, 6.546771674170436e-08 , + 2.5170413550201714e-10, 5.942565425469719e-13, 7.956149842340593e-16, + 4.744600882605979e-19, 6.212710347575392e-29, 2.1228946532983858e-106 + ], + "1000260560": [ + 0.0015996433832794234, 8.135139888332144e-6, 4.1420823134587994e-8, + 2.0697967478899172e-10, 8.354760133167253e-13, 1.0349313641570426e-15, + 2.2713129512457712e-18, 3.871036212755032e-22, 1.8597326553441085e-24 + ] + }, "Hoerandel": { "2212": [ 0.04477279875079687, 8.72999889028403e-5, 1.7021971482865675e-7, @@ -534,5 +556,27 @@ 1.5683732598520708e-10, 2.262722571198835e-13, 2.9271639177164597e-16, 5.2555897942362e-19, 1.8275958248865348e-22, 6.139956685254725e-28 ] + }, + "GlobalSplineFit_IT": { + "2212": [ + 0.04585449431423036, 7.789783588789647e-5, 1.9061451942412818e-7, + 3.913540332159069e-10, 9.171060628426788e-13, 7.933042844441958e-16, + 6.150944574738882e-19, 1.1644454409750712e-21, 3.79633922979983e-25 + ], + "1000020040": [ + 0.03338991785843059, 6.606992639645244e-5, 1.7771057235851816e-7, + 4.2287509663210625e-10, 8.159494882960421e-13, 9.87360959144915e-16, + 5.556183684902746e-19, 5.388863981471748e-22, 1.3504679127258198e-24 + ], + "1000080160": [ + 0.023110320325878606, 5.963761160902274e-05, 1.1796137928231345e-07, + 3.924962074713066e-10, 7.804815544960641e-13, 1.0230915905095007e-15, + 1.3687512894199914e-18, 4.294873263364232e-22, 1.1827498098295588e-25 + ], + "1000260560": [ + 0.0057939343800136305, 2.9715655250072123e-5, 5.976220938364594e-8, + 1.5683732598520708e-10, 2.262722571198835e-13, 2.9271639177164597e-16, + 5.2555897942362e-19, 1.8275958248865348e-22, 6.139956685254725e-28 + ] } } From f455fc60825d13625d48f7a13bed14d577b216fe Mon Sep 17 00:00:00 2001 From: "pre-commit-ci[bot]" <66853113+pre-commit-ci[bot]@users.noreply.github.com> Date: Wed, 24 Jan 2024 21:40:53 +0000 Subject: [PATCH 05/13] [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --- tests/flux_values.json | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tests/flux_values.json b/tests/flux_values.json index a80514f..f7030ab 100644 --- a/tests/flux_values.json +++ b/tests/flux_values.json @@ -141,7 +141,7 @@ 9.643447774808006e-23, 6.700436854786504e-74, 0.0 ], "1000080160": [ - 0.003971477965627884, 1.6135814814332014e-05, 6.546771674170436e-08 , + 0.003971477965627884, 1.6135814814332014e-5, 6.546771674170436e-8, 2.5170413550201714e-10, 5.942565425469719e-13, 7.956149842340593e-16, 4.744600882605979e-19, 6.212710347575392e-29, 2.1228946532983858e-106 ], @@ -569,7 +569,7 @@ 5.556183684902746e-19, 5.388863981471748e-22, 1.3504679127258198e-24 ], "1000080160": [ - 0.023110320325878606, 5.963761160902274e-05, 1.1796137928231345e-07, + 0.023110320325878606, 5.963761160902274e-5, 1.1796137928231345e-7, 3.924962074713066e-10, 7.804815544960641e-13, 1.0230915905095007e-15, 1.3687512894199914e-18, 4.294873263364232e-22, 1.1827498098295588e-25 ], From 922910428b5f6b2e9ef90bb1d20ebf474b0d8317 Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 14:49:08 -0700 Subject: [PATCH 06/13] progress on #28 --- src/simweights/_fluxes.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index 231ddc1..f2ec5ec 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -312,7 +312,9 @@ class GlobalFitGST(CosmicRayFlux): class GlobalFitGST_IT(CosmicRayFlux): r"""GlobalFitGST for four components [p, He, O, Fe]. - The Oxygen group is the sum of Nitrogen and Aluminum groups of GlobalFitGST.""" + + The Oxygen group is the sum of Nitrogen and Aluminum groups of GlobalFitGST. + """ pdgids = PDGID_4COMP _funcs = ( @@ -372,6 +374,7 @@ def __init__(self: GlobalSplineFit5Comp) -> None: class GlobalSplineFit_IT(GlobalSplineFitBase): r"""Sum of the flux of the GSF model for the standard 4 components injected by IceCube. + [(H), (He), (Li, Be, B, C, N, O, F, Ne), (Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni)] GSF is a Data-driven spline fit of the cosmic ray spectrum by Dembinski et. al. \ [#GSFDembinski]. """ From 0d4494db1a6c7cbdc609d8f71ddd9cec21ac4ca0 Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 15:50:20 -0700 Subject: [PATCH 07/13] progress on #28 --- src/simweights/__init__.py | 2 ++ src/simweights/_fluxes.py | 4 ++-- 2 files changed, 4 insertions(+), 2 deletions(-) diff --git a/src/simweights/__init__.py b/src/simweights/__init__.py index 0d8ffe8..8fcf835 100644 --- a/src/simweights/__init__.py +++ b/src/simweights/__init__.py @@ -56,6 +56,8 @@ Hoerandel5, Hoerandel_IT, Honda2004, + GlobalFitGST_IT, + GlobalSplineFit_IT, ) from ._generation_surface import GenerationSurface, generation_surface from ._genie_weighter import GenieWeighter diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index f2ec5ec..3b1ecea 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -310,7 +310,7 @@ class GlobalFitGST(CosmicRayFlux): ) -class GlobalFitGST_IT(CosmicRayFlux): +class GlobalFitGST_IT(CosmicRayFlux): #pylint: disable=invalid-name r"""GlobalFitGST for four components [p, He, O, Fe]. The Oxygen group is the sum of Nitrogen and Aluminum groups of GlobalFitGST. @@ -372,7 +372,7 @@ def __init__(self: GlobalSplineFit5Comp) -> None: super().__init__() -class GlobalSplineFit_IT(GlobalSplineFitBase): +class GlobalSplineFit_IT(GlobalSplineFitBase): #pylint: disable=invalid-name r"""Sum of the flux of the GSF model for the standard 4 components injected by IceCube. [(H), (He), (Li, Be, B, C, N, O, F, Ne), (Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni)] From 9923aae32b04ee59c8bbda527dce411ae0fda3d6 Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 15:55:23 -0700 Subject: [PATCH 08/13] progress on #28 --- src/simweights/__init__.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/src/simweights/__init__.py b/src/simweights/__init__.py index 8fcf835..a5a8dc7 100644 --- a/src/simweights/__init__.py +++ b/src/simweights/__init__.py @@ -51,7 +51,9 @@ GaisserH4a_IT, GaisserHillas, GlobalFitGST, + GlobalFitGST_IT, GlobalSplineFit5Comp, + GlobalSplineFit_IT, Hoerandel, Hoerandel5, Hoerandel_IT, From 0b15640f5dbfe1e49f0ad2d3d705e0635665a14c Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 16:05:56 -0700 Subject: [PATCH 09/13] progress on #28 --- src/simweights/__init__.py | 4 ++-- src/simweights/_fluxes.py | 4 ++-- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/src/simweights/__init__.py b/src/simweights/__init__.py index a5a8dc7..66c631d 100644 --- a/src/simweights/__init__.py +++ b/src/simweights/__init__.py @@ -27,7 +27,9 @@ "GaisserH4a_IT", "GaisserHillas", "GlobalFitGST", + "GlobalFitGST_IT", "GlobalSplineFit5Comp", + "GlobalSplineFit_IT", "Hoerandel", "Hoerandel5", "Hoerandel_IT", @@ -58,8 +60,6 @@ Hoerandel5, Hoerandel_IT, Honda2004, - GlobalFitGST_IT, - GlobalSplineFit_IT, ) from ._generation_surface import GenerationSurface, generation_surface from ._genie_weighter import GenieWeighter diff --git a/src/simweights/_fluxes.py b/src/simweights/_fluxes.py index 3b1ecea..5207960 100644 --- a/src/simweights/_fluxes.py +++ b/src/simweights/_fluxes.py @@ -310,7 +310,7 @@ class GlobalFitGST(CosmicRayFlux): ) -class GlobalFitGST_IT(CosmicRayFlux): #pylint: disable=invalid-name +class GlobalFitGST_IT(CosmicRayFlux): # pylint: disable=invalid-name r"""GlobalFitGST for four components [p, He, O, Fe]. The Oxygen group is the sum of Nitrogen and Aluminum groups of GlobalFitGST. @@ -372,7 +372,7 @@ def __init__(self: GlobalSplineFit5Comp) -> None: super().__init__() -class GlobalSplineFit_IT(GlobalSplineFitBase): #pylint: disable=invalid-name +class GlobalSplineFit_IT(GlobalSplineFitBase): # pylint: disable=invalid-name r"""Sum of the flux of the GSF model for the standard 4 components injected by IceCube. [(H), (He), (Li, Be, B, C, N, O, F, Ne), (Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni)] From c3e5cff4f363f7f158d8dada6cde94344adb80b9 Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 17:10:04 -0700 Subject: [PATCH 10/13] progress on #28 --- tests/test_fluxes.py | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/tests/test_fluxes.py b/tests/test_fluxes.py index de91aad..97eb5df 100755 --- a/tests/test_fluxes.py +++ b/tests/test_fluxes.py @@ -58,6 +58,9 @@ def test_TIG1996(self): def test_GlobalFitGST(self): self.flux_cmp("GlobalFitGST") + + def test_GlobalFitGST_IT(self): + self.flux_cmp("GlobalFitGST_IT") def test_GlobalSplineFit(self): self.flux_cmp("GlobalSplineFit") @@ -65,6 +68,9 @@ def test_GlobalSplineFit(self): def test_GlobalSplineFit5Comp(self): self.flux_cmp("GlobalSplineFit5Comp") + def test_GlobalSplineFit_IT(self): + self.flux_cmp("GlobalSplineFit_IT") + def test_FixedFractionFlux(self): self.flux_cmp("FixedFractionFlux", {2212: 0.1, 1000020040: 0.2, 1000080160: 0.3, 1000260560: 0.4}) self.flux_cmp( From 9d79e937fb1d73183246e748a0bdeb0d8e38cbf3 Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 17:15:58 -0700 Subject: [PATCH 11/13] progress on #28 --- tests/test_fluxes.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tests/test_fluxes.py b/tests/test_fluxes.py index 97eb5df..e551ed0 100755 --- a/tests/test_fluxes.py +++ b/tests/test_fluxes.py @@ -58,7 +58,7 @@ def test_TIG1996(self): def test_GlobalFitGST(self): self.flux_cmp("GlobalFitGST") - + def test_GlobalFitGST_IT(self): self.flux_cmp("GlobalFitGST_IT") @@ -70,7 +70,7 @@ def test_GlobalSplineFit5Comp(self): def test_GlobalSplineFit_IT(self): self.flux_cmp("GlobalSplineFit_IT") - + def test_FixedFractionFlux(self): self.flux_cmp("FixedFractionFlux", {2212: 0.1, 1000020040: 0.2, 1000080160: 0.3, 1000260560: 0.4}) self.flux_cmp( From 9bda715a27ecaa4baaf7fb9b4d05aee16ff0c91d Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 21:37:58 -0700 Subject: [PATCH 12/13] progress on #28 --- tests/flux_values.json | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/tests/flux_values.json b/tests/flux_values.json index f7030ab..5043bf5 100644 --- a/tests/flux_values.json +++ b/tests/flux_values.json @@ -564,19 +564,19 @@ 6.150944574738882e-19, 1.1644454409750712e-21, 3.79633922979983e-25 ], "1000020040": [ - 0.03338991785843059, 6.606992639645244e-5, 1.7771057235851816e-7, - 4.2287509663210625e-10, 8.159494882960421e-13, 9.87360959144915e-16, - 5.556183684902746e-19, 5.388863981471748e-22, 1.3504679127258198e-24 + 0.03192055738163634, 6.467091543445148e-5, 1.7585497031971908e-7, + 4.1712912655120846e-10, 8.035956231065464e-13, 9.732708568816786e-16, + 5.416659594465701e-19, 5.335274557933634e-22, 1.3445890398933065e-24 ], "1000080160": [ - 0.023110320325878606, 5.963761160902274e-5, 1.1796137928231345e-7, - 3.924962074713066e-10, 7.804815544960641e-13, 1.0230915905095007e-15, - 1.3687512894199914e-18, 4.294873263364232e-22, 1.1827498098295588e-25 + 0.018496228648837695, 4.259693993066384e-05, 8.770466269694322e-08, + 2.6164362595183446e-10, 6.022544255316352e-13, 7.62886294617275e-16, + 9.38795530909692e-19, 3.8207698885288526e-22, 1.1392141934461967e-25 ], "1000260560": [ - 0.0057939343800136305, 2.9715655250072123e-5, 5.976220938364594e-8, - 1.5683732598520708e-10, 2.262722571198835e-13, 2.9271639177164597e-16, - 5.2555897942362e-19, 1.8275958248865348e-22, 6.139956685254725e-28 + 0.012020219553048266, 4.9354088622501694e-05, 9.361411920367124e-08, + 2.9805554501291447e-10, 4.2376209855749804e-13, 5.757781836050234e-16, + 9.852668748367386e-19, 2.4213191949850687e-22, 1.0847379951951749e-26 ] } } From 069df80f711c4947574bc96a5325565cbe4da65b Mon Sep 17 00:00:00 2001 From: Amar Thakuri Date: Wed, 24 Jan 2024 21:42:21 -0700 Subject: [PATCH 13/13] progress on #28 --- tests/flux_values.json | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tests/flux_values.json b/tests/flux_values.json index 5043bf5..4ad0752 100644 --- a/tests/flux_values.json +++ b/tests/flux_values.json @@ -569,12 +569,12 @@ 5.416659594465701e-19, 5.335274557933634e-22, 1.3445890398933065e-24 ], "1000080160": [ - 0.018496228648837695, 4.259693993066384e-05, 8.770466269694322e-08, + 0.018496228648837695, 4.259693993066384e-5, 8.770466269694322e-8, 2.6164362595183446e-10, 6.022544255316352e-13, 7.62886294617275e-16, 9.38795530909692e-19, 3.8207698885288526e-22, 1.1392141934461967e-25 ], "1000260560": [ - 0.012020219553048266, 4.9354088622501694e-05, 9.361411920367124e-08, + 0.012020219553048266, 4.9354088622501694e-5, 9.361411920367124e-8, 2.9805554501291447e-10, 4.2376209855749804e-13, 5.757781836050234e-16, 9.852668748367386e-19, 2.4213191949850687e-22, 1.0847379951951749e-26 ]