add PMT response

This commit adds code to read in the PMT response from the PMTR bank from
SNOMAN. This file was used for the grey disk model in SNOMAN and was created
using a full 3D simulation of the PMT and concentrator. Since the PMT response
in SNOMAN included the quantum efficiency of the PMT, we have to divide that
out to get just the PMT response independent of the quantum efficiency.

I also updated the likelihood calculation to use the pmt response. Currently
the energy is being fit too high which I think will improve when we update the
solid angle calculation to use the radius of the concentrator instead of the
PMT.

1 files changed, 23 insertions, 0 deletions

diff --git a/src/pmt_response.h b/src/pmt_response.h
new file mode 100644
index 0000000..300f6b2
--- /dev/null
+++ b/src/pmt_response.h
@@ -0,0 +1,23 @@
+#ifndef PMT_RESPONSE_H
+#define PMT_RESPONSE_H
+
+#include "dict.h"
+
+extern char pmtr_err[256];
+
+#define NUM_ANGLES 90
+#define NUM_WAVELENGTHS 50
+
+#define KPMTR_TRTI 1
+#define KPMTR_DECA 2
+#define KPMTR_SPRD 3
+#define KPMTR_CERF 4
+#define KPMTR_RESP 5
+#define KPMTR_REFLEC 4506
+
+double get_weighted_pmt_response(double theta);
+double get_pmt_response(double wavelength, double theta);
+int pmt_response_init(dict *db);
+void pmt_response_free(void);
+
+#endif