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

2902 lines
115 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/
**************************************************************
<<< Geant4 Physics List simulation engine: FTFP_BERT
G4VModularPhysicsList::ReplacePhysics: G4EmStandard with type : 2 is replaced with G4EmStandard_opt4
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.
/tracking/verbose 0
/run/verbose 1
/process/optical/verbose 1
/run/physicsModified
/control/cout/ignoreThreadsExcept 0
/opnovice2/boxProperty RINDEX 0.000002 1.3 0.000008 1.4
The MPT for the box is now:
0: RINDEX
2e-06 1.3
8e-06 1.4
9: GROUPVEL
2e-06 218.486
8e-06 203.645
.............
/opnovice2/boxProperty ABSLENGTH 0.000002 1000000 0.000005 2000000 0.000008 3000000
The MPT for the box is now:
0: RINDEX
2e-06 1.3
8e-06 1.4
9: GROUPVEL
2e-06 218.486
8e-06 203.645
16: ABSLENGTH
2e-06 1e+06
5e-06 2e+06
8e-06 3e+06
.............
/opnovice2/worldProperty RINDEX 0.000002 1.01 0.000008 1.01
The MPT for the world is now:
0: RINDEX
2e-06 1.01
8e-06 1.01
9: GROUPVEL
2e-06 296.824
8e-06 296.824
.............
/opnovice2/worldProperty ABSLENGTH 0.000002 1000000 0.000005 2000000 0.000008 3000000
The MPT for the world is now:
0: RINDEX
2e-06 1.01
8e-06 1.01
9: GROUPVEL
2e-06 296.824
8e-06 296.824
16: ABSLENGTH
2e-06 1e+06
5e-06 2e+06
8e-06 3e+06
.............
/opnovice2/surfaceProperty SPECULARLOBECONSTANT 0.000002 0.1 0.000008 0.1
SPECULARLOBECONSTANT 0.000002 0.1 0.000008 0.1
The MPT for the surface is now:
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty SPECULARSPIKECONSTANT 0.000002 0.1 0.000008 0.1
SPECULARSPIKECONSTANT 0.000002 0.1 0.000008 0.1
The MPT for the surface is now:
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty BACKSCATTERCONSTANT 0.000002 0.1 0.000008 0.1
BACKSCATTERCONSTANT 0.000002 0.1 0.000008 0.1
The MPT for the surface is now:
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty TRANSMITTANCE 0.000002 0.1 0.000008 0.1
TRANSMITTANCE 0.000002 0.1 0.000008 0.1
The MPT for the surface is now:
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty REFLECTIVITY 0.000002 0.8 0.000008 0.8
REFLECTIVITY 0.000002 0.8 0.000008 0.8
The MPT for the surface is now:
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty EFFICIENCY 0.000002 0.1 0.000008 0.1
EFFICIENCY 0.000002 0.1 0.000008 0.1
The MPT for the surface is now:
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
.............
/opnovice2/surfaceProperty RINDEX 0.000002 1.4 0.000008 1.5
RINDEX 0.000002 1.4 0.000008 1.5
The MPT for the surface is now:
0: RINDEX
2e-06 1.4
8e-06 1.5
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
9: GROUPVEL
2e-06 203.645
8e-06 190.691
.............
/opnovice2/surfaceSigmaAlpha 0.2
/run/geometryModified
Surface sigma alpha set to: 0.2
/run/initialize
****** opticalSurface->DumpInfo:
Surface type = 1
Surface finish = 3
Surface model = 1
Surface parameter
-----------------
sigma_alpha: 0.2
****** end of opticalSurface->DumpInfo
hInelastic FTFP_BERT : threshold between BERT and FTFP is over the interval
for pions : 3 to 6 GeV
for kaons : 3 to 6 GeV
for proton : 3 to 6 GeV
for neutron : 3 to 6 GeV
### Adding tracking cuts for neutron TimeCut(ns)= 10000 KinEnergyCut(MeV)= 0
OpAbsorption is created
OpRayleigh is created
OpMieHG is created
OpBoundary is created
OpWLS is created
OpWLS2 is created
Cerenkov is created.
### Birks coefficients used in run time
G4_WATER 0.126 mm/MeV 0.0126 g/cm^2/MeV massFactor= 640.135 effCharge= 22
/opnovice2/stepping/killOnSecondSurface true
#
/gun/particle opticalphoton
/gun/energy 3 eV
/gun/position 0 0 0 cm
# will be normalized
/gun/direction .9 .44 0
/opnovice2/gun/optPhotonPolar
#
/analysis/setFileName unified
/analysis/h1/set 10 40 -1 39
/analysis/h1/set 11 100 -1.1 1.1
/analysis/h1/set 12 100 -1.1 1.1
/analysis/h1/set 13 100 -1.1 1.1
/analysis/h1/set 14 100 -1.1 1.1
/analysis/h1/set 15 100 -1.1 1.1
/analysis/h1/set 16 100 -1.1 1.1
/analysis/h1/set 17 100 -1.1 1.1
/analysis/h1/set 18 100 -1.1 1.1
/analysis/h1/set 19 100 -1.1 1.1
###############################################################################
################################# UNIFIED #####################################
###############################################################################
/opnovice2/surfaceModel unified
/run/geometryModified
/opnovice2/surfaceType dielectric_dielectric
/run/geometryModified
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ground ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish ground
/run/geometryModified
## unified, davis, lut use sigma_alpha glisur uses polish
/run/beamOn 10000
=======================================================================
====== 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 1
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 1
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 1
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 0
PIXE atomic de-excitation enabled 0
De-excitation module ignores cuts 0
Type of PIXE cross section for hadrons Empirical
Type of PIXE cross section for e+- Livermore
=======================================================================
### === Deexcitation model UAtomDeexcitation is activated for 1 region:
DefaultRegionForTheWorld 1 0 0
### === Ignore cuts flag: 0
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 ======
LowEPComptonModel : Emin= 0 eV Emax= 20 MeV Fluo
KleinNishina : Emin= 20 MeV 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 ======
BetheHeitler5D : 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 ======
PenIoni : Emin= 0 eV Emax= 100 keV
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 ======
PenIoni : Emin= 0 eV Emax= 100 keV
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
annihil: for e+ XStype:2 SubType=5 AtRestModel:Allison BuildTable=0
===== EM models for the G4Region DefaultRegionForTheWorld ======
eplusTo2or3gamma : 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
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
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
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
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
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
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
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
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
====================================================================
HADRONIC PROCESSES SUMMARY (verbose level 1)
-----------------------------------------------------------------------
Hadronic Processes for neutron
Process: hadElastic
Model: hElasticCHIPS: 0 eV ---> 100 TeV
Cr_sctns: G4NeutronElasticXS: 0 eV ---> 100 TeV
Process: neutronInelastic
Model: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: G4NeutronInelasticXS: 0 eV ---> 100 TeV
Process: nCapture
Model: nRadCapture: 0 eV ---> 100 TeV
Cr_sctns: G4NeutronCaptureXS: 0 eV ---> 100 TeV
Process: nKiller
-----------------------------------------------------------------------
Hadronic Processes for B-
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: B-Inelastic
Model: FTFP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for D-
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: D-Inelastic
Model: FTFP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for GenericIon
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: Glauber-Gribov Nucl-nucl: 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: Glauber-Gribov Nucl-nucl: 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
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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
Process: hFritiofCaptureAtRest
-------------------------------------------------------------------------
Hadronic Processes for anti_hypertriton
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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: FTFP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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
Process: hFritiofCaptureAtRest
-----------------------------------------------------------------------
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: Glauber-Gribov Nucl-nucl: 0 eV ---> 25.6 PeV
-----------------------------------------------------------------------
Hadronic Processes for e+
Process: positronNuclear
Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV
Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for e-
Process: electronNuclear
Model: G4ElectroVDNuclearModel: 0 eV ---> 1 PeV
Cr_sctns: ElectroNuclearXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for gamma
Process: photonNuclear
Model: GammaNPreco: 0 eV ---> 200 MeV
Model: BertiniCascade: 199 MeV ---> 6 GeV
Model: TheoFSGenerator: 3 GeV ---> 100 TeV
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: FTFP: 3 GeV ---> 100 TeV
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: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
Hadronic Processes for lambda
Process: hadElastic
Model: hElasticLHEP: 0 eV ---> 100 TeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: lambdaInelastic
Model: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for mu+
Process: muonNuclear
Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV
Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV
-----------------------------------------------------------------------
Hadronic Processes for mu-
Process: muonNuclear
Model: G4MuonVDNuclearModel: 0 eV ---> 1 PeV
Cr_sctns: KokoulinMuonNuclearXS: 0 eV ---> 100 TeV
Process: muMinusCaptureAtRest
-----------------------------------------------------------------------
Hadronic Processes for pi+
Process: hadElastic
Model: hElasticGlauber: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: pi+Inelastic
Model: FTFP: 3 GeV ---> 100 TeV
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: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
Hadronic Processes for proton
Process: hadElastic
Model: hElasticCHIPS: 0 eV ---> 100 TeV
Cr_sctns: BarashenkovGlauberGribov: 0 eV ---> 100 TeV
Process: protonInelastic
Model: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 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: FTFP: 3 GeV ---> 100 TeV
Model: BertiniCascade: 0 eV ---> 6 GeV
Cr_sctns: Glauber-Gribov: 0 eV ---> 100 TeV
Process: hBertiniCaptureAtRest
-----------------------------------------------------------------------
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: Glauber-Gribov Nucl-nucl: 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 : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 0 starts.
... set ntuple merging row mode : row-wise - done
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.041773s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0078
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 7
Total # of surface events: 9993
Unaccounted for: 0
Surface events by process:
Transmission: 979
Fresnel refraction: 7508
Lambertian reflection: 319
Lobe reflection: 54
Spike reflection: 52
Backscattering: 55
Absorption: 933
Detection: 93
Sum: 9993
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.687721 rms = 0.256585
Y momentum dir of backward-going photons: mean = 0.0921181 rms = 0.538817
Z momentum dir of backward-going photons: mean = -0.00593677 rms = 0.402938
X momentum dir of forward-going photons: mean = 0.818856 rms = 0.0611041
Y momentum dir of forward-going photons: mean = 0.560325 rms = 0.0805691
Z momentum dir of forward-going photons: mean = 0.00111875 rms = 0.0726935
X momentum dir of Fresnel-refracted photons: mean = 0.808486 rms = 0.0573439
Y momentum dir of Fresnel-refracted photons: mean = 0.576118 rms = 0.0719422
Z momentum dir of Fresnel-refracted photons: mean = 0.00126463 rms = 0.0772865
... write file : unified.root - done
... close file : unified.root - done
# restore default
#/opnovice2/surfaceSigmaAlpha 0.0
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ polished ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish polished
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 1 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.032070s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.008
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 11
Total # of surface events: 9989
Unaccounted for: 0
Surface events by process:
Transmission: 1003
Fresnel refraction: 7803
Fresnel reflection: 178
Absorption: 888
Detection: 117
Sum: 9989
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 2.35608e-08
Y momentum dir of backward-going photons: mean = 0.43921 rms = 2.23517e-08
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.828802 rms = 0.0249472
Y momentum dir of forward-going photons: mean = 0.557378 rms = 0.0423663
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0.819857 rms = 2.24261e-07
Y momentum dir of Fresnel-refracted photons: mean = 0.572568 rms = 2.71716e-07
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
/opnovice2/surfaceConstProperty SURFACEROUGHNESS 0.01
The MPT for the surface is now:
0: RINDEX
2e-06 1.4
8e-06 1.5
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
9: GROUPVEL
2e-06 203.645
8e-06 190.691
0: SURFACEROUGHNESS 0.01
.............
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 2 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.031449s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0075
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 6
Total # of surface events: 9994
Unaccounted for: 0
Surface events by process:
Transmission: 984
Fresnel refraction: 7855
Lambertian reflection: 167
Absorption: 883
Detection: 105
Sum: 9994
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.679312 rms = 0.236721
Y momentum dir of backward-going photons: mean = 0.0382186 rms = 0.492826
Z momentum dir of backward-going photons: mean = -0.0355269 rms = 0.486721
X momentum dir of forward-going photons: mean = 0.828599 rms = 0.0246994
Y momentum dir of forward-going photons: mean = 0.557722 rms = 0.0419455
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0.819857 rms = 2.24261e-07
Y momentum dir of Fresnel-refracted photons: mean = 0.572568 rms = 2.72023e-07
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
# restore default
/opnovice2/surfaceConstProperty SURFACEROUGHNESS 0.0
The MPT for the surface is now:
0: RINDEX
2e-06 1.4
8e-06 1.5
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
9: GROUPVEL
2e-06 203.645
8e-06 190.691
0: SURFACEROUGHNESS 0
.............
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ painted ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish polishedfrontpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 3 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.033662s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0019
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 8
Total # of surface events: 9992
Unaccounted for: 0
Surface events by process:
Transmission: 1049
Spike reflection: 7929
Absorption: 910
Detection: 104
Sum: 9992
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 6.34132e-07
Y momentum dir of backward-going photons: mean = 0.43921 rms = 1.71929e-07
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.78503e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.56007e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
/opnovice2/surfaceFinish groundfrontpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 4 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.036911s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0025
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 8
Total # of surface events: 9992
Unaccounted for: 0
Surface events by process:
Transmission: 1006
Lambertian reflection: 7974
Absorption: 907
Detection: 105
Sum: 9992
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.667769 rms = 0.234274
Y momentum dir of backward-going photons: mean = 0.001713 rms = 0.498209
Z momentum dir of backward-going photons: mean = -0.000221894 rms = 0.500984
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.77568e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.39638e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
/opnovice2/surfaceFinish polishedbackpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 5 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.070000s Real=0.069120s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0026
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 6
Total # of surface events: 9994
Unaccounted for: 0
Surface events by process:
Transmission: 974
Fresnel refraction: 8054
Lambertian reflection: 13
Lobe reflection: 3
Backscattering: 1
Absorption: 835
Detection: 114
Sum: 9994
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.896632 rms = 0.0233376
Y momentum dir of backward-going photons: mean = 0.437618 rms = 0.0482827
Z momentum dir of backward-going photons: mean = 0.000491427 rms = 0.0408103
X momentum dir of forward-going photons: mean = 0.912451 rms = 0.0288346
Y momentum dir of forward-going photons: mean = 0.405314 rms = 0.0342028
Z momentum dir of forward-going photons: mean = -0.000495469 rms = 0.0339468
X momentum dir of Fresnel-refracted photons: mean = -0.896974 rms = 0.0198106
Y momentum dir of Fresnel-refracted photons: mean = 0.438779 rms = 0.036663
Z momentum dir of Fresnel-refracted photons: mean = 0.00069122 rms = 0.034239
... write file : unified.root - done
... close file : unified.root - done
/opnovice2/surfaceFinish groundbackpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 6 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.070000s Real=0.070468s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0032
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 9
Total # of surface events: 9991
Unaccounted for: 0
Surface events by process:
Transmission: 1168
Fresnel refraction: 7605
Lambertian reflection: 19
Lobe reflection: 6
Spike reflection: 3
Backscattering: 5
Absorption: 1070
Detection: 115
Sum: 9991
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.679202 rms = 0.228788
Y momentum dir of backward-going photons: mean = 0.0104892 rms = 0.49541
Z momentum dir of backward-going photons: mean = 0.00338132 rms = 0.490701
X momentum dir of forward-going photons: mean = 0.855862 rms = 0.171656
Y momentum dir of forward-going photons: mean = 0.346258 rms = 0.25087
Z momentum dir of forward-going photons: mean = -0.00593491 rms = 0.234881
X momentum dir of Fresnel-refracted photons: mean = -0.679015 rms = 0.229939
Y momentum dir of Fresnel-refracted photons: mean = 0.0106281 rms = 0.495031
Z momentum dir of Fresnel-refracted photons: mean = 0.00286933 rms = 0.490805
... write file : unified.root - done
... close file : unified.root - done
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ dielectric_metal ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceType dielectric_metal
/run/geometryModified
/opnovice2/surfaceFinish polished
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 7 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.033668s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0022
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 7
Total # of surface events: 9993
Unaccounted for: 0
Surface events by process:
Transmission: 988
Spike reflection: 7974
Absorption: 915
Detection: 116
Sum: 9993
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 6.35094e-07
Y momentum dir of backward-going photons: mean = 0.43921 rms = 1.7201e-07
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.76627e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.34665e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
/opnovice2/surfaceFinish ground
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 8 starts.
... create file : unified.root - done
... open analysis file : unified.root - done
... open analysis file : unified.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.035939s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0024
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 8
Total # of surface events: 9992
Unaccounted for: 0
Surface events by process:
Transmission: 996
Lambertian reflection: 5618
Lobe reflection: 761
Spike reflection: 826
Backscattering: 765
Absorption: 940
Detection: 86
Sum: 9992
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.725175 rms = 0.22299
Y momentum dir of backward-going photons: mean = 0.0307967 rms = 0.486614
Z momentum dir of backward-going photons: mean = -0.00355541 rms = 0.432022
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.76941e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.34665e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : unified.root - done
... close file : unified.root - done
###############################################################################
################################# GLISUR ######################################
###############################################################################
/opnovice2/surfaceModel glisur
/run/geometryModified
/opnovice2/surfaceType dielectric_dielectric
/run/geometryModified
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ polished ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish polished
/run/geometryModified
/analysis/setFileName glisur
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 9 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.030875s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.01
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 11
Total # of surface events: 9989
Unaccounted for: 0
Surface events by process:
Transmission: 980
Fresnel refraction: 7845
Fresnel reflection: 181
Absorption: 882
Detection: 101
Sum: 9989
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 2.78775e-08
Y momentum dir of backward-going photons: mean = 0.43921 rms = 2.35608e-08
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.828578 rms = 0.0246725
Y momentum dir of forward-going photons: mean = 0.557759 rms = 0.0418998
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0.819857 rms = 2.24509e-07
Y momentum dir of Fresnel-refracted photons: mean = 0.572568 rms = 2.71921e-07
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceConstProperty SURFACEROUGHNESS 0.01
The MPT for the surface is now:
0: RINDEX
2e-06 1.4
8e-06 1.5
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
9: GROUPVEL
2e-06 203.645
8e-06 190.691
0: SURFACEROUGHNESS 0.01
.............
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 10 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.031085s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.008
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 9
Total # of surface events: 9991
Unaccounted for: 0
Surface events by process:
Transmission: 1009
Fresnel refraction: 7788
Lambertian reflection: 175
Absorption: 928
Detection: 91
Sum: 9991
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.685008 rms = 0.235191
Y momentum dir of backward-going photons: mean = -0.0672448 rms = 0.461926
Z momentum dir of backward-going photons: mean = -0.0564351 rms = 0.504348
X momentum dir of forward-going photons: mean = 0.828864 rms = 0.0250232
Y momentum dir of forward-going photons: mean = 0.557272 rms = 0.0424954
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0.819857 rms = 2.24509e-07
Y momentum dir of Fresnel-refracted photons: mean = 0.572568 rms = 2.7141e-07
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceConstProperty SURFACEROUGHNESS 0.0
The MPT for the surface is now:
0: RINDEX
2e-06 1.4
8e-06 1.5
1: REFLECTIVITY
2e-06 0.8
8e-06 0.8
4: EFFICIENCY
2e-06 0.1
8e-06 0.1
5: TRANSMITTANCE
2e-06 0.1
8e-06 0.1
6: SPECULARLOBECONSTANT
2e-06 0.1
8e-06 0.1
7: SPECULARSPIKECONSTANT
2e-06 0.1
8e-06 0.1
8: BACKSCATTERCONSTANT
2e-06 0.1
8e-06 0.1
9: GROUPVEL
2e-06 203.645
8e-06 190.691
0: SURFACEROUGHNESS 0
.............
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ground ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish ground
/run/geometryModified
/run/initialize
/opnovice2/surfacePolish 0.2
/run/geometryModified
Surface polish set to: 0.2
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 11 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.033837s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0067
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 12
Total # of surface events: 9988
Unaccounted for: 0
Surface events by process:
Transmission: 1033
Fresnel refraction: 6134
Fresnel reflection: 353
Total internal reflection: 1502
Absorption: 880
Detection: 86
Sum: 9988
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.325606 rms = 0.255169
Y momentum dir of backward-going photons: mean = 0.733737 rms = 0.277671
Z momentum dir of backward-going photons: mean = -0.0290882 rms = 0.461034
X momentum dir of forward-going photons: mean = 0.822746 rms = 0.10792
Y momentum dir of forward-going photons: mean = 0.529526 rms = 0.113763
Z momentum dir of forward-going photons: mean = 0.00308688 rms = 0.134509
X momentum dir of Fresnel-refracted photons: mean = 0.797893 rms = 0.170331
Y momentum dir of Fresnel-refracted photons: mean = 0.54611 rms = 0.117491
Z momentum dir of Fresnel-refracted photons: mean = 0.0047801 rms = 0.149302
... write file : glisur.root - done
... close file : glisur.root - done
# restore default
/opnovice2/surfacePolish 1.0
/run/geometryModified
Surface polish set to: 1
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ painted ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceFinish polishedfrontpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 12 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.035129s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0034
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 14
Total # of surface events: 9986
Unaccounted for: 0
Surface events by process:
Transmission: 1011
Spike reflection: 7969
Absorption: 889
Detection: 117
Sum: 9986
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 6.34919e-07
Y momentum dir of backward-going photons: mean = 0.43921 rms = 1.7201e-07
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.77568e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.39638e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceFinish polishedbackpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 13 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.041353s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0017
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 4
Total # of surface events: 9996
Unaccounted for: 0
Surface events by process:
Transmission: 1030
Fresnel refraction: 7900
Fresnel reflection: 6
Absorption: 952
Detection: 108
Sum: 9996
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 6.33781e-07
Y momentum dir of backward-going photons: mean = 0.43921 rms = 1.71929e-07
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.912889 rms = 2.31328e-07
Y momentum dir of forward-going photons: mean = 0.408207 rms = 2.78775e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = -0.898384 rms = 6.33694e-07
Y momentum dir of Fresnel-refracted photons: mean = 0.43921 rms = 1.71929e-07
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceFinish groundfrontpainted
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 14 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.035422s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0023
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 3
Total # of surface events: 9997
Unaccounted for: 0
Surface events by process:
Transmission: 1001
Lambertian reflection: 8020
Absorption: 875
Detection: 101
Sum: 9997
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.668223 rms = 0.233789
Y momentum dir of backward-going photons: mean = -0.00364135 rms = 0.497797
Z momentum dir of backward-going photons: mean = 0.00399356 rms = 0.500988
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.77568e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.36326e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceFinish groundbackpainted
/run/geometryModified
/opnovice2/surfacePolish 0.2
/run/geometryModified
Surface polish set to: 0.2
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 15 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.050000s Real=0.046654s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.002
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 6
Total # of surface events: 9994
Unaccounted for: 0
Surface events by process:
Transmission: 1101
Fresnel refraction: 7714
Fresnel reflection: 44
Absorption: 1025
Detection: 110
Sum: 9994
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.642445 rms = 0.26517
Y momentum dir of backward-going photons: mean = 0.00669636 rms = 0.507397
Z momentum dir of backward-going photons: mean = -0.00204362 rms = 0.509361
X momentum dir of forward-going photons: mean = 0.762545 rms = 0.290058
Y momentum dir of forward-going photons: mean = 0.29248 rms = 0.376478
Z momentum dir of forward-going photons: mean = -0.00669823 rms = 0.327209
X momentum dir of Fresnel-refracted photons: mean = -0.627929 rms = 0.28776
Y momentum dir of Fresnel-refracted photons: mean = 0.00310404 rms = 0.50949
Z momentum dir of Fresnel-refracted photons: mean = -0.00275039 rms = 0.51313
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfacePolish 1.0
/run/geometryModified
Surface polish set to: 1
#
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ dielectric_metal ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#
/opnovice2/surfaceType dielectric_metal
/run/geometryModified
/opnovice2/surfaceFinish polished
/run/geometryModified
/run/initialize
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 16 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.033437s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0021
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 8
Total # of surface events: 9992
Unaccounted for: 0
Surface events by process:
Transmission: 984
Spike reflection: 8009
Absorption: 909
Detection: 90
Sum: 9992
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.898384 rms = 6.35705e-07
Y momentum dir of backward-going photons: mean = 0.43921 rms = 1.7201e-07
Z momentum dir of backward-going photons: mean = 0 rms = 0
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.76627e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.29662e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfaceFinish ground
/run/geometryModified
/run/initialize
/opnovice2/surfacePolish 0.2
/run/geometryModified
Surface polish set to: 0.2
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 17 starts.
... create file : glisur.root - done
... open analysis file : glisur.root - done
... open analysis file : glisur.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.040000s Real=0.034942s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0022
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 7
Total # of surface events: 9993
Unaccounted for: 0
Surface events by process:
Transmission: 1013
Lobe reflection: 7950
Absorption: 928
Detection: 102
Sum: 9993
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.57728 rms = 0.27638
Y momentum dir of backward-going photons: mean = 0.308183 rms = 0.493476
Z momentum dir of backward-going photons: mean = -0.00421615 rms = 0.501845
X momentum dir of forward-going photons: mean = 0.898384 rms = 1.77568e-07
Y momentum dir of forward-going photons: mean = 0.43921 rms = 8.42937e-08
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : glisur.root - done
... close file : glisur.root - done
/opnovice2/surfacePolish 1.0
/run/geometryModified
Surface polish set to: 1
###############################################################################
################################# DAVIS LUT ###################################
###############################################################################
/opnovice2/surfaceModel DAVIS
/run/geometryModified
/opnovice2/surfaceType dielectric_LUTDAVIS
/run/geometryModified
/opnovice2/surfaceFinish RoughTeflon_LUT
G4OpticalSurface: data file /cvmfs/geant4.cern.ch/share/data/RealSurface2.2/RoughTeflon_LUT.z successfully read in.
LUT DAVIS - data file: RoughTeflon_LUT.z read in!
G4OpticalSurface: data file /cvmfs/geant4.cern.ch/share/data/RealSurface2.2/RoughTeflon_LUTR.z successfully read in.
LUT DAVIS - reflectivity data file: RoughTeflon_LUTR.z read in!
/run/geometryModified
/run/initialize
/analysis/setFileName davis_roughTeflon
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 18 starts.
... create file : davis_roughTeflon.root - done
... open analysis file : davis_roughTeflon.root - done
... open analysis file : davis_roughTeflon.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.030000s Real=0.037674s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.001
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 7
Total # of surface events: 9993
Unaccounted for: 0
Surface events by process:
Lobe reflection: 9567
Absorption: 388
Detection: 38
Sum: 9993
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.775471 rms = 0.208109
Y momentum dir of backward-going photons: mean = 0.209602 rms = 0.430271
Z momentum dir of backward-going photons: mean = 0.00265455 rms = 0.355334
X momentum dir of forward-going photons: mean = 0 rms = 0
Y momentum dir of forward-going photons: mean = 0 rms = 0
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : davis_roughTeflon.root - done
... close file : davis_roughTeflon.root - done
/opnovice2/surfaceFinish Detector_LUT
G4OpticalSurface: data file /cvmfs/geant4.cern.ch/share/data/RealSurface2.2/Detector_LUT.z successfully read in.
LUT DAVIS - data file: Detector_LUT.z read in!
G4OpticalSurface: data file /cvmfs/geant4.cern.ch/share/data/RealSurface2.2/Detector_LUTR.z successfully read in.
LUT DAVIS - reflectivity data file: Detector_LUTR.z read in!
/run/geometryModified
/analysis/setFileName davis_detector
/run/beamOn 10000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 19 starts.
... create file : davis_detector.root - done
... open analysis file : davis_detector.root - done
... open analysis file : davis_detector.root - done
Run terminated.
Run Summary
Number of events processed : 10000
User=0.020000s Real=0.020572s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 10000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.0006
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+04
OpAbsorption before surface: 6
Total # of surface events: 9994
Unaccounted for: 0
Surface events by process:
Lobe reflection: 118
Absorption: 8885
Detection: 991
Sum: 9994
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.90612 rms = 0.00201339
Y momentum dir of backward-going photons: mean = 0.422971 rms = 0.00433212
Z momentum dir of backward-going photons: mean = -0.000416827 rms = 0.00443173
X momentum dir of forward-going photons: mean = 0 rms = 0
Y momentum dir of forward-going photons: mean = 0 rms = 0
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : davis_detector.root - done
... close file : davis_detector.root - done
###############################################################################
################################# LUT (LBNL) ##################################
###############################################################################
/opnovice2/surfaceModel LUT
/run/geometryModified
/opnovice2/surfaceType dielectric_LUT
/run/geometryModified
/opnovice2/surfaceFinish polishedlumirrorair
G4OpticalSurface: data file /cvmfs/geant4.cern.ch/share/data/RealSurface2.2/PolishedLumirror.z successfully read in.
LUT - data file: PolishedLumirror.z read in!
/run/geometryModified
/analysis/setFileName lbnl_polishedLumirrorAir
/run/beamOn 1000
========= Table of registered couples ============================
Index : 0 used in the geometry : Yes
Material : G4_AIR
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
Index : 1 used in the geometry : Yes
Material : G4_WATER
Range cuts : gamma 700 um e- 700 um e+ 700 um proton 700 um
Energy thresholds : gamma 2.52521 keV e- 277.633 keV e+ 270.823 keV proton 70 keV
Region(s) which use this couple :
DefaultRegionForTheWorld
==================================================================
### Run 20 starts.
... create file : lbnl_polishedLumirrorAir.root - done
... open analysis file : lbnl_polishedLumirrorAir.root - done
... open analysis file : lbnl_polishedLumirrorAir.root - done
Run terminated.
Run Summary
Number of events processed : 1000
User=0.010000s Real=0.009039s Sys=0.000000s
Run Summary
---------------------------------
Primary particle was: opticalphoton with energy 3 eV .
Number of events: 1000
Material of world: G4_AIR
Material of tank: G4_WATER
Average number of photons absorbed by WLS per event: 0
Average number of photons created by WLS per event: 0
Average energy of WLS photons created per event: 0 eV.
Average number of photons absorbed by WLS2 per event: 0
Average number of photons created by WLS2 per event: 0
Average energy of WLS2 photons created per event: 0 eV.
Average number of OpRayleigh per event: 0
Average number of OpAbsorption per event: 0.004
Surface events (on +X surface, maximum one per photon) this run:
# of primary particles: 1e+03
OpAbsorption before surface: 1
Total # of surface events: 999
Unaccounted for: 0
Surface events by process:
Transmission: 103
Absorption: 104
Detection: 5
Polished Lumirror Air reflection: 787
Sum: 999
Unaccounted for: 0
---------------------------------
Histogram statistics for the entire run:
X momentum dir of backward-going photons: mean = -0.932034 rms = 0.0376591
Y momentum dir of backward-going photons: mean = 0.144423 rms = 0.255407
Z momentum dir of backward-going photons: mean = -0.00486819 rms = 0.209237
X momentum dir of forward-going photons: mean = 0.898384 rms = 3.65002e-08
Y momentum dir of forward-going photons: mean = 0.43921 rms = 9.12506e-09
Z momentum dir of forward-going photons: mean = 0 rms = 0
X momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Y momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
Z momentum dir of Fresnel-refracted photons: mean = 0 rms = 0
... write file : lbnl_polishedLumirrorAir.root - done
... close file : lbnl_polishedLumirrorAir.root - done
Graphics systems deleted.
Visualization Manager deleting...
================== Deleting memory pools ===================
Number of memory pools allocated: 10 of which, static: 0
Dynamic pools deleted: 10 / Total memory freed: 0.016 MB
============================================================