| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
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. 37 (2010) 4692-4708
30 // J. Comput. Phys. 274 (2014) 841-882
31 // The Geant4-DNA web site is available at htt
32 //
33 //
34 /// \file SteppingAction.cc
35 /// \brief Implementation of the SteppingActio
36
37 #include "SteppingAction.hh"
38
39 #include "EventAction.hh"
40
41 #include "G4SteppingManager.hh"
42
43 //....oooOO0OOooo........oooOO0OOooo........oo
44 SteppingAction::SteppingAction(EventAction* ev
45
46 //....oooOO0OOooo........oooOO0OOooo........oo
47 void SteppingAction::UserSteppingAction(const
48 {
49 // /////////////////////////////////////////
50 // killing primary particles (any heavy char
51 // collimator
52 auto part_name = step->GetTrack()->GetPartic
53 if (step->GetPreStepPoint()->GetPhysicalVolu
54 && step->GetTrack()->GetParticleDefiniti
55 {
56 step->GetTrack()->SetTrackStatus(fStopAndK
57 fEventAction->InvalidEvent();
58 return; // nothing more to do in such cas
59 }
60
61 // /////////////////////////////////////////
62 // recording energy of the primary particle
63 // 1. Boundary between World and Mylar wall,
64 // particle.
65 // 2. Boundary between World and silicon det
66 // particle.
67 // 3. Boundary between wall inner layer (inn
68 // interaction energy of the particle.
69
70 if (step->GetTrack()->GetParticleDefinition(
71 // 1. and 2. - particle crossing the bound
72 if (step->GetPreStepPoint()->GetPhysicalVo
73 auto physVol = step->GetPostStepPoint()-
74 // the above may be a null pointer as th
75 // has to be tested:
76 if (physVol) {
77 // 1. to Mylar wall
78 if (physVol->GetName() == "wallPhys")
79 fEventAction->RecordInitialEnergy(st
80 }
81 // 2. to silicon detector
82 if (physVol->GetName() == "enDetPhys")
83 fEventAction->RecordFinalEnergy(step
84 }
85 }
86 }
87
88 // 3. - particle crossing the boundary fro
89 if (step->GetPreStepPoint()->GetPhysicalVo
90 if (step->GetPostStepPoint()->GetPhysica
91 fEventAction->RecordInteractionEnergy(
92 }
93 }
94 }
95
96 // ////////////////////////////////////////
97 // counting individual ionization acts
98 // only ionisations in the target will be re
99 if (step->GetPreStepPoint()->GetPhysicalVolu
100 if (step->GetTotalEnergyDeposit() != 0) {
101 G4String process_name = step->GetPostSte
102
103 // counting ionisations caused only by e
104 // if(process_name=="e-_G4DNA
105 // process_name=="alpha_G4DNA
106 // fEventAction->AddIonis
107
108 // more general approach - counting ioni
109 // (including also protons, GenericIons,
110 auto pos =
111 process_name.find("Ionisation"); // l
112 if (pos != std::string::npos) // found!
113 fEventAction->AddIonisation(); // inc
114 }
115 }
116 }
117 //....oooOO0OOooo........oooOO0OOooo........oo
118