Model serialization#

This notebooks illustrates the use of the ampform_dpd.io.serialization module, which can be used to build amplitude models from the amplitude-serialization initiative.

Warning

The ampform_dpd.io.serialization module is a preview feature.

Import model#

load_workspace() loads the serialized definition and formulates its distributions, decays, and dynamics functions in one step. It also preserves kinematics, parameter points, and checksums for downstream backends and validation.

def to_latex(name: str) -> str:
    latex = {
        "Lc": R"\Lambda_c^+",
        "pi": R"\pi^+",
        "K": "K^-",
        "p": "p",
    }.get(name)
    if latex is not None:
        return latex
    mass_str = name[1:].strip("(").strip(")")
    subsystem_letter = name[0]
    subsystem = {"D": "D", "K": "K", "L": R"\Lambda"}.get(subsystem_letter)
    if subsystem is None:
        return name
    return f"{subsystem}({mass_str})"

Construct ThreeBodyDecay#

Hide code cell source

((INTENSITY_NAME, DECAY),) = WORKSPACE.decays.items()
Math(aslatex(DECAY, with_jp=True))
\[\begin{split}\displaystyle \begin{array}{c} \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(D(1232)\left[J=\frac{3}{2}\right] \to p\left[J=\frac{1}{2}\right] \pi^+\left[J=0\right]\right) K^-\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(D(1600)\left[J=\frac{3}{2}\right] \to p\left[J=\frac{1}{2}\right] \pi^+\left[J=0\right]\right) K^-\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(D(1700)\left[J=\frac{3}{2}\right] \to p\left[J=\frac{1}{2}\right] \pi^+\left[J=0\right]\right) K^-\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(K(1430)\left[J=0\right] \to \pi^+\left[J=0\right] K^-\left[J=0\right]\right) p\left[J=\frac{1}{2}\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(K(700)\left[J=0\right] \to \pi^+\left[J=0\right] K^-\left[J=0\right]\right) p\left[J=\frac{1}{2}\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(K(892)\left[J=1\right] \to \pi^+\left[J=0\right] K^-\left[J=0\right]\right) p\left[J=\frac{1}{2}\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(1405)\left[J=\frac{1}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(1520)\left[J=\frac{3}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(1600)\left[J=\frac{1}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(1670)\left[J=\frac{1}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(1690)\left[J=\frac{3}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \Lambda_c^+\left[J=\frac{1}{2}\right] \to \left(\Lambda(2000)\left[J=\frac{1}{2}\right] \to K^-\left[J=0\right] p\left[J=\frac{1}{2}\right]\right) \pi^+\left[J=0\right] \\ \end{array}\end{split}\]
WORKSPACE.kinematics[INTENSITY_NAME]
mappingproxy({'initial_state': mappingproxy({'index': 0,
                            'name': 'Lc',
                            'spin': '1/2',
                            'mass': 2.28646}),
              'final_state': (mappingproxy({'index': 1,
                             'name': 'p',
                             'spin': '1/2',
                             'mass': 0.938272046}),
               mappingproxy({'index': 2,
                             'name': 'pi',
                             'spin': '0',
                             'mass': 0.13957018}),
               mappingproxy({'index': 3,
                             'name': 'K',
                             'spin': '0',
                             'mass': 0.493677}))})

Dynamics#

CHAIN_DEFS = get_decay_chains(MODEL_DEFINITION)

Vertices#

Blatt–Weisskopf form factor#

See also

Form factor in the AmpForm documentation and TR-029

Hide code cell source

s, m1, m2, L, d = sp.symbols("s m1 m2 L R", nonnegative=True)
expr = FormFactor(s, m1, m2, L, d)
Math(aslatex(unfold_definitions(expr)))
\[\begin{split}\displaystyle \begin{aligned} \mathcal{F}_{L}\left(s, m_{1}, m_{2}\right) \;&=\; \sqrt{B_{L}^2\left(R^{2} q^2\left(s\right)\right)} \\ B_{L}^2\left(R^{2} q^2\left(s\right)\right) \;&=\; \frac{\left|{h_{L}^{(1)}\left(1\right)}\right|^{2}}{R^{2} \left|{h_{L}^{(1)}\left(R \sqrt{q^2\left(s\right)}\right)}\right|^{2} q^2\left(s\right)} \\ q^2\left(s\right) \;&=\; \frac{\left(s - \left(m_{1} - m_{2}\right)^{2}\right) \left(s - \left(m_{1} + m_{2}\right)^{2}\right)}{4 s} \\ h_{L}^{(1)}\left(z\right) \;&=\; \frac{\left(- i\right)^{L + 1} e^{i z} \sum_{k=0}^{L} \frac{\left(\frac{i}{2 z}\right)^{k} \left(k + L\right)!}{k! \left(- k + L\right)!}}{z} \\ \end{aligned}\end{split}\]
ff_L1520_prod = formulate_form_factor(
    vertex=CHAIN_DEFS[2]["vertices"][0],
    model=MODEL_DEFINITION,
)
ff_L1520_dec = formulate_form_factor(
    vertex=CHAIN_DEFS[2]["vertices"][1],
    model=MODEL_DEFINITION,
)
Math(aslatex([ff_L1520_prod, ff_L1520_dec]))
\[\begin{split}\displaystyle \begin{array}{c} \begin{array}{rcl} \frac{\sqrt{2} \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right)}{2} \\ R_{Lc} &=& 5.0 \\ \end{array} \\ \begin{array}{rcl} \frac{\sqrt{13} \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right)}{13} \\ R_\mathrm{res} &=& 1.5 \\ \end{array} \\ \end{array}\end{split}\]

Propagators#

Breit–Wigner#

Hide code cell source

s, m0, Γ0, m1, m2, L, d = sp.symbols("s m0 Gamma0 m1 m2 L R", nonnegative=True)
expr = BreitWigner(s, m0, Γ0, m1, m2, L, d)
Math(aslatex(unfold_definitions(expr)))
\[\begin{split}\displaystyle \begin{aligned} \mathcal{R}^\mathrm{BW}_{L}\left(s; m_{0}, \Gamma_{0}\right) \;&=\; \mathcal{R}^\mathrm{BW}\left(s; m_{0}, \Gamma_{0}\left(s\right)\right) \\ \mathcal{R}^\mathrm{BW}\left(s; m_{0}, \Gamma_{0}\left(s\right)\right) \;&=\; \frac{1}{m_{0}^{2} - i m_{0} \Gamma_{0}\left(s\right) - s} \\ \Gamma_{0}\left(s\right) \;&=\; \frac{\Gamma_{0} \mathcal{F}_{L}\left(s, m_{1}, m_{2}\right)^{2} \rho\left(s\right)}{\mathcal{F}_{L}\left(m_{0}^{2}, m_{1}, m_{2}\right)^{2} \rho_{0}\left(m_{0}^{2}\right)} \\ \mathcal{F}_{L}\left(s, m_{1}, m_{2}\right) \;&=\; \sqrt{B_{L}^2\left(R^{2} q^2\left(s\right)\right)} \\ \rho\left(s\right) \;&=\; \frac{\sqrt{\left(s - \left(m_{1} - m_{2}\right)^{2}\right) \left(s - \left(m_{1} + m_{2}\right)^{2}\right)}}{s} \\ B_{L}^2\left(R^{2} q^2\left(s\right)\right) \;&=\; \frac{\left|{h_{L}^{(1)}\left(1\right)}\right|^{2}}{R^{2} \left|{h_{L}^{(1)}\left(R \sqrt{q^2\left(s\right)}\right)}\right|^{2} q^2\left(s\right)} \\ q^2\left(s\right) \;&=\; \frac{\left(s - \left(m_{1} - m_{2}\right)^{2}\right) \left(s - \left(m_{1} + m_{2}\right)^{2}\right)}{4 s} \\ h_{L}^{(1)}\left(z\right) \;&=\; \frac{\left(- i\right)^{L + 1} e^{i z} \sum_{k=0}^{L} \frac{\left(\frac{i}{2 z}\right)^{k} \left(k + L\right)!}{k! \left(- k + L\right)!}}{z} \\ \end{aligned}\end{split}\]
K892_BW = formulate_breit_wigner(
    propagator=CHAIN_DEFS[20]["propagators"][0],
    resonance=to_latex(CHAIN_DEFS[20]["name"]),
    model=MODEL_DEFINITION,
)
Math(aslatex(K892_BW))
\[\begin{split}\displaystyle \begin{array}{rcl} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) &=& \mathcal{R}^\mathrm{BW}\left(\sigma_{1}; m_{K(892)}, \Gamma_{892}\left(\sigma_{1}\right)\right) \\ m_{K(892)} &=& 0.8955 \\ \Gamma_{K(892)} &=& 0.047299999999999995 \\ m_{2} &=& 0.13957018 \\ m_{3} &=& 0.493677 \\ R_\mathrm{res} &=& 1.5 \\ \end{array}\end{split}\]

Multi-channel Breit–Wigner#

Each channel in the serialized model carries a gsq value, which is a coupling squared \(g_i^2\), not an energy-dependent width. It enters the lineshape through the channel term \(\Gamma^\text{ch}\) as follows.

Hide code cell source

s, m0, m1, m2, L, d = sp.symbols("s m0 m1 m2 L R", nonnegative=True)
Γ = sp.Symbol("Gamma", nonnegative=True)
channels = tuple(
    ChannelArguments(
        s,
        m0,
        coupling_squared=sp.Symbol(f"g_{{{i}}}^2", nonnegative=True),
        m1=sp.Symbol(f"m_{{a,{i}}}", nonnegative=True),
        m2=sp.Symbol(f"m_{{b,{i}}}", nonnegative=True),
        angular_momentum=sp.Symbol(f"L_{{{i}}}", integer=True, nonnegative=True),
        meson_radius=d,
    )
    for i in [1, 2]
)
expr = MultichannelBreitWigner(s, m0, channels)
Math(aslatex(unfold_definitions(expr)))
\[\begin{split}\displaystyle \begin{aligned} \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(s; g_{1}^2, g_{2}^2\right) \;&=\; \mathcal{R}^\mathrm{BW}_{L=0}\left(s; m_{0}, \Gamma^\text{ch}\left(s; m_{0}, g_{1}^2\right) + \Gamma^\text{ch}\left(s; m_{0}, g_{2}^2\right)\right) \\ \mathcal{R}^\mathrm{BW}_{L=0}\left(s; m_{0}, \Gamma^\text{ch}\left(s; m_{0}, g_{1}^2\right) + \Gamma^\text{ch}\left(s; m_{0}, g_{2}^2\right)\right) \;&=\; \frac{1}{m_{0}^{2} - i m_{0} \left(\Gamma^\text{ch}\left(s; m_{0}, g_{1}^2\right) + \Gamma^\text{ch}\left(s; m_{0}, g_{2}^2\right)\right) - s} \\ \Gamma^\text{ch}\left(s; m_{0}, g_{1}^2\right) \;&=\; \frac{g_{1}^2 \mathcal{F}_{L_{1}}\left(s, m_{a,1}, m_{b,1}\right)^{2} \rho\left(s\right)}{m_{0}} \\ \mathcal{F}_{L_{1}}\left(s, m_{a,1}, m_{b,1}\right) \;&=\; \sqrt{B_{L_{1}}^2\left(R^{2} q^2\left(s\right)\right)} \\ \rho\left(s\right) \;&=\; \frac{\sqrt{\left(s - \left(m_{a,1} - m_{b,1}\right)^{2}\right) \left(s - \left(m_{a,1} + m_{b,1}\right)^{2}\right)}}{s} \\ B_{L_{1}}^2\left(R^{2} q^2\left(s\right)\right) \;&=\; \frac{\left|{h_{L_{1}}^{(1)}\left(1\right)}\right|^{2}}{R^{2} \left|{h_{L_{1}}^{(1)}\left(R \sqrt{q^2\left(s\right)}\right)}\right|^{2} q^2\left(s\right)} \\ q^2\left(s\right) \;&=\; \frac{\left(s - \left(m_{a,1} - m_{b,1}\right)^{2}\right) \left(s - \left(m_{a,1} + m_{b,1}\right)^{2}\right)}{4 s} \\ h_{L_{1}}^{(1)}\left(z\right) \;&=\; \frac{\left(- i\right)^{L_{1} + 1} e^{i z} \sum_{k=0}^{L_{1}} \frac{\left(\frac{i}{2 z}\right)^{k} \left(k + L_{1}\right)!}{k! \left(- k + L_{1}\right)!}}{z} \\ \end{aligned}\end{split}\]
L1405_Flatte = formulate_multichannel_breit_wigner(
    propagator=CHAIN_DEFS[0]["propagators"][0],
    resonance=to_latex(CHAIN_DEFS[0]["name"]),
    model=MODEL_DEFINITION,
)
Math(aslatex(L1405_Flatte))
\[\begin{split}\displaystyle \begin{array}{rcl} \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) &=& \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma^\text{ch}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}\right) + \Gamma^\text{ch}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right)\right) \\ m_{\Lambda(1405)} &=& 1.4051 \\ \Gamma_{\Lambda(1405)} &=& 0.328725260215546 \\ m_{3} &=& 0.938272046 \\ m_{1} &=& 0.493677 \\ R_{\Lambda(1405)} &=& 0 \\ m_{a,2} &=& 1.18937 \\ m_{b,2} &=& 0.13957018 \\ \Gamma_{\Lambda(1405)}^\text{ch. 2} &=& 0.328725260215546 \\ \end{array}\end{split}\]

