Stalled YFS corrections after hadron decay
It seems that various setups get stuck when emitting soft photons after a hadron decay (@krauss was reporting the same issue to me a few days ago). Specifically I've seen this happen for two events when YFS calculates higher order ME corrections for a Scalar_To_Scalar_Lepton_Neutrino
, see below.
It can be reproduced on a recent master (I've used 28499e91) with the run card in Examples/CI/LO_Z
, e.g. like this:
Sherpa -e0
Sherpa BATCH_MODE=4 EVENT_DISPLAY_INTERVAL=1 -e10k -R 4 -l stuck.log -o15 EVT_OUTPUT_START=1350
(if that doesn't do it, you should be able to find your own event getting stuck by generating 100k events or so, it seems to happen roughly every 5 to 10 thousand events)
Here is the output for the stuck event: stuck1.log The event does not produce any more output and event generation is stalled.
And another example: stuck2.log
Now, when I revert the recent YFS-related commit f1b4e3f5, this problem seems to disappear, all runs conclude (or crash in another unrelated way, but that's for another day