celeritas icon indicating copy to clipboard operation
celeritas copied to clipboard

Implement Celeritas-DD4hep integration plugin

Open rahmans1 opened this issue 8 months ago • 14 comments

See issue https://github.com/celeritas-project/celeritas/issues/1731

rahmans1 avatar Apr 21 '25 18:04 rahmans1

Test summary

 5 685 files   9 141 suites   7m 6s :stopwatch:  1 681 tests  1 673 :white_check_mark:   8 :zzz: 0 :x: 29 493 runs  29 345 :white_check_mark: 148 :zzz: 0 :x:

Results for commit cf3684d3.

:recycle: This comment has been updated with latest results.

github-actions[bot] avatar Apr 21 '25 18:04 github-actions[bot]

@sethrj This implementation is pretty preliminary based on streamlining the multiple test beds I was working on. I decided to use the OpenDataDetector from acts-project for validation as it should be a simpler approximation to something like the EIC ePIC geometry. Build and run instructions are inside the steeringFile.py. I can see celeritas logs but haven't validated the output yet.

/home/srahman1/celeritas-dd4hep/celeritas/src/corecel/sys/ScopedMpiInit.cc:57: debug: MPI initialization took 0.0408042s
[email protected]:86: Activating tracking manager
[email protected]:94: Preparing physics table for e+
[email protected]:94: Preparing physics table for e-
[email protected]:94: Preparing physics table for gamma
[email protected]:54: Building physics table for gamma
[email protected]:54: Building physics table for e-
[email protected]:54: Building physics table for e+
[email protected]:105: Initializing shared params
[email protected]:100: Disabling GPU support since CUDA and HIP are disabled
[email protected]:31: Activating Celeritas version 0.0.0+8fa40ce0 on CPU
[email protected]:1204: Transferring data from Geant4
[email protected]:375: Loaded 3 particles
[email protected]:434: Loaded 40 isotopes
[email protected]:487: Loaded 14 elements
[email protected]:661: Loaded 16 geo materials
[email protected]:765: Loaded 16 physics materials
[email protected]:345: Saving MSC model 'UrbanMsc' for particle e+ (-11)
[email protected]:345: Saving MSC model 'WentzelVIUni' for particle e+ (-11)
[email protected]:119: Saving process 'eIoni' for particle e+ (-11)
[email protected]:119: Saving process 'eBrem' for particle e+ (-11)
[email protected]:119: Saving process 'annihil' for particle e+ (-11)
[email protected]:119: Saving process 'CoulombScat' for particle e+ (-11)
[email protected]:345: Saving MSC model 'UrbanMsc' for particle e- (11)
[email protected]:345: Saving MSC model 'WentzelVIUni' for particle e- (11)
[email protected]:119: Saving process 'eIoni' for particle e- (11)
[email protected]:119: Saving process 'eBrem' for particle e- (11)
[email protected]:119: Saving process 'CoulombScat' for particle e- (11)
[email protected]:119: Saving process 'phot' for particle gamma (22)
[email protected]:119: Saving process 'compt' for particle gamma (22)
[email protected]:119: Saving process 'conv' for particle gamma (22)
[email protected]:119: Saving process 'Rayl' for particle gamma (22)
[email protected]:968: Loaded 11 processes
[email protected]:969: Loaded 0 optical models
[email protected]:798: Loaded 6 regions
[email protected]:1347: Loaded 1677 of 1677 volumes
[email protected]:104: ... 0.0116159 s
[email protected]:1272: Loading external elemental data
[email protected]:104: ... 0.0358131 s
[email protected]:158: No magnetic field
[email protected]:383: Initializing problem
status: Converting Geant4 geometry
[email protected]:104: ... 0.0855203 s
[email protected]:461: Initializing tracking information
[email protected]:268: Omitting coulomb_scat from physics process list
[email protected]:150: Ignored process class coulomb_scat
[email protected]:256: Set default track order none
[email protected]:368: Celeritas core setup complete
[email protected]:32: Mapped sensitive detector "PixelEndcapN" on logical volume "component0Silicon"@0x652986008ae0 (ID=0) to VecGeom volume "component0Silicon@0" (ID=0)
[email protected]:32: Mapped sensitive detector "PixelEndcapN" on logical volume "component0Silicon"@0x65298600a800 (ID=3) to VecGeom volume "component0Silicon@1" (ID=3)
[email protected]:32: Mapped sensitive detector "PixelEndcapP" on logical volume "component2Silicon"@0x652986014ff0 (ID=8) to VecGeom volume "component2Silicon@0" (ID=8)
[email protected]:32: Mapped sensitive detector "PixelEndcapP" on logical volume "component2Silicon"@0x6529860154c0 (ID=11) to VecGeom volume "component2Silicon@1" (ID=11)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x652986015660 (ID=12) to VecGeom volume "component0Silicon@2" (ID=12)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x652986015a80 (ID=15) to VecGeom volume "component0Silicon@3" (ID=15)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x652986015ed0 (ID=18) to VecGeom volume "component0Silicon@4" (ID=18)
[email protected]:32: Mapped sensitive detector "ShortStripBarrel" on logical volume "component3"@0x652986017f90 (ID=25) to VecGeom volume "component3@0" (ID=29)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x652986018180 (ID=26) to VecGeom volume "component0Silicon@5" (ID=30)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x6529860185d0 (ID=29) to VecGeom volume "component0Silicon@6" (ID=33)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x652986018a20 (ID=32) to VecGeom volume "component0Silicon@7" (ID=36)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component0Silicon"@0x652986018e70 (ID=35) to VecGeom volume "component0Silicon@8" (ID=39)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component4Silicon"@0x65298601a1d0 (ID=40) to VecGeom volume "component4Silicon@0" (ID=48)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component0Silicon"@0x65298601a3f0 (ID=41) to VecGeom volume "component0Silicon@9" (ID=49)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component4Silicon"@0x65298601ab50 (ID=46) to VecGeom volume "component4Silicon@1" (ID=58)
[email protected]:32: Mapped sensitive detector "LongStripBarrel" on logical volume "component0"@0x65298601acd0 (ID=47) to VecGeom volume "component0@2" (ID=59)
[email protected]:32: Mapped sensitive detector "LongStripBarrel" on logical volume "component9"@0x65298601be40 (ID=58) to VecGeom volume "component9" (ID=78)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component0Silicon"@0x65298601e690 (ID=80) to VecGeom volume "component0Silicon@10" (ID=100)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component4Silicon"@0x65298601edf0 (ID=85) to VecGeom volume "component4Silicon@2" (ID=109)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component0Silicon"@0x65298601ef80 (ID=86) to VecGeom volume "component0Silicon@11" (ID=110)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component4Silicon"@0x65298601f6b0 (ID=91) to VecGeom volume "component4Silicon@3" (ID=119)
[email protected]:32: Mapped sensitive detector "PixelBarrel" on logical volume "component0"@0x65298601f830 (ID=92) to VecGeom volume "component0@0" (ID=120)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986029c10 (ID=150) to VecGeom volume "slice4@48" (ID=178)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602b850 (ID=158) to VecGeom volume "slice4@49" (ID=186)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602c3d0 (ID=166) to VecGeom volume "slice4@50" (ID=194)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602ce40 (ID=174) to VecGeom volume "slice4@51" (ID=202)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602d7e0 (ID=182) to VecGeom volume "slice4@52" (ID=210)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602e370 (ID=190) to VecGeom volume "slice4@53" (ID=218)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602ed70 (ID=198) to VecGeom volume "slice4@54" (ID=226)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298602fd70 (ID=206) to VecGeom volume "slice4@55" (ID=234)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986030890 (ID=214) to VecGeom volume "slice4@56" (ID=242)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986031560 (ID=222) to VecGeom volume "slice4@57" (ID=250)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986031f60 (ID=230) to VecGeom volume "slice4@58" (ID=258)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986032960 (ID=238) to VecGeom volume "slice4@59" (ID=266)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986033360 (ID=246) to VecGeom volume "slice4@60" (ID=274)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986033d60 (ID=254) to VecGeom volume "slice4@61" (ID=282)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986034760 (ID=262) to VecGeom volume "slice4@62" (ID=290)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986035160 (ID=270) to VecGeom volume "slice4@63" (ID=298)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986036550 (ID=278) to VecGeom volume "slice4@64" (ID=306)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986037420 (ID=286) to VecGeom volume "slice4@65" (ID=314)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986037e20 (ID=294) to VecGeom volume "slice4@66" (ID=322)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986038820 (ID=302) to VecGeom volume "slice4@67" (ID=330)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986039220 (ID=310) to VecGeom volume "slice4@68" (ID=338)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986039c20 (ID=318) to VecGeom volume "slice4@69" (ID=346)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603a620 (ID=326) to VecGeom volume "slice4@70" (ID=354)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603b020 (ID=334) to VecGeom volume "slice4@71" (ID=362)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603ba20 (ID=342) to VecGeom volume "slice4@72" (ID=370)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603c420 (ID=350) to VecGeom volume "slice4@73" (ID=378)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603ce20 (ID=358) to VecGeom volume "slice4@74" (ID=386)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603d820 (ID=366) to VecGeom volume "slice4@75" (ID=394)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603e220 (ID=374) to VecGeom volume "slice4@76" (ID=402)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603ec20 (ID=382) to VecGeom volume "slice4@77" (ID=410)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298603f620 (ID=390) to VecGeom volume "slice4@78" (ID=418)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986040020 (ID=398) to VecGeom volume "slice4@79" (ID=426)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986035f80 (ID=406) to VecGeom volume "slice4@80" (ID=434)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986043610 (ID=414) to VecGeom volume "slice4@81" (ID=442)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986044010 (ID=422) to VecGeom volume "slice4@82" (ID=450)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986044a10 (ID=430) to VecGeom volume "slice4@83" (ID=458)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986045410 (ID=438) to VecGeom volume "slice4@84" (ID=466)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986045e10 (ID=446) to VecGeom volume "slice4@85" (ID=474)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986046810 (ID=454) to VecGeom volume "slice4@86" (ID=482)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986047210 (ID=462) to VecGeom volume "slice4@87" (ID=490)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986047c10 (ID=470) to VecGeom volume "slice4@88" (ID=498)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986048610 (ID=478) to VecGeom volume "slice4@89" (ID=506)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986049010 (ID=486) to VecGeom volume "slice4@90" (ID=514)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x652986049a10 (ID=494) to VecGeom volume "slice4@91" (ID=522)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298604a410 (ID=502) to VecGeom volume "slice4@92" (ID=530)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298604ae10 (ID=510) to VecGeom volume "slice4@93" (ID=538)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298600f090 (ID=518) to VecGeom volume "slice4@94" (ID=546)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x65298600fa90 (ID=526) to VecGeom volume "slice4@95" (ID=554)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986011760 (ID=533) to VecGeom volume "slice3@132" (ID=561)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986012190 (ID=537) to VecGeom volume "slice3@133" (ID=565)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986012670 (ID=541) to VecGeom volume "slice3@134" (ID=569)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986012b90 (ID=545) to VecGeom volume "slice3@135" (ID=573)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x6529860130b0 (ID=549) to VecGeom volume "slice3@136" (ID=577)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x6529860135d0 (ID=553) to VecGeom volume "slice3@137" (ID=581)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986013af0 (ID=557) to VecGeom volume "slice3@138" (ID=585)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986014010 (ID=561) to VecGeom volume "slice3@139" (ID=589)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986014530 (ID=565) to VecGeom volume "slice3@140" (ID=593)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x6529860571d0 (ID=569) to VecGeom volume "slice3@141" (ID=597)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x6529860576f0 (ID=573) to VecGeom volume "slice3@142" (ID=601)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986057c10 (ID=577) to VecGeom volume "slice3@143" (ID=605)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986058130 (ID=581) to VecGeom volume "slice3@144" (ID=609)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986058650 (ID=585) to VecGeom volume "slice3@145" (ID=613)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986058b70 (ID=589) to VecGeom volume "slice3@146" (ID=617)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986059090 (ID=593) to VecGeom volume "slice3@147" (ID=621)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x6529860599c0 (ID=597) to VecGeom volume "slice3@148" (ID=625)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605a700 (ID=601) to VecGeom volume "slice3@149" (ID=629)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605ac20 (ID=605) to VecGeom volume "slice3@150" (ID=633)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605b140 (ID=609) to VecGeom volume "slice3@151" (ID=637)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605b660 (ID=613) to VecGeom volume "slice3@152" (ID=641)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605bb80 (ID=617) to VecGeom volume "slice3@153" (ID=645)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605c0a0 (ID=621) to VecGeom volume "slice3@154" (ID=649)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605c5c0 (ID=625) to VecGeom volume "slice3@155" (ID=653)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605cae0 (ID=629) to VecGeom volume "slice3@156" (ID=657)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605d000 (ID=633) to VecGeom volume "slice3@157" (ID=661)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605d520 (ID=637) to VecGeom volume "slice3@158" (ID=665)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605da40 (ID=641) to VecGeom volume "slice3@159" (ID=669)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605df60 (ID=645) to VecGeom volume "slice3@160" (ID=673)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605e480 (ID=649) to VecGeom volume "slice3@161" (ID=677)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605e9a0 (ID=653) to VecGeom volume "slice3@162" (ID=681)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605eec0 (ID=657) to VecGeom volume "slice3@163" (ID=685)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x65298605fbf0 (ID=661) to VecGeom volume "slice3@164" (ID=689)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986060110 (ID=665) to VecGeom volume "slice3@165" (ID=693)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986060630 (ID=669) to VecGeom volume "slice3@166" (ID=697)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x652986060b50 (ID=673) to VecGeom volume "slice3@167" (ID=701)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986063a10 (ID=704) to VecGeom volume "slice4@0" (ID=736)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860643d0 (ID=712) to VecGeom volume "slice4@1" (ID=744)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986064dd0 (ID=720) to VecGeom volume "slice4@2" (ID=752)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860657d0 (ID=728) to VecGeom volume "slice4@3" (ID=760)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860661d0 (ID=736) to VecGeom volume "slice4@4" (ID=768)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986066bd0 (ID=744) to VecGeom volume "slice4@5" (ID=776)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860675d0 (ID=752) to VecGeom volume "slice4@6" (ID=784)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986067fd0 (ID=760) to VecGeom volume "slice4@7" (ID=792)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860689d0 (ID=768) to VecGeom volume "slice4@8" (ID=800)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860693d0 (ID=776) to VecGeom volume "slice4@9" (ID=808)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986069dd0 (ID=784) to VecGeom volume "slice4@10" (ID=816)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298606a7d0 (ID=792) to VecGeom volume "slice4@11" (ID=824)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298606b1d0 (ID=800) to VecGeom volume "slice4@12" (ID=832)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986040ca0 (ID=808) to VecGeom volume "slice4@13" (ID=840)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860416a0 (ID=816) to VecGeom volume "slice4@14" (ID=848)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298606e8e0 (ID=824) to VecGeom volume "slice4@15" (ID=856)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298606faf0 (ID=832) to VecGeom volume "slice4@16" (ID=864)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986070d10 (ID=840) to VecGeom volume "slice4@17" (ID=872)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986071710 (ID=848) to VecGeom volume "slice4@18" (ID=880)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986072110 (ID=856) to VecGeom volume "slice4@19" (ID=888)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986072b10 (ID=864) to VecGeom volume "slice4@20" (ID=896)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986073510 (ID=872) to VecGeom volume "slice4@21" (ID=904)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986073f10 (ID=880) to VecGeom volume "slice4@22" (ID=912)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986074910 (ID=888) to VecGeom volume "slice4@23" (ID=920)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986075310 (ID=896) to VecGeom volume "slice4@24" (ID=928)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986075d10 (ID=904) to VecGeom volume "slice4@25" (ID=936)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986076710 (ID=912) to VecGeom volume "slice4@26" (ID=944)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986077110 (ID=920) to VecGeom volume "slice4@27" (ID=952)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986077b10 (ID=928) to VecGeom volume "slice4@28" (ID=960)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986078510 (ID=936) to VecGeom volume "slice4@29" (ID=968)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986078f10 (ID=944) to VecGeom volume "slice4@30" (ID=976)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986079910 (ID=952) to VecGeom volume "slice4@31" (ID=984)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298606f2e0 (ID=960) to VecGeom volume "slice4@32" (ID=992)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607b520 (ID=968) to VecGeom volume "slice4@33" (ID=1000)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607bf20 (ID=976) to VecGeom volume "slice4@34" (ID=1008)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607c920 (ID=984) to VecGeom volume "slice4@35" (ID=1016)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607d320 (ID=992) to VecGeom volume "slice4@36" (ID=1024)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607dd20 (ID=1000) to VecGeom volume "slice4@37" (ID=1032)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607e720 (ID=1008) to VecGeom volume "slice4@38" (ID=1040)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607f120 (ID=1016) to VecGeom volume "slice4@39" (ID=1048)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298607fb20 (ID=1024) to VecGeom volume "slice4@40" (ID=1056)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298604ba90 (ID=1032) to VecGeom volume "slice4@41" (ID=1064)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298604c490 (ID=1040) to VecGeom volume "slice4@42" (ID=1072)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298604ce90 (ID=1048) to VecGeom volume "slice4@43" (ID=1080)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x65298604f0c0 (ID=1056) to VecGeom volume "slice4@44" (ID=1088)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860502d0 (ID=1064) to VecGeom volume "slice4@45" (ID=1096)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x652986050cd0 (ID=1072) to VecGeom volume "slice4@46" (ID=1104)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x6529860516d0 (ID=1080) to VecGeom volume "slice4@47" (ID=1112)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986092260 (ID=1135) to VecGeom volume "slice3@96" (ID=1167)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986092740 (ID=1139) to VecGeom volume "slice3@97" (ID=1171)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986092c60 (ID=1143) to VecGeom volume "slice3@98" (ID=1175)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986093180 (ID=1147) to VecGeom volume "slice3@99" (ID=1179)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860936a0 (ID=1151) to VecGeom volume "slice3@100" (ID=1183)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986093bc0 (ID=1155) to VecGeom volume "slice3@101" (ID=1187)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860940e0 (ID=1159) to VecGeom volume "slice3@102" (ID=1191)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986094600 (ID=1163) to VecGeom volume "slice3@103" (ID=1195)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986094b20 (ID=1167) to VecGeom volume "slice3@104" (ID=1199)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986095040 (ID=1171) to VecGeom volume "slice3@105" (ID=1203)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986095560 (ID=1175) to VecGeom volume "slice3@106" (ID=1207)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986095a80 (ID=1179) to VecGeom volume "slice3@107" (ID=1211)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986095fa0 (ID=1183) to VecGeom volume "slice3@108" (ID=1215)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860964c0 (ID=1187) to VecGeom volume "slice3@109" (ID=1219)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860969e0 (ID=1191) to VecGeom volume "slice3@110" (ID=1223)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986096f00 (ID=1195) to VecGeom volume "slice3@111" (ID=1227)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986097420 (ID=1199) to VecGeom volume "slice3@112" (ID=1231)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986097940 (ID=1203) to VecGeom volume "slice3@113" (ID=1235)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986097e60 (ID=1207) to VecGeom volume "slice3@114" (ID=1239)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986098380 (ID=1211) to VecGeom volume "slice3@115" (ID=1243)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860988a0 (ID=1215) to VecGeom volume "slice3@116" (ID=1247)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986098dc0 (ID=1219) to VecGeom volume "slice3@117" (ID=1251)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x6529860992e0 (ID=1223) to VecGeom volume "slice3@118" (ID=1255)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986099800 (ID=1227) to VecGeom volume "slice3@119" (ID=1259)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x652986099d20 (ID=1231) to VecGeom volume "slice3@120" (ID=1263)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609a240 (ID=1235) to VecGeom volume "slice3@121" (ID=1267)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609a760 (ID=1239) to VecGeom volume "slice3@122" (ID=1271)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609ac80 (ID=1243) to VecGeom volume "slice3@123" (ID=1275)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609b1a0 (ID=1247) to VecGeom volume "slice3@124" (ID=1279)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609b6c0 (ID=1251) to VecGeom volume "slice3@125" (ID=1283)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609bbe0 (ID=1255) to VecGeom volume "slice3@126" (ID=1287)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609c100 (ID=1259) to VecGeom volume "slice3@127" (ID=1291)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609ce30 (ID=1263) to VecGeom volume "slice3@128" (ID=1295)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609d350 (ID=1267) to VecGeom volume "slice3@129" (ID=1299)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609d870 (ID=1271) to VecGeom volume "slice3@130" (ID=1303)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x65298609dd90 (ID=1275) to VecGeom volume "slice3@131" (ID=1307)
[email protected]:468: Setting number of streams to 1
[email protected]:491: Wrote Geant4 diagnostic output to "trackingmanager-offload.out.json"
[email protected]:104: ... 0.355062 s
[email protected]:173: Constructing local state
[email protected]:105: Setting up thread-local hit processor for 190 sensitive detectors
[email protected]:80: Celeritas core state initialization complete
[email protected]:157: Verifying tracking manager
[email protected]:397: Reconstituted 4600543818244419859 hits for event 0
[email protected]:49: Finalizing Celeritas
[email protected]:208: Destroying local state
[email protected]:418: Finalizing Celeritas after 1936 steps from 871 offloaded tracks over 10 events, generating 4600543818244419859 hits
[email protected]:221: Deallocating hit processor
[email protected]:229: Finalizing Celeritas
[email protected]:491: Wrote Geant4 diagnostic output to "trackingmanager-offload.out.json"
[email protected]:367: Resetting shared parameters
[email protected]:94: Deallocating host core state (stream 0)
[email protected]:320: Clearing VecGeom CPU data

