From 700b8fbdb3740514d3390047a4ffc1c68761a1f2 Mon Sep 17 00:00:00 2001 From: Remco de Boer <29308176+redeboer@users.noreply.github.com> Date: Thu, 3 Oct 2024 20:02:51 +0200 Subject: [PATCH] MAINT: freeze TR-030 --- docs/030.ipynb | 256 ++++++++++++++++++++++++++++++++++++++++++++----- docs/conf.py | 1 + 2 files changed, 233 insertions(+), 24 deletions(-) diff --git a/docs/030.ipynb b/docs/030.ipynb index 8a6cd24..4508c43 100644 --- a/docs/030.ipynb +++ b/docs/030.ipynb @@ -72,6 +72,7 @@ "import re\n", "from collections import defaultdict\n", "from functools import lru_cache\n", + "from pathlib import Path\n", "from typing import Any\n", "\n", "import ampform\n", @@ -196,13 +197,26 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "remove-output" ] }, "outputs": [], "source": [ "dot = qrules.io.asdot(reaction, collapse_graphs=True)\n", - "graphviz.Source(dot)" + "graph = graphviz.Source(dot)\n", + "output_file = Path(\"030/qrules-output\")\n", + "output_file.parent.mkdir(exist_ok=True)\n", + "graph.render(output_file, format=\"svg\")\n", + "output_file.unlink()\n", + "graph" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/431c2299-2ec1-41ee-8f6f-5f3e5c036024)" ] }, { @@ -266,8 +280,26 @@ { "cell_type": "code", "execution_count": null, - "metadata": {}, - "outputs": [], + "metadata": { + "tags": [ + "keep_output" + ] + }, + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\sum_{m_{A}=-1}^{1} \\sum_{m_{1}=-1/2}^{1/2} \\sum_{m_{2}=-1/2}^{1/2}{\\left|{A^{01}_{m_{A}, 0, m_{1}, m_{2}}}\\right|^{2}}$" + ], + "text/plain": [ + "PoolSum(Abs(A^01[m_A, 0, m1, m2])**2, (m_A, (0, 1, -1)), (m1, (-1/2, 1/2)), (m2, (-1/2, 1/2)))" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "model_builder = ampform.get_builder(reaction)\n", "model_builder.adapter.permutate_registered_topologies()\n", @@ -289,10 +321,32 @@ }, "tags": [ "hide-input", - "full-width" + "full-width", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\begin{array}{rcl}\n", + " A^{01}_{0, 0, - \\frac{1}{2}, - \\frac{1}{2}} &=& - C_{J/\\psi(1S) \\to {N_1(3/2^-)}_{+1/2} \\overline{p}_{+1/2}; N_1(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,0}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{- \\frac{1}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + " &+& - C_{J/\\psi(1S) \\to {N_1(3/2^-)}_{+1/2} \\overline{p}_{-1/2}; N_1(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,1}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{\\frac{1}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + " &+& - C_{J/\\psi(1S) \\to {N_1(3/2^-)}_{+3/2} \\overline{p}_{+1/2}; N_1(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,-1}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{- \\frac{3}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + " &+& - C_{J/\\psi(1S) \\to {N_2(3/2^-)}_{+1/2} \\overline{p}_{+1/2}; N_2(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,0}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{- \\frac{1}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + " &+& - C_{J/\\psi(1S) \\to {N_2(3/2^-)}_{+1/2} \\overline{p}_{-1/2}; N_2(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,1}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{\\frac{1}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + " &+& - C_{J/\\psi(1S) \\to {N_2(3/2^-)}_{+3/2} \\overline{p}_{+1/2}; N_2(3/2^-) \\to \\eta_{0} p_{+1/2}} X_{J^P={\\frac{3}{2}}^{+}, Q=+1} D^{1}_{0,-1}\\left(- \\phi_{01},\\theta_{01},0\\right) D^{\\frac{3}{2}}_{- \\frac{3}{2},\\frac{1}{2}}\\left(- \\phi^{01}_{0},\\theta^{01}_{0},0\\right) \\\\\n", + "\\end{array}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "selected_amplitudes = {\n", " k: v for i, (k, v) in enumerate(model.amplitudes.items()) if i == 0\n", @@ -308,10 +362,29 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\begin{array}{cll}\n", + " X_{J^P={\\frac{3}{2}}^{+}, Q=+1} \\\\\n", + " N_2(3/2^-) & m=1.92\\text{ GeV} & \\Gamma=0.6\\text{ GeV} \\\\\n", + " N_1(3/2^-) & m=1.82\\text{ GeV} & \\Gamma=0.6\\text{ GeV} \\\\\n", + "\\end{array}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "src = R\"\\begin{array}{cll}\" \"\\n\"\n", "for symbol, resonances in create_dynamics_symbol.collected_symbols.items():\n", @@ -412,10 +485,29 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\begin{array}{rcl}\n", + " \\rho^\\mathrm{CM}_{m_{1},m_{2}}\\left(s\\right) &=& - 16 i \\pi \\Sigma\\left(s\\right) \\\\\n", + " \\Sigma\\left(s\\right) &=& \\frac{- \\left(m_{1}^{2} - m_{2}^{2}\\right) \\left(- \\frac{1}{\\left(m_{1} + m_{2}\\right)^{2}} + \\frac{1}{s}\\right) \\log{\\left(\\frac{m_{1}}{m_{2}} \\right)} + \\frac{2 \\log{\\left(\\frac{\\left|{m_{1}^{2} + m_{2}^{2} + 2 \\sqrt{s} q\\left(s\\right) - s}\\right|}{2 m_{1} m_{2}} \\right)} q\\left(s\\right)}{\\sqrt{s}}}{16 \\pi^{2}} \\\\\n", + " q\\left(s\\right) &=& \\frac{\\sqrt{\\lambda\\left(s, m_{1}^{2}, m_{2}^{2}\\right)}}{2 \\sqrt{s}} \\\\\n", + "\\end{array}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "s, m1, m2 = sp.symbols(\"s m1 m2\", nonnegative=True)\n", "exprs = [\n", @@ -478,10 +570,27 @@ }, "tags": [ "hide-input", - "full-width" + "full-width", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\begin{array}{rcl}\n", + " X_{J^P={\\frac{3}{2}}^{+}, Q=+1} &=& \\frac{\\Gamma_{N_1(3/2^-)} \\beta_{N_1(3/2^-)} m_{N_1(3/2^-)}}{- \\Gamma_{N_1(3/2^-)} m_{N_1(3/2^-)} \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right) - m_{01}^{2} + \\left(m_{N_1(3/2^-)}\\right)^{2}} + \\frac{\\Gamma_{N_2(3/2^-)} \\beta_{N_2(3/2^-)} m_{N_2(3/2^-)}}{- \\Gamma_{N_2(3/2^-)} m_{N_2(3/2^-)} \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right) - m_{01}^{2} + \\left(m_{N_2(3/2^-)}\\right)^{2}} \\\\\n", + "\\end{array}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "dynamics_expressions_bw = {\n", " symbol: formulate_breit_wigner(resonances)\n", @@ -582,10 +691,27 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\begin{array}{rcl}\n", + " X_{J^P={\\frac{3}{2}}^{+}, Q=+1} &=& \\frac{\\frac{\\beta_{N_1(3/2^-)} g_{N_1(3/2^-)}}{- m_{01}^{2} + \\left(m_{N_1(3/2^-)}\\right)^{2}} + \\frac{\\beta_{N_2(3/2^-)} g_{N_2(3/2^-)}}{- m_{01}^{2} + \\left(m_{N_2(3/2^-)}\\right)^{2}}}{- \\left(\\frac{\\left(g_{N_1(3/2^-)}\\right)^{2}}{- m_{01}^{2} + \\left(m_{N_1(3/2^-)}\\right)^{2}} + \\frac{\\left(g_{N_2(3/2^-)}\\right)^{2}}{- m_{01}^{2} + \\left(m_{N_2(3/2^-)}\\right)^{2}}\\right) \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right) + 1} \\\\\n", + "\\end{array}$" + ], + "text/plain": [ + "" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "dynamics_expressions_fvector = {\n", " symbol: formulate_f_vector(resonances)\n", @@ -663,10 +789,25 @@ }, "tags": [ "full-width", - "hide-input" + "hide-input", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\frac{\\Gamma_{N_1(3/2^-)} \\beta_{N_1(3/2^-)} m_{N_1(3/2^-)}}{- \\Gamma_{N_1(3/2^-)} m_{N_1(3/2^-)} \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right) - m_{01}^{2} + \\left(m_{N_1(3/2^-)}\\right)^{2}} + \\frac{\\Gamma_{N_2(3/2^-)} \\beta_{N_2(3/2^-)} m_{N_2(3/2^-)}}{- \\Gamma_{N_2(3/2^-)} m_{N_2(3/2^-)} \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right) - m_{01}^{2} + \\left(m_{N_2(3/2^-)}\\right)^{2}}$" + ], + "text/plain": [ + "\\Gamma_{N_1(3/2^-)}*\\beta_{N_1(3/2^-)}*m_{N_1(3/2^-)}/(-\\Gamma_{N_1(3/2^-)}*m_{N_1(3/2^-)}*PhaseSpaceCM(m_01**2, m_0, m_1) - m_01**2 + m_{N_1(3/2^-)}**2) + \\Gamma_{N_2(3/2^-)}*\\beta_{N_2(3/2^-)}*m_{N_2(3/2^-)}/(-\\Gamma_{N_2(3/2^-)}*m_{N_2(3/2^-)}*PhaseSpaceCM(m_01**2, m_0, m_1) - m_01**2 + m_{N_2(3/2^-)}**2)" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "dynamics_expr_bw, *_ = dynamics_expressions_bw.values()\n", "dynamics_expr_bw" @@ -681,10 +822,25 @@ }, "tags": [ "full-width", - "hide-input" + "hide-input", + "keep_output" ] }, - "outputs": [], + "outputs": [ + { + "data": { + "text/latex": [ + "$\\displaystyle \\frac{\\beta_{N_1(3/2^-)} g_{N_1(3/2^-)} \\left(m_{01}^{2} - \\left(m_{N_2(3/2^-)}\\right)^{2}\\right) + \\beta_{N_2(3/2^-)} g_{N_2(3/2^-)} \\left(m_{01}^{2} - \\left(m_{N_1(3/2^-)}\\right)^{2}\\right)}{- \\left(m_{01}^{2} - \\left(m_{N_1(3/2^-)}\\right)^{2}\\right) \\left(m_{01}^{2} - \\left(m_{N_2(3/2^-)}\\right)^{2}\\right) + \\left(- \\left(g_{N_1(3/2^-)}\\right)^{2} \\left(m_{01}^{2} - \\left(m_{N_2(3/2^-)}\\right)^{2}\\right) - \\left(g_{N_2(3/2^-)}\\right)^{2} \\left(m_{01}^{2} - \\left(m_{N_1(3/2^-)}\\right)^{2}\\right)\\right) \\rho^\\mathrm{CM}_{m_{0},m_{1}}\\left(m_{01}^{2}\\right)}$" + ], + "text/plain": [ + "(\\beta_{N_1(3/2^-)}*g_{N_1(3/2^-)}*(m_01**2 - m_{N_2(3/2^-)}**2) + \\beta_{N_2(3/2^-)}*g_{N_2(3/2^-)}*(m_01**2 - m_{N_1(3/2^-)}**2))/(-(m_01**2 - m_{N_1(3/2^-)}**2)*(m_01**2 - m_{N_2(3/2^-)}**2) + (-g_{N_1(3/2^-)}**2*(m_01**2 - m_{N_2(3/2^-)}**2) - g_{N_2(3/2^-)}**2*(m_01**2 - m_{N_1(3/2^-)}**2))*PhaseSpaceCM(m_01**2, m_0, m_1))" + ] + }, + "execution_count": null, + "metadata": {}, + "output_type": "execute_result" + } + ], "source": [ "dynamics_expr_fvector, *_ = dynamics_expressions_fvector.values()\n", "dynamics_expr_fvector.simplify(doit=False)" @@ -943,7 +1099,8 @@ }, "tags": [ "hide-input", - "scroll-input" + "scroll-input", + "remove-output" ] }, "outputs": [], @@ -1001,10 +1158,21 @@ "\n", "ax.set_ylim(0, None)\n", "fig.legend(loc=\"upper right\")\n", - "plt.tight_layout()\n", + "fig.tight_layout()\n", + "\n", + "output_file = Path(\"030/sub-intensities.svg\")\n", + "output_file.parent.mkdir(exist_ok=True)\n", + "fig.savefig(output_file, bbox_inches=\"tight\")\n", "plt.show()" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/116a418a-4394-4203-90ea-11d6c58c8aea)" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -1120,6 +1288,10 @@ " )\n", "\n", " ax1.legend(loc=\"upper left\")\n", + "\n", + " output_file = Path(f\"030/argand-{title.lower().replace(' ', '-')}.svg\")\n", + " output_file.parent.mkdir(exist_ok=True)\n", + " fig.savefig(output_file, bbox_inches=\"tight\")\n", " fig.show()" ] }, @@ -1128,7 +1300,8 @@ "execution_count": null, "metadata": { "tags": [ - "hide-input" + "hide-input", + "remove-output" ] }, "outputs": [], @@ -1140,12 +1313,20 @@ ")" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/4c4060a0-dd38-4e3a-8e4c-04ac0a26f5ac)" + ] + }, { "cell_type": "code", "execution_count": null, "metadata": { "tags": [ - "hide-input" + "hide-input", + "remove-output" ] }, "outputs": [], @@ -1157,6 +1338,13 @@ ")" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/10816903-530f-4990-9cc6-0cb35269368f)" + ] + }, { "cell_type": "markdown", "metadata": {}, @@ -1270,6 +1458,10 @@ " # Add legends\n", " fig.legend(bbox_to_anchor=(0.1, 0.9), loc=\"upper left\")\n", " fig.tight_layout()\n", + "\n", + " output_file = Path(f\"030/phase-{title.lower().replace(' ', '-')}.svg\")\n", + " output_file.parent.mkdir(exist_ok=True)\n", + " fig.savefig(output_file, bbox_inches=\"tight\")\n", " fig.show()" ] }, @@ -1281,7 +1473,8 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "remove-output" ] }, "outputs": [], @@ -1295,6 +1488,13 @@ ")" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/a98341f5-9173-4aeb-a817-e9687a0bc7ec)" + ] + }, { "cell_type": "code", "execution_count": null, @@ -1303,7 +1503,8 @@ "source_hidden": true }, "tags": [ - "hide-input" + "hide-input", + "remove-output" ] }, "outputs": [], @@ -1316,6 +1517,13 @@ " title=\"Breit-Wigner\",\n", ")" ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "![Image](https://github.com/user-attachments/assets/c718704d-2648-443e-bb44-4693f6da7e6c)" + ] } ], "metadata": { diff --git a/docs/conf.py b/docs/conf.py index 81eedde..11529c6 100644 --- a/docs/conf.py +++ b/docs/conf.py @@ -41,6 +41,7 @@ def get_nb_exclusion_patterns() -> list[str]: "021*", "022*", "028*", + "030*", "031*", "033*", }