Changeset 991 for trunk/source/processes/electromagnetic/highenergy/src/G4mplIonisationModel.cc

Timestamp:

Apr 17, 2009, 12:17:14 PM (15 years ago)

Author:

garnier

Message:

update

File:

• trunk/source/processes/electromagnetic/highenergy/src/G4mplIonisationModel.cc (modified) (4 diffs)

trunk/source/processes/electromagnetic/highenergy/src/G4mplIonisationModel.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4mplIonisationModel.cc,v 1.6 2009/02/20 16:38:33 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-02-ref-02 $
+// $Id: G4mplIonisationModel.cc,v 1.5 2007/11/13 18:36:29 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02 $
 //
 // -------------------------------------------------------------------
@@ -139 +139 @@
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
 
-G4double G4mplIonisationModel::ComputeDEDXAhlen(const G4Material* material, 
-                                                G4double bg2)
+G4double G4mplIonisationModel::ComputeDEDXAhlen(const G4Material* material, G4double bg2)
 {
   G4double eDensity = material->GetElectronDensity();
@@ -177 +176 @@
   return dedx;
 }
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-void G4mplIonisationModel::SampleSecondaries(std::vector<G4DynamicParticle*>*,
-                                             const G4MaterialCutsCouple*,
-                                             const G4DynamicParticle*,
-                                             G4double,
-                                             G4double)
-{}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -214 +204 @@
   return loss;
 }
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
-
-G4double G4mplIonisationModel::Dispersion(const G4Material* material,
-                                          const G4DynamicParticle* dp,
-                                          G4double& tmax,
-                                          G4double& length)
-{
-  G4double siga = 0.0;
-  G4double tau   = dp->GetKineticEnergy()/mass;
-  if(tau > 0.0) { 
-    G4double electronDensity = material->GetElectronDensity();
-    G4double gam   = tau + 1.0;
-    G4double invbeta2 = (gam*gam)/(tau * (tau+2.0));
-    siga  = (invbeta2 - 0.5) * twopi_mc2_rcl2 * tmax * length
-      * electronDensity * chargeSquare;
-  }
-  return siga;
-}
-
-//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....