base

Basic amplitude model

class HelicityDecayCPV(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecay

decay model for CPV

get_g_ls(charge=1)

get_ls_amp(data, data_p, **kwargs)

init_params()

model_name = 'gls-cpv'

class HelicityDecayNP(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecay

Full helicity amplitude

\[A = H_{m_1, m_2} D_{m_0, m_1-m_2}^{J_0 *}(\varphi, \theta,0)\]

fit parameters is \(H_{m_1, m_2}\).

fix_unused_h()

get_H()

get_H_zero_mask()

get_factor()

get_factor_H(data=None, data_p=None, **kwargs)

get_helicity_amp(data=None, data_p=None, **kwargs)

get_ls_amp(data, data_p, **kwargs)

get_zero_index()

init_params()

model_name = 'helicity_full'

class HelicityDecayNPbf(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecayNP

get_H_barrier_factor(data, data_p, **kwargs)

get_helicity_amp(data, data_p, **kwargs)

get_ls_amp(data, data_p, **kwargs)

init_params()

model_name = 'helicity_full-bf'

class HelicityDecayP(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecayNP

\[H_{- m1, - m2} = P_0 P_1 P_2 (-1)^{J_1 + J_2 - J_0} H_{m1, m2}\]

get_helicity_amp(data, data_p, **kwargs)

init_params()

model_name = 'helicity_parity'

class HelicityDecayReduceH0(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecay

decay model that remove helicity =0 for massless particles

get_g_ls(charge=1)

get_helicity_list2()

init_params()

model_name = 'gls_reduce_h0'

class ParticleBW(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = \frac{1}{m_0^2 - m^2 - i m_0 \Gamma_0}\]

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

get_amp(data, _data_c=None, **kwargs)

get_num_var()

get_sympy_var()

model_name = 'BW'

class ParticleBWR2(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = \frac{1}{m_0^2 - m^2 - i m_0 \Gamma(m)}\]

The difference of BWR, BWR2 is the behavior when mass is below the threshold ( \(m_0 = 0.1 < 0.1 + 0.1 = m_1 + m_2\)).

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

get_amp(data, data_c, **kwargs)

model_name = 'BWR2'

class ParticleBWRBelowThreshold(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = \frac{1}{m_0^2 - m^2 - i m_0 \Gamma(m)}\]

get_amp(data, data_c, **kwargs)

model_name = 'BWR_below'

class ParticleBWRCoupling(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

Force \(q_0=1/d\) to avoid below theshold condition for BWR model, and remove other constant parts, then the \(\Gamma_0\) is coupling parameters.

\[R(m) = \frac{1}{m_0^2 - m^2 - i m_0 \Gamma_0 \frac{q}{m} q^{2l} B_L'^2(q, 1/d, d)}\]

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

get_amp(data, data_c, **kwargs)

get_sympy_dom(m, m0, g0, m1=None, m2=None, sheet=0)

model_name = 'BWR_coupling'

class ParticleBWR_normal(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = \frac{\sqrt{m_0 \Gamma(m)}}{m_0^2 - m^2 - i m_0 \Gamma(m)}\]

get_amp(data, data_c, **kwargs)

model_name = 'BWR_normal'

class ParticleDecay(*args, has_barrier_factor=True, l_list=None, barrier_factor_mass=False, has_ql=True, has_bprime=True, aligned=False, allow_cc=True, ls_list=None, barrier_factor_norm=False, params_polar=None, below_threshold=False, force_min_l=False, params_head=None, no_q0=False, helicity_inner_full=False, ls_selector=None, **kwargs)

Bases: HelicityDecay

get_ls_amp(data, data_p, **kwargs)

model_name = 'particle-decay'

class ParticleExp(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = e^{-|a| m}\]

get_amp(data, _data_c=None, **kwargs)

init_params()

model_name = 'exp'

class ParticleExpCom(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = e^{-(a+ib) m^2}\]

lineshape when \(a=1.0, b=10.\)

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

get_amp(data, _data_c=None, **kwargs)

init_params()

model_name = 'exp_com'

class ParticleGS(*args, **kwargs)

Bases: Particle

Gounaris G.J., Sakurai J.J., Phys. Rev. Lett., 21 (1968), pp. 244-247

c_daug2Mass: mass for daughter particle 2 (\(\pi^{+}\)) 0.13957039

c_daug3Mass: mass for daughter particle 3 (\(\pi^{0}\)) 0.1349768

\[R(m) = \frac{1 + D \Gamma_0 / m_0}{(m_0^2 -m^2) + f(m) - i m_0 \Gamma(m)}\]

\[f(m) = \Gamma_0 \frac{m_0 ^2 }{q_0^3} \left[q^2 [h(m)-h(m_0)] + (m_0^2 - m^2) q_0^2 \frac{d h}{d m}|_{m0} \right]\]

\[h(m) = \frac{2}{\pi} \frac{q}{m} \ln \left(\frac{m+2q}{2m_{\pi}} \right)\]

\[\frac{d h}{d m}|_{m0} = h(m_0) [(8q_0^2)^{-1} - (2m_0^2)^{-1}] + (2\pi m_0^2)^{-1}\]

\[D = \frac{f(0)}{\Gamma_0 m_0} = \frac{3}{\pi}\frac{m_\pi^2}{q_0^2} \ln \left(\frac{m_0 + 2q_0}{2 m_\pi }\right) + \frac{m_0}{2\pi q_0} - \frac{m_\pi^2 m_0}{\pi q_0^3}\]

get_amp(data, data_c, **kwargs)

model_name = 'GS_rho'

class ParticleKmatrix(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

get_amp(data, data_c=None, **kwargs)

get_beta(m, **kwargs)

init_params()

model_name = 'Kmatrix'

class ParticleLass(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

get_amp(data, data_c=None, **kwargs)

\[R(m) = \frac{m}{q cot \delta_B - i q} + e^{2i \delta_B}\frac{m_0 \Gamma_0 \frac{m_0}{q_0}} {(m_0^2 - m^2) - i m_0\Gamma_0 \frac{q}{m}\frac{m_0}{q_0}}\]

\[cot \delta_B = \frac{1}{a q} + \frac{1}{2} r q\]

\[e^{2i\delta_B} = \cos 2 \delta_B + i \sin 2\delta_B = \frac{cot^2\delta_B -1 }{cot^2 \delta_B +1} + i \frac{2 cot \delta_B }{cot^2 \delta_B +1 }\]

init_params()

model_name = 'LASS'

class ParticleOne(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = 1\]

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

get_amp(data, _data_c=None, **kwargs)

init_params()

model_name = 'one'

class ParticlePoly(*args, running_width=True, bw_l=None, width_norm=False, params_head=None, **kwargs)

Bases: Particle

\[R(m) = \sum c_i (m-m_0)^{n-i}\]

lineshape when \(c_0=1, c_1=c_2=0\)

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

get_amp(data, _data_c=None, **kwargs)

init_params()

model_name = 'poly'

get_parity_term(j1, p1, j2, p2, j3, p3)