Download geant4.11.2

Getting Started
1

Hello!
I removed geant4.11.1 because I had problems related to zsh bus error. So, I download geant4 .11.2. I want to cope the example neuron to a new directory at the home but I want to check how to do CMake by the following lines.
1- ~%mkdir g4test

2- cd g4test

3- g4test%cp -r /Users/saraabdelaziz/Downloads/geant4-v11.2/examples/extended/medical/dna/neuron .

OR
% cp -r /usr/local/share/Geant4/examples/extended/medical/dna/neuron .
I want to know which one of them I should use, please?

4- g4test%ls

5- neuron

6- g4test % mkdir neuron-build

7- g4test%cd neuron-build

8- neuron-build% cmake … BUT in zsh terminal , it is not accepted …

SO should I replace with
% cmake …/neuron WHICH … refers to /Users/saraabdelaziz/g4test?
Could you help me, please?

2

Both are identical in terms of the code they copy, so either will do.

The correct invocation is cmake .., .i.e. two period symbols not three.

3

Oh, Yes. I just write it as to end the sentence. Thank you. I will try to use /usr/local/share/Geant4/examples/extended/medical/dna/neuron. So, when I do cmake, I will write cmake …
Because this path will be from (~) Home directory, will not it?

4

Also, in the .zshrc file, I write
export PATH=“/opt/homebrew/bin:$PATH”
When I try to write source /usr/local/bin/geant4.sh , It does not work.
ERROR: geant4.sh could NOT self-locate Geant4 installation
This is most likely because you are using ksh, dash or similar
To fix this issue, cd to the directory containing this script
and source it in that directory.

So, I try
sarahabdelaziz@Sarahs-MacBook-Pro ~ % source /usr/local/bin/geant4.sh

Is it OK? Should I start in the terminal to write source /usr/local/bin/geant4.sh
everytime?

5

If I follow where you’ve copied the neuron code to, then the operations will be:

% cd $HOME/g4test
% ls
neuron

% mkdir neuron-build
% cd neuron-build
% cmake ../neuron

I missed that you’d created your build directory alongside the source code.

With Z-shell, there’s a limitation on that script, so what you could put in your .zshrc is

pushd /usr/local/bin
source ./geant4.sh
popd

and that should be o.k. On MacOS you can probably get away with not even sourcing this script provided you installed Geant4 with -DGEANT4_INSTALL_DATA=ON.

6

Thanks. I did as you write.
and this is what I have.
sarahabdelaziz@Sarahs-MacBook-Pro ~ % nano .zshrc
sarahabdelaziz@Sarahs-MacBook-Pro ~ % source .zshrc
/usr/local/bin ~
~
sarahabdelaziz@Sarahs-MacBook-Pro ~ % cat .zshrc
export PATH=“/opt/homebrew/bin:$PATH”

pushd /usr/local/bin
source ./geant4.sh
popd

Is it OK? Should I write at everytime in the terminal ?
~% source /usr/local/bin/ geant4.sh

7

No, the lines in the .zshrc will suffice.

8

Please, find these videos, again I have the same problem about zsh: bus error.

I just try to modify in PhysicsList.cc.
Could you help me to know what is the problem, please?
Also, I compile B1 before I start. This is a video

9

Sorry, videos aren’t much help, could you please follow the guidelines: How to post code snippets

These also apply to console output like errors, and other problems.

10 
sarahabdelaziz@Sarahs-MacBook-Pro ~ % mkdir g4test
sarahabdelaziz@Sarahs-MacBook-Pro ~ % ls
Desktop		Library		Pictures	g4test
Documents	Movies		Public		try1.C
Downloads	Music		data1.txt
sarahabdelaziz@Sarahs-MacBook-Pro ~ % cd g4test 
sarahabdelaziz@Sarahs-MacBook-Pro g4test % cp -r /usr/local/share/Geant4/examples/extended/medical/dna/neuron .
sarahabdelaziz@Sarahs-MacBook-Pro g4test % ls
neuron
sarahabdelaziz@Sarahs-MacBook-Pro g4test % cd neuron 
sarahabdelaziz@Sarahs-MacBook-Pro neuron % ls
CMakeLists.txt		NeuralNETWORK.dat	neuron.in
Figs			README			neuron.out
GNUmakefile		gui.mac			plotDend.C
GranuleCell-Nr2.CNG.swc	include			src
History			neuron.cc		vis.mac
sarahabdelaziz@Sarahs-MacBook-Pro neuron % cd ..
sarahabdelaziz@Sarahs-MacBook-Pro g4test % mkdir neuron-build
sarahabdelaziz@Sarahs-MacBook-Pro g4test % cd neuron-build
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ccmake ../neuron


sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % cmake ../neuron 
-- Configuring done (0.2s)
-- Generating done (0.0s)
-- Build files have been written to: /Users/sarahabdelaziz/g4test/neuron-build
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ls
CMakeCache.txt		NeuralNETWORK.dat	neuron.in
CMakeFiles		README			plotDend.C
GranuleCell-Nr2.CNG.swc	cmake_install.cmake	vis.mac
Makefile		gui.mac
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % make
[  5%] Building CXX object CMakeFiles/neuron.dir/neuron.cc.o
[ 10%] Building CXX object CMakeFiles/neuron.dir/src/ActionInitialization.cc.o
[ 15%] Building CXX object CMakeFiles/neuron.dir/src/CommandLineParser.cc.o
[ 21%] Building CXX object CMakeFiles/neuron.dir/src/DetectorConstruction.cc.o
[ 26%] Building CXX object CMakeFiles/neuron.dir/src/EventAction.cc.o
[ 31%] Building CXX object CMakeFiles/neuron.dir/src/ITSteppingAction.cc.o
[ 36%] Building CXX object CMakeFiles/neuron.dir/src/ITTrackingAction.cc.o
[ 42%] Building CXX object CMakeFiles/neuron.dir/src/ITTrackingInteractivity.cc.o
[ 47%] Building CXX object CMakeFiles/neuron.dir/src/NeuronLoadDataFile.cc.o
[ 52%] Building CXX object CMakeFiles/neuron.dir/src/PhysicsList.cc.o
[ 57%] Building CXX object CMakeFiles/neuron.dir/src/PrimaryGeneratorAction.cc.o
[ 63%] Building CXX object CMakeFiles/neuron.dir/src/Run.cc.o
[ 68%] Building CXX object CMakeFiles/neuron.dir/src/RunAction.cc.o
[ 73%] Building CXX object CMakeFiles/neuron.dir/src/RunInitObserver.cc.o
[ 78%] Building CXX object CMakeFiles/neuron.dir/src/StackingAction.cc.o
[ 84%] Building CXX object CMakeFiles/neuron.dir/src/SteppingAction.cc.o
[ 89%] Building CXX object CMakeFiles/neuron.dir/src/TimeStepAction.cc.o
[ 94%] Building CXX object CMakeFiles/neuron.dir/src/TrackingAction.cc.o
[100%] Linking CXX executable neuron
[100%] Built target neuron
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ls
CMakeCache.txt		NeuralNETWORK.dat	neuron
CMakeFiles		README			neuron.in
GranuleCell-Nr2.CNG.swc	cmake_install.cmake	plotDend.C
Makefile		gui.mac			vis.mac
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ./neuron -mac neuron.in

**************************************************************
Geant4 version Name: geant4-11-02-patch-02 [MT]   (21-June-2024)
 << in Multi-threaded mode >> 
                      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/
