\(\Lambda_c^+ \to p \pi^+ K^-\)#

The decay chains follow the amplitude model that LHCb published for this decay. That model is distributed as a serialized JSON file, which Model serialization imports directly. Here, the qrules particle database is instead narrowed down to the twelve resonances of that model and their masses and widths are overwritten with the published values, so that the generated model stays as close to the published one as it can.

Hide code cell content

import logging
import warnings

import jax
import jax.numpy as jnp
import matplotlib.pyplot as plt
import numpy as np
import qrules
import sympy as sp
from IPython.display import Latex, Markdown
from matplotlib_inline.backend_inline import set_matplotlib_formats
from qrules.particle import create_particle
from tensorwaves.data.transform import SympyDataTransformer

from ampform_dpd import DalitzPlotDecompositionBuilder, DefinedExpression
from ampform_dpd.adapter.qrules import (
    load_particles,
    normalize_state_ids,
    to_three_body_decay,
)
from ampform_dpd.decay import State, ThreeBodyDecayChain
from ampform_dpd.dynamics import BreitWignerMinL
from ampform_dpd.dynamics.builder import create_mass_symbol, get_mandelstam_s
from ampform_dpd.io import (
    as_markdown_table,
    aslatex,
    cached,
    mute_ampform_warnings,
    simplify_latex_rendering,
    unfold_definitions,
)

jax.config.update("jax_enable_x64", True)
logging.getLogger("qrules.transition").setLevel(logging.ERROR)
logging.getLogger("absl").setLevel(logging.ERROR)  # mute JAX
mute_ampform_warnings()
set_matplotlib_formats("svg")
simplify_latex_rendering()
warnings.simplefilter("ignore", category=RuntimeWarning)

Decay definition#

Three of the published lineshapes are not relativistic Breit–Wigner functions. The \(\bar K^*_0(700)^0\) and \(\bar K^*_0(1430)^0\) use the Bugg parametrization, whose mass and width parameters are taken over as Breit–Wigner parameters here. The \(\Lambda(1405)\) is a two-channel Flatté, where the coupling of the \(\Sigma\pi\) channel reproduces a width of \(50.5\;\mathrm{MeV}\) at the pole, which is the value listed below.

Hide code cell source

PARTICLES = load_particles()
for name in [
    "K*(1410)~0",
    "K(2)*(1430)~0",
    "K*(1680)~0",
    "K(2)*(1980)~0",
    "Delta(1620)++",
    "Delta(1900)++",
    "Delta(1910)++",
    "Delta(1920)++",
    "Lambda(1800)",
    "Lambda(1810)",
    "Lambda(1890)",
]:
    PARTICLES.remove(PARTICLES[name])
LHCB_PARAMETERS = {
    "Delta(1232)++": (1.232, 0.117),
    "Delta(1600)++": (1.64, 0.3),
    "Delta(1700)++": (1.69, 0.38),
    "K(0)*(700)~0": (0.824, 0.478),
    "K*(892)~0": (0.8955, 0.0473),
    "K(0)*(1430)~0": (1.375, 0.19),
    "Lambda(1405)": (1.4051, 0.0505),
    "Lambda(1520)": (1.518, 0.0152),
    "Lambda(1600)": (1.63, 0.25),
    "Lambda(1670)": (1.67, 0.03),
    "Lambda(1690)": (1.69, 0.07),
    "Lambda(2000)": (1.988, 0.1793),
}
for name, (mass, width) in LHCB_PARAMETERS.items():
    particle = PARTICLES[name]
    PARTICLES.remove(particle)
    PARTICLES.add(create_particle(particle, mass=mass, width=width))