Breit–Wigner with exponential#

The Bugg lineshapes are serialized as generic_function expression strings. The built-in dynamics builder parses these expressions and substitutes \(i\) and \(\sigma\) with the imaginary unit and the Mandelstam variable of the propagator node, respectively.

CHAIN_DEFS[18]
mappingproxy({'propagators': (mappingproxy({'spin': '0',
                             'node': (2, 3),
                             'parametrization': 'K700_BuggBW'}),),
              'weight': '0.068908 + 2.521444i',
              'vertices': (mappingproxy({'type': 'helicity',
                             'helicities': ('0', '1/2'),
                             'node': ((2, 3), 1),
                             'formfactor': ''}),
               mappingproxy({'type': 'parity',
                             'helicities': ('0', '0'),
                             'parity_factor': '+',
                             'node': (2, 3),
                             'formfactor': ''})),
              'topology': ((2, 3), 1),
              'name': 'K700'})
get_function_definition("K700_BuggBW", MODEL_DEFINITION)
mappingproxy({'name': 'K700_BuggBW',
              'type': 'generic_function',
              'expression': '1/(0.824^2 - σ - i * 0.824 * (σ - 0.23397706275638377) / (0.824^2 - 0.23397706275638377) * 0.478 * exp(-0.941060 * σ))'})