**************************************************************

NeuronLoadDataFile::SingleNeuronSWCfile: number of strings=2922 in file GranuleCell-Nr2.CNG.swc
Total number of compartments into Neuron : 2914

Livermore + DNAphysics is activated!
---- Begin of Neuron Construction! ---- 
==========================================================
Side length of word volume is calculated : 407.347 um
Radius of homogeneous medium is calculated : 203.674 um
Overall dimensions (um) of neuron morphology : 
    width = 250.225 height = 315.035 depth = 63.8148
Volume (um3), surface (um2) and mass (ug) of Bounding Slice are calculated : 
   5.03049e+06; 	229803; 	5.03049

Volume (um3), surface (um2) and mass(ug) of Neuron are calculated : 
   12100; 	30379.7; 	0.0121
Total number of compartments into Neuron : 2914
Shift values (um) for Neuron translation are calculated : 
    shiftX = 35.13 shiftY = 99.995 shiftZ = 2.025
---- Soma for construction: ---- 
Total number of compartments into Soma : 3
---- Dendrites for construction: ---- 
Total number of compartments into Dendrites: 2775
---- Axon for construction: ---- 
Total number of compartments into Axon : 136
---- Spines not found! ---- 

---- End of Neuron Construction! ---- 
======================================
hPairProd:  for pi-  XStype:1  SubType=4
     dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
     Lambda tables from threshold to 100 TeV, 20 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

=========================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Initializing workers...
=========================================================================

Visualization Manager instantiating with verbosity "warnings (3)"...
Visualization Manager initialising...
Registering graphics systems...

You have successfully registered the following graphics systems.
Registered graphics systems are:
 ASCIITree (ATree)
 DAWNFILE (DAWNFILE)
 G4HepRepFile (HepRepFile)
 RayTracer (RayTracer)
 VRML2FILE (VRML2FILE)
 gMocrenFile (gMocrenFile)
 TOOLSSG_OFFSCREEN (TSG_OFFSCREEN)
 TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
 OpenGLImmediateQt (OGLIQt, OGLI)
 OpenGLStoredQt (OGLSQt, OGL, OGLS)
 OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK)
 OpenGLStoredX (OGLSX, OGLSQt_FALLBACK)
 Qt3D (Qt3D)
 TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_QT_GLES_FALLBACK)
 TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
 TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
Default graphics system is: TSG_OFFSCREEN (based on batch session).
Default window size hint is: 600x600-0+0 (based on G4VisManager initialisation).
Note: Parameters specified on the command line will override these defaults.
     Use "vis/open" without parameters to get these defaults.

Registering model factories...

You have successfully registered the following model factories.
Registered model factories:
 generic
 drawByAttribute
 drawByCharge
 drawByOriginVolume
 drawByParticleID
 drawByEncounteredVolume

Registered models:
 None

Registered filter factories:
 attributeFilter
 chargeFilter
 originVolumeFilter
 particleFilter
 encounteredVolumeFilter

Registered filters:
 None

You have successfully registered the following user vis actions.
Run Duration User Vis Actions: none
End of Event User Vis Actions: none
End of Run User Vis Actions: none

Some /vis commands (optionally) take a string to specify colour.
"/vis/list" to see available colours.

### ===  Deexcitation model UAtomDeexcitation is activated for 4 regions:
         DefaultRegionForTheWorld  1  0  0
         Soma  1  0  0
         Dendrites  1  0  0
         Axon  1  0  0
### ===  Auger flag: 1
### ===  Ignore cuts flag:   1
G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model.
See commands in /vis/modeling/trajectories/ for other options.

========================================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Creating 2 tasks with 1 events/task...
========================================================================================

G4WT2 > End of tracking primary particle, its final energy is :3.4748 GeV
G4WT0 > End of tracking primary particle, its final energy is :3.47513 GeV
G4WT0 > Physics stage ends
G4WT2 > Physics stage ends

======= The summary of simulation results 'neuron' ========

 Primary particle               = C12
 Kinetic energy of beam         = 3480 A*MeV
 Particle traversals the neuron = 2 of 2
 Full LET of beam as formulas   = 12.827 keV/um
 Mean LET of beam as simulation = 11.356 +- 0.45951 keV/um
 Mean track length of beam      = 418.39 +- 0.033101 um
 Particle fluence               = 383.66 particles/cm^2
 Fluence dose (full)            = 7.8848e-06 Gy
 Fluence dose ber beam          = 3.9424e-06 Gy

List of generated physical process:
               eIoni=   4653
             ionIoni=   2562
                 msc=   7083


List of generated particles outside neuron structure:
            e- :    3339  Emean =  1.6691 keV	( 242.22 eV  --> 333.82 keV) 	Eflow/event = 2.7866 MeV

Number of secondary particles inside neuron structure:
            e- :      15

Number of molecular products inside neuron structure:
   time: 1 ps - 10 ps 

Total energy (MeV) deposition :
 All volume:   9.5028
  Bounding slice:  4.2848
     Neuron:   0.0099629
       Soma:   0.0098267
   Dendrites:  0.00013624
       Axon:         0

Number of compartments traversed by particle tracks :
       Soma:        1 of total: 3
    Dendrites:       1 of total: 2775
         Axon:       0 of total: 136


Dendritic (or Axon) compartmental energy deposits 
at the distance from Soma have been written into *.out files:
Dend3DEdep.out, Axon3DEdep.out, Soma3DEdep.out
Outputs of energy deposition per event written in data file:
OutputPerEvent.out
Calculation time                      = 0.960133 s 

Graphics systems deleted.
Visualization Manager deleting...

All of this work good.
Then, I try to modify in PhysicsList.cc of g4test/neuron/src

 // a) DNAphysics and Livermore physics inside and outside neuron 
  G4cout<< "Livermore + DNAphysics is activated!"<<G4endl;
 // RegisterConstructor("emlivermore");
    RegisterConstructor("G4EmDNAPhysics_option2") ;
  // 'G4EmParameters' works together with 'G4EmDNAPhysicsActivator'
  // VI: in this example Livermore is a default in any way
  //     DNA option4 is the default configuration as well
  fDNAActivator = new G4EmDNAPhysicsActivator();

Instead of RegisterConstructor(“emlivermore”); → RegisterConstructor(“G4EmDNAPhysics_option2”) ;
Then , In the terminal

sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % make
[  5%] Building CXX object CMakeFiles/neuron.dir/src/PhysicsList.cc.o
[ 10%] Linking CXX executable neuron
[100%] Built target neuron
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % nano neuron.in

I modify in neuron.in for neuron-build to display only proton 30 MeV and comment others

/tracking/verbose 0
/control/verbose 0
/run/verbose 0
/process/eLoss/verbose 1
/process/em/fluo true
/process/em/auger true
/process/em/augerCascade true
/process/em/deexcitationIgnoreCut true
# Particle direction and position
#/gun/direction 0 0 1
#/gun/position 0 0 0 um

# Particle name and energy
/gun/particle proton
/gun/energy 30 MeV
#/gun/particle ion
#/gun/ion 6 12 6
#/gun/energy 3480 MeV

#/particle/property/dump
#/particle/list
#/process/list
#/particle/process/dump
#/run/printProgress 100

/run/initialize

/run/beamOn 2

Then, in the terminal

sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ./neuron -mac neuron.in -mt 2 -dnachemON -network NeuralNETWORK.dat 

**************************************************************
 Geant4 version Name: geant4-11-02-patch-02 [MT]   (21-June-2024)
  << in Multi-threaded mode >> 
                       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/
**************************************************************

===== neuron is started with 2 threads of MT MODE =====
NeuronLoadDataFile::NeuralNetworkDATAfile: opened NeuralNETWORK.dat
 Total number of compartments into Neuron : 18376
Livermore + DNAphysics is activated!
DNAChemistry is activated!
 ---- Begin of Neuron Construction! ---- 
 ==========================================================
 Side length of word volume is calculated : 699.017 um
 Radius of homogeneous medium is calculated : 349.509 um
 Overall dimensions (um) of neuron morphology : 
	 width = 640 height = 280 depth = 25
 Volume (um3), surface (um2) and mass (ug) of Bounding Slice are calculated : 
	4.48e+06; 	404400; 	4.48
 
 Volume (um3), surface (um2) and mass(ug) of Neuron are calculated : 
	100726; 	140835; 	0.100726
 Total number of compartments into Neuron : 18376
 Shift values (um) for Neuron translation are calculated : 
	 shiftX = 0 shiftY = 0 shiftZ = 0
 ---- Soma for construction: ---- 
 Total number of compartments into Soma : 96
 ---- Dendrites for construction: ---- 
 Total number of compartments into Dendrites: 18279
 ---- Axon not found! ---- 
 ---- Spines not found! ---- 

 ---- End of Neuron Construction! ---- 
 ========================================================== 

=======================================================================
======                 Electromagnetic Physics Parameters      ========
============================================================

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

=========================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Initializing workers...
=========================================================================

Visualization Manager instantiating with verbosity "warnings (3)"...
Visualization Manager initialising...
Registering graphics systems...

You have successfully registered the following graphics systems.
Registered graphics systems are:
  ASCIITree (ATree)
  DAWNFILE (DAWNFILE)
  G4HepRepFile (HepRepFile)
  RayTracer (RayTracer)
  VRML2FILE (VRML2FILE)
  gMocrenFile (gMocrenFile)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
  OpenGLImmediateQt (OGLIQt, OGLI)
  OpenGLStoredQt (OGLSQt, OGL, OGLS)
  OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK)
  OpenGLStoredX (OGLSX, OGLSQt_FALLBACK)
  Qt3D (Qt3D)
  TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_QT_GLES_FALLBACK)
  TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
  TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
Default graphics system is: TSG_OFFSCREEN (based on batch session).
Default window size hint is: 600x600-0+0 (based on G4VisManager initialisation).
Note: Parameters specified on the command line will override these defaults.
      Use "vis/open" without parameters to get these defaults.

Registering model factories...

You have successfully registered the following model factories.
Registered model factories:
  generic
  drawByAttribute
  drawByCharge
  drawByOriginVolume
  drawByParticleID
  drawByEncounteredVolume

Registered models:
  None

Registered filter factories:
  attributeFilter
  chargeFilter
  originVolumeFilter
  particleFilter
  encounteredVolumeFilter

Registered filters:
  None

You have successfully registered the following user vis actions.
Run Duration User Vis Actions: none
End of Event User Vis Actions: none
End of Run User Vis Actions: none

Some /vis commands (optionally) take a string to specify colour.
"/vis/list" to see available colours.
=======================================================================
======                 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                                0
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                      10 eV 
Max kinetic energy for tables                      300 MeV
Number of bins per decade of a table               20
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, 0.01 mm)
Step function for muons/hadrons                    (0.1, 0.05 mm)
Step function for light ions                       (0.1, 0.02 mm)
Step function for general ions                     (0.1, 0.001 mm)
Lowest e+e- kinetic energy                         0 eV 
Lowest muon/hadron kinetic energy                  1 keV
Use ICRU90 data                                    1
Fluctuations of dE/dx are enabled                  1
Type of fluctuation model for leptons and hadrons  Universal
Use built-in Birks satuaration                     0
Build CSDA range enabled                           0
Use cut as a final range enabled                   0
Enable angular generator interface                 1
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+-           2
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.08
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+-      3
Lambda limit for msc step limit for e+-            1 mm
Use Mott correction for e- scattering              1
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
=======================================================================
======                 Atomic Deexcitation Parameters          ========
=======================================================================
Fluorescence enabled                               1
Directory in G4LEDATA for fluorescence data files  fluor
Auger electron cascade enabled                     1
PIXE atomic de-excitation enabled                  0
De-excitation module ignores cuts                  1
Type of PIXE cross section for hadrons             Empirical
Type of PIXE cross section for e+-                 Livermore
=======================================================================
======                 DNA Physics Parameters                  ========
=======================================================================
Use fast sampling in DNA models                    1
Use Stationary option in DNA models                0
Use DNA with multiple scattering of e-             0
Use DNA e- solvation model type                    11003
=======================================================================

### ===  Deexcitation model UAtomDeexcitation is activated for 3 regions:
          DefaultRegionForTheWorld  1  1  0
          Soma  1  1  0
          Dendrites  1  1  0
### ===  Auger flag: 1
### ===  Ignore cuts flag:   1

---
proton_G4DNAElastic:  for proton SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

proton_G4DNAExcitation:  for proton SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  500 keV
DNABornExcitationModel : Emin=  500 keV Emax=  100 MeV
DNARPWBAExcitationModel : Emin=  100 MeV Emax=  300 MeV

proton_G4DNAIonisation:  for proton SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNARuddIonisationExtendedModel : Emin=    0 eV  Emax=  500 keV  deltaRudd Fluo
DNABornIonisationModel : Emin=  500 keV Emax=  100 MeV  deltaBorn Fluo
DNARPWBAIonisationModel : Emin=  100 MeV Emax=  300 MeV  deltaBorn Fluo

proton_G4DNAChargeDecrease:  for proton SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeDecreaseModel : Emin=    0 eV  Emax=  100 MeV
          DummyModel : Emin=  100 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

GenericIon_G4DNAIonisation:  for GenericIon SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

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

hIoni:  for pi-  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=297.505 keV  deltaVI
          BetheBloch : Emin=297.505 keV Emax=  300 MeV  deltaVI
G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model.
See commands in /vis/modeling/trajectories/ for other options.

========================================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Creating 2 tasks with 1 events/task...
========================================================================================

G4WT1 > End of tracking primary particle, its final energy is :28.7595 MeV
G4WT0 > End of tracking primary particle, its final energy is :28.6105 MeV
G4WT1 > Physics stage ends
G4WT1 > *** G4Scheduler starts processing 
G4WT0 > Physics stage ends
G4WT0 > *** G4Scheduler starts processing 
zsh: bus error  ./neuron -mac neuron.in -mt 2 -dnachemON -network NeuralNETWORK.dat
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ls

As at the end, there is error bus error? Also, G4Scheduler does not start.

11

Although, during work with RegisterConstructor(“emlivermore”); which in PhysicsList.cc of g4test/neuron/src. I do not have this problem!

sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % make
[  5%] Building CXX object CMakeFiles/neuron.dir/src/PhysicsList.cc.o
[ 10%] Linking CXX executable neuron
[100%] Built target neuron
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % nano neuron.in
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % ./neuron -mac neuron.in -mt 2 -dnachemON -network NeuralNETWORK.dat

**************************************************************
 Geant4 version Name: geant4-11-02-patch-02 [MT]   (21-June-2024)
  << in Multi-threaded mode >> 
                       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/
