minor updates in dec_const for vector decays

src/juobs_obs.jl

@@ -211,7 +211,7 @@ function dec_const(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64},
         aux = exp.((collect(0:T-1) .- T/2) .* [m for k = 1:T])
     end
 
-    R = corr_a0p .* aux ./ [sqrt(corr_pp[T-y0]) for k = 1:length(corr_a0p)]
+    R = corr_a0p .* aux ./ [((corr_pp[T-y0])^2)^(1/4) for k = 1:length(corr_a0p)]
     R_av = plat_av(R, plat, wpm)
 
     f = sqrt(2) * sqrt(R_av^2) / sqrt(m)
@@ -229,6 +229,7 @@ function dec_const(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64},
         fill_between(plat[1]:plat[2], v-e, v+e, color="green", alpha=0.75, label=L"$R$")
         errorbar(x, y, dy, fmt="x", color="black", label=lbl)
         legend()
+        ylim(v-10*e, v+10*e)
         ylabel(L"$R_\mathrm{av}$")
         xlabel(L"$x_0$")
     
@@ -236,6 +237,7 @@ function dec_const(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64},
             title(string(L"$\kappa_1 = $", kappa[1], L" $\kappa_2 = $", kappa[2]))
         end
         display(gcf())
+        close()
     end
     if !data
         return f