687 lines
32 KiB
Plaintext
687 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/
|
|
**************************************************************
|
|
|
|
/run/verbose 1
|
|
/tracking/verbose 0
|
|
/random/setSeeds 102 001
|
|
/mygeom/maxRange 49. um
|
|
/mygeom/material G4_WATER
|
|
/mygeom/hitSelRegZ 4. um
|
|
/mygeom/hitSelRegXY 100. um
|
|
/mygeom/siteRadius 0.5 um
|
|
/physics/addPhysics dna_opt2
|
|
PhysicsList::AddPhysicsList: <dna_opt2>
|
|
G4EmDNAPhysics_option2 physics package has been activated.
|
|
/run/setCutForAGivenParticle gamma 0.01 um
|
|
/run/setCutForAGivenParticle e- 0.01 um
|
|
/run/setCutForAGivenParticle e+ 0.01 um
|
|
/run/setCutForAGivenParticle proton 0.01 um
|
|
/run/numberOfThreads 8
|
|
*** /run/numberOfThreads command is issued in sequential mode.
|
|
Command is ignored.
|
|
/run/initialize
|
|
Material: G4_WATER
|
|
|
|
***** Table : Nb of materials = 1 *****
|
|
|
|
Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.375 cm
|
|
Imean: 78.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: H (H) Z = 1.0 N = 1 A = 1.008 g/mole
|
|
---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
|
|
---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
|
|
ElmMassFraction: 11.19 % ElmAbundance 66.67 %
|
|
|
|
---> Element: O (O) Z = 8.0 N = 16 A = 15.999 g/mole
|
|
---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
|
|
---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
|
|
---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
|
|
ElmMassFraction: 88.81 % ElmAbundance 33.33 %
|
|
|
|
|
|
|
|
|
|
================ Volume Parameters ================
|
|
Physical Volume Name: World
|
|
Position: (0,0,0) mm
|
|
Material: G4_WATER
|
|
Shape: Box
|
|
Dimensions (full): 0.102 x 0.102 x 0.104 mm
|
|
===================================================
|
|
|
|
|
|
================ Volume Parameters ================
|
|
Physical Volume Name: SDbox
|
|
Position: (0,0,0) mm
|
|
Material: G4_WATER
|
|
Shape: Box
|
|
Dimensions (full): 0.102 x 0.102 x 0.006 mm
|
|
===================================================
|
|
|
|
/gun/particle proton
|
|
/gun/energy 20.13 MeV
|
|
/beam/position/Z0 -50.99 um
|
|
/run/particle/applyCuts true e-
|
|
/run/particle/applyCuts true gamma
|
|
/run/printProgress 1
|
|
/run/beamOn 50
|
|
=======================================================================
|
|
====== 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 600 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
|
|
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 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 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 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 ========
|
|
=======================================================================
|
|
Max DNA electron energy 1 MeV
|
|
Max DNA ion energy 300 MeV
|
|
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 1 region:
|
|
DefaultRegionForTheWorld 1 1 0
|
|
### === Auger flag: 1
|
|
### === Ignore cuts flag: 1
|
|
|
|
phot: for gamma SubType=12 BuildTable=0
|
|
LambdaPrime table from 200 keV to 600 MeV in 72 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
LivermorePhElectric : Emin= 0 eV Emax= 600 MeV SauterGavrila Fluo
|
|
|
|
compt: for gamma SubType=13 BuildTable=1
|
|
Lambda table from 10 eV to 1 MeV, 21 bins/decade, spline: 1
|
|
LambdaPrime table from 1 MeV to 600 MeV in 57 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
LowEPComptonModel : Emin= 0 eV Emax= 20 MeV Fluo
|
|
KleinNishina : Emin= 20 MeV Emax= 600 MeV Fluo
|
|
|
|
conv: for gamma SubType=14 BuildTable=1
|
|
Lambda table from 1.022 MeV to 600 MeV, 28 bins/decade, spline: 1
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
BetheHeitlerLPM : Emin= 0 eV Emax= 600 MeV ModifiedTsai
|
|
|
|
Rayl: for gamma SubType=11 BuildTable=1
|
|
Lambda table from 10 eV to 150 keV, 21 bins/decade, spline: 0
|
|
LambdaPrime table from 150 keV to 600 MeV in 74 bins
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
LivermoreRayleigh : Emin= 0 eV Emax= 600 MeV CullenGenerator
|
|
|
|
msc: for e- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
GoudsmitSaunderson : Emin= 1 MeV Emax= 600 MeV Nbins=60 1 MeV - 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
eBrem: for e- XStype:4 SubType=3
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 MeV, 20 bins/decade, spline: 1
|
|
LPM flag: 1 for E > 0.6 GeV, VertexHighEnergyTh(GeV)= 100000
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eBremSB : Emin= 1 MeV Emax= 600 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= 600 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= 600 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= 600 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= 600 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= 600 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= 600 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= 600 MeV Nbins=140 100 eV - 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
eBrem: for e+ XStype:4 SubType=3
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 MeV, 20 bins/decade, spline: 1
|
|
LPM flag: 1 for E > 0.6 GeV, VertexHighEnergyTh(GeV)= 100000
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eBremSB : Emin= 0 eV Emax= 600 MeV AngularGen2BS
|
|
|
|
annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
eplus2gg : Emin= 0 eV Emax= 600 MeV
|
|
|
|
msc: for proton SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 1 MeV Emax= 600 MeV Nbins=60 1 MeV - 600 MeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm
|
|
|
|
hIoni: for proton XStype:3 SubType=2
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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 ======
|
|
BetheBloch : Emin= 300 MeV Emax= 600 MeV deltaVI
|
|
|
|
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= 600 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
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
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
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
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= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
msc: for GenericIon SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 1 MeV Emax= 600 MeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm
|
|
|
|
ionIoni: for GenericIon XStype:3 SubType=2
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 MeV, 20 bins/decade, spline: 1
|
|
StepFunction=(0.1, 0.001 mm), integ: 3, fluct: 1, linLossLim= 0.02
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
BetheBloch : Emin= 300 MeV Emax= 600 MeV deltaVI
|
|
|
|
GenericIon_G4DNAIonisation: for GenericIon SubType=53 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNAIonIonisationModel : Emin= 0 eV Emax= 300 MeV
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
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= 1 MeV Emax= 600 MeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm
|
|
|
|
ionIoni: for alpha XStype:3 SubType=2
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 MeV, 20 bins/decade, spline: 1
|
|
StepFunction=(0.1, 0.02 mm), integ: 3, fluct: 1, linLossLim= 0.02
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
BetheBloch : Emin= 300 MeV Emax= 600 MeV deltaVI
|
|
|
|
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= 600 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
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
alpha_G4DNAIonisation: for alpha SubType=53 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNARuddIonisationExtendedModel : Emin= 0 eV Emax= 300 MeV deltaRudd Fluo
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
alpha_G4DNAChargeDecrease: for alpha SubType=56 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNADingfelderChargeDecreaseModel : Emin= 0 eV Emax= 300 MeV
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
msc: for alpha+ SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
UrbanMsc : Emin= 1 MeV Emax= 600 MeV
|
|
StepLim=Minimal Rfact=0.2 Gfact=2.5 Sfact=0.6 DispFlag:0 Skin=3 Llim=1 mm
|
|
|
|
hIoni: for alpha+ XStype:3 SubType=2
|
|
dE/dx and range tables from 10 eV to 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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 ======
|
|
BetheBloch : Emin= 300 MeV Emax= 600 MeV deltaVI
|
|
|
|
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= 600 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
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
alpha+_G4DNAIonisation: for alpha+ SubType=53 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNARuddIonisationExtendedModel : Emin= 0 eV Emax= 300 MeV deltaRudd Fluo
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
alpha+_G4DNAChargeIncrease: for alpha+ SubType=57 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNADingfelderChargeIncreaseModel : Emin= 0 eV Emax= 300 MeV
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
alpha+_G4DNAChargeDecrease: for alpha+ SubType=56 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNADingfelderChargeDecreaseModel : Emin= 0 eV Emax= 300 MeV
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
msc: for anti_proton SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 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= 600 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
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
helium_G4DNAIonisation: for helium SubType=53 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNARuddIonisationExtendedModel : Emin= 0 eV Emax= 300 MeV deltaRudd Fluo
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
helium_G4DNAChargeIncrease: for helium SubType=57 BuildTable=0
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
DNADingfelderChargeIncreaseModel : Emin= 0 eV Emax= 300 MeV
|
|
DummyModel : Emin= 300 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
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= 600 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= 600 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 ======
|
|
DNARuddIonisationExtendedModel : Emin= 0 eV Emax= 100 MeV deltaRudd Fluo
|
|
DummyModel : Emin= 100 MeV Emax= 600 MeV
|
|
StepLim=UseSafety Rfact=0.04 Gfact=2.5 Sfact=0.6 DispFlag:1 Skin=1 Llim=1 mm
|
|
|
|
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= 600 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= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
msc: for kaon- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
msc: for mu+ SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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=225.219 keV deltaVI
|
|
BetheBloch : Emin=225.219 keV Emax= 600 MeV deltaVI
|
|
|
|
msc: for mu- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
msc: for pi+ SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
msc: for pi- SubType= 10
|
|
===== EM models for the G4Region DefaultRegionForTheWorld ======
|
|
WentzelVIUni : Emin= 0 eV Emax= 600 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 600 MeV in 160 bins
|
|
Lambda tables from threshold to 600 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= 600 MeV deltaVI
|
|
|
|
========= Table of registered couples ============================
|
|
|
|
Index : 0 used in the geometry : Yes
|
|
Material : G4_WATER
|
|
Range cuts : gamma 10 nm e- 10 nm e+ 10 nm proton 10 nm
|
|
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 eV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
==================================================================
|
|
|
|
### Run 0 starts.
|
|
|
|
-------- WWWW ------- G4Exception-START -------- WWWW -------
|
|
*** G4Exception : Analysis_W001
|
|
issued by : G4RootNtupleFileManager::SetNtupleMergingMode
|
|
Merging ntuples is not applicable in sequential application.
|
|
Setting was ignored.
|
|
*** This is just a warning message. ***
|
|
-------- WWWW -------- G4Exception-END --------- WWWW -------
|
|
|
|
... set ntuple merging row mode : row-wise - done
|
|
... create file : microtrack.root - done
|
|
... open analysis file : microtrack.root - done
|
|
... open analysis file : microtrack.root - done
|
|
Using
|
|
--> Event 0 starts.
|
|
--> Event 1 starts.
|
|
--> Event 2 starts.
|
|
--> Event 3 starts.
|
|
--> Event 4 starts.
|
|
--> Event 5 starts.
|
|
--> Event 6 starts.
|
|
--> Event 7 starts.
|
|
--> Event 8 starts.
|
|
--> Event 9 starts.
|
|
--> Event 10 starts.
|
|
--> Event 11 starts.
|
|
--> Event 12 starts.
|
|
--> Event 13 starts.
|
|
--> Event 14 starts.
|
|
--> Event 15 starts.
|
|
--> Event 16 starts.
|
|
--> Event 17 starts.
|
|
--> Event 18 starts.
|
|
--> Event 19 starts.
|
|
--> Event 20 starts.
|
|
--> Event 21 starts.
|
|
--> Event 22 starts.
|
|
--> Event 23 starts.
|
|
--> Event 24 starts.
|
|
--> Event 25 starts.
|
|
--> Event 26 starts.
|
|
--> Event 27 starts.
|
|
--> Event 28 starts.
|
|
--> Event 29 starts.
|
|
--> Event 30 starts.
|
|
--> Event 31 starts.
|
|
--> Event 32 starts.
|
|
--> Event 33 starts.
|
|
--> Event 34 starts.
|
|
--> Event 35 starts.
|
|
--> Event 36 starts.
|
|
--> Event 37 starts.
|
|
--> Event 38 starts.
|
|
--> Event 39 starts.
|
|
--> Event 40 starts.
|
|
--> Event 41 starts.
|
|
--> Event 42 starts.
|
|
--> Event 43 starts.
|
|
--> Event 44 starts.
|
|
--> Event 45 starts.
|
|
--> Event 46 starts.
|
|
--> Event 47 starts.
|
|
--> Event 48 starts.
|
|
--> Event 49 starts.
|
|
Run terminated.
|
|
Run Summary
|
|
Number of events processed : 50
|
|
User=33.610000s Real=33.746143s Sys=0.000000s
|
|
|
|
----> print histogram statistics for the entire run:
|
|
|
|
Single-event energy imparted:
|
|
----------------------------
|
|
|
|
Frequency-mean: \varepsilon_{1,F} = 1.58895 keV (rms = 1.53633 keV)
|
|
Dose-mean: \varepsilon_{1,D} = 3.07441 keV (rms = 1.952 keV)
|
|
Squared-weighted histogram: mean = 4.31377 keV (rms = 1.83416 keV)
|
|
|
|
Lineal energy:
|
|
-------------
|
|
|
|
Frequency-mean: y_F = 2.38343 keV/um (rms = 2.3045 keV/um)
|
|
Dose-mean: y_D = 4.61161 keV/um (rms = 2.928 keV/um)
|
|
Squared-weighted histogram: mean = 6.47066 keV/um (rms = 2.75125 keV/um)
|
|
|
|
Single-event specific energy:
|
|
----------------------------
|
|
|
|
Frequency-mean: z_{1,F} = 0.486208 Gy (rms = 0.470107 Gy)
|
|
Dose-mean: z_{1,D} = 0.940747 Gy (rms = 0.5973 Gy)
|
|
Squared-weighted histogram: mean = 1.31999 Gy (rms = 0.561242 Gy)
|
|
|
|
Number of hits per event:
|
|
------------------------
|
|
|
|
Eligible for site random placement, N_{sel}: mean = 4078.48 (rms = 1871.23)
|
|
Within the site, N_{site}: mean = 977.08 (rms = 680.095)
|
|
Within the site and eligible for site random placement, N_{int}: mean = 934.9 (rms = 667.453)
|
|
|
|
Kinetic energy of the primary particle:
|
|
--------------------------------------
|
|
|
|
At entrance of SDbox, T_{in}: mean = 20.0045 MeV (rms = 19.3827 keV)
|
|
At exit of SDbox, T_{out}: mean = 19.9891 MeV (rms = 19.3811 keV)
|
|
|
|
... write file : microtrack.root - done
|
|
... close file : microtrack.root - done
|
|
================== Deleting memory pools ===================
|
|
Number of memory pools allocated: 12 of which, static: 0
|
|
Dynamic pools deleted: 12 / Total memory freed: 2.5 MB
|
|
============================================================
|