def summarise_event_df(event_df):
summary = {}
for param in event_df.columns:
summary[f"{param}_lower"], summary[f"{param}_upper"], summary[param] = \
summarise_samples(event_df[param])
for param_type in ["_lower", "_upper", ""]:
summary[f'redshift{param_type}'] = \
float(luminosity_distance_to_redshift(
summary[f'luminosity_distance{param_type}']))
summary = {k: summary.get(k, None) for k in GWOSC_KEYS}
summary['version'] = 1
summary['reference'] = "https://github.com/gwastro/2-ogc"
summary['catalog.shortName'] = "PyCBC"
return summary
def summarise_event_df(event_df):
summary = {}
for param in event_df.columns:
summary[f"{param}_lower"], summary[f"{param}_upper"], summary[param] = \
summarise_samples(event_df[param])
for param_type in ["_lower", "_upper", ""]:
summary[f'redshift{param_type}'] = \
float(luminosity_distance_to_redshift(
summary[f'luminosity_distance{param_type}']))
for key in ['mass_1', 'mass_2', 'chirp_mass', 'total_mass']:
summary[f'{key}_source{param_type}'] = \
summary[f"{key}{param_type}"] / (1 + summary[f'redshift{param_type}'])
summary = {k: summary.get(k, None) for k in GWOSC_KEYS}
summary['version'] = 1
summary['reference'] = "https://github.com/jroulet/O2_samples/"
summary['catalog.shortName'] = "IAS"
return summary
def summarise_event(self) -> Dict:
summary = {}
for param in self.samples.columns:
summary[f"{param}_lower"], summary[f"{param}_upper"], summary[param] = \
utils.summarise_samples(self.samples[param])
for param_type in ["_lower", "_upper", ""]:
summary[f'redshift{param_type}'] = \
float(luminosity_distance_to_redshift(
summary[f'luminosity_distance{param_type}']))
for key in ['mass_1', 'mass_2', 'chirp_mass', 'total_mass']:
summary[f'{key}_source{param_type}'] = \
summary[f"{key}{param_type}"] / (
1 + summary[f'redshift{param_type}'])
summary = {k: summary.get(k, None) for k in GWOSC_KEYS}
summary['reference'] = self.datasource
summary['catalog.shortName'] = self.catalog
summary['commonName'] = self.name
return summary
def test_luminosity_redshift_with_cosmology(self):
z = conversion.luminosity_distance_to_redshift(self.distances, cosmology='WMAP9')
dl = conversion.redshift_to_luminosity_distance(z, cosmology='WMAP9')
self.assertAlmostEqual(max(abs(dl - self.distances)), 0, 4)
