update test-find-peaks to plot cerenkov rings

This commit updates the test-find-peaks script to plot Cerenkov rings for each
of the peaks. It also updates the script to use quad to find the position
instead of using the MC information. Finally, I added a -n argument to the
script to specify how many peaks to draw.

1 files changed, 1 insertions, 0 deletions

diff --git a/src/misc.h b/src/misc.h
index 7dad51c..e06e323 100644
--- a/src/misc.h
+++ b/src/misc.h
@@ -28,6 +28,7 @@
 #define LNFACT_MAX 100
 
 double trapz(const double *y, double dx, size_t n);
+int intersect_sphere(double *pos, double *dir, double R, double *l);
 void get_path_length(double *pos1, double *pos2, double R, double *l1, double *l2);
 double ln(unsigned int n);
 double lnfact(unsigned int n);