163 const G4StepPoint* preStep =
step->GetPreStepPoint();
164 const G4StepPoint* postStep =
step->GetPostStepPoint();
177 prec =
G4cout.precision(16);
180 G4cout <<
" ++G4Step Information " << G4endl;
181 G4cout <<
" Address of G4Track : " <<
track << G4endl;
182 G4cout <<
" Step Length (mm) : " <<
track->GetStepLength() << G4endl;
183 G4cout <<
" Energy Deposit (MeV) : " <<
step->GetTotalEnergyDeposit() << G4endl;
185 G4cout <<
" -------------------------------------------------------"
186 <<
"-------------------------------" << G4endl;
187 G4cout <<
" StepPoint Information " << std::setw(30) <<
"PreStep" << std::setw(30) <<
"PostStep" << G4endl;
188 G4cout <<
" -------------------------------------------------------"
189 <<
"-------------------------------" << G4endl;
190 G4cout <<
" Position - x (cm) : " << std::setw(30) << preStep->GetPosition().x() / CLHEP::cm << std::setw(30)
191 << postStep->GetPosition().x() / CLHEP::cm << G4endl;
192 G4cout <<
" Position - y (cm) : " << std::setw(30) << preStep->GetPosition().y() / CLHEP::cm << std::setw(30)
193 << postStep->GetPosition().y() / CLHEP::cm << G4endl;
194 G4cout <<
" Position - z (cm) : " << std::setw(30) << preStep->GetPosition().z() / CLHEP::cm << std::setw(30)
195 << postStep->GetPosition().z() / CLHEP::cm << G4endl;
196 G4cout <<
" Global Time (ns) : " << std::setw(30) << preStep->GetGlobalTime() / CLHEP::ns << std::setw(30)
197 << postStep->GetGlobalTime() / CLHEP::ns << G4endl;
198 G4cout <<
" Local Time (ns) : " << std::setw(30) << preStep->GetLocalTime() / CLHEP::ns << std::setw(30)
199 << postStep->GetLocalTime() / CLHEP::ns << G4endl;
200 G4cout <<
" Proper Time (ns) : " << std::setw(30) << preStep->GetProperTime() / CLHEP::ns << std::setw(30)
201 << postStep->GetProperTime() / CLHEP::ns << G4endl;
202 G4cout <<
" Momentum Direct - x : " << std::setw(30) << preStep->GetMomentumDirection().x() << std::setw(30)
203 << postStep->GetMomentumDirection().x() << G4endl;
204 G4cout <<
" Momentum Direct - y : " << std::setw(30) << preStep->GetMomentumDirection().y() << std::setw(30)
205 << postStep->GetMomentumDirection().y() << G4endl;
206 G4cout <<
" Momentum Direct - z : " << std::setw(30) << preStep->GetMomentumDirection().z() << std::setw(30)
207 << postStep->GetMomentumDirection().z() << G4endl;
208 G4cout <<
" Momentum - x (GeV/c): " << std::setw(30) << preStep->GetMomentum().x() /
CLHEP::GeV
209 << std::setw(30) << postStep->GetMomentum().x() /
CLHEP::GeV << G4endl;
210 G4cout <<
" Momentum - y (GeV/c): " << std::setw(30) << preStep->GetMomentum().y() /
CLHEP::GeV
211 << std::setw(30) << postStep->GetMomentum().y() /
CLHEP::GeV << G4endl;
212 G4cout <<
" Momentum - z (GeV/c): " << std::setw(30) << preStep->GetMomentum().z() /
CLHEP::GeV
213 << std::setw(30) << postStep->GetMomentum().z() /
CLHEP::GeV << G4endl;
214 G4cout <<
" Total Energy (GeV) : " << std::setw(30) << preStep->GetTotalEnergy() /
CLHEP::GeV << std::setw(30)
215 << postStep->GetTotalEnergy() /
CLHEP::GeV << G4endl;
216 G4cout <<
" Kinetic Energy (GeV): " << std::setw(30) << preStep->GetKineticEnergy() /
CLHEP::GeV
217 << std::setw(30) << postStep->GetKineticEnergy() /
CLHEP::GeV << G4endl;
218 G4cout <<
" Velocity (mm/ns) : " << std::setw(30) << preStep->GetVelocity() << std::setw(30)
219 << postStep->GetVelocity() << G4endl;
220 G4cout <<
" Volume Name : " << std::setw(30) << preStep->GetPhysicalVolume()->GetName();
221 G4String volName = (postStep->GetPhysicalVolume()) ? postStep->GetPhysicalVolume()->GetName() :
"OutOfWorld";
223 G4cout << std::setw(30) << volName << G4endl;
224 G4cout <<
" Safety (mm) : " << std::setw(30) << preStep->GetSafety() << std::setw(30)
225 << postStep->GetSafety() << G4endl;
226 G4cout <<
" Polarization - x : " << std::setw(30) << preStep->GetPolarization().x() << std::setw(30)
227 << postStep->GetPolarization().x() << G4endl;
228 G4cout <<
" Polarization - y : " << std::setw(30) << preStep->GetPolarization().y() << std::setw(30)
229 << postStep->GetPolarization().y() << G4endl;
230 G4cout <<
" Polarization - Z : " << std::setw(30) << preStep->GetPolarization().z() << std::setw(30)
231 << postStep->GetPolarization().z() << G4endl;
232 G4cout <<
" Weight : " << std::setw(30) << preStep->GetWeight() << std::setw(30)
233 << postStep->GetWeight() << G4endl;
234 G4cout <<
" Step Status : ";
235 G4StepStatus tStepStatus = preStep->GetStepStatus();
236 if (tStepStatus == fGeomBoundary) {
237 G4cout << std::setw(30) <<
"Geom Limit";
238 }
else if (tStepStatus == fAlongStepDoItProc) {
239 G4cout << std::setw(30) <<
"AlongStep Proc.";
240 }
else if (tStepStatus == fPostStepDoItProc) {
241 G4cout << std::setw(30) <<
"PostStep Proc";
242 }
else if (tStepStatus == fAtRestDoItProc) {
243 G4cout << std::setw(30) <<
"AtRest Proc";
244 }
else if (tStepStatus == fUndefined) {
245 G4cout << std::setw(30) <<
"Undefined";
248 tStepStatus = postStep->GetStepStatus();
249 if (tStepStatus == fGeomBoundary) {
250 G4cout << std::setw(30) <<
"Geom Limit";
251 }
else if (tStepStatus == fAlongStepDoItProc) {
252 G4cout << std::setw(30) <<
"AlongStep Proc.";
253 }
else if (tStepStatus == fPostStepDoItProc) {
254 G4cout << std::setw(30) <<
"PostStep Proc";
255 }
else if (tStepStatus == fAtRestDoItProc) {
256 G4cout << std::setw(30) <<
"AtRest Proc";
257 }
else if (tStepStatus == fUndefined) {
258 G4cout << std::setw(30) <<
"Undefined";
262 G4cout <<
" Process defined Step: ";
263 if (preStep->GetProcessDefinedStep() ==
nullptr) {
264 G4cout << std::setw(30) <<
"Undefined";
266 G4cout << std::setw(30) << preStep->GetProcessDefinedStep()->GetProcessName();
268 if (postStep->GetProcessDefinedStep() ==
nullptr) {
269 G4cout << std::setw(30) <<
"Undefined";
271 G4cout << std::setw(30) << postStep->GetProcessDefinedStep()->GetProcessName();
276 G4cout <<
" -------------------------------------------------------"
277 <<
"-------------------------------" << G4endl;
282 const G4VTouchable* tch1 = preStep->GetTouchable();
283 const G4VTouchable* tch2 = postStep->GetTouchable();
285 G4double
x = postStep->GetPosition().x();
286 G4double
y = postStep->GetPosition().y();
287 G4cout <<
"Touchable depth pre= " << tch1->GetHistoryDepth() <<
" post= " << tch2->GetHistoryDepth()
288 <<
" trans1= " << tch1->GetTranslation(tch1->GetHistoryDepth())
289 <<
" trans2= " << tch2->GetTranslation(tch2->GetHistoryDepth()) <<
" r= " <<
std::sqrt(
x *
x +
y *
y)
291 const G4VPhysicalVolume* pv1 = preStep->GetPhysicalVolume();
292 const G4VPhysicalVolume* pv2 = postStep->GetPhysicalVolume();
293 G4cout <<
"PreStepVolume: " << pv1->GetName() << G4endl;
294 G4cout <<
" Translation: " << pv1->GetObjectTranslation() << G4endl;
295 G4cout <<
" Rotation: " << pv1->GetObjectRotationValue() << G4endl;
296 const G4VSolid* sv1 = pv1->GetLogicalVolume()->GetSolid();
300 G4cout <<
"PostStepVolume: " << pv2->GetName() << G4endl;
301 G4cout <<
" Translation: " << pv2->GetObjectTranslation() << G4endl;
302 G4cout <<
" Rotation: " << pv2->GetObjectRotationValue() << G4endl;
303 const G4VSolid* sv2 = pv2->GetLogicalVolume()->GetSolid();
309 G4cout <<
"##### Geometry Depth Analysis" << G4endl;
310 for (G4int
k = 1;
k < tch1->GetHistoryDepth(); ++
k) {
311 const G4VPhysicalVolume*
pv = tch1->GetVolume(
k);
313 const G4VSolid* sol =
pv->GetLogicalVolume()->GetSolid();
314 G4cout <<
" Depth # " <<
k <<
" PhysVolume " <<
pv->GetName() << G4endl;
315 G4cout <<
" Translation: " <<
pv->GetObjectTranslation() << G4endl;
316 G4cout <<
" Rotation: " <<
pv->GetObjectRotationValue() << G4endl;
325 prec =
G4cout.precision(4);
326 G4cout << std::setw(5) <<
track->GetCurrentStepNumber() <<
" " << std::setw(8) <<
track->GetPosition().x() / CLHEP::cm
327 <<
" " << std::setw(8) <<
track->GetPosition().y() / CLHEP::cm <<
" " << std::setw(8)
328 <<
track->GetPosition().z() / CLHEP::cm <<
" " << std::setw(9) <<
track->GetKineticEnergy() /
CLHEP::GeV
331 prec =
G4cout.precision(3);
333 <<
step->GetStepLength() / CLHEP::mm <<
" " << std::setw(9) <<
track->GetTrackLength() / CLHEP::cm <<
" ";
335 G4bool endTracking =
false;
336 if (
track->GetNextVolume() !=
nullptr) {
337 G4cout << std::setw(30) <<
track->GetVolume()->GetName() <<
" ";
339 G4cout << std::setw(30) <<
"OutOfWorld"
346 }
else if (postStep->GetProcessDefinedStep() !=
nullptr) {
347 G4cout << postStep->GetProcessDefinedStep()->GetProcessName();
356 if (!endTracking && fStopAndKill !=
track->GetTrackStatus()) {
360 prec =
G4cout.precision(4);
361 const G4TrackVector* tv =
step->GetSecondary();
362 G4int
nt = tv->size();
364 G4cout <<
" ++List of " <<
nt <<
" secondaries generated; "
367 for (G4int
i = 0;
i <
nt; ++
i) {
371 G4cout <<
" (" << std::setw(9) << (*tv)[
i]->GetPosition().x() / CLHEP::cm <<
" " << std::setw(9)
372 << (*tv)[
i]->GetPosition().y() / CLHEP::cm <<
" " << std::setw(9) << (*tv)[
i]->GetPosition().z() / CLHEP::cm
373 <<
") cm, " << std::setw(9) << (*tv)[
i]->GetKineticEnergy() /
CLHEP::GeV <<
" GeV, " << std::setw(9)
374 << (*tv)[
i]->GetGlobalTime() / CLHEP::ns <<
" ns, " << std::setw(18)
375 << (*tv)[
i]->GetDefinition()->GetParticleName() << G4endl;