Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DOC: show how to collect Wigner D-functions #38

Draft
wants to merge 3 commits into
base: main
Choose a base branch
from

Conversation

redeboer
Copy link
Member

Hi @mmikhasenko, here's an example of how to 'hack' the dynamics builder so that it produces a symbol instead of an expression for the dynamics. As you said, this can then be done to collect the angular functions or one can use the same expression to insert dynamics lineshapes later on.

It's just a sketch, but you can see what it looks like in the collect-wigners.ipynb notebook.

image

@redeboer redeboer added the 📝 Docs Improvements or additions to documentation label Jan 13, 2023
@redeboer redeboer self-assigned this Jan 13, 2023
@mmikhasenko
Copy link
Collaborator

Nice! I think it might make a very popular application of the ampform-dpd. So, cool to have this example.
The expressions can be simplify-ed and trigsimp-ed.
Where are the non-digonal terms?

@redeboer
Copy link
Member Author

The expressions can be simplify-ed and trigsimp-ed.

Cool O.o

$$ \begin{eqnarray} X_{0}^{0,1/2;0,0}\left(\sigma_{1}\right) &:& 1 \\ X_{1}^{0,1/2;1,0}\left(\sigma_{1}\right) &:& \frac{1}{3} \\ X_{2}^{1,3/2;2,0}\left(\sigma_{1}\right) &:& \frac{2}{5} - \frac{3 \sin^{2}{\left(\theta_{23} \right)}}{10} \\ X_{1/2}^{0,1/2;0,1/2}\left(\sigma_{2}\right) &:& 9 \\ X_{1/2}^{0,1/2;1,1/2}\left(\sigma_{2}\right) &:& 1 \\ X_{3/2}^{1,3/2;2,1/2}\left(\sigma_{2}\right) &:& 4 - 3 \sin^{2}{\left(\theta_{31} \right)} \\ X_{3/2}^{1,3/2;1,1/2}\left(\sigma_{3}\right) &:& 4 - 3 \sin^{2}{\left(\theta_{12} \right)} \\ X_{3/2}^{1,3/2;2,1/2}\left(\sigma_{3}\right) &:& 1 - \frac{3 \sin^{2}{\left(\theta_{12} \right)}}{4} \\ \end{eqnarray} $$

@mmikhasenko
Copy link
Collaborator

Often, it is needed to give not just diagonal terms, i.e. coefficients in front of |X|^2, but also interference terms,
in front of

X_i X_j*

It matches the Hessian with a factor 1/2. So, autodiff gives a general approach to obtain the matrix.
I have been very happy discovering this simple trick for other project, recently

https://github.com/mmikhasenko/PhotoAmbiguities.jl/blob/master/src/experiments.jl#L10-L19

@redeboer redeboer modified the milestones: 0.1.6, 0.1.7 Mar 17, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
📝 Docs Improvements or additions to documentation
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants