`ddsim` error: corrupted double-linked list and infinite hang on sample with "stable" quarks

[Zehvogel]

Check duplicate issues.

• [X] Checked for duplicates

Goal

I am interested in the behavior of electrons in a certain process (e+e- to e nu qq). To save on simulation (and reconstruction) time I generated the process without running the parton shower, i.e. the two quarks show up in the file with generator status 1 (stable). My understanding was that the quarks would just not be simulated due to being in the rejection list here: https://github.com/AIDASoft/DD4hep/blob/f3304b47b3ef0450ebf2935ccb8ec1f7c398a59b/DDG4/python/DDSim/Helper/Physics.py#L21

When I try to run ddsim on the file it hangs after

LCIOFileReader   INFO  Created file reader. Try to open input data/gen/simple_whizard/enuqq.slcio
GenerationInit   INFO  +++ Initializing event 1. Within run:0 event 1.
LCIOFileReader   INFO  read collection MCParticle from event 1 in run 1 
LCIO4            INFO  +++ Particle interaction with 8 generator particles and 2 vertices ++++++++++++++++++++++++
PrimaryHandler   INFO  +++++ G4PrimaryVertex at (+0.00e+00,+0.00e+00,+0.00e+00) [mm] +0.00e+00 [ns]
PrimaryHandler   INFO  +++++ G4PrimaryVertex at (+0.00e+00,+0.00e+00,+0.00e+00) [mm] +0.00e+00 [ns]
ParticleHandler  INFO  +++ Event 0 Begin event action. Access event related information.
### Birks coefficients used in run time
corrupted double-linked list
 *** Break *** abort

when I run with edm4hep output the message is slightly different

LCIOFileReader   INFO  Created file reader. Try to open input data/gen/simple_whizard/enuqq.slcio
Geant4Output2EDM4hep INFO  Opened data/gen/simple_whizard/enuqq.10.SIM.edm4hep.root for output
GenerationInit   INFO  +++ Initializing event 1. Within run:0 event 1.
LCIOFileReader   INFO  read collection MCParticle from event 1 in run 1 
LCIO4            INFO  +++ Particle interaction with 8 generator particles and 2 vertices ++++++++++++++++++++++++
PrimaryHandler   INFO  +++++ G4PrimaryVertex at (+0.00e+00,+0.00e+00,+0.00e+00) [mm] +0.00e+00 [ns]
PrimaryHandler   INFO  +++++ G4PrimaryVertex at (+0.00e+00,+0.00e+00,+0.00e+00) [mm] +0.00e+00 [ns]
ParticleHandler  INFO  +++ Event 0 Begin event action. Access event related information.
### Birks coefficients used in run time
Geant4Output2EDM4hep INFO  +++ Saving EDM4hep event 1 run 1.
corrupted double-linked list
 *** Break *** abort

At first I suspected that my .slcio file is broken but checking it with dumpevent and anajob looks fine...

Operating System and Version

Alma9

compiler

gcc 11.4.1

ROOT Version

6.30/06

DD4hep Version

f3304b47b3ef0450ebf2935ccb8ec1f7c398a59b

Reproducer

wget https://cernbox.cern.ch/remote.php/dav/public-files/2LYvjagyMP6ytAl/enuqq.slcio
source /cvmfs/sw-nightlies.hsf.org/key4hep/setup.sh
ddsim --compactFile $K4GEO/FCCee/CLD/compact/CLD_o2_v06/CLD_o2_v06.xml \
      --outputFile enuqq.SIM.slcio \
      --inputFile enuqq.slcio \
      --numberOfEvents 10

Additional context

No response

[andresailer]

This goes via

https://github.com/AIDASoft/DD4hep/blob/f3304b47b3ef0450ebf2935ccb8ec1f7c398a59b/DDG4/src/Geant4InputHandling.cpp#L394-L400

No rejections implemented here yet.

[andresailer]

Please see #1284.

[Zehvogel]

thx!