Warning, /geant4/examples/extended/medical/dna/neuron/neuron.out is written in an unsupported language. File is not indexed.
0001 Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type...
0002
0003 ############################################
0004 !!! WARNING - FPE detection is activated !!!
0005 ############################################
0006
0007
0008 ################################
0009 !!! G4Backtrace is activated !!!
0010 ################################
0011
0012
0013 **************************************************************
0014 Geant4 version Name: geant4-11-03-ref-06 (30-June-2025)
0015 Copyright : Geant4 Collaboration
0016 References : NIM A 506 (2003), 250-303
0017 : IEEE-TNS 53 (2006), 270-278
0018 : NIM A 835 (2016), 186-225
0019 WWW : http://geant4.org/
0020 **************************************************************
0021
0022 NeuronLoadDataFile::SingleNeuronSWCfile: number of strings=2922 in file GranuleCell-Nr2.CNG.swc
0023 Total number of compartments into Neuron : 2914
0024
0025 Livermore + DNAphysics is activated!
0026 ---- Begin of Neuron Construction! ----
0027 ==========================================================
0028 Side length of word volume is calculated : 407.347 um
0029 Radius of homogeneous medium is calculated : 203.674 um
0030 Overall dimensions (um) of neuron morphology :
0031 width = 250.225 height = 315.035 depth = 63.8148
0032 Volume (um3), surface (um2) and mass (ug) of Bounding Slice are calculated :
0033 5.03049e+06 229803 5.03049
0034
0035 Volume (um3), surface (um2) and mass(ug) of Neuron are calculated :
0036 12100 30379.7 0.0121
0037 Total number of compartments into Neuron : 2914
0038 Shift values (um) for Neuron translation are calculated :
0039 shiftX = 35.13 shiftY = 99.995 shiftZ = 2.025
0040 ---- Soma for construction: ----
0041 Total number of compartments into Soma : 3
0042 ---- Dendrites for construction: ----
0043 Total number of compartments into Dendrites: 2775
0044 ---- Axon for construction: ----
0045 Total number of compartments into Axon : 136
0046 ---- Spines not found! ----
0047
0048 ---- End of Neuron Construction! ----
0049 ==========================================================
0050
0051 Visualization Manager instantiating with verbosity "warnings (3)"...
0052 Visualization Manager initialising...
0053 Registering graphics systems...
0054
0055 You have successfully registered the following graphics systems.
0056 Registered graphics systems are:
0057 ASCIITree (ATree)
0058 DAWNFILE (DAWNFILE)
0059 G4HepRepFile (HepRepFile)
0060 RayTracer (RT)
0061 VRML2FILE (VRML2FILE)
0062 gMocrenFile (gMocrenFile)
0063 TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE)
0064 OpenGLImmediateQt (OGLIQt, OGLI)
0065 OpenGLStoredQt (OGLSQt, OGLS)
0066 OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK)
0067 OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK)
0068 OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK)
0069 OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK)
0070 RayTracerX (RTX)
0071 RayTracerQt (RTQt)
0072 Qt3D (Qt3D)
0073 TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK)
0074 TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB)
0075 TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK)
0076 TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB)
0077 TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG, OGL)
0078 TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB)
0079 You may choose a graphics system (driver) with a parameter of
0080 the command "/vis/open" or "/vis/sceneHandler/create",
0081 or you may omit the driver parameter and choose at run time:
0082 - by argument in the construction of G4VisExecutive
0083 - by environment variable "G4VIS_DEFAULT_DRIVER"
0084 - by entry in "~/.g4session"
0085 - by build flags.
0086 - Note: This feature is not allowed in batch mode.
0087 For further information see "examples/basic/B1/exampleB1.cc"
0088 and "vis.mac".
0089
0090 Registering model factories...
0091
0092 You have successfully registered the following model factories.
0093 Registered model factories:
0094 generic
0095 drawByAttribute
0096 drawByCharge
0097 drawByOriginVolume
0098 drawByParticleID
0099 drawByEncounteredVolume
0100
0101 Registered models:
0102 None
0103
0104 Registered filter factories:
0105 attributeFilter
0106 chargeFilter
0107 originVolumeFilter
0108 particleFilter
0109 encounteredVolumeFilter
0110
0111 Registered filters:
0112 None
0113
0114 You have successfully registered the following user vis actions.
0115 Run Duration User Vis Actions: none
0116 End of Event User Vis Actions: none
0117 End of Run User Vis Actions: none
0118
0119 Some /vis commands (optionally) take a string to specify colour.
0120 "/vis/list" to see available colours.
0121 =======================================================================
0122 ====== Electromagnetic Physics Parameters ========
0123 =======================================================================
0124 LPM effect enabled 1
0125 Enable creation and use of sampling tables 0
0126 Apply cuts on all EM processes 0
0127 Use combined TransportationWithMsc Disabled
0128 Use general process 0
0129 Enable linear polarisation for gamma 0
0130 Enable photoeffect sampling below K-shell 1
0131 Enable sampling of quantum entanglement 0
0132 X-section factor for integral approach 0.8
0133 Min kinetic energy for tables 100 eV
0134 Max kinetic energy for tables 100 TeV
0135 Number of bins per decade of a table 20
0136 Verbose level 1
0137 Verbose level for worker thread 0
0138 Bremsstrahlung energy threshold above which
0139 primary e+- is added to the list of secondary 100 TeV
0140 Bremsstrahlung energy threshold above which primary
0141 muon/hadron is added to the list of secondary 100 TeV
0142 Positron annihilation at rest model AllisonPositronium
0143 Enable 3 gamma annihilation on fly 0
0144 Lowest triplet kinetic energy 1 MeV
0145 Enable sampling of gamma linear polarisation 0
0146 5D gamma conversion model type 0
0147 5D gamma conversion model on isolated ion 0
0148 Use RiGe 5D e+e- pair production model by muons 0
0149 Livermore data directory epics_2017
0150 =======================================================================
0151 ====== Ionisation Parameters ========
0152 =======================================================================
0153 Step function for e+- (0.2, 0.01 mm)
0154 Step function for muons/hadrons (0.1, 0.05 mm)
0155 Step function for light ions (0.1, 0.02 mm)
0156 Step function for general ions (0.1, 0.001 mm)
0157 Lowest e+e- kinetic energy 100 eV
0158 Lowest muon/hadron kinetic energy 1 keV
0159 Use ICRU90 data 1
0160 Fluctuations of dE/dx are enabled 1
0161 Type of fluctuation model for leptons and hadrons Urban
0162 Use built-in Birks satuaration 0
0163 Build CSDA range enabled 0
0164 Use cut as a final range enabled 0
0165 Enable angular generator interface 1
0166 Max kinetic energy for CSDA tables 1 GeV
0167 Max kinetic energy for NIEL computation 1 MeV
0168 Linear loss limit 0.01
0169 Read data from file for e+e- pair production by mu 0
0170 =======================================================================
0171 ====== Multiple Scattering Parameters ========
0172 =======================================================================
0173 Type of msc step limit algorithm for e+- 2
0174 Type of msc step limit algorithm for muons/hadrons 0
0175 Msc lateral displacement for e+- enabled 1
0176 Msc lateral displacement for muons and hadrons 1
0177 Urban msc model lateral displacement alg96 1
0178 Range factor for msc step limit for e+- 0.08
0179 Range factor for msc step limit for muons/hadrons 0.2
0180 Geometry factor for msc step limitation of e+- 2.5
0181 Safety factor for msc step limit for e+- 0.6
0182 Skin parameter for msc step limitation of e+- 3
0183 Lambda limit for msc step limit for e+- 1 mm
0184 Use Mott correction for e- scattering 1
0185 Factor used for dynamic computation of angular
0186 limit between single and multiple scattering 1
0187 Fixed angular limit between single
0188 and multiple scattering 3.1416 rad
0189 Upper energy limit for e+- multiple scattering 100 MeV
0190 Type of electron single scattering model 0
0191 Type of nuclear form-factor 1
0192 Screening factor 1
0193 =======================================================================
0194 ====== Atomic Deexcitation Parameters ========
0195 =======================================================================
0196 Fluorescence enabled 1
0197 Directory in G4LEDATA for fluorescence data files fluor
0198 Auger electron cascade enabled 1
0199 PIXE atomic de-excitation enabled 0
0200 De-excitation module ignores cuts 1
0201 Type of PIXE cross section for hadrons Empirical
0202 Type of PIXE cross section for e+- Livermore
0203 =======================================================================
0204 ====== DNA Physics Parameters ========
0205 =======================================================================
0206 Use fast sampling in DNA models 0
0207 Use Stationary option in DNA models 0
0208 Use DNA with multiple scattering of e- 0
0209 Use DNA e- solvation model type 11003
0210 =======================================================================
0211
0212 ### === Deexcitation model UAtomDeexcitation is activated for 4 regions:
0213 DefaultRegionForTheWorld 1 1 0
0214 Soma 1 1 0
0215 Dendrites 1 1 0
0216 Axon 1 1 0
0217 ### === Auger flag: 1
0218 ### === Ignore cuts flag: 1
0219
0220 phot: for gamma SubType=12 BuildTable=0
0221 LambdaPrime table from 200 keV to 100 TeV in 174 bins
0222 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0223 LivermorePhElectric : Emin= 0 eV Emax= 100 TeV SauterGavrila Fluo
0224
0225 compt: for gamma SubType=13 BuildTable=1
0226 Lambda table from 100 eV to 1 MeV, 20 bins/decade, spline: 1
0227 LambdaPrime table from 1 MeV to 100 TeV in 160 bins
0228 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0229 LivermoreCompton : Emin= 0 eV Emax= 1 GeV Fluo
0230 KleinNishina : Emin= 1 GeV Emax= 100 TeV Fluo
0231
0232 conv: for gamma SubType=14 BuildTable=1
0233 Lambda table from 1.022 MeV to 100 TeV, 20 bins/decade, spline: 1
0234 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0235 Livermore5DConversion : Emin= 0 eV Emax= 100 TeV ModifiedTsai
0236
0237 Rayl: for gamma SubType=11 BuildTable=1
0238 Lambda table from 100 eV to 150 keV, 20 bins/decade, spline: 0
0239 LambdaPrime table from 150 keV to 100 TeV in 176 bins
0240 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0241 LivermoreRayleigh : Emin= 0 eV Emax= 100 TeV CullenGenerator
0242
0243 msc: for e- SubType= 10
0244 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0245 GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV
0246 StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0247 WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV
0248 StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0249
0250 eIoni: for e- XStype:3 SubType=2
0251 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0252 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0253 StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
0254 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0255 LowEnergyIoni : Emin= 0 eV Emax= 100 keV deltaVI
0256 MollerBhabha : Emin= 100 keV Emax= 100 TeV deltaVI
0257
0258 eBrem: for e- XStype:4 SubType=3
0259 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0260 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0261 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0262 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0263 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0264 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0265
0266 ePairProd: for e- XStype:1 SubType=4
0267 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0268 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0269 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0270 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0271 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0272
0273 CoulombScat: for e- XStype:1 SubType=1 BuildTable=1
0274 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
0275 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0276 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0277 eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
0278
0279 msc: for e+ SubType= 10
0280 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0281 GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV
0282 StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0283 WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV
0284 StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0285
0286 eIoni: for e+ XStype:3 SubType=2
0287 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0288 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0289 StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01
0290 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0291 MollerBhabha : Emin= 0 eV Emax= 100 TeV deltaVI
0292
0293 eBrem: for e+ XStype:4 SubType=3
0294 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0295 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0296 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000
0297 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0298 eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS
0299 eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS
0300
0301 ePairProd: for e+ XStype:1 SubType=4
0302 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0303 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 0
0304 Sampling table 25x1001 from 0.1 GeV to 100 TeV
0305 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0306 ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0307
0308 annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
0309 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0310 eplus2gg : Emin= 0 eV Emax= 100 TeV
0311
0312 CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1
0313 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0
0314 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0315 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0316 eCoulombScattering : Emin= 100 MeV Emax= 100 TeV
0317
0318 msc: for proton SubType= 10
0319 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0320 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0321 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0322
0323 hIoni: for proton XStype:3 SubType=2
0324 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0325 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0326 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0327 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0328 Bragg : Emin= 0 eV Emax= 2 MeV deltaVI
0329 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0330
0331 hBrems: for proton XStype:1 SubType=3
0332 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0333 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0334 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0335 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0336 ===== Limit on energy threshold has been applied
0337
0338 hPairProd: for proton XStype:1 SubType=4
0339 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0340 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0341 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0342 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0343 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0344
0345 CoulombScat: for proton XStype:1 SubType=1 BuildTable=1
0346 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0
0347 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0348 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0349 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0350
0351 nuclearStopping: for proton SubType=8 BuildTable=0
0352 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0353 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0354
0355 msc: for GenericIon SubType= 10
0356 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0357 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0358 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0359
0360 ionIoni: for GenericIon XStype:3 SubType=2
0361 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0362 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0363 StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
0364 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0365 LindhardSorensen : Emin= 0 eV Emax= 100 TeV deltaVI
0366
0367 nuclearStopping: for GenericIon SubType=8 BuildTable=0
0368 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0369 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0370
0371 msc: for alpha SubType= 10
0372 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0373 UrbanMsc : Emin= 0 eV Emax= 100 TeV
0374 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0375
0376 ionIoni: for alpha XStype:3 SubType=2
0377 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0378 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0379 StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
0380 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0381 BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI
0382 BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI
0383
0384 nuclearStopping: for alpha SubType=8 BuildTable=0
0385 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0386 ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
0387
0388 msc: for anti_proton SubType= 10
0389 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0390 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0391 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0392
0393 hIoni: for anti_proton XStype:3 SubType=2
0394 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0395 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0396 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0397 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0398 ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI
0399 BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI
0400
0401 hBrems: for anti_proton XStype:1 SubType=3
0402 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0403 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0404 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0405 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0406 ===== Limit on energy threshold has been applied
0407
0408 hPairProd: for anti_proton XStype:1 SubType=4
0409 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0410 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0411 Sampling table 17x1001 from 7.50618 GeV to 100 TeV
0412 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0413 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0414
0415 CoulombScat: for anti_proton XStype:1 SubType=1 BuildTable=1
0416 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0
0417 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0418 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0419 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0420
0421 msc: for kaon+ SubType= 10
0422 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0423 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0424 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0425
0426 hIoni: for kaon+ XStype:3 SubType=2
0427 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0428 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0429 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0430 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0431 Bragg : Emin= 0 eV Emax=1.05231 MeV deltaVI
0432 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0433
0434 hBrems: for kaon+ XStype:1 SubType=3
0435 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0436 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0437 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0438 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0439 ===== Limit on energy threshold has been applied
0440
0441 hPairProd: for kaon+ XStype:1 SubType=4
0442 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0443 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0444 Sampling table 18x1001 from 3.94942 GeV to 100 TeV
0445 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0446 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0447
0448 CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1
0449 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0
0450 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0451 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0452 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0453
0454 msc: for kaon- SubType= 10
0455 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0456 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0457 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0458
0459 hIoni: for kaon- XStype:3 SubType=2
0460 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0461 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0462 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0463 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0464 ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI
0465 BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI
0466
0467 hBrems: for kaon- XStype:1 SubType=3
0468 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0469 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0470 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0471 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0472 ===== Limit on energy threshold has been applied
0473
0474 hPairProd: for kaon- XStype:1 SubType=4
0475 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0476 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0477 Sampling table 18x1001 from 3.94942 GeV to 100 TeV
0478 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0479 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0480
0481 CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1
0482 Used Lambda table of kaon+
0483 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0484 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0485 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0486
0487 msc: for mu+ SubType= 10
0488 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0489 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0490 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0491
0492 muIoni: for mu+ XStype:3 SubType=2
0493 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0494 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0495 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0496 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0497 Bragg : Emin= 0 eV Emax= 200 keV deltaVI
0498 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0499
0500 muBrems: for mu+ XStype:1 SubType=3
0501 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0502 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0503 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0504 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0505 ===== Limit on energy threshold has been applied
0506
0507 muPairProd: for mu+ XStype:1 SubType=4
0508 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0509 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0510 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0511 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0512 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0513
0514 CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1
0515 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0
0516 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0517 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0518 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0519
0520 msc: for mu- SubType= 10
0521 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0522 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0523 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0524
0525 muIoni: for mu- XStype:3 SubType=2
0526 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0527 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0528 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0529 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0530 ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI
0531 MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI
0532
0533 muBrems: for mu- XStype:1 SubType=3
0534 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0535 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0536 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0537 MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0538 ===== Limit on energy threshold has been applied
0539
0540 muPairProd: for mu- XStype:1 SubType=4
0541 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0542 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0543 Sampling table 21x1001 from 0.85 GeV to 100 TeV
0544 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0545 muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0546
0547 CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1
0548 Used Lambda table of mu+
0549 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0550 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0551 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0552
0553 msc: for pi+ SubType= 10
0554 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0555 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0556 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0557
0558 hIoni: for pi+ XStype:3 SubType=2
0559 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0560 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0561 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0562 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0563 Bragg : Emin= 0 eV Emax=297.505 keV deltaVI
0564 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0565
0566 hBrems: for pi+ XStype:1 SubType=3
0567 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0568 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0569 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0570 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0571 ===== Limit on energy threshold has been applied
0572
0573 hPairProd: for pi+ XStype:1 SubType=4
0574 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0575 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0576 Sampling table 20x1001 from 1.11656 GeV to 100 TeV
0577 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0578 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0579
0580 CoulombScat: for pi+ XStype:1 SubType=1 BuildTable=1
0581 Lambda table from threshold to 100 TeV, 20 bins/decade, spline: 0
0582 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0583 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0584 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0585
0586 msc: for pi- SubType= 10
0587 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0588 WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV
0589 StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm
0590
0591 hIoni: for pi- XStype:3 SubType=2
0592 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0593 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0594 StepFunction=(0.1, 0.05 mm), integ: 3, fluct: 1, linLossLim= 0.01
0595 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0596 ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI
0597 BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI
0598
0599 hBrems: for pi- XStype:1 SubType=3
0600 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0601 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0602 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0603 hBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0604 ===== Limit on energy threshold has been applied
0605
0606 hPairProd: for pi- XStype:1 SubType=4
0607 dE/dx and range tables from 100 eV to 100 TeV in 240 bins
0608 Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1
0609 Sampling table 20x1001 from 1.11656 GeV to 100 TeV
0610 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0611 hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi
0612
0613 CoulombScat: for pi- XStype:1 SubType=1 BuildTable=1
0614 Used Lambda table of pi+
0615 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531
0616 ===== EM models for the G4Region DefaultRegionForTheWorld ======
0617 eCoulombScattering : Emin= 0 eV Emax= 100 TeV
0618 End of tracking primary particle, its final energy is :3.47476 GeV
0619 Physics stage ends
0620 End of tracking primary particle, its final energy is :3.47537 GeV
0621 Physics stage ends
0622
0623 ======= The summary of simulation results 'neuron' ========
0624
0625 Primary particle = C12
0626 Kinetic energy of beam = 3480 A*MeV
0627 Particle traversals the neuron = 1 of 2
0628 Full LET of beam as formulas = 12.827 keV/um
0629 Mean LET of beam as simulation = 10.892 +- 0.31423 keV/um
0630 Mean track length of beam = 399.22 +- 0.0029698 um
0631 Particle fluence = 383.66 particles/cm^2
0632 Fluence dose (full) = 7.8848e-06 Gy
0633 Fluence dose ber beam = 3.9424e-06 Gy
0634
0635 List of generated physical process:
0636 eIoni= 3986
0637 ionIoni= 2469
0638 msc= 5330
0639
0640
0641 List of generated particles outside neuron structure:
0642 e- : 2997 Emean = 1.7843 keV ( 242.23 eV --> 434.14 keV) Eflow/event = 2.6738 MeV
0643
0644 Number of secondary particles inside neuron structure:
0645 e- : 14
0646
0647 Number of molecular products inside neuron structure:
0648 time: 1 ps - 10 ps
0649
0650 Total energy (MeV) deposition :
0651 All volume: 8.6967
0652 Bounding slice: 2.5219
0653 Neuron: 0.014207
0654 Soma: 0
0655 Dendrites: 0.014207
0656 Axon: 0
0657
0658 Number of compartments traversed by particle tracks :
0659 Soma: 0 of total: 3
0660 Dendrites: 3 of total: 2775
0661 Axon: 0 of total: 136
0662
0663
0664 Dendritic (or Axon) compartmental energy deposits
0665 at the distance from Soma have been written into *.out files:
0666 Dend3DEdep.out, Axon3DEdep.out, Soma3DEdep.out
0667 Outputs of energy deposition per event written in data file:
0668 OutputPerEvent.out
0669 Calculation time = 1.12602 s
0670
0671 Graphics systems deleted.
0672 Visualization Manager deleting...