50 {
52
53 const G4ThreeVector& myPos = aStep->GetPostStepPoint()->GetPosition();
58
59 G4TransportationManager *tm = G4TransportationManager::GetTransportationManager();
60 tm->GetNavigatorForTracking()->SetVerboseLevel(
m_config.verboseLevel);
61 tm->GetPropagatorInField()->SetVerboseLevel(
m_config.verboseLevel);
62
63 G4RunManagerKernel *rmk = G4RunManagerKernel::GetRunManagerKernel();
64 rmk->GetTrackingManager()->SetVerboseLevel(
m_config.verboseLevel);
65 rmk->GetTrackingManager()->GetSteppingManager()->SetVerboseLevel(
m_config.verboseLevel);
66 rmk->GetStackManager()->SetVerboseLevel(
m_config.verboseLevel);
67
68 G4Track *tr = aStep->GetTrack();
69 const G4ThreeVector&
mom = tr->GetMomentumDirection();
70
71 std::cout << "Moving " << tr->GetDefinition()->GetParticleName() << " at (" << myPos.x()
72 <<
", " << myPos.y() <<
", " << myPos.z() <<
") to (" <<
mom.x() <<
", " <<
mom.y()
73 <<
", " <<
mom.z() <<
") from ";
74 if (aStep->GetPreStepPoint()->GetPhysicalVolume()) std::cout << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetName();
75 else std::cout << "noVolume";
76 std::cout << " to ";
77 if (aStep->GetPostStepPoint()->GetPhysicalVolume()) std::cout << aStep->GetPostStepPoint()->GetPhysicalVolume()->GetName();
78 else std::cout << "outOfWorld";
79 std::cout << " with KE=" << tr->GetKineticEnergy() << " pT=" << tr->GetMomentum().perp()
80 << " eta=" << tr->GetMomentum().eta() << " length " << aStep->GetStepLength() << " energy "
81 << aStep->GetTotalEnergyDeposit() << " with process ";
82 if (aStep->GetPostStepPoint()->GetProcessDefinedStep()) std::cout << aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();
83 else std::cout << "Unknown";
84 std::cout << " from a ";
85 if (tr->GetCreatorProcess()) std::cout << tr->GetCreatorProcess()->GetProcessName();
86 else if (tr->GetParentID()==0) std::cout << "Generator";
87 else std::cout << "pid="<<tr->GetParentID();
88 if (tr->GetLogicalVolumeAtVertex()) std::cout << " in " << tr->GetLogicalVolumeAtVertex()->GetName() << ", ";
89 else std::cout << " nowhere, ";
90 std::cout << tr->GetTrackLength() << " mm ago" << std::endl;
91
93
94 }
95 }
96 }