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/proton.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/proton.h')
| -rw-r--r-- | src/proton.h | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/src/proton.h b/src/proton.h index 008f835..cffeb9f 100644 --- a/src/proton.h +++ b/src/proton.h @@ -3,6 +3,15 @@ #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. + * + * FIXME: Actually determine what this is. */ +#define PROTON_PHOTONS_PER_MEV 400.0 + double proton_get_range(double T, double rho); double proton_get_dEdx_rad(double T, double rho); double proton_get_dEdx(double T, double rho); |