| 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 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 // J. Comput. Phys. 274 (2014) 841-882 30 // J. Comput. Phys. 274 (2014) 841-882
31 // The Geant4-DNA web site is available at htt 31 // The Geant4-DNA web site is available at http://geant4-dna.org
32 // 32 //
33 // 33 //
34 /// \file SteppingAction.cc 34 /// \file SteppingAction.cc
35 /// \brief Implementation of the SteppingActio 35 /// \brief Implementation of the SteppingAction class
36 36
37 #include "SteppingAction.hh" 37 #include "SteppingAction.hh"
38 <<
39 #include "EventAction.hh" 38 #include "EventAction.hh"
40 <<
41 #include "G4SteppingManager.hh" 39 #include "G4SteppingManager.hh"
42 40
43 //....oooOO0OOooo........oooOO0OOooo........oo 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 SteppingAction::SteppingAction(EventAction* ev << 42 SteppingAction::SteppingAction(EventAction *event)
>> 43 : G4UserSteppingAction()
>> 44 , fEventAction(event) {}
45 45
46 //....oooOO0OOooo........oooOO0OOooo........oo 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
47 void SteppingAction::UserSteppingAction(const << 47 void SteppingAction::UserSteppingAction(const G4Step *step) {
48 { <<
49 // ///////////////////////////////////////// 48 // /////////////////////////////////////////
50 // killing primary particles (any heavy char 49 // killing primary particles (any heavy charged particle) that hit the
51 // collimator 50 // collimator
52 auto part_name = step->GetTrack()->GetPartic 51 auto part_name = step->GetTrack()->GetParticleDefinition()->GetParticleName();
53 if (step->GetPreStepPoint()->GetPhysicalVolu << 52 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "collPhys" &&
54 && step->GetTrack()->GetParticleDefiniti << 53 step->GetTrack()->GetParticleDefinition()->GetParticleName() == "alpha") {
55 { <<
56 step->GetTrack()->SetTrackStatus(fStopAndK 54 step->GetTrack()->SetTrackStatus(fStopAndKill);
57 fEventAction->InvalidEvent(); 55 fEventAction->InvalidEvent();
58 return; // nothing more to do in such cas << 56 return; // nothing more to do in such case
59 } 57 }
60 58
61 // ///////////////////////////////////////// 59 // /////////////////////////////////////////
62 // recording energy of the primary particle 60 // recording energy of the primary particle at three different locations:
63 // 1. Boundary between World and Mylar wall, 61 // 1. Boundary between World and Mylar wall, i.e. initial energy of the
64 // particle. 62 // particle.
65 // 2. Boundary between World and silicon det 63 // 2. Boundary between World and silicon detector, i.e. final energy of the
66 // particle. 64 // particle.
67 // 3. Boundary between wall inner layer (inn 65 // 3. Boundary between wall inner layer (innerWall) and Target, i.e.
68 // interaction energy of the particle. 66 // interaction energy of the particle.
69 67
70 if (step->GetTrack()->GetParticleDefinition( 68 if (step->GetTrack()->GetParticleDefinition()->GetParticleName() == "alpha") {
71 // 1. and 2. - particle crossing the bound 69 // 1. and 2. - particle crossing the boundary between World volume
72 if (step->GetPreStepPoint()->GetPhysicalVo << 70 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() ==
>> 71 "worldPhys") {
73 auto physVol = step->GetPostStepPoint()- 72 auto physVol = step->GetPostStepPoint()->GetPhysicalVolume();
74 // the above may be a null pointer as th 73 // the above may be a null pointer as the particle can leave world, so it
75 // has to be tested: 74 // has to be tested:
76 if (physVol) { 75 if (physVol) {
77 // 1. to Mylar wall 76 // 1. to Mylar wall
78 if (physVol->GetName() == "wallPhys") 77 if (physVol->GetName() == "wallPhys") {
79 fEventAction->RecordInitialEnergy(st << 78 fEventAction->RecordInitialEnergy(
>> 79 step->GetTrack()->GetKineticEnergy());
80 } 80 }
81 // 2. to silicon detector 81 // 2. to silicon detector
82 if (physVol->GetName() == "enDetPhys") 82 if (physVol->GetName() == "enDetPhys") {
83 fEventAction->RecordFinalEnergy(step 83 fEventAction->RecordFinalEnergy(step->GetTrack()->GetKineticEnergy());
84 } 84 }
85 } 85 }
86 } 86 }
87 87
88 // 3. - particle crossing the boundary fro 88 // 3. - particle crossing the boundary from wall inner layer to the target
89 if (step->GetPreStepPoint()->GetPhysicalVo << 89 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() ==
90 if (step->GetPostStepPoint()->GetPhysica << 90 "innerWallPhys") {
91 fEventAction->RecordInteractionEnergy( << 91 if (step->GetPostStepPoint()->GetPhysicalVolume()->GetName() ==
>> 92 "targetPhys") {
>> 93 fEventAction->RecordInteractionEnergy(
>> 94 step->GetTrack()->GetKineticEnergy());
92 } 95 }
93 } 96 }
94 } 97 }
95 98
96 // //////////////////////////////////////// 99 // ////////////////////////////////////////
97 // counting individual ionization acts 100 // counting individual ionization acts
98 // only ionisations in the target will be re 101 // only ionisations in the target will be recorded
99 if (step->GetPreStepPoint()->GetPhysicalVolu 102 if (step->GetPreStepPoint()->GetPhysicalVolume()->GetName() == "targetPhys") {
100 if (step->GetTotalEnergyDeposit() != 0) { 103 if (step->GetTotalEnergyDeposit() != 0) {
101 G4String process_name = step->GetPostSte << 104 G4String process_name =
>> 105 step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
102 106
103 // counting ionisations caused only by e 107 // counting ionisations caused only by electrons and alpha particles:
104 // if(process_name=="e-_G4DNA 108 // if(process_name=="e-_G4DNAIonisation" ||
105 // process_name=="alpha_G4DNA 109 // process_name=="alpha_G4DNAIonisation");
106 // fEventAction->AddIonis 110 // fEventAction->AddIonisation();
107 111
108 // more general approach - counting ioni 112 // more general approach - counting ionisations caused by any particle
109 // (including also protons, GenericIons, 113 // (including also protons, GenericIons, etc.):
110 auto pos = << 114 auto pos = process_name.find(
111 process_name.find("Ionisation"); // l << 115 "Ionisation"); // looking for "Ionisation" substring in process_name
112 if (pos != std::string::npos) // found! << 116 if (pos != std::string::npos) // found!
113 fEventAction->AddIonisation(); // inc << 117 fEventAction->AddIonisation(); // increase cluster size
114 } 118 }
115 } 119 }
116 } 120 }
117 //....oooOO0OOooo........oooOO0OOooo........oo 121 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
118 122