General propagator dynamics builder#

Hide code cell source

exprs = [
    formulate_dynamics(CHAIN_DEFS[0], MODEL_DEFINITION, to_latex),
    formulate_dynamics(CHAIN_DEFS[18], MODEL_DEFINITION, to_latex),
    formulate_dynamics(CHAIN_DEFS[20], MODEL_DEFINITION, to_latex),
]
Math(aslatex(exprs))
\[\begin{split}\displaystyle \begin{array}{c} \begin{array}{rcl} \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) &=& \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma^\text{ch}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}\right) + \Gamma^\text{ch}\left(\sigma_{2}; m_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right)\right) \\ m_{\Lambda(1405)} &=& 1.4051 \\ \Gamma_{\Lambda(1405)} &=& 0.328725260215546 \\ m_{3} &=& 0.938272046 \\ m_{1} &=& 0.493677 \\ R_{\Lambda(1405)} &=& 0 \\ m_{a,2} &=& 1.18937 \\ m_{b,2} &=& 0.13957018 \\ \Gamma_{\Lambda(1405)}^\text{ch. 2} &=& 0.328725260215546 \\ \end{array} \\ \begin{array}{rcl} \frac{1}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976} \\ \end{array} \\ \begin{array}{rcl} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) &=& \mathcal{R}^\mathrm{BW}\left(\sigma_{1}; m_{K(892)}, \Gamma_{892}\left(\sigma_{1}\right)\right) \\ m_{K(892)} &=& 0.8955 \\ \Gamma_{K(892)} &=& 0.047299999999999995 \\ m_{2} &=& 0.13957018 \\ m_{3} &=& 0.493677 \\ R_\mathrm{res} &=& 1.5 \\ \end{array} \\ \end{array}\end{split}\]