rahmans1 avatar May 29 '25 06:05 rahmans1

Running with gdb to investigate segfaults gives:

ParticleHandler  INFO  +++ Event 0 Begin event action. Access event related information.
### Birks coefficients used in run time

Thread 1 "python" received signal SIGSEGV, Segmentation fault.
dd4hep::sim::Geant4StepHandler::birkAttenuation (this=0x7fffffffaf40)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4StepHandler.cpp:89
89        const G4MaterialCutsCouple* couple = trk->GetMaterialCutsCouple();
(gdb) bt
#0  dd4hep::sim::Geant4StepHandler::birkAttenuation (this=0x7fffffffaf40)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4StepHandler.cpp:89
#1  0x00007fffe0383361 in dd4hep::sim::Geant4StepHandler::totalEnergy (this=0x7fffffffaf40)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/include/DDG4/Geant4StepHandler.h:79
#2  dd4hep::sim::Geant4HitData::extractContribution (step=0x55557e37f250, ApplyBirksLaw=true)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Data.cpp:82
#3  0x00007fffe04cf6ef in dd4hep::sim::Geant4SensitiveAction<dd4hep::sim::Geant4ScintillatorCalorimeter>::process (this=0x555569f0ac50, step=0x55557e37f250)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/plugins/Geant4SDActions.cpp:461
#4  0x00007fffe03cf402 in dd4hep::sim::Geant4SensDetActionSequence::process (this=0x55556a009930,
    step=0x55557e37f250, history=0x0)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4SensDetAction.cpp:449
