************************************************************* Geant4 version Name: geant4-08-02-patch-01-ref (23-February-2007) Copyright : Geant4 Collaboration Reference : NIM A 506 (2003), 250-303 WWW : http://cern.ch/geant4 ************************************************************* /run/verbose 1 /tracking/verbose 0 # /testem/phys/addPhysics G4standard_exp # #/testem/run/setHisto 1 100 250. 750. MeV #/testem/run/setHisto 2 100 9250. 9750. MeV # /gun/particle e- /gun/energy 10 GeV # #/random/setSavingFlag true #/random/resetEngineFrom first.rndm /testem/event/printModulo 10000 # /testem/phys/setCuts 100 mm /run/initialize ------------------------------------------------------------- ---> The calorimeter is 50 layers of: Lead: 2.3 mm liquidArgon: 5.7 mm ------------------------------------------------------------- Material: Galactic density: 0.000 kg/m3 temperature: 2.73 K pressure: 0.00 atm RadLength: 204322175.595 pc ---> Element: Galactic ( ) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: Lead density: 11.350 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 5.613 mm ---> Element: Pb (Pb) Z = 82.0 N = 207.2 A = 207.22 g/mole ---> Isotope: Pb Z = 82 N = 204 A = 203.97 g/mole abundance: 1.40 % ---> Isotope: Pb Z = 82 N = 206 A = 205.97 g/mole abundance: 24.10 % ---> Isotope: Pb Z = 82 N = 207 A = 206.98 g/mole abundance: 22.10 % ---> Isotope: Pb Z = 82 N = 208 A = 207.98 g/mole abundance: 52.40 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % Material: liquidArgon density: 1.390 g/cm3 temperature: 273.15 K pressure: 1.00 atm RadLength: 14.064 cm ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole ---> Isotope: Ar Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % ---> Isotope: Ar Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % ---> Isotope: Ar Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % ------------------------------------------------------------- PhysicsList::SetCuts:CutLength : 1 mm /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.022 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.297504 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 2.24888 MeV e- 437.915 MeV e+ 406.735 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 10 cm e- 10 cm e+ 10 cm Energy thresholds : gamma 56.6382 keV e- 24.5664 MeV e+ 22.8173 MeV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 0 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 9876, 54321 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=13036.43s Real=13060.93s Sys=0.18s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.662 GeV +- 77.23 MeV 3.187 +- 0.03187 % 5.75 m +- 6.06 cm liquidArgon: 2.1839 GeV +- 53.35 MeV 7.725 +- 0.07725 % 11 m +- 27.3 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 13.4862 GeV liquidArgon: -3.64085 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1457881611, 1784478257 ---------------------------------------- # /testem/phys/setCuts 10 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 282.47765 keV e- 15.013704 MeV e+ 13.944716 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 1 cm e- 1 cm e+ 1 cm Energy thresholds : gamma 19.052126 keV e- 2.1466594 MeV e+ 2.0435144 MeV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 1 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1457881611, 1784478257 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=14831.4s Real=16163.32s Sys=0.35s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.7536 GeV +- 71.46 MeV 2.915 +- 0.02915 % 5.52 m +- 5.29 cm liquidArgon: 2.1053 GeV +- 50.71 MeV 7.617 +- 0.07617 % 10.3 m +- 25.9 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 15.3496 GeV liquidArgon: -5.49124 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 969860042, 1003221401 ---------------------------------------- # /testem/phys/setCuts 3 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 139.61279 keV e- 3.9726571 MeV e+ 3.6898008 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 3 mm e- 3 mm e+ 3 mm Energy thresholds : gamma 10.507304 keV e- 763.25382 keV e+ 726.58019 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 2 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 969860042, 1003221401 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=15624.27s Real=15834.01s Sys=0.32s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.7699 GeV +- 67.84 MeV 2.761 +- 0.02761 % 5.51 m +- 5.02 cm liquidArgon: 2.0953 GeV +- 49.8 MeV 7.517 +- 0.07517 % 10.3 m +- 25.4 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 25.6092 GeV liquidArgon: -15.7451 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 414020352, 985459431 ---------------------------------------- # /testem/phys/setCuts 1 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= 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 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 100.51071 keV e- 1.3781428 MeV e+ 1.280018 MeV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 1 mm e- 1 mm e+ 1 mm Energy thresholds : gamma 6.1783498 keV e- 342.89066 keV e+ 334.55147 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 3 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 414020352, 985459431 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=15924.69s Real=109246.38s Sys=0.33s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.7703 GeV +- 74.95 MeV 3.05 +- 0.0305 % 5.49 m +- 5.46 cm liquidArgon: 2.0944 GeV +- 50.63 MeV 7.645 +- 0.07645 % 10.4 m +- 25.9 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 32.0543 GeV liquidArgon: -22.1916 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 374501854, 1065208692 ---------------------------------------- # /testem/phys/setCuts 0.7 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 700 um e- 700 um e+ 700 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 700 um e- 700 um e+ 700 um Energy thresholds : gamma 94.715457 keV e- 1.0006654 MeV e+ 940.92951 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 700 um e- 700 um e+ 700 um Energy thresholds : gamma 5.1445663 keV e- 271.37813 keV e+ 264.77814 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 4 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 374501854, 1065208692 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=15709.55s Real=15850.83s Sys=0.47s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.7729 GeV +- 70.51 MeV 2.868 +- 0.02868 % 5.49 m +- 5.14 cm liquidArgon: 2.0931 GeV +- 49.85 MeV 7.532 +- 0.07532 % 10.4 m +- 25.7 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 32.368 GeV liquidArgon: -22.504 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1296223024, 1247217761 ---------------------------------------- # /testem/phys/setCuts 0.3 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 300 um e- 300 um e+ 300 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 300 um e- 300 um e+ 300 um Energy thresholds : gamma 62.748006 keV e- 514.73714 keV e+ 490.0045 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 300 um e- 300 um e+ 300 um Energy thresholds : gamma 3.4073784 keV e- 157.88288 keV e+ 155.95119 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 5 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1296223024, 1247217761 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=15638.62s Real=16148.84s Sys=0.43s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.776 GeV +- 75.16 MeV 3.057 +- 0.03057 % 5.52 m +- 5.49 cm liquidArgon: 2.0906 GeV +- 51.64 MeV 7.81 +- 0.0781 % 10.5 m +- 26.8 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 31.2277 GeV liquidArgon: -21.3636 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1077804822, 1397920400 ---------------------------------------- # /testem/phys/setCuts 0.1 mm /run/initialize /run/beamOn 10000 phot: Total cross sections from Sandia parametrisation. Sampling according PhotoElectric model compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV Sampling according Klein-Nishina model tables are built for gamma Lambda tables from 100 eV to 100 GeV in 90 bins. conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; sampling secondary e+e- according Bethe-Heitler model tables are built for gamma Lambda tables from 1.0219981 MeV to 100 GeV in 100 bins. msc: Model variant of multiple scattering for e- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.02 Step limitation 1 ### The process eBrem is added to the list of collaborative processes of eIoni ### The process eBrem is added to the list of collaborative processes of eIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni ### The process muBrems is added to the list of collaborative processes of muIoni ### The process muPairProd is added to the list of collaborative processes of muIoni eIoni: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions eIoni: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Delta cross sections and sampling from MollerBhabha model Good description from 1 KeV to 100 GeV. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions eBrem: tables are built for e+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Total cross sections and sampling from StandBrem model (based on the EEDL data library) Good description from 1 KeV to 100 GeV, log scale extrapolation above 100 GeV. Subcutoff sampling in 1 regions annihil: Sampling according eplus2gg model tables are built for e+ Lambda tables from 100 eV to 100 TeV in 120 bins. msc: Model variant of multiple scattering for proton Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for GenericIon LateralDisplacementFlag= 0 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 ionIoni: tables are built for GenericIon dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from BraggIon for protons below. Stopping Power data for 8 ion/material pairs are used. Step function: finalRange(mm)= 0.1, dRoverRange= 0.1, integral: 1 Subcutoff sampling in 1 regions hIoni: tables are built for anti_proton dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for mu+ Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 muIoni: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muIoni: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Bether-Bloch model for E > 0.2 MeV, parametrisation of Bragg peak below, radiative corrections for E > 1 GeV Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions muBrems: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions muPairProd: tables are built for mu- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Parametrised model Subcutoff sampling in 1 regions hIoni: tables are built for pi+ dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions msc: Model variant of multiple scattering for pi- Lambda tables from 100 eV to 100 TeV in 120 bins. LateralDisplacementFlag= 1 Skin= 0 Boundary/stepping algorithm is active with facrange= 0.2 Step limitation 0 hIoni: tables are built for pi- dE/dx and range tables from 100 eV to 100 TeV in 120 bins. Lambda tables from threshold to 100 TeV in 120 bins. Scaling relation is used from proton dE/dx and range. Delta cross sections and sampling from BetheBloch model for scaled energy > 0.29750425 MeV Parametrisation from Bragg for protons below. Step function: finalRange(mm)= 1, dRoverRange= 0.2, integral: 1 Subcutoff sampling in 1 regions ========= Table of registered couples ============================== Index : 0 used in the geometry : Yes recalculation needed : No Material : Galactic Range cuts : gamma 100 um e- 100 um e+ 100 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV Region(s) which use this couple : DefaultRegionForTheWorld Index : 1 used in the geometry : Yes recalculation needed : No Material : Lead Range cuts : gamma 100 um e- 100 um e+ 100 um Energy thresholds : gamma 29.3406 keV e- 239.94475 keV e+ 231.24491 keV Region(s) which use this couple : DefaultRegionForTheWorld Index : 2 used in the geometry : Yes recalculation needed : No Material : liquidArgon Range cuts : gamma 100 um e- 100 um e+ 100 um Energy thresholds : gamma 1.9852965 keV e- 82.220106 keV e+ 81.214144 keV Region(s) which use this couple : DefaultRegionForTheWorld ==================================================================== ### Run 6 start. --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 1077804822, 1397920400 ---------------------------------------- Start Run processing. ---> Begin Of Event: 0 Run terminated. Run Summary Number of events processed : 10000 User=16025.42s Real=16616.74s Sys=0.59s ------------------------------------------------------------ material Total Edep sqrt(E0(GeV))*rmsE/Emean total tracklen Lead: 7.7755 GeV +- 70.94 MeV 2.885 +- 0.02885 % 5.61 m +- 5.28 cm liquidArgon: 2.0926 GeV +- 49.81 MeV 7.528 +- 0.07528 % 10.6 m +- 26.1 cm ------------------------------------------------------------ Energy deposition from Energy flow balance : material Total Edep Lead: 25.3206 GeV liquidArgon: -15.4554 GeV ------------------------------------------------------------ --------- Ranecu engine status --------- Initial seed (index) = 0 Current couple of seeds = 2045646623, 1061243697 ---------------------------------------- # G4 kernel has come to Quit state.