Index: trunk/examples/extended/medical/fanoCavity/History
===================================================================
--- trunk/examples/extended/medical/fanoCavity/History	(revision 1337)
+++ trunk/examples/extended/medical/fanoCavity/History	(revision 1342)
@@ -1,3 +1,3 @@
-$Id: History,v 1.36 2010/06/07 05:40:46 perl Exp $
+$Id: History,v 1.37 2010/10/25 13:29:24 gunter Exp $
 -------------------------------------------------------------------
 
@@ -15,4 +15,7 @@
      * Reverse chronological order (last date on top), please *
      ----------------------------------------------------------
+
+25-10-2010 G.Folger  (fano-V09-03-04)
+- Correct initialisation of directinIn in ctor; needed by clhep 2.1...
 
 06-06-10 J.Perl (fano-V09-03-03)
Index: trunk/examples/extended/medical/fanoCavity/fanoCavity.out
===================================================================
--- trunk/examples/extended/medical/fanoCavity/fanoCavity.out	(revision 1337)
+++ trunk/examples/extended/medical/fanoCavity/fanoCavity.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:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.416 cm   Imean:  75.000 eV 
+ Material:    Water     density:  1.000 g/cm3   RadL:  36.092 cm   Nucl.Int.Length:  75.537 cm   Imean:  75.000 eV 
    --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
 
- Material: Water_vapor     density:  1.000 mg/cm3  RadL: 360.924 m    Nucl.Int.Length: 754.163 m    Imean:  75.000 eV   temperature: 273.15 K  pressure:   1.00 atm
+ Material: Water_vapor     density:  1.000 mg/cm3  RadL: 360.924 m    Nucl.Int.Length: 755.365 m    Imean:  75.000 eV   temperature: 273.15 K  pressure:   1.00 atm
    --->  Element: Hydrogen (H)   Z =  1.0   N =   1.0   A =   1.01 g/mole  ElmMassFraction:  11.21 %  ElmAbundance  66.67 %
    --->  Element: Oxygen (O)   Z =  8.0   N =  16.0   A =  16.00 g/mole  ElmMassFraction:  88.79 %  ElmAbundance  33.33 %
 
- 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: Graphite     density:  2.265 g/cm3   RadL:  18.850 cm   Nucl.Int.Length:  35.387 cm   Imean:  81.000 eV 
+ Material: Graphite     density:  2.265 g/cm3   RadL:  18.850 cm   Nucl.Int.Length:  35.407 cm   Imean:  81.000 eV 
    --->  Element: Graphite ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material: Graphite_gas     density:  2.265 mg/cm3  RadL: 188.496 m    Nucl.Int.Length: 353.873 m    Imean:  81.000 eV   temperature: 273.15 K  pressure:   1.00 atm
+ Material: Graphite_gas     density:  2.265 mg/cm3  RadL: 188.496 m    Nucl.Int.Length: 354.070 m    Imean:  81.000 eV   temperature: 273.15 K  pressure:   1.00 atm
    --->  Element: Graphite_gas ( )   Z =  6.0   N =  12.0   A =  12.01 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material: Aluminium     density:  2.700 g/cm3   RadL:   8.893 cm   Nucl.Int.Length:  38.879 cm   Imean: 166.000 eV 
+ 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  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
- Material: Aluminium_gas     density:  2.700 mg/cm3  RadL:  88.925 m    Nucl.Int.Length: 388.793 m    Imean: 166.000 eV   temperature: 273.15 K  pressure:   1.00 atm
+ Material: Aluminium_gas     density:  2.700 mg/cm3  RadL:  88.925 m    Nucl.Int.Length: 388.601 m    Imean: 166.000 eV   temperature: 273.15 K  pressure:   1.00 atm
    --->  Element: Aluminium_gas ( )   Z = 13.0   N =  27.0   A =  26.98 g/mole  ElmMassFraction: 100.00 %  ElmAbundance 100.00 %
 
