add rayleigh scattering

This commit adds Rayleigh scattering to the likelihood function. The Rayleigh
scattering lengths come from rsp_rayleigh.dat from SNOMAN which only includes
photons which scattered +/- 10 ns around the prompt peak. The fraction of light
which scatters is treated the same in the likelihood as reflected light, i.e.
it is uniform across all the PMTs in the detector and the time PDF is assumed
to be a constant for a fixed amount of time after the prompt peak.

1 files changed, 9 insertions, 1 deletions

diff --git a/src/optics.h b/src/optics.h
index 06b89ca..64796c6 100644
--- a/src/optics.h
+++ b/src/optics.h
@@ -1,7 +1,15 @@
 #ifndef OPTICS_H
 #define OPTICS_H
 
-int optics_init(void);
+#include "dict.h"
+
+extern char optics_err[256];
+
+int optics_init(dict *db);
+double get_rayleigh_scattering_length_snoman_d2o(double wavelength);
+double get_rayleigh_scattering_length_snoman_h2o(double wavelength);
+double get_weighted_rayleigh_scattering_length_snoman_d2o(void);
+double get_weighted_rayleigh_scattering_length_snoman_h2o(void);
 double get_weighted_absorption_length_snoman_h2o(void);
 double get_absorption_length_snoman_h2o(double wavelength);
 double get_weighted_absorption_length_snoman_d2o(void);