PSD interpolation broken in production
The PSD interpolation needed for calculating higher mode/precessing SNR is broken in the IGWN environments. It looks like there may have been a change in the FrequencySeries API, I'm not sure if this is using gwpy or pycbc.
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/core/cli/inputs.py", line 1996, in __init__
self.samples = self.opts.samples
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/core/cli/inputs.py", line 711, in samples
self._set_samples(samples)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/cli/inputs.py", line 209, in _set_samples
super(_GWInput, self)._set_samples(*args, **kwargs)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/core/cli/inputs.py", line 755, in _set_samples
data = self.grab_data_from_input(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/core/cli/inputs.py", line 1335, in grab_data_from_input
data = self.grab_data_from_file(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/cli/inputs.py", line 190, in grab_data_from_file
data = pesummary.core.cli.inputs._Input.grab_data_from_file(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/core/cli/inputs.py", line 263, in grab_data_from_file
f.generate_all_posterior_samples(**kwargs)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/file/formats/base_read.py", line 455, in generate_all_posterior_samples
data = convert(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/__init__.py", line 147, in convert
return _Conversion(*args, resume_file=resume_file, **kwargs)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/__init__.py", line 418, in __new__
obj.__init__(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/__init__.py", line 563, in __init__
self.generate_all_posterior_samples()
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/__init__.py", line 2025, in generate_all_posterior_samples
self._rho_hm(cond)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/__init__.py", line 1174, in _rho_hm
rho, data_used = multipole_snr(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/utils/decorators.py", line 211, in wrapper_function
output = func(*new_args, **new_kwargs)
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/snr.py", line 546, in multipole_snr
psd, ANALYTIC = _setup_psd(
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/snr.py", line 175, in _setup_psd
psd = {
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pesummary/gw/conversions/snr.py", line 176, in <dictcomp>
ifo: estimate.interpolate(psd[ifo], df) for ifo in psd.keys()
File "/cvmfs/software.igwn.org/conda/envs/igwn-py310-20240410/lib/python3.10/site-packages/pycbc/psd/estimate.py", line 309, in interpolate
interpolated_series = numpy.interp(samples, series.sample_frequencies.numpy(), series.numpy())
AttributeError: 'PSD' object has no attribute 'sample_frequencies'