**************************************************************

===== neuron is started with 2 threads of MT MODE =====
NeuronLoadDataFile::NeuralNetworkDATAfile: opened NeuralNETWORK.dat
 Total number of compartments into Neuron : 18376
Livermore + DNAphysics is activated!
DNAChemistry is activated!
 ---- Begin of Neuron Construction! ---- 
 ==========================================================
 Side length of word volume is calculated : 699.017 um
 Radius of homogeneous medium is calculated : 349.509 um
 Overall dimensions (um) of neuron morphology : 
	 width = 640 height = 280 depth = 25
 Volume (um3), surface (um2) and mass (ug) of Bounding Slice are calculated : 
	4.48e+06; 	404400; 	4.48
 
 Volume (um3), surface (um2) and mass(ug) of Neuron are calculated : 
	100726; 	140835; 	0.100726
 Total number of compartments into Neuron : 18376
 Shift values (um) for Neuron translation are calculated : 
	 shiftX = 0 shiftY = 0 shiftZ = 0
 ---- Soma for construction: ---- 
 Total number of compartments into Soma : 96
 ---- Dendrites for construction: ---- 
 Total number of compartments into Dendrites: 18279
 ---- Axon not found! ---- 
 ---- Spines not found! ---- 

 ---- End of Neuron Construction! ---- 
 ========================================================== 

=======================================================================
======                 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                                0
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               20
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
---
=======================================================================
======                 Atomic Deexcitation Parameters          ========
=======================================================================
Fluorescence enabled                               1
Directory in G4LEDATA for fluorescence data files  fluor
Auger electron cascade enabled                     0
PIXE atomic de-excitation enabled                  0
De-excitation module ignores cuts                  0
Type of PIXE cross section for hadrons             Empirical
Type of PIXE cross section for e+-                 Livermore
=======================================================================
======                 DNA Physics Parameters                  ========
=======================================================================
Use fast sampling in DNA models                    0
Use Stationary option in DNA models                0
Use DNA with multiple scattering of e-             0
Use DNA e- solvation model type                    11003
=======================================================================

### ===  Deexcitation model UAtomDeexcitation is activated for 3 regions:
          DefaultRegionForTheWorld  1  0  0
          Soma  1  0  0
          Dendrites  1  0  0
### ===  Ignore cuts flag:   0

phot:  for gamma SubType=12 BuildTable=0
      LambdaPrime table from 200 keV to 100 TeV in 174 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, 20 bins/decade, spline: 1
      LambdaPrime table from 1 MeV to 100 TeV in 160 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
    LivermoreCompton : Emin=    0 eV  Emax=    1 GeV Fluo
        KleinNishina : Emin=    1 GeV Emax=  100 TeV Fluo

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

Rayl:  for gamma SubType=11 BuildTable=1
      Lambda table from 100 eV  to 150 keV, 20 bins/decade, spline: 0
      LambdaPrime table from 150 keV to 100 TeV in 176 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh : Emin=    0 eV  Emax=  100 TeV  CullenGenerator

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

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

hBrems:  for proton  XStype:1  SubType=3
      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
      Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               hBrem : Emin=    0 eV  Emax=  100 TeV  ModifiedMephi
      ===== Limit on energy threshold has been applied 

hPairProd:  for proton  XStype:1  SubType=4
      dE/dx and range tables from 100 eV  to 100 TeV in 240 bins
      Lambda tables from threshold to 100 TeV, 20 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, 20 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

nuclearStopping:  for proton SubType=8 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   ICRU49NucStopping : Emin=    0 eV  Emax=    1 MeV

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:1 Skin=3 Llim=1 mm

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

----
    ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  eCoulombScattering : Emin=    0 eV  Emax=  100 TeV

=========================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Initializing workers...
=========================================================================

Visualization Manager instantiating with verbosity "warnings (3)"...
Visualization Manager initialising...
Registering graphics systems...

You have successfully registered the following graphics systems.
Registered graphics systems are:
  ASCIITree (ATree)
  DAWNFILE (DAWNFILE)
  G4HepRepFile (HepRepFile)
  RayTracer (RayTracer)
  VRML2FILE (VRML2FILE)
  gMocrenFile (gMocrenFile)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
  OpenGLImmediateQt (OGLIQt, OGLI)
  OpenGLStoredQt (OGLSQt, OGL, OGLS)
  OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK)
  OpenGLStoredX (OGLSX, OGLSQt_FALLBACK)
  Qt3D (Qt3D)
  TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_QT_GLES_FALLBACK)
  TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
  TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
Default graphics system is: TSG_OFFSCREEN (based on batch session).
Default window size hint is: 600x600-0+0 (based on G4VisManager initialisation).
Note: Parameters specified on the command line will override these defaults.
      Use "vis/open" without parameters to get these defaults.

Registering model factories...

You have successfully registered the following model factories.
Registered model factories:
  generic
  drawByAttribute
  drawByCharge
  drawByOriginVolume
  drawByParticleID
  drawByEncounteredVolume

Registered models:
  None

Registered filter factories:
  attributeFilter
  chargeFilter
  originVolumeFilter
  particleFilter
  encounteredVolumeFilter

Registered filters:
  None

You have successfully registered the following user vis actions.
Run Duration User Vis Actions: none
End of Event User Vis Actions: none
End of Run User Vis Actions: none

Some /vis commands (optionally) take a string to specify colour.
"/vis/list" to see available colours.

### ===  Deexcitation model UAtomDeexcitation is activated for 3 regions:
          DefaultRegionForTheWorld  1  0  0
          Soma  1  0  0
          Dendrites  1  0  0
### ===  Auger flag: 1
### ===  Ignore cuts flag:   1
G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model.
See commands in /vis/modeling/trajectories/ for other options.

========================================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Creating 2 tasks with 1 events/task...
========================================================================================

G4WT0 > End of tracking primary particle, its final energy is :28.7285 MeV
G4WT1 > End of tracking primary particle, its final energy is :28.7947 MeV
G4WT0 > Physics stage ends
G4WT0 > *** G4Scheduler starts processing 
G4WT0 > *** G4Scheduler did not start because no track was found to be processed
G4WT0 > ___________________________________
G4WT1 > Physics stage ends
G4WT1 > *** G4Scheduler starts processing 
G4WT1 > *** G4Scheduler did not start because no track was found to be processed
G4WT1 > ___________________________________

 ======= The summary of simulation results 'neuron' ========

  Primary particle               = proton
  Kinetic energy of beam         = 30 A*MeV
  Particle traversals the neuron = 1 of 2
  Full LET of beam as formulas   = 1.864 keV/um
  Mean LET of beam as simulation = 1.6897 +- 0.017192 keV/um
  Mean track length of beam      = 717.46 +- 0.00057621 um
  Particle fluence               = 130.29 particles/cm^2
  Fluence dose (full)            = 3.891e-07 Gy
  Fluence dose ber beam          = 1.9455e-07 Gy

 List of generated physical process:
                eIoni=   1182
                hIoni=    942
                  msc=   1376


 List of generated particles outside neuron structure:
             e- :     961  Emean =  1.1696 keV	( 242.19 eV  --> 58.006 keV) 	Eflow/event = 562.02 keV

 Number of secondary particles inside neuron structure:
             e- :       1

 Number of molecular products inside neuron structure:
    time: 1 ps - 10 ps 

 Total energy (MeV) deposition :
  All volume:   2.4246
   Bounding slice: 0.48238
      Neuron:   0.00012086
        Soma:         0
    Dendrites:  0.00012086
        Axon:         0

 Number of compartments traversed by particle tracks :
        Soma:        0 of total: 96
     Dendrites:       1 of total: 18279
          Axon:       0 of total: 0
 

 Dendritic (or Axon) compartmental energy deposits 
 at the distance from Soma have been written into *.out files:
 Dend3DEdep.out, Axon3DEdep.out, Soma3DEdep.out
 Outputs of energy deposition per event written in data file:
 OutputPerEvent.out
 Calculation time                      = 1.33618 s 

