Commit 25ec5816 authored by LUSTIG Peter's avatar LUSTIG Peter

Added parity_threshold option

parent 261f74b9
......@@ -33,7 +33,8 @@ plt.ion()
def fake_worker(jkiter, min_threshold=2, nsources=8**2, flux=1*u.Jy,
within=(0, 1), cat_gen=pos_uniform, **kwargs):
within=(0, 1), cat_gen=pos_uniform, parity_threshold=1,
**kwargs):
"""The completness purity worker, create a fake dataset from an jackknifed
image and return catalogs
......@@ -46,7 +47,7 @@ def fake_worker(jkiter, min_threshold=2, nsources=8**2, flux=1*u.Jy,
**kwargs :
arguments for source injection
"""
img = jkiter()
img = jkiter(parity_threshold)
# img = IMAGE
# Renormalize the stddev
std = img.check_SNR()
......@@ -176,7 +177,8 @@ p = Pool(ncores)
helpfunc = partial(fake_worker, **{'min_threshold': min_detection_threshold,
'nsources': nsources, 'flux': flux,
'within': (0, 1), 'cat_gen': pos_uniform})
'within': (0, 1), 'cat_gen': pos_uniform,
'parity_threshold': 0.5})
# sources, fake_sources = helpfunc(jk_iter_list)
print('Begin Monte Carlo')
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment