From 9ffc374f6aea2c976e7ebdcfffd6d23eed03a45f Mon Sep 17 00:00:00 2001 From: John ZuHone Date: Thu, 22 Jun 2017 14:58:22 -0400 Subject: [PATCH] Use a poisson number here --- soxs/instrument.py | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/soxs/instrument.py b/soxs/instrument.py index c15f180a..62132c20 100644 --- a/soxs/instrument.py +++ b/soxs/instrument.py @@ -85,8 +85,7 @@ def detect_events(self, events, exp_time, flux, refband, prng=None): earea = self.interpolate_area(energy).value idxs = np.logical_and(energy >= refband[0], energy <= refband[1]) rate = flux/(energy[idxs].sum()*erg_per_keV)*earea[idxs].sum() - n_ph = np.modf(rate*exp_time) - n_ph = np.int64(n_ph[1]) + np.int64(n_ph[0] >= prng.uniform()) + n_ph = prng.poisson(lam=rate*exp_time) fak = float(n_ph)/energy.size if fak > 1.0: mylog.error("Number of events in sample: %d, Number of events wanted: %d" % (energy.size, n_ph))