Optimization of detector response simulations for multiple particles created by neutron-induced reactions