symmetric derivative of mpcac modified. ZV renormalization constant added

constants/juobs_const.jl

@@ -12,6 +12,8 @@ const ZM_tm_error = [42, 33, 42, 48, 55] .* 1e-4
 # 1808.09236
 const ZA_data = [0.75642, 0.76169, 0.76979, 0.78378, 0.79667]
 const ZA_err = [72, 93, 43, 47, 47] .*1e-5
+const ZV_data = [0.72259, 0.72845, 0.73886, 0.75574, 0.77074]
+const ZV_err = [74, 89, 46, 48, 49] .* 1e-5
 #2005.01352
 const ZS_over_ZP_data = [1.3497, 1.2914, 1.2317, 1.1709, 1.1343]
 const ZS_over_ZP_err = [83, 64, 48, 23, 25 ] .* 1e-4
@@ -29,6 +31,8 @@ const C3 = zeros(4, 4)
 const C4 = zeros(6, 6)
 const C5 = zeros(5, 5)
 const C6 = zeros(5, 5)
+const C7 = zeros(5, 5)
+
 for i = 1:6
     C4[i, i] = M_error[i] ^ 2
     if i<=5
@@ -36,6 +40,7 @@ for i = 1:6
         C2[i, i] = ZM_tm_error[i] ^ 2
         C5[i, i] = ZA_err[i] ^ 2
         C6[i, i] = ZS_over_ZP_err[i] ^ 2
+        C7[i, i] = ZV_err[i] ^ 2        
         if i <= 4
             C3[i, i] = t0_error[i] ^ 2
         end
@@ -49,6 +54,7 @@ const t0_ph = cobs(t0_ph_value, t0_ph_error .^ 2, "sqrt(8 t0) (fm)")
 const a_ = t0_ph ./ sqrt.(8 .* t0_)
 const M = cobs(M_values, C4, "charmed meson masses")
 const Za = cobs(ZA_data, C5, "ZA")
+const Zv = cobs(ZV_data, C7, "ZV")
 const ZS_over_ZP = cobs(ZS_over_ZP_data, C6, "ZS over ZP")
 
 zm(beta::Float64) = ZM[b_values .== beta][1]
@@ -56,4 +62,5 @@ zm_tm(beta::Float64) = ZM_tm[b_values .== beta][1]
 t0(beta::Float64) = t0_[b_values2 .== beta][1]
 a(beta::Float64) = a_[b_values2 .== beta][1]
 za(beta::Float64) = Za[b_values .== beta][1]
+zv(beta::Float64) = Zv[b_values .== beta][1]
 zs_over_zp(beta::Float64) = ZS_over_ZP[b_values .== beta][1]
\ No newline at end of file

src/juobs_obs.jl

@@ -33,6 +33,7 @@ function meff(corr::Vector{uwreal}, plat::Vector{Int64}; pl::Bool=true, data::Bo
         errorbar(x, y, dy, fmt="x", color="black")
         ylabel(L"$m_\mathrm{eff}$")
         xlabel(L"$x_0$")
+        ylim(v-10*e, v+10*e)
     
         if !isnothing(kappa)
             title(string(L"$\kappa_1 = $", kappa[1], L" $\kappa_2 = $", kappa[2]))
@@ -88,16 +89,22 @@ function mpcac(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64}; ca:
     kappa::Union{Vector{Float64}, Nothing}=nothing, mu::Union{Vector{Float64}, Nothing}=nothing, 
     wpm::Union{Dict{Int64,Vector{Float64}},Dict{String,Vector{Float64}}, Nothing}=nothing)     
     
-    corr_a0p = -a0p[2:end-1] 
-    corr_pp = pp[2:end-1]
+    corr_a0p = -a0p[1:end] 
+    corr_pp = pp[1:end]
 
     der_a0p = (corr_a0p[3:end] .- corr_a0p[1:end-2]) / 2 
     if ca != 0.0
-        der2_pp = corr_pp[1:end-2] + corr_pp[3:end] - 2 * corr_pp[2:end-1]
-        der_a0p = der_a0p + ca * der2_pp
+        der2_pp = (corr_pp[1:end-4] - 2*corr_pp[3:end-2] + corr_pp[5:end])/4
+        der_a0p = der_a0p[2:end-1] + ca * der2_pp
     end
 
-    aux = der_a0p ./ (2 .* corr_pp[2:end-1])
+    #der_a0p = (corr_a0p[3:end] .- corr_a0p[1:end-2]) / 2 
+    #if ca != 0.0
+    #    der2_pp = corr_pp[1:end-2] + corr_pp[3:end] - 2 * corr_pp[2:end-1]
+    #    der_a0p = der_a0p + ca * der2_pp
+    #end
+
+    aux = der_a0p ./ ( corr_pp[3:end-2])
     mass = plat_av(aux, plat, wpm)                                                                            
     if pl
         isnothing(wpm) ? uwerr(mass) : uwerr(mass, wpm)                       
@@ -112,6 +119,7 @@ function mpcac(a0p::Vector{uwreal}, pp::Vector{uwreal}, plat::Vector{Int64}; ca:
         errorbar(x, y, dy, fmt="x", color="black")
         ylabel(L"$m_\mathrm{PCAC}$")
         xlabel(L"$x_0$")
+        ylim(v-0.01, v+0.01)
     
         if !isnothing(kappa)
             title(string(L"$\kappa_1 = $", kappa[1], L" $\kappa_2 = $", kappa[2]))