From f93d19b96093cdd498639868700e0f722b7eb9b0 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Wed, 28 Nov 2018 09:46:52 -0600 Subject: update sno_charge.c This commit adds lots of comments to sno_charge.c and makes a couple of other changes: - use interp1d() instead of the GSL interpolation routines - increase MAX_PE to 100 I increased MAX_PE because I determined that it had a rather large impact on the likelihood function for 500 MeV electrons. This unfortunately slows down the initialization by a lot. I think I could speed this up by convolving the single PE charge distribution with a gaussian *before* convolving the charge distributions to compute the charge distributions for multiple PE. --- src/sno_charge.h | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) (limited to 'src/sno_charge.h') diff --git a/src/sno_charge.h b/src/sno_charge.h index 50b018e..6f513f9 100644 --- a/src/sno_charge.h +++ b/src/sno_charge.h @@ -2,10 +2,9 @@ #define SNO_CHARGE_H void init_charge(void); -double log_pq(double q, int n); -double pq(double q, int n); +double get_log_pq(double q, int n); +double get_pq(double q, int n); double get_log_pmiss(int n); double spe_pol2exp(double q); -void free_charge(void); #endif -- cgit