Dear @rachel.hen.shukrun,
The difference between EVENTBIN and mgdraw.f is that EVENTBIN will already impose a scoring mesh (i.e., division of your scoring volume in cells) whereas mgdraw.f does not. However, the mgdraw.f routine allows you to distinguish between continuous and discrete energy depositions as mentioned above. It also allows you to record only the energy deposition events in the energy windows you want to look at.
For both methods, the merging of the results is the responsibility of the user. Each entry you record corresponds to the energy deposition caused by the primary particle and all of its secondaries. You will need to calculate the average and standard deviation of all the entries you record (using either EVENTBIN or mgdraw.f) to get an estimate of the value per primary and error.
Cheers,
Andreas