Environment variable "G4FORCE_RUN_MANAGER_TYPE" enabled with value == Serial. Forcing G4RunManager type... ############################################ !!! WARNING - FPE detection is activated !!! ############################################ ################################ !!! G4Backtrace is activated !!! ################################ ************************************************************** Geant4 version Name: geant4-11-04-patch-01 (13-March-2026) 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! 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) RayTracer (RT) VRML2FILE (VRML2FILE) gMocrenFile (gMocrenFile) TOOLSSG_OFFSCREEN (TSG_OFFSCREEN, TSG_FILE) OpenGLImmediateQt (OGLIQt, OGLI) OpenGLStoredQt (OGLSQt, OGLS) OpenGLImmediateXm (OGLIXm, OGLIQt_FALLBACK) OpenGLStoredXm (OGLSXm, OGLSQt_FALLBACK) OpenGLImmediateX (OGLIX, OGLIQt_FALLBACK, OGLIXm_FALLBACK) OpenGLStoredX (OGLSX, OGLSQt_FALLBACK, OGLSXm_FALLBACK) RayTracerX (RTX) RayTracerQt (RTQt) TOOLSSG_X11_GLES (TSG_X11_GLES, TSGX11, TSG_XT_GLES_FALLBACK) TOOLSSG_X11_ZB (TSG_X11_ZB, TSGX11ZB) TOOLSSG_XT_GLES (TSG_XT_GLES, TSGXt, TSG_QT_GLES_FALLBACK) TOOLSSG_XT_ZB (TSG_XT_ZB, TSGXtZB) TOOLSSG_QT_GLES (TSG_QT_GLES, TSGQt, TSG, OGL) TOOLSSG_QT_ZB (TSG_QT_ZB, TSGQtZB, TSGZB) You may choose a graphics system (driver) with a parameter of the command "/vis/open" or "/vis/sceneHandler/create", or you may omit the driver parameter and choose at run time: - by argument in the construction of G4VisExecutive - by environment variable "G4VIS_DEFAULT_DRIVER" - by entry in "~/.g4session" - by build flags. - Note: This feature is not allowed in batch mode. For further information see "examples/basic/B1/exampleB1.cc" and "vis.mac". 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. ---- 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! ---- ========================================================== ======================================================================= ====== 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 Positron annihilation at rest model AllisonPositronium Enable 3 gamma annihilation on fly 0 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 Use RiGe 5D e+e- pair production model by muons 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 100 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 Urban Use built-in Birks saturation 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 1 MeV 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 1 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 ======== ======================================================================= Max DNA electron energy 1 MeV Max DNA ion energy 300 MeV 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 4 regions: DefaultRegionForTheWorld 1 1 0 Soma 1 1 0 Dendrites 1 1 0 Axon 1 1 0 ### === Auger flag: 1 ### === Ignore cuts flag: 1 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 e- SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV 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 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 20 bins/decade, spline: 1 StepFunction=(0.2, 0.01 mm), integ: 3, fluct: 1, linLossLim= 0.01 ===== EM models for the G4Region DefaultRegionForTheWorld ====== LowEnergyIoni : Emin= 0 eV Emax= 100 keV deltaVI MollerBhabha : Emin= 100 keV Emax= 100 TeV deltaVI eBrem: for e- XStype:4 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 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e- 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: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for e- XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV 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= 100 MeV Emax= 100 TeV msc: for e+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== GoudsmitSaunderson : Emin= 0 eV Emax= 100 MeV Nbins=120 100 eV - 100 MeV StepLim=SafetyPlus Rfact=0.08 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm WentzelVIUni : Emin= 100 MeV Emax= 100 TeV Nbins=120 100 MeV - 100 TeV 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 100 eV to 100 TeV in 240 bins Lambda tables from threshold to 100 TeV, 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= 100 TeV deltaVI eBrem: for e+ XStype:4 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 LPM flag: 1 for E > 1 GeV, VertexHighEnergyTh(GeV)= 100000 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eBremSB : Emin= 0 eV Emax= 1 GeV AngularGen2BS eBremLPM : Emin= 1 GeV Emax= 100 TeV AngularGen2BS ePairProd: for e+ 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: 0 Sampling table 25x1001 from 0.1 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== ePairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eplus2gg : Emin= 0 eV Emax= 100 TeV CoulombScat: for e+ XStype:1 SubType=1 BuildTable=1 Lambda table from 100 MeV to 100 TeV, 20 bins/decade, spline: 0 ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 100 MeV Emax= 100 TeV msc: for proton SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=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 nuclearStopping: for GenericIon SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for alpha SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== UrbanMsc : Emin= 0 eV Emax= 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm ionIoni: for alpha 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.02 mm), integ: 3, fluct: 1, linLossLim= 0.02 ===== EM models for the G4Region DefaultRegionForTheWorld ====== BraggIon : Emin= 0 eV Emax=7.9452 MeV deltaVI BetheBloch : Emin=7.9452 MeV Emax= 100 TeV deltaVI nuclearStopping: for alpha SubType=8 BuildTable=0 ===== EM models for the G4Region DefaultRegionForTheWorld ====== ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV msc: for anti_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 anti_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 ====== ICRU73QO : Emin= 0 eV Emax= 2 MeV deltaVI BetheBloch : Emin= 2 MeV Emax= 100 TeV deltaVI hBrems: for anti_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 anti_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 anti_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 msc: for kaon+ 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 kaon+ 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=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon+ 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 kaon+ 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 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for kaon+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 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 msc: for kaon- 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 kaon- 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 ====== ICRU73QO : Emin= 0 eV Emax=1.05231 MeV deltaVI BetheBloch : Emin=1.05231 MeV Emax= 100 TeV deltaVI hBrems: for kaon- 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 kaon- 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 18x1001, from 3.94942 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== hPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for kaon- XStype:1 SubType=1 BuildTable=1 Used Lambda table of kaon+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for mu+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=240 100 eV - 100 TeV StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=3 Llim=1 mm muIoni: for mu+ 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= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu+ 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 ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied muPairProd: for mu+ 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 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for mu+ XStype:1 SubType=1 BuildTable=1 Lambda table from threshold to 100 TeV, 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 msc: for mu- 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 muIoni: for mu- 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 ====== ICRU73QO : Emin= 0 eV Emax= 200 keV deltaVI MuBetheBloch : Emin= 200 keV Emax= 100 TeV deltaVI muBrems: for mu- 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 ====== MuBrem : Emin= 0 eV Emax= 100 TeV ModifiedMephi ===== Limit on energy threshold has been applied muPairProd: for mu- 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 21x1001, from 0.85 GeV to 100 TeV ===== EM models for the G4Region DefaultRegionForTheWorld ====== muPairProd : Emin= 0 eV Emax= 100 TeV ModifiedMephi CoulombScat: for mu- XStype:1 SubType=1 BuildTable=1 Used Lambda table of mu+ ThetaMin(p) < Theta(degree) < 180, pLimit(GeV^1)= 0.139531 ===== EM models for the G4Region DefaultRegionForTheWorld ====== eCoulombScattering : Emin= 0 eV Emax= 100 TeV msc: for pi+ SubType= 10 ===== EM models for the G4Region DefaultRegionForTheWorld ====== WentzelVIUni : Emin= 0 eV Emax= 100 TeV Nbins=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 pi+ 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=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi+ 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 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 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 msc: for pi- 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 pi- 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 ====== ICRU73QO : Emin= 0 eV Emax=297.505 keV deltaVI BetheBloch : Emin=297.505 keV Emax= 100 TeV deltaVI hBrems: for pi- 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 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 End of tracking primary particle, its final energy is :3.47476 GeV Physics stage ends End of tracking primary particle, its final energy is :3.47513 GeV Physics stage ends ======= The summary of simulation results 'neuron' ======== Primary particle = C12 Kinetic energy of beam = 3480 A*MeV Particle traversals the neuron = 1 of 2 Full LET of beam as formulas = 12.827 keV/um Mean LET of beam as simulation = 11.42 +- 0.45547 keV/um Mean track length of beam = 406.78 +- 0.010536 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= 4366 ionIoni= 2662 msc= 723 List of generated particles outside neuron structure: e- : 3245 Emean = 1.6075 keV ( 242.31 eV --> 424.07 keV) Eflow/event = 2.6082 MeV 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: 9.2907 Bounding slice: 3.5625 Neuron: 0.00036408 Soma: 0 Dendrites: 0.00036408 Axon: 0 Number of compartments traversed by particle tracks : Soma: 0 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 = 1.07753 s Graphics systems deleted. Visualization Manager deleting...