| 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 // The Geant4-DNA web site is available at htt 30 // The Geant4-DNA web site is available at http://geant4-dna.org
31 // 31 //
32 // If you use this example, please cite the fo 32 // If you use this example, please cite the following publication:
33 // Rad. Prot. Dos. 133 (2009) 2-11 33 // Rad. Prot. Dos. 133 (2009) 2-11
34 // << 34
>> 35 #include "G4SystemOfUnits.hh"
>> 36 #include "G4SteppingManager.hh"
>> 37
35 #include "SteppingAction.hh" 38 #include "SteppingAction.hh"
36 #include "RunAction.hh" 39 #include "RunAction.hh"
37 #include "DetectorConstruction.hh" 40 #include "DetectorConstruction.hh"
38 #include "G4AnalysisManager.hh" << 41 #include "Analysis.hh"
39 #include "G4SystemOfUnits.hh" << 42
40 #include "G4SteppingManager.hh" <<
41 #include "G4Alpha.hh" 43 #include "G4Alpha.hh"
42 #include "G4Electron.hh" 44 #include "G4Electron.hh"
43 45
>> 46 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 47
44 SteppingAction::SteppingAction(RunAction* run, 48 SteppingAction::SteppingAction(RunAction* run,const DetectorConstruction* det)
45 :fRun(run),fDetector(det) 49 :fRun(run),fDetector(det)
46 {} 50 {}
47 51
>> 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
>> 53
48 SteppingAction::~SteppingAction() 54 SteppingAction::~SteppingAction()
49 {} 55 {}
>> 56
>> 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50 58
51 void SteppingAction::UserSteppingAction(const 59 void SteppingAction::UserSteppingAction(const G4Step* aStep)
52 60
53 { 61 {
54 // Analysis manager 62 // Analysis manager
55 63
56 G4AnalysisManager* man = G4AnalysisManager:: 64 G4AnalysisManager* man = G4AnalysisManager::Instance();
57 65
58 // Read phantom - Singleton 66 // Read phantom - Singleton
59 67
60 fMyCellParameterisation = CellParameterisati 68 fMyCellParameterisation = CellParameterisation::Instance();
61 69
62 // Material : 1 is cytoplasm, 2 is nucleus 70 // Material : 1 is cytoplasm, 2 is nucleus
63 71
64 G4int matVoxelPRE = -1; 72 G4int matVoxelPRE = -1;
65 G4int matVoxelPOST = -1; 73 G4int matVoxelPOST = -1;
66 74
67 const G4StepPoint* preStep = aStep->GetPreSt 75 const G4StepPoint* preStep = aStep->GetPreStepPoint();
68 const G4StepPoint* postStep = aStep->GetPost 76 const G4StepPoint* postStep = aStep->GetPostStepPoint();
69 const G4Track* track = aStep->GetTrack(); 77 const G4Track* track = aStep->GetTrack();
70 78
71 const G4LogicalVolume* preVolume = 79 const G4LogicalVolume* preVolume =
72 preStep->GetPhysicalVolume()->GetLogicalVo 80 preStep->GetPhysicalVolume()->GetLogicalVolume();
73 81
74 const G4LogicalVolume* postVolume = nullptr; 82 const G4LogicalVolume* postVolume = nullptr;
75 if(postStep->GetPhysicalVolume()) 83 if(postStep->GetPhysicalVolume())
76 { 84 {
77 postVolume = postStep->GetPhysicalVolume() 85 postVolume = postStep->GetPhysicalVolume()->GetLogicalVolume();
78 } 86 }
79 const G4ParticleDefinition* particle = 87 const G4ParticleDefinition* particle =
80 track->GetDynamicParticle()->GetDefinition 88 track->GetDynamicParticle()->GetDefinition();
81 89
82 G4int preReplicaNumber = preStep->GetTouchab 90 G4int preReplicaNumber = preStep->GetTouchableHandle()->GetReplicaNumber();
83 G4double edep = aStep->GetTotalEnergyDeposit 91 G4double edep = aStep->GetTotalEnergyDeposit();
84 92
85 if (preReplicaNumber>0) 93 if (preReplicaNumber>0)
86 { 94 {
87 matVoxelPRE = fMyCellParameterisation->Ge 95 matVoxelPRE = fMyCellParameterisation->GetTissueType(preReplicaNumber);
88 } 96 }
89 97
90 if(postVolume) 98 if(postVolume)
91 { 99 {
92 G4int postReplicaNumber = postStep->GetTou 100 G4int postReplicaNumber = postStep->GetTouchableHandle()->GetReplicaNumber();
93 if (postReplicaNumber>0) 101 if (postReplicaNumber>0)
94 { 102 {
95 matVoxelPOST = fMyCellParameterisation-> 103 matVoxelPOST = fMyCellParameterisation->GetTissueType(postReplicaNumber);
96 } 104 }
97 } 105 }
98 106
99 // COUNT GAS DETECTOR HITS 107 // COUNT GAS DETECTOR HITS
100 108
101 if (particle == G4Alpha::AlphaDefinition()) 109 if (particle == G4Alpha::AlphaDefinition())
102 { 110 {
103 if(postVolume == fDetector->GetLogicalIsob 111 if(postVolume == fDetector->GetLogicalIsobutane() &&
104 ((preVolume == fDetector->GetLogicalCol 112 ((preVolume == fDetector->GetLogicalCollDetYoke())
105 || 113 ||
106 (preVolume == fDetector->GetLogicalCollDetGa 114 (preVolume == fDetector->GetLogicalCollDetGap4())
107 || 115 ||
108 (preVolume == fDetector->GetLogicalCollDetGa 116 (preVolume == fDetector->GetLogicalCollDetGap4())))
109 { 117 {
110 fRun->AddNbOfHitsGas(); 118 fRun->AddNbOfHitsGas();
111 } 119 }
112 120
113 // STOPPING POWER AND BEAM SPOT SIZE AT CE 121 // STOPPING POWER AND BEAM SPOT SIZE AT CELL ENTRANCE
114 if(preVolume == fDetector->GetLogicalPolyp 122 if(preVolume == fDetector->GetLogicalPolyprop() &&
115 ( (postVolume == fDetector->GetLogicalK 123 ( (postVolume == fDetector->GetLogicalKgm()) ||
116 (matVoxelPOST == 1)) ) 124 (matVoxelPOST == 1)) )
117 { 125 {
118 G4double deltaE = preStep->GetKineticEne 126 G4double deltaE = preStep->GetKineticEnergy()
119 - postStep->GetKineticEnergy(); 127 - postStep->GetKineticEnergy();
120 if(deltaE > 0.0) 128 if(deltaE > 0.0)
121 { 129 {
122 //Fill ntupleid=1 130 //Fill ntupleid=1
123 man->FillNtupleDColumn(1,0,preStep->GetKinet 131 man->FillNtupleDColumn(1,0,preStep->GetKineticEnergy()/keV);
124 man->FillNtupleDColumn(1,1,deltaE*micrometer 132 man->FillNtupleDColumn(1,1,deltaE*micrometer/(keV*aStep->GetStepLength()));
125 man->AddNtupleRow(1); 133 man->AddNtupleRow(1);
126 } 134 }
127 135
128 // Average dE over step suggested by Mic 136 // Average dE over step suggested by Michel Maire
129 G4ThreeVector coord1 = preStep->GetPosit 137 G4ThreeVector coord1 = preStep->GetPosition();
130 const G4AffineTransform transformation1 138 const G4AffineTransform transformation1 =
131 preStep->GetTouchable()->GetHistory()->GetTo 139 preStep->GetTouchable()->GetHistory()->GetTopTransform();
132 G4ThreeVector localPosition1 = transform 140 G4ThreeVector localPosition1 = transformation1.TransformPoint(coord1);
133 141
134 G4ThreeVector coord2 = postStep->GetPosi 142 G4ThreeVector coord2 = postStep->GetPosition();
135 const G4AffineTransform transformation2 143 const G4AffineTransform transformation2 =
136 postStep->GetTouchable()->GetHistory()->GetT 144 postStep->GetTouchable()->GetHistory()->GetTopTransform();
137 G4ThreeVector localPosition2 = transform 145 G4ThreeVector localPosition2 = transformation2.TransformPoint(coord2);
138 146
139 G4ThreeVector localPosition = 147 G4ThreeVector localPosition =
140 localPosition1 + G4UniformRand()*(localPosit 148 localPosition1 + G4UniformRand()*(localPosition2-localPosition1);
141 149
142 //Fill ntupleid=2 150 //Fill ntupleid=2
143 man->FillNtupleDColumn(2,0,localPosition 151 man->FillNtupleDColumn(2,0,localPosition.x()/micrometer);
144 man->FillNtupleDColumn(2,1,localPosition 152 man->FillNtupleDColumn(2,1,localPosition.y()/micrometer);
145 man->AddNtupleRow(2); 153 man->AddNtupleRow(2);
146 } 154 }
147 155
148 // ALPHA RANGE 156 // ALPHA RANGE
149 if (postStep->GetKineticEnergy() < eV && 157 if (postStep->GetKineticEnergy() < eV &&
150 ( (matVoxelPOST==1) || 158 ( (matVoxelPOST==1) ||
151 (postVolume == fDetector->GetLogicalKgm()) 159 (postVolume == fDetector->GetLogicalKgm()) ||
152 (matVoxelPOST==2) ) ) 160 (matVoxelPOST==2) ) )
153 { 161 {
154 //Fill ntupleid=3 162 //Fill ntupleid=3
155 man->FillNtupleDColumn(3,0,postStep->Get 163 man->FillNtupleDColumn(3,0,postStep->GetPosition().x()/micrometer);
156 man->FillNtupleDColumn(3,1,postStep->Get 164 man->FillNtupleDColumn(3,1,postStep->GetPosition().y()/micrometer);
157 man->FillNtupleDColumn(3,2,postStep->Get 165 man->FillNtupleDColumn(3,2,postStep->GetPosition().z()/micrometer);
158 man->AddNtupleRow(3); 166 man->AddNtupleRow(3);
159 } 167 }
160 168
161 // TOTAL DOSE DEPOSIT AND DOSE DEPOSIT WIT 169 // TOTAL DOSE DEPOSIT AND DOSE DEPOSIT WITHIN A PHANTOM VOXEL
162 // FOR ALL PARTICLES 170 // FOR ALL PARTICLES
163 } 171 }
164 172
165 if (matVoxelPRE == 2) 173 if (matVoxelPRE == 2)
166 { 174 {
167 G4double dose = (edep/joule)/(fRun->GetMas 175 G4double dose = (edep/joule)/(fRun->GetMassNucleus()/kg);
168 fRun->AddDoseN(dose); 176 fRun->AddDoseN(dose);
169 fRun->AddDoseBox(preReplicaNumber, edep/eV 177 fRun->AddDoseBox(preReplicaNumber, edep/eV);
170 } 178 }
171 else if (matVoxelPRE == 1) 179 else if (matVoxelPRE == 1)
172 { 180 {
173 G4double dose = (edep/joule)/(fRun->GetMas 181 G4double dose = (edep/joule)/(fRun->GetMassCytoplasm()/kg);
174 fRun->AddDoseC(dose); 182 fRun->AddDoseC(dose);
175 fRun->AddDoseBox(preReplicaNumber, edep/eV 183 fRun->AddDoseBox(preReplicaNumber, edep/eV);
176 } 184 }
177 185
178 // PROTECTION AGAINST POSSIBLE MSC LOOPS FOR 186 // PROTECTION AGAINST POSSIBLE MSC LOOPS FOR e-
179 187
180 // if ( edep/MeV<1e-25 && particle == G4Ele 188 // if ( edep/MeV<1e-25 && particle == G4Electron::Electron())
181 // { 189 // {
182 //aStep->GetTrack()->SetTrackStatus(fSto 190 //aStep->GetTrack()->SetTrackStatus(fStopAndKill);
183 /* 191 /*
184 G4cout << "*** Warning *** : msc loop for " 192 G4cout << "*** Warning *** : msc loop for "
185 << track->GetDefinition()->GetParticleName() 193 << track->GetDefinition()->GetParticleName()
186 << " in " << 194 << " in " <<
187 postPoint->GetTouchableHandle()->GetVolume() 195 postPoint->GetTouchableHandle()->GetVolume()->GetName() << G4endl;
188 */ 196 */
189 // } 197 // }
190 } 198 }
191 199