@@ -65,15 +65,15 @@
 phot:   for  gamma    SubType= 12
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-       PhotoElectric :     Emin=          0 eV         Emax=   10 GeV
+       PhotoElectric :   Emin=        0 eV       Emax=         10 GeV
 
 compt:   for  gamma    SubType= 13
       Lambda tables from 100 eV  to 10 GeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-     myKlein-Nishina :     Emin=          0 eV         Emax=   10 GeV
+     myKlein-Nishina :   Emin=        0 eV       Emax=         10 GeV
 
 conv:   for  gamma    SubType= 14
       Lambda tables from 1.022 MeV to 10 GeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-       Bethe-Heitler :     Emin=          0 eV         Emax=   10 GeV
+        BetheHeitler :   Emin=        0 eV       Emax=         10 GeV
 
 msc:   for e-    SubType= 10
@@ -81,5 +81,5 @@
       RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-          UrbanMsc93 :     Emin=          0 eV         Emax=   10 GeV
+          UrbanMsc93 :   Emin=        0 eV       Emax=         10 GeV
 
 eIoni:   for  e-    SubType= 2
@@ -88,5 +88,5 @@
       finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-      myMollerBhabha :     Emin=          0 eV         Emax=   10 GeV
+      myMollerBhabha :   Emin=        0 eV       Emax=         10 GeV
       CSDA range table up to 10 GeV in 160 bins
 
@@ -96,5 +96,5 @@
       finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-      myMollerBhabha :     Emin=          0 eV         Emax=   10 GeV
+      myMollerBhabha :   Emin=        0 eV       Emax=         10 GeV
       CSDA range table up to 10 GeV in 160 bins
 
@@ -102,5 +102,5 @@
       Lambda tables from 100 eV  to 10 GeV in 160 bins, spline: 1
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-            eplus2gg :     Emin=          0 eV         Emax=   10 GeV
+            eplus2gg :   Emin=        0 eV       Emax=         10 GeV
 
 msc:   for proton    SubType= 10
@@ -108,5 +108,5 @@
       RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5
       ===== EM models for the G4Region  DefaultRegionForTheWorld ======
-          UrbanMsc90 :     Emin=          0 eV         Emax=   10 GeV
+          UrbanMsc90 :   Emin=        0 eV       Emax=         10 GeV
 
 hIoni:   for  proton    SubType= 2
@@ -115,6 +115,6 @@
       finalRange(mm)= 0.01, 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 GeV
+               Bragg :   Emin=        0 eV       Emax=          2 MeV
+          BetheBloch :   Emin=        2 MeV      Emax=         10 GeV
       CSDA range table up to 10 GeV in 160 bins
 
@@ -156,13 +156,13 @@
 Start Run processing.
 
- ---> NbofEvents= 1000   NbOfelectr= 70538   Tkin= 590.669 keV (0.000 %)   NbOfelec in cav= 4726   Dose/EnFluence= 30.091 cm2/g (0.000 %)
-
- ---> NbofEvents= 2000   NbOfelectr= 141090   Tkin= 589.502 keV (-0.198 %)   NbOfelec in cav= 9358   Dose/EnFluence= 28.953 cm2/g (-3.782 %)
-
- ---> NbofEvents= 3000   NbOfelectr= 211481   Tkin= 589.178 keV (-0.055 %)   NbOfelec in cav= 14156   Dose/EnFluence= 28.931 cm2/g (-0.076 %)
+ ---> NbofEvents= 1000   NbOfelectr= 70332   Tkin= 589.082 keV (0.000 %)   NbOfelec in cav= 4589   Dose/EnFluence= 26.638 cm2/g (0.000 %)
+
+ ---> NbofEvents= 2000   NbOfelectr= 141252   Tkin= 588.179 keV (-0.153 %)   NbOfelec in cav= 9374   Dose/EnFluence= 28.803 cm2/g (8.129 %)
+
+ ---> NbofEvents= 3000   NbOfelectr= 211798   Tkin= 588.153 keV (-0.004 %)   NbOfelec in cav= 13975   Dose/EnFluence= 28.118 cm2/g (-2.379 %)
 Run terminated.
 Run Summary
   Number of events processed : 4000
