change mass distributions' variates to "massX_source" instead of just "mass"

This way, we don't have to worry about the annoying fact that currently we draw in "mass" but then rename those attributes when we compute detector-frame masses. This complicates downstream processing and would be better to be avoided if at all possible.