From b0c6382a80ef9704858c3a77511fa58f7949d7a3 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Wed, 22 May 2019 12:44:15 -0400 Subject: 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. --- src/muon.h | 1 + 1 file changed, 1 insertion(+) (limited to 'src/muon.h') 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); -- cgit