diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2018-10-17 09:33:29 -0500 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2018-10-17 09:33:29 -0500 |
| commit | 643204e807d5e78f883fc30dc7383a209e86dbc5 (patch) | |
| tree | e9e09dec819cd9b7c45bf8081b0acb7a329d5611 /src/scattering.h | |
| parent | 09f7f3ec8bbff5102d0447ee664df3f3a404c9bc (diff) | |
| download | sddm-643204e807d5e78f883fc30dc7383a209e86dbc5.tar.gz sddm-643204e807d5e78f883fc30dc7383a209e86dbc5.tar.bz2 sddm-643204e807d5e78f883fc30dc7383a209e86dbc5.zip | |
fix a bug in the theta0 calculation for a path
This commit fixes a bug in the calculation of the average rms width of the
angular distribution for a path with a KL expansion. I also made a lot of
updates to the test-path program:
- plot the distribution of the KL expansion coefficients
- plot the standard deviation of the angular distribution as a function of
distance along with the prediction
- plot the simulated and reconstructed path in 3D
Diffstat (limited to 'src/scattering.h')
| -rw-r--r-- | src/scattering.h | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/src/scattering.h b/src/scattering.h index 308f8ef..65765b0 100644 --- a/src/scattering.h +++ b/src/scattering.h @@ -1,6 +1,11 @@ #ifndef SCATTERING_H #define SCATTERING_H +/* Maximum value for the scattering RMS `theta0`. + * + * This is to prevent the interpolation from going out of bounds. */ +#define MAX_THETA0 1.0 + void init_interpolation(void); double get_probability(double beta, double cos_theta, double theta0); double get_probability2(double beta); |