I’ve recently started using openMC and am blown away by the simplicity and power of the interface. I have an application for which I’d usually use G4/MCNP that requires me to run two simulations.
One with no restrictions - OpenMC is more than capable.
The second, I need to be able to do one of two things:
Either kill the particle/event immediately after the first reaction in my “cell of interest” has occurred.
Or, use the simulation’s output data to examine the type and order of reactions in events that actually had at least one of a specific reaction type.
I’ve examined the documentation and haven’t found a conclusive answer - and so I was wondering if anyone here would be able to provide one!
Essentially, does OpenMC afford me this level of control over each event, or can I configure the settings to provide a sufficiently verbose output such that I can examine which reactions occurred in which event, and in what order?
the standard track information gives x,y,z postion, u,v,w vectors, energy, time, particle weight, cell ID, cell instance, and material ID (see particle track documentation). This is easy to extract, however,
The source code for collision track shows that some additional information that seems relevant to your question is stored, such as the nuclide, parent and progeny ID involved in the collision, but I’m not sure if/how this information can be extracted through the tracks.h5 file.
Furthermore, if you use settings.verbosity = 10 in your simulation, it will output the full sequence of interactions/reactions to the terminal, but again, i don’t know how this information can be extracted.
if the information above can be extracted it would be fairly simple to filter data to look for your reaction of interest.
I hope some of this is useful - hopefully someone with deeper knowledge can provide more info/guidance, I would be very interested in any developments in extracting the extra particle track information.