def getShotNoise(myeff, effbins, samples, samplesize):
effs = numpy.zeros(effbins)
for _ in range(samples):
rnds = numpy.random.random(samplesize)
accs = (rnds < myeff).sum()
eff = efficiencyDeltaFix((accs / samplesize,))
ndx = genEffIndex(eff[0], effbins)
effs[ndx] += 1 / samples
return effs
def populateMatrixColumnWithBurst(self, rbinindex, reff, burst, weight):
effval = genRandomBurstEffs(reff, burst)
effval = efficiencyDeltaFix((effval, ))
effndx = genEffIndex(effval, self.EffBins)
self.tm[rbinindex, effndx] += weight
def populateMatrixColumnWithBurst(self, rbinindex, reff, burst, weight):
effval = genRandomBurstEffs(reff, burst)
effval = efficiencyDeltaFix((effval,))
effndx = genEffIndex(effval, self.EffBins)
self.tm[rbinindex, effndx] += weight
def test_deltaFix(self):
arr = numpy.random.rand(1000)
aavg = arr.mean()
self.assertAlmostEqual(aavg,
efficiencyDeltaFix(arr).mean(),
delta=0.01)
def test_deltaFix(self):
arr = numpy.random.rand(1000)
aavg = arr.mean()
self.assertAlmostEqual(aavg, efficiencyDeltaFix(arr).mean(), delta = 0.01)
