Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/splitting/src/G4ElectronCapture.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/extended/medical/dna/splitting/src/G4ElectronCapture.cc (Version 11.3.0) and /examples/extended/medical/dna/splitting/src/G4ElectronCapture.cc (Version 11.0)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 // This example is provided by the Geant4-DNA      26 // This example is provided by the Geant4-DNA collaboration
 27 // Any report or published results obtained us <<  27 // Any report or published results obtained using the Geant4-DNA software 
 28 // shall cite the following Geant4-DNA collabo     28 // shall cite the following Geant4-DNA collaboration publication:
 29 // Med. Phys. 37 (2010) 4692-4708                  29 // Med. Phys. 37 (2010) 4692-4708
 30 // The Geant4-DNA web site is available at htt     30 // The Geant4-DNA web site is available at http://geant4-dna.org
 31 //                                                 31 //
 32 //--------------------------------------------     32 //---------------------------------------------------------------------------
 33 //                                                 33 //
 34 // ClassName:   G4ElectronCapture                  34 // ClassName:   G4ElectronCapture
 35 //                                                 35 //
 36 // Description: The process to kill particles      36 // Description: The process to kill particles to save CPU
 37 //                                                 37 //
 38 // Author:      V.Ivanchenko 31 August 2010        38 // Author:      V.Ivanchenko 31 August 2010
 39 //                                                 39 //
 40 //--------------------------------------------     40 //----------------------------------------------------------------------------
 41 //                                                 41 //
 42 /// \file G4ElectronCapture.cc                     42 /// \file G4ElectronCapture.cc
 43 /// \brief Implementation of the G4ElectronCap     43 /// \brief Implementation of the G4ElectronCapture class
 44                                                    44 
 45 #include "G4ElectronCapture.hh"                    45 #include "G4ElectronCapture.hh"
 46                                                <<  46 #include "G4SystemOfUnits.hh"
 47 #include "G4Electron.hh"                       << 
 48 #include "G4ParticleDefinition.hh"                 47 #include "G4ParticleDefinition.hh"
 49 #include "G4Region.hh"                         << 
 50 #include "G4RegionStore.hh"                    << 
 51 #include "G4Step.hh"                               48 #include "G4Step.hh"
 52 #include "G4SystemOfUnits.hh"                  << 
 53 #include "G4Track.hh"                              49 #include "G4Track.hh"
                                                   >>  50 #include "G4Region.hh"
                                                   >>  51 #include "G4RegionStore.hh"
                                                   >>  52 #include "G4Electron.hh"
 54                                                    53 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 56                                                    55 
 57 G4ElectronCapture::G4ElectronCapture(const G4S     56 G4ElectronCapture::G4ElectronCapture(const G4String& regName, G4double ekinlim)
 58   : G4VDiscreteProcess("eCapture", fElectromag <<  57 : G4VDiscreteProcess("eCapture", fElectromagnetic),
 59     fKinEnergyThreshold(ekinlim),              <<  58   fKinEnergyThreshold(ekinlim),
 60     fRegionName(regName),                      <<  59   fRegionName(regName), fpRegion(0)
 61     fpRegion(0)                                << 
 62 {                                                  60 {
 63   if (regName == "" || regName == "world") {   <<  61   if(regName == "" || regName == "world") { 
 64     fRegionName = "DefaultRegionForTheWorld";      62     fRegionName = "DefaultRegionForTheWorld";
 65   }                                                63   }
 66   pParticleChange = &fParticleChange;              64   pParticleChange = &fParticleChange;
 67 }                                                  65 }
 68                                                    66 
 69 //....oooOO0OOooo........oooOO0OOooo........oo     67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 70                                                    68 
 71 G4ElectronCapture::~G4ElectronCapture() {}     <<  69 G4ElectronCapture::~G4ElectronCapture() 
                                                   >>  70 {}
 72                                                    71 
 73 //....oooOO0OOooo........oooOO0OOooo........oo     72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 74                                                    73 
 75 void G4ElectronCapture::SetKinEnergyLimit(G4do     74 void G4ElectronCapture::SetKinEnergyLimit(G4double val)
 76 {                                                  75 {
 77   fKinEnergyThreshold = val;                       76   fKinEnergyThreshold = val;
 78   if (verboseLevel > 0) {                      <<  77   if(verboseLevel > 0) {
 79     G4cout << "### G4ElectronCapture: Tracking <<  78     G4cout << "### G4ElectronCapture: Tracking cut E(MeV) = " 
 80            << G4endl;                          <<  79         << fKinEnergyThreshold/MeV << G4endl;
 81   }                                                80   }
 82 }                                                  81 }
 83                                                    82 
 84 //....oooOO0OOooo........oooOO0OOooo........oo     83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 85                                                    84 
 86 void G4ElectronCapture::BuildPhysicsTable(cons     85 void G4ElectronCapture::BuildPhysicsTable(const G4ParticleDefinition&)
 87 {                                                  86 {
 88   fpRegion = (G4RegionStore::GetInstance())->G     87   fpRegion = (G4RegionStore::GetInstance())->GetRegion(fRegionName);
 89   if (fpRegion && verboseLevel > 0) {          <<  88   if(fpRegion && verboseLevel > 0) {
 90     G4cout << "### G4ElectronCapture: Tracking <<  89     G4cout << "### G4ElectronCapture: Tracking cut E(MeV) = " 
 91            << " is assigned to " << fRegionNam <<  90         << fKinEnergyThreshold/MeV << " is assigned to " << fRegionName
                                                   >>  91         << G4endl;
 92   }                                                92   }
 93 }                                                  93 }
 94                                                    94 
 95 //....oooOO0OOooo........oooOO0OOooo........oo     95 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 96                                                    96 
 97 G4bool G4ElectronCapture::IsApplicable(const G     97 G4bool G4ElectronCapture::IsApplicable(const G4ParticleDefinition&)
 98 {                                                  98 {
 99   return true;                                     99   return true;
100 }                                                 100 }
101                                                   101 
102 //....oooOO0OOooo........oooOO0OOooo........oo    102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
103                                                   103 
104 G4double G4ElectronCapture::PostStepGetPhysica << 104 G4double 
105                                                << 105 G4ElectronCapture::PostStepGetPhysicalInteractionLength(const G4Track& aTrack,
                                                   >> 106     G4double,
                                                   >> 107     G4ForceCondition* condition)
106 {                                                 108 {
107   // condition is set to "Not Forced"             109   // condition is set to "Not Forced"
108   *condition = NotForced;                         110   *condition = NotForced;
109                                                   111 
110   G4double limit = DBL_MAX;                    << 112   G4double limit = DBL_MAX; 
111   if (fpRegion) {                              << 113   if(fpRegion) {
112     if (aTrack.GetVolume()->GetLogicalVolume() << 114     if(aTrack.GetVolume()->GetLogicalVolume()->GetRegion() == fpRegion && 
113         && aTrack.GetKineticEnergy() < fKinEne << 115         aTrack.GetKineticEnergy() < fKinEnergyThreshold) { limit = 0.0; }
114     {                                          << 
115       limit = 0.0;                             << 
116     }                                          << 
117   }                                               116   }
118   return limit;                                   117   return limit;
119 }                                                 118 }
120                                                   119 
121 //....oooOO0OOooo........oooOO0OOooo........oo    120 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
122                                                   121 
123 G4VParticleChange* G4ElectronCapture::PostStep << 122 G4VParticleChange* G4ElectronCapture::PostStepDoIt(const G4Track& aTrack, 
                                                   >> 123     const G4Step&)
124 {                                                 124 {
125   pParticleChange->Initialize(aTrack);            125   pParticleChange->Initialize(aTrack);
126   pParticleChange->ProposeTrackStatus(fStopAnd    126   pParticleChange->ProposeTrackStatus(fStopAndKill);
127   pParticleChange->ProposeLocalEnergyDeposit(a    127   pParticleChange->ProposeLocalEnergyDeposit(aTrack.GetKineticEnergy());
128   fParticleChange.SetProposedKineticEnergy(0.0    128   fParticleChange.SetProposedKineticEnergy(0.0);
129   return pParticleChange;                         129   return pParticleChange;
130 }                                                 130 }
131                                                   131 
132 //....oooOO0OOooo........oooOO0OOooo........oo    132 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
133                                                   133 
134 G4double G4ElectronCapture::GetMeanFreePath(co << 134 G4double G4ElectronCapture::GetMeanFreePath(const G4Track&,G4double,
                                                   >> 135     G4ForceCondition*)
135 {                                                 136 {
136   return DBL_MAX;                                 137   return DBL_MAX;
137 }                                                 138 }
                                                   >> 139 
138                                                   140