Validation#

The propagators are compared with the reference values that the model file provides under misc.amplitude_model_checksums. The marks 🟢, 🟡, and 🔴 indicate an accuracy of \(<10^{-10}\), \(<10^{-2}\), and \(\geq10^{-2}\), respectively, on the difference between the reference and the computed value.

See also

The same table for the reference implementation in Julia

Hide code cell source

def to_complex(value: complex | str) -> complex:
    if isinstance(value, str):
        return complex(value.replace(" ", "").replace("i", "j"))
    return complex(value)


def label_diff(difference: complex) -> str:
    absolute_difference = abs(difference)
    if absolute_difference < 1e-10:
        return "🟢"
    if absolute_difference < 1e-2:
        return "🟡"
    return "🔴"


CHAINS_BY_PROPAGATOR = {
    propagator["parametrization"]: chain
    for chain in CHAIN_DEFS
    for propagator in chain["propagators"]
}
PARAMETER_POINTS = {
    point["name"]: {par["name"]: par["value"] for par in point["parameters"]}
    for point in WORKSPACE.reference_points
}

propagator_checks = []
for checksum in WORKSPACE.checksums:
    parametrization = checksum["distribution"]
    chain = CHAINS_BY_PROPAGATOR.get(parametrization)
    if chain is None:  # the checksum is for the full intensity, not a propagator
        continue
    dynamics = WORKSPACE.functions[parametrization]
    expr = dynamics.expression.doit().xreplace(dynamics.parameters)
    (variable,) = expr.free_symbols
    (value,) = PARAMETER_POINTS[checksum["point"]].values()
    computed = complex(expr.subs(variable, value).evalf())
    expected = to_complex(checksum["value"])
    difference = expected - computed
    propagator_checks.append({
        "Propagator": parametrization,
        "Point": checksum["point"],
        "Computed": computed,
        "Expected": expected,
        "Difference": abs(difference),
        "Status": label_diff(difference),
    })
propagator_validation = pd.DataFrame(propagator_checks)
propagator_validation[["Propagator", "Point", "Difference", "Status"]]
Propagator Point Difference Status
0 L1405_Flatte validation_point_m31sq 0.000000e+00 🟢
1 L1690_BW validation_point_m31sq 1.404333e-15 🟢
2 D1232_BW validation_point_m12sq 0.000000e+00 🟢
3 L1520_BW validation_point_m31sq 1.570092e-16 🟢
4 L1600_BW validation_point_m31sq 1.110223e-16 🟢
5 L2000_BW validation_point_m31sq 0.000000e+00 🟢
6 D1600_BW validation_point_m12sq 1.053250e-15 🟢
7 D1700_BW validation_point_m12sq 1.241267e-16 🟢
8 K892_BW validation_point_m23sq 3.554448e-16 🟢
9 L1670_BW validation_point_m31sq 4.440892e-16 🟢

Hide code cell source

np.testing.assert_allclose(
    propagator_validation["Computed"].to_numpy(),
    propagator_validation["Expected"].to_numpy(),
    atol=1e-10,
    rtol=0,
)

Construct AmplitudeModel#

Unpolarized intensity#

