begin
cd(joinpath(@__DIR__, ".."))
using Pkg
Pkg.activate(".")
Pkg.instantiate()
#
using Plots
using RecipesBase
#
using LinearAlgebra
using Parameters
#
using ThreeBodyDecay
end
 ❔❔
  Activating project at `/builds/polarimetry/Lc2pKpi/julia`
7.1 sconst ms
const ms = ThreeBodyMasses(0.938, 0.140, 0.493, m0=2.28)
14.2 msconst ms²
const ms² = ms^2
35.7 ms
@x_str (macro with 1 method)
begin
struct minusone end
import Base: ^
^(x::minusone, n::Int) = isodd(n) ? -1 : 1
macro x_str(s::String)
minusone()
end
end
1.6 msconst two_jΛ
1
const two_jΛ = 1
104 μs
wignerd_angle_doublearg (generic function with 1 method)
function wignerd_angle_doublearg(two_j, two_λ1, two_λ2, θ)
s, c = sincos(θ)
wd = wignerd_angle_doublearg(two_j, two_λ1, two_λ2, c)
return s > 0 ? wd : wd * x"-1"^(two_λ1 - two_λ2)
end
3.1 ms
cosαΔ (generic function with 1 method)
function cosαΔ(σs)
γp_Λc = (ms²[4] + ms²[1] - σs[1]) / (2 * ms[4] * ms[1])
γΔ_Λc = (ms²[4] + σs[3] - ms²[3]) / (2 * ms[4] * sqrt(σs[3]))
γp_Δ = (σs[3] + ms²[1] - ms²[2]) / (2 * ms[1] * sqrt(σs[3]))
(1 + γp_Λc + γΔ_Λc + γp_Δ)^2 / ((1 + γp_Λc) * (1 + γΔ_Λc) * (1 + γp_Δ)) - 1
end
1.3 ms
cosαΛ (generic function with 1 method)
function cosαΛ(σs)
γp_Λc = (ms²[4] + ms²[1] - σs[1]) / (2 * ms[4] * ms[1])
γΛ_Λc = (ms²[4] + σs[2] - ms²[2]) / (2 * ms[4] * sqrt(σs[2]))
γp_Λ = (σs[2] + ms²[1] - ms²[3]) / (2 * ms[1] * sqrt(σs[2]))
(1 + γp_Λc + γΛ_Λc + γp_Λ)^2 / ((1 + γp_Λc) * (1 + γΛ_Λc) * (1 + γp_Λ)) - 1
end
1.3 ms
plot(layout=grid(1, 2), size=(700, 250),
plot(ms, cosαΔ, colorbar=true),
plot(ms, cosαΛ, colorbar=true))
5.3 s
ξΛ (generic function with 1 method)
function ξΛ(σs)
θp = -acos(-cosθ31(σs, ms²)) # - (π-θ)
θΛ = acos(-cosζ12_for0(σs, ms²)) # π-θ
α = acos(cosαΛ(σs))
return θp + θΛ + α
end
752 μs
plot(ms, ξΛ, colorbar=true)
131 ms
ξΔ (generic function with 1 method)
function ξΔ(σs)
θp = acos(cosθ12(σs, ms²)) # just θ
θΔ = -acos(-cosζ31_for0(σs, ms²)) # -(π-θ)
α = acos(cosαΔ(σs))
return θp + θΔ - α
end
705 μs
plot(layout=grid(1, 2),
plot(ms, σs -> acos(cosζ13_for1(σs, ms²)) / ξΔ(σs) - 1,
colorbar=true, title="ξΔ = ζ13_for1"),
plot(ms, σs -> -acos(cosζ21_for1(σs, ms²)) / ξΛ(σs) - 1,
colorbar=true, title="ξΛ = - ζ21_for1"))
273 ms