| 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 // This example is provided by the Geant4-DNA <<
27 // Any report or published results obtained us <<
28 // shall cite the following Geant4-DNA collabo <<
29 // Med. Phys. 45 (2018) e722-e739 <<
30 // Phys. Med. 31 (2015) 861-874 <<
31 // Med. Phys. 37 (2010) 4692-4708 <<
32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) <<
33 // <<
34 // The Geant4-DNA web site is available at htt <<
35 // <<
36 /// \file SteppingAction.cc 26 /// \file SteppingAction.cc
37 /// \brief Implementation of the SteppingActio 27 /// \brief Implementation of the SteppingAction class
38 28
39 #include "SteppingAction.hh" 29 #include "SteppingAction.hh"
>> 30 #include "Run.hh"
40 31
41 #include "G4DNAGenericIonsManager.hh" << 32 #include "G4RunManager.hh"
42 #include "G4Electron.hh" 33 #include "G4Electron.hh"
43 #include "G4Gamma.hh" 34 #include "G4Gamma.hh"
44 #include "G4Proton.hh" 35 #include "G4Proton.hh"
45 #include "G4RunManager.hh" << 36 #include "G4DNAGenericIonsManager.hh"
46 37
47 //....oooOO0OOooo........oooOO0OOooo........oo 38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
48 39
49 SteppingAction::SteppingAction() : G4UserStepp << 40 SteppingAction::SteppingAction()
>> 41 :G4UserSteppingAction()
50 { 42 {
51 fSumOfStepLength = 0; << 43 fSumOfStepLength = 0;
52 fLength = 0; << 44 fLength = 0;
53 fDeltaE = -7; << 45 fDeltaE = -7;
54 fTotalStoppingPower = -8; 46 fTotalStoppingPower = -8;
55 fTotalNumberOfSteps = 0; 47 fTotalNumberOfSteps = 0;
56 fNumberOfSteps = 0; << 48 fNumberOfSteps = 0;
57 fDepositedEnergy = 0; << 49 fDepositedEnergy = 0;
58 } 50 }
59 51
60 //....oooOO0OOooo........oooOO0OOooo........oo 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
61 53
62 SteppingAction::~SteppingAction() {} << 54 SteppingAction::~SteppingAction()
>> 55 {}
63 56
64 //....oooOO0OOooo........oooOO0OOooo........oo 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65 58
66 void SteppingAction::UserSteppingAction(const 59 void SteppingAction::UserSteppingAction(const G4Step* step)
67 { 60 {
68 //*** WARNING: this line will kill all secon <<
69 61
70 if (step->GetTrack()->GetDefinition() == G4E << 62 //*** WARNING: this line will kill all secondary e- and gammas ***
71 && step->GetTrack()->GetTrackID() != 1) <<
72 step->GetTrack()->SetTrackStatus(fStopAndK <<
73 63
74 if (step->GetTrack()->GetDefinition() == G4G << 64 if ( step->GetTrack()->GetDefinition()==G4Electron::ElectronDefinition()
>> 65 && step->GetTrack()->GetTrackID() != 1 )
75 step->GetTrack()->SetTrackStatus(fStopAndK 66 step->GetTrack()->SetTrackStatus(fStopAndKill);
76 67
77 // << 68 if ( step->GetTrack()->GetDefinition()==G4Gamma::GammaDefinition() )
78 << 69 step->GetTrack()->SetTrackStatus(fStopAndKill);
79 G4double charge = step->GetTrack()->GetDefin <<
80 G4double steplen = step->GetStepLength(); <<
81 G4double edep = step->GetTotalEnergyDeposit( <<
82 G4int maxNumberOfSteps = -1; <<
83 <<
84 G4DNAGenericIonsManager* instance; <<
85 instance = G4DNAGenericIonsManager::Instance <<
86 <<
87 if ( <<
88 // ELECTRONS <<
89 (step->GetTrack()->GetTrackID() == 1 <<
90 && step->GetTrack()->GetParticleDefinitio <<
91 <<
92 || <<
93 <<
94 // PROTONS, HYDROGEN, ALPHA PARTICLES AND <<
95 (charge != -1. && step->GetTrack()->GetPar <<
96 <<
97 ) <<
98 { <<
99 fSumOfStepLength = fSumOfStepLength + step <<
100 fDepositedEnergy = fDepositedEnergy + edep <<
101 fNumberOfSteps = fNumberOfSteps + 1; <<
102 <<
103 // G4cout << fSumOfStepLength << " " << fD <<
104 <<
105 // ELECTRONS <<
106 <<
107 if (step->GetTrack()->GetParticleDefinitio <<
108 maxNumberOfSteps = 1000; <<
109 <<
110 // PROTONS, HYDROGEN <<
111 70
112 if (step->GetTrack()->GetParticleDefinitio << 71 //
113 || step->GetTrack()->GetParticleDefini << 72
>> 73 G4double charge = step->GetTrack()->GetDefinition()->GetPDGCharge();
>> 74 G4double steplen = step->GetStepLength();
>> 75 G4double edep = step->GetTotalEnergyDeposit();
>> 76 G4int maxNumberOfSteps = -1;
>> 77
>> 78 G4DNAGenericIonsManager *instance;
>> 79 instance = G4DNAGenericIonsManager::Instance();
>> 80
>> 81 if (
>> 82
>> 83 // ELECTRONS
>> 84 (
>> 85 step->GetTrack()->GetTrackID() == 1 &&
>> 86 step->GetTrack()->GetParticleDefinition()==G4Electron::ElectronDefinition()
>> 87 )
>> 88
>> 89 ||
>> 90
>> 91 // PROTONS, HYDROGEN, ALPHA PARTICLES AND CHARGED STATES
>> 92 (
>> 93 charge != -1. &&
>> 94 step->GetTrack()->GetParticleDefinition()!=G4Gamma::GammaDefinition()
>> 95 )
>> 96
>> 97 )
>> 98 {
>> 99 fSumOfStepLength = fSumOfStepLength + steplen;
>> 100 fDepositedEnergy = fDepositedEnergy + edep;
>> 101 fNumberOfSteps = fNumberOfSteps + 1;
>> 102
>> 103 //G4cout << fSumOfStepLength << " " << fDepositedEnergy << " " << fNumberOfSteps << G4endl;
>> 104
>> 105 // ELECTRONS
>> 106
>> 107 if (step->GetTrack()->GetParticleDefinition()==G4Electron::ElectronDefinition())
114 maxNumberOfSteps = 1000; 108 maxNumberOfSteps = 1000;
115 109
116 // He0, He+, He2+ << 110 // PROTONS, HYDROGEN
>> 111
>> 112 if
>> 113 (
>> 114 step->GetTrack()->GetParticleDefinition()==G4Proton::ProtonDefinition()
>> 115 ||
>> 116 step->GetTrack()->GetParticleDefinition()==instance->GetIon("hydrogen")
>> 117 )
>> 118 maxNumberOfSteps = 1000;
>> 119
>> 120 // He0, He+, He2+
>> 121
>> 122 if
>> 123 (
>> 124 step->GetTrack()->GetParticleDefinition()==instance->GetIon("alpha++")
>> 125 ||
>> 126 step->GetTrack()->GetParticleDefinition()==instance->GetIon("alpha+")
>> 127 ||
>> 128 step->GetTrack()->GetParticleDefinition()==instance->GetIon("helium")
>> 129 )
>> 130 maxNumberOfSteps = 10000;
>> 131
>> 132
>> 133 // All
>> 134
>> 135 if (fNumberOfSteps > maxNumberOfSteps)
>> 136 {
>> 137 fLength = fSumOfStepLength;
>> 138
>> 139 fTotalStoppingPower = fDepositedEnergy/fLength;
>> 140 //G4cout << fTotalStoppingPower/(MeV/cm) << G4endl;
>> 141
>> 142 G4RunManager::GetRunManager()->AbortEvent();
>> 143 fSumOfStepLength=0.;
>> 144 fDepositedEnergy=0;
>> 145 fNumberOfSteps = 0;
>> 146
>> 147 };
>> 148 };
117 149
118 if (step->GetTrack()->GetParticleDefinitio <<
119 || step->GetTrack()->GetParticleDefini <<
120 || step->GetTrack()->GetParticleDefini <<
121 maxNumberOfSteps = 10000; <<
122 <<
123 // IONS (above helium) <<
124 if (charge > 2) maxNumberOfSteps = 10000; <<
125 <<
126 // ALL <<
127 <<
128 if (fNumberOfSteps > maxNumberOfSteps) { <<
129 fLength = fSumOfStepLength; <<
130 <<
131 fTotalStoppingPower = fDepositedEnergy / <<
132 // G4cout << fTotalStoppingPower/(MeV/cm <<
133 <<
134 G4RunManager::GetRunManager()->AbortEven <<
135 fSumOfStepLength = 0.; <<
136 fDepositedEnergy = 0; <<
137 fNumberOfSteps = 0; <<
138 }; <<
139 }; <<
140 } 150 }
141 151