don't need to divide out the extra four bits

lib/imaging.py

@@ -199,7 +199,7 @@ def compute_mot_number(image, QE, gain, exposure_time, solid_angle,
     H, W = image.shape
     im = image[H//3:2*H//3, W//3:2*W//3]
     pd.Plotter().imshow(im).savefig("img.png").close()
-    electron_rate = im.astype(np.float64).sum() / 16 # get to 12-bit values
+    electron_rate = im.astype(np.float64).sum()
     photon_rate \
             = electron_rate \
             / 10**(gain / 10) \