diff options
Diffstat (limited to 'src/muon.c')
-rw-r--r-- | src/muon.c | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -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; |