Gabriele Cosmo
673 lines
32 KiB
Plaintext
673 lines
32 KiB
Plaintext
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...
|