sizheng-ma/qnm_filter

Ringdown analysis with rational filters.

qnm_filter

Ringdown analysis with rational filters. The properties of the rational filters can be found in Ma, et al. (2023a), Ma, et al. (2023b), and Ma, et al. (2022).

Installation

Install this Python package in "editable mode":

git clone [email protected]:Sizheng-Ma/qnm_filter.git
cd qnm_filter
pip install -e path/to/this/repository

Examples

1. The analysis of GW150914 can be found in code_demo/Filter.ipynb.
2. Two injection systems with a numerical-relativity waveform (GW150914-like) are available at code_demo/NRInjection_bilby.ipynb and code_demo/NRInjection_customized_noise.ipynb.
3. The evidence curve is a useful tool for determining the statistical significance of a quasinormal mode and for limiting its start time. An example is at code_demo/Evidence.ipynb.

Cite this code

Please cite these three papers if you use this code for academic research:

• Black Hole Spectroscopy by Mode Cleaning, Ma et al., Phys. Rev. Lett. 130, 141401 (2023).
• Using rational filters to uncover the first ringdown overtone in GW150914, Ma et al., Phys. Rev. D 107, 084010 (2023).
• Quasinormal-mode filters: A new approach to analyze the gravitational-wave ringdown of binary black-hole mergers, Ma et al., Phys. Rev. D 106, 084036 (2022).

Documentation

The documentation can be viewed at here.

qnm_filter's Issues

detector_alignment shouldn't take t_init

t_init is already an attribute. The current infrastructure may cause issues.

A bug in noise

While reading a noise file, Noise assumes the frequency is equally sampled.

Fix parallel_compute

1. parallel_compute takes three arguments instead of 2.
2. Custom n_jobs, and test it on clusters.

Move a function to source files

Move likelihood_vs_mass_spin to Network class.

import_data_array

import_data_array might be redundant

pad_data_for_fft

The comment for len_pow in pad_data_for_fft is wrong

A bug in Noise

The noise injections within bilby are not as expected. The sampling_rate of the acf is affecting the SNR

Add comments

add nperseg to get_acf

NRInjection_without_bilby bug

When running the NRInjection_without_bilby.ipynb and seeking to generate a new realisation with a different SNR (i.e. SNR=33 -> SNR=13) by changing luminous_intensity (which is in MPc). If you don't clear the variables in memory errors will occur and the generated waveform will be incorrect.

Not sure why this error occurs and will look into it further