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

Timestamp:

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

Author:

garnier

Message:

update processes

File:

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

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4PEEffectModel.cc,v 1.6 2007/05/22 17:34:36 vnivanch Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4PEEffectModel.cc,v 1.7 2009/02/20 12:06:37 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -41 +41 @@
 //
 // 04.12.05 : SetProposedKineticEnergy(0.) for the killed photon (mma)
+// 20.02.09 : Added initialisation of deexcitation flag and method
+//            CrossSectionPerVolume instead of mfp (V.Ivanchenko)
 //
 // Class Description:
@@ -66 +68 @@
   theGamma    = G4Gamma::Gamma();
   theElectron = G4Electron::Electron();
+  fminimalEnergy = 1.0*eV;
 }
 
@@ -71 +74 @@
 
 G4PEEffectModel::~G4PEEffectModel()
-{
-}
+{}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
@@ -79 +81 @@
                                  const G4DataVector&)
 {
- if (isInitialized) return;
- if (pParticleChange)
-   fParticleChange =
-                  reinterpret_cast<G4ParticleChangeForGamma*>(pParticleChange);
-  else
-   fParticleChange = new G4ParticleChangeForGamma();
-
- fminimalEnergy = 1.0*eV;
+  // always false before the run
+  SetDeexcitationFlag(false);
+
+  if (isInitialized) return;
+  if (pParticleChange) {
+    fParticleChange = reinterpret_cast<G4ParticleChangeForGamma*>(pParticleChange);
+  } else {
+    fParticleChange = new G4ParticleChangeForGamma();
+  }
+  isInitialized = true;
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
+
+G4double G4PEEffectModel::ComputeCrossSectionPerAtom(const G4ParticleDefinition*,
+                                                     G4double energy,
+                                                     G4double Z, G4double,
+                                                     G4double, G4double)
+{
+  G4double* SandiaCof = G4SandiaTable::GetSandiaCofPerAtom((G4int)Z, energy);
+
+  G4double energy2 = energy*energy;
+  G4double energy3 = energy*energy2;
+  G4double energy4 = energy2*energy2;
+
+  return SandiaCof[0]/energy  + SandiaCof[1]/energy2 +
+    SandiaCof[2]/energy3 + SandiaCof[3]/energy4;
+}
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
+
+G4double G4PEEffectModel::CrossSectionPerVolume(const G4Material* material,
+                                                const G4ParticleDefinition*,
+                                                G4double energy,
+                                                G4double, G4double)
+{
+  G4double* SandiaCof = 
+    material->GetSandiaTable()->GetSandiaCofForMaterial(energy);
+                                
+  G4double energy2 = energy*energy;
+  G4double energy3 = energy*energy2;
+  G4double energy4 = energy2*energy2;
+          
+  return SandiaCof[0]/energy  + SandiaCof[1]/energy2 +
+    SandiaCof[2]/energy3 + SandiaCof[3]/energy4; 
 }