acts: Bug? in ActsExampleFatrasTGeo event processing with Pythia input

I’ve spotted an issue of Fatras simulation when processing Pythia events.

I generated 50 pp 14 TeV events as:

./ActsExamplePythia8 -n 50 \                                                                                
--output-dir=data/pythia8/pp14_n50 \                                                
--output-csv \                                               
--gen-cms-energy-gev=14000

When doing their propagation with ActsExampleFatrasTGeo most events were processed within seconds or faster but then for a few events it took 15/10 mins to process and after that it gets stuck until I stop it manually:

10:34:47    Sequencer      INFO      finished event 36
10:34:49    Sequencer      INFO      finished event 42
10:34:51    Sequencer      INFO      finished event 15
10:34:52    Sequencer      INFO      finished event 43
10:34:53    Simulation     ERROR     Propagation reached the step count limit of 1000 (did 1000 steps)
10:34:54    FatrasSimula   ERROR     event 44 particle id=(62|0|74|0|0)|pdg=K+|q=1|m=0.49368|p=12.7496 failed to simulate with error PropagatorError:3: Propagation reached the configured maximum number of steps
10:34:54    Sequencer      INFO      finished event 16
10:34:55    Sequencer      INFO      finished event 44
10:34:59    Sequencer      INFO      finished event 45
10:35:03    Sequencer      INFO      finished event 46
10:35:06    Sequencer      INFO      finished event 47
10:35:10    Sequencer      INFO      finished event 48
10:35:14    Sequencer      INFO      finished event 49
10:35:17    Sequencer      INFO      finished event 21
10:35:20    Sequencer      INFO      finished event 22
10:35:24    Sequencer      INFO      finished event 23
10:35:28    Sequencer      INFO      finished event 24
10:50:41    Sequencer      INFO      finished event 13
10:59:40    Sequencer      INFO      finished event 14

Events 13 and 14 look similar to other in terms of size: Screenshot 2021-10-22 at 11 52 02

While processing events 13 or 14 my hard drive lost 10 GB of space. The issue was also reproduced with a similar dataset – when around 5 events are left to be processed it gets slow/stuck.

About this issue

  • Original URL
  • State: closed
  • Created 3 years ago
  • Comments: 15 (13 by maintainers)

Most upvoted comments

So after looking at the code in a bit more detail the only explanation possible is that something goes wrong with the simulation actor : #include “ActsFatras/Kernel/detail/SimulationActor.hpp” I will let you have a look, maybe you could pin point where the issue is happening.

To give you additional info : in Fatras/include/ActsFatras/Physics/ they are a few while loop. To check if the issue is coming from there or not you could just put a cout statement in each of those and see if you get stuck in one.