diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2019-01-27 21:08:25 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2019-01-27 21:08:25 -0600 |
| commit | 1d77bacaae25d40d160f2bcd14ba3a355921213e (patch) | |
| tree | ad283c9c0d7bfd5326690bc43af7c82c5a1bb40d /src/electron.h | |
| parent | b9491718282f86b77c2594f161b096903706edc1 (diff) | |
| download | sddm-1d77bacaae25d40d160f2bcd14ba3a355921213e.tar.gz sddm-1d77bacaae25d40d160f2bcd14ba3a355921213e.tar.bz2 sddm-1d77bacaae25d40d160f2bcd14ba3a355921213e.zip | |
add photons from delta rays to likelihood calculation
This commit updates the likelihood function to take into account Cerenkov light
produced from delta rays produced by muons. The angular distribution of this
light is currently assumed to be constant along the track and parameterized in
the same way as the Cerenkov light from an electromagnetic shower. Currently I
assume the light is produced uniformly along the track which isn't exactly
correct, but should be good enough.
Diffstat (limited to 'src/electron.h')
| -rw-r--r-- | src/electron.h | 8 |
1 files changed, 8 insertions, 0 deletions
diff --git a/src/electron.h b/src/electron.h index cb3f24e..20e5ea2 100644 --- a/src/electron.h +++ b/src/electron.h @@ -3,9 +3,17 @@ #include <stddef.h> /* for size_t */ +/* Number of photons in the range 200 nm - 800 nm generated per MeV of energy + * lost to radiation for electrons. + * + * FIXME: This is just a rough estimate, should use an energy dependent + * quantity from simulation. */ +#define ELECTRON_PHOTONS_PER_MEV 400.0 + void electron_get_position_distribution_parameters(double T0, double *a, double *b); double electron_get_angular_distribution_alpha(double T0); double electron_get_angular_distribution_beta(double T0); +double electron_get_angular_pdf_delta_ray(double cos_theta, double alpha, double beta, double mu); double electron_get_angular_pdf(double cos_theta, double alpha, double beta, double mu); double electron_get_range(double T, double rho); double electron_get_dEdx_rad(double T, double rho); |