diff --git a/regimes.py b/regimes.py
index 7392662..d9b3730 100644
--- a/regimes.py
+++ b/regimes.py
@@ -28,7 +28,7 @@ sims_path = data_path / 'simulations'
 
 trials = 1000
 spec_trials = 100   # set to zero to only recompute firng rates
-sigma = 0.002
+sigma = 0.001
 nfft = 2**18
 
 recompute = False
@@ -268,18 +268,18 @@ def plot_psd(ax, s, path, contrast, spikes, nfft, dt, beatf1, beatf2, eodf):
         ax.plot(beatf2 + 2*beatf1, decibel(peak_ampl(freqs, psd, beatf2 + 2*beatf1)) + offsm,
                 label=r'$\Delta f_2 \pm k\Delta f_1$', zorder=40, **s.psF012m)
         ax.plot(beatf2 + 3*beatf1, decibel(peak_ampl(freqs, psd, beatf2 + 3*beatf1)) + offsm, zorder=40, **s.psF012m)
-        ax.plot(beatf2 - 2*beatf1, decibel(peak_ampl(freqs, psd, beatf2 - 2*beatf1)) + offsm, zorder=40, **s.psF012m)
-        ax.plot(beatf2 - 3*beatf1, decibel(peak_ampl(freqs, psd, beatf2 - 3*beatf1)) + offsm, zorder=40, **s.psF012m)
+        #ax.plot(beatf2 - 2*beatf1, decibel(peak_ampl(freqs, psd, beatf2 - 2*beatf1)) + offsm, zorder=40, **s.psF012m)
+        #ax.plot(beatf2 - 3*beatf1, decibel(peak_ampl(freqs, psd, beatf2 - 3*beatf1)) + offsm, zorder=40, **s.psF012m)
         ax.plot(3*beatf1, decibel(peak_ampl(freqs, psd, 3*beatf1)) + offsm, zorder=40, **s.psF01m)
         ax.plot(4*beatf1, decibel(peak_ampl(freqs, psd, 4*beatf1)) + offsm, zorder=40, **s.psF01m)
         ax.plot(eodf - 2*beatf1, decibel(peak_ampl(freqs, psd, eodf - 2*beatf1)) + offsm, zorder=40, **s.psFEODm)
         ax.plot(eodf - 3*beatf1, decibel(peak_ampl(freqs, psd, eodf - 3*beatf1)) + offsm, zorder=40, **s.psFEODm)
-        ax.plot(eodf - 4*beatf1, decibel(peak_ampl(freqs, psd, eodf - 4*beatf1)) + offsm, zorder=40, **s.psFEODm)
+        #ax.plot(eodf - 4*beatf1, decibel(peak_ampl(freqs, psd, eodf - 4*beatf1)) + offsm, zorder=40, **s.psFEODm)
         ax.plot(eodf - 2*beatf2, decibel(peak_ampl(freqs, psd, eodf - 2*beatf2)) + offsm, zorder=40, **s.psF0m)
-        ax.plot(eodf - beatf2 + 2*beatf1, decibel(peak_ampl(freqs, psd, eodf - beatf2 + 2*beatf1)) + offsm,
-                zorder=40, **s.psF02m)
-        ax.plot(eodf - beatf2 + 3*beatf1, decibel(peak_ampl(freqs, psd, eodf - beatf2 + 2*beatf1)) + offsm,
-                zorder=40, **s.psF02m)
+        #ax.plot(eodf - beatf2 + 2*beatf1, decibel(peak_ampl(freqs, psd, eodf - beatf2 + 2*beatf1)) + offsm,
+        #        zorder=40, **s.psF02m)
+        #ax.plot(eodf - beatf2 + 3*beatf1, decibel(peak_ampl(freqs, psd, eodf - beatf2 + 2*beatf1)) + offsm,
+        #        zorder=40, **s.psF02m)
     ax.set_xlim(0, 750)
     ax.set_ylim(-60, 0)
     ax.set_xticks_delta(200)
@@ -290,8 +290,8 @@ def plot_psd(ax, s, path, contrast, spikes, nfft, dt, beatf1, beatf2, eodf):
 
 def plot_example(axs, axr, axf, axp, s, path, cell, alpha, beatf1, beatf2, eodf,
                  nfft, trials):
-    spec_path = path.with_name(path.stem + f'-contrastspectrum-{1000*alpha:03.0f}.npz')
-    rate_path = path.with_name(path.stem + f'-contrastrates-{1000*alpha:03.0f}.npz')
+    spec_path = path.with_name(path.stem + f'-contrastspectrum-c{1000*alpha:03.0f}.npz')
+    rate_path = path.with_name(path.stem + f'-contrastrates-c{1000*alpha:03.0f}-s{1000*sigma:.0f}.npz')
     dt = 0.0001
     tmax = nfft*dt
     t0 = 0.112