Changeset 1055 for trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc

Timestamp:

May 28, 2009, 4:26:57 PM (15 years ago)

Author:

garnier

Message:

maj sur la beta de geant 4.9.3

File:

• trunk/source/processes/electromagnetic/highenergy/src/G4hhIonisation.cc (modified) (2 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4hhIonisation.cc,v 1.8 2008/10/16 14:29:48 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-02 $
+// $Id: G4hhIonisation.cc,v 1.9 2009/02/20 16:38:33 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
 //
 // -------------------------------------------------------------------
@@ -75 +75 @@
 G4hhIonisation::~G4hhIonisation()
 {}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4bool G4hhIonisation::IsApplicable(const G4ParticleDefinition& p)
+{
+  return (p.GetPDGCharge() != 0.0 && p.GetPDGMass() > 100.0*MeV &&
+         !p.IsShortLived());
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4hhIonisation::MinPrimaryEnergy(const G4ParticleDefinition*,
+                                          const G4Material*,
+                                          G4double cut)
+{
+  G4double x = 0.5*cut/electron_mass_c2;
+  G4double y = electron_mass_c2/mass;
+  G4double g = x*y + std::sqrt((1. + x)*(1. + x*y*y));
+  return mass*(g - 1.0);
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....