Graphics systems deleted.
Visualization Manager deleting...
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build %

Also, I noticed that DNA Physics Parameters as shown above is the same for two cases. Dose it mean there is a problem related to EMDNAphysics model process?

12

My Mac is Sonoma 14.6.1, 500 GB, chip M2. Does it not work well with geant4? Should I download VMware with Geant4 on mac? But, I am afraid because VMware has a limit space to storage data and output files. Also, I want to work for energy range from 30 10 300 MeV. So, could anyone advice me if I work with VMware, what the space I need? . I try to work on geant4.11.1 and geant4.11.2 but I have the same problem related to EMDNAPhysics processes.

13

I’m not an expert on this particular example (perhaps @ivana and @guatelli can help), but this is probably indicative of a memory or similar issue caused by the modification you’ve made. The easiest way to start to get a handle on that is to run through a debugger, which in macOS is lldb. Try running as:

lldb -- ./neuron -mac neuron.in -mt 2 -dnachemON -network NeuralNETWORK.dat

then type run when it’s loaded and you get a (lldb) prompt. It should then execute and stop if something goes wrong. You can usually type bt at that point to get a traceback of where in the code the problem occurred.

Geant4 works fine to the best of our testing and experience on Sonoma. A VM in this case won’t help in any way and will simply complicate matters.

14

[/quote]

Thanks for your help. I did as you say and this is what I have …

/usr/local/bin ~/g4test/neuron-build
~/g4test/neuron-build
sarahabdelaziz@Sarahs-MacBook-Pro neuron-build % lldb -- ./neuron -mac neuron.in -mt 2 -dnachemON -network NeuralNETWORK.dat
(lldb) target create "./neuron"
Current executable set to '/Users/sarahabdelaziz/g4test/neuron-build/neuron' (arm64).
(lldb) settings set -- target.run-args  "-mac" "neuron.in" "-mt" "2" "-dnachemON" "-network" "NeuralNETWORK.dat"
(lldb) run
Process 9415 launched: '/Users/sarahabdelaziz/g4test/neuron-build/neuron' (arm64)

**************************************************************
 Geant4 version Name: geant4-11-02-patch-02 [MT]   (21-June-2024)
  << in Multi-threaded mode >> 
                       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/
**************************************************************

===== neuron is started with 2 threads of MT MODE =====
NeuronLoadDataFile::NeuralNetworkDATAfile: opened NeuralNETWORK.dat
 Total number of compartments into Neuron : 18376
Livermore + DNAphysics is activated!
DNAChemistry is activated!
 ---- Begin of Neuron Construction! ---- 
 ==========================================================
 Side length of word volume is calculated : 699.017 um
 Radius of homogeneous medium is calculated : 349.509 um
 Overall dimensions (um) of neuron morphology : 
	 width = 640 height = 280 depth = 25
 Volume (um3), surface (um2) and mass (ug) of Bounding Slice are calculated : 
	4.48e+06; 	404400; 	4.48
 
 Volume (um3), surface (um2) and mass(ug) of Neuron are calculated : 
	100726; 	140835; 	0.100726
 Total number of compartments into Neuron : 18376
 Shift values (um) for Neuron translation are calculated : 
	 shiftX = 0 shiftY = 0 shiftZ = 0
 ---- Soma for construction: ---- 
 Total number of compartments into Soma : 96
 ---- Dendrites for construction: ---- 
 Total number of compartments into Dendrites: 18279
 ---- Axon not found! ---- 
 ---- Spines not found! ---- 

 ---- End of Neuron Construction! ---- 
 ========================================================== 

=======================================================================
======                 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                                0
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                      10 eV 
Max kinetic energy for tables                      300 MeV
Number of bins per decade of a table               20
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, 0.01 mm)
Step function for muons/hadrons                    (0.1, 0.05 mm)
Step function for light ions                       (0.1, 0.02 mm)
Step function for general ions                     (0.1, 0.001 mm)
Lowest e+e- kinetic energy                         0 eV 
Lowest muon/hadron kinetic energy                  1 keV
Use ICRU90 data                                    1
Fluctuations of dE/dx are enabled                  1
Type of fluctuation model for leptons and hadrons  Universal
Use built-in Birks satuaration                     0
Build CSDA range enabled                           0
Use cut as a final range enabled                   0
Enable angular generator interface                 1
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+-           2
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.08
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+-      3
Lambda limit for msc step limit for e+-            1 mm
Use Mott correction for e- scattering              1
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
=======================================================================
======                 Atomic Deexcitation Parameters          ========
=======================================================================
Fluorescence enabled                               1
Directory in G4LEDATA for fluorescence data files  fluor
Auger electron cascade enabled                     1
PIXE atomic de-excitation enabled                  0
De-excitation module ignores cuts                  1
Type of PIXE cross section for hadrons             Empirical
Type of PIXE cross section for e+-                 Livermore
=======================================================================
======                 DNA Physics Parameters                  ========
=======================================================================
Use fast sampling in DNA models                    1
Use Stationary option in DNA models                0
Use DNA with multiple scattering of e-             0
Use DNA e- solvation model type                    11003
=======================================================================

### ===  Deexcitation model UAtomDeexcitation is activated for 3 regions:
          DefaultRegionForTheWorld  1  1  0
          Soma  1  1  0
          Dendrites  1  1  0
### ===  Auger flag: 1
### ===  Ignore cuts flag:   1

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

compt:  for gamma SubType=13 BuildTable=1
      Lambda table from 10 eV  to 1 MeV, 19 bins/decade, spline: 1
      LambdaPrime table from 1 MeV to 300 MeV in 46 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LowEPComptonModel : Emin=    0 eV  Emax=   20 MeV Fluo
        KleinNishina : Emin=   20 MeV Emax=  300 MeV Fluo

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

Rayl:  for gamma SubType=11 BuildTable=1
      Lambda table from 10 eV  to 150 keV, 19 bins/decade, spline: 0
      LambdaPrime table from 150 keV to 300 MeV in 62 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh : Emin=    0 eV  Emax=  300 MeV  CullenGenerator

msc:  for e-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  GoudsmitSaunderson : Emin=    1 MeV Emax=  300 MeV Nbins=40   1 MeV - 300 MeV
          StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm

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

eBrem:  for e-  XStype:4  SubType=3
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      LPM flag: 1 for E > 0.3 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    1 MeV Emax=  300 MeV  AngularGen2BS

e-_G4DNAElectronSolvation:  for e- SubType=58 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAOneStepThermalizationModel_Meesungnoen2002 : Emin=    0 eV  Emax=  7.4 eV 
          DummyModel : Emin=  7.4 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAElastic:  for e- SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAChampionElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAExcitation:  for e- SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNABornExcitationModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAIonisation:  for e- SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNABornIonisationModel : Emin=    0 eV  Emax=    1 MeV  deltaBorn Fluo
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAVibExcitation:  for e- SubType=54 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNASancheExcitationModel : Emin=    0 eV  Emax=  100 eV 
          DummyModel : Emin=  100 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAAttachment:  for e- SubType=55 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMeltonAttachmentModel : Emin=    0 eV  Emax=   13 eV 
          DummyModel : Emin=   13 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

