| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // This example is provided by the Geant4-DNA collaboration
27 // Any report or published results obtained using the Geant4-DNA software
28 // shall cite the following Geant4-DNA collaboration publication:
29 // Med. Phys. 37 (2010) 4692-4708
30 // J. Comput. Phys. 274 (2014) 841-882
31 // The Geant4-DNA web site is available at http://geant4-dna.org
32 //
33 //
34 /// \file SteppingAction.cc
35 /// \brief Implementation of the SteppingAction class
36
37 #include "SteppingAction.hh"
38
39 #include "EventAction.hh"
40
41 #include "G4SteppingManager.hh"
42
43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 SteppingAction::SteppingAction(EventAction* event) : G4UserSteppingAction(), fEventAction(event) {}
45
46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 void SteppingAction::UserSteppingAction(const G4Step* step)
48 {
49 // /////////////////////////////////////////
50 // killing primary particles (any heavy charged particle) that hit the
51 // collimator
52 auto part_name = step->GetTrack()->GetParticleDefinition()->GetParticleName();
53 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "collPhys"
54 && step->GetTrack()->GetParticleDefinition()->GetParticleName() == "alpha")
55 {
56 step->GetTrack()->SetTrackStatus(fStopAndKill);
57 fEventAction->InvalidEvent();
58 return; // nothing more to do in such case
59 }
60
61 // /////////////////////////////////////////
62 // recording energy of the primary particle at three different locations:
63 // 1. Boundary between World and Mylar wall, i.e. initial energy of the
64 // particle.
65 // 2. Boundary between World and silicon detector, i.e. final energy of the
66 // particle.
67 // 3. Boundary between wall inner layer (innerWall) and Target, i.e.
68 // interaction energy of the particle.
69
70 if (step->GetTrack()->GetParticleDefinition()->GetParticleName() == "alpha") {
71 // 1. and 2. - particle crossing the boundary between World volume
72 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "worldPhys") {
73 auto physVol = step->GetPostStepPoint()->GetPhysicalVolume();
74 // the above may be a null pointer as the particle can leave world, so it
75 // has to be tested:
76 if (physVol) {
77 // 1. to Mylar wall
78 if (physVol->GetName() == "wallPhys") {
79 fEventAction->RecordInitialEnergy(step->GetTrack()->GetKineticEnergy());
80 }
81 // 2. to silicon detector
82 if (physVol->GetName() == "enDetPhys") {
83 fEventAction->RecordFinalEnergy(step->GetTrack()->GetKineticEnergy());
84 }
85 }
86 }
87
88 // 3. - particle crossing the boundary from wall inner layer to the target
89 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "innerWallPhys") {
90 if (step->GetPostStepPoint()->GetPhysicalVolume()->GetName() == "targetPhys") {
91 fEventAction->RecordInteractionEnergy(step->GetTrack()->GetKineticEnergy());
92 }
93 }
94 }
95
96 // ////////////////////////////////////////
97 // counting individual ionization acts
98 // only ionisations in the target will be recorded
99 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "targetPhys") {
100 if (step->GetTotalEnergyDeposit() != 0) {
101 G4String process_name = step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
102
103 // counting ionisations caused only by electrons and alpha particles:
104 // if(process_name=="e-_G4DNAIonisation" ||
105 // process_name=="alpha_G4DNAIonisation");
106 // fEventAction->AddIonisation();
107
108 // more general approach - counting ionisations caused by any particle
109 // (including also protons, GenericIons, etc.):
110 auto pos =
111 process_name.find("Ionisation"); // looking for "Ionisation" substring in process_name
112 if (pos != std::string::npos) // found!
113 fEventAction->AddIonisation(); // increase cluster size
114 }
115 }
116 }
117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
118