Changeset 961 for trunk/source/processes/electromagnetic/standard/src/G4eIonisation.cc

Timestamp:

Apr 6, 2009, 12:21:12 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/standard/src/G4eIonisation.cc (modified) (5 diffs)

trunk/source/processes/electromagnetic/standard/src/G4eIonisation.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4eIonisation.cc,v 1.53 2007/05/22 17:34:36 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4eIonisation.cc,v 1.57 2009/02/20 12:06:37 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -90 +90 @@
     isInitialised(false)
 {
-  SetStepFunction(0.2, 1*mm);
-  SetIntegral(true);
-  SetVerboseLevel(1);
+  //  SetStepFunction(0.2, 1*mm);
+  // SetIntegral(true);
+  // SetVerboseLevel(1);
+  SetProcessSubType(fIonisation);
 }
 
@@ -99 +100 @@
 G4eIonisation::~G4eIonisation()
 {}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4eIonisation::MinPrimaryEnergy(const G4ParticleDefinition*,
+                                         const G4Material*,
+                                         G4double cut)
+{
+  G4double x = cut;
+  if(isElectron) x += cut;
+  return x;
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4bool G4eIonisation::IsApplicable(const G4ParticleDefinition& p)
+{
+  return (&p == G4Electron::Electron() || &p == G4Positron::Positron());
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -110 +129 @@
     SetSecondaryParticle(theElectron);
     if (!EmModel()) SetEmModel(new G4MollerBhabhaModel());
-    EmModel()->SetLowEnergyLimit (100*eV);
-    EmModel()->SetHighEnergyLimit(100*TeV);
+    EmModel()->SetLowEnergyLimit (MinKinEnergy());
+    EmModel()->SetHighEnergyLimit(MaxKinEnergy());
     if (!FluctModel()) SetFluctModel(new G4UniversalFluctuation());
                 
@@ -122 +141 @@
 
 void G4eIonisation::PrintInfo()
-{
-  if(EmModel())
-    G4cout << "      Delta cross sections and sampling from " 
-           << EmModel()->GetName() << " model"
-           << "\n      Good description from 1 KeV to 100 GeV."
-           << G4endl;
-}
+{}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....