STM = qrules.StateTransitionManager(
    initial_state=["Lambda(c)+"],
    final_state=["p", "K-", "pi+"],
    mass_conservation_factor=3,
    allowed_intermediate_particles=["K", "Delta", "Lambda"],
    particle_db=PARTICLES,
    max_angular_momentum=2,
    formalism="canonical-helicity",
)
STM.set_allowed_interaction_types([qrules.InteractionType.STRONG], node_id=1)
problem_sets = STM.create_problem_sets()
REACTION = STM.find_solutions(problem_sets)
REACTION = normalize_state_ids(REACTION)
src = qrules.io.asmermaid(REACTION, collapse_graphs=True, markdown=True)
Markdown(src)
        flowchart LR
    T0_N0["$$\Lambda_{c}^{+}$$"]
    T0_1["$$1: p$$"]
    T0_2["$$2: K^{-}$$"]
    T0_3["$$3: \pi^{+}$$"]
    T0_N1@{ shape: text, label: " " }
    T0_4("$$\begin{gathered} \Lambda(1405) \\\ \Lambda(1520) \\\ \Lambda(1600) \\\ \Lambda(1670) \\\ \Lambda(1690) \\\ \Lambda(2000) \end{gathered}$$")
    T0_N0 --- T0_4
    T0_4 --- T0_N1
    T0_N0 --- T0_3
    T0_N1 --- T0_1
    T0_N1 --- T0_2
    T1_N0["$$\Lambda_{c}^{+}$$"]
    T1_1["$$1: p$$"]
    T1_2["$$2: K^{-}$$"]
    T1_3["$$3: \pi^{+}$$"]
    T1_N1@{ shape: text, label: " " }
    T1_4("$$\begin{gathered} \Delta(1232)^{++} \\\ \Delta(1600)^{++} \\\ \Delta(1700)^{++} \end{gathered}$$")
    T1_N0 --- T1_4
    T1_4 --- T1_N1
    T1_N0 --- T1_2
    T1_N1 --- T1_1
    T1_N1 --- T1_3
    T2_N0["$$\Lambda_{c}^{+}$$"]
    T2_1["$$1: p$$"]
    T2_2["$$2: K^{-}$$"]
    T2_3["$$3: \pi^{+}$$"]
    T2_N1@{ shape: text, label: " " }
    T2_4("$$\begin{gathered} \overline{K}_{0}^{*}(700)^{0} \\\ \overline{K}^{*}(892)^{0} \\\ \overline{K}_{0}^{*}(1430)^{0} \end{gathered}$$")
    T2_N0 --- T2_4
    T2_4 --- T2_N1
    T2_N0 --- T2_1
    T2_N1 --- T2_2
    T2_N1 --- T2_3
    

Hide code cell source

DECAY = to_three_body_decay(REACTION.transitions, min_ls=True)
Markdown(as_markdown_table([DECAY.initial_state, *DECAY.final_state.values()]))

index

name

LaTeX

\(J^P\)

mass (MeV)

width (MeV)

0

Lambda(c)+

\(\Lambda_{c}^{+}\)

\(\frac{1}{2}^+\)

2,286

0

1

p

\(p\)

\(\frac{1}{2}^+\)

938

0

2

K-

\(K^{-}\)

\(0^-\)

493

0

3

pi+

\(\pi^{+}\)

\(0^-\)

139

0

Hide code cell source

resonances = sorted(
    {t.resonance for t in DECAY.chains},
    key=lambda p: (p.name[0], p.mass),
)
resonance_names = [p.name for p in resonances]
Markdown(as_markdown_table(resonances))

name

LaTeX

\(J^P\)

mass (MeV)

width (MeV)

Delta(1232)++

\(\Delta(1232)^{++}\)

\(\frac{3}{2}^+\)

1,232

117

Delta(1600)++

\(\Delta(1600)^{++}\)

\(\frac{3}{2}^+\)

1,640

300

Delta(1700)++

\(\Delta(1700)^{++}\)

\(\frac{3}{2}^-\)

1,690

380

K(0)*(700)~0

\(\overline{K}_{0}^{*}(700)^{0}\)

\(0^+\)

824

478

K*(892)~0

\(\overline{K}^{*}(892)^{0}\)

\(1^-\)

895

47

K(0)*(1430)~0

\(\overline{K}_{0}^{*}(1430)^{0}\)

\(0^+\)

1,375

190

Lambda(1405)

\(\Lambda(1405)\)

\(\frac{1}{2}^-\)

1,405

50

Lambda(1520)

\(\Lambda(1520)\)

\(\frac{3}{2}^-\)

1,518

15

Lambda(1600)

\(\Lambda(1600)\)

\(\frac{1}{2}^+\)

1,630

250

Lambda(1670)

\(\Lambda(1670)\)

\(\frac{1}{2}^-\)

1,670

30

Lambda(1690)

\(\Lambda(1690)\)

\(\frac{3}{2}^-\)

1,690

70

Lambda(2000)

\(\Lambda(2000)\)

\(\frac{1}{2}^-\)

1,988

179

Hide code cell source

Latex(aslatex(DECAY, with_jp=True))
\[\begin{split}\begin{array}{c} \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=3/2]{L=1} \left(\Delta(1232)^{++}\left[\frac{3}{2}^+\right] \xrightarrow[S=1/2]{L=1} p\left[\frac{1}{2}^+\right] \pi^{+}\left[0^-\right]\right) K^{-}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=3/2]{L=1} \left(\Delta(1600)^{++}\left[\frac{3}{2}^+\right] \xrightarrow[S=1/2]{L=1} p\left[\frac{1}{2}^+\right] \pi^{+}\left[0^-\right]\right) K^{-}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=3/2]{L=1} \left(\Delta(1700)^{++}\left[\frac{3}{2}^-\right] \xrightarrow[S=1/2]{L=2} p\left[\frac{1}{2}^+\right] \pi^{+}\left[0^-\right]\right) K^{-}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\overline{K}_{0}^{*}(1430)^{0}\left[0^+\right] \xrightarrow[S=0]{L=0} K^{-}\left[0^-\right] \pi^{+}\left[0^-\right]\right) p\left[\frac{1}{2}^+\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\overline{K}_{0}^{*}(700)^{0}\left[0^+\right] \xrightarrow[S=0]{L=0} K^{-}\left[0^-\right] \pi^{+}\left[0^-\right]\right) p\left[\frac{1}{2}^+\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\overline{K}^{*}(892)^{0}\left[1^-\right] \xrightarrow[S=0]{L=1} K^{-}\left[0^-\right] \pi^{+}\left[0^-\right]\right) p\left[\frac{1}{2}^+\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\Lambda(1405)\left[\frac{1}{2}^-\right] \xrightarrow[S=1/2]{L=0} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=3/2]{L=1} \left(\Lambda(1520)\left[\frac{3}{2}^-\right] \xrightarrow[S=1/2]{L=2} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\Lambda(1600)\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=1} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\Lambda(1670)\left[\frac{1}{2}^-\right] \xrightarrow[S=1/2]{L=0} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=3/2]{L=1} \left(\Lambda(1690)\left[\frac{3}{2}^-\right] \xrightarrow[S=1/2]{L=2} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \Lambda_{c}^{+}\left[\frac{1}{2}^+\right] \xrightarrow[S=1/2]{L=0} \left(\Lambda(2000)\left[\frac{1}{2}^-\right] \xrightarrow[S=1/2]{L=0} p\left[\frac{1}{2}^+\right] K^{-}\left[0^-\right]\right) \pi^{+}\left[0^-\right] \\ \end{array}\end{split}\]

Lineshapes for dynamics#

Hide code cell source

