square_dalitz_plot

class SDPGenerator(m0, mi, legacy=True)

Bases: BaseGenerator

generate(N)

>>> from tf_pwa.generator.square_dalitz_plot import SDPGenerator
>>> gen = SDPGenerator(3.0, [1.0, 0.5, 0.1])
>>> p1, p2, p3 = gen.generate(100)

generate_SDP(m0, mi, N=1000, legacy=True)

generate square dalitz plot ditribution for 1,2

The legacy mode will include a cut off in the threshold.

(Source code, png, hires.png, pdf)

square_dalitz_cut(p)

Copy from EvtGen old version

\[|J| = 4 p q m_{12} \frac{\partial m_{12}}{\partial m'} \frac{\partial \cos\theta_{12}}{\partial \theta'}\]

\[\frac{\partial m_{12}}{\partial m'} = -\frac{\pi}{2} \sin (\pi m') (m_{12}^{max} - m_{12}^{min})\]

\[\frac{\partial \cos\theta_{12}}{\partial \theta'} = -\pi \sin (\pi \theta')\]

square_dalitz_variables(p)

Variables used of square dalitz plot, the first 2 is \(m'\) and \(\theta'\).

\[m' = \frac{1}{\pi}\cos^{-1} \left(2 \frac{m_{12}-m^{min}_{12}}{m^{max}_{12}-m^{min}_{12}} - 1\right)\]

\[\theta' = \frac{1}{\pi} \theta_{12}\]