diff --git a/reconstruction.py b/reconstruction.py
index 30cb87a790c7445af81fe0d0256da5c52886b40b..cd41c2a160780d17bd67e8a317fa5c759d12f283 100644
--- a/reconstruction.py
+++ b/reconstruction.py
@@ -29,23 +29,6 @@ from config import doms, eps, min_timestamps_per_bin, nthreads
 from config import sky_movie_mf as sky
 from config import startt, dt, npixt
 
-def stat_plotting(pos, KL):
-    if master:
-        sc = ift.StatCalculator()
-        sc_mean = ift.StatCalculator()
-        sc_spectral = ift.StatCalculator()
-        skysamples = {'hi': [], 'lo': []}
-
-    # CAUTION: the loop and the if-clause cannot be interchanged!
-    for ss in KL.samples:
-        if master:
-            samp = sky(pos+ss)
-            skysamples['hi'].append(samp['hi'])
-            skysamples['lo'].append(samp['lo'])
-            sc.add(samp)
-            sc_mean.add(0.5*(samp['hi']+samp['lo']).val)
-            sc_spectral.add(2 * (samp['lo']-samp['hi']).val / (samp['lo']+samp['hi']).val )
-
 def optimization_heuristic(ii, likelihoods):
     lh_full, lh_full_amp, lh_full_ph, lh_cut, lh_cut_amp, lh_cut_ph = likelihoods 
     N_samples = 10 * (1 + ii // 8)