I am interested in implemeting my own (low-energy) nuclear reaction model for a specific projectile-target combination. Let’s say 3He+7Li. This reaction has several open channels. Let’s pretend that I’m only interested in modelling two of these, say, (1) p+9Be and (2) 4He+6Li. Finally, let’s assume that I know the cross sections of these two channels, cs1 and cs2, and their angular distributions.
Based on my (somewhat limited) understanding of GEANT4, I believe this could be implemented as follows,
i) Create a new physical process, e.g. called G4He3Li7NuclearReaction, which inherits from G4HadronicProcess
ii) Create a new cross section dataset with the total (cs1+cs2) cross section. This class should inherit from G4VCrossSectionDataSet.
iii) Add the dataset from ii) to the physics process from i)
iv) Create a new interaction model that inherits from G4HadronicInteraction. This interaction model should model both channels. In particular, the interaction model would use the cross section ratio cs1/cs2 to determine on an event-by-event basis through monte-carlo samping if the reaction proceeds via channel (1) or (2).
Does this seem like a sensible approach to achieving the desired goal?