From c92657e33b81e37eb5851217fcdea8767c183e62 Mon Sep 17 00:00:00 2001
From: tlatorre <tlatorre@uchicago.edu>
Date: Tue, 26 Mar 2019 16:11:40 -0500
Subject: add energy dependent number of shower photons

This commit updates the code to calculate the number of Cerenkov photons from
secondary particles produced in an electromagnetic shower from electrons to use
an energy dependent formula I fit to data simulated with RAT-PAC.
---
 src/likelihood.c | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

(limited to 'src/likelihood.c')

diff --git a/src/likelihood.c b/src/likelihood.c
index f0430cf..c53e1fc 100644
--- a/src/likelihood.c
+++ b/src/likelihood.c
@@ -107,7 +107,7 @@ particle *particle_init(int id, double T0, size_t n)
          * range here using the dE/dx table instead of reading in the range. */
         p->T[n-1] = 0;
 
-        p->shower_photons = rad*ELECTRON_PHOTONS_PER_MEV;
+        p->shower_photons = electron_get_shower_photons(T0, rad);
 
         break;
     case IDP_MU_MINUS:
-- 
cgit