jiangyi15 / tf-pwa

Partial Wave Analysis using TensorFlow.
MIT License
27 stars 13 forks source link

Problem with the TF-PWA angulare distributions #42

Open Leongrim opened 3 years ago

Leongrim commented 3 years ago

I continued the comparison and started to look into the channel D1->K+K-K0 with an a1 in the intermediate state.

There we have the problem that the model of TF-PWA the angulares distributions are modeled flat.

D(1)(2420)0_R_K~0K+_cos(beta) R_K~0K+_K~0_alpha (Black is data from Tensorwaves and rad the model of TF-PWA)

In our opinion the distribution shouldn't be flat. I used the following config file: config.txt

I hope you have an idea what is going on.

jiangyi15 commented 3 years ago

L46.

 spins:
 - 0

R_K~0K+: There is a typo.

This simple script can calculate the formula of angular distribution https://github.com/jiangyi15/tf-pwa/blob/dev/tutorials/examples/get_angle_distribution.py The projection of cos theta^cms can have such shape , but the phi^K0 is flat.

jiangyi15 commented 3 years ago

Here is the outputs

[D(1)(2420)0->a(1)(1260)++K-, a(1)(1260)+->K~0+K+] ls: ((0, 1.0), (1, 0.0)) lambda: (0, -0.0, -0.0, -0.0) : (sin(theta0)sin(theta1)cos(phi1) - cos(theta0)cos(theta1))2/3 ls: ((1, 1.0), (1, 0.0)) lambda: (0, -0.0, -0.0, -0.0) : sin(phi1)2sin(theta0)*2sin(theta1)2/2 ls: ((2, 1.0), (1, 0.0)) lambda: (0, -0.0, -0.0, -0.0) : (sin(theta0)sin(theta1)cos(phi1) + 2cos(theta0)cos(theta1))2/6 There is 3 possible LS components.

Leongrim commented 3 years ago

Thank you very much I will look into it in more detail.

jiangyi15 commented 3 years ago

@Leongrim . I find the problem of Tensorwaves, this is the formula from Tensorwaves

result = q.generate_transitions(        
     initial_state="D(1)(2420)0",        
     final_state=["K+","K-","K~0"],      
     allowed_intermediate_particles=["a(1)(1260)+"],
     formalism_type="canonical-helicity")

modol_builder = get_builder(result)
model_no_dynamics = modol_builder.generate()
sym.print_latex(model_no_dynamics.expression)

Screenshot_20210519_232545

The summation should be LS instead of helicity. As Screenshot_20210519_234509