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

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4FinalStateChargeDecrease.cc (modified) (12 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-//
-// $Id: G4FinalStateChargeDecrease.cc,v 1.2 2007/11/09 20:11:04 pia Exp $
-// GEANT4 tag $Name:  $
-// 
-// Contact Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
-//
-// Reference: TNS Geant4-DNA paper
-// Reference for implementation model: NIM. 155, pp. 145-156, 1978
-
-// History:
-// -----------
-// Date         Name              Modification
-// 28 Apr 2007  M.G. Pia          Created in compliance with design described in TNS paper
-//
-// -------------------------------------------------------------------
-
-// Class description:
-// Reference: TNS Geant4-DNA paper
-// S. Chauvie et al., Geant4 physics processes for microdosimetry simulation:
-// design foundation and implementation of the first set of models,
-// IEEE Trans. Nucl. Sci., vol. 54, no. 6, Dec. 2007.
-// Further documentation available from http://www.ge.infn.it/geant4/dna
-
-// -------------------------------------------------------------------
-
+// $Id: G4FinalStateChargeDecrease.cc,v 1.4 2009/01/20 07:50:28 sincerti Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 
 #include "G4FinalStateChargeDecrease.hh"
-#include "G4Track.hh"
-#include "G4Step.hh"
-#include "G4DynamicParticle.hh"
-//#include "Randomize.hh"
 
-#include "G4ParticleTypes.hh"
-#include "G4ParticleDefinition.hh"
-#include "G4SystemOfUnits.hh"
-//#include "G4ParticleMomentum.hh"
-#include "G4DNAGenericIonsManager.hh"
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4FinalStateChargeDecrease::G4FinalStateChargeDecrease()
 {
-  name = "ChargeDecrease";
   lowEnergyLimit = 1 * keV;
   highEnergyLimit = 10 * MeV;
 }
 
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4FinalStateChargeDecrease::~G4FinalStateChargeDecrease()
 {}
  
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 const G4FinalStateProduct& G4FinalStateChargeDecrease::GenerateFinalState(const G4Track& track, const G4Step& /* step */)
 {
-  // Clear previous secondaries, energy deposit and particle kill status
   product.Clear();
 
@@ -95 +64 @@
   
   if (outK<0)
-    {
-      G4String message;
-      message="ChargeDecreaseDingfelder::GenerateFinalState - Final kinetic energy is below 0! Process ";
-      G4Exception(message);
-    }
+  {
+    G4String message;
+    message="ChargeDecreaseDingfelder::GenerateFinalState - Final kinetic energy is below 0! Process ";
+    G4Exception(message);
+  }
   
-  // Primary particle
+  //SI - Added protection against total energy deposit
+  product.DoNotDepositEnergy();
+  //
   product.KillPrimaryParticle();
+
   product.AddEnergyDeposit(waterBindingEnergy);
-  
-  //Secondary particle
+
   G4DynamicParticle* aSecondary = new G4DynamicParticle(OutgoingParticleDefinition(definition, finalStateIndex), 
                                                         track.GetDynamicParticle()->GetMomentumDirection(), 
@@ -114 +85 @@
   return product;
 }
+
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4int G4FinalStateChargeDecrease::NumberOfFinalStates(G4ParticleDefinition* particleDefinition, 
@@ -125 +98 @@
   
   if (particleDefinition == instance->GetIon("alpha++") ) 
-    {
-      if (finalStateIndex==0)  return 1;
-      return 2;   
-    }
+  {
+    if (finalStateIndex==0)  return 1;
+    return 2;   
+  }
   
   if (particleDefinition == instance->GetIon("alpha+") ) return 1;
@@ -135 +108 @@
 }
 
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4ParticleDefinition* G4FinalStateChargeDecrease::OutgoingParticleDefinition (G4ParticleDefinition* particleDefinition, 
@@ -144 +118 @@
   
   if (particleDefinition == instance->GetIon("alpha++") ) 
-    {
-      if (finalStateIndex == 0) return instance->GetIon("alpha+"); 
-
-      return instance->GetIon("helium");    
-    }
+  {
+    if (finalStateIndex == 0) return instance->GetIon("alpha+"); 
+    return instance->GetIon("helium");    
+  }
   
   if (particleDefinition == instance->GetIon("alpha+") ) return instance->GetIon("helium");    
@@ -155 +128 @@
 }
 
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4double G4FinalStateChargeDecrease::WaterBindingEnergyConstant(G4ParticleDefinition* particleDefinition, 
@@ -168 +142 @@
 
   if (particleDefinition == instance->GetIon("alpha++") ) 
-    {
+  {
       // Binding energy for    W+ -> W++ + e-    10.79 eV
       // Binding energy for    W  -> W+  + e-    10.79 eV
@@ -178 +152 @@
   
       return 10.79*2*eV;
-    }
+  }
 
   if (particleDefinition == instance->GetIon("alpha+") ) 
-    {
+  {
       // Binding energy for    W+ -> W++ + e-    10.79 eV
       // Binding energy for    W  -> W+  + e-    10.79 eV
@@ -189 +163 @@
 
       return 10.79*eV;
-    }
+  }
   
   return 0;
 }
 
+//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 
 G4double G4FinalStateChargeDecrease::OutgoingParticleBindingEnergyConstant(G4ParticleDefinition* particleDefinition, 
@@ -203 +178 @@
 
   if (particleDefinition == instance->GetIon("alpha++") ) 
-    {
+  {
       // Binding energy for    He+ -> He++ + e-    54.509 eV
       // Binding energy for    He  -> He+  + e-    24.587 eV
@@ -210 +185 @@
       
       return (54.509 + 24.587)*eV;
-    }
+  }
   
   if (particleDefinition == instance->GetIon("alpha+") ) 
-    {
+  {
       // Binding energy for    He+ -> He++ + e-    54.509 eV
       // Binding energy for    He  -> He+  + e-    24.587 eV
       
       return 24.587*eV;
-    }  
+  }  
   
   return 0;