msc:  for e+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  GoudsmitSaunderson : Emin=    0 eV  Emax=  300 MeV Nbins=120 100 eV  - 300 MeV
          StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm

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

eBrem:  for e+  XStype:4  SubType=3
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      LPM flag: 1 for E > 0.3 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 eV  Emax=  300 MeV  AngularGen2BS

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

proton_G4DNAElastic:  for proton SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

proton_G4DNAExcitation:  for proton SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  500 keV
DNABornExcitationModel : Emin=  500 keV Emax=  100 MeV
DNARPWBAExcitationModel : Emin=  100 MeV Emax=  300 MeV

proton_G4DNAIonisation:  for proton SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNARuddIonisationExtendedModel : Emin=    0 eV  Emax=  500 keV  deltaRudd Fluo
DNABornIonisationModel : Emin=  500 keV Emax=  100 MeV  deltaBorn Fluo
DNARPWBAIonisationModel : Emin=  100 MeV Emax=  300 MeV  deltaBorn Fluo

proton_G4DNAChargeDecrease:  for proton SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeDecreaseModel : Emin=    0 eV  Emax=  300 MeV

GenericIon_G4DNAIonisation:  for GenericIon SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha_G4DNAElastic:  for alpha SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

alpha_G4DNAExcitation:  for alpha SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  300 MeV

alpha_G4DNAIonisation:  for alpha SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha_G4DNAChargeDecrease:  for alpha SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeDecreaseModel : Emin=    0 eV  Emax=  300 MeV

alpha+_G4DNAElastic:  for alpha+ SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

alpha+_G4DNAExcitation:  for alpha+ SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  300 MeV

alpha+_G4DNAIonisation:  for alpha+ SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha+_G4DNAChargeIncrease:  for alpha+ SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeIncreaseModel : Emin=    0 eV  Emax=  300 MeV

alpha+_G4DNAChargeDecrease:  for alpha+ SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeDecreaseModel : Emin=    0 eV  Emax=  300 MeV

msc:  for anti_proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

helium_G4DNAElastic:  for helium SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

helium_G4DNAExcitation:  for helium SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  300 MeV

helium_G4DNAIonisation:  for helium SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

helium_G4DNAChargeIncrease:  for helium SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeIncreaseModel : Emin=    0 eV  Emax=  300 MeV

hydrogen_G4DNAElastic:  for hydrogen SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

hydrogen_G4DNAExcitation:  for hydrogen SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  300 MeV

hydrogen_G4DNAIonisation:  for hydrogen SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

hydrogen_G4DNAChargeIncrease:  for hydrogen SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeIncreaseModel : Emin=    0 eV  Emax=  300 MeV

msc:  for kaon+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

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

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

msc:  for mu+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

hIoni:  for mu+  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=225.219 keV  deltaVI
          BetheBloch : Emin=225.219 keV Emax=  300 MeV  deltaVI

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

hIoni:  for mu-  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=225.219 keV  deltaVI
          BetheBloch : Emin=225.219 keV Emax=  300 MeV  deltaVI

msc:  for pi+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

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

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

=========================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Initializing workers...
=========================================================================

Visualization Manager instantiating with verbosity "warnings (3)"...
Visualization Manager initialising...
Registering graphics systems...

You have successfully registered the following graphics systems.
Registered graphics systems are:
  ASCIITree (ATree)
  DAWNFILE (DAWNFILE)
  G4HepRepFile (HepRepFile)
  RayTracer (RayTracer)
  VRML2FILE (VRML2FILE)
  gMocrenFile (gMocrenFile)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN)
  TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
  OpenGLImmediateQt (OGLIQt, OGLI)
  OpenGLStoredQt (OGLSQt, OGL, OGLS)
  OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK)
  OpenGLStoredX (OGLSX, OGLSQt_FALLBACK)
  Qt3D (Qt3D)
  TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_QT_GLES_FALLBACK)
  TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG)
  TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
Default graphics system is: TSG_OFFSCREEN (based on batch session).
Default window size hint is: 600x600-0+0 (based on G4VisManager initialisation).
Note: Parameters specified on the command line will override these defaults.
      Use "vis/open" without parameters to get these defaults.

Registering model factories...

You have successfully registered the following model factories.
Registered model factories:
  generic
  drawByAttribute
  drawByCharge
  drawByOriginVolume
  drawByParticleID
  drawByEncounteredVolume

Registered models:
  None

Registered filter factories:
  attributeFilter
  chargeFilter
  originVolumeFilter
  particleFilter
  encounteredVolumeFilter

Registered filters:
  None

You have successfully registered the following user vis actions.
Run Duration User Vis Actions: none
End of Event User Vis Actions: none
End of Run User Vis Actions: none

Some /vis commands (optionally) take a string to specify colour.
"/vis/list" to see available colours.
=======================================================================
======                 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                                0
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                      10 eV 
Max kinetic energy for tables                      300 MeV
Number of bins per decade of a table               20
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, 0.01 mm)
Step function for muons/hadrons                    (0.1, 0.05 mm)
Step function for light ions                       (0.1, 0.02 mm)
Step function for general ions                     (0.1, 0.001 mm)
Lowest e+e- kinetic energy                         0 eV 
Lowest muon/hadron kinetic energy                  1 keV
Use ICRU90 data                                    1
Fluctuations of dE/dx are enabled                  1
Type of fluctuation model for leptons and hadrons  Universal
Use built-in Birks satuaration                     0
Build CSDA range enabled                           0
Use cut as a final range enabled                   0
Enable angular generator interface                 1
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+-           2
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.08
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+-      3
Lambda limit for msc step limit for e+-            1 mm
Use Mott correction for e- scattering              1
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
=======================================================================
======                 Atomic Deexcitation Parameters          ========
=======================================================================
Fluorescence enabled                               1
Directory in G4LEDATA for fluorescence data files  fluor
Auger electron cascade enabled                     1
PIXE atomic de-excitation enabled                  0
De-excitation module ignores cuts                  1
Type of PIXE cross section for hadrons             Empirical
Type of PIXE cross section for e+-                 Livermore
=======================================================================
======                 DNA Physics Parameters                  ========
=======================================================================
Use fast sampling in DNA models                    1
Use Stationary option in DNA models                0
Use DNA with multiple scattering of e-             0
Use DNA e- solvation model type                    11003
=======================================================================
15

The following code …

### ===  Deexcitation model UAtomDeexcitation is activated for 3 regions:
          DefaultRegionForTheWorld  1  1  0
          Soma  1  1  0
          Dendrites  1  1  0
### ===  Auger flag: 1
### ===  Ignore cuts flag:   1

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

compt:  for gamma SubType=13 BuildTable=1
      Lambda table from 10 eV  to 1 MeV, 19 bins/decade, spline: 1
      LambdaPrime table from 1 MeV to 300 MeV in 46 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LowEPComptonModel : Emin=    0 eV  Emax=   20 MeV Fluo
        KleinNishina : Emin=   20 MeV Emax=  300 MeV Fluo

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