λ0, λ1, λ2, λ3 = sp.symbols("lambda(:4)", rational=True)
amplitude_expr, _ = formulate_aligned_amplitude(MODEL_DEFINITION, λ0, λ1, λ2, λ3)
amplitude_expr.cleanup()
\[\displaystyle \sum_{\lambda_0^{\prime}=-1/2}^{1/2} \sum_{\lambda_1^{\prime}=-1/2}^{1/2}{A^{1}_{\lambda_0^{\prime}, \lambda_1^{\prime}, 0, 0} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{1(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{1(1)}\right) + A^{2}_{\lambda_0^{\prime}, \lambda_1^{\prime}, 0, 0} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{2(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{2(1)}\right) + A^{3}_{\lambda_0^{\prime}, \lambda_1^{\prime}, 0, 0} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{3(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{3(1)}\right)}\]

Amplitude for the decay chain#

Helicity recouplings#

Hide code cell source

λa = sp.Symbol(R"\lambda_a", rational=True)
λb = sp.Symbol(R"\lambda_b", rational=True)
λa0 = sp.Symbol(R"\lambda_a^0", rational=True)
λb0 = sp.Symbol(R"\lambda_b^0", rational=True)
f = sp.Symbol("f", integer=True)
l = sp.Symbol("l", integer=True, nonnegative=True)
s = sp.Symbol("s", nonnegative=True, rational=True)
ja = sp.Symbol("j_a", nonnegative=True, rational=True)
jb = sp.Symbol("j_b", nonnegative=True, rational=True)
j = sp.Symbol("j", nonnegative=True, rational=True)
exprs = [
    HelicityRecoupling(λa, λb, λa0, λb0),
    ParityRecoupling(λa, λb, λa0, λb0, f),
    LSRecoupling(λa, λb, l, s, ja, jb, j),
]
Math(aslatex({e: e.doit(deep=False) for e in exprs}))
\[\begin{split}\displaystyle \begin{aligned} \mathcal{H}^\text{helicity}\left(\lambda_{a},\lambda_{b}\middle|\lambda^{0}_{a},\lambda^{0}_{b}\right) \;&=\; \delta_{\lambda_{a} \lambda^{0}_{a}} \delta_{\lambda_{b} \lambda^{0}_{b}} \\ \mathcal{H}^\text{parity}\left(\lambda_{a},\lambda_{b}\middle|\lambda^{0}_{a},\lambda^{0}_{b},f\right) \;&=\; f \delta_{\lambda_{a}, - \lambda^{0}_{a}} \delta_{\lambda_{b}, - \lambda^{0}_{b}} + \delta_{\lambda_{a} \lambda^{0}_{a}} \delta_{\lambda_{b} \lambda^{0}_{b}} \\ \mathcal{H}^\text{parity}\left(\lambda_{a},\lambda_{b}\middle|l,s,j_{a},j_{b},j\right) \;&=\; \frac{\sqrt{2 l + 1} C^{s,\lambda_{a} - \lambda_{b}}_{j_{a},\lambda_{a},j_{b},- \lambda_{b}} C^{j,\lambda_{a} - \lambda_{b}}_{l,0,s,\lambda_{a} - \lambda_{b}}}{\sqrt{2 j + 1}} \\ \end{aligned}\end{split}\]

Recoupling deserialization#

Hide code cell source

recouplings = [
    formulate_recoupling(MODEL_DEFINITION, chain_idx=0, vertex_idx=i) for i in range(2)
]
Math(aslatex({e: e.doit(deep=False) for e in recouplings}))
\[\begin{split}\displaystyle \begin{aligned} \mathcal{H}^\text{helicity}\left(\lambda_{R},\lambda_{2}\middle|\frac{1}{2},0\right) \;&=\; \delta_{0 \lambda_{2}} \delta_{\frac{1}{2} \lambda_{R}} \\ \mathcal{H}^\text{parity}\left(\lambda_{3},\lambda_{1}\middle|0,\frac{1}{2},1\right) \;&=\; \delta_{- \frac{1}{2} \lambda_{1}} \delta_{0 \lambda_{3}} + \delta_{0 \lambda_{3}} \delta_{\frac{1}{2} \lambda_{1}} \\ \end{aligned}\end{split}\]

Chain amplitudes#

definitions = formulate_chain_amplitude(λ0, λ1, λ2, λ3, MODEL_DEFINITION, chain_idx=0)
Math(aslatex(definitions))
\[\begin{split}\displaystyle \begin{aligned} A^{2}_{\lambda_{0}, \lambda_{1}, \lambda_{2}, \lambda_{3}} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\left(-1\right)^{- \lambda_{2}} \left(-1\right)^{\frac{1}{2} - \lambda_{1}} \sqrt{2} c^{L1405[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\lambda_{0}, \lambda_{R} - \lambda_{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\lambda_{2}\middle|\frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{L1405}, \Gamma_{L1405}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(\lambda_{3},\lambda_{1}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \lambda_{1} + \lambda_{3}}\left(\theta_{31}\right)} \\ c^{L1405[1/2]}_{\frac{1}{2}, 0, 0} \;&=\; 7.38649400481717+1.971018433257411i \\ m_{L1405} \;&=\; 1.4051 \\ \Gamma_{L1405} \;&=\; 0.328725260215546 \\ m_{3} \;&=\; 0.938272046 \\ m_{1} \;&=\; 0.493677 \\ R_{L1405} \;&=\; 0 \\ m_{a,2} \;&=\; 1.18937 \\ m_{b,2} \;&=\; 0.13957018 \\ \Gamma_{L1405}^\text{ch. 2} \;&=\; 0.328725260215546 \\ \theta_{31} \;&=\; \operatorname{acos}{\left(\frac{2 \sigma_{2} \left(- m_{2}^{2} - m_{3}^{2} + \sigma_{1}\right) - \left(m_{0}^{2} - m_{2}^{2} - \sigma_{2}\right) \left(- m_{1}^{2} + m_{3}^{2} + \sigma_{2}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{2}^{2}, \sigma_{2}\right)} \sqrt{\lambda\left(\sigma_{2}, m_{3}^{2}, m_{1}^{2}\right)}} \right)} \\ \end{aligned}\end{split}\]

Full amplitude model#

MODEL = WORKSPACE.distributions[INTENSITY_NAME]
MODEL.intensity
\[\displaystyle \sum_{\lambda_{0}=-1/2}^{1/2} \sum_{\lambda_{1}=-1/2}^{1/2} \sum_{\lambda_{2}=0} \sum_{\lambda_{3}=0}{\left|{\sum_{\lambda_0^{\prime}=-1/2}^{1/2} \sum_{\lambda_1^{\prime}=-1/2}^{1/2} \sum_{\lambda_2^{\prime}=0} \sum_{\lambda_3^{\prime}=0}{A^{1}_{\lambda_0^{\prime}, \lambda_1^{\prime}, \lambda_2^{\prime}, \lambda_3^{\prime}} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{1(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{1(1)}\right) + A^{2}_{\lambda_0^{\prime}, \lambda_1^{\prime}, \lambda_2^{\prime}, \lambda_3^{\prime}} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{2(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{2(1)}\right) + A^{3}_{\lambda_0^{\prime}, \lambda_1^{\prime}, \lambda_2^{\prime}, \lambda_3^{\prime}} d^{\frac{1}{2}}_{\lambda_1^{\prime},\lambda_{1}}\left(\zeta^1_{3(1)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{3(1)}\right)}}\right|^{2}}\]

Hide code cell source

if "EXECUTE_NB" in os.environ:
    selected_amplitudes = MODEL.amplitudes
else:
    selected_amplitudes = {
        k: v for i, (k, v) in enumerate(MODEL.amplitudes.items()) if i < 2
    }
Math(aslatex(selected_amplitudes, terms_per_line=1))
\[\begin{split}\displaystyle \begin{aligned} A^{2}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1405)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1405)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1670)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1670)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(2000)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(2000)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- c^{\Lambda(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- c^{\Lambda(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1520)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1520)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1690)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1690)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ A^{3}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ A^{1}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=0}{- \frac{c^{K(1430)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(1430)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(700)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(700)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[-1]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|-1,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[1]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|1,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1405)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1405)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1670)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1670)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(2000)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(2000)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{c^{\Lambda(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{c^{\Lambda(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1520)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1520)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1690)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1690)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ A^{3}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ A^{1}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=0}{\frac{c^{K(1430)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(1430)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(700)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(700)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[-1]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|-1,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[0]}_{- \frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[0]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[1]}_{\frac{1}{2}, 0, 0} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|1,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1405)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1405)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1670)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(1670)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(2000)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- \sqrt{2} c^{\Lambda(2000)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- c^{\Lambda(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{- c^{\Lambda(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1520)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1520)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1690)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{26} c^{\Lambda(1690)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,- \frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ A^{3}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(- \frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ A^{1}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=0}{- \frac{c^{K(1430)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(1430)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(700)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{c^{K(700)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[-1]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|-1,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{6} c^{K(892)[1]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},- \frac{1}{2}\middle|1,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1405)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1405)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{R}^\mathrm{BW}_\mathrm{multi}\left(\sigma_{2}; \Gamma_{\Lambda(1405)}, \Gamma_{\Lambda(1405)}^\text{ch. 2}\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1670)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(1670)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(1670)}, \Gamma_{\Lambda(1670)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(2000)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\sqrt{2} c^{\Lambda(2000)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=0}\left(\sigma_{2}; m_{\Lambda(2000)}, \Gamma_{\Lambda(2000)}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{c^{\Lambda(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{c^{\Lambda(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{2}; m_{\Lambda(1600)}, \Gamma_{\Lambda(1600)}\right) \mathcal{F}_{1}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1520)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1520)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1520)}, \Gamma_{\Lambda(1520)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1690)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{\Lambda(1690)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{2}; m_{\Lambda(1690)}, \Gamma_{\Lambda(1690)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{2}}, m_{2}\right) \mathcal{F}_{2}\left(\sigma_{2}, m_{3}, m_{1}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(0,\frac{1}{2}\middle|0,\frac{1}{2},-1\right) d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{13}} \\ A^{3}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1232)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1232)}, \Gamma_{D(1232)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{c^{D(1600)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{3}; m_{D(1600)}, \Gamma_{D(1600)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{1}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[-1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|- \frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{26} c^{D(1700)[1/2]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{L=2}\left(\sigma_{3}; m_{D(1700)}, \Gamma_{D(1700)}\right) \mathcal{F}_{1}\left(m_{0}^{2}, \sqrt{\sigma_{3}}, m_{3}\right) \mathcal{F}_{2}\left(\sigma_{3}, m_{1}, m_{2}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},0\middle|\frac{1}{2},0\right) \mathcal{H}^\text{parity}\left(\frac{1}{2},0\middle|\frac{1}{2},0,-1\right) d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{13}} \\ A^{1}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=0}{\frac{c^{K(1430)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(1430)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.15769795991457068 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.020981 \sigma_{1}} + 1.890625}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(700)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{c^{K(700)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{- \sigma_{1} - 0.88510773180649964 i \left(\sigma_{1} - 0.23397706275638377\right) e^{- 0.94106 \sigma_{1}} + 0.678976}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[-1]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|-1,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[0]}_{- \frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,- \frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[0]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|0,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{6} c^{K(892)[1]}_{\frac{1}{2}, 0, 0} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{L=1}\left(\sigma_{1}; m_{K(892)}, \Gamma_{K(892)}\right) \mathcal{F}_{1}\left(\sigma_{1}, m_{2}, m_{3}\right) \mathcal{H}^\text{helicity}\left(\lambda_{R},\frac{1}{2}\middle|1,\frac{1}{2}\right) \mathcal{H}^\text{parity}\left(0,0\middle|0,0,1\right) d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \end{aligned}\end{split}\]

Hide code cell source

Math(aslatex(MODEL.variables))
\[\begin{split}\displaystyle \begin{aligned} \theta_{31} \;&=\; \operatorname{acos}{\left(\frac{2 \sigma_{2} \left(- m_{2}^{2} - m_{3}^{2} + \sigma_{1}\right) - \left(m_{0}^{2} - m_{2}^{2} - \sigma_{2}\right) \left(- m_{1}^{2} + m_{3}^{2} + \sigma_{2}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{2}^{2}, \sigma_{2}\right)} \sqrt{\lambda\left(\sigma_{2}, m_{3}^{2}, m_{1}^{2}\right)}} \right)} \\ \theta_{12} \;&=\; \operatorname{acos}{\left(\frac{2 \sigma_{3} \left(- m_{1}^{2} - m_{3}^{2} + \sigma_{2}\right) - \left(m_{0}^{2} - m_{3}^{2} - \sigma_{3}\right) \left(m_{1}^{2} - m_{2}^{2} + \sigma_{3}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{3}^{2}, \sigma_{3}\right)} \sqrt{\lambda\left(\sigma_{3}, m_{1}^{2}, m_{2}^{2}\right)}} \right)} \\ \theta_{23} \;&=\; \operatorname{acos}{\left(\frac{2 \sigma_{1} \left(- m_{1}^{2} - m_{2}^{2} + \sigma_{3}\right) - \left(m_{0}^{2} - m_{1}^{2} - \sigma_{1}\right) \left(m_{2}^{2} - m_{3}^{2} + \sigma_{1}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{1}^{2}, \sigma_{1}\right)} \sqrt{\lambda\left(\sigma_{1}, m_{2}^{2}, m_{3}^{2}\right)}} \right)} \\ \zeta^0_{1(1)} \;&=\; 0 \\ \zeta^1_{1(1)} \;&=\; 0 \\ \zeta^0_{2(1)} \;&=\; - \operatorname{acos}{\left(\frac{- 2 m_{0}^{2} \left(- m_{1}^{2} - m_{2}^{2} + \sigma_{3}\right) + \left(m_{0}^{2} + m_{1}^{2} - \sigma_{1}\right) \left(m_{0}^{2} + m_{2}^{2} - \sigma_{2}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{2}^{2}, \sigma_{2}\right)} \sqrt{\lambda\left(m_{0}^{2}, \sigma_{1}, m_{1}^{2}\right)}} \right)} \\ \zeta^1_{2(1)} \;&=\; \operatorname{acos}{\left(\frac{2 m_{1}^{2} \left(- m_{0}^{2} - m_{3}^{2} + \sigma_{3}\right) + \left(m_{0}^{2} + m_{1}^{2} - \sigma_{1}\right) \left(- m_{1}^{2} - m_{3}^{2} + \sigma_{2}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{1}^{2}, \sigma_{1}\right)} \sqrt{\lambda\left(\sigma_{2}, m_{1}^{2}, m_{3}^{2}\right)}} \right)} \\ \zeta^0_{3(1)} \;&=\; \operatorname{acos}{\left(\frac{- 2 m_{0}^{2} \left(- m_{1}^{2} - m_{3}^{2} + \sigma_{2}\right) + \left(m_{0}^{2} + m_{1}^{2} - \sigma_{1}\right) \left(m_{0}^{2} + m_{3}^{2} - \sigma_{3}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{1}^{2}, \sigma_{1}\right)} \sqrt{\lambda\left(m_{0}^{2}, \sigma_{3}, m_{3}^{2}\right)}} \right)} \\ \zeta^1_{3(1)} \;&=\; - \operatorname{acos}{\left(\frac{2 m_{1}^{2} \left(- m_{0}^{2} - m_{2}^{2} + \sigma_{2}\right) + \left(m_{0}^{2} + m_{1}^{2} - \sigma_{1}\right) \left(- m_{1}^{2} - m_{2}^{2} + \sigma_{3}\right)}{\sqrt{\lambda\left(m_{0}^{2}, m_{1}^{2}, \sigma_{1}\right)} \sqrt{\lambda\left(\sigma_{3}, m_{1}^{2}, m_{2}^{2}\right)}} \right)} \\ \end{aligned}\end{split}\]

Hide code cell source

Math(aslatex({**MODEL.invariants, **MODEL.masses}))
\[\begin{split}\displaystyle \begin{aligned} \sigma_{1} \;&=\; m_{0}^{2} + m_{1}^{2} + m_{2}^{2} + m_{3}^{2} - \sigma_{2} - \sigma_{3} \\ \sigma_{2} \;&=\; m_{0}^{2} + m_{1}^{2} + m_{2}^{2} + m_{3}^{2} - \sigma_{1} - \sigma_{3} \\ \sigma_{3} \;&=\; m_{0}^{2} + m_{1}^{2} + m_{2}^{2} + m_{3}^{2} - \sigma_{1} - \sigma_{2} \\ m_{0} \;&=\; 2.28646 \\ m_{1} \;&=\; 0.938272046 \\ m_{2} \;&=\; 0.13957018 \\ m_{3} \;&=\; 0.493677 \\ \end{aligned}\end{split}\]

Numeric results#

intensity_expr = cached.xreplace(cached.unfold(MODEL), MODEL.variables)
intensity_expr = cached.xreplace(intensity_expr, MODEL.parameter_defaults)

Hide code cell source

free_symbols = intensity_expr.free_symbols
assert len(free_symbols) == 3
assert str(sorted(free_symbols, key=str)) == "[sigma1, sigma2, sigma3]"

Hide code cell source

intensity_funcs = {}
for s, s_expr in tqdm(MODEL.invariants.items()):
    k = int(str(s)[-1])
    s_expr = s_expr.xreplace(MODEL.masses).doit()
    expr = cached.doit(intensity_expr.xreplace({s: s_expr}))
    func = cached.lambdify(expr, backend="jax")
    assert len(func.argument_order) == 2, func.argument_order
    intensity_funcs[k] = func

Validation#

The following checks whether the full amplitude model has been deserialized correctly.

checksums = {
    checksum["point"]: checksum["value"]
    for checksum in WORKSPACE.checksums
    if checksum["distribution"] == INTENSITY_NAME
}
checksums
{'validation_point': 9345.853380852355}
checksum_points = {
    point["name"]: {par["name"]: par["value"] for par in point["parameters"]}
    for point in WORKSPACE.reference_points
}
checksum_points
{'validation_point': {'cos_theta_31': -0.2309352648098208,
  'phi_31': 0.0,
  'm_31': 1.9101377207489973,
  'cos_theta_31_2': 0.0,
  'phi_31_2': 0.0,
  'm_31_2': 2.28646},
 'validation_point_m12sq': {'m_12_sq': 3.2},
 'validation_point_m23sq': {'m_23_sq': 1.4},
 'validation_point_m31sq': {'m_31_sq': 3.2}}

Hide code cell source

theta31 = next(symbol for symbol in MODEL.variables if str(symbol) == "theta_31")
cos_theta31_expr = sp.cos(MODEL.variables[theta31]).xreplace(MODEL.parameter_defaults)
sigma1 = next(s for s in cos_theta31_expr.free_symbols if str(s) == "sigma1")
sigma2 = next(s for s in cos_theta31_expr.free_symbols if str(s) == "sigma2")
z = sp.Symbol("z", real=True)
s1_expr = sp.solve(sp.Eq(cos_theta31_expr, z), sigma1)[0]
array = []
for point_name, expected in checksums.items():
    parameters = checksum_points[point_name]
    s2 = parameters["m_31"] ** 2
    s1 = s1_expr.xreplace({z: parameters["cos_theta_31"], sigma2: s2})
    computed = intensity_funcs[3]({"sigma1": float(s1), "sigma2": s2})
    status = label_diff(expected - computed)
    array.append((INTENSITY_NAME, point_name, computed, expected, status))
intensity_validation = pd.DataFrame(
    array,
    columns=["Distribution", "Point", "Computed", "Expected", "Status"],
)
intensity_validation
Distribution Point Computed Expected Status
0 default_model validation_point 9345.853380852332 9345.853381 🟢

Hide code cell source

np.testing.assert_allclose(
    intensity_validation["Computed"].astype(float).to_numpy(),
    intensity_validation["Expected"].to_numpy(),
    atol=1e-10,
    rtol=0,
)

Dalitz plot#

Hide code cell source

i, j = (2, 1)
k, *_ = {1, 2, 3} - {i, j}
σk, σk_expr = list(MODEL.invariants.items())[k - 1]
Math(aslatex({σk: σk_expr}))
\[\begin{split}\displaystyle \begin{aligned} \sigma_{3} \;&=\; m_{0}^{2} + m_{1}^{2} + m_{2}^{2} + m_{3}^{2} - \sigma_{1} - \sigma_{2} \\ \end{aligned}\end{split}\]

Hide code cell source

resolution = 1_000
m = sorted(MODEL.masses, key=str)
x_min = float(((m[j] + m[k]) ** 2).xreplace(MODEL.masses))
x_max = float(((m[0] - m[i]) ** 2).xreplace(MODEL.masses))
y_min = float(((m[i] + m[k]) ** 2).xreplace(MODEL.masses))
y_max = float(((m[0] - m[j]) ** 2).xreplace(MODEL.masses))
x_diff = x_max - x_min
y_diff = y_max - y_min
x_min -= 0.05 * x_diff
x_max += 0.05 * x_diff
y_min -= 0.05 * y_diff
y_max += 0.05 * y_diff
X, Y = jnp.meshgrid(
    jnp.linspace(x_min, x_max, num=resolution),
    jnp.linspace(y_min, y_max, num=resolution),
)
dalitz_data = {
    f"sigma{i}": X,
    f"sigma{j}": Y,
}
intensities = intensity_funcs[k](dalitz_data)

Hide code cell source

assert not jnp.all(jnp.isnan(intensities)), "All intensities are NaN"

Hide code cell source

def get_decay_products(
    decay: ThreeBodyDecay, subsystem_id: FinalStateID
) -> tuple[State, State]:
    if subsystem_id not in decay.final_state:
        msg = f"Subsystem ID {subsystem_id} is not a valid final state ID"
        raise ValueError(msg)
    return tuple(s for s in decay.final_state.values() if s.index != subsystem_id)


plt.rc("font", size=18)
I_tot = jnp.nansum(intensities)
normalized_intensities = intensities / I_tot

fig, ax = plt.subplots(figsize=(14, 10))
mesh = ax.pcolormesh(X, Y, normalized_intensities, rasterized=True)
ax.set_aspect("equal")
c_bar = plt.colorbar(mesh, ax=ax, pad=0.01)
c_bar.ax.set_ylabel("Normalized intensity (a.u.)")
sigma_labels = {
    i: Rf"$\sigma_{i} = M^2\left({' '.join(p.latex for p in get_decay_products(DECAY, i))}\right)$"
    for i in (1, 2, 3)
}
ax.set_xlabel(sigma_labels[i])
ax.set_ylabel(sigma_labels[j])
plt.show()
_images/b47c955081ad151ff131f12b075d47b50a107724c5eb4aae027b782c8bf19ab6.svg