From a0c1f4560325311558eedf39bd191ddaad182367 Mon Sep 17 00:00:00 2001 From: tlatorre Date: Tue, 27 Nov 2018 16:07:56 -0600 Subject: 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. --- src/optics.h | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) (limited to 'src/optics.h') 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); -- cgit