-  User=32.81s Real=48.36s Sys=0s
+  User=30.34s Real=30.35s Sys=0s
 
  ======================== run summary ======================
@@ -174,30 +174,30 @@
  ============================================================
 
- Process calls frequency --->  compt= 282107  Transportation= 77642  msc= 2087713  eIoni= 1755986
+ Process calls frequency --->  compt= 281849  Transportation= 77349  msc= 2094229  eIoni= 1736138
 
  Gamma crossSections in wall material :  compt= 63.447 cm2/g   --> total= 63.447 cm2/g
 
- Mean energy of secondary e- = 588.9126 keV +- 0.1029 %  (--> range in wall material = 2.2146 mm )
- Mass_energy_transfer coef: 29.8918 cm2/g
-
- StoppingPower in wall   = 1.9720 MeV*cm2/g
-               in cavity = 1.9720 MeV*cm2/g
+ Mean energy of secondary e- = 588.2741 keV +- 0.1028 %  (--> range in wall material = 2.2114 mm )
+ Mass_energy_transfer coef: 29.8594 cm2/g
+
+ StoppingPower in wall   = 1.9724 MeV*cm2/g
+               in cavity = 1.9724 MeV*cm2/g
 
  Charged particle flow in cavity :
-      Enter --> nbParticles = 18882	 Energy = 9.5906 GeV
-      Exit  --> nbParticles = 18879	 Energy = 9.5913 GeV
-
- Total edep in cavity = 13.0468 MeV +- 3.3776 %	 Total charged trackLength = 51.4975 m     (mean value = 2.7273 mm )
- Total dose in cavity = 20.7646 MeV/mg
- Dose/EnergyFluence   = 29.3552 cm2/g
-
- (Dose/EnergyFluence)/Mass_energy_transfer = 0.98205 +- 0.03318
-
- StepSize of ch. tracks in wall   = 122.8120 um  +- 121.5387 um 	 (nbSteps/track = 13.6004)
- StepSize of ch. tracks in cavity = 668.5550 um  +- 853.6674 um 	 (nbSteps/track = 4.0794)
+      Enter --> nbParticles = 18697	 Energy = 9.4723 GeV
+      Exit  --> nbParticles = 18694	 Energy = 9.4730 GeV
+
+ Total edep in cavity = 12.4835 MeV +- 2.3625 %	 Total charged trackLength = 51.1392 m     (mean value = 2.7352 mm )
+ Total dose in cavity = 19.8681 MeV/mg
+ Dose/EnergyFluence   = 28.0878 cm2/g
+
+ (Dose/EnergyFluence)/Mass_energy_transfer = 0.94067 +- 0.02224
+
+ StepSize of ch. tracks in wall   = 122.6753 um  +- 120.8436 um 	 (nbSteps/track = 13.5663)
+ StepSize of ch. tracks in cavity = 668.6433 um  +- 847.1076 um 	 (nbSteps/track = 4.0906)
 
 --------- Ranecu engine status ---------
  Initial seed (index) = 0
- Current couple of seeds = 1335550696, 1317159614
+ Current couple of seeds = 1568745475, 2068190730
 ----------------------------------------
 UserDetectorConstruction deleted.
Index: trunk/examples/extended/medical/fanoCavity/src/SteppingAction.cc
===================================================================
--- trunk/examples/extended/medical/fanoCavity/src/SteppingAction.cc	(revision 1337)
+++ trunk/examples/extended/medical/fanoCavity/src/SteppingAction.cc	(revision 1342)
@@ -24,6 +24,6 @@
 // ********************************************************************
 //
-// $Id: SteppingAction.cc,v 1.4 2007/10/29 17:09:53 maire Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: SteppingAction.cc,v 1.5 2010/10/25 13:29:08 gunter Exp $
+// GEANT4 tag $Name: examples-V09-03-09 $
 //
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
@@ -51,5 +51,5 @@
   first = true;
   trackSegm = 0.;
-  directionIn = 0.;
+  directionIn = G4ThreeVector(0.,0.,0.);
 }
 
