405 lines
20 KiB
Plaintext
405 lines
20 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/
|
|
**************************************************************
|
|
|
|
/control/cout/ignoreThreadsExcept 0
|
|
/run/numberOfThreads 2
|
|
*** /run/numberOfThreads command is issued in sequential mode.
|
|
Command is ignored.
|
|
/run/verbose 1
|
|
#
|
|
/testhadr/det/setMat Molybdenum98
|
|
#
|
|
/process/had/verbose 1
|
|
#
|
|
/run/initialize
|
|
|
|
The Box is 10 m of Molybdenum98
|
|
|
|
Material: Molybdenum98 density: 10.280 g/cm3 RadL: 9.731 mm Nucl.Int.Length: 15.692 cm
|
|
Imean: 424.000 eV temperature: 293.15 K pressure: 1.00 atm
|
|
|
|
---> Element: Molybdenum98 (Mo98) Z = 42.0 N = 98 A = 97.905 g/mole
|
|
---> Isotope: Mo98 Z = 42 N = 98 A = 97.91 g/mole abundance: 100.000 %
|
|
ElmMassFraction: 100.00 % ElmAbundance 100.00 %
|
|
|
|
|
|
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
|
|
#
|
|
/gun/particle proton
|
|
/gun/energy 10 MeV
|
|
#
|
|
/run/printProgress 1000
|
|
/run/beamOn 10000
|
|
|
|
=======================================================
|
|
====== ParticleHP Physics Parameters ========
|
|
=======================================================
|
|
Use only photo-evaporation 1
|
|
Skip missing isotopes 0
|
|
Neglect Doppler 0
|
|
Do not adjust final state 1
|
|
Produce fission fragments 1
|
|
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
|
|
-----------------------------------------------------------------------
|
|
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: 0 eV ---> 100 TeV
|
|
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 : Molybdenum98
|
|
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
|
|
Energy thresholds : gamma 42.9626 keV e- 1.44668 MeV e+ 1.35953 MeV proton 0 eV
|
|
Region(s) which use this couple :
|
|
DefaultRegionForTheWorld
|
|
|
|
==================================================================
|
|
|
|
### Run 0 starts.
|
|
|
|
------- MixMaxRng engine status -------
|
|
Current state vector is:
|
|
mixmax state, file version 1.0
|
|
N=17 V[N]={906770732717044781, 629165745432651234, 1235682547346241386, 68420008233404568, 2262190991329242458, 2266470399991071809, 1976726662926872232, 245458862506414172, 1955974201201518530, 2155248512522080758, 604170912935414061, 1116171330120743511, 1861018313684488333, 1296715403254578286, 1549011045957234151, 370819759640195970, 2230139271784837643} counter= 17sumtot= 1977567618660788324
|
|
---------------------------------------
|
|
--> 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=0.050000s Real=0.053883s Sys=0.000000s
|
|
|
|
The run is 10000 proton of 10 MeV through 10 m of Molybdenum98 (density: 10.28 g/cm3 )
|
|
|
|
Process calls frequency:
|
|
hadElastic= 8041 protonInelastic= 1959
|
|
|
|
|
|
MeanFreePath: 4.983 cm +- 4.8835 cm massic: 51.225 g/cm2
|
|
CrossSection: 0.20068 cm^-1 massic: 1.9522 mm2/g
|
|
crossSection per atom: 3.1738 barn
|
|
|
|
Verification: crossSections from G4HadronicProcessStore
|
|
|
|
hadElastic = 1.5946 mm2/g 2.5924 barn
|
|
protonInelastic = 371.67 um2/mg 604.25 mbarn
|
|
total = 1.9663 mm2/g 3.1967 barn
|
|
|
|
List of nuclear reactions:
|
|
|
|
proton + Mo98 --> N gamma or e- + Tc99: 22 Q = 6.4945 MeV
|
|
proton + Mo98 --> N gamma or e- + Tc99[142.68360]: 1 Q = 6.3518 MeV
|
|
proton + Mo98 --> N gamma or e- + neutron + Tc98: 1776 Q = -2.4725 MeV
|
|
proton + Mo98 --> N gamma or e- + neutron + Tc98[138.59000]: 120 Q = -2.6111 MeV
|
|
proton + Mo98 --> N gamma or e- + neutron + Tc98[90.77000]: 10 Q = -2.5633 MeV
|
|
proton + Mo98 --> N gamma or e- + proton + Mo98: 27 Q = -0.042596 meV
|
|
proton + Mo98 --> N gamma or e- + proton + Mo98[734.75000]: 3 Q = -734.75 keV
|
|
proton + Mo98 --> proton + Mo98: 8041 Q = -0.0057252 meV
|
|
|
|
number of gamma or e- (ic): N = 1 --> 12
|
|
|
|
List of generated particles:
|
|
Mo98: 8068 Emean = 33.183 keV ( 1.9562 eV --> 401.7 keV)
|
|
Mo98[734.75000]: 3 Emean = 91.773 keV ( 57.063 keV --> 144.47 keV)
|
|
Tc98: 1776 Emean = 118.67 keV ( 12.565 keV --> 295.23 keV)
|
|
Tc98[138.59000]: 120 Emean = 117.01 keV ( 23.062 keV --> 276.46 keV)
|
|
Tc98[90.77000]: 10 Emean = 143.27 keV ( 47.602 keV --> 214.71 keV)
|
|
Tc99: 22 Emean = 105.24 keV ( 91.697 keV --> 118.68 keV)
|
|
Tc99[142.68360]: 1 Emean = 96.941 keV ( 96.941 keV --> 96.941 keV)
|
|
gamma: 11380 Emean = 990.61 keV ( 7.9307 keV --> 12.073 MeV)
|
|
neutron: 1906 Emean = 1.7323 MeV ( 104.2 keV --> 7.007 MeV)
|
|
proton: 30 Emean = 6.643 MeV ( 3.6611 MeV --> 9.5463 MeV)
|
|
|
|
Momentum balance: Pmean = 933.25 meV ( 0.00078225 meV --> 60.587 eV )
|
|
|
|
|
|
------- MixMaxRng engine status -------
|
|
Current state vector is:
|
|
mixmax state, file version 1.0
|
|
N=17 V[N]={1754264211422920542, 1082955125318936279, 402599475922890425, 146407241558192287, 1225671683646015187, 2098918383306110564, 960188014190822303, 1565495293656637349, 1956609809941606200, 1116121548833418844, 2046196586284439102, 1313911472640651208, 125255053698066180, 229358084532735297, 1889518433785165640, 349423770074452278, 1929443708649348895} counter= 6sumtot= 1745593823752856972
|
|
---------------------------------------
|
|
================== Deleting memory pools ===================
|
|
Number of memory pools allocated: 11 of which, static: 0
|
|
Dynamic pools deleted: 11 / Total memory freed: 0.021 MB
|
|
============================================================
|