Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4SteppingVerboseWithUnits.cc
Go to the documentation of this file.
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/// \file
27/// \brief Implementation of the G4SteppingVerboseWithUnits class
28//
29//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
30//
31// Stepping Verbose with units for all the applicable double values
32// This class is ported from TestEm2 extended example
33// Original author : Michel Maire (LAPP)
34// Porting with addition of UI command : Makoto Asai (SLAC) Feb.23.2021
35//
36//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
37
38#include "G4SteppingVerboseWithUnits.hh"
39
40#include "G4SteppingManager.hh"
41#include "G4UnitsTable.hh"
42#include "G4GenericMessenger.hh"
43
44//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
45
46G4SteppingVerboseWithUnits::G4SteppingVerboseWithUnits(G4int prec)
47 : G4SteppingVerbose(), fprec(prec)
48{
49}
50
51//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52
53G4SteppingVerboseWithUnits::~G4SteppingVerboseWithUnits()
54{
55 delete fmessenger;
56}
57
58//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59
60void G4SteppingVerboseWithUnits::SetManager(G4SteppingManager* const fMan)
61{
62 fManager = fMan;
63 fmessenger = new G4GenericMessenger(this,"/tracking/",
64 "precision of verbose output");
65 auto& cmd = fmessenger->DeclareProperty("setVerbosePrecision",fprec,
66 "set precision of verbose output");
67 cmd.SetStates(G4State_PreInit,G4State_Idle);
68}
69
70//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
71
72void G4SteppingVerboseWithUnits::TrackingStarted()
73{
74 CopyState();
75 G4long oldprec = G4cout.precision(fprec);
76
77 // Step zero
78 //
79 if( verboseLevel > 0 )
80 {
81 G4cout << std::setw( 5) << "Step#" << " "
82 << std::setw(fprec+3) << "X" << " "
83 << std::setw(fprec+3) << "Y" << " "
84 << std::setw(fprec+3) << "Z" << " "
85 << std::setw(fprec+6) << "KineE" << " "
86 << std::setw(fprec+10) << "dEStep" << " "
87 << std::setw(fprec+7) << "StepLeng"
88 << std::setw(fprec+7) << "TrakLeng"
89 << std::setw( 10) << "Volume" << " "
90 << std::setw( 10) << "Process" << G4endl;
91
92 G4cout << std::setw(5) << fTrack->GetCurrentStepNumber() << " "
93 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().x(),"Length")
94 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().y(),"Length")
95 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().z(),"Length")
96 << std::setw(fprec+3) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
97 << std::setw(fprec+7) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
98 << std::setw(fprec+3) << G4BestUnit(fStep->GetStepLength(),"Length")
99 << std::setw(fprec+3) << G4BestUnit(fTrack->GetTrackLength(),"Length")
100 << std::setw(10) << fTrack->GetVolume()->GetName()
101 << std::setw( 9) << " initStep" << G4endl;
102 }
103 G4cout.precision(oldprec);
104}
105
106//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
107
108void G4SteppingVerboseWithUnits::StepInfo()
109{
110 CopyState();
111 G4long oldprec = G4cout.precision(fprec);
112
113 if( verboseLevel >= 1 )
114 {
115 if( verboseLevel >= 4 ) VerboseTrack();
116 if( verboseLevel >= 3 )
117 {
118 G4cout << G4endl;
119 G4cout << std::setw(5) << "#Step#" << " "
120 << std::setw(fprec+3) << "X" << " "
121 << std::setw(fprec+3) << "Y" << " "
122 << std::setw(fprec+3) << "Z" << " "
123 << std::setw(fprec+6) << "KineE" << " "
124 << std::setw(fprec+10) << "dEStep" << " "
125 << std::setw(fprec+7) << "StepLeng"
126 << std::setw(fprec+7) << "TrakLeng"
127 << std::setw(10) << "Volume" << " "
128 << std::setw(10) << "Process" << G4endl;
129 }
130
131 G4cout << std::setw(5) << fTrack->GetCurrentStepNumber() << " "
132 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().x(),"Length")
133 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().y(),"Length")
134 << std::setw(fprec+3) << G4BestUnit(fTrack->GetPosition().z(),"Length")
135 << std::setw(fprec+3) << G4BestUnit(fTrack->GetKineticEnergy(),"Energy")
136 << std::setw(fprec+7) << G4BestUnit(fStep->GetTotalEnergyDeposit(),"Energy")
137 << std::setw(fprec+3) << G4BestUnit(fStep->GetStepLength(),"Length")
138 << std::setw(fprec+3) << G4BestUnit(fTrack->GetTrackLength(),"Length")
139 << std::setw(10) << fTrack->GetVolume()->GetName();
140
141 const G4VProcess* process
142 = fStep->GetPostStepPoint()->GetProcessDefinedStep();
143 G4String procName = " UserLimit";
144 if (process != nullptr) procName = process->GetProcessName();
145 if (fStepStatus == fWorldBoundary) procName = "OutOfWorld";
146 G4cout << " " << std::setw( 9) << procName;
147 G4cout << G4endl;
148
149 if (verboseLevel == 2)
150 {
151 const std::vector<const G4Track*>* secondary
152 = fStep->GetSecondaryInCurrentStep();
153 std::size_t nbtrk = (*secondary).size();
154 if (nbtrk)
155 {
156 G4cout << "\n :----- List of secondaries ----------------" << G4endl;
157 G4cout.precision(4);
158 for (std::size_t lp=0; lp<(*secondary).size(); ++lp)
159 {
160 G4cout << " "
161 << std::setw(13)
162 << (*secondary)[lp]->GetDefinition()->GetParticleName()
163 << ": energy ="
164 << std::setw(6)
165 << G4BestUnit((*secondary)[lp]->GetKineticEnergy(),"Energy")
166 << " time ="
167 << std::setw(6)
168 << G4BestUnit((*secondary)[lp]->GetGlobalTime(),"Time");
169 G4cout << G4endl;
170 }
171
172 G4cout << " :------------------------------------------\n" << G4endl;
173 }
174 }
175 }
176 G4cout.precision(oldprec);
177}
178
179//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
180
181void G4SteppingVerboseWithUnits::AtRestDoItInvoked()
182{
183 G4VProcess* ptProcManager;
184 CopyState();
185
186 if(verboseLevel >= 3)
187 {
188 G4int npt=0;
189 G4cout << " **List of AtRestDoIt invoked:" << G4endl;
190 for(std::size_t np=0; np<MAXofAtRestLoops; ++np)
191 {
192 std::size_t npGPIL = MAXofAtRestLoops-np-1;
193 if( (*fSelectedAtRestDoItVector)[npGPIL] == G4ForceCondition::Forced )
194 {
195 ++npt;
196 ptProcManager = (*fAtRestDoItVector)[(G4int)np];
197 G4cout << " # " << npt << " : "
198 << ptProcManager->GetProcessName()
199 << " (Forced)" << G4endl;
200 }
201 else if ( (*fSelectedAtRestDoItVector)[npGPIL] == G4ForceCondition::NotForced )
202 {
203 ++npt;
204 ptProcManager = (*fAtRestDoItVector)[(G4int)np];
205 G4cout << " # " << npt << " : " << ptProcManager->GetProcessName()
206 << G4endl;
207 }
208 }
209
210 G4cout << " Generated secondaries = " << fN2ndariesAtRestDoIt << G4endl;
211
212 if( fN2ndariesAtRestDoIt > 0 )
213 {
214 G4cout << " -- List of secondaries generated : "
215 << "(x,y,z,kE,t,PID) --" << G4endl;
216 for( std::size_t lp1=(*fSecondary).size()-fN2ndariesAtRestDoIt;
217 lp1<(*fSecondary).size(); ++lp1)
218 {
219 G4cout << " "
220 << std::setw( 9)
221 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(),"Length")
222 << " " << std::setw( 9)
223 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(),"Length")
224 << " " << std::setw( 9)
225 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(),"Length")
226 << " " << std::setw( 9)
227 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(),"Energy")
228 << " " << std::setw( 9)
229 << G4BestUnit((*fSecondary)[lp1]->GetGlobalTime(),"Time")
230 << " " << std::setw(18)
231 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName()
232 << G4endl;
233 }
234 }
235 }
236
237 if( verboseLevel >= 4 )
238 {
239 ShowStep();
240 G4cout << G4endl;
241 }
242}
243
244//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
245
246void G4SteppingVerboseWithUnits::AlongStepDoItAllDone()
247{
248 G4VProcess* ptProcManager;
249
250 CopyState();
251
252 if(verboseLevel >= 3)
253 {
254 G4cout << G4endl;
255 G4cout << " >>AlongStepDoIt (after all invocations):" << G4endl;
256 G4cout << " ++List of invoked processes " << G4endl;
257
258 for(std::size_t ci=0; ci<MAXofAlongStepLoops; ++ci)
259 {
260 ptProcManager = (*fAlongStepDoItVector)((G4int)ci);
261 G4cout << " " << ci+1 << ") ";
262 if(ptProcManager != nullptr)
263 {
264 G4cout << ptProcManager->GetProcessName() << G4endl;
265 }
266 }
267
268 ShowStep();
269 G4cout << G4endl;
270 G4cout << " ++List of secondaries generated "
271 << "(x,y,z,kE,t,PID):"
272 << " No. of secondaries = "
273 << (*fSecondary).size() << G4endl;
274
275 if((*fSecondary).size()>0)
276 {
277 for(std::size_t lp1=0; lp1<(*fSecondary).size(); ++lp1)
278 {
279 G4cout << " "
280 << std::setw( 9)
281 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(),"Length")
282 << " " << std::setw( 9)
283 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(),"Length")
284 << " " << std::setw( 9)
285 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(),"Length")
286 << " " << std::setw( 9)
287 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(),"Energy")
288 << " " << std::setw( 9)
289 << G4BestUnit((*fSecondary)[lp1]->GetGlobalTime(),"Time")
290 << " " << std::setw(18)
291 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName()
292 << G4endl;
293 }
294 }
295 }
296}
297
298//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
299
300void G4SteppingVerboseWithUnits::PostStepDoItAllDone()
301{
302 G4VProcess* ptProcManager;
303
304 CopyState();
305
306 if( (fStepStatus == fPostStepDoItProc) |
307 (fCondition == Forced) |
308 (fCondition == Conditionally) |
309 (fCondition == ExclusivelyForced) |
310 (fCondition == StronglyForced) )
311 {
312 if(verboseLevel >= 3)
313 {
314 G4int npt=0;
315 G4cout << G4endl;
316 G4cout << " **PostStepDoIt (after all invocations):" << G4endl;
317 G4cout << " ++List of invoked processes " << G4endl;
318
319 for(std::size_t np=0; np<MAXofPostStepLoops; ++np)
320 {
321 std::size_t npGPIL = MAXofPostStepLoops-np-1;
322 if( (*fSelectedPostStepDoItVector)[npGPIL] == G4ForceCondition::Forced )
323 {
324 ++npt;
325 ptProcManager = (*fPostStepDoItVector)[(G4int)np];
326 G4cout << " " << npt << ") "
327 << ptProcManager->GetProcessName()
328 << " (Forced)" << G4endl;
329 }
330 else if ( (*fSelectedPostStepDoItVector)[npGPIL] == G4ForceCondition::NotForced )
331 {
332 ++npt;
333 ptProcManager = (*fPostStepDoItVector)[(G4int)np];
334 G4cout << " " << npt << ") " << ptProcManager->GetProcessName()
335 << G4endl;
336 }
337 }
338
339 ShowStep();
340 G4cout << G4endl;
341 G4cout << " ++List of secondaries generated "
342 << "(x,y,z,kE,t,PID):"
343 << " No. of secondaries = "
344 << (*fSecondary).size() << G4endl;
345 G4cout << " [Note]Secondaries from AlongStepDoIt included."
346 << G4endl;
347
348 if((*fSecondary).size()>0)
349 {
350 for(std::size_t lp1=0; lp1<(*fSecondary).size(); ++lp1)
351 {
352 G4cout << " "
353 << std::setw( 9)
354 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(), "Length")
355 << " " << std::setw( 9)
356 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(), "Length")
357 << " " << std::setw( 9)
358 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(), "Length")
359 << " " << std::setw( 9)
360 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(), "Energy")
361 << " " << std::setw( 9)
362 << G4BestUnit((*fSecondary)[lp1]->GetGlobalTime(), "Time")
363 << " " << std::setw(18)
364 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName()
365 << G4endl;
366 }
367 }
368 }
369 }
370}
371
372//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
373
374void G4SteppingVerboseWithUnits::DPSLStarted()
375{
376 CopyState();
377
378 if( verboseLevel > 5 )
379 {
380 G4cout << G4endl
381 << " >>DefinePhysicalStepLength (List of proposed StepLengths): "
382 << G4endl;
383 }
384}
385
386//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
387
388void G4SteppingVerboseWithUnits::DPSLUserLimit()
389{
390 CopyState();
391
392 if( verboseLevel > 5 )
393 {
394 G4cout << G4endl << G4endl;
395 G4cout << "=== Defined Physical Step Length (DPSL)" << G4endl;
396 G4cout << " ++ProposedStep(UserLimit) = "
397 << std::setw( 9) << G4BestUnit(physIntLength , "Length")
398 << " : ProcName = User defined maximum allowed Step" << G4endl;
399 }
400}
401
402//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
403
404void G4SteppingVerboseWithUnits::DPSLPostStep()
405{
406 CopyState();
407
408 if( verboseLevel > 5 )
409 {
410 G4cout << " ++ProposedStep(PostStep ) = "
411 << std::setw( 9) << G4BestUnit(physIntLength , "Length")
412 << " : ProcName = " << fCurrentProcess->GetProcessName() << " (";
413 if(fCondition==ExclusivelyForced)
414 {
415 G4cout << "ExclusivelyForced)" << G4endl;
416 }
417 else if(fCondition==StronglyForced)
418 {
419 G4cout << "StronglyForced)" << G4endl;
420 }
421 else if(fCondition==Conditionally)
422 {
423 G4cout << "Conditionally)" << G4endl;
424 }
425 else if(fCondition==Forced)
426 {
427 G4cout << "Forced)" << G4endl;
428 }
429 else
430 {
431 G4cout << "No ForceCondition)" << G4endl;
432 }
433 }
434}
435
436//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
437
438void G4SteppingVerboseWithUnits::DPSLAlongStep()
439{
440 CopyState();
441
442 if( verboseLevel > 5 )
443 {
444 G4cout << " ++ProposedStep(AlongStep) = "
445 << std::setw( 9) << G4BestUnit(physIntLength , "Length")
446 << " : ProcName = "
447 << fCurrentProcess->GetProcessName()
448 << " (";
449 if(fGPILSelection==CandidateForSelection)
450 {
451 G4cout << "CandidateForSelection)" << G4endl;
452 }
453 else if(fGPILSelection==NotCandidateForSelection)
454 {
455 G4cout << "NotCandidateForSelection)" << G4endl;
456 }
457 else
458 {
459 G4cout << "?!?)" << G4endl;
460 }
461 }
462}
463
464//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
465
466void G4SteppingVerboseWithUnits::AlongStepDoItOneByOne()
467{
468 CopyState();
469
470 if(verboseLevel >= 4)
471 {
472 G4cout << G4endl;
473 G4cout << " >>AlongStepDoIt (process by process): "
474 << " Process Name = "
475 << fCurrentProcess->GetProcessName() << G4endl;
476
477 ShowStep();
478 G4cout << " "
479 << "!Note! Safety of PostStep is only valid "
480 << "after all DoIt invocations."
481 << G4endl;
482
483 VerboseParticleChange();
484 G4cout << G4endl;
485
486 G4cout << " ++List of secondaries generated "
487 << "(x,y,z,kE,t,PID):"
488 << " No. of secondaries = "
489 << fN2ndariesAlongStepDoIt << G4endl;
490
491 if(fN2ndariesAlongStepDoIt>0)
492 {
493 for(std::size_t lp1=(*fSecondary).size()-fN2ndariesAlongStepDoIt;
494 lp1<(*fSecondary).size(); ++lp1)
495 {
496 G4cout << " "
497 << std::setw( 9)
498 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(), "Length")
499 << " " << std::setw( 9)
500 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(), "Length")
501 << " " << std::setw( 9)
502 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(), "Length")
503 << " " << std::setw( 9)
504 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(), "Energy")
505 << " " << std::setw( 9)
506 << G4BestUnit((*fSecondary)[lp1]->GetGlobalTime(), "Time")
507 << " " << std::setw(18)
508 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName()
509 << G4endl;
510 }
511 }
512 }
513}
514
515//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
516
517void G4SteppingVerboseWithUnits::PostStepDoItOneByOne()
518{
519 CopyState();
520
521 if(verboseLevel >= 4)
522 {
523 G4cout << G4endl;
524 G4cout << " >>PostStepDoIt (process by process): "
525 << " Process Name = "
526 << fCurrentProcess->GetProcessName() << G4endl;
527
528 ShowStep();
529 G4cout << G4endl;
530 VerboseParticleChange();
531 G4cout << G4endl;
532
533 G4cout << " ++List of secondaries generated "
534 << "(x,y,z,kE,t,PID):"
535 << " No. of secondaries = "
536 << fN2ndariesPostStepDoIt << G4endl;
537
538 if(fN2ndariesPostStepDoIt>0)
539 {
540 for(std::size_t lp1=(*fSecondary).size()-fN2ndariesPostStepDoIt;
541 lp1<(*fSecondary).size(); ++lp1)
542 {
543 G4cout << " "
544 << std::setw( 9)
545 << G4BestUnit((*fSecondary)[lp1]->GetPosition().x(), "Length")
546 << " " << std::setw( 9)
547 << G4BestUnit((*fSecondary)[lp1]->GetPosition().y(), "Length")
548 << " " << std::setw( 9)
549 << G4BestUnit((*fSecondary)[lp1]->GetPosition().z(), "Length")
550 << " " << std::setw( 9)
551 << G4BestUnit((*fSecondary)[lp1]->GetKineticEnergy(), "Energy")
552 << " " << std::setw( 9)
553 << G4BestUnit((*fSecondary)[lp1]->GetGlobalTime(), "Time")
554 << " " << std::setw(18)
555 << (*fSecondary)[lp1]->GetDefinition()->GetParticleName()
556 << G4endl;
557 }
558 }
559 }
560}
561
562//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
563
564void G4SteppingVerboseWithUnits::VerboseTrack()
565{
566 CopyState();
567
568 G4cout << G4endl;
569 G4cout << " ++G4Track Information " << G4endl;
570 G4long oldprec = G4cout.precision(fprec);
571
572
573 G4cout << " -----------------------------------------------" << G4endl;
574 G4cout << " G4Track Information " << std::setw(20) << G4endl;
575 G4cout << " -----------------------------------------------" << G4endl;
576
577 G4cout << " Step number : "
578 << std::setw(20) << fTrack->GetCurrentStepNumber() << G4endl;
579 G4cout << " Position - x : "
580 << std::setw(20) << G4BestUnit(fTrack->GetPosition().x(), "Length")
581 << G4endl;
582 G4cout << " Position - y : "
583 << std::setw(20) << G4BestUnit(fTrack->GetPosition().y(), "Length")
584 << G4endl;
585 G4cout << " Position - z : "
586 << std::setw(20) << G4BestUnit(fTrack->GetPosition().z(), "Length")
587 << G4endl;
588 G4cout << " Global Time : "
589 << std::setw(20) << G4BestUnit(fTrack->GetGlobalTime(), "Time")
590 << G4endl;
591 G4cout << " Local Time : "
592 << std::setw(20) << G4BestUnit(fTrack->GetLocalTime(), "Time")
593 << G4endl;
594 G4cout << " Momentum Direct - x : "
595 << std::setw(20) << fTrack->GetMomentumDirection().x()
596 << G4endl;
597 G4cout << " Momentum Direct - y : "
598 << std::setw(20) << fTrack->GetMomentumDirection().y()
599 << G4endl;
600 G4cout << " Momentum Direct - z : "
601 << std::setw(20) << fTrack->GetMomentumDirection().z()
602 << G4endl;
603 G4cout << " Kinetic Energy : "
604 << std::setw(20) << G4BestUnit(fTrack->GetKineticEnergy(), "Energy")
605 << G4endl;
606 G4cout << " Polarization - x : "
607 << std::setw(20) << fTrack->GetPolarization().x()
608 << G4endl;
609 G4cout << " Polarization - y : "
610 << std::setw(20) << fTrack->GetPolarization().y()
611 << G4endl;
612 G4cout << " Polarization - z : "
613 << std::setw(20) << fTrack->GetPolarization().z()
614 << G4endl;
615 G4cout << " Track Length : "
616 << std::setw(20) << G4BestUnit(fTrack->GetTrackLength(), "Length")
617 << G4endl;
618 G4cout << " Track ID # : "
619 << std::setw(20) << fTrack->GetTrackID()
620 << G4endl;
621 G4cout << " Parent Track ID # : "
622 << std::setw(20) << fTrack->GetParentID()
623 << G4endl;
624 G4cout << " Next Volume : "
625 << std::setw(20);
626 if (fTrack->GetNextVolume() != 0)
627 G4cout << fTrack->GetNextVolume()->GetName() << " ";
628 else G4cout << "OutOfWorld" << " ";
629 G4cout << G4endl;
630 G4cout << " Track Status : "
631 << std::setw(20);
632 if( fTrack->GetTrackStatus() == fAlive )
633 G4cout << " Alive";
634 else if( fTrack->GetTrackStatus() == fStopButAlive )
635 G4cout << " StopButAlive";
636 else if( fTrack->GetTrackStatus() == fStopAndKill )
637 G4cout << " StopAndKill";
638 else if( fTrack->GetTrackStatus() == fKillTrackAndSecondaries )
639 G4cout << " KillTrackAndSecondaries";
640 else if( fTrack->GetTrackStatus() == fSuspend )
641 G4cout << " Suspend";
642 else if( fTrack->GetTrackStatus() == fPostponeToNextEvent )
643 G4cout << " PostponeToNextEvent";
644 G4cout << G4endl;
645 G4cout << " Vertex - x : "
646 << std::setw(20)
647 << G4BestUnit(fTrack->GetVertexPosition().x(),"Length")
648 << G4endl;
649 G4cout << " Vertex - y : "
650 << std::setw(20)
651 << G4BestUnit(fTrack->GetVertexPosition().y(),"Length")
652 << G4endl;
653 G4cout << " Vertex - z : "
654 << std::setw(20)
655 << G4BestUnit(fTrack->GetVertexPosition().z(),"Length")
656 << G4endl;
657 G4cout << " Vertex - Px (MomDir): "
658 << std::setw(20) << fTrack->GetVertexMomentumDirection().x()
659 << G4endl;
660 G4cout << " Vertex - Py (MomDir): "
661 << std::setw(20) << fTrack->GetVertexMomentumDirection().y()
662 << G4endl;
663 G4cout << " Vertex - Pz (MomDir): "
664 << std::setw(20) << fTrack->GetVertexMomentumDirection().z()
665 << G4endl;
666 G4cout << " Vertex - KineE : "
667 << std::setw(20)
668 << G4BestUnit(fTrack->GetVertexKineticEnergy(),"Energy")
669 << G4endl;
670
671 G4cout << " Creator Process : "
672 << std::setw(20);
673 if (fTrack->GetCreatorProcess() == 0)
674 G4cout << " Event Generator" << G4endl;
675 else
676 G4cout << fTrack->GetCreatorProcess()->GetProcessName() << G4endl;
677
678 G4cout << " -----------------------------------------------" << G4endl;
679 G4cout.precision(oldprec);
680}
681
682//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
683
684void G4SteppingVerboseWithUnits::VerboseParticleChange()
685{
686 G4cout << G4endl;
687 G4cout << " ++G4ParticleChange Information " << G4endl;
688 fParticleChange->DumpInfo();
689}
690
691//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
692
693void G4SteppingVerboseWithUnits::ShowStep() const
694{
695 // Show header
696 G4cout << G4endl;
697 G4cout << " ++G4Step Information " << G4endl;
698 G4long oldprc = G4cout.precision(fprec);
699
700 // Show G4Step specific information
701 G4cout << " Address of G4Track : "
702 << fStep->GetTrack() << G4endl;
703 G4cout << " Step Length : "
704 << G4BestUnit(fStep->GetTrack()->GetStepLength(), "Length") << G4endl;
705 G4cout << " Energy Deposit : "
706 << G4BestUnit(fStep->GetTotalEnergyDeposit(), "Energy") << G4endl;
707
708 // Show G4StepPoint specific information
709 G4cout << " -------------------------------------------------------"
710 << "----------------" << G4endl;
711 G4cout << " StepPoint Information "
712 << std::setw(20) << "PreStep" << std::setw(20) << "PostStep" << G4endl;
713 G4cout << " -------------------------------------------------------"
714 << "----------------" << G4endl;
715 G4cout << " Position - x : "
716 << std::setw(17)
717 << G4BestUnit(fStep->GetPreStepPoint()->GetPosition().x(), "Length")
718 << std::setw(17)
719 << G4BestUnit(fStep->GetPostStepPoint()->GetPosition().x(), "Length")
720 << G4endl;
721 G4cout << " Position - y : "
722 << std::setw(17)
723 << G4BestUnit(fStep->GetPreStepPoint()->GetPosition().y(), "Length")
724 << std::setw(17)
725 << G4BestUnit(fStep->GetPostStepPoint()->GetPosition().y(), "Length")
726 << G4endl;
727 G4cout << " Position - z : "
728 << std::setw(17)
729 << G4BestUnit(fStep->GetPreStepPoint()->GetPosition().z(), "Length")
730 << std::setw(17)
731 << G4BestUnit(fStep->GetPostStepPoint()->GetPosition().z(), "Length")
732 << G4endl;
733 G4cout << " Global Time : "
734 << std::setw(17)
735 << G4BestUnit(fStep->GetPreStepPoint()->GetGlobalTime(), "Time")
736 << std::setw(17)
737 << G4BestUnit(fStep->GetPostStepPoint()->GetGlobalTime(), "Time")
738 << G4endl;
739 G4cout << " Local Time : "
740 << std::setw(17)
741 << G4BestUnit(fStep->GetPreStepPoint()->GetLocalTime(), "Time")
742 << std::setw(17)
743 << G4BestUnit(fStep->GetPostStepPoint()->GetLocalTime(), "Time")
744 << G4endl;
745 G4cout << " Proper Time : "
746 << std::setw(17)
747 << G4BestUnit(fStep->GetPreStepPoint()->GetProperTime(), "Time")
748 << std::setw(17)
749 << G4BestUnit(fStep->GetPostStepPoint()->GetProperTime(), "Time")
750 << G4endl;
751 G4cout << " Momentum Direct - x : "
752 << std::setw(20) << fStep->GetPreStepPoint()->GetMomentumDirection().x()
753 << std::setw(20) << fStep->GetPostStepPoint()->GetMomentumDirection().x()
754 << G4endl;
755 G4cout << " Momentum Direct - y : "
756 << std::setw(20) << fStep->GetPreStepPoint()->GetMomentumDirection().y()
757 << std::setw(20) << fStep->GetPostStepPoint()->GetMomentumDirection().y()
758 << G4endl;
759 G4cout << " Momentum Direct - z : "
760 << std::setw(20) << fStep->GetPreStepPoint()->GetMomentumDirection().z()
761 << std::setw(20) << fStep->GetPostStepPoint()->GetMomentumDirection().z()
762 << G4endl;
763 G4cout << " Momentum - x : "
764 << std::setw(14)
765 << G4BestUnit(fStep->GetPreStepPoint()->GetMomentum().x(),"Momentum")
766 << std::setw(14)
767 << G4BestUnit(fStep->GetPostStepPoint()->GetMomentum().x(),"Momentum")
768 << G4endl;
769 G4cout << " Momentum - y : "
770 << std::setw(14)
771 << G4BestUnit(fStep->GetPreStepPoint()->GetMomentum().y(),"Momentum")
772 << std::setw(14)
773 << G4BestUnit(fStep->GetPostStepPoint()->GetMomentum().y(),"Momentum")
774 << G4endl;
775 G4cout << " Momentum - z : "
776 << std::setw(14)
777 << G4BestUnit(fStep->GetPreStepPoint()->GetMomentum().z(),"Momentum")
778 << std::setw(14)
779 << G4BestUnit(fStep->GetPostStepPoint()->GetMomentum().z(),"Momentum")
780 << G4endl;
781 G4cout << " Total Energy : "
782 << std::setw(16)
783 << G4BestUnit(fStep->GetPreStepPoint()->GetTotalEnergy(),"Energy")
784 << std::setw(16)
785 << G4BestUnit(fStep->GetPostStepPoint()->GetTotalEnergy(),"Energy")
786 << G4endl;
787 G4cout << " Kinetic Energy : "
788 << std::setw(16)
789 << G4BestUnit(fStep->GetPreStepPoint()->GetKineticEnergy(),"Energy")
790 << std::setw(16)
791 << G4BestUnit(fStep->GetPostStepPoint()->GetKineticEnergy(),"Energy")
792 << G4endl;
793 G4cout << " Velocity : "
794 << std::setw(14)
795 << G4BestUnit(fStep->GetPreStepPoint()->GetVelocity(),"Velocity")
796 << std::setw(14)
797 << G4BestUnit(fStep->GetPostStepPoint()->GetVelocity(),"Velocity")
798 << G4endl;
799 G4cout << " Volume Name : "
800 << std::setw(20)
801 << fStep->GetPreStepPoint()->GetPhysicalVolume()->GetName();
802 G4String volName = "OutOfWorld";
803 if (fStep->GetPostStepPoint()->GetPhysicalVolume())
804 volName = fStep->GetPostStepPoint()->GetPhysicalVolume()->GetName();
805 G4cout << std::setw(20) << volName
806 << G4endl;
807 G4cout << " Safety : "
808 << std::setw(17)
809 << G4BestUnit(fStep->GetPreStepPoint()->GetSafety(),"Length")
810 << std::setw(17)
811 << G4BestUnit(fStep->GetPostStepPoint()->GetSafety(),"Length")
812 << G4endl;
813 G4cout << " Polarization - x : "
814 << std::setw(20) << fStep->GetPreStepPoint()->GetPolarization().x()
815 << std::setw(20) << fStep->GetPostStepPoint()->GetPolarization().x()
816 << G4endl;
817 G4cout << " Polarization - y : "
818 << std::setw(20) << fStep->GetPreStepPoint()->GetPolarization().y()
819 << std::setw(20) << fStep->GetPostStepPoint()->GetPolarization().y()
820 << G4endl;
821 G4cout << " Polarization - Z : "
822 << std::setw(20) << fStep->GetPreStepPoint()->GetPolarization().z()
823 << std::setw(20) << fStep->GetPostStepPoint()->GetPolarization().z()
824 << G4endl;
825 G4cout << " Weight : "
826 << std::setw(20) << fStep->GetPreStepPoint()->GetWeight()
827 << std::setw(20) << fStep->GetPostStepPoint()->GetWeight()
828 << G4endl;
829 G4cout << " Step Status : " ;
830 G4StepStatus tStepStatus = fStep->GetPreStepPoint()->GetStepStatus();
831 if( tStepStatus == fGeomBoundary )
832 {
833 G4cout << std::setw(20) << "Geom Limit";
834 }
835 else if ( tStepStatus == fAlongStepDoItProc )
836 {
837 G4cout << std::setw(20) << "AlongStep Proc.";
838 }
839 else if ( tStepStatus == fPostStepDoItProc )
840 {
841 G4cout << std::setw(20) << "PostStep Proc";
842 }
843 else if ( tStepStatus == fAtRestDoItProc )
844 {
845 G4cout << std::setw(20) << "AtRest Proc";
846 }
847 else if ( tStepStatus == fUndefined )
848 {
849 G4cout << std::setw(20) << "Undefined";
850 }
851
852 tStepStatus = fStep->GetPostStepPoint()->GetStepStatus();
853 if( tStepStatus == fGeomBoundary )
854 {
855 G4cout << std::setw(20) << "Geom Limit";
856 }
857 else if ( tStepStatus == fAlongStepDoItProc )
858 {
859 G4cout << std::setw(20) << "AlongStep Proc.";
860 }
861 else if ( tStepStatus == fPostStepDoItProc )
862 {
863 G4cout << std::setw(20) << "PostStep Proc";
864 }
865 else if ( tStepStatus == fAtRestDoItProc )
866 {
867 G4cout << std::setw(20) << "AtRest Proc";
868 }
869 else if ( tStepStatus == fUndefined )
870 {
871 G4cout << std::setw(20) << "Undefined";
872 }
873
874 G4cout << G4endl;
875 G4cout << " Process defined Step: " ;
876 if( fStep->GetPreStepPoint()->GetProcessDefinedStep() == 0 )
877 {
878 G4cout << std::setw(20) << "Undefined";
879 }
880 else
881 {
882 G4cout << std::setw(20)
883 << fStep->GetPreStepPoint()
884 ->GetProcessDefinedStep()->GetProcessName();
885 }
886 if( fStep->GetPostStepPoint()->GetProcessDefinedStep() == 0)
887 {
888 G4cout << std::setw(20) << "Undefined";
889 }
890 else
891 {
892 G4cout << std::setw(20)
893 << fStep->GetPostStepPoint()
894 ->GetProcessDefinedStep()->GetProcessName();
895 }
896 G4cout.precision(oldprc);
897
898 G4cout << G4endl;
899 G4cout << " -------------------------------------------------------"
900 << "----------------" << G4endl;
901}
902
903//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4State_Idle
@ G4State_Idle
Definition: G4ApplicationState.hh:86
G4State_PreInit
@ G4State_PreInit
Definition: G4ApplicationState.hh:84
StronglyForced
@ StronglyForced
Definition: G4ForceCondition.hh:55
Conditionally
@ Conditionally
Definition: G4ForceCondition.hh:49
ExclusivelyForced
@ ExclusivelyForced
Definition: G4ForceCondition.hh:52
Forced
@ Forced
Definition: G4ForceCondition.hh:44
CandidateForSelection
@ CandidateForSelection
Definition: G4GPILSelection.hh:42
NotCandidateForSelection
@ NotCandidateForSelection
Definition: G4GPILSelection.hh:47
G4StepStatus
G4StepStatus
Definition: G4StepStatus.hh:40
fGeomBoundary
@ fGeomBoundary
Definition: G4StepStatus.hh:43
fWorldBoundary
@ fWorldBoundary
Definition: G4StepStatus.hh:41
fUndefined
@ fUndefined
Definition: G4StepStatus.hh:55
fPostStepDoItProc
@ fPostStepDoItProc
Definition: G4StepStatus.hh:49
fAtRestDoItProc
@ fAtRestDoItProc
Definition: G4StepStatus.hh:45
fAlongStepDoItProc
@ fAlongStepDoItProc
Definition: G4StepStatus.hh:47
G4SteppingVerboseWithUnits.hh
Definition of the G4SteppingVerboseWithUnits class.
G4BestUnit
#define G4BestUnit(a, b)
Definition: G4SteppingVerbose.cc:45
fKillTrackAndSecondaries
@ fKillTrackAndSecondaries
Definition: G4TrackStatus.hh:48
fSuspend
@ fSuspend
Definition: G4TrackStatus.hh:50
fAlive
@ fAlive
Definition: G4TrackStatus.hh:42
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
fStopButAlive
@ fStopButAlive
Definition: G4TrackStatus.hh:44
fPostponeToNextEvent
@ fPostponeToNextEvent
Definition: G4TrackStatus.hh:52
G4long
long G4long
Definition: G4Types.hh:87
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
CLHEP::Hep3Vector::z
double z() const
CLHEP::Hep3Vector::x
double x() const
CLHEP::Hep3Vector::y
double y() const
G4GenericMessenger::DeclareProperty
Command & DeclareProperty(const G4String &name, const G4AnyType &variable, const G4String &doc="")
Definition: G4GenericMessenger.cc:84
G4StepPoint::GetTotalEnergy
G4double GetTotalEnergy() const
G4StepPoint::GetStepStatus
G4StepStatus GetStepStatus() const
G4StepPoint::GetVelocity
G4double GetVelocity() const
G4StepPoint::GetProperTime
G4double GetProperTime() const
G4StepPoint::GetGlobalTime
G4double GetGlobalTime() const
G4StepPoint::GetSafety
G4double GetSafety() const
G4StepPoint::GetProcessDefinedStep
const G4VProcess * GetProcessDefinedStep() const
G4StepPoint::GetMomentum
G4ThreeVector GetMomentum() const
G4StepPoint::GetPosition
const G4ThreeVector & GetPosition() const
G4StepPoint::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4StepPoint::GetLocalTime
G4double GetLocalTime() const
G4StepPoint::GetPhysicalVolume
G4VPhysicalVolume * GetPhysicalVolume() const
G4StepPoint::GetPolarization
const G4ThreeVector & GetPolarization() const
G4StepPoint::GetKineticEnergy
G4double GetKineticEnergy() const
G4StepPoint::GetWeight
G4double GetWeight() const
G4Step::GetTrack
G4Track * GetTrack() const
G4Step::GetSecondaryInCurrentStep
const std::vector< const G4Track * > * GetSecondaryInCurrentStep() const
Definition: G4Step.cc:202
G4Step::GetPreStepPoint
G4StepPoint * GetPreStepPoint() const
G4Step::GetStepLength
G4double GetStepLength() const
G4Step::GetTotalEnergyDeposit
G4double GetTotalEnergyDeposit() const
G4Step::GetPostStepPoint
G4StepPoint * GetPostStepPoint() const
G4SteppingVerboseWithUnits::SetManager
void SetManager(G4SteppingManager *const) override
Definition: G4SteppingVerboseWithUnits.cc:60
G4Track::GetTrackStatus
G4TrackStatus GetTrackStatus() const
G4Track::GetVertexKineticEnergy
G4double GetVertexKineticEnergy() const
G4Track::GetTrackID
G4int GetTrackID() const
G4Track::GetVolume
G4VPhysicalVolume * GetVolume() const
G4Track::GetCreatorProcess
const G4VProcess * GetCreatorProcess() const
G4Track::GetNextVolume
G4VPhysicalVolume * GetNextVolume() const
G4Track::GetPosition
const G4ThreeVector & GetPosition() const
G4Track::GetTrackLength
G4double GetTrackLength() const
G4Track::GetVertexMomentumDirection
const G4ThreeVector & GetVertexMomentumDirection() const
G4Track::GetGlobalTime
G4double GetGlobalTime() const
G4Track::GetCurrentStepNumber
G4int GetCurrentStepNumber() const
G4Track::GetVertexPosition
const G4ThreeVector & GetVertexPosition() const
G4Track::GetLocalTime
G4double GetLocalTime() const
G4Track::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4Track::GetPolarization
const G4ThreeVector & GetPolarization() const
G4Track::GetStepLength
G4double GetStepLength() const
G4Track::GetParentID
G4int GetParentID() const
G4VParticleChange::DumpInfo
virtual void DumpInfo() const
Definition: G4VParticleChange.cc:133
G4VPhysicalVolume::GetName
const G4String & GetName() const
G4VProcess::GetProcessName
const G4String & GetProcessName() const
Definition: G4VProcess.hh:386
G4VSteppingVerbose::MAXofAlongStepLoops
std::size_t MAXofAlongStepLoops
Definition: G4VSteppingVerbose.hh:152
G4VSteppingVerbose::fCondition
G4ForceCondition fCondition
Definition: G4VSteppingVerbose.hh:187
G4VSteppingVerbose::fGPILSelection
G4GPILSelection fGPILSelection
Definition: G4VSteppingVerbose.hh:188
G4VSteppingVerbose::fParticleChange
G4VParticleChange * fParticleChange
Definition: G4VSteppingVerbose.hh:130
G4VSteppingVerbose::fSelectedPostStepDoItVector
G4SelectedPostStepDoItVector * fSelectedPostStepDoItVector
Definition: G4VSteppingVerbose.hh:178
G4VSteppingVerbose::fSelectedAtRestDoItVector
G4SelectedAtRestDoItVector * fSelectedAtRestDoItVector
Definition: G4VSteppingVerbose.hh:176
G4VSteppingVerbose::fCurrentProcess
G4VProcess * fCurrentProcess
Definition: G4VSteppingVerbose.hh:139
G4VSteppingVerbose::fManager
G4SteppingManager * fManager
Definition: G4VSteppingVerbose.hh:114
G4VSteppingVerbose::MAXofPostStepLoops
std::size_t MAXofPostStepLoops
Definition: G4VSteppingVerbose.hh:153
G4VSteppingVerbose::fSecondary
G4TrackVector * fSecondary
Definition: G4VSteppingVerbose.hh:132
G4GenericMessenger::Command::SetStates
Command & SetStates(G4ApplicationState s0)
Definition: G4GenericMessenger.hh:81