add a function to compute the number of shower photons for muons

Similarly to electrons, I fit an analytic form to the ratio of the number of
photons produced via shower particles over the radiative energy loss. In this
case, I chose the functional form:

    ratio = a*(1-exp(-T/b))

since the ratio seemed to reach a constant value after a certain energy. I then
simulated a 10 GeV muon and it appears that the ratio might actually decrease
after that, so for higher energies I may have to come up with a different fit
function.

1 files changed, 1 insertions, 0 deletions

diff --git a/src/muon.h b/src/muon.h
index e259999..a8d2ea1 100644
--- a/src/muon.h
+++ b/src/muon.h
@@ -30,6 +30,7 @@
  * FIXME: Actually determine what this is. */
 #define MUON_PHOTONS_PER_MEV 7368.0
 
+double muon_get_shower_photons(double T0, double rad);
 void muon_get_position_distribution_parameters(double T0, double *a, double *b);
 double muon_get_angular_distribution_alpha(double T0);
 double muon_get_angular_distribution_beta(double T0);