Index: trunk/examples/extended/electromagnetic/TestEm17/TestEm17.out
===================================================================
--- trunk/examples/extended/electromagnetic/TestEm17/TestEm17.out	(revision 1337)
+++ trunk/examples/extended/electromagnetic/TestEm17/TestEm17.out	(revision 1342)
@@ -5,5 +5,5 @@
 
 *************************************************************
- Geant4 version Name: geant4-09-03-ref-06    (25-June-2010)
+ Geant4 version Name: geant4-09-03-ref-09    (25-June-2010)
                       Copyright : Geant4 Collaboration
                       Reference : NIM A 506 (2003), 250-303
@@ -14,26 +14,26 @@
 ***** Table : Nb of materials = 7 *****
 
- Material: galactic     density:  0.000 kg/m3   RadL: 204727512.315 pc   Nucl.Int.Length: 113804112.800 pc   Imean:  19.200 eV   temperature:   2.73 K  pressure:   0.00 atm
+ 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  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Nucl.Int.Length: 662.680 m    Imean:  85.703 eV   temperature: 273.15 K  pressure:   1.00 atm
+ Material:      Air     density:  1.290 mg/cm3  RadL: 285.161 m    Nucl.Int.Length: 662.904 m    Imean:  85.703 eV   temperature: 273.15 K  pressure:   1.00 atm
    --->  Element: Nitrogen (N)   Z =  7.0   N =  14.0   A =  14.01 g/mole  ElmMassFraction:  70.00 %  ElmAbundance  72.71 %
    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  30.00 %  ElmAbundance  27.29 %
 
- Material:    CaCO3     density:  2.800 g/cm3   RadL:   8.581 cm   Nucl.Int.Length:  34.749 cm   Imean: 123.236 eV 
+ Material:    CaCO3     density:  2.800 g/cm3   RadL:   8.581 cm   Nucl.Int.Length:  34.767 cm   Imean: 123.236 eV 
    --->  Element: Calcium (Ca)   Z = 20.0   N =  40.1   A =  40.08 g/mole  ElmMassFraction:  40.04 %  ElmAbundance  20.00 %
    --->  Element: Carbon (C)   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction:  12.00 %  ElmAbundance  20.00 %
    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  47.96 %  ElmAbundance  60.00 %
 
- Material:   Carbon     density:  2.265 g/cm3   RadL:  18.850 cm   Nucl.Int.Length:  35.387 cm   Imean:  81.000 eV 
+ Material:   Carbon     density:  2.265 g/cm3   RadL:  18.850 cm   Nucl.Int.Length:  35.407 cm   Imean:  81.000 eV 
    --->  Element: Carbon ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material:     Iron     density:  7.870 g/cm3   RadL:   1.759 cm   Nucl.Int.Length:  16.999 cm   Imean: 286.000 eV 
+ 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  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material:      Tin     density:  7.310 g/cm3   RadL:   1.206 cm   Nucl.Int.Length:  23.531 cm   Imean: 488.000 eV 
+ Material:      Tin     density:  7.310 g/cm3   RadL:   1.206 cm   Nucl.Int.Length:  23.491 cm   Imean: 488.000 eV 
    --->  Element: Tin ( )   Z = 50.0   N = 118.7   A = 118.70 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material:       HI     density:  2.800 g/cm3   RadL:   3.049 cm   Nucl.Int.Length:  60.888 cm   Imean: 462.393 eV 
+ Material:       HI     density:  2.800 g/cm3   RadL:   3.049 cm   Nucl.Int.Length:  60.873 cm   Imean: 462.393 eV 
    --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:   0.79 %  ElmAbundance  50.00 %
    --->  Element: Iodine (I)   Z = 53.0   N = 126.9   A = 126.90 g/mole  ElmMassFraction:  99.21 %  ElmAbundance  50.00 %
@@ -71,6 +71,6 @@
       finalRange(mm)= 1, 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=   1000 PeV
+               Bragg :   Emin=        0 eV       Emax=          2 MeV
+          BetheBloch :   Emin=        2 MeV      Emax=       1000 PeV
 
 hBrems:   for  proton    SubType= 3
@@ -78,5 +78,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-               hBrem :     Emin=          0 eV         Emax=   100 TeV
+               hBrem :   Emin=        0 eV       Emax=        100 TeV
 
 hPairProd:   for  proton    SubType= 4
