diff --git a/monte_carlo.py b/monte_carlo.py
index b87acda5390dfb169446b25d9bad4c25c3dfab39..f62af7ee5d14ffe9dce7d810dc0d4733ed9e9caf 100644
--- a/monte_carlo.py
+++ b/monte_carlo.py
@@ -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
-nsources = 2
 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
 band = '1mm'
 
 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/"
-                       "HLS091828_common_mode_one_block/v_1")
+# DATA_DIR_SERVER = Path("/data/NIKA/Reduced/"
+#                        "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")
 
 
@@ -223,19 +230,6 @@ for _flux in flux:
         jackknifes.append(jk_iter())
     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
 
     helpfunc = partial(fake_worker, **{'min_threshold':
@@ -252,7 +246,7 @@ for _flux in flux:
     else:
         print('simulation without source injection')
 
-    if(0):
+    if(1):
         res = p.map(helpfunc, jackknifes)
         res = list(zip(*res))