Merge pull request #106 from lsst-camera-dh/LSSTTD-1558_pd_integratio…

…n_bug treat each picoammeter point as the average of the pd current over the preceding interval
lsst-camera-dh · May 24, 2021 · 091d930 · 091d930
2 parents 60499a5 + 794a63d
commit 091d930
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/python/bot_eo_analyses.py b/python/bot_eo_analyses.py
@@ -128,7 +128,9 @@ def mondiode_value(flat_file, exptime, factor=5,
                                      pd_filename))[0]
     x, y = np.recfromtxt(pd_file).transpose()
 
-    integral = sum((y[1:] + y[:-1])/2*(x[1:] - x[:-1]))
+    # Each picoammeter measurement point is the average of the pd
+    # current over the preceding interval.
+    integral = sum(y[1:]*(x[1:] - x[:-1]))
     return integral/exptime