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

[Zehvogel]

Check duplicate issues.

• 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:

DD4hep/DDG4/python/DDSim/Helper/Physics.py

Line 21 in f3304b4

self._rejectPDGs = {1, 2, 3, 4, 5, 6, # quarks

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

f3304b4

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

DD4hep/DDG4/src/Geant4InputHandling.cpp

Lines 394 to 400 in f3304b4

	if ( status.isSet(G4PARTICLE_GEN_STABLE) ) {
	if ( prim.find(p->id) == prim.end() ) {
	G4PrimaryParticle* p4 = createG4Primary(p);
	prim[p->id] = p4;
	res.emplace_back(p,p4);
	}
	}

No rejections implemented here yet.

[andresailer]

Please see #1284.

[Zehvogel]

thx!