DOC: manual formulation of helicity amplitudes for pgamma->LambdaKPi - #104
DOC: manual formulation of helicity amplitudes for pgamma->LambdaKPi#104shenvitor wants to merge 13 commits into
Conversation
shenvitor
added
✨ Feature
redeboer
left a comment
redeboer
left a comment
There was a problem hiding this comment.
Looks like going in the right direction 💪
See comments and be careful to rename the PR to a conventional commit message
| "zeta_1_0, zeta_1_1, zeta_2_0, zeta_2_1, zeta_3_0, zeta_3_1 = sp.symbols(\n", | ||
| " r\"\\zeta_{1(1)}^0 \\zeta_{1(1)}^1 \\zeta_{2(1)}^0 \\zeta_{2(1)}^1 \\zeta_{3(1)}^0 \\zeta_{3(1)}^1\"\n", | ||
| ")\n", |
There was a problem hiding this comment.
In this case I would use sp.Symbol and define each zeta separately. It's more scalable and helps finding bugs if each symbol definition is below the other.
| "A1 = (\n", | ||
| " formualte_A1(\"K^{*}(1410)^+\", 1)\n", | ||
| " + formualte_A1(\"K^{*}(1680)^+\", 1)\n", | ||
| " + formualte_A1(\"K^{*}(892)^+\", 1)\n", | ||
| " + formualte_A1(\"K^{*}_0(1430)^+\", 0)\n", | ||
| " + formualte_A1(\"K^{*}_0(700)\", 0)\n", | ||
| " + formualte_A1(\"K^{*}_2(1430)\", 2)\n", | ||
| " + formualte_A1(\"K^{*}_3(1780)\", 3)\n", | ||
| " + formualte_A1(\"K^{*}_4(2045)\", 4)\n", | ||
| ")\n", | ||
| "A1" |
There was a problem hiding this comment.
This will need to be automated, i.e., some loop over the formulate_A1() is more something like formulate_chain_amplitude_a1() and formulate_a1() creates the coherent sum of those chain amplitudes.
| " r\"s_{31} m_{\\Sigma^*} \\Gamma_{\\Sigma^*} C_{\\Sigma^*}\"\n", | ||
| ")\n", | ||
| "theta31, phi31, delta = sp.symbols(r\"theta_31 phi_31 \\delta_{\\lambda_\\Sigma^*1/2}\")\n", | ||
| "d_s = sp.Function(r\"d_{\\lambda_\\Sigma^*1/2}^{1/2}\")\n", |
There was a problem hiding this comment.
Replace with Wigner.d and use directly inline in the expression definition
| "metadata": {}, | ||
| "outputs": [], | ||
| "source": [ | ||
| "A1 = (\n", |
There was a problem hiding this comment.
Be careful not to overwrite the previous A1 definition. This should be the A1_expr, or even something like a mapping of sp.Symbol to sp.Expr as dictionary. Something like:
amplitudes = {}
amplitudes[A1[-half, +half]] = formulate_A1(...) + ...
shenvitor
changed the title
2 participants
Uh oh!
There was an error while loading. Please reload this page.