@@ -84,5 +84,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-           hPairProd :     Emin=          0 eV         Emax=   100 TeV
+           hPairProd :   Emin=        0 eV       Emax=        100 TeV
 
 muIoni:   for  mu+    SubType= 2
@@ -91,7 +91,7 @@
       finalRange(mm)= 1, 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=   1000 PeV
+               Bragg :   Emin=        0 eV       Emax=        200 keV
+          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
+        MuBetheBloch :   Emin=        1 GeV      Emax=       1000 PeV
 
 muBrems:   for  mu+    SubType= 3
@@ -99,5 +99,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-              MuBrem :     Emin=          0 eV         Emax=   1000 PeV
+              MuBrem :   Emin=        0 eV       Emax=       1000 PeV
 
 muPairProd:   for  mu+    SubType= 4
@@ -105,5 +105,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-          muPairProd :     Emin=          0 eV         Emax=   1000 PeV
+          muPairProd :   Emin=        0 eV       Emax=       1000 PeV
 
 muIoni:   for  mu-    SubType= 2
@@ -112,7 +112,7 @@
       finalRange(mm)= 1, 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=   1000 PeV
+            ICRU73QO :   Emin=        0 eV       Emax=        200 keV
+          BetheBloch :   Emin=      200 keV      Emax=          1 GeV
+        MuBetheBloch :   Emin=        1 GeV      Emax=       1000 PeV
 
 muBrems:   for  mu-    SubType= 3
@@ -120,5 +120,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-              MuBrem :     Emin=          0 eV         Emax=   1000 PeV
+              MuBrem :   Emin=        0 eV       Emax=       1000 PeV
 
 muPairProd:   for  mu-    SubType= 4
@@ -126,5 +126,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-          muPairProd :     Emin=          0 eV         Emax=   1000 PeV
+          muPairProd :   Emin=        0 eV       Emax=       1000 PeV
 
 hIoni:   for  pi+    SubType= 2
@@ -133,6 +133,6 @@
       finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-               Bragg :     Emin=          0 eV         Emax=   297.504 keV
-          BetheBloch :     Emin=    297.504 keV        Emax=   1000 PeV
+               Bragg :   Emin=        0 eV       Emax=    297.505 keV
+          BetheBloch :   Emin=  297.505 keV      Emax=       1000 PeV
 
 hBrems:   for  pi+    SubType= 3
@@ -140,5 +140,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-               hBrem :     Emin=          0 eV         Emax=   100 TeV
+               hBrem :   Emin=        0 eV       Emax=        100 TeV
 
 hPairProd:   for  pi+    SubType= 4
@@ -146,5 +146,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-           hPairProd :     Emin=          0 eV         Emax=   100 TeV
+           hPairProd :   Emin=        0 eV       Emax=        100 TeV
 
 hIoni:   for  pi-    SubType= 2
@@ -153,6 +153,6 @@
       finalRange(mm)= 1, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-            ICRU73QO :     Emin=          0 eV         Emax=   297.504 keV
-          BetheBloch :     Emin=    297.504 keV        Emax=   1000 PeV
+            ICRU73QO :   Emin=        0 eV       Emax=    297.505 keV
+          BetheBloch :   Emin=  297.505 keV      Emax=       1000 PeV
 
 hBrems:   for  pi-    SubType= 3
@@ -160,5 +160,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-               hBrem :     Emin=          0 eV         Emax=   100 TeV
+               hBrem :   Emin=        0 eV       Emax=        100 TeV
 
 hPairProd:   for  pi-    SubType= 4
@@ -166,5 +166,5 @@
       Lambda tables from threshold to 1000 PeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-           hPairProd :     Emin=          0 eV         Emax=   100 TeV
+           hPairProd :   Emin=        0 eV       Emax=        100 TeV
 
 Region <DefaultRegionForTheWorld> -- appears in <Iron> world volume
@@ -220,16 +220,16 @@
 Run Summary
   Number of events processed : 10000
-  User=4.58s Real=9.39s Sys=0.25s
+  User=4.08s Real=4.33s Sys=0.25s
 
  The run consists of 10000 mu+ of 10 TeV through 1 m   of Iron (density: 7.9 g/cm3 )
 
