Crash when calculating normalisation for model with two nearby resonances
Bug in LauIsobarDynamics::calcDPNormalisation(), around line 829 in
https://laura.hepforge.org/doc/doxygen/v2r1/LauIsobarDynamics_8cc_source.html#l00829
When you have two resonances, one close to threshold, you assume that res1Max is smaller than res2Min, while if the two resonances are close together, this may not be true. I had this problem running with Lc(2880) and Lc(2940) in Dp spectrum. I had to add another "if" which checks that if res1Max>res2Min, one just needs to divide the DP into two regions, (m13min, res2Max) and (res2Max, m13max). I suspect similar thing can affect other cases.
Change History (6)
Description: |
modified (diff)
|
Owner: |
set to Thomas Latham
|
Status: |
new →
assigned
|
Status: |
assigned →
testing
|
Resolution: |
→ fixed
|
Status: |
testing →
closed
|
(In [119]) Tidy up code for determining the integration scheme once and using for all reclaculations thereafter, addresses #11. Also fix (in this v3 branch) the bug in the determination of the scheme, addresses #35.