Rayl:  for gamma SubType=11 BuildTable=1
      Lambda table from 10 eV  to 150 keV, 19 bins/decade, spline: 0
      LambdaPrime table from 150 keV to 300 MeV in 62 bins 
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
   LivermoreRayleigh : Emin=    0 eV  Emax=  300 MeV  CullenGenerator

msc:  for e-  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  GoudsmitSaunderson : Emin=    1 MeV Emax=  300 MeV Nbins=40   1 MeV - 300 MeV
          StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm

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

eBrem:  for e-  XStype:4  SubType=3
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      LPM flag: 1 for E > 0.3 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    1 MeV Emax=  300 MeV  AngularGen2BS

e-_G4DNAElectronSolvation:  for e- SubType=58 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAOneStepThermalizationModel_Meesungnoen2002 : Emin=    0 eV  Emax=  7.4 eV 
          DummyModel : Emin=  7.4 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAElastic:  for e- SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAChampionElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAExcitation:  for e- SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNABornExcitationModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAIonisation:  for e- SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNABornIonisationModel : Emin=    0 eV  Emax=    1 MeV  deltaBorn Fluo
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAVibExcitation:  for e- SubType=54 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNASancheExcitationModel : Emin=    0 eV  Emax=  100 eV 
          DummyModel : Emin=  100 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

e-_G4DNAAttachment:  for e- SubType=55 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMeltonAttachmentModel : Emin=    0 eV  Emax=   13 eV 
          DummyModel : Emin=   13 eV  Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

msc:  for e+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  GoudsmitSaunderson : Emin=    0 eV  Emax=  300 MeV Nbins=120 100 eV  - 300 MeV
          StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm

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

eBrem:  for e+  XStype:4  SubType=3
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      LPM flag: 1 for E > 0.3 GeV,  VertexHighEnergyTh(GeV)= 100000
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
             eBremSB : Emin=    0 eV  Emax=  300 MeV  AngularGen2BS

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

proton_G4DNAElastic:  for proton SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

proton_G4DNAExcitation:  for proton SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  500 keV
DNABornExcitationModel : Emin=  500 keV Emax=  100 MeV
DNARPWBAExcitationModel : Emin=  100 MeV Emax=  300 MeV

proton_G4DNAIonisation:  for proton SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNARuddIonisationExtendedModel : Emin=    0 eV  Emax=  500 keV  deltaRudd Fluo
DNABornIonisationModel : Emin=  500 keV Emax=  100 MeV  deltaBorn Fluo
DNARPWBAIonisationModel : Emin=  100 MeV Emax=  300 MeV  deltaBorn Fluo

proton_G4DNAChargeDecrease:  for proton SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeDecreaseModel : Emin=    0 eV  Emax=  100 MeV
          DummyModel : Emin=  100 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

GenericIon_G4DNAIonisation:  for GenericIon SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha_G4DNAElastic:  for alpha SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

alpha_G4DNAExcitation:  for alpha SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNAMillerGreenExcitationModel : Emin=    1 keV Emax=  300 MeV

alpha_G4DNAIonisation:  for alpha SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha_G4DNAChargeDecrease:  for alpha SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNADingfelderChargeDecreaseModel : Emin=    1 keV Emax=  300 MeV

alpha+_G4DNAElastic:  for alpha+ SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

alpha+_G4DNAExcitation:  for alpha+ SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNAMillerGreenExcitationModel : Emin=    1 keV Emax=  300 MeV

alpha+_G4DNAIonisation:  for alpha+ SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

alpha+_G4DNAChargeIncrease:  for alpha+ SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNADingfelderChargeIncreaseModel : Emin=    1 keV Emax=  300 MeV

alpha+_G4DNAChargeDecrease:  for alpha+ SubType=56 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNADingfelderChargeDecreaseModel : Emin=    1 keV Emax=  300 MeV

msc:  for anti_proton  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

helium_G4DNAElastic:  for helium SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

helium_G4DNAExcitation:  for helium SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNAMillerGreenExcitationModel : Emin=    1 keV Emax=  300 MeV

helium_G4DNAIonisation:  for helium SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

helium_G4DNAChargeIncrease:  for helium SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=    1 keV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNADingfelderChargeIncreaseModel : Emin=    1 keV Emax=  300 MeV

hydrogen_G4DNAElastic:  for hydrogen SubType=51 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
  DNAIonElasticModel : Emin=    0 eV  Emax=    1 MeV
          DummyModel : Emin=    1 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

hydrogen_G4DNAExcitation:  for hydrogen SubType=52 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNAMillerGreenExcitationModel : Emin=    0 eV  Emax=  500 keV
          DummyModel : Emin=  500 keV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

hydrogen_G4DNAIonisation:  for hydrogen SubType=53 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
          DummyModel : Emin=    0 eV  Emax=  100 eV 
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
DNARuddIonisationExtendedModel : Emin=  100 eV  Emax=  300 MeV  deltaRudd Fluo

hydrogen_G4DNAChargeIncrease:  for hydrogen SubType=57 BuildTable=0
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
DNADingfelderChargeIncreaseModel : Emin=    0 eV  Emax=  100 MeV
          DummyModel : Emin=  100 MeV Emax=  300 MeV
          StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm

msc:  for kaon+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

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

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

msc:  for mu+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

hIoni:  for mu+  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
               Bragg : Emin=    0 eV  Emax=225.219 keV  deltaVI
          BetheBloch : Emin=225.219 keV Emax=  300 MeV  deltaVI

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

hIoni:  for mu-  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=225.219 keV  deltaVI
          BetheBloch : Emin=225.219 keV Emax=  300 MeV  deltaVI

msc:  for pi+  SubType= 10
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
        WentzelVIUni : Emin=    0 eV  Emax=  300 MeV
          StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm

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

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

hIoni:  for pi-  XStype:3  SubType=2
      dE/dx and range tables from 10 eV  to 300 MeV in 140 bins
      Lambda tables from threshold to 300 MeV, 20 bins/decade, spline: 1
      StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
      ===== EM models for the G4Region  DefaultRegionForTheWorld ======
            ICRU73QO : Emin=    0 eV  Emax=297.505 keV  deltaVI
          BetheBloch : Emin=297.505 keV Emax=  300 MeV  deltaVI
G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model.
See commands in /vis/modeling/trajectories/ for other options.
16

The following code …

========================================================================================
--> G4TaskRunManager::CreateAndStartWorkers() --> Creating 2 tasks with 1 events/task...
========================================================================================

G4WT0 > End of tracking primary particle, its final energy is :28.7595 MeV
G4WT1 > End of tracking primary particle, its final energy is :28.6105 MeV
G4WT0 > Physics stage ends
G4WT0 > *** G4Scheduler starts processing 
G4WT1 > Physics stage ends
G4WT1 > *** G4Scheduler starts processing 
Process 9415 stopped
* thread #3, stop reason = EXC_BAD_ACCESS (code=2, address=0x16ff1bf90)
    frame #0: 0x0000000102b7b468 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 4
