From 79bcbbfb17a12de5ddd10873e846723a74a559cc Mon Sep 17 00:00:00 2001
From: marcocle <marcocle@4525493e-7705-40b1-a816-d608a930855b>
Date: Wed, 9 Sep 2015 15:43:34 +0000
Subject: [PATCH] move_package: moved the Geant4/G4externalslcg package to the
obsolete project
git-svn-id: svn+ssh://svn.cern.ch/reps/lhcb/Geant4/trunk@194757 4525493e-7705-40b1-a816-d608a930855b
---
CMakeLists.txt | 2 +-
Geant4/G4config/CMakeLists.txt | 1 +
Geant4/G4config/doc/release.notes | 8 +
Geant4/G4config/patches/geant4.cmake.patch | 210 +-
Geant4/G4examples/cmt/requirements | 3 +-
Geant4/G4examples/doc/release.notes | 4 +-
.../G4TestEm3/Test/opt1NoApplyCuts/run.sh | 27 +-
.../G4TestEm3/cmt/requirements | 6 +-
.../G4TestEm3/scripts/Makefile | 3 -
.../G4TestEm3/scripts/OUTPUT_3.000 | Bin 166 -> 0 bytes
.../electromagnetic/G4TestEm3/scripts/Plot.C | 439 +++
.../G4TestEm3/scripts/TestEm3.in | 36 -
.../G4TestEm3/scripts/atlashec.mac | 32 -
.../G4TestEm3/scripts/dedx.mac | 31 -
.../G4TestEm3/scripts/emtutor.mac | 20 -
.../G4TestEm3/scripts/ionC12.mac | 36 -
.../G4TestEm3/scripts/lhcb.mac | 131 -
.../G4TestEm3/scripts/linac.mac | 31 -
.../G4TestEm3/scripts/lockwood.mac | 35 -
.../G4TestEm3/scripts/opt1noapplycuts.mac | 307 +++
.../G4TestEm3/scripts/plotHisto.C | 31 -
.../G4TestEm3/scripts/retrieveTables.mac | 25 -
.../electromagnetic/G4TestEm3/scripts/run.sh | 27 +
.../G4TestEm3/scripts/run01.mac | 30 -
.../G4TestEm3/scripts/run02.mac | 32 -
.../G4TestEm3/scripts/storeTables.mac | 28 -
.../G4TestEm3/scripts/tileCal.mac | 30 -
.../electromagnetic/G4TestEm3/scripts/vis.mac | 79 -
.../G4TestEm5/cmt/requirements | 4 +
.../G4TestEm5/scripts/Em5run.sh | 30 +
.../G4TestEm5/scripts/MacGen.py | 13 +
.../G4TestEm5/scripts/TestEm5.in | 63 -
.../G4TestEm5/scripts/TestEm5.out | 2368 -----------------
.../G4TestEm5/scripts/_CMakeLists.txt | 29 -
.../G4TestEm5/scripts/acosta.mac | 29 -
.../G4TestEm5/scripts/anthony.mac | 31 -
.../G4TestEm5/scripts/berger.mac | 36 -
.../G4TestEm5/scripts/bichsel.mac | 31 -
.../G4TestEm5/scripts/dedx1.mac | 29 -
.../G4TestEm5/scripts/dedx2.mac | 29 -
.../G4TestEm5/scripts/fluo.mac | 35 -
.../G4TestEm5/scripts/gammaSpectrum.mac | 31 -
.../G4TestEm5/scripts/gottsch.mac | 31 -
.../G4TestEm5/scripts/hanson.mac | 35 -
.../G4TestEm5/scripts/hunger.mac | 30 -
.../electromagnetic/G4TestEm5/scripts/ion.mac | 26 -
.../G4TestEm5/scripts/kulchi.mac | 31 -
.../G4TestEm5/scripts/mumsc.mac | 34 -
.../G4TestEm5/scripts/mutev.mac | 39 -
.../electromagnetic/G4TestEm5/scripts/pixe.C | 23 -
.../G4TestEm5/scripts/pixe.mac | 40 -
.../G4TestEm5/scripts/plotHisto.C | 14 -
.../G4TestEm5/scripts/plotHisto.kumac | 14 -
.../G4TestEm5/scripts/shen1.mac | 31 -
.../G4TestEm5/scripts/shen2.mac | 32 -
.../G4TestEm5/scripts/tavora.mac | 29 -
.../G4TestEm5/scripts/tramu.mac | 34 -
.../G4TestEm5/scripts/vincour.mac | 31 -
.../electromagnetic/G4TestEm5/scripts/vis.mac | 77 -
cmt/project.cmt | 2 +-
toolchain.cmake | 2 +-
61 files changed, 995 insertions(+), 3962 deletions(-)
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Makefile
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/OUTPUT_3.000
create mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Plot.C
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/TestEm3.in
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/atlashec.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/dedx.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/emtutor.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/ionC12.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lhcb.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/linac.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lockwood.mac
create mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/opt1noapplycuts.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/plotHisto.C
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/retrieveTables.mac
create mode 100755 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run.sh
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run01.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run02.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/storeTables.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/tileCal.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/vis.mac
create mode 100755 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/Em5run.sh
create mode 100755 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/MacGen.py
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.in
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.out
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/_CMakeLists.txt
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/acosta.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/anthony.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/berger.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/bichsel.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx1.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx2.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/fluo.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gammaSpectrum.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gottsch.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hanson.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hunger.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/ion.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/kulchi.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mumsc.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mutev.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.C
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.C
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.kumac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen1.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen2.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tavora.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tramu.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vincour.mac
delete mode 100644 Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vis.mac
diff --git a/CMakeLists.txt b/CMakeLists.txt
index 6eafb2e..c919055 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -13,7 +13,7 @@ if(DEFINED ENV{G4_NATIVE_VERSION})
endif()
find_package(GaudiProject)
-gaudi_project(Geant4 v96r4p2
+gaudi_project(Geant4 v96r4p2g1
DATA Geant4Files VERSION v96r*)
# FIXME: hack for compatibility with CMT build
diff --git a/Geant4/G4config/CMakeLists.txt b/Geant4/G4config/CMakeLists.txt
index eb08f3b..c70c6d3 100644
--- a/Geant4/G4config/CMakeLists.txt
+++ b/Geant4/G4config/CMakeLists.txt
@@ -27,6 +27,7 @@ ExternalProject_Add(Geant4
URL http://geant4.cern.ch/support/source/geant${Geant4_version}.tar.gz
PATCH_COMMAND ${PYTHON_EXECUTABLE} ${CMAKE_CURRENT_SOURCE_DIR}/patches/patch_g4_sources.py ${CMAKE_SOURCE_DIR} ${Geant4_version}
CMAKE_ARGS -DCMAKE_TOOLCHAIN_FILE=${CMAKE_SOURCE_DIR}/toolchain.cmake
+ -DGEANT4_BUILD_CXXSTD=${GAUDI_CXX_STANDARD}
-DGEANT4_INSTALL_DATA=OFF
-DGEANT4_USE_SYSTEM_CLHEP=ON
-DGEANT4_USE_XM=ON
diff --git a/Geant4/G4config/doc/release.notes b/Geant4/G4config/doc/release.notes
index 601bb37..325affe 100755
--- a/Geant4/G4config/doc/release.notes
+++ b/Geant4/G4config/doc/release.notes
@@ -4,6 +4,14 @@
! Purpose : Configuration package for Geant4 build
!-----------------------------------------------------------------------------
+! 2015-09-09 - Marco Clemencic
+ - Fixed ccache (and custom launchers) support in CMake build.
+ - Use the same C++ standard as Gaudi in CMake builds.
+ - Build and install a few examples
+ - enable G4analysis
+ - fixed Geant4Config.cmake
+ - enabled the builds of selected examples
+
! ======================= G4config v96r4p2 2015-07-24 ========================
! 2015-07-24 - Gloria Corti
- Only change local version to match the fact that we use Geant4 9.6.p04
diff --git a/Geant4/G4config/patches/geant4.cmake.patch b/Geant4/G4config/patches/geant4.cmake.patch
index 57b8985..896a6b3 100644
--- a/Geant4/G4config/patches/geant4.cmake.patch
+++ b/Geant4/G4config/patches/geant4.cmake.patch
@@ -1,6 +1,6 @@
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4BuildProjectConfig.cmake b/geant4.9.6.p03/cmake/Modules/Geant4BuildProjectConfig.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4BuildProjectConfig.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4BuildProjectConfig.cmake 2014-11-18 13:48:32.512838961 +0100
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4BuildProjectConfig.cmake b/geant4.9.6.p04/cmake/Modules/Geant4BuildProjectConfig.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4BuildProjectConfig.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4BuildProjectConfig.cmake 2015-09-09 17:03:22.291420708 +0200
@@ -173,13 +173,13 @@
# Use ABSOLUTE paths...
set(GEANT4_INCLUDE_DIR_SETUP "
@@ -17,9 +17,9 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4BuildProjectConfig.cmake b/geant4
)
set(GEANT4_INCLUDE_DIR_SETUP "
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4ConfigureConfigScript.cmake b/geant4.9.6.p03/cmake/Modules/Geant4ConfigureConfigScript.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4ConfigureConfigScript.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4ConfigureConfigScript.cmake 2014-11-18 13:48:32.512838961 +0100
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4ConfigureConfigScript.cmake b/geant4.9.6.p04/cmake/Modules/Geant4ConfigureConfigScript.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4ConfigureConfigScript.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4ConfigureConfigScript.cmake 2015-09-09 17:03:22.291420708 +0200
@@ -230,7 +230,7 @@
set(GEANT4_CONFIG_INSTALL_PREFIX "${CMAKE_INSTALL_PREFIX}")
set(GEANT4_CONFIG_INSTALL_EXECPREFIX \"\")
@@ -38,9 +38,9 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4ConfigureConfigScript.cmake b/gea
endif()
# Configure the install tree script
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4MacroLibraryTargets.cmake b/geant4.9.6.p03/cmake/Modules/Geant4MacroLibraryTargets.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4MacroLibraryTargets.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4MacroLibraryTargets.cmake 2014-11-18 13:48:32.512838961 +0100
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4MacroLibraryTargets.cmake b/geant4.9.6.p04/cmake/Modules/Geant4MacroLibraryTargets.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4MacroLibraryTargets.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4MacroLibraryTargets.cmake 2015-09-09 17:03:22.292420731 +0200
@@ -166,7 +166,7 @@
# Header install?
@@ -68,10 +68,10 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4MacroLibraryTargets.cmake b/geant
COMPONENT Development)
# Store the include path of the component so that the build tree
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4MakeRules_cxx.cmake b/geant4.9.6.p03/cmake/Modules/Geant4MakeRules_cxx.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4MakeRules_cxx.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4MakeRules_cxx.cmake 2014-11-18 13:50:59.456948558 +0100
-@@ -240,3 +240,64 @@
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4MakeRules_cxx.cmake b/geant4.9.6.p04/cmake/Modules/Geant4MakeRules_cxx.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4MakeRules_cxx.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4MakeRules_cxx.cmake 2015-09-09 17:03:22.297420846 +0200
+@@ -240,3 +240,73 @@
# - CMake may do a reasonable job on its own here...
endif()
@@ -136,9 +136,18 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4MakeRules_cxx.cmake b/geant4.9.6.
+ set(GAUDI_RULE_LAUNCH_LINK "env ${GAUDI_RULE_LAUNCH_LINK}")
+ endif()
+endif()
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4OptionalComponents.cmake b/geant4.9.6.p03/cmake/Modules/Geant4OptionalComponents.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4OptionalComponents.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4OptionalComponents.cmake 2014-11-18 13:55:24.512555325 +0100
++
++# apply launch rules
++foreach(_rule COMPILE LINK CUSTOM)
++ if(GAUDI_RULE_LAUNCH_${_rule})
++ string(STRIP "${GAUDI_RULE_LAUNCH_${_rule}}" GAUDI_RULE_LAUNCH_${_rule})
++ set_property(GLOBAL PROPERTY RULE_LAUNCH_${_rule} "${GAUDI_RULE_LAUNCH_${_rule}}")
++ message(STATUS "Prefix ${_rule} commands with '${GAUDI_RULE_LAUNCH_${_rule}}'")
++ endif()
++endforeach()
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4OptionalComponents.cmake b/geant4.9.6.p04/cmake/Modules/Geant4OptionalComponents.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4OptionalComponents.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4OptionalComponents.cmake 2015-09-09 17:03:22.297420846 +0200
@@ -34,7 +34,7 @@
if(GEANT4_USE_SYSTEM_CLHEP)
# We keep this as required, because if the user chooses to use a
@@ -148,9 +157,9 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4OptionalComponents.cmake b/geant4
set(GEANT4_USE_SYSTEM_CLHEP TRUE)
else()
set(CLHEP_FOUND TRUE)
-diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake b/geant4.9.6.p03/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake
---- a/geant4.9.6.p03/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake 2014-03-17 12:47:01.000000000 +0100
-+++ b/geant4.9.6.p03/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake 2014-11-18 13:48:32.512838961 +0100
+diff -urN a/geant4.9.6.p04/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake b/geant4.9.6.p04/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake
+--- a/geant4.9.6.p04/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake 2015-02-04 09:53:22.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Modules/Geant4ToolchainBackwardCompatibility.cmake 2015-09-09 17:03:22.297420846 +0200
@@ -531,7 +531,7 @@
# If we are relocatable, then the structure we will have is
# +- CMAKE_INSTALL_PREFIX
@@ -169,10 +178,40 @@ diff -urN a/geant4.9.6.p03/cmake/Modules/Geant4ToolchainBackwardCompatibility.cm
)
set(G4INCLUDE "\"`cd \$geant4make_root/${G4MAKE_TO_INCLUDEDIR} > /dev/null \; pwd`\"")
-diff -urN a/geant4.9.6.p03/CMakeLists.txt b/geant4.9.6.p03/CMakeLists.txt
---- a/geant4.9.6.p03/CMakeLists.txt 2014-03-17 12:47:39.000000000 +0100
-+++ b/geant4.9.6.p03/CMakeLists.txt 2014-11-18 13:48:32.512838961 +0100
-@@ -113,8 +113,8 @@
+diff -urN a/geant4.9.6.p04/cmake/Templates/Geant4Config.cmake.in b/geant4.9.6.p04/cmake/Templates/Geant4Config.cmake.in
+--- a/geant4.9.6.p04/cmake/Templates/Geant4Config.cmake.in 2015-02-04 09:53:23.000000000 +0100
++++ b/geant4.9.6.p04/cmake/Templates/Geant4Config.cmake.in 2015-09-09 16:11:28.000000000 +0200
+@@ -382,8 +382,6 @@
+ set(_geant4_internal_libraries
+ G4Tree${_geant4_lib_use_suffix}
+ G4FR${_geant4_lib_use_suffix}
+- G4GMocren${_geant4_lib_use_suffix}
+- G4visHepRep${_geant4_lib_use_suffix}
+ G4RayTracer${_geant4_lib_use_suffix}
+ G4VRML${_geant4_lib_use_suffix}
+ )
+@@ -413,13 +411,14 @@
+ list(APPEND _geant4_internal_libraries
+ G4vis_management${_geant4_lib_use_suffix}
+ G4modeling${_geant4_lib_use_suffix}
+- G4interfaces${_geant4_lib_use_suffix}
++ G4UIGAG${_geant4_lib_use_suffix}
++ G4UIbasic${_geant4_lib_use_suffix}
++ G4UIcommon${_geant4_lib_use_suffix}
+ G4persistency${_geant4_lib_use_suffix}
+ ${_geant4_g3tog4_library}
+ G4analysis${_geant4_lib_use_suffix}
+- G4error_propagation${_geant4_lib_use_suffix}
+ G4readout${_geant4_lib_use_suffix}
+- G4physicslists${_geant4_lib_use_suffix}
++ G4physics_lists${_geant4_lib_use_suffix}
+ G4run${_geant4_lib_use_suffix}
+ G4event${_geant4_lib_use_suffix}
+ G4tracking${_geant4_lib_use_suffix}
+diff -urN a/geant4.9.6.p04/CMakeLists.txt b/geant4.9.6.p04/CMakeLists.txt
+--- a/geant4.9.6.p04/CMakeLists.txt 2015-02-04 10:03:00.000000000 +0100
++++ b/geant4.9.6.p04/CMakeLists.txt 2015-09-09 17:03:22.298420869 +0200
+@@ -136,8 +136,8 @@
# Generate any Use/Config files here once everything else has been processed
# e.g. "UseGeant4.cmake", "Geant4Config.cmake" - library dependencies etc
#
@@ -183,14 +222,35 @@ diff -urN a/geant4.9.6.p03/CMakeLists.txt b/geant4.9.6.p03/CMakeLists.txt
#----------------------------------------------------------------------------
# Create the Geant4Config files and supporting modules
-diff -urN a/geant4.9.6.p03/source/CMakeLists.txt b/geant4.9.6.p03/source/CMakeLists.txt
---- a/geant4.9.6.p03/source/CMakeLists.txt 2014-03-17 12:46:55.000000000 +0100
-+++ b/geant4.9.6.p03/source/CMakeLists.txt 2014-11-18 13:48:32.509838896 +0100
-@@ -11,11 +11,8 @@
- #
- #------------------------------------------------------------------------------
+@@ -154,10 +154,18 @@
+ add_subdirectory(benchmarks)
+ endif()
+ endif()
+-if(GEANT4_BUILD_EXAMPLES)
+- set(Geant4_DIR ${CMAKE_BINARY_DIR} CACHE PATH "Current build directory")
+- add_subdirectory(examples)
+-endif()
++
++# build examples required by LHCb (see LHCBGAUSS-505)
++set(Geant4_DIR ${CMAKE_BINARY_DIR} CACHE PATH "Current build directory")
++option(WITH_GEANT4_UIVIS "Build example with Geant4 UI and Vis drivers" OFF)
++set(G4UI_NONE ON CACHE BOOL "disable UI for examples")
++set(G4VIS_NONE ON CACHE BOOL "disable Vis for examples")
++add_subdirectory(examples/extended/electromagnetic/TestEm3)
++add_subdirectory(examples/extended/electromagnetic/TestEm5)
++add_subdirectory(examples/extended/electromagnetic/TestEm7)
++add_subdirectory(examples/extended/electromagnetic/TestEm9)
++add_subdirectory(examples/extended/hadronic/Hadr00)
++
+ #-----------------------------------------------------------------------------
+ # Provide packaging with CPack.
+ include(Geant4CPackBase)
+diff -urN a/geant4.9.6.p04/source/CMakeLists.txt b/geant4.9.6.p04/source/CMakeLists.txt
+--- a/geant4.9.6.p04/source/CMakeLists.txt 2015-02-04 09:57:54.000000000 +0100
++++ b/geant4.9.6.p04/source/CMakeLists.txt 2015-09-09 17:03:22.298420869 +0200
+@@ -13,9 +13,7 @@
--add_subdirectory(analysis)
+ add_subdirectory(analysis)
add_subdirectory(digits_hits)
-add_subdirectory(error_propagation)
add_subdirectory(event)
@@ -198,9 +258,9 @@ diff -urN a/geant4.9.6.p03/source/CMakeLists.txt b/geant4.9.6.p03/source/CMakeLi
# Target G4g3tog4 is optional based on user selection...
if(GEANT4_USE_G3TOG4)
-diff -urN a/geant4.9.6.p03/source/interfaces/basic/sources.cmake b/geant4.9.6.p03/source/interfaces/basic/sources.cmake
---- a/geant4.9.6.p03/source/interfaces/basic/sources.cmake 2014-03-17 12:46:13.000000000 +0100
-+++ b/geant4.9.6.p03/source/interfaces/basic/sources.cmake 2014-11-18 13:48:32.509838896 +0100
+diff -urN a/geant4.9.6.p04/source/interfaces/basic/sources.cmake b/geant4.9.6.p04/source/interfaces/basic/sources.cmake
+--- a/geant4.9.6.p04/source/interfaces/basic/sources.cmake 2015-02-04 09:53:34.000000000 +0100
++++ b/geant4.9.6.p04/source/interfaces/basic/sources.cmake 2015-09-09 17:03:22.298420869 +0200
@@ -145,6 +145,7 @@
G4globman
G4intercoms
@@ -209,9 +269,9 @@ diff -urN a/geant4.9.6.p03/source/interfaces/basic/sources.cmake b/geant4.9.6.p0
G4global
G4intercoms
LINK_LIBRARIES
-diff -urN a/geant4.9.6.p03/source/interfaces/CMakeLists.txt b/geant4.9.6.p03/source/interfaces/CMakeLists.txt
---- a/geant4.9.6.p03/source/interfaces/CMakeLists.txt 2014-03-17 12:46:14.000000000 +0100
-+++ b/geant4.9.6.p03/source/interfaces/CMakeLists.txt 2014-11-18 13:48:32.510838918 +0100
+diff -urN a/geant4.9.6.p04/source/interfaces/CMakeLists.txt b/geant4.9.6.p04/source/interfaces/CMakeLists.txt
+--- a/geant4.9.6.p04/source/interfaces/CMakeLists.txt 2015-02-04 09:53:37.000000000 +0100
++++ b/geant4.9.6.p04/source/interfaces/CMakeLists.txt 2015-09-09 17:03:22.299420892 +0200
@@ -18,11 +18,27 @@
if(NOT GEANT4_BUILD_GRANULAR_LIBS)
@@ -244,9 +304,9 @@ diff -urN a/geant4.9.6.p03/source/interfaces/CMakeLists.txt b/geant4.9.6.p03/sou
+
endif()
-diff -urN a/geant4.9.6.p03/source/interfaces/GAG/sources.cmake b/geant4.9.6.p03/source/interfaces/GAG/sources.cmake
---- a/geant4.9.6.p03/source/interfaces/GAG/sources.cmake 2014-03-17 12:46:14.000000000 +0100
-+++ b/geant4.9.6.p03/source/interfaces/GAG/sources.cmake 2014-11-18 13:48:32.510838918 +0100
+diff -urN a/geant4.9.6.p04/source/interfaces/GAG/sources.cmake b/geant4.9.6.p04/source/interfaces/GAG/sources.cmake
+--- a/geant4.9.6.p04/source/interfaces/GAG/sources.cmake 2015-02-04 09:53:37.000000000 +0100
++++ b/geant4.9.6.p04/source/interfaces/GAG/sources.cmake 2015-09-09 17:03:22.299420892 +0200
@@ -39,6 +39,7 @@
G4globman
G4intercoms
@@ -255,9 +315,9 @@ diff -urN a/geant4.9.6.p03/source/interfaces/GAG/sources.cmake b/geant4.9.6.p03/
G4global
G4intercoms
LINK_LIBRARIES
-diff -urN a/geant4.9.6.p03/source/persistency/CMakeLists.txt b/geant4.9.6.p03/source/persistency/CMakeLists.txt
---- a/geant4.9.6.p03/source/persistency/CMakeLists.txt 2014-03-17 12:46:15.000000000 +0100
-+++ b/geant4.9.6.p03/source/persistency/CMakeLists.txt 2014-11-18 13:48:32.510838918 +0100
+diff -urN a/geant4.9.6.p04/source/persistency/CMakeLists.txt b/geant4.9.6.p04/source/persistency/CMakeLists.txt
+--- a/geant4.9.6.p04/source/persistency/CMakeLists.txt 2015-02-04 09:57:44.000000000 +0100
++++ b/geant4.9.6.p04/source/persistency/CMakeLists.txt 2015-09-09 17:03:22.299420892 +0200
@@ -12,30 +12,17 @@
#
#------------------------------------------------------------------------------
@@ -289,9 +349,9 @@ diff -urN a/geant4.9.6.p03/source/persistency/CMakeLists.txt b/geant4.9.6.p03/so
)
endif()
endif()
-diff -urN a/geant4.9.6.p03/source/physics_lists/builders/sources.cmake b/geant4.9.6.p03/source/physics_lists/builders/sources.cmake
---- a/geant4.9.6.p03/source/physics_lists/builders/sources.cmake 2014-03-17 12:46:22.000000000 +0100
-+++ b/geant4.9.6.p03/source/physics_lists/builders/sources.cmake 2014-11-18 13:54:07.010915908 +0100
+diff -urN a/geant4.9.6.p04/source/physics_lists/builders/sources.cmake b/geant4.9.6.p04/source/physics_lists/builders/sources.cmake
+--- a/geant4.9.6.p04/source/physics_lists/builders/sources.cmake 2015-02-04 09:54:21.000000000 +0100
++++ b/geant4.9.6.p04/source/physics_lists/builders/sources.cmake 2015-09-09 17:03:22.299420892 +0200
@@ -40,12 +40,6 @@
include_directories(${CMAKE_SOURCE_DIR}/source/particles/shortlived/include)
include_directories(${CMAKE_SOURCE_DIR}/source/processes/cuts/include)
@@ -333,9 +393,9 @@ diff -urN a/geant4.9.6.p03/source/physics_lists/builders/sources.cmake b/geant4.
G4emhighenergy
G4emlowenergy
G4emstandard
-diff -urN a/geant4.9.6.p03/source/physics_lists/CMakeLists.txt b/geant4.9.6.p03/source/physics_lists/CMakeLists.txt
---- a/geant4.9.6.p03/source/physics_lists/CMakeLists.txt 2014-03-17 12:46:21.000000000 +0100
-+++ b/geant4.9.6.p03/source/physics_lists/CMakeLists.txt 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/physics_lists/CMakeLists.txt b/geant4.9.6.p04/source/physics_lists/CMakeLists.txt
+--- a/geant4.9.6.p04/source/physics_lists/CMakeLists.txt 2015-02-04 09:54:23.000000000 +0100
++++ b/geant4.9.6.p04/source/physics_lists/CMakeLists.txt 2015-09-09 17:03:22.300420915 +0200
@@ -17,10 +17,11 @@
if(NOT GEANT4_BUILD_GRANULAR_LIBS)
@@ -349,9 +409,9 @@ diff -urN a/geant4.9.6.p03/source/physics_lists/CMakeLists.txt b/geant4.9.6.p03/
endif()
+add_subdirectory(LHCb)
-diff -urN a/geant4.9.6.p03/source/physics_lists/LHCb/CMakeLists.txt b/geant4.9.6.p03/source/physics_lists/LHCb/CMakeLists.txt
---- a/geant4.9.6.p03/source/physics_lists/LHCb/CMakeLists.txt 1970-01-01 01:00:00.000000000 +0100
-+++ b/geant4.9.6.p03/source/physics_lists/LHCb/CMakeLists.txt 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/physics_lists/LHCb/CMakeLists.txt b/geant4.9.6.p04/source/physics_lists/LHCb/CMakeLists.txt
+--- a/geant4.9.6.p04/source/physics_lists/LHCb/CMakeLists.txt 1970-01-01 01:00:00.000000000 +0100
++++ b/geant4.9.6.p04/source/physics_lists/LHCb/CMakeLists.txt 2015-09-09 17:03:22.300420915 +0200
@@ -0,0 +1,20 @@
+#------------------------------------------------------------------------------
+# CMakeLists.txt
@@ -373,9 +433,9 @@ diff -urN a/geant4.9.6.p03/source/physics_lists/LHCb/CMakeLists.txt b/geant4.9.6
+ GEANT4_GLOBAL_LIBRARY_TARGET(COMPONENTS sources.cmake)
+endif()
+
-diff -urN a/geant4.9.6.p03/source/physics_lists/LHCb/sources.cmake b/geant4.9.6.p03/source/physics_lists/LHCb/sources.cmake
---- a/geant4.9.6.p03/source/physics_lists/LHCb/sources.cmake 1970-01-01 01:00:00.000000000 +0100
-+++ b/geant4.9.6.p03/source/physics_lists/LHCb/sources.cmake 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/physics_lists/LHCb/sources.cmake b/geant4.9.6.p04/source/physics_lists/LHCb/sources.cmake
+--- a/geant4.9.6.p04/source/physics_lists/LHCb/sources.cmake 1970-01-01 01:00:00.000000000 +0100
++++ b/geant4.9.6.p04/source/physics_lists/LHCb/sources.cmake 2015-09-09 17:03:22.300420915 +0200
@@ -0,0 +1,210 @@
+#------------------------------------------------------------------------------
+# sources.cmake
@@ -587,9 +647,9 @@ diff -urN a/geant4.9.6.p03/source/physics_lists/LHCb/sources.cmake b/geant4.9.6.
+
+# List any source specific properties here
+
-diff -urN a/geant4.9.6.p03/source/physics_lists/lists/sources.cmake b/geant4.9.6.p03/source/physics_lists/lists/sources.cmake
---- a/geant4.9.6.p03/source/physics_lists/lists/sources.cmake 2014-03-17 12:46:21.000000000 +0100
-+++ b/geant4.9.6.p03/source/physics_lists/lists/sources.cmake 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/physics_lists/lists/sources.cmake b/geant4.9.6.p04/source/physics_lists/lists/sources.cmake
+--- a/geant4.9.6.p04/source/physics_lists/lists/sources.cmake 2015-02-04 09:54:23.000000000 +0100
++++ b/geant4.9.6.p04/source/physics_lists/lists/sources.cmake 2015-09-09 17:03:22.301420938 +0200
@@ -43,7 +43,6 @@
include_directories(${CMAKE_SOURCE_DIR}/source/physics_lists/builders/include)
include_directories(${CMAKE_SOURCE_DIR}/source/processes/cuts/include)
@@ -606,9 +666,9 @@ diff -urN a/geant4.9.6.p03/source/physics_lists/lists/sources.cmake b/geant4.9.6
G4emhighenergy
G4emlowenergy
G4emstandard
-diff -urN a/geant4.9.6.p03/source/processes/CMakeLists.txt b/geant4.9.6.p03/source/processes/CMakeLists.txt
---- a/geant4.9.6.p03/source/processes/CMakeLists.txt 2014-03-17 12:46:27.000000000 +0100
-+++ b/geant4.9.6.p03/source/processes/CMakeLists.txt 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/processes/CMakeLists.txt b/geant4.9.6.p04/source/processes/CMakeLists.txt
+--- a/geant4.9.6.p04/source/processes/CMakeLists.txt 2015-02-04 09:57:35.000000000 +0100
++++ b/geant4.9.6.p04/source/processes/CMakeLists.txt 2015-09-09 17:03:22.301420938 +0200
@@ -37,12 +37,6 @@
cuts/sources.cmake
decay/sources.cmake
@@ -622,9 +682,9 @@ diff -urN a/geant4.9.6.p03/source/processes/CMakeLists.txt b/geant4.9.6.p03/sour
electromagnetic/highenergy/sources.cmake
electromagnetic/lowenergy/sources.cmake
electromagnetic/muons/sources.cmake
-diff -urN a/geant4.9.6.p03/source/processes/electromagnetic/CMakeLists.txt b/geant4.9.6.p03/source/processes/electromagnetic/CMakeLists.txt
---- a/geant4.9.6.p03/source/processes/electromagnetic/CMakeLists.txt 2014-03-17 12:46:53.000000000 +0100
-+++ b/geant4.9.6.p03/source/processes/electromagnetic/CMakeLists.txt 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/processes/electromagnetic/CMakeLists.txt b/geant4.9.6.p04/source/processes/electromagnetic/CMakeLists.txt
+--- a/geant4.9.6.p04/source/processes/electromagnetic/CMakeLists.txt 2015-02-04 09:55:49.000000000 +0100
++++ b/geant4.9.6.p04/source/processes/electromagnetic/CMakeLists.txt 2015-09-09 17:03:22.301420938 +0200
@@ -12,7 +12,6 @@
#------------------------------------------------------------------------------
@@ -633,9 +693,9 @@ diff -urN a/geant4.9.6.p03/source/processes/electromagnetic/CMakeLists.txt b/gea
add_subdirectory(highenergy)
add_subdirectory(lowenergy)
add_subdirectory(muons)
-diff -urN a/geant4.9.6.p03/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc b/geant4.9.6.p03/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc
---- a/geant4.9.6.p03/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc 2014-03-17 12:46:37.000000000 +0100
-+++ b/geant4.9.6.p03/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc 2014-11-18 14:09:28.495401650 +0100
+diff -urN a/geant4.9.6.p04/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc b/geant4.9.6.p04/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc
+--- a/geant4.9.6.p04/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc 2015-02-04 09:57:06.000000000 +0100
++++ b/geant4.9.6.p04/source/processes/hadronic/models/im_r_matrix/src/G4CollisionInitialState.cc 2015-09-09 17:03:22.301420938 +0200
@@ -100,7 +100,7 @@
return theFSGenerator->GetFinalState(thePrimary, theTs);
}
@@ -645,9 +705,9 @@ diff -urN a/geant4.9.6.p03/source/processes/hadronic/models/im_r_matrix/src/G4Co
void G4CollisionInitialState::Print() const
{
-diff -urN a/geant4.9.6.p03/source/visualization/CMakeLists.txt b/geant4.9.6.p03/source/visualization/CMakeLists.txt
---- a/geant4.9.6.p03/source/visualization/CMakeLists.txt 2014-03-17 12:47:00.000000000 +0100
-+++ b/geant4.9.6.p03/source/visualization/CMakeLists.txt 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/visualization/CMakeLists.txt b/geant4.9.6.p04/source/visualization/CMakeLists.txt
+--- a/geant4.9.6.p04/source/visualization/CMakeLists.txt 2015-02-04 09:54:58.000000000 +0100
++++ b/geant4.9.6.p04/source/visualization/CMakeLists.txt 2015-09-09 17:03:22.302420961 +0200
@@ -12,13 +12,10 @@
#------------------------------------------------------------------------------
@@ -662,9 +722,9 @@ diff -urN a/geant4.9.6.p03/source/visualization/CMakeLists.txt b/geant4.9.6.p03/
add_subdirectory(management)
add_subdirectory(modeling)
-diff -urN a/geant4.9.6.p03/source/visualization/OpenGL/sources.cmake b/geant4.9.6.p03/source/visualization/OpenGL/sources.cmake
---- a/geant4.9.6.p03/source/visualization/OpenGL/sources.cmake 2014-03-17 12:46:56.000000000 +0100
-+++ b/geant4.9.6.p03/source/visualization/OpenGL/sources.cmake 2014-11-18 13:48:32.510838919 +0100
+diff -urN a/geant4.9.6.p04/source/visualization/OpenGL/sources.cmake b/geant4.9.6.p04/source/visualization/OpenGL/sources.cmake
+--- a/geant4.9.6.p04/source/visualization/OpenGL/sources.cmake 2015-02-04 09:54:39.000000000 +0100
++++ b/geant4.9.6.p04/source/visualization/OpenGL/sources.cmake 2015-09-09 17:03:22.302420961 +0200
@@ -312,7 +312,8 @@
G4global
G4graphics_reps
@@ -675,9 +735,9 @@ diff -urN a/geant4.9.6.p03/source/visualization/OpenGL/sources.cmake b/geant4.9.
G4modeling
G4tracking
G4vis_management
-diff -urN a/geant4.9.6.p03/source/visualization/OpenInventor/sources.cmake b/geant4.9.6.p03/source/visualization/OpenInventor/sources.cmake
---- a/geant4.9.6.p03/source/visualization/OpenInventor/sources.cmake 2014-03-17 12:47:00.000000000 +0100
-+++ b/geant4.9.6.p03/source/visualization/OpenInventor/sources.cmake 2014-11-18 13:48:32.512838961 +0100
+diff -urN a/geant4.9.6.p04/source/visualization/OpenInventor/sources.cmake b/geant4.9.6.p04/source/visualization/OpenInventor/sources.cmake
+--- a/geant4.9.6.p04/source/visualization/OpenInventor/sources.cmake 2015-02-04 09:54:47.000000000 +0100
++++ b/geant4.9.6.p04/source/visualization/OpenInventor/sources.cmake 2015-09-09 17:03:22.302420961 +0200
@@ -195,7 +195,7 @@
G4global
G4graphics_reps
diff --git a/Geant4/G4examples/cmt/requirements b/Geant4/G4examples/cmt/requirements
index f48513e..3b3e742 100644
--- a/Geant4/G4examples/cmt/requirements
+++ b/Geant4/G4examples/cmt/requirements
@@ -16,7 +16,8 @@ branches cmt doc extended
use G4Hadr00 v* Geant4/G4examples/extended/hadronic
use G4TestEm7 v* Geant4/G4examples/extended/electromagnetic
use G4TestEm9 v* Geant4/G4examples/extended/electromagnetic
-
+use G4TestEm3 v* Geant4/G4examples/extended/electromagnetic
+use G4TestEm5 v* Geant4/G4examples/extended/electromagnetic
include_path none
#============================================================================
diff --git a/Geant4/G4examples/doc/release.notes b/Geant4/G4examples/doc/release.notes
index e58fb91..0200630 100644
--- a/Geant4/G4examples/doc/release.notes
+++ b/Geant4/G4examples/doc/release.notes
@@ -3,7 +3,9 @@
! Responsible : James McCarthy, Luca Pescatore
! Purpose : Examples packages provided my Geant4, used as standalone tests
!-----------------------------------------------------------------------------
-
+!2015-07-29- Tim Williams
+ - Installed scripts to run TestEm3 and TestEm5, can be run with 'Calorimeter_test' and 'MultipleScattering_test' respectively
+ - Removed unused and unnecessary scripts from TestEm3 and TestEm5
!===================== Geant4/G4examples v6r1 2015-07-23 =====================
! 2015-07-23- Tim Williams
- Updated scripts to run testEm3 and testEm5 to be compiler independent
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/Test/opt1NoApplyCuts/run.sh b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/Test/opt1NoApplyCuts/run.sh
index c7159e2..4722bd0 100755
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/Test/opt1NoApplyCuts/run.sh
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/Test/opt1NoApplyCuts/run.sh
@@ -1,18 +1,27 @@
-#run with . run.sh to keep environment for SetupProject
+#!/bin/bash
-if [ "$1" = "sim" ]
-then
- Binary=$(find ../.. -name "testEm3.exe")
- echo $Binary
- $Binary opt1noapplycuts.mac
+CurrentDir="${PWD##*/}"
+echo $CurrentDir
+TestDir=CaloTest
+if [[ ! -d "CaloTest" && $CurrentDir != $TestDir ]]; then
+mkdir CaloTest
+cd CaloTest
fi
+
+if [[ -d "CaloTest" && $CurrentDir != $TestDir ]]; then
+cd CaloTest
+fi
+
+testEm3.exe $G4TESTEM3ROOT/scripts/opt1noapplycuts.mac
+
#There is a bug with ROOT 6 opening Geant4 created histograms, have to use old root for now
-SetupProject ROOT --version 72a
+source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.30/$CMTCONFIG/root/bin/thisroot.sh
-g++ -c `root-config --cflags` Plot.C
+g++ -c `root-config --cflags` $G4TESTEM3ROOT/scripts/Plot.C
g++ -o Plot `root-config --glibs` Plot.o
./Plot
-#creates full file of plots Save.root and Key results in Selectedresults.root
+echo "Results Produced in Save.root, Selectedresults.root and selectedresults.txt"
+
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/cmt/requirements b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/cmt/requirements
index 7ae8d36..389269b 100755
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/cmt/requirements
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/cmt/requirements
@@ -13,7 +13,7 @@ use G4LHCblists v* Geant4
include_dirs ../include
-#------ Example TestEm7- -------
+#------ Example TestEm3- -------
#------ Building the library -------
@@ -26,6 +26,10 @@ application testEm3 ../TestEm3.cc
#========================================================================
+#-------Install scripts -------
+apply_pattern install_scripts
+alias 'Calorimeter_test' 'run.sh'
+
private
apply_pattern G4_copy_examples_source
macro_append testEm3_dependencies " G4TestEm3 "
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Makefile b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Makefile
deleted file mode 100644
index ca71bc1..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Makefile
+++ /dev/null
@@ -1,3 +0,0 @@
-all:
- g++ -c -pthread -std=c++11 -Wno-deprecated-declarations -m64 -I/afs/cern.ch/sw/lcg/releases/LCG_70/ROOT/5.34.22/x86_64-slc6-gcc48-opt/include Plot.C
- g++ -o Plot -L/afs/cern.ch/sw/lcg/releases/LCG_70/ROOT/5.34.22/x86_64-slc6-gcc48-opt/lib -lGui -lCore -lCint -lRIO -lNet -lHist -lGraf -lGraf3d -lGpad -lTree -lRint -lPostscript -lMatrix -lPhysics -lMathCore -lThread -pthread -lm -ldl -rdynamic Plot.o
\ No newline at end of file
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/OUTPUT_3.000 b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/OUTPUT_3.000
deleted file mode 100644
index b9f997dd90a6e5c19fa6984ab39f4e972bf32181..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001
literal 166
zcmXTQ&o5zMnB%~}zyQL_fEXFrumE|C$Re2RGzKQ1paob$53JY$#AXd~%gjmT@ed6N
i2n~rh)-y0LU_j+DzYLfx{vL$EdKnmOptJ*&hEo7@DiVYM
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Plot.C b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Plot.C
new file mode 100644
index 0000000..768bace
--- /dev/null
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/Plot.C
@@ -0,0 +1,439 @@
+#include "TH1.h"
+#include<vector>
+#include "TCanvas.h"
+#include "TFile.h"
+#include "TStyle.h"
+#include "TF1.h"
+#include "TGraphErrors.h"
+#include "TMinuit.h"
+#include"TVirtualFitter.h"
+#include"TMatrixD.h"
+#include"TFitResult.h"
+#include<fstream>
+#include<iomanip>
+#include"TLegend.h"
+#include"TText.h"
+#include"TLatex.h"
+#include"TPaveStats.h"
+using namespace std;
+
+struct M{
+ int Name;
+ double Energy;
+ int rebin;
+};
+
+struct SeperateProfiles
+{
+ TH1D* PhotonsLead;
+ TH1D* ElectronsLead;
+ TH1D* PositronsLead;
+ TH1D* PhotonsScintillator;
+ TH1D* ElectronsScintillator;
+ TH1D* PositronsScintillator;
+};
+
+vector<SeperateProfiles> Profiles(15);
+vector<M> Control(15);
+double nEvents=10000;
+int main()
+{
+
+
+ TCanvas *c[15];
+ ofstream outt;
+ outt.open ("results.txt");
+ ofstream selected;
+ selected.open("selectedresults.txt");
+ ofstream ErrorInv("Errors.txt");
+
+ for(int i=0; i<15;i++)
+ {
+ Control[i].Name=i+1;
+
+ }
+ Control[0].Energy=100.0;
+ Control[1].Energy=44.44;
+ Control[2].Energy=25.00;
+ Control[3].Energy=16.00;
+ Control[4].Energy=11.11;
+ Control[5].Energy=8.16;
+ Control[6].Energy=6.25;
+ Control[7].Energy=4.94;
+ Control[8].Energy=4.00;
+ Control[9].Energy=3.31;
+ Control[10].Energy=2.78;
+ Control[11].Energy=2.37;
+ Control[12].Energy=2.04;
+ Control[13].Energy=1.78;
+ Control[14].Energy=1.56;
+
+ Control[14].rebin=6;
+ Control[13].rebin=5;
+ Control[12].rebin=5;
+ Control[11].rebin=5;
+ Control[10].rebin=5;
+ Control[9].rebin=4;
+ Control[8].rebin=4;
+ Control[7].rebin=3;
+ Control[6].rebin=2;
+ Control[5].rebin=2;
+ Control[4].rebin=2;
+ Control[3].rebin=2;
+ Control[2].rebin=2;
+ Control[1].rebin=1;
+ Control[0].rebin=1;
+
+
+
+ TH1D* h1[15];
+ TH1D* h2[15];
+ TH1D* l1[15];
+ TH1D* l2[15];
+ //TH1D* N[15];
+ TH1D* hnew[15];
+ TH1D* H1;
+ vector<double> PMean;
+ vector<double> AMean;
+ vector<double> PMeanError;
+ vector<double> AMeanError;
+ vector<double> PSigma;
+ vector<double> PSigmaError;
+ vector<double> ASigmaError;
+ vector<double> ASigma;
+ vector<double> CovElement;
+ //vector<double> SF;
+ //vector<double>SFError;
+
+
+ //c1->Divide(2,2);
+ for(int i=1;i<14;i++) //loop over all simulation output files
+ {
+
+ Char_t input[50];
+ sprintf(input,"NOUTPUT_%d.root",Control[i].Name);
+ TFile* In= new TFile(input);
+ //get histograms from files.
+ h1[i]=(TH1D*)In->Get("1");
+ h2[i]=(TH1D*)In->Get("2");
+ l1[i]=(TH1D*)In->Get("11");
+ l2[i]=(TH1D*)In->Get("12");
+ hnew[i]= (TH1D*)h2[i]->Clone("hnew");
+ hnew[i]->Rebin(Control[i].rebin);
+ cout<<"nbins after cloning= "<<hnew[i]->GetNbinsX()<<endl;
+
+ Profiles[i].PhotonsLead=(TH1D*)In->Get("23");
+ Profiles[i].ElectronsLead=(TH1D*)In->Get("24");
+ Profiles[i].PositronsLead=(TH1D*)In->Get("25");
+ Profiles[i].PhotonsScintillator=(TH1D*)In->Get("26");
+ Profiles[i].ElectronsScintillator=(TH1D*)In->Get("27");
+ Profiles[i].PositronsScintillator=(TH1D*)In->Get("28");
+ Char_t Phl[50];
+ Char_t El[50];
+ Char_t Pol[50];
+ Char_t Phs[50];
+ Char_t Es[50];
+ Char_t Pos[50];
+ sprintf(Phl,"v9.6.p04 Photons in Lead Shower Profile at %.2f GeV",Control[i].Energy);
+ sprintf(El,"v9.6.p04 Electrons in Lead Shower Profile at %.2f GeV",Control[i].Energy);
+ sprintf(Pol,"v9.6.p04 Positrons in Lead Shower Profile at %.2f GeV",Control[i].Energy);
+ sprintf(Phs,"v9.6.p04 Photons in Scintillator Shower Profile at %.2f GeV",Control[i].Energy);
+ sprintf(Es,"v9.6.p04 Electrons in Scintillator Shower Profile at %.2f GeV",Control[i].Energy);
+ sprintf(Pos,"v9.6.p04 Positrons in Scintillator Shower Profile at %.2f GeV;",Control[i].Energy);
+
+ Profiles[i].PhotonsLead->SetTitle(Phl);
+ Profiles[i].ElectronsLead->SetTitle(El);
+ Profiles[i].PositronsLead->SetTitle(Pol);
+ Profiles[i].PhotonsScintillator->SetTitle(Phs);
+ Profiles[i].ElectronsScintillator->SetTitle(Es);
+ Profiles[i].PositronsScintillator->SetTitle(Pos);
+
+ //set titles of histograms
+ Char_t titleActive[50];
+ Char_t titlePassive[50];
+ sprintf(titleActive,"Scintillator E- %f GeV",Control[i].Energy);
+ sprintf(titlePassive,"Lead E - %f GeV",Control[i].Energy);
+ //fit gaussians to histograms
+ TFitResultPtr g= h1[i]->Fit("gaus","S");
+ gStyle->SetOptFit();
+ h1[i]->SetTitle(titlePassive);
+ TF1 *fit1=h1[i]->GetFunction("gaus");
+
+ //get mean, sigma and uncertainties from fit to gaussians for lead layers
+ PMean.push_back(fit1->GetParameter(1));
+ PMeanError.push_back(fit1->GetParError(1));
+ PSigma.push_back(fit1->GetParameter(2));
+ PSigmaError.push_back(fit1->GetParError(2));
+ TFitResultPtr g2;
+ double element;
+ TF1* range = new TF1("range","gaus",0.1,0.3);
+ g2=hnew[i]->Fit("gaus","S");
+ TVirtualFitter *fitter = TVirtualFitter::GetFitter();
+ element=fitter->GetCovarianceMatrixElement(2,1); //get covariance matrix element
+ CovElement.push_back(element);
+ gStyle->SetOptFit();
+ hnew[i]->SetTitle(titleActive);
+ TF1 *fit2=hnew[i]->GetFunction("gaus");
+
+ //get mean and sigma for scintillator layers
+ AMean.push_back(fit2->GetParameter(1));
+ AMeanError.push_back(fit2->GetParError(1));
+ ASigma.push_back(fit2->GetParameter(2));
+ ASigmaError.push_back(fit2->GetParError(2));
+
+ }
+
+
+ double SigmaOverE[15]; //sigma over E
+ double Xaxis[15]; //1 / sqrt(E)
+ double SigmaOverEError[15];
+ double XaxisError[15]; //=0
+
+ double XB[15]; //Fille with Energy specified in .mac file
+ double XBE[15];
+
+ for(int k=0;k<13;k++) //arrays run from 0 to 12.
+ {
+
+ SigmaOverE[k]=ASigma[k]/AMean[k];
+ Xaxis[k]=1/sqrt(Control[k+1].Energy); //change when 100GeV point removed
+ SigmaOverEError[k]=(1/AMean[k])*sqrt(pow(ASigmaError[k],2)+pow((ASigma[k]*AMeanError[k]/(AMean[k])),2)+((2*ASigma[k]*
+CovElement[k])/AMean[k]));
+ cout<<"CovElement= "<<CovElement[k]<<endl;
+
+
+ XaxisError[k]=0;
+ XB[k]=Control[k+1].Energy;
+ XBE[k]=0;
+
+ }
+ TFile* out= new TFile("Save.root","RECREATE");
+ TCanvas*CStraightRes = new TCanvas("CStraightRes","CStraightRes",1800,1000);
+
+ for(int l=1;l<14;++l)
+ {
+ Profiles[l].PhotonsLead->Write();
+ Profiles[l].ElectronsLead->Write();
+ Profiles[l].PositronsLead->Write();
+ Profiles[l].PhotonsScintillator->Write();
+ Profiles[l].ElectronsScintillator->Write();
+ Profiles[l].PositronsScintillator->Write();
+ }
+
+
+ //------------1/sqrt(E)------------------------------------------------------------------------
+ TF1* fitstraight= new TF1("fitstraight","sqrt(([0]*x)**2+(([1])**2))",0,2);
+ fitstraight->SetParName(0,"a");
+ fitstraight->SetParName(1,"c");
+ fitstraight->SetParameter(0,0.01);
+ fitstraight->SetParameter(1,0.01);
+ TGraphErrors*Resolution1=new TGraphErrors(13,Xaxis,SigmaOverE,XaxisError,SigmaOverEError);
+ Resolution1->Draw("AP");
+ Resolution1->SetTitle("10000 events,Physics List=emstandard_opt1noapplycuts;");
+ Resolution1->Fit("fitstraight"); //perform fit
+ Resolution1->GetXaxis()->SetTitle("#frac{1}{#sqrt{E_{Incident}}}[GeV]");
+ Resolution1->GetYaxis()->SetTitle("#frac{#sigma}{E}");
+ Resolution1->GetYaxis()->SetTitleOffset(1.1);
+ /*TText* Res1Text= new TText(0.1,0.1," #Sigma=Resolution of"
+ "Energy Distribution in Scintillator layers E= Total Energy Deposited in"
+ "Active Layers");*/
+ TLatex* Res1Text= new TLatex();
+ CStraightRes->cd();
+ Resolution1->Draw("AP");
+ //Add Labels
+ Res1Text->SetTextSize(0.03);
+ Res1Text->DrawLatex(0.13,0.064,"#sigma=Resolution of Energy Distribution in Active Layers");
+ TText * Res1Line2= new TText();
+ Res1Line2->SetTextSize(0.03);
+ Res1Line2->SetTextSize(0.03);
+ Res1Line2->DrawText(0.13,0.059,"E= Total Energy Deposited in Active Layers");
+ TLatex* FitEquation= new TLatex();
+ FitEquation->DrawLatex(0.13,0.052,"Fit: #frac{#sigma}{E}=#frac{a}{#sqrt{E}} #oplus C");
+ FitEquation->SetTextSize(0.03);
+ gStyle->SetOptStat();
+ TPaveStats*ps= (TPaveStats*)CStraightRes->GetPrimitive("stats");
+ ps->SetX1NDC(0.13);
+ ps->SetX2NDC(0.47);
+ ps->SetY1NDC(0.86);
+ ps->SetY2NDC(0.76);
+ // fitstraight->Draw();
+ CStraightRes->Write();
+ //write to root file
+ Resolution1->Write("Resolution1");
+ //output to text file
+ selected<<"Term,Value,Uncertainty;"<<endl;
+ selected<<"a, "<<setprecision(3)<<fitstraight->GetParameter(0)<<" , "<<setprecision(1)<<fitstraight->GetParError(0)<<";"<<endl;
+ selected<<"c, "<<setprecision(1)<<fitstraight->GetParameter(1)<<" , "<<setprecision(1)<<fitstraight->GetParError(1)<<";"<<endl;
+ selected<<"\n"<<endl;
+ selected<<"Energy,Sampling Fraction,Uncertainty"<<endl;
+
+ //-------------------------------Curved Fit----------------------------------
+ TGraphErrors*Resolution2=new TGraphErrors(13,XB,SigmaOverE,XBE,SigmaOverEError);
+ TCanvas*CResolution2=new TCanvas("CResolution2","CResolution2",1800,1000);
+
+ TF1* fit1= new TF1("fit1","sqrt(([0]/sqrt(x))**2+(([1])**2))",0,2);
+ fit1->SetParName(0,"a");
+ fit1->SetParName(1,"c");
+
+ fit1->SetParameter(0,0.01);
+ fit1->SetParameter(1,0.01);
+ fit1->SetLineColor(kBlue);
+ Resolution2->Fit("fit1");
+
+ Resolution2->SetTitle("10000 events, Physics List=emstandard_opt1noapplycuts");
+ Resolution2->GetXaxis()->SetTitle("E_{Incident}[GeV]");
+ Resolution2->GetYaxis()->SetTitle("#frac{#sigma}{E}");
+ Resolution2->Draw("AP");
+ gStyle->SetOptStat();
+ TLatex* Res2Text= new TLatex();
+ Res1Text->SetTextSize(0.03);
+ Res1Text->DrawLatex(23.5,0.062,"#sigma=Resolution of Energy Distribution in Active Layers");
+ TText * Res2Line2= new TText();
+ Res2Line2->SetTextSize(0.03);
+ Res2Line2->SetTextSize(0.03);
+ Res2Line2->DrawText(23.5,0.057,"E= Total Energy Deposited Per Event in Active Layers");
+ TLatex* FitEquation2= new TLatex();
+ FitEquation2->DrawLatex(23.5,0.052,"Fit: #frac{#sigma}{E}=#frac{a}{#sqrt{E}} #oplus C");
+ FitEquation2->SetTextSize(0.03);
+ //write to root file
+ Resolution2->Write("Resolution2");
+ CResolution2->Write();
+ //numerical results are inferred from straight line plot but are consistent with cureved plot.
+
+ //----------------------------shower profiles-----------------------------------------
+ double En[15];
+ double SFError[15];
+ double SF[15];
+ double SFErrorProp[15];
+ double LeadSF[15];
+ double ScintSF[15];
+ double LeadSFError[15];
+ double ScintSFError[15];
+ double Energ[15];
+ double EnergyEr[15];for(int n=1;n<14;++n){
+ int k=n-1;
+ Double_t xmin=l2[n]->GetXaxis()->GetXmin();
+ Double_t xmax=l2[n]->GetXaxis()->GetXmax();
+ Double_t integral=l2[n]->Integral(xmin,xmax);
+ Double_t xminlead=l1[n]->GetXaxis()->GetXmin();
+ Double_t xmaxlead=l1[n]->GetXaxis()->GetXmax();
+ Double_t integrallead=l1[n]->Integral(xminlead,xmaxlead);
+ Double_t leadGeV=integrallead/1000.0;
+ Double_t ScGeV=integral/1000.0;
+
+
+ gStyle->SetOptStat(0);
+
+ // l1[n]->GetYaxis()->SetTitle("#frac{dE}{dX}");
+ // l1[n]->GetXaxis()->SetTitle("Distance of Material Traversed/Layers");
+ // l1[n]->Draw();
+ // l2[n]->Draw("Same");
+
+
+ //Leg[n]->AddEntry(l1[n],"Lead","P");
+ //Leg[n]->AddEntry(l2[n],"Scintillator","P");
+
+ En[k]=Control[n].Energy;
+ //Sampling->Fit("pol1")`
+
+ SF[k]=integral/integrallead;
+ LeadSF[k]=leadGeV/(Control[n].Energy);
+ ScintSF[k]=ScGeV/(Control[n].Energy);
+ LeadSFError[k]=(PSigma[k]/sqrt(nEvents))*(1/Control[n].Energy);
+ ScintSFError[k]=(ASigma[k]/sqrt(nEvents))*(1/Control[n].Energy);
+ //SFError[k]=SF[k]*sqrt(pow(((ASigma[k]/sqrt(nEvents))/ScGeV),2)+pow(((PSigma[k]/sqrt(nEvents)))/leadGeV,2));
+ SFError[k]=(ASigma[k]/sqrt(nEvents))*(1/leadGeV);
+ ErrorInv<<"SFError= "<<SFError[k]<<"k= "<<k<<endl;
+ ErrorInv<<"Asigma= "<<ASigma[k]<<"for k= "<<k<<endl;
+
+
+ cout<<"Active Integral for Histogram "<<n<<"= "<<ScGeV<<endl;
+ cout<<"Absorber Integral for Histogram "<<n<<"= "<<leadGeV<<endl;
+ cout<<"sampling fraction for energy "<<Control[n].Energy<<"GeV = "<<ScGeV/leadGeV<<endl;
+ cout<<"--------------------------------"<<endl;
+ outt<<"Active Integral for Histogram "<<n<<"= "<<ScGeV<<endl;
+ outt<<"Absorber Integral for Histogram "<<n<<"= "<<leadGeV<<endl;
+ outt<<"sampling fraction for energy "<<Control[n].Energy<<"GeV = "<<ScGeV/leadGeV<<endl;
+ outt<<"--------------------------------"<<endl;
+ selected<<setprecision(4)<<Control[n].Energy<<" , "<<setprecision(4)<<SF[k]<<" , "<<setprecision(2)<<SFError[k]<<";"<<endl;
+ }
+ TGraphErrors* ScintSampling =new TGraphErrors(13,En,ScintSF,0,ScintSFError);
+ ScintSampling->Write("ScintSampling");
+ TGraphErrors* LeadSampling = new TGraphErrors(13,En,LeadSF,0,LeadSFError);
+ LeadSampling->Write("LeadSampling");
+ //create sampling fraction plot and write it to file.
+ TGraphErrors*Sampling=new TGraphErrors(13,En,SF,0,SFError);
+ Sampling->SetTitle("Sampling Fraction as a Funciton of Energy, Phys List=emstandard_opt1noapplycuts,10000 events");
+ Sampling->GetXaxis()->SetTitle("E/GeV");
+ Sampling->GetYaxis()->SetTitle("#frac{Scintillator E}{Lead E}");
+ // Sampling->Fit("pol1");
+ Sampling->Write("Sampling");
+ TCanvas* CSampling = new TCanvas("CSampling","CSampling",1800,1000);
+ Sampling->GetYaxis()->SetTitleOffset(1.2);
+ Sampling->Draw("AP");
+ gPad->SetLogx();
+ CSampling->Write();
+
+
+ TLegend*Leg[15];
+ //write longitudinal showers to file and sort formatting for shower profiles
+ for(int j=1;j<14;j++)
+ {
+ Char_t canvasname[20];
+ sprintf(canvasname,"c_%d",j);
+ c[j]=new TCanvas(canvasname,canvasname,1800,1000);
+
+ Char_t title1[50];
+ Char_t title2[50];
+ Char_t longitudinaltitle[50];
+ sprintf(title1,"1_%d",Control[j].Name);
+ sprintf(title2,"2_%d",Control[j].Name);
+ sprintf(longitudinaltitle,"Energy deposited per event, %5.2f GeV incident energy",Control[j].Energy);
+ h1[j]->Write(title1);
+ h2[j]->Write(title2);
+ hnew[j]->Write("Hnew");
+ l2[j]->SetLineColor(kRed);
+ l1[j]->SetLineColor(kBlue);
+ l1[j]->SetTitle(longitudinaltitle);
+ l2[j]->SetTitle(longitudinaltitle);
+ l1[j]->GetXaxis()->SetTitle("Layer number");
+ l2[j]->GetXaxis()->SetTitle("Layer number");
+ l1[j]->GetYaxis()->SetTitle("#frac{MeV}{Layer}");
+ l2[j]->GetYaxis()->SetTitle("#frac{MeV}{Layer}");
+ l1[j]->Write();
+ l2[j]->Write();
+ c[j]->cd();
+ l1[j]->Draw();
+ l2[j]->Draw("SAME");
+ Leg[j]=new TLegend(0.75,0.60,0.95,0.75);
+ Leg[j]->AddEntry(l1[j],"Lead","L");
+ Leg[j]->AddEntry(l2[j],"Scintillator","L");
+ Leg[j]->Draw();
+ c[j]->SetLogy();
+ c[j]->Write();
+
+ //N[j]->Write();
+ // hnew[j]->Write();
+ }
+
+
+
+
+ for(int E=0;E<13;++E)
+ {
+ ErrorInv<<"Psigma for E = "<<E<<" "<<PSigma.at(E)<<endl;
+ }
+
+ //output selected results to file
+ TFile* Selectedroot= new TFile("Selectedresults.root","RECREATE");
+ CStraightRes->Write();
+ CResolution2->Write();
+ CSampling->Write();
+ for(int n=1;n<14;++n){
+ c[n]->Write();
+ }
+ selected.close();
+ outt.close();
+ //cout<<"this was compiled with my script"<<endl;
+
+}
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/TestEm3.in b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/TestEm3.in
deleted file mode 100644
index 64cba5f..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/TestEm3.in
+++ /dev/null
@@ -1,36 +0,0 @@
-# $Id: TestEm3.in,v 1.27 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "exampleN03.cc"
-# (can be run in batch, without graphic)
-#
-# Lead-liquidArgon 50 layers; electron 1 GeV
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/phys/addPhysics emstandard_opt0
-#
-/run/initialize
-/process/em/verbose 1
-/process/msc/Skin 1
-#
-/gun/particle e-
-/gun/energy 1 GeV
-#
-###/analysis/setFileName testem3
-###/analysis/h1/set 1 100 700. 1000. MeV #edep in absor1
-###/analysis/h1/set 2 100 0. 300. MeV #edep in absor2
-###/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-###/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-###/analysis/h1/set 21 102 0. 102. none #energy flow
-###/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-#
-/testem/event/printModulo 50
-#
-/run/beamOn 100
-#
-/process/em/setSecBiasing eIoni World 0.3 2 MeV
-/process/em/setSecBiasing eBrem World 0.2 1 MeV
-/run/physicsModified
-/run/beamOn 100
-#
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/atlashec.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/atlashec.mac
deleted file mode 100644
index 5476a54..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/atlashec.mac
+++ /dev/null
@@ -1,32 +0,0 @@
-# $Id: atlashec.mac,v 1.8 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-#
-# Atlas HEC model;
-#
-/control/verbose 2
-/run/verbose 1
-#
-/testem/det/setNbOfLayers 25
-/testem/det/setNbOfAbsor 2
-/testem/det/setAbsor 1 liquidArgon 0.85 cm
-/testem/det/setAbsor 2 Copper 2.5 cm
-/testem/det/setSizeYZ 17 cm
-#
-/testem/phys/addPhysics emstandard_opt0
-/testem/phys/setCuts 0.7 mm
-#
-/run/initialize
-#
-/analysis/setFileName 10gev
-/analysis/h1/set 1 100 0. 1000. MeV
-/analysis/h1/set 2 100 9. 10. GeV
-#
-/gun/particle e-
-/gun/particle pi+
-/gun/energy 100 GeV
-/gun/position -40.635 0 0
-#
-/testem/event/printModulo 100
-/run/beamOn 100
-
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/dedx.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/dedx.mac
deleted file mode 100644
index 31f6b5e..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/dedx.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: dedx.mac,v 1.5 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-# (can be run in batch, without graphic)
-#
-# To control dE/dx calculation:
-# 1 layer; minimum ionizing particle
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setNbOfLayers 1
-/testem/det/setNbOfAbsor 1
-/testem/det/setAbsor 1 liquidArgon 1 cm
-/testem/det/setSizeYZ 1 cm
-#
-/testem/phys/addPhysics local
-/testem/phys/setCuts 1 km
-#
-/run/initialize
-#
-/process/inactivate msc
-/process/eLoss/fluct false
-#
-/testem/gun/setDefault
-/gun/particle mu-
-/gun/energy 50 GeV
-#
-/tracking/verbose 1
-#
-/run/beamOn 2
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/emtutor.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/emtutor.mac
deleted file mode 100644
index b21e545..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/emtutor.mac
+++ /dev/null
@@ -1,20 +0,0 @@
-# $Id: emtutor.mac,v 1.1 2007-04-06 17:41:54 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-#
-# water 1 layer;
-#
-/control/verbose 2
-#
-/testem/det/setNbOfLayers 1
-/testem/det/setNbOfAbsor 1
-/testem/det/setAbsor 1 Aluminium 10 cm
-/testem/det/setSizeYZ 10 cm
-#
-/testem/phys/setCuts 2 cm
-#
-/run/initialize
-#
-/control/execute vis.mac
-#
-/testem/gun/setDefault
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/ionC12.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/ionC12.mac
deleted file mode 100644
index de7d605..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/ionC12.mac
+++ /dev/null
@@ -1,36 +0,0 @@
-# $Id: ionC12.mac,v 1.11 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-# (can be run in batch, without graphic)
-#
-# CsI 1 layer;
-#
-/control/verbose 2
-/run/verbose 2
-/tracking/verbose 1
-#
-# select EM PhysicsList
-#
-/testem/phys/addPhysics local
-#
-/testem/phys/setCuts 1 mm
-#
-/testem/det/setNbOfLayers 1
-/testem/det/setNbOfAbsor 1
-/testem/det/setAbsor 1 CsI 1.99 cm
-/testem/det/setSizeYZ 3 cm
-#
-/gun/particle ion
-/gun/ion 6 12
-/run/initialize
-#
-/testem/gun/setDefault
-#/gun/particle ion
-#/gun/ion 6 12
-/gun/energy 108 GeV
-#
-/analysis/h1/set 1 200 350. 550. MeV
-#
-/testem/event/printModulo 1000
-/run/beamOn 5000
-
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lhcb.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lhcb.mac
deleted file mode 100644
index b9300a8..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lhcb.mac
+++ /dev/null
@@ -1,131 +0,0 @@
-# $Id: lhcb.mac,v 1.4 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-#
-# LHCB ECAL model
-#
-/control/verbose 2
-/run/verbose 1
-#
-/testem/det/setNbOfLayers 66
-/testem/det/setNbOfAbsor 2
-/testem/det/setAbsor 1 Lead 2 mm
-/testem/det/setAbsor 2 Scintillator 4 mm
-/testem/det/setSizeYZ 20 cm
-#
-/testem/phys/addPhysics emstandard_opt0
-#
-/gun/particle e-
-/gun/energy 10 GeV
-#
-/testem/event/printModulo 50
-#
-/testem/phys/setCuts 3 mm
-/run/initialize
-/analysis/setFileName OUTPUT_3.000
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 1 mm
-/run/initialize
-/analysis/setFileName OUTPUT_1.000
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.5 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.500
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.1 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.100
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.07 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.070
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.030 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.030
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.01 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.010
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.007 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.007
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.003 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.003
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
-#
-/testem/phys/setCuts 0.001 mm
-/run/initialize
-/analysis/setFileName OUTPUT_0.001
-/analysis/h1/set 1 500 5. 10. GeV
-/analysis/h1/set 2 500 0. 5. GeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-/run/beamOn 5000
\ No newline at end of file
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/linac.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/linac.mac
deleted file mode 100644
index 8cdb04d..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/linac.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: linac.mac,v 1.6 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3"
-#
-# Linac/Ecal from Graham Wilson;
-#
-/control/verbose 2
-/run/verbose 1
-#
-/testem/det/setNbOfLayers 42
-/testem/det/setNbOfAbsor 4
-/testem/det/setAbsor 1 Tungsten 2.50 mm
-/testem/det/setAbsor 2 NemaG10 0.68 mm
-/testem/det/setAbsor 3 Silicon 0.32 mm
-/testem/det/setAbsor 4 Air20 0.25 mm
-/testem/det/setSizeYZ 30 cm
-#
-/testem/phys/addPhysics local
-/testem/phys/setCuts 10 um
-#
-/run/initialize
-#
-/analysis/setFileName linac
-/analysis/h1/set 3 100 0. 20. MeV #edep in absor3
-#
-/gun/particle gamma
-/gun/energy 1 GeV
-#
-/testem/event/printModulo 100
-#
-/run/beamOn 1000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lockwood.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lockwood.mac
deleted file mode 100644
index 7a37b09..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/lockwood.mac
+++ /dev/null
@@ -1,35 +0,0 @@
-# $Id: lockwood.mac,v 1.6 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-# (can be run in batch, without graphic)
-#
-# Al-Au-Al 1 layer;
-#
-# G.L.Lockwood et al. SAND79-0414 (1980)
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/phys/addPhysics local
-#
-/testem/det/setNbOfLayers 1
-/testem/det/setNbOfAbsor 3
-/testem/det/setAbsor 1 Aluminium 170 um
-/testem/det/setAbsor 2 Gold 23 um
-/testem/det/setAbsor 3 Aluminium 850 um
-/testem/det/setSizeYZ 1050 um
-#
-/testem/phys/setCuts 1 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 1 MeV
-/testem/gun/rndm 0.8
-#
-/testem/stepMax/absorber 2 5 um
-#
-/testem/event/printModulo 1000
-#
-/run/beamOn 10
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/opt1noapplycuts.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/opt1noapplycuts.mac
new file mode 100644
index 0000000..cdd282b
--- /dev/null
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/opt1noapplycuts.mac
@@ -0,0 +1,307 @@
+# $Id: lhcb.mac,v 1.4 2009-09-17 20:06:26 maire Exp $
+#
+# Macro file for "TestEm3.cc"
+#
+# LHCB ECAL model
+#
+/control/verbose 2
+/run/verbose 1
+#
+/testem/det/setNbOfLayers 66
+/testem/det/setNbOfAbsor 2
+/testem/det/setAbsor 1 Lead 2 mm
+/testem/det/setAbsor 2 Scintillator 4 mm
+/testem/det/setSizeYZ 20 cm
+#
+/testem/phys/addPhysics emstandard_opt1nocuts
+#
+#/gun/particle e-
+#/gun/energy 100 GeV
+#
+#/testem/event/printModulo 100
+#
+#/testem/phys/setCuts 5.0 mm
+#/run/initialize
+#/analysis/setFileName NOUTPUT_1
+#/analysis/h1/set 1 500 60. 100. GeV
+#/analysis/h1/set 2 500 5. 35. GeV
+#/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+#/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+#/analysis/h1/set 21 102 0. 102. none #energy flow
+#/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+#/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 44.44 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_2
+/analysis/h1/set 1 500 26. 45. GeV
+/analysis/h1/set 2 500 2.22 16. GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 25.0 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_3
+/analysis/h1/set 1 500 15.00 25.00 GeV
+/analysis/h1/set 2 500 1.25 8.75 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 16.00 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_4
+/analysis/h1/set 1 500 9.60 16.00 GeV
+/analysis/h1/set 2 500 0.80 5.60 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 11.11 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_5
+/analysis/h1/set 1 500 6.67 11.11 GeV
+/analysis/h1/set 2 500 0.56 3.89 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 8.16 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_6
+/analysis/h1/set 1 500 4.90 8.16 GeV
+/analysis/h1/set 2 500 0.41 2.86 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 6.25 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_7
+/analysis/h1/set 1 500 3.75 6.25 GeV
+/analysis/h1/set 2 500 0.31 2.19 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 4.94 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_8
+/analysis/h1/set 1 500 2.96 4.94 GeV
+/analysis/h1/set 2 500 0.25 1.73 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 4.00 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_9
+/analysis/h1/set 1 500 2.4 4.0 GeV
+/analysis/h1/set 2 500 0.2 1.4 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 3.31 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_10
+/analysis/h1/set 1 500 1.98 3.31 GeV
+/analysis/h1/set 2 500 0.17 1.16 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 2.78 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_11
+/analysis/h1/set 1 500 1.67 2.78 GeV
+/analysis/h1/set 2 500 0.14 0.97 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 2.37 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_12
+/analysis/h1/set 1 500 1.42 2.37 GeV
+/analysis/h1/set 2 500 0.12 0.83 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 2.04 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_13
+/analysis/h1/set 1 500 1.22 2.04 GeV
+/analysis/h1/set 2 500 0.1 0.71 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 1.78 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_14
+/analysis/h1/set 1 500 1.07 1.78 GeV
+/analysis/h1/set 2 500 0.09 0.62 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
+#
+/gun/particle e-
+/gun/energy 1.56 GeV
+#
+/testem/phys/setCuts 5.0 mm
+/run/initialize
+/analysis/setFileName NOUTPUT_15
+/analysis/h1/set 1 500 0.94 1.56 GeV
+/analysis/h1/set 2 500 0.08 0.55 GeV
+/analysis/h1/set 11 66 0. 66. none #long. profile in absor1
+/analysis/h1/set 12 66 0. 66. none #long. profile in absor2
+/analysis/h1/set 21 102 0. 102. none #energy flow
+/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
+/analysis/h1/setActivation 23 true
+/analysis/h1/setActivation 24 true
+/analysis/h1/setActivation 25 true
+/analysis/h1/setActivation 26 true
+/analysis/h1/setActivation 27 true
+/analysis/h1/setActivation 28 true
+/run/beamOn 10000
\ No newline at end of file
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/plotHisto.C b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/plotHisto.C
deleted file mode 100644
index 46991c2..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/plotHisto.C
+++ /dev/null
@@ -1,31 +0,0 @@
-{
- gROOT->Reset();
-
- // Draw histos filled by Geant4 simulation
- //
- ///TFile f = TFile("testem3.root");
- TFile f = TFile("run01.root");
- TCanvas* c1 = new TCanvas("c1", " ");
-
- TH1D* hist1 = (TH1D*)f.Get("1");
- hist1->Draw("HIST");
-
- TH1D* hist2 = (TH1D*)f.Get("2");
- hist2->Draw("HIST");
-
- TH1D* hist3 = (TH1D*)f.Get("11");
- c1->SetLogy(1);
- c1->cd();
- c1->Update();
- hist3->Draw("HIST");
-
- TH1D* hist4 = (TH1D*)f.Get("12");
- hist4->Draw("HIST");
-
- TH1D* hist5 = (TH1D*)f.Get("21");
- hist5->Draw("HIST");
-
- TH1D* hist6 = (TH1D*)f.Get("22");
- hist6->Draw("HIST");
-
-}
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/retrieveTables.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/retrieveTables.mac
deleted file mode 100644
index 5fea5b1..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/retrieveTables.mac
+++ /dev/null
@@ -1,25 +0,0 @@
-# $Id: retrieveTables.mac,v 1.6 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "exampleN03.cc"
-#
-# can be run in batch, without graphic
-# or interactively: Idle> /control/execute run1.mac
-#
-/control/verbose 2
-#
-# select EM PhysicsList
-#
-/testem/phys/addPhysics emstandard_opt0
-#
-/run/particle/retrievePhysicsTable physdata
-/run/initialize
-/process/eLoss/verbose 1
-/process/eLoss/CSDARange true
-#
-/gun/particle e-
-/gun/energy 1 GeV
-#
-/run/beamOn 3
-#
-/run/particle/storePhysicsTable physdata1
-#
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run.sh b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run.sh
new file mode 100755
index 0000000..4722bd0
--- /dev/null
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run.sh
@@ -0,0 +1,27 @@
+#!/bin/bash
+
+CurrentDir="${PWD##*/}"
+echo $CurrentDir
+TestDir=CaloTest
+if [[ ! -d "CaloTest" && $CurrentDir != $TestDir ]]; then
+mkdir CaloTest
+cd CaloTest
+fi
+
+if [[ -d "CaloTest" && $CurrentDir != $TestDir ]]; then
+cd CaloTest
+fi
+
+testEm3.exe $G4TESTEM3ROOT/scripts/opt1noapplycuts.mac
+
+#There is a bug with ROOT 6 opening Geant4 created histograms, have to use old root for now
+source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.30/$CMTCONFIG/root/bin/thisroot.sh
+
+g++ -c `root-config --cflags` $G4TESTEM3ROOT/scripts/Plot.C
+g++ -o Plot `root-config --glibs` Plot.o
+
+./Plot
+
+echo "Results Produced in Save.root, Selectedresults.root and selectedresults.txt"
+
+
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run01.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run01.mac
deleted file mode 100644
index 4c69c0e..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run01.mac
+++ /dev/null
@@ -1,30 +0,0 @@
-# $Id: run01.mac,v 1.19 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "exampleN03.cc"
-# (can be run in batch, without graphic)
-#
-# Lead-liquidArgon 50 layers; electron 1 GeV
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/phys/addPhysics local
-#
-/testem/phys/setCuts 100 um
-#
-/run/initialize
-#
-/gun/particle e-
-/gun/energy 1 GeV
-#
-/analysis/setFileName run01
-/analysis/h1/set 1 100 700. 1000. MeV
-/analysis/h1/set 2 100 0. 300. MeV
-/analysis/h1/set 11 52 0. 52. none #long. profile in absor1
-/analysis/h1/set 12 52 0. 52. none #long. profile in absor2
-/analysis/h1/set 21 102 0. 102. none #energy flow
-/analysis/h1/set 22 102 0. 102. none #lateral energy leakage
-#
-/testem/event/printModulo 100
-#
-/run/beamOn 1000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run02.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run02.mac
deleted file mode 100644
index adaf939..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/run02.mac
+++ /dev/null
@@ -1,32 +0,0 @@
-# $Id: run02.mac,v 1.16 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-# (can be run in batch, without graphic)
-#
-# Tungsten-Silicon 50 layers; electron 1 GeV
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/phys/addPhysics local
-#
-/testem/det/setNbOfLayers 50
-/testem/det/setNbOfAbsor 2
-/testem/det/setAbsor 1 Tungsten 2 mm
-/testem/det/setAbsor 2 Silicon 0.3 mm
-#
-/testem/phys/setCuts 100 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 1 GeV
-#
-/analysis/setFileName run02
-/analysis/h1/set 1 100 900. 1000. MeV
-/analysis/h1/set 2 100 0. 50. MeV
-#
-/testem/event/printModulo 50
-#
-/run/beamOn 200
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/storeTables.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/storeTables.mac
deleted file mode 100644
index 4a55c22..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/storeTables.mac
+++ /dev/null
@@ -1,28 +0,0 @@
-# $Id: storeTables.mac,v 1.10 2009-09-17 20:06:26 maire Exp $
-#
-# Macro file for "testEm3.cc"
-#
-# can be run in batch, without graphic
-# or interactively: Idle> /control/execute run1.mac
-#
-/control/verbose 2
-#
-/testem/phys/addPhysics emstandard_opt0
-#
-/run/initialize
-/process/eLoss/verbose 1
-/process/eLoss/CSDARange true
-/run/beamOn 3
-#
-/run/particle/storePhysicsTable physdata
-#
-/gun/particle e-
-/gun/energy 1 GeV
-#
-#/run/beamOn 3
-#
-/run/particle/retrievePhysicsTable physdata
-#
-/random/resetEngineFrom currentRun.rndm
-/run/beamOn 3
-#
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/tileCal.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/tileCal.mac
deleted file mode 100644
index 77496fb..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/tileCal.mac
+++ /dev/null
@@ -1,30 +0,0 @@
-# $Id: tileCal.mac,v 1.1 2011/01/03 13:54:15 maire Exp $
-#
-# Macro file for "TestEm3.cc"
-#
-# ATLAS tileCal
-#
-/control/verbose 2
-/run/verbose 1
-#
-/testem/det/setNbOfLayers 50
-/testem/det/setNbOfAbsor 2
-/testem/det/setAbsor 1 Iron 14 mm
-/testem/det/setAbsor 2 Scintillator 3 mm
-/testem/det/setSizeYZ 20 cm
-#
-/testem/phys/addPhysics emstandard_opt0
-/testem/phys/setCuts 1 um
-#
-/run/initialize
-#
-/analysis/setFileName mu.1mm
-/analysis/h1/set 1 100 500. 1000. MeV
-/analysis/h1/set 2 100 0. 50. MeV
-#
-/gun/particle mu+
-/gun/energy 20 GeV
-#
-/testem/event/printModulo 4000
-#
-/run/beamOn 20000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/vis.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/vis.mac
deleted file mode 100644
index 3cd77e0..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm3/scripts/vis.mac
+++ /dev/null
@@ -1,79 +0,0 @@
-#
-# Macro file for the initialization phase of "TestEm3.cc"
-#
-# Sets some default verbose
-# and initializes the graphic.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/run/initialize
-#
-# Use this open statement to create an OpenGL view:
-/vis/open OGL 600x600-0+0
-#
-# Use this open statement to create a .prim file suitable for
-# viewing in DAWN:
-#/vis/open DAWNFILE
-#
-# Use this open statement to create a .heprep file suitable for
-# viewing in HepRApp:
-#/vis/open HepRepFile
-#
-# Use this open statement to create a .wrl file suitable for
-# viewing in a VRML viewer:
-#/vis/open VRML2FILE
-#
-# Disable auto refresh and quieten vis messages whilst scene and
-# trajectories are established:
-/vis/viewer/set/autoRefresh false
-/vis/verbose errors
-#
-# Draw geometry:
-/vis/drawVolume
-#
-# Specify view angle:
-#/vis/viewer/set/viewpointThetaPhi 90. 0.
-#
-# Specify zoom value:
-/vis/viewer/zoom 1.4
-#
-# Specify style (surface or wireframe):
-#/vis/viewer/set/style wireframe
-#
-# Draw coordinate axes:
-#/vis/scene/add/axes 0 0 0 1 m
-#
-# Draw smooth trajectories at end of event, showing trajectory points
-# as markers 2 pixels wide:
-/vis/scene/add/trajectories smooth
-/vis/modeling/trajectories/create/drawByCharge
-#/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
-#/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 1
-# (if too many tracks cause core dump => /tracking/storeTrajectory 0)
-#
-# Draw hits at end of event:
-#/vis/scene/add/hits
-#
-# To draw only gammas:
-#/vis/filtering/trajectories/create/particleFilter
-#/vis/filtering/trajectories/particleFilter-0/add gamma
-#
-# To invert the above, drawing all particles except gammas,
-# keep the above two lines but also add:
-#/vis/filtering/trajectories/particleFilter-0/invert true
-#
-# Many other options are available with /vis/modeling and /vis/filtering.
-# For example, to select colour by particle ID:
-#/vis/modeling/trajectories/create/drawByParticleID
-#/vis/modeling/trajectories/drawByParticleID-0/set e- blue
-#
-# To superimpose all of the events from a given run:
-/vis/scene/endOfEventAction accumulate
-#
-# Re-establish auto refreshing and verbosity:
-/vis/viewer/set/autoRefresh true
-/vis/verbose warnings
-#
-# For file-based drivers, use this to create an empty detector view:
-#/vis/viewer/flush
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/cmt/requirements b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/cmt/requirements
index 1df27ba..fc1abdf 100755
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/cmt/requirements
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/cmt/requirements
@@ -29,6 +29,10 @@ application testEm5 ../TestEm5.cc
#========================================================================
+#--------Install scripts-------
+apply_pattern install_scripts
+alias 'MultipleScattering_test' 'Em5run.sh'
+
private
apply_pattern G4_copy_examples_source
macro_append testEm5_dependencies " G4TestEm5 "
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/Em5run.sh b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/Em5run.sh
new file mode 100755
index 0000000..c017106
--- /dev/null
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/Em5run.sh
@@ -0,0 +1,30 @@
+#!/bin/bash
+if [[ $# -eq 0 ]] ; then
+ echo 'Running Multiple Scattering Test'
+ echo 'Number of runs for uncertainties not specified, using defauly value of 1000'
+ n=1000
+fi
+
+if [[ $# -eq 1 ]] ; then
+ n=$1
+ echo 'Running Multiple Scattering test'
+fi
+
+CurrentDir="${PWD##*/}"
+echo $CurrentDir
+TestDir=MSCTest
+if [[ ! -d "MSCTest" && $CurrentDir != $TestDir ]]; then
+mkdir $TestDir
+cd $TestDir
+fi
+
+if [[ -d "MSCTest" && $CurrentDir != $TestDir ]]; then
+cd $TestDir
+fi
+
+MacGen.py $n
+testEm5.exe BruteForce.mac
+
+
+
+
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/MacGen.py b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/MacGen.py
new file mode 100755
index 0000000..3291495
--- /dev/null
+++ b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/MacGen.py
@@ -0,0 +1,13 @@
+#!/usr/bin/env python
+
+import sys
+n=int(sys.argv[1])
+macout=file('BruteForce.mac','w')
+macout.write('/control/verbose 1 \n/run/verbose 1 \n/testem/det/setAbsMat Silicon \n/testem/det/setAbsThick 300 um \n/testem/det/setAbsYZ 42 mm \n/testem/phys/addPhysics emstandard_opt1nocuts \n/testem/phys/setCuts 5 mm \n')
+Energies=[1,2,3,4,5,7,9,12,15,20,25,30,40]
+Limits=[4.0,3.0,2.0,1.5,1.0,0.6,0.5,0.5,0.4,0.2,0.2,0.15,0.15]
+for E in Energies:
+ j=0
+ for i in range(0,n):
+ macout.write('/run/initialize \n/testem/gun/setDefault \n/gun/particle e- \n/gun/energy %.2f GeV \n/analysis/setFileName Output_%.2f \n/analysis/h1/set 1 100 0 100 keV #energy depostied in absorber \n/analysis/h1/set 10 100 0 100 keV #K.E at exit of world \n/analysis/h1/set 11 100 0 100 keV #Energy fluence dE/dOmega \n/analysis/h1/set 12 100 0.0 0.09 mrad #space angle dN/dOmega \n/analysis/h1/set 13 100 -%.2f %.2f mrad #projected angle at exit of world \n/analysis/h1/set 14 100 -14 14 nm #projected positon at exit of world \n/analysis/h1/set 15 100 0 44 mm #radius at exit of world \n/run/beamOn 10000 \n'%(E,E,Limits[j],Limits[j]))
+ j=j+1
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.in b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.in
deleted file mode 100644
index 40c3bc2..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.in
+++ /dev/null
@@ -1,63 +0,0 @@
-# $Id: TestEm5.in,v 1.23 2009-01-23 13:26:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silicon
-/testem/det/setAbsThick 20 um
-/testem/det/setAbsYZ 50 um
-#
-/testem/phys/addPhysics emstandard_opt3
-#
-/testem/phys/setCuts 7 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle pi+
-/gun/energy 5 GeV
-#
-/run/beamOn 100000
-#
-/testem/det/setAbsMat Gold
-/testem/det/setAbsThick 9.658 um
-/testem/det/setAbsYZ 1 mm
-/testem/det/update
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 15.7 MeV
-/run/beamOn 50000
-#
-/testem/phys/setCuts 10 um
-#
-/testem/det/setAbsMat G4_POLYETHYLENE
-/testem/det/setAbsThick 4.74 mm
-/testem/det/setAbsYZ 50 mm
-/testem/det/update
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 96.2 MeV
-/run/beamOn 50000
-#
-/testem/phys/setCuts 1 mm
-#
-/testem/stack/killSecondaries
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 8.004 mm
-/testem/det/update
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 174.664 GeV
-/run/beamOn 50000
-#
-/testem/det/setAbsMat Iron
-/testem/det/setAbsThick 1 m
-/testem/det/update
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 100 GeV
-/run/beamOn 50000
-#
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.out b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.out
deleted file mode 100644
index 44346ba..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/TestEm5.out
+++ /dev/null
@@ -1,2368 +0,0 @@
-
- ############################################
- !!! WARNING - FPE detection is activated !!!
- ############################################
-
-*************************************************************
- Geant4 version Name: geant4-09-06-patch-04 (30-January-2015)
- Copyright : Geant4 Collaboration
- Reference : NIM A 506 (2003), 250-303
- WWW : http://cern.ch/geant4
-*************************************************************
-
-/run/verbose 2
-#
-/testem/det/setAbsMat Silicon
-/testem/det/setAbsThick 20 um
-/testem/det/setAbsYZ 50 um
-#
-/testem/phys/addPhysics emstandard_opt3
-PhysicsList::AddPhysicsList: <emstandard_opt3>
-#
-/testem/phys/setCuts 7 um
-#
-/run/initialize
-userDetector->Construct() start.
-
- Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
- ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: 1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: 2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- Material: Silicon density: 2.330 g/cm3 RadL: 9.368 cm Nucl.Int.Length: 45.761 cm Imean: 173.000 eV
- ---> Element: Silicon ( ) Z = 14.0 N = 28.1 A = 28.09 g/mole
- ---> Isotope: 28 Z = 14 N = 28 A = 27.98 g/mole abundance: 92.23 %
- ---> Isotope: 29 Z = 14 N = 29 A = 28.98 g/mole abundance: 4.68 %
- ---> Isotope: 30 Z = 14 N = 30 A = 29.97 g/mole abundance: 3.09 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- The WORLD is made of 30 um of Galactic. The transverse size (YZ) of the world is 60 um
- The ABSORBER is made of 20 um of Silicon. The transverse size (YZ) is 50 um
- X position of the middle of the absorber 0 fm
-World is registered to the default region.
-physicsList->Construct() start.
-physicsList->Construct() start.
-physicsList->CheckParticleList() start.
-physicsList->setCut() start.
-PhysicsList::SetCuts:CutLength : 1 mm
-#
-/testem/gun/setDefault
-/gun/particle pi+
-/gun/energy 5 GeV
-#
-/run/beamOn 100000
-
-phot: for gamma SubType= 12
- LambdaPrime table from 200 keV to 10 TeV in 154 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila FluoActive
-
-compt: for gamma SubType= 13
- Lambda table from 10 eV to 1 MeV in 100 bins, spline: 1
- LambdaPrime table from 1 MeV to 10 TeV in 140 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- KleinNishina : Emin= 0 eV Emax= 10 TeV FluoActive
-
-conv: for gamma SubType= 14
- Lambda table from 1.022 MeV to 10 TeV in 140 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BetheHeitler : Emin= 0 eV Emax= 80 GeV
- BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV
-
-Rayl: for gamma SubType= 11
- Lambda table from 10 eV to 100 keV in 80 bins, spline: 0
- LambdaPrime table from 100 keV to 10 TeV in 160 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- LivermoreRayleigh : Emin= 0 eV Emax= 10 TeV CullenGenerator
-
-msc: for e- SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-### === Deexcitation model UAtomDeexcitation is activated for 1 region:
- DefaultRegionForTheWorld
-
-eIoni: for e- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-msc: for e+ SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-eIoni: for e+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-annihil: for e+ SubType= 5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eplus2gg : Emin= 0 eV Emax= 10 TeV
-
-msc: for proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for GenericIon SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV
-
-ionIoni: for GenericIon SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- Stopping Power data for 17 ion/material pairs
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ParamICRU73 : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for GenericIon SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for alpha SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-ionIoni: for alpha SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BraggIon : Emin= 0 eV Emax= 7.9452 MeV
- BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV
-
-nuclearStopping: for alpha SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for anti_proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for anti_proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for anti_proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for anti_proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for anti_proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for kaon+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for kaon- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu+ SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu+ SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu- SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu- SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
- This region is in the mass world.
- Root logical volume(s) : World
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials : Galactic Silicon
- Production cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
-
-Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
- Root logical volume(s) :
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials :
- Production cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
-
-========= Table of registered couples ==============================
-
-Index : 0 used in the geometry : Yes recalculation needed : No
- Material : Galactic
- Range cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-Index : 1 used in the geometry : Yes recalculation needed : No
- Material : Silicon
- Range cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 keV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-====================================================================
-
-Start closing geometry.
-G4GeometryManager::ReportVoxelStats -- Voxel Statistics
-
- Total memory consumed for geometry optimisation: 0 kByte
- Total CPU time elapsed for geometry optimisation: 0 seconds
-### Run 0 start.
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 9876, 54321
-----------------------------------------
-Start Run processing.
-
----> Begin of Event: 0
-
----> Begin of Event: 10000
-
----> Begin of Event: 20000
-
----> Begin of Event: 30000
-
----> Begin of Event: 40000
-
----> Begin of Event: 50000
-
----> Begin of Event: 60000
-
----> Begin of Event: 70000
-
----> Begin of Event: 80000
-
----> Begin of Event: 90000
-Run terminated.
-Run Summary
- Number of events processed : 100000
- User=2.6s Real=6.47s Sys=0.06s
-
- ======================== run summary ======================
-
- The run was 100000 pi+ of 5 GeV through 20 um of Silicon (density: 2.33 g/cm3 )
-
- Total energy deposit in absorber per event = 6.321 keV +- 16.01 eV
-
- -----> Mean dE/dx = 3.16 MeV/cm (1.356 MeV*cm2/g)
-
- From formulas :
- restricted dEdx = 2.963 MeV/cm (1.272 MeV*cm2/g)
- full dEdx = 4.725 MeV/cm (2.028 MeV*cm2/g)
-
- Leakage : primary = 5 GeV +- 303.7 eV secondaries = 2.098 keV +- 299.2 eV
- Energy balance : edep + eleak = 5 GeV
-
- Total track length (charged) in absorber per event = 20.25 um +- 7.542 nm
- Total track length (neutral) in absorber per event = 4.992 Ang +- 3.311 Ang
-
- Number of steps (charged) in absorber per event = 1.584 +- 8.684e-06
- Number of steps (neutral) in absorber per event = 4e-05 +- 8.684e-06
-
- Number of secondaries per event : Gammas = 4e-05 electrons = 0.0164 positrons = 1e-05
-
- Number of events with the primary particle transmitted = 100 %
- Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
-
- Number of events with the primary particle reflected = 0 %
- Number of events with at least 1 particle reflected (same charge as primary) = 0 %
-
- MultipleScattering:
- rms proj angle of transmit primary particle = 0.02527 mrad (central part only)
- computed theta0 (Highland formula) = 0.02627 mrad
- central part defined as +- 0.0788 mrad Tail ratio = 2.96 %
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 1107637770, 1688293817
-----------------------------------------
-#
-/testem/det/setAbsMat Gold
-/testem/det/setAbsThick 9.658 um
-/testem/det/setAbsYZ 1 mm
-/testem/det/update
-
- Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
- ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: 1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: 2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- Material: Gold density: 19.320 g/cm3 RadL: 3.344 mm Nucl.Int.Length: 10.539 cm Imean: 790.000 eV
- ---> Element: Gold ( ) Z = 79.0 N = 197.0 A = 196.97 g/mole
- ---> Isotope: 197 Z = 79 N = 197 A = 196.97 g/mole abundance: 100.00 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- The WORLD is made of 14.487 um of Galactic. The transverse size (YZ) of the world is 1.2 mm
- The ABSORBER is made of 9.658 um of Gold. The transverse size (YZ) is 1 mm
- X position of the middle of the absorber 0 fm
-Obsolete world logical volume is removed from the default region.
-World is registered to the default region.
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 15.7 MeV
-/run/beamOn 50000
-
-phot: for gamma SubType= 12
- LambdaPrime table from 200 keV to 10 TeV in 154 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila FluoActive
-
-compt: for gamma SubType= 13
- Lambda table from 10 eV to 1 MeV in 100 bins, spline: 1
- LambdaPrime table from 1 MeV to 10 TeV in 140 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- KleinNishina : Emin= 0 eV Emax= 10 TeV FluoActive
-
-conv: for gamma SubType= 14
- Lambda table from 1.022 MeV to 10 TeV in 140 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BetheHeitler : Emin= 0 eV Emax= 80 GeV
- BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV
-
-Rayl: for gamma SubType= 11
- Lambda table from 10 eV to 100 keV in 80 bins, spline: 0
- LambdaPrime table from 100 keV to 10 TeV in 160 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- LivermoreRayleigh : Emin= 0 eV Emax= 10 TeV CullenGenerator
-
-msc: for e- SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-### === Deexcitation model UAtomDeexcitation is activated for 1 region:
- DefaultRegionForTheWorld
-
-eIoni: for e- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-msc: for e+ SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-eIoni: for e+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-annihil: for e+ SubType= 5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eplus2gg : Emin= 0 eV Emax= 10 TeV
-
-msc: for proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for GenericIon SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV
-
-ionIoni: for GenericIon SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- Stopping Power data for 17 ion/material pairs
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ParamICRU73 : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for GenericIon SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for alpha SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-ionIoni: for alpha SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BraggIon : Emin= 0 eV Emax= 7.9452 MeV
- BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV
-
-nuclearStopping: for alpha SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for anti_proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for anti_proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for anti_proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for anti_proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for anti_proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for kaon+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for kaon- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu+ SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu+ SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu- SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu- SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
- This region is in the mass world.
- Root logical volume(s) : World
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials : Galactic Gold
- Production cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
-
-Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
- Root logical volume(s) :
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials :
- Production cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
-
-========= Table of registered couples ==============================
-
-Index : 0 used in the geometry : Yes recalculation needed : No
- Material : Galactic
- Range cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-Index : 1 used in the geometry : No recalculation needed : No
- Material : Silicon
- Range cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 21.4886 keV e+ 20.8954 keV proton 100 keV
-
-Index : 2 used in the geometry : Yes recalculation needed : No
- Material : Gold
- Range cuts : gamma 7 um e- 7 um e+ 7 um proton 1 mm
- Energy thresholds : gamma 6.54229 keV e- 65.1406 keV e+ 63.8405 keV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-====================================================================
-
-Start closing geometry.
-G4GeometryManager::ReportVoxelStats -- Voxel Statistics
-
- Total memory consumed for geometry optimisation: 0 kByte
- Total CPU time elapsed for geometry optimisation: 0 seconds
-### Run 1 start.
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 1107637770, 1688293817
-----------------------------------------
-Start Run processing.
-
----> Begin of Event: 0
-
----> Begin of Event: 10000
-
----> Begin of Event: 20000
-
----> Begin of Event: 30000
-
----> Begin of Event: 40000
-Run terminated.
-Run Summary
- Number of events processed : 50000
- User=2.64s Real=9.21s Sys=0.11s
-
- ======================== run summary ======================
-
- The run was 50000 e- of 15.7 MeV through 9.66 um of Gold (density: 19.3 g/cm3 )
-
- Total energy deposit in absorber per event = 18.41 keV +- 60.7 eV
-
- -----> Mean dE/dx = 19.06 MeV/cm (0.9865 MeV*cm2/g)
-
- From formulas :
- restricted dEdx = 17.91 MeV/cm (0.9269 MeV*cm2/g)
- full dEdx = 61.41 MeV/cm (3.179 MeV*cm2/g)
-
- Leakage : primary = 15.64 MeV +- 2.102 keV secondaries = 39.42 keV +- 2.099 keV
- Energy balance : edep + eleak = 15.7 MeV
-
- Total track length (charged) in absorber per event = 9.905 um +- 10.26 nm
- Total track length (neutral) in absorber per event = 120.7 nm +- 4.013 nm
-
- Number of steps (charged) in absorber per event = 6.69 +- 1.433e-05
- Number of steps (neutral) in absorber per event = 0.02738 +- 1.433e-05
-
- Number of secondaries per event : Gammas = 0.02728 electrons = 0.02166 positrons = 0
-
- Number of events with the primary particle transmitted = 99.98 %
- Number of events with at least 1 particle transmitted (same charge as primary) = 99.98 %
-
- Number of events with the primary particle reflected = 0.022 %
- Number of events with at least 1 particle reflected (same charge as primary) = 0.332 %
-
- MultipleScattering:
- rms proj angle of transmit primary particle = 34.66 mrad (central part only)
- computed theta0 (Highland formula) = 35.1 mrad
- central part defined as +- 105.3 mrad Tail ratio = 2.608 %
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 500007974, 1226617912
-----------------------------------------
-#
-/testem/phys/setCuts 10 um
-#
-/testem/det/setAbsMat G4_POLYETHYLENE
-/testem/det/setAbsThick 4.74 mm
-/testem/det/setAbsYZ 50 mm
-/testem/det/update
-
- Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
- ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: 1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: 2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- Material: G4_POLYETHYLENE (C_2H_4)_N-Polyethylene density: 940.000 mg/cm3 RadL: 47.632 cm Nucl.Int.Length: 72.122 cm Imean: 57.400 eV
- ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 14.37 % ElmAbundance 66.67 %
-
- ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole
- ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 %
- ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 %
- ElmMassFraction: 85.63 % ElmAbundance 33.33 %
-
-
- The WORLD is made of 7.11 mm of Galactic. The transverse size (YZ) of the world is 6 cm
- The ABSORBER is made of 4.74 mm of G4_POLYETHYLENE. The transverse size (YZ) is 5 cm
- X position of the middle of the absorber 0 fm
-Obsolete world logical volume is removed from the default region.
-World is registered to the default region.
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 96.2 MeV
-/run/beamOn 50000
-
-phot: for gamma SubType= 12
- LambdaPrime table from 200 keV to 10 TeV in 154 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila FluoActive
-
-compt: for gamma SubType= 13
- Lambda table from 10 eV to 1 MeV in 100 bins, spline: 1
- LambdaPrime table from 1 MeV to 10 TeV in 140 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- KleinNishina : Emin= 0 eV Emax= 10 TeV FluoActive
-
-conv: for gamma SubType= 14
- Lambda table from 1.022 MeV to 10 TeV in 140 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BetheHeitler : Emin= 0 eV Emax= 80 GeV
- BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV
-
-Rayl: for gamma SubType= 11
- Lambda table from 10 eV to 100 keV in 80 bins, spline: 0
- LambdaPrime table from 100 keV to 10 TeV in 160 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- LivermoreRayleigh : Emin= 0 eV Emax= 10 TeV CullenGenerator
-
-msc: for e- SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-### === Deexcitation model UAtomDeexcitation is activated for 1 region:
- DefaultRegionForTheWorld
-
-eIoni: for e- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-msc: for e+ SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-eIoni: for e+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-annihil: for e+ SubType= 5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eplus2gg : Emin= 0 eV Emax= 10 TeV
-
-msc: for proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for GenericIon SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV
-
-ionIoni: for GenericIon SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- Stopping Power data for 17 ion/material pairs
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ParamICRU73 : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for GenericIon SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for alpha SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-ionIoni: for alpha SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BraggIon : Emin= 0 eV Emax= 7.9452 MeV
- BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV
-
-nuclearStopping: for alpha SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for anti_proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for anti_proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for anti_proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for anti_proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for anti_proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for kaon+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for kaon- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu+ SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu+ SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu- SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu- SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
- This region is in the mass world.
- Root logical volume(s) : World
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials : Galactic G4_POLYETHYLENE
- Production cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
-
-Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
- Root logical volume(s) :
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials :
- Production cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
-
-========= Table of registered couples ==============================
-
-Index : 0 used in the geometry : Yes recalculation needed : No
- Material : Galactic
- Range cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-Index : 1 used in the geometry : No recalculation needed : No
- Material : Silicon
- Range cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 31.9516 keV e+ 31.7376 keV proton 100 keV
-
-Index : 2 used in the geometry : No recalculation needed : No
- Material : Gold
- Range cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
- Energy thresholds : gamma 8.27707 keV e- 80.7687 keV e+ 79.1567 keV proton 100 keV
-
-Index : 3 used in the geometry : Yes recalculation needed : No
- Material : G4_POLYETHYLENE
- Range cuts : gamma 10 um e- 10 um e+ 10 um proton 1 mm
- Energy thresholds : gamma 990 eV e- 14.1748 keV e+ 13.7548 keV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-====================================================================
-
-Start closing geometry.
-G4GeometryManager::ReportVoxelStats -- Voxel Statistics
-
- Total memory consumed for geometry optimisation: 0 kByte
- Total CPU time elapsed for geometry optimisation: 0 seconds
-### Run 2 start.
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 500007974, 1226617912
-----------------------------------------
-Start Run processing.
-
----> Begin of Event: 0
-
----> Begin of Event: 10000
-
----> Begin of Event: 20000
-
----> Begin of Event: 30000
-
----> Begin of Event: 40000
-Run terminated.
-Run Summary
- Number of events processed : 50000
- User=5.83s Real=15.19s Sys=0.22s
-
- ======================== run summary ======================
-
- The run was 50000 mu+ of 96.2 MeV through 4.74 mm of G4_POLYETHYLENE (density: 940 mg/cm3)
-
- Total energy deposit in absorber per event = 1.05 MeV +- 1.005 keV
-
- -----> Mean dE/dx = 2.216 MeV/cm (2.357 MeV*cm2/g)
-
- From formulas :
- restricted dEdx = 1.754 MeV/cm (1.866 MeV*cm2/g)
- full dEdx = 2.266 MeV/cm (2.41 MeV*cm2/g)
-
- Leakage : primary = 95.12 MeV +- 1.344 keV secondaries = 25.79 keV +- 696.5 eV
- Energy balance : edep + eleak = 96.2 MeV
-
- Total track length (charged) in absorber per event = 5.176 mm +- 3.805 um
- Total track length (neutral) in absorber per event = 23.46 um +- 2.722 um
-
- Number of steps (charged) in absorber per event = 16.98 +- 0.0002759
- Number of steps (neutral) in absorber per event = 0.00896 +- 0.0002759
-
- Number of secondaries per event : Gammas = 0.00856 electrons = 4.057 positrons = 0
-
- Number of events with the primary particle transmitted = 100 %
- Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
-
- Number of events with the primary particle reflected = 0 %
- Number of events with at least 1 particle reflected (same charge as primary) = 0 %
-
- MultipleScattering:
- rms proj angle of transmit primary particle = 7.31 mrad (central part only)
- computed theta0 (Highland formula) = 7.636 mrad
- central part defined as +- 22.91 mrad Tail ratio = 1.256 %
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 1071691336, 1171911329
-----------------------------------------
-#
-/testem/phys/setCuts 1 mm
-#
-/testem/stack/killSecondaries
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 8.004 mm
-/testem/det/update
-
- Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
- ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: 1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: 2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.860 cm Imean: 166.000 eV
- ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole
- ---> Isotope: 27 Z = 13 N = 27 A = 26.98 g/mole abundance: 100.00 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- The WORLD is made of 1.2006 cm of Galactic. The transverse size (YZ) of the world is 6 cm
- The ABSORBER is made of 8.004 mm of Aluminium. The transverse size (YZ) is 5 cm
- X position of the middle of the absorber 0 fm
-Obsolete world logical volume is removed from the default region.
-World is registered to the default region.
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 174.664 GeV
-/run/beamOn 50000
-
-phot: for gamma SubType= 12
- LambdaPrime table from 200 keV to 10 TeV in 154 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila FluoActive
-
-compt: for gamma SubType= 13
- Lambda table from 10 eV to 1 MeV in 100 bins, spline: 1
- LambdaPrime table from 1 MeV to 10 TeV in 140 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- KleinNishina : Emin= 0 eV Emax= 10 TeV FluoActive
-
-conv: for gamma SubType= 14
- Lambda table from 1.022 MeV to 10 TeV in 140 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BetheHeitler : Emin= 0 eV Emax= 80 GeV
- BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV
-
-Rayl: for gamma SubType= 11
- Lambda table from 10 eV to 100 keV in 80 bins, spline: 0
- LambdaPrime table from 100 keV to 10 TeV in 160 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- LivermoreRayleigh : Emin= 0 eV Emax= 10 TeV CullenGenerator
-
-msc: for e- SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-### === Deexcitation model UAtomDeexcitation is activated for 1 region:
- DefaultRegionForTheWorld
-
-eIoni: for e- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-msc: for e+ SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-eIoni: for e+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-annihil: for e+ SubType= 5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eplus2gg : Emin= 0 eV Emax= 10 TeV
-
-msc: for proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for GenericIon SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV
-
-ionIoni: for GenericIon SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- Stopping Power data for 17 ion/material pairs
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ParamICRU73 : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for GenericIon SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for alpha SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-ionIoni: for alpha SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BraggIon : Emin= 0 eV Emax= 7.9452 MeV
- BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV
-
-nuclearStopping: for alpha SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for anti_proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for anti_proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for anti_proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for anti_proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for anti_proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for kaon+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for kaon- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu+ SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu+ SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu- SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu- SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
- This region is in the mass world.
- Root logical volume(s) : World
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials : Galactic Aluminium
- Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
-
-Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
- Root logical volume(s) :
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials :
- Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
-
-========= Table of registered couples ==============================
-
-Index : 0 used in the geometry : Yes recalculation needed : No
- Material : Galactic
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-Index : 1 used in the geometry : No recalculation needed : No
- Material : Silicon
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 100 keV
-
-Index : 2 used in the geometry : No recalculation needed : No
- Material : Gold
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 100 keV
-
-Index : 3 used in the geometry : No recalculation needed : No
- Material : G4_POLYETHYLENE
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 100 keV
-
-Index : 4 used in the geometry : Yes recalculation needed : No
- Material : Aluminium
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-====================================================================
-
-Start closing geometry.
-G4GeometryManager::ReportVoxelStats -- Voxel Statistics
-
- Total memory consumed for geometry optimisation: 0 kByte
- Total CPU time elapsed for geometry optimisation: 0 seconds
-### Run 3 start.
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 1071691336, 1171911329
-----------------------------------------
-Start Run processing.
-
----> Begin of Event: 0
-
----> Begin of Event: 10000
-
----> Begin of Event: 20000
-
----> Begin of Event: 30000
-
----> Begin of Event: 40000
-Run terminated.
-Run Summary
- Number of events processed : 50000
- User=1.27s Real=2.87s Sys=0s
-
- ======================== run summary ======================
-
- The run was 50000 proton of 175 GeV through 8 mm of Aluminium (density: 2.7 g/cm3 )
-
- Total energy deposit in absorber per event = 4.794 MeV +- 188 keV
-
- -----> Mean dE/dx = 5.99 MeV/cm (2.219 MeV*cm2/g)
-
- From formulas :
- restricted dEdx = 4.024 MeV/cm (1.49 MeV*cm2/g)
- full dEdx = 6.049 MeV/cm (2.24 MeV*cm2/g)
-
- Leakage : primary = 174.7 GeV +- 188 keV secondaries = 0 eV +- 0 eV
- Energy balance : edep + eleak = 174.7 GeV
-
- Total track length (charged) in absorber per event = 8.004 mm +- 0 fm
- Total track length (neutral) in absorber per event = 0 fm +- 0 fm
-
- Number of steps (charged) in absorber per event = 1.269 +- 0
- Number of steps (neutral) in absorber per event = 0 +- 0
-
- Number of secondaries per event : Gammas = 0 electrons = 0.269 positrons = 0.00078
-
- Number of events with the primary particle transmitted = 100 %
- Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
-
- Number of events with the primary particle reflected = 0 %
- Number of events with at least 1 particle reflected (same charge as primary) = 0 %
-
- MultipleScattering:
- rms proj angle of transmit primary particle = 0.02072 mrad (central part only)
- computed theta0 (Highland formula) = 0.02111 mrad
- central part defined as +- 0.06333 mrad Tail ratio = 1.915 %
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 855224844, 1534153149
-----------------------------------------
-#
-/testem/det/setAbsMat Iron
-/testem/det/setAbsThick 1 m
-/testem/det/update
-
- Material: Galactic density: 0.000 kg/m3 RadL: 204727512.315 pc Nucl.Int.Length: 114561548.020 pc Imean: 19.200 eV temperature: 2.73 K pressure: 0.00 atm
- ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole
- ---> Isotope: 1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %
- ---> Isotope: 2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- Material: Iron density: 7.870 g/cm3 RadL: 1.759 cm Nucl.Int.Length: 16.969 cm Imean: 286.000 eV
- ---> Element: Iron ( ) Z = 26.0 N = 55.8 A = 55.85 g/mole
- ---> Isotope: 54 Z = 26 N = 54 A = 53.94 g/mole abundance: 5.84 %
- ---> Isotope: 56 Z = 26 N = 56 A = 55.93 g/mole abundance: 91.75 %
- ---> Isotope: 57 Z = 26 N = 57 A = 56.94 g/mole abundance: 2.12 %
- ---> Isotope: 58 Z = 26 N = 58 A = 57.93 g/mole abundance: 0.28 %
- ElmMassFraction: 100.00 % ElmAbundance 100.00 %
-
-
- The WORLD is made of 1.5 m of Galactic. The transverse size (YZ) of the world is 6 cm
- The ABSORBER is made of 1 m of Iron. The transverse size (YZ) is 5 cm
- X position of the middle of the absorber 0 fm
-Obsolete world logical volume is removed from the default region.
-World is registered to the default region.
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 100 GeV
-/run/beamOn 50000
-
-phot: for gamma SubType= 12
- LambdaPrime table from 200 keV to 10 TeV in 154 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- PhotoElectric : Emin= 0 eV Emax= 10 TeV AngularGenSauterGavrila FluoActive
-
-compt: for gamma SubType= 13
- Lambda table from 10 eV to 1 MeV in 100 bins, spline: 1
- LambdaPrime table from 1 MeV to 10 TeV in 140 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- KleinNishina : Emin= 0 eV Emax= 10 TeV FluoActive
-
-conv: for gamma SubType= 14
- Lambda table from 1.022 MeV to 10 TeV in 140 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BetheHeitler : Emin= 0 eV Emax= 80 GeV
- BetheHeitlerLPM : Emin= 80 GeV Emax= 10 TeV
-
-Rayl: for gamma SubType= 11
- Lambda table from 10 eV to 100 keV in 80 bins, spline: 0
- LambdaPrime table from 100 keV to 10 TeV in 160 bins
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- LivermoreRayleigh : Emin= 0 eV Emax= 10 TeV CullenGenerator
-
-msc: for e- SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-### === Deexcitation model UAtomDeexcitation is activated for 1 region:
- DefaultRegionForTheWorld
-
-eIoni: for e- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-msc: for e+ SubType= 10
- RangeFactor= 0.01, stepLimitType: 2, latDisplacement: 1, skin= 1, geomFactor= 2.5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-eIoni: for e+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MollerBhabha : Emin= 0 eV Emax= 10 TeV
-
-eBrem: for e+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- LPM flag: 1 for E > 1 GeV
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eBremSB : Emin= 0 eV Emax= 1 GeV DipBustGen
- eBremLPM : Emin= 1 GeV Emax= 10 TeV DipBustGen
-
-annihil: for e+ SubType= 5
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eplus2gg : Emin= 0 eV Emax= 10 TeV
-
-msc: for proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for GenericIon SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV
-
-ionIoni: for GenericIon SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.001, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- Stopping Power data for 17 ion/material pairs
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ParamICRU73 : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for GenericIon SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for alpha SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 0
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-ionIoni: for alpha SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.01, dRoverRange= 0.1, integral: 1, fluct: 1, linLossLimit= 0.02
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- BraggIon : Emin= 0 eV Emax= 7.9452 MeV
- BetheBloch : Emin= 7.9452 MeV Emax= 10 TeV
-
-nuclearStopping: for alpha SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 10 TeV
-
-msc: for anti_proton SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for anti_proton SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 2 MeV
- BetheBloch : Emin= 2 MeV Emax= 10 TeV
-
-hBrems: for anti_proton SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for anti_proton SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-nuclearStopping: for anti_proton SubType= 8
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU49NucStopping : Emin= 0 eV Emax= 1 MeV
-
-msc: for kaon+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for kaon- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for kaon- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 1.05231 MeV
- BetheBloch : Emin= 1.05231 MeV Emax= 10 TeV
-
-hBrems: for kaon- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for kaon- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu+ SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu+ SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for mu- SubType= 10
- RangeFactor= 0.2, step limit type: 0, lateralDisplacement: 1, polarAngleLimit(deg)= 180
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- WentzelVIUni : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-muIoni: for mu- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 200 keV
- BetheBloch : Emin= 200 keV Emax= 1 GeV
- MuBetheBloch : Emin= 1 GeV Emax= 10 TeV
-
-muBrems: for mu- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- MuBrem : Emin= 0 eV Emax= 10 TeV
-
-muPairProd: for mu- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- muPairProd : Emin= 0 eV Emax= 10 TeV
-
-CoulombScat: for mu- SubType= 1
- Lambda table from 10 eV to 10 TeV in 108 bins, spline: 1
- 180 < Theta(degree) < 180 pLimit(GeV^1)= 0.139531
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- eCoulombScattering : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi+ SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi+ SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- Bragg : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi+ SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi+ SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-msc: for pi- SubType= 10
- RangeFactor= 0.2, stepLimitType: 0, latDisplacement: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- UrbanMsc95 : Emin= 0 eV Emax= 10 TeV Table with 220 bins Emin= 100 eV Emax= 10 TeV
-
-hIoni: for pi- SubType= 2
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- finalRange(mm)= 0.05, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- ICRU73QO : Emin= 0 eV Emax= 297.505 keV
- BetheBloch : Emin= 297.505 keV Emax= 10 TeV
-
-hBrems: for pi- SubType= 3
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hBrem : Emin= 0 eV Emax= 10 TeV
-
-hPairProd: for pi- SubType= 4
- dE/dx and range tables from 10 eV to 10 TeV in 240 bins
- Lambda tables from threshold to 10 TeV in 240 bins, spline: 1
- ===== EM models for the G4Region DefaultRegionForTheWorld ======
- hPairProd : Emin= 0 eV Emax= 10 TeV
-
-Region <DefaultRegionForTheWorld> -- -- appears in <World> world volume
- This region is in the mass world.
- Root logical volume(s) : World
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials : Galactic Iron
- Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
-
-Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
- Root logical volume(s) :
- Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
- Materials :
- Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
-
-========= Table of registered couples ==============================
-
-Index : 0 used in the geometry : Yes recalculation needed : No
- Material : Galactic
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-Index : 1 used in the geometry : No recalculation needed : No
- Material : Silicon
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 6.95018 keV e- 548.291 keV e+ 526.624 keV proton 100 keV
-
-Index : 2 used in the geometry : No recalculation needed : No
- Material : Gold
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 112.644 keV e- 2.30976 MeV e+ 2.13081 MeV proton 100 keV
-
-Index : 3 used in the geometry : No recalculation needed : No
- Material : G4_POLYETHYLENE
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 2.26261 keV e- 351.877 keV e+ 342.545 keV proton 100 keV
-
-Index : 4 used in the geometry : No recalculation needed : No
- Material : Aluminium
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 6.90363 keV e- 598.345 keV e+ 570.85 keV proton 100 keV
-
-Index : 5 used in the geometry : Yes recalculation needed : No
- Material : Iron
- Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 1 mm
- Energy thresholds : gamma 20.6438 keV e- 1.29592 MeV e+ 1.21169 MeV proton 100 keV
- Region(s) which use this couple :
- DefaultRegionForTheWorld
-
-====================================================================
-
-Start closing geometry.
-G4GeometryManager::ReportVoxelStats -- Voxel Statistics
-
- Total memory consumed for geometry optimisation: 0 kByte
- Total CPU time elapsed for geometry optimisation: 0 seconds
-### Run 4 start.
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 855224844, 1534153149
-----------------------------------------
-Start Run processing.
-
----> Begin of Event: 0
-
----> Begin of Event: 10000
-
----> Begin of Event: 20000
-
----> Begin of Event: 30000
-
----> Begin of Event: 40000
-Run terminated.
-Run Summary
- Number of events processed : 50000
- User=16.36s Real=50.59s Sys=0.01s
-
- ======================== run summary ======================
-
- The run was 50000 mu+ of 100 GeV through 1 m of Iron (density: 7.87 g/cm3 )
-
- Total energy deposit in absorber per event = 2.153 GeV +- 15.31 MeV
-
- -----> Mean dE/dx = 21.53 MeV/cm (2.735 MeV*cm2/g)
-
- From formulas :
- restricted dEdx = 11.22 MeV/cm (1.425 MeV*cm2/g)
- full dEdx = 21.6 MeV/cm (2.745 MeV*cm2/g)
-
- Leakage : primary = 97.85 GeV +- 15.31 MeV secondaries = 0 eV +- 0 eV
- Energy balance : edep + eleak = 100 GeV
-
- Total track length (charged) in absorber per event = 100 cm +- 6.052 um
- Total track length (neutral) in absorber per event = 0 fm +- 0 fm
-
- Number of steps (charged) in absorber per event = 47.93 +- 0.0003479
- Number of steps (neutral) in absorber per event = 0 +- 0
-
- Number of secondaries per event : Gammas = 0.04914 electrons = 46.33 positrons = 1.256
-
- Number of events with the primary particle transmitted = 100 %
- Number of events with at least 1 particle transmitted (same charge as primary) = 100 %
-
- Number of events with the primary particle reflected = 0 %
- Number of events with at least 1 particle reflected (same charge as primary) = 0 %
-
- MultipleScattering:
- rms proj angle of transmit primary particle = 1.127 mrad (central part only)
- computed theta0 (Highland formula) = 1.182 mrad
- central part defined as +- 3.545 mrad Tail ratio = 0.312 %
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 2037953065, 756024662
-----------------------------------------
-#
-G4 kernel has come to Quit state.
-UserDetectorConstruction deleted.
-UserPhysicsList deleted.
-UserRunAction deleted.
-UserPrimaryGenerator deleted.
-EventManager deleted.
-UImanager deleted.
-Units table cleared.
-StateManager deleted.
-RunManagerKernel is deleted.
-RunManager is deleted.
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/_CMakeLists.txt b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/_CMakeLists.txt
deleted file mode 100644
index a22b19a..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/_CMakeLists.txt
+++ /dev/null
@@ -1,29 +0,0 @@
-# $Id: GNUmakefile,v 1.15 2008-06-11 22:15:20 maire Exp $
-# --------------------------------------------------------------
-# GNUmakefile for examples module. Gabriele Cosmo, 06/04/98.
-# --------------------------------------------------------------
-
-name := TestEm5
-G4TARGET := $(name)
-G4EXLIB := true
-
-ifndef G4INSTALL
- G4INSTALL = ../../../..
-endif
-
-.PHONY: all
-all: hbook lib bin
-
-#### G4_USE_HBOOK := true
-include GNUmakefile.tools_hbook
-
-include $(G4INSTALL)/config/architecture.gmk
-
-include $(G4INSTALL)/config/binmake.gmk
-
-histclean:
- rm -f $(G4WORKDIR)/tmp/$(G4SYSTEM)/$(G4TARGET)/HistoManager.o
-
-visclean:
- rm -f g4*.prim g4*.eps g4*.wrl
- rm -f .DAWN_*
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/acosta.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/acosta.mac
deleted file mode 100644
index 8bcb9f8..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/acosta.mac
+++ /dev/null
@@ -1,29 +0,0 @@
-# $Id: acosta.mac,v 1.2 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# Back x-ray emission by 20 keV electrons in Silver.
-# E. Acosta et al. Journal of Applied Physics 83(11) 1998 page 6038
-# Fig. 4-5-6
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silver
-/testem/det/setAbsThick 1 cm
-/testem/det/setAbsYZ 1 cm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 20 keV
-#
-/analysis/setFileName acosta
-/analysis/h1/set 40 40 0 20 keV #energy at exit. (reflect,neutral)
-#
-/testem/phys/setGCut 1 um
-#
-/testem/event/printModulo 100000
-#
-/run/beamOn 5000000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/anthony.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/anthony.mac
deleted file mode 100644
index 90f61a8..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/anthony.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: anthony.mac,v 1.4 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# LPM and dielectric effect measurement:
-# 25 GeV electrons through thin foils.
-# P.L. Anthony et al. Phys.Rev. D 56 (1997) page 1373.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 5.30 mm
-###/testem/det/setAbsMat Lead
-###/testem/det/setAbsThick 0.15 mm
-/testem/det/setAbsYZ 10 cm
-#
-/testem/phys/setGCut 1 mm
-/testem/phys/setECut 5 um
-/testem/phys/setPCut 5 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 25 GeV
-#
-/analysis/setFileName anthony
-/analysis/h1/set 20 80 1.e-4 1. GeV #gamma energy at exit
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/berger.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/berger.mac
deleted file mode 100644
index 390408a..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/berger.mac
+++ /dev/null
@@ -1,36 +0,0 @@
-# $Id: berger.mac,v 1.8 2010-05-19 09:35:37 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# energy deposit by 1 MeV electrons in silicon counters.
-# M.J.Berger et al. NIM 69 (1969) page 181
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silicon
-###/testem/det/setAbsThick 61 um
-###/testem/det/setAbsThick 105 um
-###/testem/det/setAbsThick 191 um
-/testem/det/setAbsThick 530 um
-###/testem/det/setAbsThick 1 mm
-###/testem/det/setAbsThick 3 mm
-/testem/det/setAbsYZ 1 cm
-#
-###/testem/phys/addPhysics emstandard_opt0
-/testem/phys/addPhysics local
-#
-/testem/phys/setCuts 10 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 1 MeV
-#
-/analysis/setFileName berger/local
-/analysis/h1/set 1 105 0 1.05 MeV #energy deposit
-#
-/testem/event/printModulo 100000
-#
-/run/beamOn 1000000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/bichsel.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/bichsel.mac
deleted file mode 100644
index d23e38e..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/bichsel.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: bichsel.mac,v 1.4 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# 0.766 MeV protons, transmitted through 1.37 mg/cm2 Al
-# H.Bichsel Phys.Rev. 112 (1958) page 182.
-#
-# proj. angle distr. exp(-(th/th0)**2) th0 = 4.44+- 0.20 deg
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 5.0741 um
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setCuts 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 766 keV
-#
-/analysis/setFileName bichsel
-/analysis/h1/set 12 100 0 9 deg #space angle
-/analysis/h1/set 13 100 -9 +9 deg #projected angle
-#
-/testem/stack/killSecondaries
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx1.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx1.mac
deleted file mode 100644
index a4547bd..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx1.mac
+++ /dev/null
@@ -1,29 +0,0 @@
-# $Id: dedx1.mac,v 1.6 2009-09-18 14:03:24 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# to control dE/dx calculation
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Air
-/testem/det/setAbsThick 1 cm
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/addPhysics local
-#
-/testem/phys/setCuts 1 km
-#
-/run/initialize
-#
-/process/inactivate msc
-/process/eLoss/fluct false
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 10 MeV
-#
-/tracking/verbose 1
-#
-/run/beamOn 5
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx2.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx2.mac
deleted file mode 100644
index 813f516..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/dedx2.mac
+++ /dev/null
@@ -1,29 +0,0 @@
-# $Id: dedx2.mac,v 1.6 2009-09-18 14:03:24 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# to control dE/dx calculation
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Graphite
-/testem/det/setAbsThick 10 um
-/testem/det/setAbsYZ 10 um
-#
-/testem/phys/addPhysics local
-###/testem/phys/addPhysics emlivermore
-###/testem/phys/addPhysics empenelope
-#
-/testem/phys/setCuts 1 km
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 1 MeV
-#
-/analysis/setFileName dedx2
-/analysis/h1/set 1 100 0. 10. keV #energy in absorber
-#
-/run/beamOn 10000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/fluo.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/fluo.mac
deleted file mode 100644
index 719e812..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/fluo.mac
+++ /dev/null
@@ -1,35 +0,0 @@
-# macro file for TestEm5.cc
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Gold
-/testem/det/setAbsThick 50 um
-/testem/det/setAbsYZ 50 um
-#
-/testem/phys/addPhysics emlivermore
-#
-/testem/phys/setGCut 10 nm
-/testem/phys/setECut 10 nm
-/testem/phys/setPCut 10 nm
-#
-/cuts/setLowEdge 250 eV
-#
-/run/initialize
-#
-/process/em/deexcitation World true true true
-/process/em/fluo true
-/process/em/auger true
-/process/em/pixe true
-#
-/testem/gun/setDefault
-/gun/particle gamma
-/gun/energy 1 MeV
-#
-/analysis/setFileName fluo
-/analysis/h1/set 2 100 0 1000 keV #e-: energy at vertex
-/analysis/h1/set 3 100 1 1000 keV #gamma: energy at vertex
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gammaSpectrum.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gammaSpectrum.mac
deleted file mode 100644
index 06ced36..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gammaSpectrum.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# macro file for TestEm5.cc
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Tungsten
-/testem/det/setAbsThick 3 mm
-/testem/det/setAbsYZ 10 mm
-#
-/testem/phys/addPhysics local
-###/testem/phys/addPhysics emstandard_opt3
-/testem/phys/setCuts 100 um
-###/run/setCut 100 um
-#
-/run/initialize
-#
-###/process/em/fluo true
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 6 MeV
-#
-/analysis/setFileName gammaSpectrum
-/analysis/h1/set 3 200 0.01 10 MeV #gamma: energy at vertex
-/analysis/h1/set 5 200 0.01 10 MeV log10 #gamma: energy at vertex (log10)
-/analysis/h1/set 20 200 0 6 MeV #gamma: energy at exit
-/analysis/h1/set 40 200 0 6 MeV #gamma: energy at back
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gottsch.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gottsch.mac
deleted file mode 100644
index 27fcab2..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/gottsch.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: gottsch.mac,v 1.3 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# 158.6 MeV protons, transmitted through 0.2160 g/cm2 Al
-# B.Gottschalk et al. NIM B74 (1993) page 467.
-#
-# proj. angle distr. exp(-0.5*(th/th0)**2) th0 = 3.534 +- 0.491 mrad
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 0.8 mm
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setCuts 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 158.6 MeV
-#
-/analysis/setFileName gottsch
-/analysis/h1/set 12 100 0.0 0.6 deg #space angle
-/analysis/h1/set 13 100 -0.6 +0.6 deg #projected angle
-#
-/testem/stack/killSecondaries
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hanson.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hanson.mac
deleted file mode 100644
index c96ad57..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hanson.mac
+++ /dev/null
@@ -1,35 +0,0 @@
-# $Id: hanson.mac,v 1.5 2009-06-19 12:42:04 vnivanch Exp $
-#
-# macro file for TestEm5.cc
-#
-# angle distribution of 15.7 MeV electrons
-# transmitted through thin gold foils.
-# A.O.Hanson et al. Phys.Rev.84 (1951) page 634.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Gold
-###/testem/det/setAbsThick 9.658 um
-/testem/det/setAbsThick 19.296 um
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setGCut 13.0 um
-/testem/phys/setECut 1.12 um
-/testem/phys/setPCut 1.12 um
-#
-###/testem/phys/addPhysics emstandard_opt0
-/testem/phys/addPhysics local
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 15.7 MeV
-#
-/analysis/setFileName hanson/19um.local
-/anlysis/h1/set 12 120 0 30 deg #space angle
-#
-/testem/event/printModulo 100000
-#
-/run/beamOn 1000000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hunger.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hunger.mac
deleted file mode 100644
index c44d2e4..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/hunger.mac
+++ /dev/null
@@ -1,30 +0,0 @@
-# $Id: hunger.mac,v 1.3 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# Back scattering of 41 keV electrons.
-# H.J. Hunger and L. Kuchler Phys. Stat. Sol.(a) 56, K45 (1979)
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silver
-/testem/det/setAbsThick 50 um
-/testem/det/setAbsYZ 60 um
-#
-/run/initialize
-#
-/testem/phys/setCuts 1 mm
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 41 keV
-#
-/analysis/setFileName hunger
-/analysis/h1/set 30 50 0.0 50 keV #energy reflected_charged
-#
-###/testem/stepMax 1 um
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/ion.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/ion.mac
deleted file mode 100644
index 1ff3aad..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/ion.mac
+++ /dev/null
@@ -1,26 +0,0 @@
-# $Id: ion.mac,v 1.1 2009-06-18 19:08:27 vnivanch Exp $
-#
-# macro file for TestEm5.cc
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat G4_Fe
-#
-/testem/det/setAbsThick 0.1 mm
-/testem/det/setAbsYZ 1 m
-#
-/testem/phys/addPhysics emstandard
-#
-/testem/phys/setCuts 1 km
-#
-/run/initialize
-#
-#/testem/gun/setDefault
-/gun/particle ion
-/gun/ion 6 12
-/gun/energy 1200 MeV
-#
-/tracking/verbose 0
-#
-/run/beamOn 1000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/kulchi.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/kulchi.mac
deleted file mode 100644
index 23b8c75..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/kulchi.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: kulchi.mac,v 1.4 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# 2.25 MeV e-, transmitted through 26.60 mg/cm2 Al
-# L.Kulchitsky Phys.Rev. 61 (1941) page 254.
-#
-# Gaussian th1/e = 9.50 +- 0.40 deg
-# proj. distr. fit: exp(-0.5*(th/th0)**2) ====> th0 = th1/e/sqrt(2)
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 98.5185 um
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setCuts 1 mm
-###/run/setCut 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 2.25 MeV
-#
-/analysis/setFileName kulchi
-/analysis/h1/set 12 100 0. 25. deg #space angle
-/analysis/h1/set 13 100 -25. +25. deg #projected angle
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mumsc.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mumsc.mac
deleted file mode 100644
index fd8ab7d..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mumsc.mac
+++ /dev/null
@@ -1,34 +0,0 @@
-# $Id: mumsc.mac,v 1.21 2010-05-19 09:35:37 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# mu+, transmitted through 1 m of iron
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Iron
-/testem/det/setAbsThick 1 m
-/testem/det/setAbsYZ 1 m
-/testem/det/setWorldX 1 m
-/testem/det/setWorldYZ 1 m
-#
-/testem/phys/addPhysics local
-###/testem/phys/addPhysics emstandard_opt0
-#
-/testem/phys/setCuts 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 100 GeV
-#
-/analysis/setFileName mumsc/local
-/analysis/h1/set 14 100 -5.0 +5.0 mm #x-position at exit
-/analysis/h1/set 15 100 0. +10.0 mm #radius at exit
-#
-/testem/stack/killSecondaries
-#
-/testem/event/printModulo 4000
-/run/beamOn 40000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mutev.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mutev.mac
deleted file mode 100644
index 1501893..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/mutev.mac
+++ /dev/null
@@ -1,39 +0,0 @@
-# $Id: mumsc.mac,v 1.21 2010-05-19 09:35:37 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# mu+, transmitted through 1 m of iron
-#
-/control/verbose 2
-/run/verbose 2
-#/tracking/verbose 1
-#
-/testem/det/setAbsMat Iron
-/testem/det/setAbsMat G4_WATER
-/testem/det/setAbsThick 10 m
-/testem/det/setAbsYZ 10 m
-/testem/det/setWorldX 10 m
-/testem/det/setWorldYZ 10 m
-#
-/testem/phys/addPhysics local
-###/testem/phys/addPhysics emstandard_opt0
-#
-/testem/phys/setCuts 1 km
-#
-/run/initialize
-#
-#/process/msc/ThetaLimit 0. rad
-#/process/msc/LateralDisplacement false
-#
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 1000 GeV
-#
-#/analysis/setFileName mumsc/local
-/analysis/h1/set 14 100 -5.0 +5.0 mm #x-position at exit
-/analysis/h1/set 15 100 0. +10.0 mm #radius at exit
-#
-/testem/stack/killSecondaries 2
-#
-/testem/event/printModulo 100
-/run/beamOn 1000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.C b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.C
deleted file mode 100644
index ba2edb9..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.C
+++ /dev/null
@@ -1,23 +0,0 @@
-{
- gROOT->Reset();
-
- // Draw histos filled by Geant4 simulation
- //
- TFile f = TFile("pixe.root");
- TCanvas* c1 = new TCanvas("c1", " ");
- c1->SetLogy(1);
- c1->cd();
- c1->Update();
-
- TH1D* hist1 = (TH1D*)f.Get("3");
- hist1->Draw("HIST");
-
- TH1D* hist2 = (TH1D*)f.Get("5");
- hist2->Draw("HIST");
-
- TH1D* hist3 = (TH1D*)f.Get("20");
- hist3->Draw("HIST");
-
- TH1D* hist4 = (TH1D*)f.Get("40");
- hist4->Draw("HIST");
-}
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.mac
deleted file mode 100644
index 5651eff..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/pixe.mac
+++ /dev/null
@@ -1,40 +0,0 @@
-# macro file for TestEm5.cc
-#
-/control/verbose 2
-/run/verbose 2
-#
-#/testem/det/setAbsMat Gold
-/testem/det/setAbsMat G4_Au
-/testem/det/setAbsThick 50 um
-/testem/det/setAbsYZ 50 um
-#
-#/testem/phys/addPhysics local
-#/testem/phys/addPhysics emstandard_opt0
-/testem/phys/addPhysics emlivermore
-#/testem/phys/addPhysics empenelope
-###/testem/phys/addPhysics emstandard_opt3
-/testem/phys/setGCut 0.5 um
-/testem/phys/setECut 1 mm
-/testem/phys/setPCut 1 mm
-#
-/run/initialize
-#
-#/process/em/deexcitation world true false true
-/process/em/fluo true
-/process/em/pixe true
-/process/em/auger false
-#
-/testem/gun/setDefault
-/gun/particle e-
-#/gun/particle proton
-/gun/energy 20 MeV
-#
-/analysis/setFileName pixe
-/analysis/h1/set 3 100 0 100 keV #gamma: energy at vertex
-/analysis/h1/set 5 100 1 100 keV log10 #gamma: energy at vertex (log scale)
-/analysis/h1/set 20 100 0 100 keV #gamma: energy at exit
-/analysis/h1/set 40 100 0 100 keV #gamma: energy at back
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 10000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.C b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.C
deleted file mode 100644
index 22f82b4..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.C
+++ /dev/null
@@ -1,14 +0,0 @@
-{
- gROOT->Reset();
-
- // Draw histos filled by Geant4 simulation
- //
- TFile f = TFile("kulchi.root");
- TCanvas* c1 = new TCanvas("c1", " ");
-
- TH1D* hist12 = (TH1D*)f.Get("12");
- hist12->Draw("HIST");
-
- TH1D* hist13 = (TH1D*)f.Get("13");
- hist13->Draw("HIST");
-}
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.kumac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.kumac
deleted file mode 100644
index 4174307..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/plotHisto.kumac
+++ /dev/null
@@ -1,14 +0,0 @@
-macro plotHisto.kumac
-*
-h/file 1 testem5.hbook
-*
-h/list
-*
-opt grid
-opt stat
-h/pl 13
-*
-opt logy
-h/pl 13
-*
-return
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen1.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen1.mac
deleted file mode 100644
index 75027ed..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen1.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: shen1.mac,v 1.4 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# angle distribution of high energy (50-200 GeV/c) protons
-# transmitted through different targets.
-# G. Shen et al. Phys.Rev. D20 (1979) page 1584.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 4.005 mm
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setGCut 2.25 mm
-/testem/phys/setECut 6.05 um
-/testem/phys/setPCut 6.05 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 69.068 GeV #momentum 70 GeV/c
-#
-/analysis/setFileName shen1
-/analysis/h1/set 12 100 0.00 0.01 deg #space angle
-#
-/testem/stack/killSecondaries
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen2.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen2.mac
deleted file mode 100644
index 76dea72..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/shen2.mac
+++ /dev/null
@@ -1,32 +0,0 @@
-# $Id: shen2.mac,v 1.3 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# proton 175 GeV/c, transmitted through 8.004 mm Al
-# G. Shen et al. Phys.Rev. D20 (1979) page 1584.
-#
-# data theta0 = (1.93+-0.39)*e-5 rad ???
-# proj. distr. exp(-0.5*(th/th0)**2)
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Aluminium
-/testem/det/setAbsThick 8.004 mm
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setCuts 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 174.664 GeV #momentum 175 GeV/c
-#
-/analysis/setFileName shen2
-/analysis/h1/set 12 100 0.00 0.004 deg #space angle
-/analysis/h1/set 13 100 -0.004 +0.004 deg #projected angle
-#
-/testem/stack/killSecondaries
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tavora.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tavora.mac
deleted file mode 100644
index 49d8b66..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tavora.mac
+++ /dev/null
@@ -1,29 +0,0 @@
-# $Id: tavora.mac,v 1.2 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# Back scattering of 35 keV electrons in Silver.
-# L.M. Tavora et al. J.Phys.D: Appl. Phys. 33 (2000) page 2497
-# Fig. 7
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silver
-/testem/det/setAbsThick 50 um
-/testem/det/setAbsYZ 60 um
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle e-
-/gun/energy 35 keV
-#
-/analysis/setFileName tavora
-/analysis/h1/set 30 40 0.0 40 keV #energy reflected_charged
-#
-/testem/stepMax 50 nm
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tramu.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tramu.mac
deleted file mode 100644
index 6079aef..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/tramu.mac
+++ /dev/null
@@ -1,34 +0,0 @@
-# $Id: tramu.mac,v 1.14 2010-05-19 09:35:37 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# 1 TeV mu+, transmitted through 3 m of iron
-# Rev. of Particle Physics Eur. Phys. Jour. C (2000) page 172.
-# Rev. of Particle Physics Letters B 592 (2004) page 251.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Iron
-/testem/det/setAbsThick 3 m
-/testem/det/setAbsYZ 3 m
-#
-###/testem/phys/addPhysics emstandard_opt0
-/testem/phys/addPhysics local
-#
-/testem/phys/setCuts 1 mm
-#
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle mu+
-/gun/energy 1 TeV
-#
-/analysis/setFileName tramu/local
-/analysis/h1/set 10 100 950 1000 GeV #energy at exit
-#
-/testem/stack/killSecondaries
-#
-/testem/event/printModulo 10000
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vincour.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vincour.mac
deleted file mode 100644
index 49b7183..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vincour.mac
+++ /dev/null
@@ -1,31 +0,0 @@
-# $Id: vincour.mac,v 1.4 2007-07-31 16:51:29 maire Exp $
-#
-# macro file for TestEm5.cc
-#
-# angle distribution of 6.56 MeV protons
-# transmitted through thin silicon targets.
-# J.Vincour,P.Bem NIM 148 (1978) page 396.
-#
-/control/verbose 2
-/run/verbose 2
-#
-/testem/det/setAbsMat Silicon
-/testem/det/setAbsThick 50 um
-###/testem/det/setAbsThick 92.6 um
-/testem/det/setAbsYZ 1 cm
-#
-/testem/phys/setGCut 2.13 mm
-/testem/phys/setECut 7 um
-/testem/phys/setPCut 7 um
-/run/initialize
-#
-/testem/gun/setDefault
-/gun/particle proton
-/gun/energy 6.56 MeV
-#
-/analysis/setFileName vincour
-/analysis/h1/set 12 60 0. 6. deg #space angle
-#
-/testem/stack/killSecondaries
-#
-/run/beamOn 100000
diff --git a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vis.mac b/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vis.mac
deleted file mode 100644
index a6968d9..0000000
--- a/Geant4/G4examples/extended/electromagnetic/G4TestEm5/scripts/vis.mac
+++ /dev/null
@@ -1,77 +0,0 @@
-#
-# Macro file for the initialization phase of "TestEm5.cc"
-# Sets some default verbose
-# and initializes the graphic.
-#
-/control/verbose 2
-/run/verbose 2
-/run/initialize
-#
-# Use this open statement to create an OpenGL view:
-#/vis/open OGL
-#
-# Use this open statement to create a .prim file suitable for
-# viewing in DAWN:
-/vis/open DAWNFILE
-#
-# Use this open statement to create a .heprep file suitable for
-# viewing in HepRApp:
-#/vis/open HepRepFile
-#
-# Use this open statement to create a .wrl file suitable for
-# viewing in a VRML viewer:
-#/vis/open VRML2FILE
-#
-# Disable auto refresh and quieten vis messages whilst scene and
-# trajectories are established:
-/vis/viewer/set/autoRefresh false
-/vis/verbose errors
-#
-# Draw geometry:
-/vis/drawVolume
-#
-# Specify view angle:
-#/vis/viewer/set/viewpointThetaPhi 90. 0.
-#
-# Specify zoom value:
-/vis/viewer/zoom 1.4
-#
-# Specify style (surface or wireframe):
-#/vis/viewer/set/style wireframe
-#
-# Draw coordinate axes:
-#/vis/scene/add/axes 0 0 0 1 m
-#
-# Draw smooth trajectories at end of event, showing trajectory points
-# as markers 2 pixels wide:
-/vis/scene/add/trajectories smooth
-/vis/modeling/trajectories/create/drawByCharge
-/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
-/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 1
-# (if too many tracks cause core dump => /tracking/storeTrajectory 0)
-#
-# Draw hits at end of event:
-#/vis/scene/add/hits
-#
-# To draw only gammas:
-#/vis/filtering/trajectories/create/particleFilter
-#/vis/filtering/trajectories/particleFilter-0/add gamma
-#
-# To invert the above, drawing all particles except gammas,
-# keep the above two lines but also add:
-#/vis/filtering/trajectories/particleFilter-0/invert true
-#
-# Many other options are available with /vis/modeling and /vis/filtering.
-# For example, to select colour by particle ID:
-#/vis/modeling/trajectories/create/drawByParticleID
-#/vis/modeling/trajectories/drawByParticleID-0/set e- blue
-#
-# To superimpose all of the events from a given run:
-/vis/scene/endOfEventAction accumulate
-#
-# Re-establish auto refreshing and verbosity:
-/vis/viewer/set/autoRefresh true
-/vis/verbose warnings
-#
-# For file-based drivers, use this to create an empty detector view:
-#/vis/viewer/flush
diff --git a/cmt/project.cmt b/cmt/project.cmt
index 0ffbddc..e8466e9 100755
--- a/cmt/project.cmt
+++ b/cmt/project.cmt
@@ -1,6 +1,6 @@
project GEANT4
-use GAUDI GAUDI_v26r2
+use GAUDI GAUDI_v26r3
use PARAM
build_strategy with_installarea
diff --git a/toolchain.cmake b/toolchain.cmake
index 5d6a2eb..6d47517 100644
--- a/toolchain.cmake
+++ b/toolchain.cmake
@@ -1,5 +1,5 @@
# Special wrapper to load the declared version of the heptools toolchain.
-set(heptools_version 76root6)
+set(heptools_version 79)
find_file(use_heptools_module UseHEPTools.cmake
HINTS ${CMAKE_CURRENT_SOURCE_DIR}/cmake)
--
GitLab