Geant4 Cross Reference |
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 //-------------------------------------------- 26 //--------------------------------------------------------------------------- 27 // 27 // 28 // ClassName: G4ElectronCapture 28 // ClassName: G4ElectronCapture 29 // 29 // 30 // Description: The process to kill particles 30 // Description: The process to kill particles to save CPU 31 // 31 // 32 // Author: V.Ivanchenko 31 August 2010 32 // Author: V.Ivanchenko 31 August 2010 33 // 33 // 34 //-------------------------------------------- 34 //---------------------------------------------------------------------------- 35 // 35 // 36 // This example is provided by the Geant4-DNA 36 // This example is provided by the Geant4-DNA collaboration 37 // Any report or published results obtained us << 37 // Any report or published results obtained using the Geant4-DNA software 38 // and the DNA geometry given in the Geom_DNA << 38 // and the DNA geometry given in the Geom_DNA example 39 // shall cite the following Geant4-DNA collabo 39 // shall cite the following Geant4-DNA collaboration publications: 40 // [1] NIM B 298 (2013) 47-54 40 // [1] NIM B 298 (2013) 47-54 41 // [2] Med. Phys. 37 (2010) 4692-4708 41 // [2] Med. Phys. 37 (2010) 4692-4708 42 // The Geant4-DNA web site is available at htt 42 // The Geant4-DNA web site is available at http://geant4-dna.org 43 // 43 // 44 //....oooOO0OOooo........oooOO0OOooo........oo 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 // 45 // 46 /// \file G4ElectronCapture.cc 46 /// \file G4ElectronCapture.cc 47 /// \brief Implementation of the G4ElectronCap 47 /// \brief Implementation of the G4ElectronCapture class 48 48 49 #include "G4ElectronCapture.hh" 49 #include "G4ElectronCapture.hh" 50 << 50 #include "G4SystemOfUnits.hh" 51 #include "G4Electron.hh" << 52 #include "G4ParticleDefinition.hh" 51 #include "G4ParticleDefinition.hh" 53 #include "G4Region.hh" << 54 #include "G4RegionStore.hh" << 55 #include "G4Step.hh" 52 #include "G4Step.hh" 56 #include "G4SystemOfUnits.hh" << 57 #include "G4Track.hh" 53 #include "G4Track.hh" >> 54 #include "G4Region.hh" >> 55 #include "G4RegionStore.hh" >> 56 #include "G4Electron.hh" 58 57 59 //....oooOO0OOooo........oooOO0OOooo........oo 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 60 59 61 G4ElectronCapture::G4ElectronCapture(const G4S 60 G4ElectronCapture::G4ElectronCapture(const G4String& regName, G4double ekinlim) 62 : G4VDiscreteProcess("eCapture", fElectromag << 61 : G4VDiscreteProcess("eCapture", fElectromagnetic), 63 fKinEnergyThreshold(ekinlim), << 62 fKinEnergyThreshold(ekinlim), 64 fRegionName(regName), << 63 fRegionName(regName), 65 region(0) << 64 region(0) 66 { 65 { 67 if (regName == "" || regName == "world") { << 66 if (regName == "" || regName == "world") >> 67 { 68 fRegionName = "DefaultRegionForTheWorld"; 68 fRegionName = "DefaultRegionForTheWorld"; 69 } 69 } 70 pParticleChange = &fParticleChange; 70 pParticleChange = &fParticleChange; 71 } 71 } 72 72 73 //....oooOO0OOooo........oooOO0OOooo........oo 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 74 74 75 G4ElectronCapture::~G4ElectronCapture() {} << 75 G4ElectronCapture::~G4ElectronCapture() >> 76 { >> 77 } 76 78 77 //....oooOO0OOooo........oooOO0OOooo........oo 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 78 80 79 void G4ElectronCapture::SetKinEnergyLimit(G4do << 81 void >> 82 G4ElectronCapture::SetKinEnergyLimit(G4double val) 80 { 83 { 81 fKinEnergyThreshold = val; 84 fKinEnergyThreshold = val; 82 if (verboseLevel > 0) { << 85 if (verboseLevel > 0) 83 G4cout << "### G4ElectronCapture: Tracking << 86 { 84 << G4endl; << 87 G4cout << "### G4ElectronCapture: Tracking cut E(MeV) = " >> 88 << fKinEnergyThreshold / MeV<< G4endl; 85 } 89 } 86 } 90 } 87 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 88 92 89 void G4ElectronCapture::BuildPhysicsTable(cons << 93 void >> 94 G4ElectronCapture::BuildPhysicsTable(const G4ParticleDefinition&) 90 { 95 { 91 region = (G4RegionStore::GetInstance())->Get 96 region = (G4RegionStore::GetInstance())->GetRegion(fRegionName); 92 if (region && verboseLevel > 0) { << 97 if (region && verboseLevel > 0) 93 G4cout << "### G4ElectronCapture: Tracking << 98 { 94 << " is assigned to " << fRegionNam << 99 G4cout << "### G4ElectronCapture: Tracking cut E(MeV) = " >> 100 << fKinEnergyThreshold / MeV<< " is assigned to " << fRegionName >> 101 << G4endl; 95 } 102 } 96 } 103 } 97 104 98 //....oooOO0OOooo........oooOO0OOooo........oo 105 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 99 106 100 G4bool G4ElectronCapture::IsApplicable(const G << 107 G4bool >> 108 G4ElectronCapture::IsApplicable(const G4ParticleDefinition&) 101 { 109 { 102 return true; 110 return true; 103 } 111 } 104 112 105 G4double G4ElectronCapture::PostStepGetPhysica << 113 G4double 106 << 114 G4ElectronCapture::PostStepGetPhysicalInteractionLength(const G4Track& aTrack, >> 115 G4double, >> 116 G4ForceCondition* condition) 107 { 117 { 108 // condition is set to "Not Forced" 118 // condition is set to "Not Forced" 109 *condition = NotForced; 119 *condition = NotForced; 110 120 111 G4double limit = DBL_MAX; 121 G4double limit = DBL_MAX; 112 if (region) { << 122 if (region) >> 123 { 113 if (aTrack.GetVolume()->GetLogicalVolume() 124 if (aTrack.GetVolume()->GetLogicalVolume()->GetRegion() == region 114 && aTrack.GetKineticEnergy() < fKinEne 125 && aTrack.GetKineticEnergy() < fKinEnergyThreshold) 115 { 126 { 116 limit = 0.0; 127 limit = 0.0; 117 } 128 } 118 } 129 } 119 return limit; 130 return limit; 120 } 131 } 121 132 122 //....oooOO0OOooo........oooOO0OOooo........oo 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 123 134 124 G4VParticleChange* G4ElectronCapture::PostStep << 135 G4VParticleChange* >> 136 G4ElectronCapture::PostStepDoIt(const G4Track& aTrack, const G4Step&) 125 { 137 { 126 pParticleChange->Initialize(aTrack); 138 pParticleChange->Initialize(aTrack); 127 pParticleChange->ProposeTrackStatus(fStopAnd 139 pParticleChange->ProposeTrackStatus(fStopAndKill); 128 pParticleChange->ProposeLocalEnergyDeposit(a 140 pParticleChange->ProposeLocalEnergyDeposit(aTrack.GetKineticEnergy()); 129 fParticleChange.SetProposedKineticEnergy(0.0 141 fParticleChange.SetProposedKineticEnergy(0.0); 130 return pParticleChange; 142 return pParticleChange; 131 } 143 } 132 144 133 //....oooOO0OOooo........oooOO0OOooo........oo 145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 146 135 G4double G4ElectronCapture::GetMeanFreePath(co << 147 G4double >> 148 G4ElectronCapture::GetMeanFreePath(const G4Track&, G4double, G4ForceCondition*) 136 { 149 { 137 return DBL_MAX; 150 return DBL_MAX; 138 } 151 } 139 152 140 //....oooOO0OOooo........oooOO0OOooo........oo 153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 154 141 155