When running a slightly edited version of the molecular_dna’s human_cell.mac example, I have been consistently running into a problem whereby the application stops running when it generates some of the molecular_dna threads (just gets stuck at generating one of them). This is happening regardless of the number of threads or the number of events being run. The executable file being built is the exact same as the one in the molecular_dna example and works for all macro examples provided out of the box.
I am currently running this Geant4 11.2.2 on a Windows 11 System with 8GB of RAM and
Intel Core i5-1035G1 CPU.
My macro file that I am attempting to run is:
/process/dna/e-SolvationSubType Meesungnoen2002
# print macro commands on screen and save run history for all runs
/control/verbose 1
/control/saveHistory
/chem/activate true
# setting up world with cells
/world/worldSize 1000 um
/cell/radiusSize 500 500 500 um
/scheduler/endTime 5.0 ns
/scheduler/maxNullTimeSteps 10000000
/dnageom/radicalKillDistance 9 nm
/dnageom/interactionDirectRange 3.5 angstrom
# DNA Geometry Model
/dnageom/definitionFile geometries/cube-centred-X-8.txt
/dnageom/placementVolume turn geometries/turned_solenoid_750_withHistone.txt
/dnageom/placementVolume turntwist geometries/turned_twisted_solenoid_750_withHistone.txt true
/dnageom/placementVolume straight geometries/straight_solenoid_750_withHistone.txt
# DNA Damage Model
/dnadamage/directDamageLower 5 eV
/dnadamage/directDamageUpper 37.5 eV
/dnadamage/indirectOHBaseChance 1.0
/dnadamage/indirectOHStrandChance 0.405
/dnadamage/inductionOHChance 0.00
/dnadamage/indirectHBaseChance 1.0
/dnadamage/indirectHStrandChance 0.0
/dnadamage/inductionHChance 0.00
/dnadamage/indirectEaqBaseChance 1.0
/dnadamage/indirectEaqStrandChance 0.0
/dnadamage/inductionEaqChance 0.00
# Add the dna to the cell
/chromosome/add cell ellipse 7005 2500 5300 112061.6 0 0 nm 0 0 0
/run/initialize
# source conditions, spherical output:
/gps/particle alpha
/gps/pos/type Volume
/gps/pos/shape Sphere
/gps/pos/centre 0 0 0 um
/gps/pos/radius 0.1 um
/gps/ang/type iso
# Gaussian energy distribution
/gps/ene/type Gauss
/gps/ene/mono 7.1 MeV
/gps/ene/sigma 0.43333 MeV
/run/beamOn 100
Is there anything that I am missing?