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

Timestamp:

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

Author:

garnier

Message:

update processes

File:

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

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4BetheHeitlerModel.cc,v 1.11 2007/05/22 17:34:36 vnivanch Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4BetheHeitlerModel.cc,v 1.12 2008/10/15 15:54:57 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -61 +61 @@
 #include "G4PhysicsLogVector.hh"
 #include "G4ParticleChangeForGamma.hh"
+#include "G4LossTableManager.hh"
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -70 +71 @@
   : G4VEmModel(nam),
     theCrossSectionTable(0),
-    nbins(200)
-{
+    nbins(10)
+{
+  fParticleChange = 0;
   theGamma    = G4Gamma::Gamma();
   thePositron = G4Positron::Positron();
@@ -92 +94 @@
                                      const G4DataVector&)
 {
-  if(pParticleChange)
-    fParticleChange = reinterpret_cast<G4ParticleChangeForGamma*>(pParticleChange);
-  else
-    fParticleChange = new G4ParticleChangeForGamma();
+  if(!fParticleChange) {
+    if(pParticleChange) {
+      fParticleChange = reinterpret_cast<G4ParticleChangeForGamma*>(pParticleChange);
+    } else {
+      fParticleChange = new G4ParticleChangeForGamma();
+    }
+  }
 
   if(theCrossSectionTable) {
@@ -108 +113 @@
   G4double emin = LowEnergyLimit();
   G4double emax = HighEnergyLimit();
+  G4int n = nbins*G4int(log10(emax/emin));
+  G4bool spline = G4LossTableManager::Instance()->SplineFlag(); 
   G4double e, value;
 
   for(size_t j=0; j<nvect ; j++) { 
 
-    ptrVector  = new G4PhysicsLogVector(emin, emax, nbins);
+    ptrVector  = new G4PhysicsLogVector(emin, emax, n);
+    ptrVector->SetSpline(spline);
     G4double Z = (*theElementTable)[j]->GetZ();
     G4int   iz = G4int(Z);