Put number of detected photons into charge value for channel.

1 files changed, 2 insertions, 0 deletions

diff --git a/src/daq.cu b/src/daq.cu
index 9824d7d..2b95560 100644
--- a/src/daq.cu
+++ b/src/daq.cu
@@ -35,6 +35,7 @@ __global__ void run_daq(curandState *s, unsigned int detection_state,
 			unsigned int *photon_histories,
 			int *last_hit_triangles, int *solid_map,
 			int nsolids, unsigned int *earliest_time_int,
+			unsigned int *channel_q,
 			unsigned int *channel_histories)
 {
 
@@ -57,6 +58,7 @@ __global__ void run_daq(curandState *s, unsigned int detection_state,
 			  unsigned int time_int = float_to_sortable_int(time);
 			  
 			  atomicMin(earliest_time_int + solid_id, time_int);
+			  atomicAdd(channel_q + solid_id, 1);
 			  atomicOr(channel_histories + solid_id, history);
 			}