figure 4 addapted regime borders

regimes.py

@@ -10,11 +10,12 @@ from plotstyle import plot_style, lighter, darker
 
 
 model_cell = '2018-05-08-ad-invivo-1'      # 228Hz, CV=0.67
-alphas = [0.002, 0.01, 0.03, 0.06]
+#alphas = [0.002, 0.01, 0.03, 0.06]
+alphas = [0.002, 0.01, 0.025, 0.05]
 rmax = 500
 amax = 60
 cthresh1 = 1.2
-cthresh2 = 3.05
+cthresh2 = 2.6
 
 #model_cell = '2018-05-08-ab-invivo-1'      # 116, CV=0.68
 #alphas = [0.002, 0.008, 0.025, 0.05]
@@ -351,8 +352,7 @@ def amplitude_squarefit(contrast, power, max_contrast):
 
 
 def plot_peaks(ax, s, alphas, contrasts, powerf1, powerf2, powerfsum,
-               powerfdiff):
-    cmax = 10
+               powerfdiff, cmax=8):
     contrasts *= 100
     ax.plot(contrasts, amplitude_linearfit(contrasts, powerf1, cthresh1),
             **s.lsF01m)
@@ -391,7 +391,7 @@ def plot_peaks(ax, s, alphas, contrasts, powerf1, powerf2, powerfsum,
                 ha='center', va='bottom')
     ax.set_xlim(0, cmax)
     ax.set_ylim(0, amax)
-    ax.set_xticks_delta(2)
+    ax.set_xticks_delta(1)
     ax.set_yticks_delta(20)
     ax.set_xlabel('Contrast', r'\%')
     ax.set_ylabel('Amplitude', 'Hz')