diff options
| author | tlatorre <tlatorre@uchicago.edu> | 2020-02-10 09:37:06 -0600 |
|---|---|---|
| committer | tlatorre <tlatorre@uchicago.edu> | 2020-02-10 09:37:06 -0600 |
| commit | f41d2cfab0ebe9d8a41333b6a24cf44550ce68d2 (patch) | |
| tree | 33b06a6c0864655a9cbf8c05cff9cee297665ab2 /src | |
| parent | a2b3a8097c5ac3e2e74af2e10e75f6aa571990a8 (diff) | |
| download | sddm-f41d2cfab0ebe9d8a41333b6a24cf44550ce68d2.tar.gz sddm-f41d2cfab0ebe9d8a41333b6a24cf44550ce68d2.tar.bz2 sddm-f41d2cfab0ebe9d8a41333b6a24cf44550ce68d2.zip | |
fix formula for calculating b for position distribution of light from electron showers
Also update the a parameter based on a simple 0 degree polynomial fit to the
shower profiles above 100 MeV.
Diffstat (limited to 'src')
| -rw-r--r-- | src/electron.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/electron.c b/src/electron.c index 3f5d4d8..967026c 100644 --- a/src/electron.c +++ b/src/electron.c @@ -118,9 +118,9 @@ void electron_get_position_distribution_parameters(double T0, double *a, double * See http://pdg.lbl.gov/2014/reviews/rpp2014-rev-passage-particles-matter.pdf. */ double tmax; - *b = RADIATION_LENGTH; + *b = 43.51; tmax = log(T0/ELECTRON_CRITICAL_ENERGY_D2O) - 0.5; - *a = fmax(1.1,tmax*(*b)/RADIATION_LENGTH + 1); + *a = fmax(1.1,tmax*RADIATION_LENGTH/(*b) + 1); } double electron_get_angular_distribution_alpha(double T0) |