#5  0x00007fffd3a958b5 in celeritas::detail::HitProcessor::operator()(celeritas::DetectorStepOutput const&, unsigned long) const () from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libceleritas.so
#6  0x00007fffd3a9532e in celeritas::detail::HitProcessor::operator()(celeritas::StepStateData<(celeritas::Ownership)1, (celeritas::MemSpace)0> const&) ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libceleritas.so
#7  0x00007fffd3a78050 in celeritas::detail::StepGatherAction<(celeritas::StepPoint)1>::step(celeritas::CoreParams const&, celeritas::CoreState<(celeritas::MemSpace)0>&) const ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libceleritas.so
#8  0x00007fffd396bf5a in void celeritas::ActionSequence::step<(celeritas::MemSpace)0>(celeritas::CoreParams const&, celeritas::CoreState<(celeritas::MemSpace)0>&) ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libceleritas.so
#9  0x00007fffd397fbdf in celeritas::Stepper<(celeritas::MemSpace)0>::operator()() ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libceleritas.so
#10 0x00007fffd3c5d703 in celeritas::LocalTransporter::Flush() ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libaccel.so
#11 0x00007fffd3c76508 in celeritas::TrackingManager::FlushEvent() ()
   from /home/srahman1/celeritas-dd4hep/celeritas/install/lib/libaccel.so
#12 0x00007fffe13f558a in G4EventManager::DoProcessing(G4Event*) () from /opt/local/lib/libG4event.so
#13 0x00007fffdbb7396f in G4RunManager::DoEventLoop(int, char const*, int) () from /opt/local/lib/libG4run.so
#14 0x00007fffdbb70a7e in G4RunManager::BeamOn(int, char const*, int) () from /opt/local/lib/libG4run.so
#15 0x00007fffe03db978 in dd4hep::sim::Geant4UIManager::start (this=this@entry=0x555569800f90)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4UIManager.cpp:253
--Type <RET> for more, q to quit, c to continue without paging--
#16 0x00007fffe03dbafb in dd4hep::sim::Geant4UIManager::operator() (this=0x555569800f90)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4UIManager.cpp:176
#17 0x00007fffe03dbafb in virtual thunk to dd4hep::sim::Geant4UIManager::operator()(void*) ()
   from /home/srahman1/celeritas-dd4hep/DD4hep/install/lib/libDDG4.so.1.32
#18 0x00007fffe0387fb2 in dd4hep::sim::Geant4Exec::run (kernel=...)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Exec.cpp:645
#19 0x00007fffe03b2f4a in dd4hep::sim::Geant4Kernel::run (
    this=0x7fffe041f248 <dd4hep::sim::Geant4StepHandler::birkAttenuation() const::s_emSaturation>)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Kernel.cpp:409

rahmans1 avatar Jun 20 '25 14:06 rahmans1

Building Geant4

cmake -DCMAKE_BUILD_TYPE=Debug -DCMAKE_CXX_STANDARD=20 -DCMAKE_INSTALL_PREFIX=install -DGEANT4_USE_GDML=ON -DGEANT4_BUILD_TLS_MODEL=global-dynamic -DGEANT4_BUILD_MULTITHREADED=ON -S . -B build
cmake --build build -j 8
cmake --install build
source $PWD/install/bin/geant4.sh
export LD_LIBRARY_PATH=$PWD/install/lib:${LD_LIBRARY_PATH}

Building DD4hep

cmake -DGeant4_DIR=$PWD/../geant4/install/lib/cmake/Geant4 -DCMAKE_CXX_STANDARD=20 -DDD4HEP_USE_GEANT4=ON -DDD4HEP_USE_HEPMC3=ON -DDD4HEP_USE_EDM4HEP=ON -DCMAKE_BUILD_TYPE=Debug -DCMAKE_INSTALL_PREFIX=install -S . -B build
cmake --build build -j 8
cmake --install build
source $PWD/install/bin/thisdd4hep.sh
export LD_LIBRARY_PATH=$PWD/install/lib:${LD_LIBRARY_PATH}

Building celeritas

cmake -DCMAKE_BUILD_TYPE=Debug -DCELERITAS_USE_DD4hep=ON -DCMAKE_INSTALL_PREFIX=install -S . -B build
cmake --build build -j 8
cmake --install build
export LD_LIBRARY_PATH=$PWD/install/lib:${LD_LIBRARY_PATH}

Running with gdb

CELER_LOG=debug CELER_LOG_LOCAL=debug gdb --args python ~/celeritas-dd4hep/DD4hep/install/bin/ddsim --action.run DDcelerRunAction --compactFile=$PWD/install-debug/share/OpenDataDetector/xml/OpenDataDetector.xml --steering steeringFile.py -N 10 -G --gun.particle e- --gun.etaMax 4 --gun.etaMin 1 --gun.distribution uniform
(gdb) run
Starting program: /usr/bin/python /home/srahman1/celeritas-dd4hep/DD4hep/install/bin/ddsim --action.run DDcelerRunAction --compactFile=/home/srahman1/celeritas-dd4hep/OpenDataDetector/install/share/OpenDataDetector/xml/OpenDataDetector.xml --steering steeringFile.py -N 10 -G --gun.particle e- --gun.etaMax 4 --gun.etaMin 1 --gun.distribution uniform
[Thread debugging using libthread_db enabled]
Using host libthread_db library "/lib/x86_64-linux-gnu/libthread_db.so.1".
[Detaching after vfork from child process 247698]
[Detaching after vfork from child process 247702]
[Detaching after vfork from child process 247704]
[New Thread 0x7fffe26166c0 (LWP 247715)]
PersistencyIO    INFO  +++ Set Streamer to dd4hep::OpaqueDataBlock
Info in <TGeoManager::TGeoManager>: Geometry default, Detector Geometry created
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
DDG4             INFO +++ Imported 0 global values to namespace:DDG4
Geant4Kernel           OutputLevel:  2
Geant4Kernel           UI:           UI
Geant4Kernel           NumEvents:    10
Geant4Kernel           NumThreads:   0
DDG4             INFO +++  List of sensitive detectors:
Geant4UI         INFO  +++ PhysicsList> Install Geant4 control directory:/ddg4/PhysicsList/
PhysicsList            +++ Dump of physics list component(s)
PhysicsList            +++ Extension name       FTFP_BERT
PhysicsList            +++ Transportation flag: 0
PhysicsList            +++ Program decays:      1
PhysicsList            +++ RangeCut:            0.700000
PhysicsList            +++ Verbosity:           1
PhysicsList            +++ Dump of physics list component(s)
PhysicsList            +++ Extension name       FTFP_BERT
PhysicsList            +++ Transportation flag: 0
PhysicsList            +++ Program decays:      1
PhysicsList            +++ RangeCut:            0.700000
PhysicsList            +++ Verbosity:           1
XMLLoader        INFO  +++ Processing XML file: file:/home/srahman1/celeritas-dd4hep/OpenDataDetector/install/share/OpenDataDetector/xml/OpenDataDetector.xml
..........................................................................................................................................................................................

**************************************************************
 Geant4 version Name: geant4-11-02-patch-02 [MT]   (21-June-2024)
                       Copyright : Geant4 Collaboration
                      References : NIM A 506 (2003), 250-303
                                 : IEEE-TNS 53 (2006), 270-278
                                 : NIM A 835 (2016), 186-225
                             WWW : http://geant4.org/
**************************************************************

Geant4UI         INFO  +++ Geant4RunManager> Install Geant4 control directory:/ddg4/Geant4RunManager/
Geant4RunManager WARN  +++ Configured run manager of type: G4RunManager.
Geant4Kernel     WARN  +++ Multi-threaded mode requested, but not supported by this compilation of Geant4.
Geant4Kernel     WARN  +++ Falling back to single threaded mode.
G4PhysListFactory::GetReferencePhysList <FTFP_BERT>  EMoption= 0
<<< Geant4 Physics List simulation engine: FTFP_BERT


DDcelerTMI       INFO  Using FTFP_BERT physics list with Celeritas tracking for e-/e+/gamma.
[New Thread 0x7fffd20a96c0 (LWP 247718)]
/home/srahman1/celeritas-dd4hep/celeritas/src/corecel/sys/ScopedMpiInit.cc:57: debug: MPI initialization took 0.0793693s
DDcelerTMI       INFO  Celeritas TrackingManager registered.
Myphysics        INFO  Registered Geant4 physics constructor G4StepLimiterPhysics to physics list
UserInitialization INFO  +++ Executing Geant4UserActionInitialization::Build. Context:0xccac7b0 Kernel:0xcc481b0 [-1]
SignalHandler    INFO  ++ Re-apply signal handlers
Geant4Converter  INFO  +++  Successfully converted geometry to Geant4. [  0.370 seconds]
Geant4UI         INFO  +++ ConstructGeo> Install Geant4 control directory:/ddg4/ConstructGeo/

 hInelastic FTFP_BERT : threshold between BERT and FTFP is over the interval
 for pions :   3 to 6 GeV
 for kaons :   3 to 6 GeV
 for proton :  3 to 6 GeV
 for neutron : 3 to 6 GeV

### Adding tracking cuts for neutron  TimeCut(ns)= 10000  KinEnergyCut(MeV)= 0
[email protected]:86: Activating tracking manager
GlobalRangeCut   INFO  +++ Set PhysicsList RangeCut to 0.700000 mm
UI               INFO  ++ Start run with 10 events.
[email protected]:94: Preparing physics table for e+
[email protected]:94: Preparing physics table for e-
[email protected]:94: Preparing physics table for gamma
[email protected]:54: Building physics table for gamma
=======================================================================
======                 Electromagnetic Physics Parameters      ========
=======================================================================
LPM effect enabled                                 1
Enable creation and use of sampling tables         0
Apply cuts on all EM processes                     0
Use combined TransportationWithMsc                 Disabled
Use general process                                1
Enable linear polarisation for gamma               0
Enable photoeffect sampling below K-shell          1
Enable sampling of quantum entanglement            0
X-section factor for integral approach             0.8
Min kinetic energy for tables                      100 eV
Max kinetic energy for tables                      100 TeV
Number of bins per decade of a table               7
Verbose level                                      1
Verbose level for worker thread                    0
Bremsstrahlung energy threshold above which
  primary e+- is added to the list of secondary    100 TeV
Bremsstrahlung energy threshold above which primary
  muon/hadron is added to the list of secondary    100 TeV