libG4processes.dylib`G4KDTree::__NearestToNode<G4KDNode_Base>:
->  0x102b7b468 <+4>:  stp    d11, d10, [sp, #0x30]
    0x102b7b46c <+8>:  stp    d9, d8, [sp, #0x40]
    0x102b7b470 <+12>: stp    x28, x27, [sp, #0x50]
    0x102b7b474 <+16>: stp    x26, x25, [sp, #0x60]
Target 0: (neuron) stopped.
(lldb) bt
* thread #3, stop reason = EXC_BAD_ACCESS (code=2, address=0x16ff1bf90)
  * frame #0: 0x0000000102b7b468 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 4
    frame #1: 0x0000000102b7b850 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 1004
    frame #2: 0x0000000102b7b850 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 1004
    frame #3: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #4: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #5: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #6: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #7: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #8: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #9: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #10: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #11: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #12: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #13: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #14: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #15: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #16: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #17: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #18: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #19: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #20: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #21: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #22: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #23: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #24: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #25: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #26: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #27: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #28: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #29: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #30: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #31: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #32: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #33: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #34: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #35: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #36: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #37: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #38: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #39: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #40: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #41: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #42: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #43: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #44: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #45: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #46: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #47: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #48: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #49: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #50: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
17

The following code …

    frame #2894: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2895: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2896: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2897: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2898: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2899: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2900: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2901: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2902: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2903: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2904: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2905: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2906: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2907: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2908: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2909: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2910: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2911: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2912: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2913: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2914: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2915: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2916: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2917: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2918: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2919: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2920: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2921: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2922: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2923: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2924: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2925: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2926: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2927: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2928: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2929: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2930: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2931: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2932: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2933: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2934: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2935: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2936: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2937: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2938: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2939: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2940: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2941: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2942: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2943: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2944: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2945: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2946: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2947: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2948: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2949: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2950: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2951: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2952: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2953: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2954: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2955: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2956: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2957: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2958: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2959: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2960: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2961: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2962: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2963: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2964: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2965: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2966: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2967: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2968: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2969: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2970: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2971: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2972: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2973: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2974: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2975: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2976: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2977: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2978: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2979: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2980: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2981: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2982: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2983: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2984: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2985: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2986: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2987: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2988: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
18

The following code …

    frame #2989: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2990: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2991: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2992: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2993: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2994: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2995: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2996: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2997: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2998: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #2999: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3000: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3001: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3002: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3003: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3004: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3005: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3006: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3007: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3008: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3009: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3010: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3011: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3012: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3013: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3014: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3015: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3016: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3017: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3018: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3019: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3020: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3021: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3022: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3023: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3024: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3025: 0x0000000102b7b568 libG4processes.dylib`void G4KDTree::__NearestToNode<G4KDNode_Base>(G4KDNode_Base*, G4KDNode_Base*, G4KDNode_Base const&, std::__1::vector<G4KDNode_Base*, std::__1::allocator<G4KDNode_Base*>>&, double*, G4KDTree::HyperRect*, int&) + 260
    frame #3026: 0x0000000102b7b22c libG4processes.dylib`G4KDTree::Nearest(G4KDNode_Base*) + 236
    frame #3027: 0x0000000102a78878 libG4processes.dylib`G4ITFinder<G4Molecule>::FindNearest(G4Molecule const*, int) + 136
    frame #3028: 0x0000000102a80a08 libG4processes.dylib`G4DNAMoleculeEncounterStepper::CalculateStep(G4Track const&, double const&) + 948
    frame #3029: 0x0000000102a827b8 libG4processes.dylib`G4DNAMoleculeEncounterStepper::CalculateMinTimeStep(double, double) + 480
    frame #3030: 0x0000000102b3a8d0 libG4processes.dylib`G4ITModelProcessor::CalculateMinTimeStep(double, double) + 248
    frame #3031: 0x0000000102b7fc5c libG4processes.dylib`G4Scheduler::Stepping() + 748
    frame #3032: 0x0000000102b7f45c libG4processes.dylib`G4Scheduler::DoProcess() + 104
    frame #3033: 0x0000000102b7ee98 libG4processes.dylib`G4Scheduler::SynchronizeTracks() + 496
    frame #3034: 0x0000000102b7e858 libG4processes.dylib`G4Scheduler::Process() + 568
    frame #3035: 0x0000000102adc3ac libG4processes.dylib`G4DNAChemistryManager::Run() + 640
    frame #3036: 0x0000000100046658 neuron`StackingAction::NewStage() + 84
    frame #3037: 0x0000000100d94088 libG4event.dylib`G4StackManager::PopNextTrack(G4VTrajectory**) + 776
    frame #3038: 0x0000000100d604dc libG4event.dylib`G4EventManager::DoProcessing(G4Event*) + 1580
    frame #3039: 0x0000000101000bd8 libG4run.dylib`G4WorkerTaskRunManager::ProcessOneEvent(int) + 60
    frame #3040: 0x0000000101000b48 libG4run.dylib`G4WorkerTaskRunManager::DoEventLoop(int, char const*, int) + 232
    frame #3041: 0x00000001010024cc libG4run.dylib`G4WorkerTaskRunManager::DoWork() + 520
    frame #3042: 0x0000000100fec840 libG4run.dylib`std::__1::__packaged_task_func<std::__1::enable_if<std::is_void<void>::value, void>::type PTL::TaskGroup<void, void, 0l>::exec<G4TaskRunManager::AddEventTask(int)::$_4, void>(G4TaskRunManager::AddEventTask(int)::$_4)::'lambda'(), std::__1::allocator<std::__1::enable_if<std::is_void<void>::value, void>::type PTL::TaskGroup<void, void, 0l>::exec<G4TaskRunManager::AddEventTask(int)::$_4, void>(G4TaskRunManager::AddEventTask(int)::$_4)::'lambda'()>, void ()>::operator()() + 44
    frame #3043: 0x0000000100fec91c libG4run.dylib`std::__1::packaged_task<void ()>::operator()() + 80
    frame #3044: 0x00000001010bd484 libG4ptl.2.dylib`PTL::ThreadPool::execute_thread(PTL::VUserTaskQueue*) + 996
    frame #3045: 0x00000001010bcbc0 libG4ptl.2.dylib`PTL::ThreadPool::start_thread(PTL::ThreadPool*, std::__1::vector<std::__1::shared_ptr<PTL::ThreadData>, std::__1::allocator<std::__1::shared_ptr<PTL::ThreadData>>>*, long) + 688
    frame #3046: 0x00000001010c3b38 libG4ptl.2.dylib`void* std::__1::__thread_proxy[abi:ue170006]<std::__1::tuple<std::__1::unique_ptr<std::__1::__thread_struct, std::__1::default_delete<std::__1::__thread_struct>>, void (*)(PTL::ThreadPool*, std::__1::vector<std::__1::shared_ptr<PTL::ThreadData>, std::__1::allocator<std::__1::shared_ptr<PTL::ThreadData>>>*, long), PTL::ThreadPool*, std::__1::vector<std::__1::shared_ptr<PTL::ThreadData>, std::__1::allocator<std::__1::shared_ptr<PTL::ThreadData>>>*, unsigned long>>(void*) + 56
    frame #3047: 0x000000019fd2df94 libsystem_pthread.dylib`_pthread_start + 136
(lldb)

What should I do, please? @ivana and @guatelli

19

That ugly traceback looks to me like you have a literal stack overflow – the Geant4 function G4KDTree::__NearestToNode<G4KDNode_Base>() is calling itself recursively until your program runs out of program stack space.

I’ve run into this occasionally in my own code, where I was doing a recursive evaluation and didn’t put in a safety check (“abort job if too many calls”). But this is internal to Geant4-DNA, so you’ll need one of those experts to help.

20

Hello,

I am adding in the thread the developers who contributed to this example: @ngoc , @civanch can you, please, take a look ?

Thank you,