Conversation
redeboer
force-pushed
the
FEAT/final-state-polarization
branch
4 times, most recently
from
July 31, 2026 08:11
d6d12d4 to
6703753
Compare
Add a polarized_final_states argument to DalitzPlotDecompositionBuilder.formulate() that attaches one additional Wigner-D function per amplitude for each selected final-state particle, connecting its helicity (summed coherently) to the helicity of its decay product (summed incoherently), weighted by final-state decay couplings H^fs. With the default couplings of 1, the intensity is identical to the model without the additional Wigner-D functions; for weak decays such as Sigma+ -> p pi0, the couplings follow from the decay-asymmetry parameter via the new formulate_weak_decay_couplings() helper. This makes production-coupling phases visible that are flat directions in a Dalitz-plot-only fit. The new polarization example notebook demonstrates this on J/psi -> K0 Sigma+ pbar: the proton direction in the Sigma+ rest frame becomes sensitive to the Sigma+ polarization, while the Dalitz distribution is unchanged.
The decay angles (phi_k, theta_k) of a polarized final-state particle cannot be computed from the Mandelstam variables, since the decay product lies outside the three-body decay plane. Formulate them with AmpForm's compute_helicity_angles() over an extended four-body topology in which final-state particle k decays further, and include the resulting boost-chain expressions in AmplitudeModel.variables. The expressions take four-momentum arrays p1, p2, p3 and q<k> (the analyzed decay product) in the center-of-momentum frame as input, following the boost chain of the reference subsystem.
redeboer
force-pushed
the
FEAT/final-state-polarization
branch
from
July 31, 2026 08:13
ea66d3b to
cee6354
Compare
redeboer
changed the title
redeboer
changed the title
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
3 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
✨ New features
polarized_final_statestoDalitzPlotDecompositionBuilder.formulate()to measure final-state polarization through subsequent decays. The model sums each selected particle's helicity coherently with a Wigner-$D$ function and introduces decay couplings; unit default couplings recover the unpolarized intensity.ampform_dpd.polarizationhelpers for decay-angle and coupling symbols, weak two-body decay couplings from the asymmetry parameter, and decay angles expressed through four-momenta. The angle definitions follow the selected reference subsystem's boost chain for both spectators and isobar decay products and are included in the model variables.📝 Documentation
Squash commit messages