s, m0, Γ0, m1, m2 = sp.symbols("s m0 Gamma0 m1 m2", nonnegative=True)
m_top, m_spec = sp.symbols(R"m_\mathrm{top} m_\mathrm{spectator}", nonnegative=True)
R_dec, R_prod = sp.symbols(R"R_\mathrm{res} R_{\Lambda_c}", nonnegative=True)
l_Λc, l_R = sp.symbols(R"l_{\Lambda_c} l_R", integer=True, nonnegative=True)
bw = BreitWignerMinL(s, m_top, m_spec, m0, Γ0, m1, m2, l_R, l_Λc, R_dec, R_prod)
Latex(aslatex(unfold_definitions(bw)))
\[\begin{split}\begin{aligned} \mathcal{R}^\mathrm{BW}_{l_{R},l_{\Lambda_c}}\left(s\right) \;&=\; \frac{\frac{\mathcal{F}_{l_{R}}\left(s, m_{1}, m_{2}\right)}{\mathcal{F}_{l_{R}}\left(m_{0}^{2}, m_{1}, m_{2}\right)} \frac{\mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, \sqrt{s}, m_\mathrm{spectator}\right)}{\mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, m_{0}, m_\mathrm{spectator}\right)}}{m_{0}^{2} - i m_{0} \Gamma_{0}\left(s\right) - s} \\ \frac{\frac{\mathcal{F}_{l_{R}}\left(s, m_{1}, m_{2}\right)}{\mathcal{F}_{l_{R}}\left(m_{0}^{2}, m_{1}, m_{2}\right)} \frac{\mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, \sqrt{s}, m_\mathrm{spectator}\right)}{\mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, m_{0}, m_\mathrm{spectator}\right)}}{m_{0}^{2} - i m_{0} \Gamma_{0}\left(s\right) - s} \;&=\; \frac{\mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, \sqrt{s}, m_\mathrm{spectator}\right) \mathcal{F}_{l_{R}}\left(s, m_{1}, m_{2}\right)}{\left(m_{0}^{2} - i m_{0} \Gamma_{0}\left(s\right) - s\right) \mathcal{F}_{l_{R}}\left(m_{0}^{2}, m_{1}, m_{2}\right) \mathcal{F}_{l_{\Lambda_c}}\left(m_\mathrm{top}^{2}, m_{0}, m_\mathrm{spectator}\right)} \\ \mathcal{F}_{l_{\Lambda_c}}\left(s, m_{1}, m_\mathrm{spectator}\right) \;&=\; \sqrt{B_{l_{\Lambda_c}}^2\left(R_{\Lambda_c}^{2} q^2\left(s\right)\right)} \\ \Gamma_{0}\left(s\right) \;&=\; \frac{\Gamma_{0} \mathcal{F}_{l_{R}}\left(s, m_{1}, m_{2}\right)^{2} \rho\left(s\right)}{\mathcal{F}_{l_{R}}\left(m_{0}^{2}, m_{1}, m_{2}\right)^{2} \rho_{0}\left(m_{0}^{2}\right)} \\ B_{l_{\Lambda_c}}^2\left(R_{\Lambda_c}^{2} q^2\left(s\right)\right) \;&=\; \frac{\left|{h_{l_{\Lambda_c}}^{(1)}\left(1\right)}\right|^{2}}{R_{\Lambda_c}^{2} \left|{h_{l_{\Lambda_c}}^{(1)}\left(R_{\Lambda_c} \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_\mathrm{spectator}\right)^{2}\right) \left(s - \left(m_{1} + m_\mathrm{spectator}\right)^{2}\right)}{4 s} \\ \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} \\ h_{l_{\Lambda_c}}^{(1)}\left(z\right) \;&=\; \frac{\left(- i\right)^{l_{\Lambda_c} + 1} e^{i z} \sum_{k=0}^{l_{\Lambda_c}} \frac{\left(\frac{i}{2 z}\right)^{k} \left(k + l_{\Lambda_c}\right)!}{k! \left(- k + l_{\Lambda_c}\right)!}}{z} \\ \end{aligned}\end{split}\]

Hide code cell source

def formulate_breit_wigner(
    decay_chain: ThreeBodyDecayChain,
) -> DefinedExpression:
    s = get_mandelstam_s(decay_chain.decay_node)
    child1_mass, child2_mass = map(create_mass_symbol, decay_chain.decay_products)
    assert decay_chain.incoming_ls is not None
    assert decay_chain.outgoing_ls is not None
    l_prod = sp.Rational(decay_chain.incoming_ls.L)
    l_dec = sp.Rational(decay_chain.outgoing_ls.L)
    parent_mass = sp.Symbol(f"m_{{{decay_chain.parent.latex}}}", nonnegative=True)
    spectator_mass = sp.Symbol(f"m_{{{decay_chain.spectator.latex}}}", nonnegative=True)
    resonance_mass = sp.Symbol(f"m_{{{decay_chain.resonance.latex}}}", nonnegative=True)
    resonance_width = sp.Symbol(
        Rf"\Gamma_{{{decay_chain.resonance.latex}}}", nonnegative=True
    )
    R_dec = sp.Symbol(R"R_\mathrm{res}", nonnegative=True)
    R_prod = sp.Symbol(R"R_{\Lambda_c}", nonnegative=True)
    parameter_defaults = {
        parent_mass: decay_chain.parent.mass,
        spectator_mass: decay_chain.spectator.mass,
        resonance_mass: decay_chain.resonance.mass,
        resonance_width: decay_chain.resonance.width,
        child1_mass: decay_chain.decay_products[0].mass,
        child2_mass: decay_chain.decay_products[1].mass,
        # https://github.com/ComPWA/polarimetry/pull/11#issuecomment-1128784376
        R_dec: 1.5,
        R_prod: 5,
    }
    dynamics = BreitWignerMinL(
        s,
        parent_mass,
        spectator_mass,
        resonance_mass,
        resonance_width,
        child1_mass,
        child2_mass,
        l_dec,
        l_prod,
        R_dec,
        R_prod,
    )
    return DefinedExpression(dynamics, parameter_defaults)

