Can high dose rates be simulated in GEANT4?
One example to look at would be the eFLASH_radiotherapy advanced example (see also Advanced Examples - Geant4).
Geant4 doesn’t do coupled simulation, so if you need things like, e.g. material changes in response to preceding events that influence future things, that’s not directly possible (I believe there may be ways to tackle this partially, but other should comment).
Ben is correct about this. If you need a coupled simulation (in particular, if you need to modify material composition to add long-lived radiocontamination) you need to do that as a series of consecutive runs:
- Run a sample and collect the radioisotope production rate (i.e., vs. beamOn argument). You may need a fair amount of detail, collecting production per volume, for example.
- Use that output to compute what the new material composition should be, and pass that back into your DetectorConstruction class.
- Create the new materials, and rebuild your geometry.
- Run a new sample; rinse and repeat until you’re happy.
It’s certainly possible to do all of that in a single G4 user application. You could do the intermediate calculations in EndOfRunAction, then set the RunManager “geometry has been modified” flag so that when you start the next run, everything will update.
Alternatively, you can put the “analysis” and “configuration” stuff into a Python script or ROOT macros or something, and do separate executions in a shell script. Which method you choose depends on your own workflow.
Thank you all for your replies, I have benefited greatly.
or you can try the UHDR example (examples/extended/medical/dna/UHDR) for water radiolysis.
This topic was automatically closed 30 days after the last reply. New replies are no longer allowed.