49 = G4SDManager::GetSDMpointer()->GetCollectionID(
"crystal/edep");
55 = G4SDManager::GetSDMpointer()->GetCollectionID(
"patient/dose");
59 G4int evtNb =
event->GetEventID();
62 G4cout << G4endl <<
"---> end of event: " << evtNb << G4endl;
72 const G4double eThreshold = 500*keV;
77 std::map<G4int,G4double*>::iterator itr;
78 for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
79 G4double edep = *(itr->second);
80 if (edep > eThreshold) nbOfFired++;
92 for (itr = evtMap->GetMap()->begin(); itr != evtMap->GetMap()->end(); itr++) {
94 dose = *(itr->second);
99 G4Run::RecordEvent(event);