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

Timestamp:

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

Author:

garnier

Message:

update

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.12 2009/01/29 08:13:34 pandola Exp $
-// GEANT4 tag $Name: geant4-09-02-ref-02 $
+// $Id: G4eIonisationCrossSectionHandler.cc,v 1.11 2006/06/29 19:42:00 gunter Exp $
+// GEANT4 tag $Name: geant4-09-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
 //
 // -------------------------------------------------------------------
@@ -100 +98 @@
     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];
@@ -132 +129 @@
             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
@@ -143 +139 @@
                      << " Z= " << Z
                      << G4endl;
-              }
+            }
 
           }
@@ -159 +155 @@
 }
 
-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;
-}
+