I’m currently modifying the detector geometry in my simulation and encounter an unexpected exception using /gps/confine. In this new version I include a radioactive spot to which I will try to confine the event generation. This radioactive spot is defined in the DetectorConstruction class as a G4tubs with a height of 1nm and a 1mm radius centered on (0., 0., -1. mm).
In the macro I defined:
/gps/pos/centre 0. 0. -1. mm
/gps/pos/radius 1. mm
Which works as expected (for charged geantinos and 10 events). Keeping everything the same except for a bigger value of /gps/pos/radius, e.g. 1.1 mm, results in G4Exception: G4GPS001.
LoopCount = 100000
Either the source distribution >> confinement or any confining volume may not overlap with the source distribution or any confining volumes may not exist
This exception was unexpected for me, especially when looking at the output when using /gps/verbose 1 where there are several events which according to me fit the requirements to be excepted:
G4WT0 > Rotated and Translated position (0.8435785306149446,-0.5855505638464635,-1)
G4WT0 > Rotated and Translated position (0.04277877150439346,0.1528173779970043,-1) Should be excepted?
G4WT0 > Rotated and Translated position (0.4877390072176273,0.6246413596337022,-1) Should be excepted?
G4WT0 > Rotated and Translated position (0.5520109266627489,-0.9381969957858898,-1)
G4WT0 > Rotated and Translated position (-0.5594148641537221,0.5836322704609065,-1) Should be excepted?
Can anyone help me with this issue? It’s the first time that I use /gps/confine, don’t hesitate to let me know if I need to include additional information.
Thanks in advance!