Geant4 Cross Reference |
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 // This example is provided by the Geant4-DNA collaboration 27 // Any report or published results obtained using the Geant4-DNA software 28 // shall cite the following Geant4-DNA collaboration publications: 29 // Med. Phys. 45 (2018) e722-e739 30 // Phys. Med. 31 (2015) 861-874 31 // Med. Phys. 37 (2010) 4692-4708 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178 33 // 34 // The Geant4-DNA web site is available at http://geant4-dna.org 35 // 36 /// \file SteppingVerbose.cc 37 /// \brief Implementation of the SteppingVerbose class 38 39 #include "SteppingVerbose.hh" 40 41 #include "G4UnitsTable.hh" 42 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 44 45 SteppingVerbose::SteppingVerbose() : G4SteppingVerbose() {} 46 47 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 48 49 SteppingVerbose::~SteppingVerbose() {} 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 53 void SteppingVerbose::TrackingStarted() 54 { 55 CopyState(); 56 57 G4int prec = G4cout.precision(3); 58 59 // Step zero 60 if (verboseLevel > 0) { 61 G4cout << std::setw(5) << "Step#" 62 << " " << std::setw(6) << "X" 63 << " " << std::setw(6) << "Y" 64 << " " << std::setw(6) << "Z" 65 << " " << std::setw(9) << "KineE" 66 << " " << std::setw(9) << "dEStep" 67 << " " << std::setw(10) << "StepLeng" << std::setw(10) << "TrakLeng" << std::setw(10) 68 << "Volume" 69 << " " << std::setw(10) << "Process" << G4endl; 70 71 G4cout << std::setw(5) << fTrack->GetCurrentStepNumber() << " " << std::setw(6) 72 << G4BestUnit(fTrack->GetPosition().x(), "Length") << std::setw(6) 73 << G4BestUnit(fTrack->GetPosition().y(), "Length") << std::setw(6) 74 << G4BestUnit(fTrack->GetPosition().z(), "Length") << std::setw(6) 75 << G4BestUnit(fTrack->GetKineticEnergy(), "Energy") << std::setw(6) 76 << G4BestUnit(fStep->GetTotalEnergyDeposit(), "Energy") << std::setw(6) 77 << G4BestUnit(fStep->GetStepLength(), "Length") << std::setw(6) 78 << G4BestUnit(fTrack->GetTrackLength(), "Length") << std::setw(10) 79 << fTrack->GetVolume()->GetName() << " initStep" << G4endl; 80 } 81 G4cout.precision(prec); 82 } 83 84 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 85 86 void SteppingVerbose::StepInfo() 87 { 88 CopyState(); 89 90 G4int prec = G4cout.precision(3); 91 92 if (verboseLevel >= 1) { 93 if (verboseLevel >= 4) VerboseTrack(); 94 if (verboseLevel >= 3) { 95 G4cout << G4endl; 96 G4cout << std::setw(5) << "#Step#" 97 << " " << std::setw(6) << "X" 98 << " " << std::setw(6) << "Y" 99 << " " << std::setw(6) << "Z" 100 << " " << std::setw(9) << "KineE" 101 << " " << std::setw(9) << "dEStep" 102 << " " << std::setw(10) << "StepLeng" << std::setw(10) << "TrakLeng" << std::setw(10) 103 << "Volume" 104 << " " << std::setw(10) << "Process" << G4endl; 105 } 106 107 G4cout << std::setw(5) << fTrack->GetCurrentStepNumber() << " " << std::setw(6) 108 << G4BestUnit(fTrack->GetPosition().x(), "Length") << std::setw(6) 109 << G4BestUnit(fTrack->GetPosition().y(), "Length") << std::setw(6) 110 << G4BestUnit(fTrack->GetPosition().z(), "Length") << std::setw(6) 111 << G4BestUnit(fTrack->GetKineticEnergy(), "Energy") << std::setw(6) 112 << G4BestUnit(fStep->GetTotalEnergyDeposit(), "Energy") << std::setw(6) 113 << G4BestUnit(fStep->GetStepLength(), "Length") << std::setw(6) 114 << G4BestUnit(fTrack->GetTrackLength(), "Length") << std::setw(10) 115 << fTrack->GetVolume()->GetName(); 116 117 const G4VProcess* process = fStep->GetPostStepPoint()->GetProcessDefinedStep(); 118 G4String procName = " UserLimit"; 119 if (process) procName = process->GetProcessName(); 120 if (fStepStatus == fWorldBoundary) procName = "OutOfWorld"; 121 G4cout << " " << std::setw(10) << procName; 122 G4cout << G4endl; 123 124 if (verboseLevel == 2) { 125 const std::vector<const G4Track*>* secondary = fStep->GetSecondaryInCurrentStep(); 126 size_t nbtrk = (*secondary).size(); 127 if (nbtrk) { 128 G4cout << "\n :----- List of secondaries ----------------" << G4endl; 129 G4cout.precision(4); 130 for (size_t lp = 0; lp < (*secondary).size(); lp++) { 131 G4cout << " " << std::setw(13) << (*secondary)[lp]->GetDefinition()->GetParticleName() 132 << ": energy =" << std::setw(6) 133 << G4BestUnit((*secondary)[lp]->GetKineticEnergy(), "Energy") 134 << " time =" << std::setw(6) 135 << G4BestUnit((*secondary)[lp]->GetGlobalTime(), "Time"); 136 G4cout << G4endl; 137 } 138 139 G4cout << " :------------------------------------------\n" << G4endl; 140 } 141 } 142 } 143 G4cout.precision(prec); 144 } 145