Index: trunk/source/processes/electromagnetic/standard/include/G4ComptonScattering52.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4ComptonScattering52.hh	(revision 1316)
+++ 	(revision )
@@ -1,162 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-//
-// $Id: G4ComptonScattering52.hh,v 1.4 2007/05/16 14:00:56 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-//------------------ G4ComptonScattering52 physics process -----------------------
-//                   by Michel Maire, April 1996
-//
-// 10-06-96, updated by M.Maire
-// 21-06-96, SetCuts implementation, M.Maire
-// 06-01-97, crossection table + meanfreepath table, M.Maire
-// 17-02-97, New Physics scheme
-// 25-02-97, GetMeanFreePath() now is public function
-// 12-03-97, new physics scheme again
-// 13-08-98, new methods SetBining()  PrintInfo()
-// 03-08-01, new methods Store/Retrieve PhysicsTable (mma)
-// 06-08-01, BuildThePhysicsTable() called from constructor (mma)
-// 19-09-01, come back to previous ProcessName "compt"
-// 20-09-01, DoIt: fminimalEnergy = 1*eV (mma)
-// 01-10-01, come back to BuildPhysicsTable(const G4ParticleDefinition&)
-// 13-08-04, suppress icc file; make public ComputeCrossSectionPerAtom()   (mma)
-// 09-11-04, Remove Retrieve tables (V.Ivantchenko)
-// 04-05-05, Add 52 to class name (V.Ivanchenko)
-// -----------------------------------------------------------------------------
-
-// class description
-//
-// inherit from G4VDiscreteProcess
-//
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#ifndef G4ComptonScattering52_h
-#define G4ComptonScattering52_h 1
-
-#include "G4ios.hh" 
-#include "globals.hh"
-#include "Randomize.hh"
-#include "G4VDiscreteProcess.hh"
-#include "G4PhysicsTable.hh"
-#include "G4PhysicsLogVector.hh"
-#include "G4Element.hh"
-#include "G4Gamma.hh"
-#include "G4Electron.hh"
-#include "G4Step.hh"
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4ComptonScattering52 : public G4VDiscreteProcess
-
-{
-  public:  // with description
-
-  G4ComptonScattering52(const G4String& processName ="compt",
-		            G4ProcessType type = fElectromagnetic);
-
-    ~G4ComptonScattering52();
-
-     G4bool IsApplicable(const G4ParticleDefinition&);
-       // true for Gamma only.
-
-     void SetPhysicsTableBining(G4double lowE, G4double highE, G4int nBins);
-       // Allows to define the binning of the PhysicsTables,
-       // before to build them.
-
-     void BuildPhysicsTable(const G4ParticleDefinition&);
-       // It builds the total CrossSectionPerAtom table, for Gamma,
-       // and for every element contained in the elementTable.
-       // It builds the MeanFreePath table, for Gamma,
-       // and for every material contained in the materialTable.
-
-     G4bool StorePhysicsTable(const G4ParticleDefinition* ,
-			      const G4String& directory, G4bool);
-       // store CrossSection and MeanFreePath tables into an external file
-       // specified by 'directory' (must exist before invokation)
-
-       //G4bool RetrievePhysicsTable(const G4ParticleDefinition* ,
-       //			 const G4String& directory, G4bool);
-       // retrieve CrossSection and MeanFreePath tables from an external file
-       // specified by 'directory'
-
-     void PrintInfoDefinition();
-       // Print few lines of informations about the process: validity range,
-       // origine ..etc..
-       // Invoked by BuildThePhysicsTable().
-
-     G4double GetMeanFreePath(const G4Track& aTrack,
-                              G4double previousStepSize,
-                              G4ForceCondition* condition);
-       // It returns the MeanFreePath of the process for the current track :
-       // (energy, material)
-       // The previousStepSize and G4ForceCondition* are not used.
-       // This function overloads a virtual function of the base class.
-       // It is invoked by the ProcessManager of the Particle.
-
-     G4double GetCrossSectionPerAtom(G4DynamicParticle* aDynamicGamma,
-                                         G4Element*         anElement);
-       // It returns the total CrossSectionPerAtom of the process, 
-       // for the current DynamicGamma (energy), in anElement.					 
-					  
-     G4VParticleChange* PostStepDoIt(const G4Track& aTrack,
-                                    const G4Step& aStep);
-       // It computes the final state of the process (at end of step),
-       // returned as a ParticleChange object.			    
-       // This function overloads a virtual function of the base class.
-       // It is invoked by the ProcessManager of the Particle.
-
-     virtual
-     G4double ComputeCrossSectionPerAtom(G4double GammaEnergy, 
-                                         G4double AtomicNumber);
-
-     G4double ComputeMeanFreePath(G4double GammaEnergy, 
-                                  G4Material* aMaterial);
-  private:
-  
-     // hide assignment operator as private 
-     G4ComptonScattering52& operator=(const G4ComptonScattering52 &right);
-     G4ComptonScattering52(const G4ComptonScattering52& );
-                                          
-  private:
-     
-     G4PhysicsTable* theCrossSectionTable;    // table for crosssection
-     G4PhysicsTable* theMeanFreePathTable;    // table for mean free path
-       
-     G4double LowestEnergyLimit;      // low  energy limit of the tables
-     G4double HighestEnergyLimit;     // high energy limit of the tables
-     G4int NumbBinTable;              // number of bins in the tables
-     
-  protected:   
-     
-     G4double fminimalEnergy;         // minimalEnergy of produced particles
-};
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-  
-#endif
- 
Index: trunk/source/processes/electromagnetic/standard/include/G4MscModel71.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4MscModel71.hh	(revision 1316)
+++ 	(revision )
@@ -1,171 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-// $Id: G4MscModel71.hh,v 1.6 2009/11/01 13:04:12 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-// -------------------------------------------------------------------
-//
-//
-// GEANT4 Class header file
-//
-//
-// File name:     G4MscMOdel71
-//
-// Author:        Laszlo Urban
-//
-// Creation date: 01.03.2001
-//
-// Modifications:
-//
-// 27-03-03 Move model part from G4MultipleScattering (V.Ivanchenko)
-// 27-03-03 Rename (V.Ivanchenko)
-//
-// 05-08-03 angle distribution has been modified (L.Urban)
-// 26-11-03 new data member currentRange (L.Urban)
-// 01-03-04 changes in data members + signature changed in SampleCosineTheta
-// 11-03-04 changes in data members (L.Urban)
-// 23-04-04 changes in data members and in signature of SampleCosineTheta
-//          (L.Urban)
-// 17-08-04 name of data member facxsi changed to factail (L.Urban)
-// 08-04-05 Major optimisation of internal interfaces (V.Ivantchenko)
-// 15-04-05 optimize internal interface - add SampleSecondaries method (V.Ivanchenko)
-// 03-10-05 Model is freezed with the name McsModel71 (V.Ivanchenko)
-// 17-02-06 Save table of transport cross sections not mfp (V.Ivanchenko)
-// 07-03-06 Create G4UrbanMscModel and move there step limit calculation (V.Ivanchenko)
-//
-
-//
-// Class Description:
-//
-// Implementation of the model of multiple scattering based on
-// H.W.Lewis Phys Rev 78 (1950) 526 and L.Urban model
-
-// -------------------------------------------------------------------
-//
-
-#ifndef G4MscModel71_h
-#define G4MscModel71_h 1
-
-#include "G4VEmModel.hh"
-
-class G4ParticleChangeForMSC;
-class G4Navigator;
-
-class G4MscModel71 : public G4VEmModel
-{
-
-public:
-
-  G4MscModel71(G4double&, G4double&, G4double&, G4double&, G4bool&,
-	       const G4String& nam = "UrbanMsc71");
-
-  virtual ~G4MscModel71();
-
-  virtual void Initialise(const G4ParticleDefinition*, const G4DataVector&);
-
-  virtual G4double ComputeCrossSectionPerAtom(
-                             const G4ParticleDefinition* particle,
-                                   G4double KineticEnergy,
-                                   G4double AtomicNumber,
-                                   G4double AtomicWeight=0., 
-				   G4double cut =0.,
-				   G4double emax=DBL_MAX);
-
-  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
-				 const G4MaterialCutsCouple*,
-				 const G4DynamicParticle*,
-				 G4double length,
-				 G4double safety);
-
-  G4double GeomPathLength(G4PhysicsTable* theLambdaTable,
-                    const G4MaterialCutsCouple* couple,
-                    const G4ParticleDefinition* particle,
-                          G4double& kineticEnergy,
-                          G4double lambda,
-                          G4double range,
-                          G4double truePathLength);
-
-  G4double TrueStepLength(G4double geomStepLength);
-
-  G4double SampleCosineTheta(G4double trueStepLength,G4double KineticEnergy);
-
-  G4double SampleDisplacement();
-
-  void SetLateralDisplasmentFlag(G4bool val);
-
-private:
-
-  // hide assignment operator
-  G4MscModel71 & operator=(const  G4MscModel71 &right);
-  G4MscModel71(const  G4MscModel71&);
-
-  const G4ParticleDefinition* particle;
-  G4ParticleChangeForMSC*     fParticleChange;
-  G4Navigator*                navigator;
-
-  G4double mass;
-  G4double charge;
-  G4double massRate;
-
-  G4double taubig;
-  G4double tausmall;
-  G4double taulim;
-  G4double currentTau;
-  G4double dtrl;
-  G4double NuclCorrPar;
-  G4double FactPar;
-  G4double factail ;
-
-  G4double sigmafactor;
-  G4double b;
-  G4double xsi;
-
-  G4double lambda0;
-  G4double tPathLength;
-  G4double par1,par2,par3 ;
-
-  G4double stepmin ;
-
-  G4double currentKinEnergy;
-  G4double currentRange ; 
-  G4double currentRadLength;
-
-  G4bool   samplez;
-  G4bool   latDisplasment;
-  G4bool   isInitialized;
-};
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-inline void G4MscModel71::SetLateralDisplasmentFlag(G4bool val)
-{
-  latDisplasment = val;
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-#endif
-
Index: trunk/source/processes/electromagnetic/standard/include/G4MultipleScattering.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4MultipleScattering.hh	(revision 1316)
+++ 	(revision )
@@ -1,147 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-// $Id: G4MultipleScattering.hh,v 1.36 2008/03/10 10:39:21 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-// -----------------------------------------------------------------------------
-//
-// GEANT4 Class header file
-//
-// File name:     G4MultipleScattering
-//
-// Author:        Laszlo Urban
-//
-// Creation date: March 2001
-// 
-// Modifications:
-//
-// 07-08-01 new methods Store/Retrieve PhysicsTable
-// 23-08-01 new angle and z distribution,energy dependence reduced,
-//          Store,Retrieve methods commented out temporarily, L.Urban
-// 11-09-01 G4MultipleScatteringx put as default: G4MultipleScattering
-//          Store,Retrieve methods reactived (mma)
-// 13-09-01 Unused TrueToGeomTransformation method deleted,
-//          class description (L.Urban)
-// 19-09-01 come back to previous process name msc
-// 17-04-02 NEW angle distribution + boundary algorithm modified, L.Urban
-// 22-04-02 boundary algorithm modified -> important improvement in timing !!!!
-//          (L.Urban)
-// 24-05-02 changes in data members, L.Urban
-// 30-10-02 changes in data members, L.Urban
-// 20-01-03 Migrade to cut per region (V.Ivanchenko)
-// 05-02-03 changes in data members, L.Urban
-// 28-03-03 Move to model design (V.Ivanchenko)
-// 18-04-03 Change name (V.Ivanchenko)
-// 16-06-03: ShortLived are not applicable any more (V.Ivanchenko)
-// 17-08-04 name of data member facxsi changed to factail together
-//          with the corresponding set function (L.Urban)
-// 08-11-04 Migration to new interface of Store/Retrieve tables (V.Ivantchenko)
-// 15-04-05 optimize internal interfaces (V.Ivanchenko)
-// 02-10-05 new algorithm for step limitation, new data members (L.Urban)
-// 05-10-05 value of data member tlimitmin has been changed (L.Urban)
-// 23-10-05 new Boolean data member prec (false ~ 7.1 like, true new step
-//          limit in TruePathLengthLimit, L.Urban)
-// 25-10-05 prec renamed to steppingAlgorithm, set function triggers
-//          'default' facrange too, true - 0.02, false - 0.2 (L.Urban)
-// 26-10-05 the above is put in the function MscStepLimitation() (mma)
-// 05-11-05 new data member rangecut (L.Urban)
-// 13-11-05 some code cleaning (L.Urban)
-// 07-12-05 GeomLimit is protected instead of public
-// 11-12-05 data menber rangecut removed (L.Urban)
-// 19-01-07 tlimitmin = facrange*50*micrometer, i.e. it depends on the
-//          value of facrange (L.Urban)
-// 16-02-06 set function for data member factail (L.Urban)
-// 13-10-06 data member factail removed, new data member skin
-//          together with set function, data member tkinlimit
-//          changed to lambdalimit (L.Urban)
-// 07-03-07 remove method SetSkin, because it is now in the base class (VI)
-//
-//------------------------------------------------------------------------------
-//
-
-// class description
-//
-//  The class simulates the multiple scattering for any kind
-//  of charged particle.
-//
-// class description - end
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#ifndef G4MultipleScattering_h
-#define G4MultipleScattering_h 1
-
-#include "G4VMultipleScattering.hh"
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4VMscModel;
-
-class G4MultipleScattering : public G4VMultipleScattering
-
-{
-public:    // with description
-
-  G4MultipleScattering(const G4String& processName="msc");
-
-  virtual ~G4MultipleScattering();
-
-  // returns true for charged particles, false otherwise
-  G4bool IsApplicable (const G4ParticleDefinition& p);
-
-  // Print few lines of informations about the process: validity range,
-  void PrintInfo();
-
-  // geom. step length distribution should be sampled or not
-  void Setsamplez(G4bool value) { samplez = value;};
-
-  // to reduce the energy/step dependence
-  void Setdtrl(G4double value) { dtrl = value;};
-
-  // 'soften' step limitation above lambdalimit
-  void SetLambdalimit(G4double value) { lambdalimit = value;};
-
-protected:
-
-  // This function initialise models
-  void InitialiseProcess(const G4ParticleDefinition*);
-
-private:        // data members
-
-  G4VMscModel* mscUrban;
-
-  G4double lambdalimit;
-  G4double dtrl;
-
-  G4bool   samplez;
-  G4bool   isInitialized;
-
-};
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#endif
Index: trunk/source/processes/electromagnetic/standard/include/G4MultipleScattering71.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4MultipleScattering71.hh	(revision 1316)
+++ 	(revision )
@@ -1,252 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-// $Id: G4MultipleScattering71.hh,v 1.6 2008/07/16 11:27:41 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-//
-//------------- G4MultipleScattering71 physics process --------------------------
-//               by Laszlo Urban, March 2001
-//
-// 07-08-01 new methods Store/Retrieve PhysicsTable
-// 23-08-01 new angle and z distribution,energy dependence reduced,
-//          Store,Retrieve methods commented out temporarily, L.Urban
-// 11-09-01 G4MultipleScatteringx put as default: G4MultipleScattering
-//          Store,Retrieve methods reactived (mma)
-// 13-09-01 Unused TrueToGeomTransformation method deleted,
-//          class description (L.Urban)
-// 19-09-01 come back to previous process name msc
-// 17-04-02 NEW angle distribution + boundary algorithm modified, L.Urban
-// 22-04-02 boundary algorithm modified -> important improvement in timing !!!!
-//          (L.Urban)
-// 24-05-02 changes in data members, L.Urban
-// 30-10-02 changes in data members, L.Urban
-// 20-01-03 Migrade to cut per region (V.Ivanchenko)
-// 05-02-03 changes in data members, L.Urban
-// 28-03-03 Move to model design (V.Ivanchenko)
-// 18-04-03 Change name (V.Ivanchenko)
-// 16-06-03: ShortLived are not applicable any more (V.Ivanchenko)
-// 17-08-04 name of data member facxsi changed to factail together
-//          with the corresponding set function (L.Urban)
-// 08-11-04 Migration to new interface of Store/Retrieve tables (V.Ivantchenko)
-// 15-04-05 optimize internal interfaces (V.Ivanchenko)
-// 03-10-05 Process is freezed with the name 71 (V.Ivanchenko)
-// 07-03-06 Create G4UrbanMscModel and move there step limit calculation (V.Ivanchenko)
-//
-//------------------------------------------------------------------------------
-//
-// $Id: G4MultipleScattering71.hh,v 1.6 2008/07/16 11:27:41 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-
-// class description
-//
-//  The class simulates the multiple scattering for any kind
-//  of charged particle.
-//
-// class description - end
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#ifndef G4MultipleScattering71_h
-#define G4MultipleScattering71_h 1
-
-#include "G4VMultipleScattering.hh"
-#include "G4MscModel71.hh"
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4MultipleScattering71 : public G4VMultipleScattering
-
-{
-public:    // with description
-
-  G4MultipleScattering71(const G4String& processName="msc");
-
-  virtual ~G4MultipleScattering71();
-
-  // returns true for charged particles, false otherwise
-  G4bool IsApplicable (const G4ParticleDefinition& p);
-
-  virtual G4VParticleChange* AlongStepDoIt(const G4Track&, const G4Step&);
-
-  virtual G4VParticleChange* PostStepDoIt(const G4Track&, const G4Step&);
-
-  G4double TruePathLengthLimit(const G4Track&  track,
-			       G4double& lambda,
-			       G4double  currentMinimalStep);
-
-  // Print few lines of informations about the process: validity range,
-  void PrintInfo();
-
-  // geom. step length distribution should be sampled or not
-  void Setsamplez(G4bool value);
-
-  // activate boundary algorithm
-  void SetBoundary(G4bool value);
-
-  // to reduce the energy/step dependence
-  void Setdtrl(G4double value);
-
-  void Setfactail(G4double value);
-
-  // Steplimit after boundary crossing = facrange*range
-  // estimated nb of steps at boundary nsmallstep = 1/facrange
-  void SetFacrange(G4double val);
-
-  // corrs to transport cross section for high energy
-  void SetNuclCorrPar(G4double val);
-  void SetFactPar(G4double val);
-
-protected:
-
-  // This function initialise models
-  void InitialiseProcess(const G4ParticleDefinition*);
-
-  // This method is used for tracking, it returns step limit
-  virtual G4double GetContinuousStepLimit(const G4Track& track,
-                                        G4double previousStepSize,
-                                        G4double currentMinimalStep,
-                                        G4double& currentSafety);
-
-private:
-
-  //  hide assignment operator as  private
-  G4MultipleScattering71 & operator = (const G4MultipleScattering71 &right);
-  G4MultipleScattering71 ( const G4MultipleScattering71 &);
-
-private:        // data members
-
-  G4MscModel71* model;
-
-  G4double lowKineticEnergy;
-  G4double highKineticEnergy;
-  G4int    totBins;
-
-  G4double truePathLength;
-  G4double geomPathLength;
-  G4double trueStepLength;
-  G4double range;
-
-  G4double facrange;
-  G4double tlimit;
-  G4double tlimitmin;
-  G4double dtrl;
-  G4double NuclCorrPar;
-  G4double FactPar;
-  G4double factail;
-  G4double cf;
-
-  G4int    stepnolastmsc;
-  G4int    nsmallstep;
-
-  G4bool   samplez;
-  G4bool   boundary;
-  G4bool   isInitialized;
-
-};
-
-//--------------------------------------------------------------------
-//  inline methods
-//--------------------------------------------------------------------
-
-inline G4bool G4MultipleScattering71::IsApplicable (const G4ParticleDefinition& p)
-{
-  return (p.GetPDGCharge() != 0.0 && !p.IsShortLived());
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-inline G4VParticleChange* G4MultipleScattering71::AlongStepDoIt(
-                                                        const G4Track&,
-                                                        const G4Step& step)
-{
-  G4double geomStepLength = step.GetStepLength();
-  if((geomStepLength == geomPathLength) && (truePathLength <= range))
-     trueStepLength = truePathLength;
-  else
-     trueStepLength = model->TrueStepLength(geomStepLength);
-  fParticleChange.ProposeTrueStepLength(trueStepLength);
-  return &fParticleChange;
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-inline G4VParticleChange* G4MultipleScattering71::PostStepDoIt(const G4Track& track,
-							       const G4Step& step)
-{
-  fParticleChange.Initialize(track);
-  std::vector<G4DynamicParticle*>* p = 0;
-  model->SampleSecondaries(p, CurrentMaterialCutsCouple(),track.GetDynamicParticle(),
-		    step.GetStepLength(),step.GetPostStepPoint()->GetSafety());
-  return &fParticleChange;
-}
-
-// geom. step length distribution should be sampled or not
-inline void G4MultipleScattering71::Setsamplez(G4bool value)
-{
-  samplez = value;
-}
-
-// activate boundary algorithm
-inline void G4MultipleScattering71::SetBoundary(G4bool value)
-{
-  boundary = value;
-}
-
-// to reduce the energy/step dependence
-inline void G4MultipleScattering71::Setdtrl(G4double value)                
-{
-  dtrl = value;
-}
-
-inline void G4MultipleScattering71::Setfactail(G4double value)              
-{
-  factail = value;
-}
-
-// Steplimit after boundary crossing = facrange*range
-// estimated nb of steps at boundary nsmallstep = 1/facrange
-inline void G4MultipleScattering71::SetFacrange(G4double val)              
-{
-  facrange = val;
-  nsmallstep = G4int(std::log((cf+facrange-1.)/facrange)/std::log(cf))+1;
-}
-
-// corrs to transport cross section for high energy
-inline  void G4MultipleScattering71::SetNuclCorrPar(G4double val)            
-{
-  NuclCorrPar = val;
-}
-
-inline  void G4MultipleScattering71::SetFactPar(G4double val)                
-{
-  FactPar = val;
-}
-
-#endif
-
-
-
Index: trunk/source/processes/electromagnetic/standard/include/G4PolarizedComptonScattering.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4PolarizedComptonScattering.hh	(revision 1316)
+++ 	(revision )
@@ -1,78 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-//
-// $Id: G4PolarizedComptonScattering.hh,v 1.9 2006/06/29 19:51:14 gunter Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-// --------- G4PolarizedComptonScattering physics process ----------------------
-//                   by Vicente Lara, March 1998
-//
-// -----------------------------------------------------------------------------
-
-// class description
-//
-// inherit from G4ComptonScattering
-//
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#ifndef G4PolarizedComptonScattering_h
-#define G4PolarizedComptonScattering_h 1
-
-#include "G4ComptonScattering52.hh"
-#include "G4EnergyLossTables.hh"
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4PolarizedComptonScattering : public G4ComptonScattering52
-{
- public:  // with description
-
-  G4PolarizedComptonScattering(const G4String& processName = "polarCompt");
-
-  G4VParticleChange* PostStepDoIt(const G4Track& aTrack, const G4Step& aStep);
-
- private:
-  
-  G4ThreeVector SetNewPolarization(G4double, G4double,
-  				   G4double, G4double,
-  				   G4ThreeVector& );
-
-  G4double SetPhi(G4double, G4double, G4double, G4double);
-
-
-  void SystemOfRefChange(G4ThreeVector&, G4ThreeVector&, G4ThreeVector&, 
-                         G4ThreeVector&);
-
-  // hide assignment operator as private 
-  G4PolarizedComptonScattering& operator=(const G4PolarizedComptonScattering &right);
-  G4PolarizedComptonScattering(const G4PolarizedComptonScattering& );
-
-};
-  
-#endif
- 
Index: trunk/source/processes/electromagnetic/standard/include/G4UrbanMscModel.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4UrbanMscModel.hh	(revision 1316)
+++ 	(revision )
@@ -1,229 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-// $Id: G4UrbanMscModel.hh,v 1.35 2009/04/29 13:30:22 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-// -------------------------------------------------------------------
-//
-//
-// GEANT4 Class header file
-//
-//
-// File name:     G4UrbanMscModel
-//
-// Author:        Laszlo Urban
-//
-// Creation date: 06.03.2008
-//
-// Modifications:
-//
-// 28-04-09  move G4UrbanMscModel2 from the g49.2 to G4UrbanMscModel. 
-//           now it is frozen (V.Ivanchenk0)
-//
-// Class Description:
-//
-// Implementation of the model of multiple scattering based on
-// H.W.Lewis Phys Rev 78 (1950) 526 and L.Urban model
-
-// -------------------------------------------------------------------
-//
-
-#ifndef G4UrbanMscModel_h
-#define G4UrbanMscModel_h 1
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#include "G4VMscModel.hh"
-#include "G4PhysicsTable.hh"
-#include "G4MscStepLimitType.hh"
-
-class G4ParticleChangeForMSC;
-class G4SafetyHelper;
-class G4LossTableManager;
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4UrbanMscModel : public G4VMscModel
-{
-
-public:
-
-  G4UrbanMscModel(const G4String& nam = "UrbanMscUni");
-
-  virtual ~G4UrbanMscModel();
-
-  void Initialise(const G4ParticleDefinition*, const G4DataVector&);
-
-  G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition* particle,
-				      G4double KineticEnergy,
-				      G4double AtomicNumber,
-				      G4double AtomicWeight=0., 
-				      G4double cut =0.,
-				      G4double emax=DBL_MAX);
-
-  void SampleScattering(const G4DynamicParticle*,
-			G4double safety);
-
-  G4double ComputeTruePathLengthLimit(const G4Track& track,
-				      G4PhysicsTable* theLambdaTable,
-				      G4double currentMinimalStep);
-
-  G4double ComputeGeomPathLength(G4double truePathLength);
-
-  G4double ComputeTrueStepLength(G4double geomStepLength);
-
-  G4double ComputeTheta0(G4double truePathLength,
-                         G4double KineticEnergy);
-
-private:
-
-  G4double SimpleScattering(G4double xmeanth, G4double x2meanth);
-
-  G4double SampleCosineTheta(G4double trueStepLength, G4double KineticEnergy);
-
-  G4double SampleDisplacement();
-
-  G4double LatCorrelation();
-
-  inline G4double GetLambda(G4double kinEnergy);
-
-  inline void SetParticle(const G4ParticleDefinition*);
-
-  inline void UpdateCache();
-
-  //  hide assignment operator
-  G4UrbanMscModel & operator=(const  G4UrbanMscModel &right);
-  G4UrbanMscModel(const  G4UrbanMscModel&);
-
-  const G4ParticleDefinition* particle;
-  G4ParticleChangeForMSC*     fParticleChange;
-
-  G4PhysicsTable*             theLambdaTable;
-  const G4MaterialCutsCouple* couple;
-  G4LossTableManager*         theManager;
-
-  G4double mass;
-  G4double charge,ChargeSquare;
-  G4double masslimite,lambdalimit,fr;
-
-  G4double taubig;
-  G4double tausmall;
-  G4double taulim;
-  G4double currentTau;
-  G4double tlimit;
-  G4double tlimitmin;
-  G4double tlimitminfix;
-  G4double tgeom;
-
-  G4double geombig;
-  G4double geommin;
-  G4double geomlimit;
-  G4double skindepth;
-  G4double smallstep;
-
-  G4double presafety;
-
-  G4double lambda0;
-  G4double lambdaeff;
-  G4double tPathLength;
-  G4double zPathLength;
-  G4double par1,par2,par3;
-
-  G4double stepmin;
-
-  G4double currentKinEnergy;
-  G4double currentRange; 
-  G4double rangeinit;
-  G4double currentRadLength;
-
-  G4double theta0max,rellossmax;
-  G4double third;
-
-  G4int    currentMaterialIndex;
-
-  G4double y;
-  G4double Zold;
-  G4double Zeff,Z2,Z23,lnZ;
-  G4double coeffth1,coeffth2;
-  G4double coeffc1,coeffc2;
-  G4double scr1ini,scr2ini,scr1,scr2;
-
-  G4bool   isInitialized;
-  G4bool   inside;
-  G4bool   insideskin;
-};
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-G4double G4UrbanMscModel::GetLambda(G4double e)
-{
-  G4double x;
-  if(theLambdaTable) {
-    G4bool b;
-    x = ((*theLambdaTable)[currentMaterialIndex])->GetValue(e, b);
-  } else {
-    x = CrossSection(couple,particle,e);
-  }
-  if(x > DBL_MIN) x = 1./x;
-  else            x = DBL_MAX;
-  return x;
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-void G4UrbanMscModel::SetParticle(const G4ParticleDefinition* p)
-{
-  if (p != particle) {
-    particle = p;
-    mass = p->GetPDGMass();
-    charge = p->GetPDGCharge()/eplus;
-    ChargeSquare = charge*charge;
-  }
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-void G4UrbanMscModel::UpdateCache()                                   
-{
-    lnZ = std::log(Zeff);
-    coeffth1 = 0.885+lnZ*(0.104-0.0170*lnZ);
-    coeffth2 = 0.028+lnZ*(0.012-0.00125*lnZ);
-    coeffc1  = 2.134-lnZ*(0.1045-0.00602*lnZ);
-    coeffc2  = 0.001126-lnZ*(0.0001089+0.0000247*lnZ);
-    Z2 = Zeff*Zeff;
-    Z23 = std::exp(2.*lnZ/3.);
-    scr1     = scr1ini*Z23;
-    scr2     = scr2ini*Z2*ChargeSquare;
-  //  lastMaterial = couple->GetMaterial();
-    Zold = Zeff;
-}
-
-#endif
-
Index: trunk/source/processes/electromagnetic/standard/include/G4UrbanMscModel2.hh
===================================================================
--- trunk/source/processes/electromagnetic/standard/include/G4UrbanMscModel2.hh	(revision 1316)
+++ 	(revision )
@@ -1,232 +1,0 @@
-//
-// ********************************************************************
-// * License and Disclaimer                                           *
-// *                                                                  *
-// * The  Geant4 software  is  copyright of the Copyright Holders  of *
-// * the Geant4 Collaboration.  It is provided  under  the terms  and *
-// * conditions of the Geant4 Software License,  included in the file *
-// * LICENSE and available at  http://cern.ch/geant4/license .  These *
-// * include a list of copyright holders.                             *
-// *                                                                  *
-// * Neither the authors of this software system, nor their employing *
-// * institutes,nor the agencies providing financial support for this *
-// * work  make  any representation or  warranty, express or implied, *
-// * regarding  this  software system or assume any liability for its *
-// * use.  Please see the license in the file  LICENSE  and URL above *
-// * for the full disclaimer and the limitation of liability.         *
-// *                                                                  *
-// * This  code  implementation is the result of  the  scientific and *
-// * technical work of the GEANT4 collaboration.                      *
-// * By using,  copying,  modifying or  distributing the software (or *
-// * any work based  on the software)  you  agree  to acknowledge its *
-// * use  in  resulting  scientific  publications,  and indicate your *
-// * acceptance of all terms of the Geant4 Software license.          *
-// ********************************************************************
-//
-// $Id: G4UrbanMscModel2.hh,v 1.17 2009/05/15 09:26:42 urban Exp $
-// GEANT4 tag $Name: geant4-09-03 $
-//
-// -------------------------------------------------------------------
-//
-//
-// GEANT4 Class header file
-//
-//
-// File name:     G4UrbanMscModel2
-//
-// Author:        Laszlo Urban
-//
-// Creation date: 06.03.2008
-//
-// Modifications:
-// 23.04.2009 L.Urban updated parameterization in UpdateCache method
-//
-//
-// Class Description:
-//
-// Implementation of the model of multiple scattering based on
-// H.W.Lewis Phys Rev 78 (1950) 526 and L.Urban model
-
-// -------------------------------------------------------------------
-//
-
-#ifndef G4UrbanMscModel2_h
-#define G4UrbanMscModel2_h 1
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-#include "G4VMscModel.hh"
-#include "G4PhysicsTable.hh"
-#include "G4MscStepLimitType.hh"
-
-class G4ParticleChangeForMSC;
-class G4SafetyHelper;
-class G4LossTableManager;
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-class G4UrbanMscModel2 : public G4VMscModel
-{
-
-public:
-
-  G4UrbanMscModel2(const G4String& nam = "UrbanMscUni2");
-
-  virtual ~G4UrbanMscModel2();
-
-  void Initialise(const G4ParticleDefinition*, const G4DataVector&);
-
-  G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition* particle,
-				      G4double KineticEnergy,
-				      G4double AtomicNumber,
-				      G4double AtomicWeight=0., 
-				      G4double cut =0.,
-				      G4double emax=DBL_MAX);
-
-  void SampleScattering(const G4DynamicParticle*,
-			G4double safety);
-
-  G4double ComputeTruePathLengthLimit(const G4Track& track,
-				      G4PhysicsTable* theLambdaTable,
-				      G4double currentMinimalStep);
-
-  G4double ComputeGeomPathLength(G4double truePathLength);
-
-  G4double ComputeTrueStepLength(G4double geomStepLength);
-
-  G4double ComputeTheta0(G4double truePathLength,
-                         G4double KineticEnergy);
-
-private:
-
-  G4double SimpleScattering(G4double xmeanth, G4double x2meanth);
-
-  G4double SampleCosineTheta(G4double trueStepLength, G4double KineticEnergy);
-
-  G4double SampleDisplacement();
-
-  G4double LatCorrelation();
-
-  inline G4double GetLambda(G4double kinEnergy);
-
-  inline void SetParticle(const G4ParticleDefinition*);
-
-  inline void UpdateCache();
-
-  //  hide assignment operator
-  G4UrbanMscModel2 & operator=(const  G4UrbanMscModel2 &right);
-  G4UrbanMscModel2(const  G4UrbanMscModel2&);
-
-  const G4ParticleDefinition* particle;
-  G4ParticleChangeForMSC*     fParticleChange;
-
-  G4PhysicsTable*             theLambdaTable;
-  const G4MaterialCutsCouple* couple;
-  G4LossTableManager*         theManager;
-
-  G4double mass;
-  G4double charge,ChargeSquare;
-  G4double masslimite,lambdalimit,fr;
-
-  G4double taubig;
-  G4double tausmall;
-  G4double taulim;
-  G4double currentTau;
-  G4double tlimit;
-  G4double tlimitmin;
-  G4double tlimitminfix;
-  G4double tgeom;
-
-  G4double geombig;
-  G4double geommin;
-  G4double geomlimit;
-  G4double skindepth;
-  G4double smallstep;
-
-  G4double presafety;
-
-  G4double lambda0;
-  G4double lambdaeff;
-  G4double tPathLength;
-  G4double zPathLength;
-  G4double par1,par2,par3;
-
-  G4double stepmin;
-
-  G4double currentKinEnergy;
-  G4double currentRange; 
-  G4double rangeinit;
-  G4double currentRadLength;
-
-  G4double theta0max,rellossmax;
-  G4double third;
-
-  G4int    currentMaterialIndex;
-
-  G4double y;
-  G4double Zold;
-  G4double Zeff,Z2,Z23,lnZ;
-  G4double coeffth1,coeffth2;
-  G4double coeffc1,coeffc2;
-  G4double scr1ini,scr2ini,scr1,scr2;
-
-  G4bool   isInitialized;
-  G4bool   inside;
-  G4bool   insideskin;
-};
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-G4double G4UrbanMscModel2::GetLambda(G4double e)
-{
-  G4double x;
-  if(theLambdaTable) {
-    G4bool b;
-    x = ((*theLambdaTable)[currentMaterialIndex])->GetValue(e, b);
-  } else {
-    x = CrossSection(couple,particle,e);
-  }
-  if(x > DBL_MIN) x = 1./x;
-  else            x = DBL_MAX;
-  return x;
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-void G4UrbanMscModel2::SetParticle(const G4ParticleDefinition* p)
-{
-  if (p != particle) {
-    particle = p;
-    mass = p->GetPDGMass();
-    charge = p->GetPDGCharge()/eplus;
-    ChargeSquare = charge*charge;
-  }
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
-
-inline
-void G4UrbanMscModel2::UpdateCache()                                   
-{
-    lnZ = std::log(Zeff);
-    //new correction in theta0 formula
-    coeffth1 = (1.-8.7780e-2/Zeff)*(0.87+0.03*lnZ);                   
-    coeffth2 = (4.0780e-2+1.7315e-4*Zeff)*(0.87+0.03*lnZ);              
-    // tail parameters
-    G4double lnZ1 = std::log(Zeff+1.);
-    coeffc1  = 2.943-0.197*lnZ1;                  
-    coeffc2  = 0.0987-0.0143*lnZ1;                              
-    // for single scattering
-    Z2 = Zeff*Zeff;
-    Z23 = std::exp(2.*lnZ/3.);
-    scr1     = scr1ini*Z23;
-    scr2     = scr2ini*Z2*ChargeSquare;
-                                              
-    Zold = Zeff;
-}
-
-#endif
-