Lowest triplet kinetic energy                      1 MeV
Enable sampling of gamma linear polarisation       0
5D gamma conversion model type                     0
5D gamma conversion model on isolated ion          0
Livermore data directory                           epics_2017
=======================================================================
======                 Ionisation Parameters                   ========
=======================================================================
Step function for e+-                              (0.2, 1 mm)
Step function for muons/hadrons                    (0.2, 0.1 mm)
Step function for light ions                       (0.2, 0.1 mm)
Step function for general ions                     (0.2, 0.1 mm)
Lowest e+e- kinetic energy                         1 keV
Lowest muon/hadron kinetic energy                  1 keV
Use ICRU90 data                                    0
Fluctuations of dE/dx are enabled                  1
Type of fluctuation model for leptons and hadrons  Urban
Use built-in Birks satuaration                     0
Build CSDA range enabled                           0
Use cut as a final range enabled                   0
Enable angular generator interface                 0
Max kinetic energy for CSDA tables                 1 GeV
Max kinetic energy for NIEL computation            0 eV
Linear loss limit                                  0.01
Read data from file for e+e- pair production by mu 0
=======================================================================
======                 Multiple Scattering Parameters          ========
=======================================================================
Type of msc step limit algorithm for e+-           1
Type of msc step limit algorithm for muons/hadrons 0
Msc lateral displacement for e+- enabled           1
Msc lateral displacement for muons and hadrons     0
Urban msc model lateral displacement alg96         1
Range factor for msc step limit for e+-            0.04
Range factor for msc step limit for muons/hadrons  0.2
Geometry factor for msc step limitation of e+-     2.5
Safety factor for msc step limit for e+-           0.6
Skin parameter for msc step limitation of e+-      1
Lambda limit for msc step limit for e+-            1 mm
Use Mott correction for e- scattering              0
Factor used for dynamic computation of angular
  limit between single and multiple scattering     1
Fixed angular limit between single
  and multiple scattering                          3.1416 rad
Upper energy limit for e+- multiple scattering     100 MeV
Type of electron single scattering model           0
Type of nuclear form-factor                        1
Screening factor                                   1
=======================================================================

phot:  for gamma SubType=12 BuildTable=0
      LambdaPrime table from 200 keV to 100 TeV in 61 bins
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
 LivermorePhElectric : Emin=    0 eV  Emax=  100 TeV  SauterGavrila Fluo

compt:  for gamma SubType=13 BuildTable=1
      Lambda table from 100 eV  to 1 MeV, 7 bins/decade, spline: 1
      LambdaPrime table from 1 MeV to 100 TeV in 56 bins
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
       Klein-Nishina : Emin=    0 eV  Emax=  100 TeV

conv:  for gamma SubType=14 BuildTable=1
      Lambda table from 1.022 MeV to 100 TeV, 18 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
     BetheHeitlerLPM : Emin=    0 eV  Emax=  100 TeV  ModifiedTsai

Rayl:  for gamma SubType=11 BuildTable=1
      Lambda table from 100 eV  to 150 keV, 7 bins/decade, spline: 0
      LambdaPrime table from 150 keV to 100 TeV in 62 bins
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh : Emin=    0 eV  Emax=  100 TeV  CullenGenerator
[email protected]:54: Building physics table for e-

msc:  for e-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 eV  Emax=  100 MeV Nbins=42 100 eV  - 100 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
        WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=42 100 MeV - 100 TeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

eIoni:  for e-  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        MollerBhabha : Emin=    0 eV  Emax=  100 TeV

eBrem:  for e-  XStype:4  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 eV  Emax=    1 GeV  ModifiedTsai
            eBremLPM : Emin=    1 GeV Emax=  100 TeV  ModifiedTsai

CoulombScat:  for e- XStype:1 SubType=1 BuildTable=1
      Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=  100 MeV Emax=  100 TeV
[email protected]:54: Building physics table for e+

msc:  for e+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 eV  Emax=  100 MeV Nbins=42 100 eV  - 100 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
        WentzelVIUni : Emin=  100 MeV Emax=  100 TeV Nbins=42 100 MeV - 100 TeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

eIoni:  for e+  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        MollerBhabha : Emin=    0 eV  Emax=  100 TeV

eBrem:  for e+  XStype:4  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      LPM flag: 1 for E > 1 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 eV  Emax=    1 GeV  ModifiedTsai
            eBremLPM : Emin=    1 GeV Emax=  100 TeV  ModifiedTsai

annihil:  for e+ XStype:2 SubType=5 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            eplus2gg : Emin=    0 eV  Emax=  100 TeV

CoulombScat:  for e+ XStype:1 SubType=1 BuildTable=1
      Lambda table from 100 MeV to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=  100 MeV Emax=  100 TeV

msc:  for proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for proton  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

hBrems:  for proton  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for proton  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 17x1001; from 7.50618 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for proton XStype:1 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for GenericIon  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 eV  Emax=  100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

ionIoni:  for GenericIon  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            BraggIon : Emin=    0 eV  Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

msc:  for alpha  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            UrbanMsc : Emin=    0 eV  Emax=  100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

ionIoni:  for alpha  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.02
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            BraggIon : Emin=    0 eV  Emax=7.9452 MeV
          BetheBloch : Emin=7.9452 MeV Emax=  100 TeV

msc:  for anti_proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for anti_proton  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=    2 MeV
          BetheBloch : Emin=    2 MeV Emax=  100 TeV

hBrems:  for anti_proton  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for anti_proton  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 17x1001; from 7.50618 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for anti_proton XStype:1 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for kaon+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for kaon+  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=1.05231 MeV
          BetheBloch : Emin=1.05231 MeV Emax=  100 TeV

hBrems:  for kaon+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for kaon+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 18x1001; from 3.94942 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for kaon+ XStype:1 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for kaon-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for kaon-  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=1.05231 MeV
          BetheBloch : Emin=1.05231 MeV Emax=  100 TeV

hBrems:  for kaon-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for kaon-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 18x1001; from 3.94942 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for kaon- XStype:1 SubType=1 BuildTable=1
      Used Lambda table of kaon+
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for mu+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

muIoni:  for mu+  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=  200 keV
        MuBetheBloch : Emin=  200 keV Emax=  100 TeV

muBrems:  for mu+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

muPairProd:  for mu+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 21x1001; from 0.85 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for mu+ XStype:1 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for mu-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

muIoni:  for mu-  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=  200 keV
        MuBetheBloch : Emin=  200 keV Emax=  100 TeV

muBrems:  for mu-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
              MuBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

muPairProd:  for mu-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 21x1001; from 0.85 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          muPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for mu- XStype:1 SubType=1 BuildTable=1
      Used Lambda table of mu+
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for pi+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for pi+  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=297.505 keV
          BetheBloch : Emin=297.505 keV Emax=  100 TeV

hBrems:  for pi+  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for pi+  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 20x1001; from 1.11656 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for pi+ XStype:1 SubType=1 BuildTable=1
      Lambda table from threshold  to 100 TeV, 7 bins/decade, spline: 0
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

msc:  for pi-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  100 TeV Nbins=84 100 eV  - 100 TeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=1 Llim=1 mm

