def chirp_volume_montecarlo(found_d, missed_d, found_mchirp, missed_mchirp, distribution_param, distribution, limits_param, min_param, max_param): assert distribution_param == 'chirp_distance' assert limits_param == 'chirp_distance' found_dchirp = chirp_distance(found_d, found_mchirp) missed_dchirp = chirp_distance(missed_d, missed_mchirp) # treat chirp distances in MC volume estimate as physical distances return volume_montecarlo(found_dchirp, missed_dchirp, found_mchirp, missed_mchirp, 'distance', distribution, 'distance', min_param, max_param)
def chirp_distance(dist, mchirp, ref_mass=1.4): return conversions.chirp_distance(dist, mchirp, ref_mass=ref_mass)