search

gwastro / pycbc

Core package to analyze gravitational-wave data, find signals, and study their parameters. This package was used in the first direct detection of gravitational waves (GW150914), and is used in the ongoing analysis of LIGO/Virgo data.
http://pycbc.org
GNU General Public License v3.0
312 stars 347 forks source link

Upload PSDs to GraceDB #550

Closed lpsinger closed 8 years ago

lpsinger commented 8 years ago

Some automatic analyses of CBC triggers that are uploaded to GraceDB, including rapid sky localization with BAYESTAR, depend on noise PSDs from the time of the candidate. In order to enable rapid sky localization for pycbc events, pycbc must upload a PSD to GraceDB along with each coinc.

The best place to do this is probably in pycbc_upload_xml_to_gracedb. The PSD should be in the LIGO-LW XML format that is produced by lal.series.make_psd_xmldoc. The file should be named psd.xml.gz and should be uploaded with the psd tag.

Note that the disk usage of one PSD file per GraceDB event is not a major consideration, because the BAYESTAR sky maps files typically have about the same size as the PSD files.

lpsinger commented 8 years ago

I take back what I originally said about the PSD file size: gstlal and MBTA produce PSD files with sizes of ~100 KB, whereas this patch is generating PSD files of ~10 MB. However, from looking at a handful of pycbc runs, it looks like there are about as many entries in the loudest events table as there are PSD segments. So, my conclusion is the same: storing the PSDs separately would produce no significant savings in disk space in GraceDB.

spxiwh commented 8 years ago

Great! Thanks a lot for doing this.