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