- Number of process calls --->	muIoni : 461894	muPairProd : 63886	muBrems : 665
-
- Simulation: total CrossSection = 0.52644 /cm	 MeanFreePath = 1.8995 cm 	 massicCrossSection = 0.066893 cm2/g
+ Number of process calls --->	muIoni : 461338	muPairProd : 63718	muBrems : 700
+
+ Simulation: total CrossSection = 0.52576 /cm	 MeanFreePath = 1.902 cm 	 massicCrossSection = 0.066805 cm2/g
  Theory:     total CrossSection = 0.53938 /cm	 MeanFreePath = 1.854 cm 	 massicCrossSection = 0.068536 cm2/g
 
 --------- Ranecu engine status ---------
  Initial seed (index) = 0
- Current couple of seeds = 1599898295, 1410323926
+ Current couple of seeds = 776914426, 1041080034
 ----------------------------------------
 #
@@ -248,5 +248,5 @@
  Material : Iron
  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm 
- Energy thresholds :  gamma  20.643844 keV    e-  1.295918 MeV    e+  1.2116932 MeV proton 0 eV 
+ Energy thresholds :  gamma  20.6438 keV    e-  1.29592 MeV    e+  1.21169 MeV proton 0 eV 
  Region(s) which use this couple : 
     DefaultRegionForTheWorld
@@ -258,5 +258,5 @@
 --------- Ranecu engine status ---------
  Initial seed (index) = 0
- Current couple of seeds = 1599898295, 1410323926
+ Current couple of seeds = 776914426, 1041080034
 ----------------------------------------
 Start Run processing.
@@ -284,15 +284,15 @@
 Run Summary
   Number of events processed : 10000
-  User=4.12s Real=5.27s Sys=0.24s
+  User=3.9s Real=4.13s Sys=0.23s
 
  The run consists of 10000 pi+ of 10 TeV through 1 m   of Iron (density: 7.9 g/cm3 )
 
- Number of process calls --->	hIoni : 433928	hPairProd : 59452	hBrems : 347
-
- Simulation: total CrossSection = 0.49373 /cm	 MeanFreePath = 2.0254 cm 	 massicCrossSection = 0.062735 cm2/g
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 1194972054, 916833196
+ Number of process calls --->	hIoni : 433540	hPairProd : 59301	hBrems : 363
+
+ Simulation: total CrossSection = 0.4932 /cm	 MeanFreePath = 2.0276 cm 	 massicCrossSection = 0.062669 cm2/g
+
+--------- Ranecu engine status ---------
+ Initial seed (index) = 0
+ Current couple of seeds = 1649005378, 938556311
 ----------------------------------------
 #
@@ -311,5 +311,5 @@
  Material : Iron
  Range cuts        :  gamma  1 mm     e-  1 mm     e+  1 mm  proton 0 fm 
- Energy thresholds :  gamma  20.643844 keV    e-  1.295918 MeV    e+  1.2116932 MeV proton 0 eV 
+ Energy thresholds :  gamma  20.6438 keV    e-  1.29592 MeV    e+  1.21169 MeV proton 0 eV 
  Region(s) which use this couple : 
     DefaultRegionForTheWorld
@@ -321,5 +321,5 @@
 --------- Ranecu engine status ---------
  Initial seed (index) = 0
- Current couple of seeds = 1194972054, 916833196
+ Current couple of seeds = 1649005378, 938556311
 ----------------------------------------
 Start Run processing.
@@ -347,15 +347,15 @@
 Run Summary
   Number of events processed : 10000
-  User=3.58s Real=3.91s Sys=0.24s
+  User=3.38s Real=3.61s Sys=0.23s
 
  The run consists of 10000 proton of 10 TeV through 1 m   of Iron (density: 7.9 g/cm3 )
 
- Number of process calls --->	hIoni : 433160	hPairProd : 34070	hBrems : 10
-
- Simulation: total CrossSection = 0.46724 /cm	 MeanFreePath = 2.1402 cm 	 massicCrossSection = 0.05937 cm2/g
-
---------- Ranecu engine status ---------
- Initial seed (index) = 0
- Current couple of seeds = 250921008, 1691860124
+ Number of process calls --->	hIoni : 434451	hPairProd : 33980	hBrems : 6
+
+ Simulation: total CrossSection = 0.46844 /cm	 MeanFreePath = 2.1348 cm 	 massicCrossSection = 0.059522 cm2/g
+
+--------- Ranecu engine status ---------
+ Initial seed (index) = 0
+ Current couple of seeds = 669649146, 967029133
 ----------------------------------------
 #
