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.

1 files changed, 0 insertions, 2 deletions

diff --git a/src/fit.c b/src/fit.c
index 23eab74..0634116 100644
--- a/src/fit.c
+++ b/src/fit.c
@@ -5619,7 +5619,6 @@ int main(int argc, char **argv)
 end:
 
     free_interpolation();
-    free_charge();
     pmt_response_free();
 
     db_free(db);
@@ -5632,7 +5631,6 @@ end:
 
 err:
     free_interpolation();
-    free_charge();
     pmt_response_free();
 
     db_free(db);