diff --git a/teaser_critical_filter.py b/teaser_critical_filter.py
index 1c13b7772aa6252cb25ed0b310e1c87bd4fd4e2f..0b57738e8f3eb97b73ea3d1e5f5a4ea7be786679 100644
--- a/teaser_critical_filter.py
+++ b/teaser_critical_filter.py
@@ -96,10 +96,8 @@ plot_WF('teaser_unknown_power', ground_truth, data, m=mean,
         samples=posterior_samples)
 
 # Plot the reconstruction of the power spectrum
-mysterious_spectrum = lambda k: 5/((7**2 - k**2)**2 + 3**2*k**2)
-# FIXME There is a simpler way, isn't it?
-ground_truth_spectrum = ift.makeField(power_space,
-                                      mysterious_spectrum(power_space.k_lengths))
+ground_truth_spectrum = ift.PS_field(power_space,
+                                     lambda k: 5/((7**2 - k**2)**2 + 3**2*k**2))
 posterior_power_samples = [pspec.force(KL.position+samp) for samp in KL.samples]
 power_mean = reduce(add, posterior_power_samples)
 power_plot('teaser_power_reconstruction', ground_truth_spectrum, power_mean,