Decay modes #1
Description
The simulation currently only allows a maximum of 6 decay modes for a given state. Although this will have a marginal affect on the simulated BGO detection efficiency (because the systematic uncertainty of the BGO detection efficiency is nominally ~10%), neglected decay modes can have a drastic affect on the separator transmission, especially when dealing with reactions with a large recoil cone angle.
Currently, decays are handled in the user subroutine ureact.f which then calls the geant3 subroutine gsdk. gsdk contains the arrays BRATIO and MODE each containing 6 elements. This routine should be modified (in the form of a new user routine udecay.f) mimicking gsdk but that handles 10 decay modes and properly addresses the particle pointer.