| 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 // -------------------------------------------
27 // MONTE CARLO SIMULATION OF REALISTIC G
28 //
29 // Authors and contributors:
30 // P. Barberet, S. Incerti, N. H. Tran, L. Mor
31 //
32 // University of Bordeaux, CNRS, LP2i, UMR5797
33 //
34 // If you use this code, please cite the follo
35 // P. Barberet et al.,
36 // "Monte-Carlo dosimetry on a realistic cell
37 // geometry exposed to alpha particles."
38 // Ph. Barberet et al 2012 Phys. Med. Biol. 57
39 // doi: 110.1088/0031-9155/57/8/2189
40 // -------------------------------------------
41
42 #include "SteppingAction.hh"
43
44 #include "G4SteppingManager.hh"
45
46 //....oooOO0OOooo........oooOO0OOooo........oo
47
48 SteppingAction::SteppingAction(RunAction* runA
49 :G4UserSteppingAction(), fRunAction(runAction)
50 {}
51
52
53 //....oooOO0OOooo........oooOO0OOooo........oo
54
55 void SteppingAction::UserSteppingAction(const
56 {
57 // *****************************************
58 // Avoid string comparison to extract materi
59 // *****************************************
60
61 fMyPhantomParam = CellParameterisation::Inst
62 const G4StepPoint* preStep = aStep->GetPreSt
63 G4int preReplicaNumber = preStep->GetTouchab
64 G4int voxelMaterial = fMyPhantomParam->GetM
65
66 // The absorbed energy is added to the "voxe
67 // Added protection to make sure Replica Num
68
69 if (aStep->GetTotalEnergyDeposit()>0. && pre
70 {
71 if (voxelMaterial == 1)
72 {
73 fRunAction->AddDoseBox(preReplicaNumber,
74 }
75 else if (voxelMaterial == 2)
76 {
77 fRunAction->AddDoseBox(preReplicaNumber,
78 }
79 else if (voxelMaterial == 3)
80 {
81 fRunAction->AddDoseBox(preReplicaNumber,
82 }
83 }
84 }
85