found a bug in drawing theta_inj from a unifor distribution over the sky within pareto_hrss
fixed a bug in which work was duplicated in multiprocessing set_* algorithms if num_proc==1
tweaked a bit to get better information and to get a statistical set from which we can draw conclusions
added a reference line to pp plots