Changeset 961 for trunk/source/processes/electromagnetic/lowenergy/src/G4eIonisationCrossSectionHandler.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4eIonisationCrossSectionHandler.cc (modified) (6 diffs)

trunk/source/processes/electromagnetic/lowenergy/src/G4eIonisationCrossSectionHandler.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4eIonisationCrossSectionHandler.cc,v 1.11 2006/06/29 19:42:00 gunter Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4eIonisationCrossSectionHandler.cc,v 1.12 2009/01/29 08:13:34 pandola Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -42 +42 @@
 // 19 Jul 2002   VI          Create composite data set for material
 // 21 Jan 2003  V.Ivanchenko Cut per region
+// 28 Jan 2009  L.Pandola    Added public method to make a easier migration of 
+//                           G4LowEnergyIonisation to G4LivermoreIonisationModel
 //
 // -------------------------------------------------------------------
@@ -98 +100 @@
     G4int nElements = material->GetNumberOfElements();
 
-    if(verbose > 0) {
-      G4cout << "eIonisation CS for " << m << "th material "
-             << material->GetName()
-             << "  eEl= " << nElements << G4endl;
-    }
+    if(verbose > 0) 
+      {
+        G4cout << "eIonisation CS for " << m << "th material "
+               << material->GetName()
+               << "  eEl= " << nElements << G4endl;
+      }
 
     G4double tcut  = (*energyCuts)[m];
@@ -129 +132 @@
             value += cross * p * density;
 
-            if(verbose>0 && m == 0 && e>=1. && e<=0.) {
+            if(verbose>0 && m == 0 && e>=1. && e<=0.) 
+            {
               G4cout << "G4eIonCrossSH: e(MeV)= " << e/MeV
                      << " n= " << n
@@ -139 +143 @@
                      << " Z= " << Z
                      << G4endl;
-            }
+              }
 
           }
@@ -155 +159 @@
 }
 
-
+G4double G4eIonisationCrossSectionHandler::GetCrossSectionAboveThresholdForElement(G4double energy,
+                                                                                   G4double cutEnergy,
+                                                                                   G4int Z)
+{
+  G4int nShells = NumberOfComponents(Z);
+  G4double value = 0.;
+  if(energy > cutEnergy) 
+    {
+      for (G4int n=0; n<nShells; n++) {
+        G4double cross = FindValue(Z, energy, n);
+        G4double p = theParam->Probability(Z, cutEnergy, energy, energy, n);
+        value += cross * p;
+      }
+    }
+  return value;
+}