hIoni:  for pi-  XStype:3  SubType=2
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      StepFunction=(0.2, 0.1 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=297.505 keV
          BetheBloch : Emin=297.505 keV Emax=  100 TeV

hBrems:  for pi-  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

hPairProd:  for pi-  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 84 bins
      Lambda tables from threshold to 100 TeV, 7 bins/decade, spline: 1
      Sampling table 20x1001; from 1.11656 GeV to 100 TeV
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
           hPairProd : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi

CoulombScat:  for pi- XStype:1 SubType=1 BuildTable=1
      Used Lambda table of pi+
      ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

====================================================================
                  HADRONIC PROCESSES SUMMARY (verbose level 1)
-----------------------------------------------------------------------
                           Hadronic Processes for neutron
  Process: hadElastic
        Model:             hElasticCHIPS: 0 eV  ---> 100 TeV
     Cr_sctns:        G4NeutronElasticXS: 0 eV  ---> 100 TeV
  Process: neutronInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:      G4NeutronInelasticXS: 0 eV  ---> 100 TeV
  Process: nCapture
        Model:               nRadCapture: 0 eV  ---> 100 TeV
     Cr_sctns:        G4NeutronCaptureXS: 0 eV  ---> 100 TeV
  Process: nKiller
-----------------------------------------------------------------------
                           Hadronic Processes for B-
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: B-Inelastic
        Model:                      FTFP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for D-
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: D-Inelastic
        Model:                      FTFP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for GenericIon
  Process: ionInelastic
        Model:  Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
-----------------------------------------------------------------------
                           Hadronic Processes for He3
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
  Process: He3Inelastic
        Model:  Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
-----------------------------------------------------------------------
                           Hadronic Processes for alpha
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
  Process: alphaInelastic
        Model:  Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
-----------------------------------------------------------------------
                           Hadronic Processes for anti_He3
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_He3Inelastic
        Model:                      FTFP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_alpha
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_alphaInelastic
        Model:                      FTFP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_deuteron
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_deuteronInelastic
        Model:                      FTFP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-------------------------------------------------------------------------
                           Hadronic Processes for anti_hypertriton
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_lambda
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: anti_lambdaInelastic
        Model:                      FTFP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_neutron
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100.1 MeV
        Model:              AntiAElastic: 100 MeV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_neutronInelastic
        Model:                      FTFP: 0 eV  ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_proton
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100.1 MeV
        Model:              AntiAElastic: 100 MeV ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_protonInelastic
        Model:                      FTFP: 0 eV  ---> 100 TeV
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for anti_triton
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
        Model:              AntiAElastic: 100 MeV/n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: anti_tritonInelastic
        Model:                      FTFP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:              AntiAGlauber: 0 eV  ---> 25.6 PeV
  Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for deuteron
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
  Process: dInelastic
        Model:  Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
-----------------------------------------------------------------------
                           Hadronic Processes for e+
  Process: positronNuclear
        Model:   G4ElectroVDNuclearModel: 0 eV  ---> 1 PeV
     Cr_sctns:          ElectroNuclearXS: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for e-
  Process: electronNuclear
        Model:   G4ElectroVDNuclearModel: 0 eV  ---> 1 PeV
     Cr_sctns:          ElectroNuclearXS: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for gamma
  Process: photonNuclear
        Model:               GammaNPreco: 0 eV  ---> 200 MeV
        Model:            BertiniCascade: 199 MeV ---> 6 GeV
        Model:           TheoFSGenerator: 3 GeV ---> 100 TeV
     Cr_sctns:            GammaNuclearXS: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for kaon+
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: kaon+Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for kaon-
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: kaon-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for lambda
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: lambdaInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for mu+
  Process: muonNuclear
        Model:      G4MuonVDNuclearModel: 0 eV  ---> 1 PeV
     Cr_sctns:     KokoulinMuonNuclearXS: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for mu-
  Process: muonNuclear
        Model:      G4MuonVDNuclearModel: 0 eV  ---> 1 PeV
     Cr_sctns:     KokoulinMuonNuclearXS: 0 eV  ---> 100 TeV
  Process: muMinusCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for pi+
  Process: hadElastic
        Model:           hElasticGlauber: 0 eV  ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
  Process: pi+Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for pi-
  Process: hadElastic
        Model:           hElasticGlauber: 0 eV  ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
  Process: pi-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
  Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for proton
  Process: hadElastic
        Model:             hElasticCHIPS: 0 eV  ---> 100 TeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
  Process: protonInelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:  BarashenkovGlauberGribov: 0 eV  ---> 100 TeV
-----------------------------------------------------------------------
                           Hadronic Processes for sigma-
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV  ---> 100 TeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: sigma-Inelastic
        Model:                      FTFP: 3 GeV ---> 100 TeV
        Model:            BertiniCascade: 0 eV  ---> 6 GeV
     Cr_sctns:            Glauber-Gribov: 0 eV  ---> 100 TeV
  Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
                           Hadronic Processes for triton
  Process: hadElastic
        Model:              hElasticLHEP: 0 eV /n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
  Process: tInelastic
        Model:  Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
        Model:                      FTFP: 3 GeV/n ---> 100 TeV/n
     Cr_sctns:  Glauber-Gribov Nucl-nucl: 0 eV  ---> 25.6 PeV
=======================================================================
======       Geant4 Native Pre-compound Model Parameters       ========
=======================================================================
Type of pre-compound inverse x-section              3
Pre-compound model active                           1
Pre-compound excitation low energy                  100 keV
Pre-compound excitation high energy                 30 MeV
Angular generator for pre-compound model            1
Use NeverGoBack option for pre-compound model       0
Use SoftCutOff option for pre-compound model        0
Use CEM transitions for pre-compound model          1
Use GNASH transitions for pre-compound model        0
Use HETC submodel for pre-compound model            0
=======================================================================
======       Nuclear De-excitation Module Parameters           ========
=======================================================================
Type of de-excitation inverse x-section             3
Type of de-excitation factory                       Evaporation+GEM
Number of de-excitation channels                    68
Min excitation energy                               10 eV
Min energy per nucleon for multifragmentation       200 GeV
Limit excitation energy for Fermi BreakUp           20 MeV
Level density (1/MeV)                               0.075
Use simple level density model                      1
Use discrete excitation energy of the residual      1
Time limit for long lived isomeres                  1 ns
Isomer production flag                              1
Internal e- conversion flag                         1
Store e- internal conversion data                   0
Correlated gamma emission flag                      0
Max 2J for sampling of angular correlations         10
=======================================================================
DDcelerRunAction INFO  Begin of run
[email protected]:105: Initializing shared params
[email protected]:100: Disabling GPU support since CUDA and HIP are disabled
[email protected]:31: Activating Celeritas version 0.0.0+45675e37 on CPU
[email protected]:1204: Transferring data from Geant4
[email protected]:375: Loaded 3 particles
[email protected]:434: Loaded 40 isotopes
[email protected]:487: Loaded 14 elements
[email protected]:661: Loaded 16 geo materials
[email protected]:765: Loaded 16 physics materials
[email protected]:345: Saving MSC model 'UrbanMsc' for particle e+ (-11)
[email protected]:345: Saving MSC model 'WentzelVIUni' for particle e+ (-11)
[email protected]:119: Saving process 'eIoni' for particle e+ (-11)
[email protected]:119: Saving process 'eBrem' for particle e+ (-11)
[email protected]:119: Saving process 'annihil' for particle e+ (-11)
[email protected]:119: Saving process 'CoulombScat' for particle e+ (-11)
[email protected]:345: Saving MSC model 'UrbanMsc' for particle e- (11)
[email protected]:345: Saving MSC model 'WentzelVIUni' for particle e- (11)
[email protected]:119: Saving process 'eIoni' for particle e- (11)
[email protected]:119: Saving process 'eBrem' for particle e- (11)
[email protected]:119: Saving process 'CoulombScat' for particle e- (11)
[email protected]:119: Saving process 'phot' for particle gamma (22)
[email protected]:119: Saving process 'compt' for particle gamma (22)
[email protected]:119: Saving process 'conv' for particle gamma (22)
[email protected]:119: Saving process 'Rayl' for particle gamma (22)
[email protected]:968: Loaded 11 processes
[email protected]:969: Loaded 0 optical models
[email protected]:798: Loaded 6 regions
[email protected]:1347: Loaded 1677 of 1677 volumes
[email protected]:104: ... 0.039671 s
[email protected]:1272: Loading external elemental data
[email protected]:104: ... 0.0423031 s
[email protected]:143: Using a uniform field: 3 [T] in all volumes
[email protected]:383: Initializing problem
status: Converting Geant4 geometry
[email protected]:104: ... 0.139454 s
[email protected]:461: Initializing tracking information
[email protected]:268: Omitting coulomb_scat from physics process list
[email protected]:150: Ignored process class coulomb_scat
[email protected]:256: Set default track order none
[email protected]:368: Celeritas core setup complete
[email protected]:32: Mapped sensitive detector "PixelEndcapN" on logical volume "component0Silicon"@0x16200d30 (ID=0) to VecGeom volume "component0Silicon@0" (ID=0)
[email protected]:32: Mapped sensitive detector "PixelEndcapN" on logical volume "component0Silicon"@0x1620cbe0 (ID=3) to VecGeom volume "component0Silicon@1" (ID=3)
[email protected]:32: Mapped sensitive detector "PixelEndcapP" on logical volume "component2Silicon"@0x1620d390 (ID=8) to VecGeom volume "component2Silicon@0" (ID=8)
[email protected]:32: Mapped sensitive detector "PixelEndcapP" on logical volume "component2Silicon"@0x1620d860 (ID=11) to VecGeom volume "component2Silicon@1" (ID=11)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x1620da00 (ID=12) to VecGeom volume "component0Silicon@2" (ID=12)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x1620de20 (ID=15) to VecGeom volume "component0Silicon@3" (ID=15)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapN" on logical volume "component0Silicon"@0x1620e270 (ID=18) to VecGeom volume "component0Silicon@4" (ID=18)
[email protected]:32: Mapped sensitive detector "ShortStripBarrel" on logical volume "component3"@0x16210330 (ID=25) to VecGeom volume "component3@0" (ID=29)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x16210520 (ID=26) to VecGeom volume "component0Silicon@5" (ID=30)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x16210970 (ID=29) to VecGeom volume "component0Silicon@6" (ID=33)
[email protected]:32: Mapped sensitive detector "ShortStripEndcapP" on logical volume "component0Silicon"@0x16210dc0 (ID=32) to VecGeom volume "component0Silicon@7" (ID=36)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component0Silicon"@0x16211210 (ID=35) to VecGeom volume "component0Silicon@8" (ID=39)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component4Silicon"@0x16212570 (ID=40) to VecGeom volume "component4Silicon@0" (ID=48)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component0Silicon"@0x16212790 (ID=41) to VecGeom volume "component0Silicon@9" (ID=49)
[email protected]:32: Mapped sensitive detector "LongStripEndcapN" on logical volume "component4Silicon"@0x16212ef0 (ID=46) to VecGeom volume "component4Silicon@1" (ID=58)
[email protected]:32: Mapped sensitive detector "LongStripBarrel" on logical volume "component0"@0x16213070 (ID=47) to VecGeom volume "component0@2" (ID=59)
[email protected]:32: Mapped sensitive detector "LongStripBarrel" on logical volume "component9"@0x162141e0 (ID=58) to VecGeom volume "component9" (ID=78)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component0Silicon"@0x16216a00 (ID=80) to VecGeom volume "component0Silicon@10" (ID=100)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component4Silicon"@0x16217190 (ID=85) to VecGeom volume "component4Silicon@2" (ID=109)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component0Silicon"@0x16217320 (ID=86) to VecGeom volume "component0Silicon@11" (ID=110)
[email protected]:32: Mapped sensitive detector "LongStripEndcapP" on logical volume "component4Silicon"@0x16217a50 (ID=91) to VecGeom volume "component4Silicon@3" (ID=119)
[email protected]:32: Mapped sensitive detector "PixelBarrel" on logical volume "component0"@0x16217bd0 (ID=92) to VecGeom volume "component0@0" (ID=120)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16221f80 (ID=150) to VecGeom volume "slice4@48" (ID=178)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16223be0 (ID=158) to VecGeom volume "slice4@49" (ID=186)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16224740 (ID=166) to VecGeom volume "slice4@50" (ID=194)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162251b0 (ID=174) to VecGeom volume "slice4@51" (ID=202)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16225b70 (ID=182) to VecGeom volume "slice4@52" (ID=210)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16226700 (ID=190) to VecGeom volume "slice4@53" (ID=218)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16227100 (ID=198) to VecGeom volume "slice4@54" (ID=226)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16228140 (ID=206) to VecGeom volume "slice4@55" (ID=234)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16228c60 (ID=214) to VecGeom volume "slice4@56" (ID=242)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16229930 (ID=222) to VecGeom volume "slice4@57" (ID=250)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622a330 (ID=230) to VecGeom volume "slice4@58" (ID=258)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622ad30 (ID=238) to VecGeom volume "slice4@59" (ID=266)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622b730 (ID=246) to VecGeom volume "slice4@60" (ID=274)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622c130 (ID=254) to VecGeom volume "slice4@61" (ID=282)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622cb30 (ID=262) to VecGeom volume "slice4@62" (ID=290)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622d530 (ID=270) to VecGeom volume "slice4@63" (ID=298)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622e920 (ID=278) to VecGeom volume "slice4@64" (ID=306)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622f7f0 (ID=286) to VecGeom volume "slice4@65" (ID=314)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162301f0 (ID=294) to VecGeom volume "slice4@66" (ID=322)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16230bf0 (ID=302) to VecGeom volume "slice4@67" (ID=330)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162315f0 (ID=310) to VecGeom volume "slice4@68" (ID=338)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16231ff0 (ID=318) to VecGeom volume "slice4@69" (ID=346)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162329f0 (ID=326) to VecGeom volume "slice4@70" (ID=354)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162333f0 (ID=334) to VecGeom volume "slice4@71" (ID=362)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16233df0 (ID=342) to VecGeom volume "slice4@72" (ID=370)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162347f0 (ID=350) to VecGeom volume "slice4@73" (ID=378)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162351f0 (ID=358) to VecGeom volume "slice4@74" (ID=386)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16235bf0 (ID=366) to VecGeom volume "slice4@75" (ID=394)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162365f0 (ID=374) to VecGeom volume "slice4@76" (ID=402)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16236ff0 (ID=382) to VecGeom volume "slice4@77" (ID=410)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162379f0 (ID=390) to VecGeom volume "slice4@78" (ID=418)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162383f0 (ID=398) to VecGeom volume "slice4@79" (ID=426)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1622e350 (ID=406) to VecGeom volume "slice4@80" (ID=434)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623b9e0 (ID=414) to VecGeom volume "slice4@81" (ID=442)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623c3e0 (ID=422) to VecGeom volume "slice4@82" (ID=450)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623cde0 (ID=430) to VecGeom volume "slice4@83" (ID=458)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623d7e0 (ID=438) to VecGeom volume "slice4@84" (ID=466)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623e1e0 (ID=446) to VecGeom volume "slice4@85" (ID=474)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623ebe0 (ID=454) to VecGeom volume "slice4@86" (ID=482)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623f5e0 (ID=462) to VecGeom volume "slice4@87" (ID=490)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x1623ffe0 (ID=470) to VecGeom volume "slice4@88" (ID=498)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162409e0 (ID=478) to VecGeom volume "slice4@89" (ID=506)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162413e0 (ID=486) to VecGeom volume "slice4@90" (ID=514)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16241de0 (ID=494) to VecGeom volume "slice4@91" (ID=522)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162427e0 (ID=502) to VecGeom volume "slice4@92" (ID=530)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x162431e0 (ID=510) to VecGeom volume "slice4@93" (ID=538)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16207070 (ID=518) to VecGeom volume "slice4@94" (ID=546)
[email protected]:32: Mapped sensitive detector "ECalEndcap" on logical volume "slice4"@0x16207a70 (ID=526) to VecGeom volume "slice4@95" (ID=554)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162092f0 (ID=533) to VecGeom volume "slice3@132" (ID=561)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16209cf0 (ID=537) to VecGeom volume "slice3@133" (ID=565)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620a1f0 (ID=541) to VecGeom volume "slice3@134" (ID=569)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620a710 (ID=545) to VecGeom volume "slice3@135" (ID=573)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620ac30 (ID=549) to VecGeom volume "slice3@136" (ID=577)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620b150 (ID=553) to VecGeom volume "slice3@137" (ID=581)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620b670 (ID=557) to VecGeom volume "slice3@138" (ID=585)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620bb90 (ID=561) to VecGeom volume "slice3@139" (ID=589)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620c0b0 (ID=565) to VecGeom volume "slice3@140" (ID=593)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1620c5d0 (ID=569) to VecGeom volume "slice3@141" (ID=597)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1624f690 (ID=573) to VecGeom volume "slice3@142" (ID=601)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x1624fbb0 (ID=577) to VecGeom volume "slice3@143" (ID=605)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162500d0 (ID=581) to VecGeom volume "slice3@144" (ID=609)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162505f0 (ID=585) to VecGeom volume "slice3@145" (ID=613)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16250b10 (ID=589) to VecGeom volume "slice3@146" (ID=617)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16251030 (ID=593) to VecGeom volume "slice3@147" (ID=621)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16251960 (ID=597) to VecGeom volume "slice3@148" (ID=625)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162526a0 (ID=601) to VecGeom volume "slice3@149" (ID=629)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16252bc0 (ID=605) to VecGeom volume "slice3@150" (ID=633)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162530e0 (ID=609) to VecGeom volume "slice3@151" (ID=637)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16253600 (ID=613) to VecGeom volume "slice3@152" (ID=641)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16253b20 (ID=617) to VecGeom volume "slice3@153" (ID=645)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16254040 (ID=621) to VecGeom volume "slice3@154" (ID=649)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16254560 (ID=625) to VecGeom volume "slice3@155" (ID=653)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16254a80 (ID=629) to VecGeom volume "slice3@156" (ID=657)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16254fa0 (ID=633) to VecGeom volume "slice3@157" (ID=661)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162554c0 (ID=637) to VecGeom volume "slice3@158" (ID=665)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162559e0 (ID=641) to VecGeom volume "slice3@159" (ID=669)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16255f00 (ID=645) to VecGeom volume "slice3@160" (ID=673)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16256420 (ID=649) to VecGeom volume "slice3@161" (ID=677)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16256940 (ID=653) to VecGeom volume "slice3@162" (ID=681)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16256e60 (ID=657) to VecGeom volume "slice3@163" (ID=685)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16257b90 (ID=661) to VecGeom volume "slice3@164" (ID=689)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162580b0 (ID=665) to VecGeom volume "slice3@165" (ID=693)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x162585d0 (ID=669) to VecGeom volume "slice3@166" (ID=697)
[email protected]:32: Mapped sensitive detector "HCalEndcap" on logical volume "slice3"@0x16258af0 (ID=673) to VecGeom volume "slice3@167" (ID=701)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625b960 (ID=704) to VecGeom volume "slice4@0" (ID=736)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625c340 (ID=712) to VecGeom volume "slice4@1" (ID=744)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625cd40 (ID=720) to VecGeom volume "slice4@2" (ID=752)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625d740 (ID=728) to VecGeom volume "slice4@3" (ID=760)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625e140 (ID=736) to VecGeom volume "slice4@4" (ID=768)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625eb40 (ID=744) to VecGeom volume "slice4@5" (ID=776)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625f540 (ID=752) to VecGeom volume "slice4@6" (ID=784)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1625ff40 (ID=760) to VecGeom volume "slice4@7" (ID=792)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16260940 (ID=768) to VecGeom volume "slice4@8" (ID=800)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16261340 (ID=776) to VecGeom volume "slice4@9" (ID=808)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16261d40 (ID=784) to VecGeom volume "slice4@10" (ID=816)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16262740 (ID=792) to VecGeom volume "slice4@11" (ID=824)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16263140 (ID=800) to VecGeom volume "slice4@12" (ID=832)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16239070 (ID=808) to VecGeom volume "slice4@13" (ID=840)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16239a70 (ID=816) to VecGeom volume "slice4@14" (ID=848)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16266850 (ID=824) to VecGeom volume "slice4@15" (ID=856)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16267a60 (ID=832) to VecGeom volume "slice4@16" (ID=864)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16268c80 (ID=840) to VecGeom volume "slice4@17" (ID=872)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16269680 (ID=848) to VecGeom volume "slice4@18" (ID=880)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626a080 (ID=856) to VecGeom volume "slice4@19" (ID=888)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626aa80 (ID=864) to VecGeom volume "slice4@20" (ID=896)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626b480 (ID=872) to VecGeom volume "slice4@21" (ID=904)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626be80 (ID=880) to VecGeom volume "slice4@22" (ID=912)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626c880 (ID=888) to VecGeom volume "slice4@23" (ID=920)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626d280 (ID=896) to VecGeom volume "slice4@24" (ID=928)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626dc80 (ID=904) to VecGeom volume "slice4@25" (ID=936)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626e680 (ID=912) to VecGeom volume "slice4@26" (ID=944)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626f080 (ID=920) to VecGeom volume "slice4@27" (ID=952)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x1626fa80 (ID=928) to VecGeom volume "slice4@28" (ID=960)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16270480 (ID=936) to VecGeom volume "slice4@29" (ID=968)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16270e80 (ID=944) to VecGeom volume "slice4@30" (ID=976)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16271880 (ID=952) to VecGeom volume "slice4@31" (ID=984)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16267250 (ID=960) to VecGeom volume "slice4@32" (ID=992)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162734b0 (ID=968) to VecGeom volume "slice4@33" (ID=1000)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16273eb0 (ID=976) to VecGeom volume "slice4@34" (ID=1008)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162748b0 (ID=984) to VecGeom volume "slice4@35" (ID=1016)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162752b0 (ID=992) to VecGeom volume "slice4@36" (ID=1024)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16275cb0 (ID=1000) to VecGeom volume "slice4@37" (ID=1032)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162766b0 (ID=1008) to VecGeom volume "slice4@38" (ID=1040)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162770b0 (ID=1016) to VecGeom volume "slice4@39" (ID=1048)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16277ab0 (ID=1024) to VecGeom volume "slice4@40" (ID=1056)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16243e60 (ID=1032) to VecGeom volume "slice4@41" (ID=1064)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16244860 (ID=1040) to VecGeom volume "slice4@42" (ID=1072)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16245260 (ID=1048) to VecGeom volume "slice4@43" (ID=1080)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16247490 (ID=1056) to VecGeom volume "slice4@44" (ID=1088)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162486a0 (ID=1064) to VecGeom volume "slice4@45" (ID=1096)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x162490a0 (ID=1072) to VecGeom volume "slice4@46" (ID=1104)
[email protected]:32: Mapped sensitive detector "ECalBarrel" on logical volume "slice4"@0x16249aa0 (ID=1080) to VecGeom volume "slice4@47" (ID=1112)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628a1d0 (ID=1135) to VecGeom volume "slice3@96" (ID=1167)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628a6b0 (ID=1139) to VecGeom volume "slice3@97" (ID=1171)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628abd0 (ID=1143) to VecGeom volume "slice3@98" (ID=1175)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628b0f0 (ID=1147) to VecGeom volume "slice3@99" (ID=1179)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628b610 (ID=1151) to VecGeom volume "slice3@100" (ID=1183)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628bb30 (ID=1155) to VecGeom volume "slice3@101" (ID=1187)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628c050 (ID=1159) to VecGeom volume "slice3@102" (ID=1191)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628c570 (ID=1163) to VecGeom volume "slice3@103" (ID=1195)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628ca90 (ID=1167) to VecGeom volume "slice3@104" (ID=1199)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628cfb0 (ID=1171) to VecGeom volume "slice3@105" (ID=1203)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628d4d0 (ID=1175) to VecGeom volume "slice3@106" (ID=1207)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628d9f0 (ID=1179) to VecGeom volume "slice3@107" (ID=1211)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628df10 (ID=1183) to VecGeom volume "slice3@108" (ID=1215)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628e430 (ID=1187) to VecGeom volume "slice3@109" (ID=1219)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628e950 (ID=1191) to VecGeom volume "slice3@110" (ID=1223)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628ee70 (ID=1195) to VecGeom volume "slice3@111" (ID=1227)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628f390 (ID=1199) to VecGeom volume "slice3@112" (ID=1231)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628f8b0 (ID=1203) to VecGeom volume "slice3@113" (ID=1235)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x1628fdd0 (ID=1207) to VecGeom volume "slice3@114" (ID=1239)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x162902f0 (ID=1211) to VecGeom volume "slice3@115" (ID=1243)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16290810 (ID=1215) to VecGeom volume "slice3@116" (ID=1247)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16290d30 (ID=1219) to VecGeom volume "slice3@117" (ID=1251)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16291250 (ID=1223) to VecGeom volume "slice3@118" (ID=1255)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16291770 (ID=1227) to VecGeom volume "slice3@119" (ID=1259)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16291c90 (ID=1231) to VecGeom volume "slice3@120" (ID=1263)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x162921b0 (ID=1235) to VecGeom volume "slice3@121" (ID=1267)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x162926d0 (ID=1239) to VecGeom volume "slice3@122" (ID=1271)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16292bf0 (ID=1243) to VecGeom volume "slice3@123" (ID=1275)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16293110 (ID=1247) to VecGeom volume "slice3@124" (ID=1279)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16293630 (ID=1251) to VecGeom volume "slice3@125" (ID=1283)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16293b50 (ID=1255) to VecGeom volume "slice3@126" (ID=1287)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16294070 (ID=1259) to VecGeom volume "slice3@127" (ID=1291)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16294da0 (ID=1263) to VecGeom volume "slice3@128" (ID=1295)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x162952c0 (ID=1267) to VecGeom volume "slice3@129" (ID=1299)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x162957e0 (ID=1271) to VecGeom volume "slice3@130" (ID=1303)
[email protected]:32: Mapped sensitive detector "HCalBarrel" on logical volume "slice3"@0x16295d00 (ID=1275) to VecGeom volume "slice3@131" (ID=1307)
[email protected]:468: Setting number of streams to 1
[email protected]:491: Wrote Geant4 diagnostic output to "trackingmanager-offload.out.json"
[email protected]:104: ... 0.712984 s
[email protected]:173: Constructing local state
[email protected]:105: Setting up thread-local hit processor for 190 sensitive detectors
[email protected]:80: Celeritas core state initialization complete
[email protected]:157: Verifying tracking manager
Geant4Output2EDM4hep INFO  Opened ddsimOutput.root for output
SignalHandler    INFO  ++ Re-apply signal handlers
GenerationInit   INFO  +++ Initializing event 1. Within run:0 event 1.
Gun              INFO  Shoot [0] [0.000 , 10.000] GeV e- pos:(0.000 0.000 0.000)[mm] dir:( 0.000  0.000  1.000)
Gun              INFO  Particle [0] e-           Mom:1.953 GeV vertex:( 0.000  0.000  0.000)[mm] direction:( 0.067 -0.042  0.997)
Gun              INFO  +-> Interaction [0] [0.000 , 10.000] GeV e- pos:(0.000 0.000 0.000)[mm]
Gun              INFO  +++   +-> ID:  0 e-           status:00000002 PDG:    11 Vtx:(+0.00e+00,+0.00e+00,+0.00e+00)[mm] time: +0.00e+00 [ns] #Dau:  0 #Par:0
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

Thread 1 "python" received signal SIGSEGV, Segmentation fault.
G4Step::GetPreStepPoint (this=0x0) at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4Step.icc:37
37        return fpPreStepPoint;

Complete backtrace

(gdb) bt
#0  G4Step::GetPreStepPoint (this=0x0) at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4Step.icc:37
#1  0x00007fffd9cba099 in G4Track::GetMaterialCutsCouple (this=0x29b062b0)
    at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4Track.icc:192
#2  0x00007fffd9cb9ff1 in dd4hep::sim::Geant4StepHandler::birkAttenuation (this=0x7fffffff9c10)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4StepHandler.cpp:89
#3  0x00007fffd9c1bc6e in dd4hep::sim::Geant4StepHandler::totalEnergy (this=0x7fffffff9c10)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/include/DDG4/Geant4StepHandler.h:79
#4  0x00007fffd9c1a5d2 in dd4hep::sim::Geant4HitData::extractContribution (step=0x29b19f40, ApplyBirksLaw=true)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Data.cpp:82
#5  0x00007fffe03d4794 in dd4hep::sim::Geant4SensitiveAction<dd4hep::sim::Geant4ScintillatorCalorimeter>::process (
    this=0x156acac0, step=0x29b19f40) at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/plugins/Geant4SDActions.cpp:461
