Commit cc062c4a by Alessandro

### plot labels updated

parent 1ac1cf56
 ... ... @@ -476,7 +476,8 @@ uwerr(phi2_ph) x = [1 ./(8 .* getfield.(ensobs,:t0)) phi2] #x[1] = a^2 / (8t0), x[2] = 8t0 mpi^2 uwerr.(x) #modify CL+chiral fit model @. cl_chir_model(x, p) = (p[1] + p[2] * x[:, 1]) + (p[3]) * x[:, 2] #linear fits @. cl_chir_model(x, p) = (p[1] + p[2] * x[:, 1]) + (p[3]) * x[:, 2] #linear fits #@. cl_chir_model(x, p) = (p[1] + p[2] * x[:, 1]) * ((p[3]) * x[:, 2]) #non-linear fits obs_t0 = Vector{uwreal}(undef, length(obs)) #sqrt(8t0)obs @ CL & phi2_phys xarr = Float64.(range(0.0, stop=0.05, length=100)) ... ... @@ -492,7 +493,7 @@ for k = 1:length(obs) figure() for b in unique(phi2) #select point with same beta nn = findall(x-> x == b, phi2) lgnd = string(L"$m_\pi =$", Int32(round(sqrt(value(b))*hc/t0_ph_value[1], sigdigits=3)), " MeV") lgnd = string(L"$m_\pi =$", Int32(round(sqrt(value(b))*hc/t0_ph_value[1], sigdigits=3)), " MeV") errorbar(value.(x[nn,1]), value.(obs[k][nn]), err.(obs[k][nn]), ms=5, fmt="s", mfc="none", mew=0.8, elinewidth=0.8, capsize=2, label=lgnd) end fill_between(xarr, value.(y) .+ err.(y),value.(y) .- err.(y) , color="tomato", alpha=0.4) ... ...
 ... ... @@ -27,7 +27,7 @@ ylbl_db = Dict( "lh" => [L"$m_D \sqrt{8t_0}$", L"$m_{D^*} \sqrt{8t_0}$", L"$f_D \sqrt{8t_0}$", L"$f_{D^*} \sqrt{8t_0}$"], "ss" => [L"$m_{\eta'} \sqrt{8t_0}$", L"$m_{\phi} \sqrt{8t_0}$", L"$f_{η'} \sqrt{8t_0}$", L"$f_{\phi} \sqrt{8t_0}$"], "sh" => [L"$m_{D_s} \sqrt{8t_0}$", L"$m_{D_s^*} \sqrt{8t_0}$", L"$f_{D_s} \sqrt{8t_0}$", L"$f_{D_s^*} \sqrt{8t_0}$"], "hh" => [L"$m_etac \sqrt{8t_0}$", L"$m_{J/ψ} \sqrt{8t_0}$", L"$f_{η_c} \sqrt{8t_0}$", L"$f_{J/ψ} \sqrt{8t_0}$"] "hh" => [L"$m_{η_c} \sqrt{8t_0}$", L"$m_{J/ψ} \sqrt{8t_0}$", L"$f_{η_c} \sqrt{8t_0}$", L"$f_{J/ψ} \sqrt{8t_0}$"] ) #PDG ... ...
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment