GWSkyNet : a machine learning algorithm used for estimating the probabiliy that a gravitational wave candidate is of astrophysical origin
Project description
GWSkyNet
GWSkyNet, developed by Cabero et al 2020, is a machine learning classifier capable of distinguishing between astrophysical events and instrumental artifacts. Designed to facilitate potential EM follow-up observations, the classifier can be operated in low-latency and provide information complementary to what has previously been released in an Open Public Alert(OPA) seconds after an alert is published. A unique feature of the classifier is that it needs only the publicly available information from the OPA system. The classifier can also be expanded easily to intake other information such as the SNRs of a GW candidate in each detector in the network.
GWSkyNet is now being developed as a low-latency annotation pipeline for the fourth observing run of the LIGO Scientific and Virgo Collaboration (O4).
For more detail, please refer to the following document
Getting started
The current version of GWSkyNet uses the following libaries:
- numpy 1.20.3
- astropy 4.3.1
- reproject 0.7.1
- ligo.skymap 0.5.3
- tensorflow 2.9.1
- scipy 1.8.0
Using GWSkyNet is easy.
To apply GWSkyNet to a graviational wave candidate to determine the origin of the gravitational wave candidate (glitch or astrophysical) , all is needed is a FITS file corresponding to the gravitational wave alert. Such a FITS file may be downloaded from GraceDB
In the following example, we will use the FITS files corresponding to the Open Public Alert S190421ar.
To apply GWSkyNet to S190421ar, perform the following:
import GWSkyNet.GWSkyNet as G
model = G.load_GWSkyNet_model()
data = G.prepare_data('S190421ar.fits')
class_score = G.predict(model, data)
FAR, FNR = G.get_rates(class_score[0])
Running the above code will output a class score equal to 0.995. A class_score closer to 1 indicates the origin is more likely to be astrophysical. FAR and FNR stand for False Alarm Rate and False Negative Rate respectively.
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.