From 6e3af3bb59c729bbcfa1126a6c05694c52e21616 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Mon, 1 Oct 2018 16:07:33 -0500 Subject: use the PMT response table to calculate the amount of reflected light To calculate the expected number of photons from reflected light we now integrate over the track and use the PMT response table to calculate what fraction of the light is reflected. Previously we were just using a constant fraction of the total detected light which was faster since we only had to integrate over the track once, but this should be more accurate. --- src/likelihood.h | 9 ++------- 1 file changed, 2 insertions(+), 7 deletions(-) (limited to 'src/likelihood.h') diff --git a/src/likelihood.h b/src/likelihood.h index d53b777..b4b3be9 100644 --- a/src/likelihood.h +++ b/src/likelihood.h @@ -18,13 +18,8 @@ * value for n. */ #define MIN_RATIO -10 -/* Roughly equal to the fraction of light which is reflected at each PMT. - * - * Note: The proper thing to do here would be to separately integrate over the - * track to calculate the amount of reflected light using the data from the - * PMTR bank. However, this would make the likelihood function almost twice as - * long, so we just assume a constant fraction of the light is reflected now. */ -#define CHARGE_FRACTION 0.3 +/* The fraction of reflected light which is detected. */ +#define CHARGE_FRACTION 1.0 /* Dark rate of the PMTs (Hz). * -- cgit