Handle cases where more than 2 billion photons enter a pixel. (#188)

spacetelescope · Dec 6, 2024 · 72f0084 · 72f0084
1 parent 09f9d75
commit 72f0084
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/romanisim/image.py b/romanisim/image.py
@@ -470,7 +470,7 @@ def simulate_counts_generic(image, exptime, objlist=None, psf=None,
         rng_numpy_seed = rng.raw()
         rng_numpy = np.random.default_rng(rng_numpy_seed)
         image.array[:, :] = rng_numpy.binomial(
-            image.array.astype('i4'), flat / maxflat)
+            np.clip(image.array, 0, 2**31 - 1).astype('i4'), flat / maxflat)
 
     if dark is not None:
         workim = image * 0