Model formulation#

model_builder = DalitzPlotDecompositionBuilder(DECAY, min_ls=(False, True))
for chain in model_builder.decay.chains:
    model_builder.dynamics_choices.register_builder(chain, formulate_breit_wigner)
model = model_builder.formulate()
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(3)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{1(3)}\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(3)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{2(3)}\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(3)}\right) d^{\frac{1}{2}}_{\lambda_{0},\lambda_0^{\prime}}\left(\zeta^0_{3(3)}\right)}}\right|^{2}}\]

Notice that subsystem 3 (with the \(\varLambda^*\) resonances and recoil \(\pi^+\)) has been selected as reference subsystem for the Wigner‑\(d\) functions. This is computationally more efficient, because that subsystem contains most resonances.

Latex(aslatex(model.amplitudes, terms_per_line=1))
\[\begin{split}\begin{aligned} A^{1}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{1,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}^{*}(892)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}^{*}(892)^{0}, 0, 0} d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(1430)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(1430)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(700)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(700)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1232)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1232)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1600)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1600)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1700)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1700)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ A^{3}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1405), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1405), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1600), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1600), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1670), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1670), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(2000), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(2000), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1520), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1520), - \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1690), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1690), - \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ A^{1}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{1,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}^{*}(892)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}^{*}(892)^{0}, 0, 0} d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(1430)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(1430)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(700)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(700)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1232)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1232)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1600)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1600)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1700)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1700)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ A^{3}_{- \frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1405), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1405), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1600), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1600), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1670), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1670), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(2000), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(2000), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1520), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1520), \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{- \frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1690), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1690), \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ A^{1}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{1,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}^{*}(892)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}^{*}(892)^{0}, 0, 0} d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(1430)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(1430)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(700)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(700)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1232)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1232)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1600)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1600)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{- \frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1700)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1700)^{++}, 0, - \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ A^{3}_{\frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1405), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1405), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1600), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1600), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1670), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1670), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(2000), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(2000), - \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1520), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1520), - \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1690), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1690), - \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ A^{1}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1}^{1}{\frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{1,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}^{*}(892)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}^{*}(892)^{0}, 0, 0} d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(1430)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(1430)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{\frac{\sqrt{2} \delta_{\frac{1}{2}, \lambda_{R} - \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} - \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} - \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},- \frac{1}{2}} \mathcal{H}^\mathrm{LS,production}_{\overline{K}_{0}^{*}(700)^{0}, 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\overline{K}_{0}^{*}(700)^{0}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ A^{2}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1232)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1232)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1600)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1600)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{2}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Delta(1700)^{++}, 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Delta(1700)^{++}, 0, \frac{1}{2}} d^{\frac{3}{2}}_{\lambda_{R},- \frac{1}{2}}\left(\theta_{31}\right)}{2}} \\ A^{3}_{\frac{1}{2}, \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1405), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1405), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1600), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1600), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1670), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1670), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-1/2}^{1/2}{\frac{\sqrt{2} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{0,0,\frac{1}{2},\lambda_{R}} C^{\frac{1}{2},\lambda_{R}}_{\frac{1}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(2000), 0, \frac{1}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(2000), \frac{1}{2}, 0} d^{\frac{1}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1520), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1520), \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=-3/2}^{3/2}{\frac{\sqrt{6} \delta_{\frac{1}{2} \lambda_{R}} \mathcal{R}^\mathrm{BW}_{2,1}\left(\sigma_{3}\right) C^{\frac{1}{2},\lambda_{R}}_{1,0,\frac{3}{2},\lambda_{R}} C^{\frac{3}{2},\lambda_{R}}_{\frac{3}{2},\lambda_{R},0,0} \mathcal{H}^\mathrm{LS,production}_{\Lambda(1690), 1, \frac{3}{2}} \mathcal{H}^\mathrm{decay}_{\Lambda(1690), \frac{1}{2}, 0} d^{\frac{3}{2}}_{\lambda_{R},\frac{1}{2}}\left(\theta_{12}\right)}{2}} \\ \end{aligned}\end{split}\]

