Commit eae82af5 authored by LUSTIG Peter's avatar LUSTIG Peter

new data directory

parent 5a6b3e6a
...@@ -139,27 +139,34 @@ def CreateParallelJackkifes(jkiter, parity_threshold): ...@@ -139,27 +139,34 @@ def CreateParallelJackkifes(jkiter, parity_threshold):
############################################################## ##############################################################
np.random.seed(0)
flux = np.array([10]) * u.mJy
nsim = 3
min_detection_threshold = 4 min_detection_threshold = 4
nsources = 2
photometry = True photometry = True
detection = False detection = True
ncores = 2 flux = np.geomspace(0.1, 10, 50) * u.mJy
nsim = 1000
ncores = 24
nsources = 50
parity_threshold = 1 parity_threshold = 1
band = '1mm' band = '1mm'
outdir = Path("montecarlo_results/") outdir = Path("montecarlo_results/")
outdir = Path("testdir")
filenamecomment = '' if 1:
np.random.seed(0)
flux = np.array([10]) * u.mJy
nsim = 1
ncores = 2
nsources = 2
outdir = Path("testdir")
filenamecomment = 'NewRealisation_'
############################################################## ##############################################################
DATA_DIR_SERVER = Path("/data/NIKA/Reduced/" # DATA_DIR_SERVER = Path("/data/NIKA/Reduced/"
"HLS091828_common_mode_one_block/v_1") # "HLS091828_common_mode_one_block/v_1")
DATA_DIR_SERVER = ("/data/NIKA/Reduced/May2018/"
"HLS091828_common_mode_one_block/v_1/")
DATA_DIR_MYPC = Path("/home/peter/Dokumente/Uni/Paris/Stage/data/v_1") DATA_DIR_MYPC = Path("/home/peter/Dokumente/Uni/Paris/Stage/data/v_1")
...@@ -223,19 +230,6 @@ for _flux in flux: ...@@ -223,19 +230,6 @@ for _flux in flux:
jackknifes.append(jk_iter()) jackknifes.append(jk_iter())
print('Done in {:.2f}s'.format(clock()-t0)) print('Done in {:.2f}s'.format(clock()-t0))
'''
# %%
# %matplotlib tk
plt.figure()
jackknifes[0].plot()
plt.figure()
jackknifes[0].match_filter(jackknifes[0].beam).plot()
plt.show(block=True)
# %%
sys.exit()
# set forced minimal distance of two fake sources
'''
min_source_dist = 2 * jackknifes[0].beam.fwhm_pix.value min_source_dist = 2 * jackknifes[0].beam.fwhm_pix.value
helpfunc = partial(fake_worker, **{'min_threshold': helpfunc = partial(fake_worker, **{'min_threshold':
...@@ -252,7 +246,7 @@ for _flux in flux: ...@@ -252,7 +246,7 @@ for _flux in flux:
else: else:
print('simulation without source injection') print('simulation without source injection')
if(0): if(1):
res = p.map(helpfunc, jackknifes) res = p.map(helpfunc, jackknifes)
res = list(zip(*res)) res = list(zip(*res))
......
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