A histogram-based objective function of the form described above has the computational advantage that if an single beam weight or position is varied, one only need look at the voxels intersecting the old and new beam positions in order to update the value of the objective function. This may provide a significant computational advantage when using optimization algorithms that adjust only a few beams at a time while exploring the space of possible solutions.