By default, the aligned amplitudes are built up of Wigner-\(d\) functions, Clebsch–Gordan coefficients (\(C\)), a resonance parametrization (\(\mathcal{R}(\sigma)\)), and two coupling symbols \(\mathcal{H}^\text{prod}_{\dots}, \mathcal{H}^\text{dec}_{\dots}\). In some cases, you want to combine the couplings into one scaling factor. That can be done with the use_coefficients flag. The following example also explicitly specifies the reference subsystem (which is otherwise automatically chosen).

coefficient_model = model_builder.formulate(
    reference_subsystem=1, use_coefficients=True
)

Hide code cell source

(symbol, expr), *_ = coefficient_model.amplitudes.items()
Latex(aslatex({symbol: expr}, terms_per_line=1))
\[\begin{split}\begin{aligned} A^{1}_{- \frac{1}{2}, - \frac{1}{2}, 0, 0} \;&=\; \sum_{\lambda_{R}=-1}^{1}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{1,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{1,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{\lambda,LS,\overline{K}^{*}(892)^{0}}_{0, \frac{1}{2}, 0, 0} d^{1}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{\lambda,LS,\overline{K}_{0}^{*}(1430)^{0}}_{0, \frac{1}{2}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \;&+\; \sum_{\lambda_{R}=0}{- \frac{\sqrt{2} \delta_{- \frac{1}{2}, \lambda_{R} + \frac{1}{2}} \mathcal{R}^\mathrm{BW}_{0,0}\left(\sigma_{1}\right) C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,0,\frac{1}{2},\lambda_{R} + \frac{1}{2}} C^{\frac{1}{2},\lambda_{R} + \frac{1}{2}}_{0,\lambda_{R},\frac{1}{2},\frac{1}{2}} \mathcal{H}^\mathrm{\lambda,LS,\overline{K}_{0}^{*}(700)^{0}}_{0, \frac{1}{2}, 0, 0} d^{0}_{\lambda_{R},0}\left(\theta_{23}\right)}{2}} \\ \end{aligned}\end{split}\]

Dalitz plot#

The intensity can now be rendered over the Dalitz plane, just as in Model serialization, with the difference that the model here was formulated from the ThreeBodyDecay above instead of imported from a serialized JSON file. There is therefore no compiler that infers the kinematics from serialized metadata: the amplitudes are written in terms of helicity angles, so the two Mandelstam variables that span the plane first have to be converted with a SympyDataTransformer.

Note that the couplings all carry their default value of \(1\), so the relative strength of the resonances is not the one that LHCb fitted. It is the lineshapes and the phase space that carry the published values.

The third Mandelstam variable is fixed by the other two.

Hide code cell source

i, j = (3, 1)
k, *_ = {1, 2, 3} - {i, j}
σk, σk_expr = list(model.invariants.items())[k - 1]
Latex(aslatex({σk: σk_expr}))
\[\begin{split}\begin{aligned} \sigma_{2} \;&=\; m_{0}^{2} + m_{1}^{2} + m_{2}^{2} + m_{3}^{2} - \sigma_{1} - \sigma_{3} \\ \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),
)

Hide code cell source

definitions = dict(model.variables)
definitions[σk] = σk_expr
definitions = {
    symbol: expr.xreplace(definitions).xreplace(model.masses)
    for symbol, expr in definitions.items()
}
data_transformer = SympyDataTransformer.from_sympy(definitions, backend="jax")
dalitz_data = {
    f"sigma{i}": X,
    f"sigma{j}": Y,
}
dalitz_data.update(data_transformer(dalitz_data))
intensity_expr = cached.xreplace(
    cached.unfold(model),
    {**model.parameter_defaults, **model.masses},
)
intensity_func = cached.lambdify(intensity_expr, backend="jax")
intensities = intensity_func(dalitz_data)

Hide code cell source

def get_decay_products(subsystem_id: int) -> tuple[State, State]:
    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, aspect=25, ax=ax, pad=0.01, shrink=0.66)
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(i))}\right)$"
    for i in (1, 2, 3)
}
ax.set_xlabel(sigma_labels[i])
ax.set_ylabel(sigma_labels[j])
plt.show()
_images/b1dd7e5ff773a44cf9cca2b520a6dbd2a7d8f2c4b3368cc01985add7a9b0b47a.svg

Tip

Compare this Dalitz plot, built directly from the reaction, with the one produced from the deserialized LHCb model in Model serialization.