| 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 // ------------------------------------------- 26 // -------------------------------------------------------------------
27 // ------------------------------------------- 27 // -------------------------------------------------------------------
28 28
29 #include "RunAction.hh" 29 #include "RunAction.hh"
30 #include "G4Run.hh" 30 #include "G4Run.hh"
31 #include "TrackingAction.hh" 31 #include "TrackingAction.hh"
32 #include "G4ParticleDefinition.hh" 32 #include "G4ParticleDefinition.hh"
33 #include "G4RunManager.hh" 33 #include "G4RunManager.hh"
34 #include "G4AnalysisManager.hh" 34 #include "G4AnalysisManager.hh"
35 #include "G4Threading.hh" 35 #include "G4Threading.hh"
36 36
37 void PrintNParticles(std::map<const G4Particle 37 void PrintNParticles(std::map<const G4ParticleDefinition*, int>& container);
38 38
39 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
40 40
41 RunAction::RunAction() 41 RunAction::RunAction()
42 { 42 {
43 fFileName = "microelectronics"; 43 fFileName = "microelectronics";
44 fpTrackingAction = 0; 44 fpTrackingAction = 0;
45 fInitialized = 0; 45 fInitialized = 0;
46 fDebug = false; 46 fDebug = false;
47 } 47 }
48 48
49 //....oooOO0OOooo........oooOO0OOooo........oo 49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
50 50
51 RunAction::~RunAction() 51 RunAction::~RunAction()
52 {} 52 {}
53 53
54 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
55 55
56 void RunAction::BeginOfRunAction(const G4Run* 56 void RunAction::BeginOfRunAction(const G4Run* run)
57 { 57 {
58 // In this example, we considered that the s 58 // In this example, we considered that the same class was
59 // used for both master and worker threads. 59 // used for both master and worker threads.
60 // However, in case the run action is long, 60 // However, in case the run action is long,
61 // for better code review, this practice is 61 // for better code review, this practice is not recommanded.
62 // 62 //
63 63
64 if(isMaster) // WARNING : in sequential mode 64 if(isMaster) // WARNING : in sequential mode, isMaster == true
65 BeginMaster(run); 65 BeginMaster(run);
66 else 66 else
67 BeginWorker(run); 67 BeginWorker(run);
68 } 68 }
69 69
70 //....oooOO0OOooo........oooOO0OOooo........oo 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
71 71
72 void RunAction::EndOfRunAction(const G4Run* ru 72 void RunAction::EndOfRunAction(const G4Run* run)
73 { 73 {
74 if(isMaster) 74 if(isMaster)
75 EndMaster(run); 75 EndMaster(run);
76 else 76 else
77 EndWorker(run); 77 EndWorker(run);
78 } 78 }
79 79
80 80
81 //....oooOO0OOooo........oooOO0OOooo........oo 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
82 82
83 void RunAction::BeginMaster(const G4Run* run) 83 void RunAction::BeginMaster(const G4Run* run)
84 { 84 {
85 bool sequential = (G4RunManager::GetRunManag 85 bool sequential = (G4RunManager::GetRunManager()->GetRunManagerType() == G4RunManager::sequentialRM);
86 86
87 if(fDebug) 87 if(fDebug)
88 { 88 {
89 G4cout << "°°°°°°°°°°°°°° 89 G4cout << "°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°" << G4endl;
90 if(!sequential) 90 if(!sequential)
91 G4cout << "°°°°°°°°°°°°°°°° 91 G4cout << "°°°°°°°°°°°°°°°° RunAction::BeginMaster" << G4endl;
92 PrintRunInfo(run); 92 PrintRunInfo(run);
93 G4cout << "°°°°°°°°°°°°°° 93 G4cout << "°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°" << G4endl;
94 } 94 }
95 95
96 if(sequential) 96 if(sequential)
97 { 97 {
98 if(!fInitialized) 98 if(!fInitialized)
99 InitializeWorker(run); 99 InitializeWorker(run);
100 // Note: fpTrackingAction could be used 100 // Note: fpTrackingAction could be used as a flag for initialization instead
101 101
102 CreateHistogram(); 102 CreateHistogram();
103 } 103 }
104 } 104 }
105 105
106 //....oooOO0OOooo........oooOO0OOooo........oo 106 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
107 107
108 void RunAction::BeginWorker(const G4Run* run) 108 void RunAction::BeginWorker(const G4Run* run)
109 { 109 {
110 if(fDebug) 110 if(fDebug)
111 { 111 {
112 G4cout << "°°°°°°°°°°°°°° 112 G4cout << "°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°" << G4endl;
113 G4cout << "°°°°°°°°°°°°°° 113 G4cout << "°°°°°°°°°°°°°°°° RunAction::BeginWorker" << G4endl;
114 PrintRunInfo(run); 114 PrintRunInfo(run);
115 G4cout << "°°°°°°°°°°°°°° 115 G4cout << "°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°°" << G4endl;
116 } 116 }
117 if(!fInitialized) 117 if(!fInitialized)
118 InitializeWorker(run); 118 InitializeWorker(run);
119 119
120 CreateHistogram(); 120 CreateHistogram();
121 } 121 }
122 122
123 //....oooOO0OOooo........oooOO0OOooo........oo 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
124 124
125 void RunAction::EndMaster(const G4Run* run) 125 void RunAction::EndMaster(const G4Run* run)
126 { 126 {
127 bool sequential = (G4RunManager::GetRunManag 127 bool sequential = (G4RunManager::GetRunManager()->GetRunManagerType()
128 == G4RunManager::sequentialRM); 128 == G4RunManager::sequentialRM);
129 if(sequential) 129 if(sequential)
130 EndWorker(run); 130 EndWorker(run);
131 } 131 }
132 132
133 //....oooOO0OOooo........oooOO0OOooo........oo 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
134 134
135 void RunAction::EndWorker(const G4Run* run) 135 void RunAction::EndWorker(const G4Run* run)
136 { 136 {
137 if(fDebug) 137 if(fDebug)
138 { 138 {
139 PrintRunInfo(run); 139 PrintRunInfo(run);
140 } 140 }
141 141
142 G4int nofEvents = run->GetNumberOfEvent(); 142 G4int nofEvents = run->GetNumberOfEvent();
143 if ( nofEvents == 0 ) 143 if ( nofEvents == 0 )
144 { 144 {
145 if(fDebug) 145 if(fDebug)
146 { 146 {
147 G4cout << "°°°°°°°°°°°°°°° 147 G4cout << "°°°°°°°°°°°°°°°° NO EVENTS TREATED IN THIS RUN ==> LEAVING RunAction::EndOfRunAction "<< G4endl;
148 } 148 }
149 return; 149 return;
150 } 150 }
151 151
152 /////////////// 152 ///////////////
153 // Write Histo 153 // Write Histo
154 // 154 //
155 WriteHistogram(); 155 WriteHistogram();
156 156
157 /////////////// 157 ///////////////
158 // Complete cleanup 158 // Complete cleanup
159 // 159 //
160 G4AnalysisManager* analysisManager = G4Analy 160 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
161 analysisManager->Clear(); 161 analysisManager->Clear();
162 162
163 /////////////// 163 ///////////////
164 // Printouts 164 // Printouts
165 // 165 //
166 std::map<const G4ParticleDefinition*, int>& 166 std::map<const G4ParticleDefinition*, int>&
167 particlesCreatedInWorld = fpTrackingAction 167 particlesCreatedInWorld = fpTrackingAction->GetNParticlesCreatedInWorld();
168 168
169 G4cout << "Number and type of particles crea 169 G4cout << "Number and type of particles created outside region \"Target\" :" << G4endl;
170 170
171 PrintNParticles(particlesCreatedInWorld); 171 PrintNParticles(particlesCreatedInWorld);
172 172
173 G4cout << "_______________________" << G4end 173 G4cout << "_______________________" << G4endl;
174 174
175 std::map<const G4ParticleDefinition*, int>& 175 std::map<const G4ParticleDefinition*, int>&
176 particlesCreatedInTarget = fpTrackingActio 176 particlesCreatedInTarget = fpTrackingAction->GetNParticlesCreatedInTarget();
177 177
178 G4cout << "Number and type of particles crea 178 G4cout << "Number and type of particles created in region \"Target\" :" << G4endl;
179 179
180 PrintNParticles(particlesCreatedInTarget); 180 PrintNParticles(particlesCreatedInTarget);
181 181
182 } 182 }
183 183
184 //....oooOO0OOooo........oooOO0OOooo........oo 184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
185 185
186 void RunAction::InitializeWorker(const G4Run*) 186 void RunAction::InitializeWorker(const G4Run*)
187 { 187 {
188 if (fpTrackingAction == 0) 188 if (fpTrackingAction == 0)
189 { 189 {
190 fpTrackingAction = (TrackingAction*) G4Run 190 fpTrackingAction = (TrackingAction*) G4RunManager::GetRunManager()->GetUserTrackingAction();
191 191
192 if(fpTrackingAction == 0 && isMaster == fa 192 if(fpTrackingAction == 0 && isMaster == false)
193 { 193 {
194 G4ExceptionDescription exDescrption ; 194 G4ExceptionDescription exDescrption ;
195 exDescrption << "fpTrackingAction is a n 195 exDescrption << "fpTrackingAction is a null pointer. Has it been correctly initialized ?";
196 G4Exception("RunAction::BeginOfRunAction 196 G4Exception("RunAction::BeginOfRunAction","RunAction001",FatalException, exDescrption);
197 } 197 }
198 } 198 }
199 199
200 fInitialized = true; 200 fInitialized = true;
201 } 201 }
202 202
203 //....oooOO0OOooo........oooOO0OOooo........oo 203 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
204 204
205 void RunAction::CreateHistogram() 205 void RunAction::CreateHistogram()
206 { 206 {
207 // Book histograms, ntuple 207 // Book histograms, ntuple
208 208
209 // Create analysis manager 209 // Create analysis manager
210 // The choice of analysis technology is done 210 // The choice of analysis technology is done via selection of a namespace
211 // in Analysis.hh 211 // in Analysis.hh
212 212
213 G4cout << "##### Create analysis manager " < 213 G4cout << "##### Create analysis manager " << " " << this << G4endl;
214 G4AnalysisManager* analysisManager = G4Analy 214 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
215 analysisManager->SetDefaultFileType("root"); 215 analysisManager->SetDefaultFileType("root");
216 216
217 G4cout << "Using " << analysisManager->GetTy 217 G4cout << "Using " << analysisManager->GetType() << " analysis manager" << G4endl;
218 218
219 // Create directories 219 // Create directories
220 220
221 //analysisManager->SetHistoDirectoryName("hi 221 //analysisManager->SetHistoDirectoryName("histograms");
222 //analysisManager->SetNtupleDirectoryName("n 222 //analysisManager->SetNtupleDirectoryName("ntuple");
223 analysisManager->SetVerboseLevel(1); 223 analysisManager->SetVerboseLevel(1);
224 224
225 // Open an output file 225 // Open an output file
226 226
227 analysisManager->OpenFile(fFileName); 227 analysisManager->OpenFile(fFileName);
228 228
229 // Creating ntuple 229 // Creating ntuple
230 230
231 analysisManager->CreateNtuple("microelectron 231 analysisManager->CreateNtuple("microelectronics", "physics");
232 analysisManager->CreateNtupleDColumn("flagPa 232 analysisManager->CreateNtupleDColumn("flagParticle");
233 analysisManager->CreateNtupleDColumn("flagPr 233 analysisManager->CreateNtupleDColumn("flagProcess");
234 analysisManager->CreateNtupleDColumn("x"); 234 analysisManager->CreateNtupleDColumn("x");
235 analysisManager->CreateNtupleDColumn("y"); 235 analysisManager->CreateNtupleDColumn("y");
236 analysisManager->CreateNtupleDColumn("z"); 236 analysisManager->CreateNtupleDColumn("z");
237 analysisManager->CreateNtupleDColumn("totalE 237 analysisManager->CreateNtupleDColumn("totalEnergyDeposit");
238 analysisManager->CreateNtupleDColumn("stepLe 238 analysisManager->CreateNtupleDColumn("stepLength");
239 analysisManager->CreateNtupleDColumn("kineti 239 analysisManager->CreateNtupleDColumn("kineticEnergyDifference");
240 analysisManager->FinishNtuple(); 240 analysisManager->FinishNtuple();
241 } 241 }
242 242
243 //....oooOO0OOooo........oooOO0OOooo........oo 243 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
244 244
245 void RunAction::WriteHistogram() 245 void RunAction::WriteHistogram()
246 { 246 {
247 // print histogram statistics 247 // print histogram statistics
248 // 248 //
249 G4AnalysisManager* analysisManager = G4Analy 249 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
250 250
251 // save histograms 251 // save histograms
252 // 252 //
253 analysisManager->Write(); 253 analysisManager->Write();
254 analysisManager->CloseFile(); 254 analysisManager->CloseFile();
255 } 255 }
256 256
257 //....oooOO0OOooo........oooOO0OOooo........oo 257 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
258 258
259 void RunAction::PrintRunInfo(const G4Run* run) 259 void RunAction::PrintRunInfo(const G4Run* run)
260 { 260 {
261 G4cout << "°°°°°°°°°°°°°°°° 261 G4cout << "°°°°°°°°°°°°°°°° Run is = " << run->GetRunID() << G4endl;
262 G4cout << "°°°°°°°°°°°°°°°° 262 G4cout << "°°°°°°°°°°°°°°°° Run type is = " << G4RunManager::GetRunManager()->GetRunManagerType() << G4endl;
263 G4cout << "°°°°°°°°°°°°°°°° 263 G4cout << "°°°°°°°°°°°°°°°° Event processed = " << run->GetNumberOfEventToBeProcessed() << G4endl;
264 G4cout << "°°°°°°°°°°°°°°°° 264 G4cout << "°°°°°°°°°°°°°°°° N° Event = " << run->GetNumberOfEvent() << G4endl;
265 } 265 }
266 266
267 //....oooOO0OOooo........oooOO0OOooo........oo 267 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
268 void PrintNParticles(std::map<const G4Particle 268 void PrintNParticles(std::map<const G4ParticleDefinition*, int>& container)
269 { 269 {
270 std::map<const G4ParticleDefinition*, int> 270 std::map<const G4ParticleDefinition*, int>::iterator it;
271 for(it = container.begin() ; 271 for(it = container.begin() ;
272 it != container.end(); it ++) 272 it != container.end(); it ++)
273 { 273 {
274 G4cout << "N " << it->first->GetPartic 274 G4cout << "N " << it->first->GetParticleName() << " : " << it->second << G4endl;
275 } 275 }
276 } 276 }
277 277
278 //....oooOO0OOooo........oooOO0OOooo........oo 278 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
279 279