#6  0x00007fffd9caaa7c in dd4hep::sim::Geant4SensDetActionSequence::process (this=0x15692910, step=0x29b19f40,
    history=0x0) at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4SensDetAction.cpp:449
#7  0x00007fffe03df602 in dd4hep::sim::Geant4SensDet::ProcessHits (this=0xd279760, step=0x29b19f40, hist=0x0)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/plugins/Geant4SensDet.cpp:155
#8  0x00007fffd5fad3c6 in G4VSensitiveDetector::Hit (this=0xd279770, aStep=0x29b19f40)
    at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4VSensitiveDetector.hh:92
#9  0x00007fffd5fab5af in celeritas::detail::HitProcessor::operator() (this=0x299a4c20, out=..., i=0)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/ext/detail/HitProcessor.cc:374
#10 0x00007fffd5faa8aa in celeritas::detail::HitProcessor::operator() (this=0x299a4c20, out=...)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/ext/detail/HitProcessor.cc:271
#11 0x00007fffd5faa66a in celeritas::detail::HitProcessor::operator() (this=0x299a4c20, states=...)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/ext/detail/HitProcessor.cc:239
#12 0x00007fffd5f88c7a in celeritas::GeantSd::process_steps (this=0x2c9afa20, state=...)
--Type <RET> for more, q to quit, c to continue without paging--
    at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/ext/GeantSd.cc:176
#13 0x00007fffd5f5a6ce in celeritas::detail::StepGatherAction<(celeritas::StepPoint)1>::step (this=0x2c978770,
    params=..., state=...) at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/user/detail/StepGatherAction.cc:75
#14 0x00007fffd5c8d108 in celeritas::ActionSequence::step<(celeritas::MemSpace)0> (this=0x2c97aa00, params=...,
    state=...) at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/global/ActionSequence.cc:132
#15 0x00007fffd5ccf94d in celeritas::Stepper<(celeritas::MemSpace)0>::operator() (this=0x29691300)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/celeritas/global/Stepper.cc:135
#16 0x00007fffd6566291 in celeritas::LocalTransporter::Flush (this=0x192b6300)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/accel/LocalTransporter.cc:384

#17 0x00007fffd65a84cc in celeritas::TrackingManager::FlushEvent (this=0x15d6f4d0)
    at /home/srahman1/celeritas-dd4hep/celeritas/src/accel/TrackingManager.cc:164
#18 0x00007fffe096a4aa in G4EventManager::DoProcessing (this=0x156c9eb0, anEvent=<optimized out>)
    at /home/srahman1/celeritas-dd4hep/geant4/source/event/src/G4EventManager.cc:288
#19 0x00007fffe0a0f594 in G4RunManager::ProcessOneEvent (this=0x151ac610, i_event=0)
    at /home/srahman1/celeritas-dd4hep/geant4/source/run/src/G4RunManager.cc:404

#20 0x00007fffe0a0f3ac in G4RunManager::DoEventLoop (this=0x151ac610, n_event=<optimized out>,
    macroFile=<optimized out>, n_select=<optimized out>)
    at /home/srahman1/celeritas-dd4hep/geant4/source/run/src/G4RunManager.cc:370
#21 0x00007fffe0a0e94c in G4RunManager::BeamOn (this=0x151ac610, n_event=10, macroFile=0x0, n_select=-1)
    at /home/srahman1/celeritas-dd4hep/geant4/source/run/src/G4RunManager.cc:261
#22 0x00007fffd9cc591f in dd4hep::sim::Geant4UIManager::start (this=0x1507e400)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4UIManager.cpp:253

#23 0x00007fffd9cc502d in dd4hep::sim::Geant4UIManager::operator() (this=0x1507e400)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4UIManager.cpp:176
--Type <RET> for more, q to quit, c to continue without paging--

#24 0x00007fffd9c25dee in dd4hep::sim::Geant4Exec::run (kernel=...)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Exec.cpp:645

#25 0x00007fffd9c6b395 in dd4hep::sim::Geant4Kernel::run (this=0xcc481b0)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4Kernel.cpp:409
#26 0x00007fffd08a5012 in ?? ()
(gdb) break G4Step.icc:137

....................................................................

### Birks coefficients used in run time

Thread 1 "python" received signal SIGSEGV, Segmentation fault.
G4Step::GetPreStepPoint (this=0x0) at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4Step.icc:37
37        return fpPreStepPoint;
(gdb) print this
$1 = (const G4Step *) 0x0
(gdb) print *this
Cannot access memory at address 0x0
(gdb) print *step
$2 = {int (const char *, const char *)} 0x7ffff7dffe00 <step>
(gdb) print step
$3 = {int (const char *, const char *)} 0x7ffff7dffe00 <step>
(gdb)

rahmans1 avatar Jun 26 '25 15:06 rahmans1

(gdb) frame 1
#1  0x00007fffd9cba099 in G4Track::GetMaterialCutsCouple (this=0x29b062b0)
    at /home/srahman1/celeritas-dd4hep/geant4/install/include/Geant4/G4Track.icc:192
192       return fpStep->GetPreStepPoint()->GetMaterialCutsCouple();
(gdb) print step
$14 = {int (const char *, const char *)} 0x7ffff7dffe00 <step>
(gdb) print *step
$15 = {int (const char *, const char *)} 0x7ffff7dffe00 <step>
(gdb) print track
No symbol "track" in current context.
(gdb) print *track
No symbol "track" in current context.
(gdb) print this
$16 = (const G4Track *) 0x29b062b0
(gdb) print *this
$17 = {fPosition = {static ToleranceTicks = 100, data = {315.30225090937984, 44.617875305131129, 3288.4999999999995},
    static tolerance = 2.22045e-14}, fGlobalTime = 11.021136586226234, fLocalTime = 0, fTrackLength = 0,
  fVelocity = 299.79245800000001, fpTouchable = {fObj = 0x0}, fpNextTouchable = {fObj = 0x0}, fpOriginTouchable = {
    fObj = 0x0}, fpDynamicParticle = 0x199c74f0, fTrackStatus = fAlive, fStepLength = 0, fWeight = 1, fpStep = 0x0,
  fVtxPosition = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxMomentumDirection = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxKineticEnergy = 0, fpLVAtVertex = 0x0, fpCreatorProcess = 0x0, fpUserInformation = 0x0, prev_mat = 0x0,
  groupvel = 0x0, prev_velocity = 0, prev_momentum = 0, fpAuxiliaryTrackInformationMap = 0x0, fCurrentStepNumber = 0,
  fCreatorModelID = -1, fParentResonanceDef = 0x0, fParentResonanceID = 0, fParentID = -1, fTrackID = -1,
  fBelowThreshold = false, fGoodForTracking = false, is_OpticalPhoton = false, useGivenVelocity = false}
(gdb) print trk
No symbol "trk" in current context.
(gdb) print *trk
No symbol "trk" in current context.
(gdb)
(gdb) frame 2
#2  0x00007fffd9cb9ff1 in dd4hep::sim::Geant4StepHandler::birkAttenuation (this=0x7fffffff9c10)
    at /home/srahman1/celeritas-dd4hep/DD4hep/DDG4/src/Geant4StepHandler.cpp:89
89        const G4MaterialCutsCouple* couple = trk->GetMaterialCutsCouple();
(gdb) print *this
$4 = {<dd4hep::sim::Geant4HitHandler> = {track = 0x29b062b0, touchable_ptr = 0x29a545b0}, step = 0x295e59e0,
  pre = 0x29a5bb50, post = 0x2968ee60, applyBirksLaw = true}
(gdb) print this
$5 = (const dd4hep::sim::Geant4StepHandler *) 0x7fffffff9c10
(gdb) print step
$6 = (const G4Step *) 0x295e59e0
(gdb) print *step
$7 = {fTotalEnergyDeposit = 0.25599308920274544, fNonIonizingEnergyDeposit = -inf, fpPreStepPoint = 0x29a5bb50,
  fpPostStepPoint = 0x2968ee60, fStepLength = 0.50200461920594952, fpTrack = 0x29b062b0,
  fpSteppingControlFlag = NormalCondition, fFirstStepInVolume = false, fLastStepInVolume = false,
  fSecondary = 0x2b682ca0, nSecondaryByLastStep = 0, secondaryInCurrentStep = 0x2b682ae0,
  fpVectorOfAuxiliaryPointsPointer = 0x0}
(gdb) print trk
$8 = (const G4Track *) 0x29b062b0
(gdb) print *trk
$9 = {fPosition = {static ToleranceTicks = 100, data = {315.30225090937984, 44.617875305131129, 3288.4999999999995},
    static tolerance = 2.22045e-14}, fGlobalTime = 11.021136586226234, fLocalTime = 0, fTrackLength = 0,
  fVelocity = 299.79245800000001, fpTouchable = {fObj = 0x0}, fpNextTouchable = {fObj = 0x0}, fpOriginTouchable = {
    fObj = 0x0}, fpDynamicParticle = 0x199c74f0, fTrackStatus = fAlive, fStepLength = 0, fWeight = 1, fpStep = 0x0,
  fVtxPosition = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxMomentumDirection = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxKineticEnergy = 0, fpLVAtVertex = 0x0, fpCreatorProcess = 0x0, fpUserInformation = 0x0, prev_mat = 0x0,
  groupvel = 0x0, prev_velocity = 0, prev_momentum = 0, fpAuxiliaryTrackInformationMap = 0x0, fCurrentStepNumber = 0,
  fCreatorModelID = -1, fParentResonanceDef = 0x0, fParentResonanceID = 0, fParentID = -1, fTrackID = -1,
  fBelowThreshold = false, fGoodForTracking = false, is_OpticalPhoton = false, useGivenVelocity = false}

(gdb) print track
$10 = (const G4Track *) 0x29b062b0
(gdb) print *track
$11 = {fPosition = {static ToleranceTicks = 100, data = {315.30225090937984, 44.617875305131129, 3288.4999999999995},
    static tolerance = 2.22045e-14}, fGlobalTime = 11.021136586226234, fLocalTime = 0, fTrackLength = 0,
  fVelocity = 299.79245800000001, fpTouchable = {fObj = 0x0}, fpNextTouchable = {fObj = 0x0}, fpOriginTouchable = {
    fObj = 0x0}, fpDynamicParticle = 0x199c74f0, fTrackStatus = fAlive, fStepLength = 0, fWeight = 1, fpStep = 0x0,
  fVtxPosition = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxMomentumDirection = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxKineticEnergy = 0, fpLVAtVertex = 0x0, fpCreatorProcess = 0x0, fpUserInformation = 0x0, prev_mat = 0x0,
  groupvel = 0x0, prev_velocity = 0, prev_momentum = 0, fpAuxiliaryTrackInformationMap = 0x0, fCurrentStepNumber = 0,
  fCreatorModelID = -1, fParentResonanceDef = 0x0, fParentResonanceID = 0, fParentID = -1, fTrackID = -1,
  fBelowThreshold = false, fGoodForTracking = false, is_OpticalPhoton = false, useGivenVelocity = false}

(gdb) print step->GetTrack()
$12 = (G4Track *) 0x29b062b0
(gdb) print *(step->GetTrack())
$13 = {fPosition = {static ToleranceTicks = 100, data = {315.30225090937984, 44.617875305131129, 3288.4999999999995},
    static tolerance = 2.22045e-14}, fGlobalTime = 11.021136586226234, fLocalTime = 0, fTrackLength = 0,
  fVelocity = 299.79245800000001, fpTouchable = {fObj = 0x0}, fpNextTouchable = {fObj = 0x0}, fpOriginTouchable = {
    fObj = 0x0}, fpDynamicParticle = 0x199c74f0, fTrackStatus = fAlive, fStepLength = 0, fWeight = 1, fpStep = 0x0,
  fVtxPosition = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxMomentumDirection = {static ToleranceTicks = 100, data = {0, 0, 0}, static tolerance = 2.22045e-14},
  fVtxKineticEnergy = 0, fpLVAtVertex = 0x0, fpCreatorProcess = 0x0, fpUserInformation = 0x0, prev_mat = 0x0,
  groupvel = 0x0, prev_velocity = 0, prev_momentum = 0, fpAuxiliaryTrackInformationMap = 0x0, fCurrentStepNumber = 0,
  fCreatorModelID = -1, fParentResonanceDef = 0x0, fParentResonanceID = 0, fParentID = -1, fTrackID = -1,
  fBelowThreshold = false, fGoodForTracking = false, is_OpticalPhoton = false, useGivenVelocity = false}

rahmans1 avatar Jun 26 '25 17:06 rahmans1

With changes to the hit processor (https://github.com/celeritas-project/celeritas/pull/1835), we can see some hits in the ODD ecal if we ignore saving MC truth particle info in edm4hep output.

Screenshot 2025-06-27 123202

Now we need to consider the following:

  1. How do we reconcile celeritas hit processor with the dd4hep edm4hep writer that looks to retrieve a realistic track ID? Celeritas by default seems to be always returning -1.
  2. The tracker hit collections in ODD are returning empty. This needs to be understood.
  3. Try the EIC ePIC detector and see if the behavior can be reproduced.

rahmans1 avatar Jun 27 '25 16:06 rahmans1

Running with the track id returns the same error of non-matching truth record when trying to write a edm4hep output. There may be a truth record that needs to be propagated through for the edm4hep plugin to not complain. Might be something to discuss with the dd4hep expert. There is 3 plugins available based on different output configuration. The failure we are getting so far is with the EDM4HEP plugin. I haven't tried LCIO.

--outputConfig.forceDD4HEP Use the DD4HEP output plugin regardless of outputfilename. --outputConfig.forceEDM4HEP Use the EDM4HEP output plugin regardless of outputfilename. --outputConfig.forceLCIO Use the LCIO output plugin regardless of outputfilename.

Geant4Output2EDM4hep INFO  Opened output.edm4hep.root for output
SignalHandler    INFO  ++ Re-apply signal handlers
GenerationInit   INFO  +++ Initializing event 1. Within run:0 event 1.
Gun              INFO  Shoot [0] [0.000 , 10.000] GeV e- pos:(0.000 0.000 0.000)[mm] dir:( 0.000  0.000  1.000)
Gun              INFO  Particle [0] e-           Mom:9.240 GeV vertex:( 0.000  0.000  0.000)[mm] direction:(-0.011  0.492  0.870)
Gun              INFO  +-> Interaction [0] [0.000 , 10.000] GeV e- pos:(0.000 0.000 0.000)[mm]
Gun              INFO  +++   +-> ID:  0 e-           status:00000002 PDG:    11 Vtx:(+0.00e+00,+0.00e+00,+0.00e+00)[mm] time: +0.00e+00 [ns] #Dau:  0 #Par:0
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
[email protected]:397: Reconstituted 4600543818244423333 hits for event 0
EDM4hepOutput    WARN  +++ [Event:0] No valid MC truth info present. Is a Particle handler installed ?
Geant4Output2EDM4hep INFO  +++ Saving EDM4hep event 0 run 0.
Geant4ParticleMap ERROR +++ No Equivalent particle for track:0. Monte Carlo truth record looks broken!
Particle map:
 [0   :0x2952b9f0]
Equivalents:

EDM4hepOutput    ERROR +++ [Event:0] Exception while saving event:deque::_M_range_check: __n (which is 18446744073709551615)>= this->size() (which is 1)
EDM4hepOutput    ERROR +++ [Event:0] Exception while saving event:deque::_M_range_check: __n (which is 18446744073709551615)>= this->size() (which is 1)
Geant4Kernel     FATAL +++ Exception while simulating:deque::_M_range_check: __n (which is 18446744073709551615)>= this->size() (which is 1)
DDSim            ERROR Simulation failed!
Geant4Kernel     INFO  ++ Terminate Geant4 and delete associated actions.
DDSim            INFO Total Time:   29.29 s (User), 1.44 s (System)
DDSim            INFO StartUp Time: 7.32 s, Processing and Init: 21.97 s (~21.97 s/Event)
DDcelerRunAction INFO  Destructing Geant4RunAction
[Thread 0x7fffe26136c0 (LWP 373684) exited]
Fatal Python error: none_dealloc: deallocating None: bug likely caused by a refcount error in a C extension
Python runtime state: finalizing (tstate=0x0000000000a840f8)

Current thread 0x00007ffff7c85740 (most recent call first):
  Garbage-collecting
  <no Python frame>

Thread 1 "python" received signal SIGABRT, Aborted.
__pthread_kill_implementation (threadid=<optimized out>, signo=signo@entry=6, no_tid=no_tid@entry=0) at ./nptl/pthread_kill.c:44
44      ./nptl/pthread_kill.c: No such file or directory.
(gdb)

rahmans1 avatar Jul 01 '25 03:07 rahmans1

@rahmans1 I think a next step toward getting this ready is finding a simple calorimeter-type detector example that uses DD4HEP and integrate it into the celeritas "examples" directory.

sethrj avatar Jul 01 '25 10:07 sethrj

@sethrj Do you want a separate implementation of detector or may be download the ODD as a submodule if the ddceler plugin is enabled?

CELER_LOG=debug CELER_LOG_LOCAL=debug gdb --args python ~/celeritas-dd4hep/DD4hep/install/bin/ddsim --action.run DDcelerRunAction --compactFile=$PWD/install/share/OpenDataDetector/xml/OpenDataDetector.xml --steering steeringFile.py -N 1000 -G --gun.particle e- --gun.etaMax 4 --gun.etaMin 1 --gun.distribution uniform --outputFile output.root --outputConfig.forceDD4HEP

Btw I also did a little more work with the standard DD4hep output plugin on the ODD and wrote a small analysis script to make some plots. Should that go into benchmarks?

--outputConfig.forceDD4HEP Use the DD4HEP output plugin regardless of outputfilename.

image

image

image

rahmans1 avatar Jul 02 '25 15:07 rahmans1

A little confusing. When running with CELER_DISABLE=1, it seems to makes dd4hep ignore its internal flags.

CELER_DISABLE=1 CELER_LOG=debug CELER_LOG_LOCAL=debug ddsim --action.run DDcelerRunAction --compactFile=$PWD/install/share/CeleritasDD4hep/celeritas-dd4hep.xml --steering steeringFile.py -N 500 -G --gun.particle e- --gun.etaMax 4 --gun.etaMin 1 --gun.energy 5*GeV --gun.distribution uniform --outputFile output2.root --outputConfig.forceDD4HEP  --action.tracker "Geant4TrackerAction" --action.trackerSDTypes "tracker" --action.calorimeter "Geant4CalorimeterAction" --action.calorimeterSDTypes "calorimeter"
CELER_LOG=debug CELER_LOG_LOCAL=debug ddsim --action.run DDcelerRunAction --compactFile=$PWD/install/share/CeleritasDD4hep/celeritas-dd4hep.xml --steering steeringFile.py -N 500 -G --gun.particle e- --gun.etaMax 4 --gun.etaMin 1 --gun.energy 5*GeV --gun.distribution uniform --outputFile output1.root --outputConfig.forceDD4HEP  --action.tracker "Geant4TrackerAction" --action.trackerSDTypes "tracker" --action.calorimeter "Geant4CalorimeterAction" --action.calorimeterSDTypes "calorimeter"
image For example, configuration 1 (celer+G4) eta looks range limited as expected but the configuration 2 (no celer) looks like eta is unrestricted.

However, the energy deposit in calorimeter looks similar Calorimeter_Energy_Deposit

rahmans1 avatar Jul 21 '25 18:07 rahmans1

Ok. I see the issue now. It's not ignoring the internal flags. But when running without celeritas, there are a lot of secondaries that end up in the MCparticle collection. If i cut those out, i get the same distribution. Now the question is why these don't show up in the MCparticle distribution in the celeritas offload case. I am guessing that none of the secondary MCParticle information is propagated through celeritas.

image

rahmans1 avatar Jul 25 '25 15:07 rahmans1

Calorimeter edep comparison using logarithmically increasing bin-width along the x-axis.

We seem to missing a lot of stuff that deposits less than 90 MeV energy in the calorimeters.

image

rahmans1 avatar Jul 31 '25 17:07 rahmans1

As discussed this morning, we're going to reduce the scope of this PR so that it does not include full field support: we hopefully can query DD4HEP for the field types and see if they're uniform and constant.

sethrj avatar Oct 07 '25 18:10 sethrj

Codecov Report

:white_check_mark: All modified and coverable lines are covered by tests. :white_check_mark: Project coverage is 83.43%. Comparing base (068a4ea) to head (cf3684d).

Additional details and impacted files

Impacted file tree graph

@@             Coverage Diff             @@
##           develop    #1756      +/-   ##
===========================================
- Coverage    84.87%   83.43%   -1.45%     
===========================================
  Files         1255     1235      -20     
  Lines        44185    41895    -2290     
  Branches     16501    13194    -3307     
===========================================
- Hits         37503    34953    -2550     
+ Misses        4725     4634      -91     
- Partials      1957     2308     +351     

see 685 files with indirect coverage changes

:rocket: New features to boost your workflow:
  • :snowflake: Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

codecov[bot] avatar Oct 13 '25 18:10 codecov[bot]