Hello @logan12 ,
In case you are simply interested in getting the energy deposition distribution, irregardless of what primary or secondary particle triggered it, it suffices to use the DETECT card and set the energy limits to the signal range you indicated (10 to 200keV). Using the same simulation, simply count the amount of inelastic events by indeed looking at the star density (and taking into account the volume of the regions which you need to calculate yourself).
The total number of collisions can indeed be obtained by adding all the events, also considering different simulations runs. Of course here you need to take into account the amount of primary particles you simulated, hence you need the neutron fluence that you want to scale your results to.
For the inelastic collisions and subsequent processes which you would like information on, I would advise you to look at this and this on the use of mgdraw.f, in addition to the manual, particularly the entries ENDRAW (for energy deposition events) and USDRAW (selecting inelastic collisions by using ICODE=101). In case you want to select only events caused by neutrons you need to select JTRACK.EQ.8. It is up to the user to then customize the routine to ask for more desired information.
Hope this helps,
Andreas