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

Timestamp:

Apr 20, 2009, 4:53:50 PM (15 years ago)

Author:

garnier

Message:

fichier oublies

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.5 2007/11/13 18:36:29 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-02 $
+// $Id: G4mplIonisationModel.cc,v 1.6 2009/02/20 16:38:33 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-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();
@@ -176 +177 @@
   return dedx;
 }
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+void G4mplIonisationModel::SampleSecondaries(std::vector<G4DynamicParticle*>*,
+                                             const G4MaterialCutsCouple*,
+                                             const G4DynamicParticle*,
+                                             G4double,
+                                             G4double)
+{}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -204 +214 @@
   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....