Suppose I make a cylindrical crystal and assign its material as 60Co (or any radioactive ion). I make this my scoring volume to deposit the energy absorbed.
Just a clarification, this is not my source, as in not defined in the event generator but rather in detector construction. The relevant G4Radioactive and G4Decay classes have been included in the physics list.
When I run the simulation will this material undergo radioactive decay and also contribute to the resulting spectra?
Can someone please share their thoughts on the same?
You found the correct answer. It is possible to use GPS to specify a region of the world (a “volume”, but not a true G4LogicalVolume) within which to generate events. You could set up your GPS so that the source region corresponds to your “detector crystal”, and throw Co-60 decays in that region.
In our experiment, we are very concerned about radiocontaminants in both the support structures and the detectors themselves, and this is the main way we generate simulated events. We wrote code where you can pass in a specific G4 volume by name (either LV or PV), and it picks random positions in the corresponding G4VSolid for each event.
Yes, this approach is quite helpful of “overlapping” source and detector regions is quite helpful.
Just a follow up though are there any examples/codes that achieve the tasks of simulating decays uniformly inside a given region (cylindrical or other geometry)?