aboutsummaryrefslogtreecommitdiff
path: root/src/muon.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/muon.c')
-rw-r--r--src/muon.c4
1 files changed, 2 insertions, 2 deletions
diff --git a/src/muon.c b/src/muon.c
index 6ab0aa7..452dfba 100644
--- a/src/muon.c
+++ b/src/muon.c
@@ -98,7 +98,7 @@ double muon_get_angular_distribution_alpha(double T0)
*
* where T0 is the initial energy of the muon in MeV and c0, c1, c2,
* and c3 are constants which I fit out. */
- return 8.238633e-01 + 3.896665e-03/log(1.581060e-05*T0 + 9.991576e-01);
+ return fmin(100.0,8.238633e-01 + 3.896665e-03/log(1.581060e-05*T0 + 9.991576e-01));
}
double muon_get_angular_distribution_beta(double T0)
@@ -144,7 +144,7 @@ void muon_get_delta_ray_distribution_parameters(double T0, double *a, double *b)
*
* where T0 is the initial energy of the muon in MeV and c0, c1, c2,
* and c3 are constants which I fit out. */
- *a = 3.463093e-01 + 1.110835e-02/log(5.662948e-06*T0 + 1.009215e+00);
+ *a = fmin(100.0,3.463093e-01 + 1.110835e-02/log(5.662948e-06*T0 + 1.009215e+00));
*b = 2.297358e-01 + 4.085721e-03/log(8.218774e-06*T0 + 1.007135e+00);
return;