| Message: Making multiple sensitive Detectors - a probelm | Not Logged In (login) |
|
Hi, I am simulating a compton camera arrangement and so have two detectors with gamma interactions, (currently set at <1MeV incoming). I have tried to set the two detectors to be sensitive, which I believe works, but when I output the hit information I find something strange:
--->Hits Collection: in this event there are 0 hits in the detector : BUT SHOWS TWO HITS trackID: 1 time: 325.428 ps SPDetector No: 0 energy dep: 11.067 keV trackID: 2 time: 325.492 ps SPDetector No: 0 energy dep: 134.739 keV
--->Hits Collection: in this event there are 2 hits in the detector : trackID: 1 time: 325.428 ps SPDetector No: 0 energy dep: 11.067 keV trackID: 2 time: 325.492 ps SPDetector No: 0 energy dep: 134.739 keV
--->Hits Collection: in this event there are 0 hits in the detector : BUT SHOWS TWO HITS trackID: 1 time: 162.511 ps SPDetector No: 0 energy dep: 11.067 keV trackID: 2 time: 162.581 ps SPDetector No: 0 energy dep: 145.671 keV
--->Hits Collection: in this event there are 2 hits in the detector : trackID: 1 time: 162.511 ps SPDetector No: 0 energy dep: 11.067 keV trackID: 2 time: 162.581 ps SPDetector No: 0 energy dep: 145.671 keV
--->Hits Collection: in this event there are 0 hits in the detector : BUT SHOWS THREE HITS trackID: 1 time: 183.959 ps SPDetector No: 0 energy dep: 129.166 keV trackID: 1 time: 205.561 ps SPDetector No: 0 energy dep: 61.9996 keV trackID: 1 time: 208.869 ps SPDetector No: 0 energy dep: 14.6158 keV
--->Hits Collection: in this event there are 3 hits in the detector : trackID: 1 time: 183.959 ps SPDetector No: 0 energy dep: 129.166 keV trackID: 1 time: 205.561 ps SPDetector No: 0 energy dep: 61.9996 keV trackID: 1 time: 208.869 ps SPDetector No: 0 energy dep: 14.6158 keV So I appear to be duplicating the hits collection rather than looking at the detectors independently. I don't understand if this is a problem in my detectorHits file or my sensitivedetector file or if I haven't declared the sensitive detectors correctly:
G4SDManager *SDman = G4SDManager::GetSDMpointer(); SPDetectorSD* aDetectorSD = new SPDetectorSD( detName[i] ); SDman->AddNewDetector( aDetectorSD ); detector_log[i]->SetSensitiveDetector(aDetectorSD); (within a large looped geometry) Any advice would be greatly appreciated, thanks in advance,
Laura
|
| Inline Depth: | Outline Depth: | Add message: |
|
to: |