Files
2026-03-19 16:51:22 +01:00

443 lines
23 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
#
/testhadr/det/setIsotopeMat Li7 3 7 1.85 g/cm3
#
/testhadr/det/setNbOfAbsor 1
/run/reinitializeGeometry
/testhadr/det/setAbsor 1 Li7 30 cm
/run/reinitializeGeometry
/testhadr/det/setSizeYZ 30 cm
/run/reinitializeGeometry
#
/process/em/verbose 0
/process/had/verbose 1
#
/run/initialize
-------------------------------------------------------------
---> The Absorber is 1 layers of:
Li7: 30 cm
-------------------------------------------------------------
hInelastic QGSP_BIC_HP Thresholds:
0) between BIC and BERT for p, n over the interval 1 to 1.5 GeV.
1) between BERT and FTF/P over the interval 3 to 6 GeV.
2) between FTF/P and QGS/P over the interval 12 to 25 GeV.
-- quasiElastic: 1 for QGS and 0 for FTF
@@@ G4ParticleHPInelasticData instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1
NeutronHP: /Capture file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Capture/CrossSection/18_37_Argon
NeutronHP: /Elastic file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_37_Argon
NeutronHP: /Inelastic file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Inelastic/CrossSection/18_37_Argon
NeutronHP: /Capture file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Capture/CrossSection/18_39_Argon
NeutronHP: /Elastic file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_39_Argon
NeutronHP: /Inelastic file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Inelastic/CrossSection/18_39_Argon
#
/testhadr/gun/setDefault
/gun/particle neutron
/gun/energy 14.1 MeV
#
/run/printProgress 1000
/run/beamOn 10000
======================================================================
====== Radioactive Decay Physics Parameters =======
======================================================================
min MeanLife (from G4NuclideTable) 1 ns
Max life time (from G4DeexPrecoParameters) 1 ns
Internal e- conversion flag 1
Stored internal conversion coefficients 1
Enabled atomic relaxation mode 0
Enable correlated gamma emission 0
Max 2J for sampling of angular correlations 10
Atomic de-excitation enabled 0
Auger electron emission enabled 0
Check EM cuts disabled for atomic de-excitation 0
Use Bearden atomic level energies 0
Use ANSTO fluorescence model 0
Threshold for very long decay time at rest 1 y
======================================================================
NeutronHP: /Capture file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Capture/CrossSection/18_37_Argon
NeutronHP: /Elastic file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_37_Argon
NeutronHP: /Inelastic file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Inelastic/CrossSection/18_37_Argon
NeutronHP: /Capture file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Capture/CrossSection/18_39_Argon
NeutronHP: /Elastic file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_39_Argon
NeutronHP: /Inelastic file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Inelastic/CrossSection/18_39_Argon
NeutronHP: /Elastic file for Z = 18, A = 36 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_37_Argon
NeutronHP: /Elastic file for Z = 18, A = 38 is not found and NeutronHP will use /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Elastic/CrossSection/18_39_Argon
=======================================================
====== ParticleHP Physics Parameters ========
=======================================================
Use only photo-evaporation 1
Skip missing isotopes 0
Neglect Doppler 0
Do not adjust final state 1
Produce fission fragments 0
Use WendtFissionModel 0
Use NRESP71Model 1
Use DBRC 0
PHP use Poisson 0
PHP check 1
CHECK HP NAMES 0
Enable DEBUG 0
Use probability tables from
=======================================================
@@@ G4ParticleHPInelastic instantiated for particle neutron/n data directory is /cvmfs/geant4.cern.ch/share/data/G4NDL4.7.1/Inelastic
====================================================================
HADRONIC PROCESSES SUMMARY (verbose level 1)
-----------------------------------------------------------------------
Hadronic Processes for GenericIon
Process: ionElastic
Model: NNDiffuseElastic: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: ionInelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: RadioactiveDecay
-----------------------------------------------------------------------
Hadronic Processes for He3
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: He3Inelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: He3ParticleXS: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for alpha
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: alphaInelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: alphaParticleXS: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_He3
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_He3Inelastic
Model: FTFP: 0 eV /n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_alpha
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_alphaInelastic
Model: FTFP: 0 eV /n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_deuteron
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_deuteronInelastic
Model: FTFP: 0 eV /n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_lambda
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: anti_lambdaInelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 0 eV ---> 25 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for anti_neutron
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100.1 MeV
Model: AntiAElastic: 100 MeV ---> 100 TeV
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_neutronInelastic
Model: FTFP: 0 eV ---> 100 TeV
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_proton
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100.1 MeV
Model: AntiAElastic: 100 MeV ---> 100 TeV
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_protonInelastic
Model: FTFP: 0 eV ---> 100 TeV
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for anti_triton
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100.1 MeV/n
Model: AntiAElastic: 100 MeV/n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
Process: anti_tritonInelastic
Model: FTFP: 0 eV /n ---> 100 TeV/n
Cr_sctns: AntiAGlauber: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for deuteron
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: dInelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: deuteronParticleXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for gamma
Process: photonNuclear
Model: GammaNPreco: 0 eV ---> 200 MeV
Model: BertiniCascade: 199 MeV ---> 10 GeV
Cr_sctns: GammaNuclearXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for kaon+
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: kaon+Inelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for kaon-
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: kaon-Inelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for lambda
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: lambdaInelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for neutron
Process: hadElastic
Model: hElasticCHIPS: 19.5 MeV ---> 100 TeV
Model: NeutronHPElastic: 0 eV ---> 20 MeV
Cr_sctns: NeutronHPElasticXS: 0 eV ---> 20 MeV
Cr_sctns: G4NeutronElasticXS: 0 eV ---> 100 TeV
Process: neutronInelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 1 GeV ---> 6 GeV
Model: Binary Cascade: 19.9 MeV ---> 1.5 GeV
Model: NeutronHPInelastic: 0 eV ---> 20 MeV
Cr_sctns: NeutronHPInelasticXS: 0 eV ---> 20 MeV
Cr_sctns: G4NeutronInelasticXS: 0 eV ---> 100 TeV
Process: nCaptureHP
Model: nRadCaptureHP: 0 eV ---> 100 TeV
Cr_sctns: NeutronHPCaptureXS: 0 eV ---> 100 TeV
Cr_sctns: G4NeutronCaptureXS: 0 eV ---> 100 TeV
Process: nFissionHP
Model: NeutronHPFission: 0 eV ---> 20 MeV
Cr_sctns: NeutronHPFissionXS: 0 eV ---> 20 MeV
Cr_sctns: ZeroXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for pi+
Process: hadElastic
Model: hElasticGlauber: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: pi+Inelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for pi-
Process: hadElastic
Model: hElasticGlauber: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: pi-Inelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for proton
Process: hadElastic
Model: hElasticCHIPS: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: protonInelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 1 GeV ---> 6 GeV
Model: Binary Cascade: 0 eV ---> 1.5 GeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for sigma-
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: sigma-Inelastic
Model: QGSP: 12 GeV ---> 100 TeV
Model: FTFP: 3 GeV ---> 25 GeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for triton
Process: hadElastic
Model: hElasticLHEP: 0 eV /n ---> 100 TeV/n
Cr_sctns: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
Process: tInelastic
Model: Binary Light Ion Cascade: 0 eV /n ---> 6 GeV/n
Model: FTFP: 3 GeV/n ---> 100 TeV/n
Cr_sctns: tritonParticleXS: 0 eV ---> 25.6 PeV
=======================================================================
====== Geant4 Native Pre-compound Model Parameters ========
=======================================================================
Type of pre-compound model 0
Type of pre-compound inverse x-section 1
Pre-compound model active 1
Pre-compound excitation low energy 0.1 MeV
Pre-compound excitation high energy 15 MeV
Angular generator for pre-compound model 1
Use NeverGoBack option for pre-compound model 0
Use SoftCutOff option for pre-compound model 0
Use CEM transitions for pre-compound model 1
Use GNASH transitions for pre-compound model 0
Use HETC submodel for pre-compound model 0
=======================================================================
====== Nuclear De-excitation Module Parameters ========
=======================================================================
Type of de-excitation inverse x-section 3
Type of de-excitation factory Evaporation+GEM
Number of de-excitation channels 68
Type of Fermi BreakUp model ModelVI
Min excitation energy 0.01 keV
Min energy per nucleon for multifragmentation 2e+05 MeV
Level density (1/MeV) 0.075
Use simple level density model 1
Use discrete excitation energy of the residual 0
Time limit for long lived isomeres 1 ns
Isomer production flag 1
Internal e- conversion flag 1
Store e- internal conversion data 1
Correlated gamma emission flag 0
Max 2J for sampling of angular correlations 10
=======================================================================
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : Galactic
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : Li7
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
Energy thresholds : gamma 1.50862 keV e- 466.359 keV e+ 452.491 keV proton 0 eV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 0 starts.
--------- Ranecu engine status ---------
Initial seed (index) = 0
Current couple of seeds = 9876, 54321
----------------------------------------
--> Event 0 starts.
--> Event 1000 starts.
--> Event 2000 starts.
--> Event 3000 starts.
--> Event 4000 starts.
--> Event 5000 starts.
--> Event 6000 starts.
--> Event 7000 starts.
--> Event 8000 starts.
--> Event 9000 starts.
Run terminated.
Run Summary
Number of events processed : 10000
User=2.350000s Real=2.799897s Sys=0.000000s
======================== run summary =====================
The run is 10000 neutron of 14.1 MeV through 1 absorbers:
1 30 cm of Li7 (density: 1.85 g/cm3 )
Process calls frequency :
Decay= 20 NoProcess= 99239 RadioactiveDecay= 171
Rayl= 1570 Transportation= 58522 annihil= 138
compt= 52807 conv= 143 eBrem= 734
eIoni= 108323 hIoni= 197 hadElastic= 85875
ionIoni= 171674 msc= 75 nCaptureHP= 4
neutronInelastic= 14215 phot= 783 photonNuclear= 20
Edep in absorber 1 = 6.93 MeV (315 eV -->22.2 MeV)
Energy leakage = 4.65 MeV (157 eV -->15.7 MeV)
Energy total = 11.6 MeV (347 keV-->32.2 MeV)
List of created particles in absorber 1:
Be8[3030.00000]: 4 Emean = 7.63 keV ( 3.18 keV --> 10.8 keV) mean life = 0 ps
He6: 163 Emean = 1.27 MeV ( 866 keV --> 1.34 MeV) mean life = 1.16 s
Li6: 163 Emean = 777 eV ( 27.1 eV --> 1.4 keV) stable
Li7: 99076 Emean = 564 keV ( 1.79 meV --> 6.22 MeV) stable
Li8: 4 Emean = 245 keV ( 3.85 keV --> 488 keV) mean life = 1.21 s
alpha: 339 Emean = 1.18 MeV ( 28.1 keV --> 4.02 MeV) stable
anti_nu_e: 167 Emean = 2.11 MeV ( 233 keV --> 9.81 MeV) stable
deuteron: 163 Emean = 423 keV ( 290 keV --> 447 keV) stable
e+: 143 Emean = 1.85 MeV ( 153 keV --> 3.53 MeV) stable
e-: 53513 Emean = 245 keV ( 100 eV --> 9.4 MeV) stable
gamma: 14190 Emean = 2.48 MeV ( 1.52 keV --> 4.89 MeV) stable
neutron: 14903 Emean = 3.66 MeV ( 5.65 keV --> 13.6 MeV) mean life = 14.7 min
triton: 20 Emean = 1.24 MeV ( 1.18 MeV --> 1.29 MeV) mean life = 17.8 y
List of particles emerging from absorbers :
anti_nu_e: 167 Emean = 2.11 MeV ( 233 keV --> 9.81 MeV)
e+: 6 Emean = 1.93 MeV ( 794 keV --> 3.25 MeV)
e-: 159 Emean = 1.84 MeV ( 13.3 keV --> 4.31 MeV)
gamma: 13244 Emean = 1.68 MeV ( 4.43 keV --> 4.63 MeV)
neutron: 10684 Emean = 2.21 MeV ( 129 meV --> 14.1 MeV)
Nb of events with primary absorbed = 87 %, transmit = 7 %, reflected = 6.4 %
--------- Ranecu engine status ---------
Initial seed (index) = 0
Current couple of seeds = 974953601, 2012462375
----------------------------------------
================== Deleting memory pools ===================
Number of memory pools allocated: 11 of which, static: 0
Dynamic pools deleted: 11 / Total memory freed: 0.041 MB
============================================================