All-sky search for short gravitational-wave bursts in the first Advanced LIGO run
We present the results from an all-sky search for short-duration gravitational waves in the data of the first run of the Advanced LIGO detectors between September 2015 and January 2016. The search algorithms use minimal assumptions on the signal morphology, so they are sensitive to a wide range of sources emitting gravitational waves. The analyses target transient signals with duration ranging from milliseconds to seconds over the frequency band of 32 to 4096 Hz. The first observed gravitational-wave event, GW150914, has been detected with high confidence in this search; other known gravitational-wave events fall below the search’s sensitivity. Besides GW150914, all of the search results are consistent with the expected rate of accidental noise coincidences. Finally, we estimate rate-density limits for a broad range of non-BBH transient gravitational-wave sources as a function of their gravitational radiation emission energy and their characteristic frequency. These rate-density upper-limits are stricter than those previously published by an order-of-magnitude.
The first observing period of the Advanced LIGO detectors ;  has been completed recently with the most sensitive gravitational-wave (GW) detectors ever built. The two LIGO observatories in Hanford, WA and Livingston, LA achieved a major milestone in gravitational wave astronomy: the first direct detection of gravitational waves on September 14, 2015, referred as GW150914 . Advanced LIGO is the first of a new generation of instruments, including GEO 600 , Advanced Virgo , KAGRA  and LIGO-India .
This paper reports on a search for short-duration transient gravitational-wave events, commonly referred to as GW bursts, during the first observing run (O1) of the Advanced LIGO detectors, from September 2015 to January 2016. The first 16 days of coincident data have been already analyzed, resulting in a high-significance detection statement for the GW150914 event . GW bursts can be generated by a wide variety of astrophysical sources, such as merging compact binary systems ; , core-collapse supernovae of massive stars , neutron stars collapsing to form black holes, pulsar glitches, and cosmic string cusps . To search broadly for these phenomena, we employ searches with minimal assumptions regarding the expected waveform characteristics and the source direction. The search we report here is more sensitive than the previous burst searches  because of both the increased sensitivity of the Advanced detectors  and improvements in the search algorithms in rejecting transient non-Gaussian noise artifacts (glitches) ; ; ; .
The described un-modeled all-sky search for GW bursts consists of three different algorithms. This paper shows the result of these algorithms, and gives limits on the rate-density of transient GW events. All of these algorithms have independently claimed high-significance detections of GW150914 . The lower-mass GW event, GW151226 , and the LVT151012 candidate ;  were not detected by these searches.
Ii Observing Run 1
Our data set extends over 130 calendar days from September 12, 2015 to January 19, 2016. This first observing period (called O1) of Advanced LIGO began after a series of major upgrades to both the Hanford and Livingston detectors .
In the most sensitive frequency band, 100-300 Hz, the O1 LIGO detectors are 3 to 5 times more sensitive than the initial LIGO detectors . Future observing runs are expected to increase sensitivity by an additional factor of 3 .
As in the previous LIGO/Virgo searches ; ; , intervals of poor data quality are identified and excluded from the analysis. To monitor environmental disturbances and their influence on the detectors, each observatory is equipped with an array of sensors: seismometers, accelerometers, microphones, magnetometers, radio receivers, weather sensors, ac-power line monitors, and a cosmic-ray detector. Hundreds of thousands of auxiliary channels within the instrument are also monitored. Characterization of the relationship of the strain data to this additional information allows many non-GW transients to be removed with high statistical confidence ; .
The livetime in which the two detectors were individually locked is about 79 days for H1 and 67 days for L1. After data quality flags have been applied, the total analyzable time is about 75 days for H1 and 65 days for L1. The coincident livetime between H1 and L1 is about 48 days. This livetime includes the 16 days of this coincident data that has already been analyzed in . Finally, the estimated calibration uncertainty () below 2 kHz is less than 10% in amplitude and 10 degrees in phase . The calibration uncertainty above 2 kHz is less certain, although the limited data obtained at these frequencies suggests upper bounds of 20% in amplitude and 10 degrees in phase. These estimates will be further refined through future measurements and analyses .
This search covers the most sensitive frequency band of the involved detectors, i.e. 32 - 4096 Hz, and it consists of the same three burst algorithms used to measure the significance of GW150914 . They consist of two end-to-end algorithms: coherent Waveburst (cWB) ;  and omicron-LIB (oLIB) ; and a follow-up algorithm applied to cWB events: BayesWave (BW) ; . Using multiple search algorithms has two advantages: it can provide independent validation of results, and it can also improve the search sensitivity in regions of parameter space where a single algorithm outperforms the others.
The three algorithms ran over the 48 days of coincident data. However, due to internal segmentation
The three algorithms also ran in low-latency mode during O1. In this mode, both cWB and oLIB produced independent alerts of the GW150914 event and the result was validated by a BW follow-up .
To characterize the statistical rate of transient noise glitches occurring simultaneously at the two LIGO sites by chance, this analysis uses the time-shift method: data from one interferometer is shifted in time with respect to the other interferometer by multiple delays much larger than the maximum GW travel time between the interferometers. In this way, we can accumulate a significant duration of estimated background that we use to estimate the false-alarm rate (FAR) for each algorithm.
We set a FAR threshold of 1 in 100 years for identifying a detection candidate, which roughly corresponds to a 3 sigma detection statement for the duration of our observation. If an event in this search were to have a FAR less than this threshold, a refined analysis (i.e., more time-shifts) would be performed to assign the appropriate significance in the detection statement for this event.
iii.1 Coherent WaveBurst
Coherent WaveBurst (cWB) has been used in multiple searches for transient GWs ; . It calculates a maximum-likelihood-ratio statistic for power excesses identified in the time-frequency domain. A primary selection cut is applied to the network correlation coefficient , which measures the degree of correlation between the detectors. Events with are discarded from the analysis. Events are ranked according to their coherent network signal-to-noise ratio (SNR) , which is related to the matched-filter SNR, favoring GW signals correlated in both detectors and suppressing uncorrelated glitches. A detailed explanation of the algorithm and the definition of these statistics are given in .
The cWB analysis is divided in two frequency bands, where the splitting frequency is 1024 Hz. For the low-frequency band, the data is downsampled to reduce the computational cost of the analysis.
Low-frequency cWB events are divided into three search classes according to their morphology, as described in . The class is based on cuts which primarily select so-called “blip” glitches and non-stationary power-spectrum lines. The former are non-Gaussian noise transients of unknown origin consisting of a few cycles around 100 Hz. The class is based on cuts that select events whose frequency increases with time, i.e. those similar in morphology to the merger of compact objects. The class is composed of all the remaining events.
The FAR of each identified event is estimated using the time-slide background distribution of similar class. Since there are three independent classes, we apply a trials factor of 3 to estimate the final significance. The high-frequency analysis consists of only a single class.
Omicron-LIB (oLIB) is a hierarchical search algorithm that first analyzes the data streams of individual detectors, which we refer to as an incoherent analysis. It then follows up stretches of data that are potentially correlated across the detector network, which we refer to as a coherent analysis. The incoherent analysis (“Omicron”)  flags stretches of coincident excess power. The coherent follow-up (“LIB”)  models gravitational wave signals and noise transients with a single sine-Gaussian, and it produces two different Bayes factors. Each of these Bayes factors is expressed as the natural logarithm of the evidence ratio of two hypotheses: a GW signal vs Gaussian noise (BSN) and a coherent GW signal vs incoherent noise transients (BCI). The joint likelihood ratio of these two Bayes factors is used as ranking statistic to assign a significance to each event. See  for further technical details on the implementation of these steps.
For this analysis, oLIB events are divided into two classes, based on the inferred parameters of the best-fit sine-Gaussian. The exact parameter ranges of these search classes are chosen in order to group noise transients of similar morphology together. Particularly noisy regions of the parameter space are excluded from the analysis entirely (e.g., events with median quality factor > 108). Both classes contain only events whose median frequency , as estimated by LIB, lies within the range of 48 - 1024 Hz. The first, analogous to cWB’s C1 class, is a “low-Q” class that contains only events whose median quality factor , lies within the range 0.1 - 2. The second, analogous to the union of cWB’s C2 and C3 classes, is a “high-Q” class that contains only events whose median lies within the range 2 - 108. In both classes, event candidates were also required to have positive Bayes factors, i.e., BSN > 0 and BCI > 0, meaning the evidence for the signal model was greater than the evidences for the noise models. A trials factor of 2 accounts for these independent search classes.
The oLIB background analysis is performed using 456 years of background data. We select single-detector events with SNR > 5.0. This is lower than the threshold of 6.5 adopted in  and it is chosen to allow us to make a significance estimation of low-SNR events. For this reason, we cannot directly compare the two set of results reported in  and in this study using the likelihood ratios , but we have to consider the reported FAR. The results are presented in Fig. c.
iii.3 BayesWave Follow-up
BayesWave (BW) tests if the data in multiple detectors are best explained by coincident glitches or a signal, and it is used as a follow-up to events produced by cWB. It has been shown that BW is able to increase the detection confidence for GW signals of complex morphology.
The BW algorithm uses a variable number of sine-Gaussian wavelets to reconstruct the data independently for the signal and glitch models, then computes the natural logarithm of the Bayes factor between these two models, . The number of wavelets used is determined by using a reversible jump Markov chain Monte Carlo, with more complex signals requiring more wavelets . The Bayes factor scales as , where is number of wavelets used. This means the detection statistic depends on waveform complexity in addition to SNR. Full details of the algorithm can be found in .
In this search, BW followed up events produced by cWB in any of the three low-frequency search classes with a coherent network SNR of and correlation coefficient of . There are no additional cuts performed on the data, and all of these events (C1+C2+C3) are analyzed as a single class. The cumulative FAR as a function of is shown in Fig. d.
The detection efficiency of the search is measured by adding simulated signals into the detectors’ data and evaluating whether or not they pass the selection cuts explained in the Section III for the different search algorithms. A variety of GW signal morphologies were tested, spanning a wide range of amplitudes and duration, and with characteristic frequencies within the sensitive bandwidth of the detectors. We identify two different waveform sets: a set of generic bursts, and a set of simulated astrophysical signals coming from the coalescence and merging of binary black holes (BBH). All of the results in this section refer to a FAR detection threshold of 1 in 100 years.
iv.1 Generic bursts
This family includes the waveform types described in , all with elliptical polarization: gaussian pulses (GA), parametrized by their duration parameter ; sine-Gaussian wavelets (SG), sinusoids within a Gaussian envelope, characterized by the frequency of the sinusoid and a quality factor ; white-noise bursts (WNB), white noise bounded in frequency over a bandwidth and with a Gaussian envelope, described by the lower frequency , , and the duration . Table 1 lists the waveforms that have been considered for this work.
The amplitudes of the test signals are chosen to cover a wide range of values and are expressed in terms of the root-mean-square strain amplitude at Earth (before accounting for the detection response patterns), denoted . For this search, we injected signals according to the distance distribution where is a constant. The constant is chosen to produce at least several test events with large .
|Hz, Hz, s||2.0||N/A||3.0|
|Hz, Hz, s||2.2||N/A||9.2|
Table 1, shows the value at which 50% of the injections are detected for each signal morphology and algorithm. There are some morphology-dependent features that affect each of the different algorithms at the FAR threshold of 1 in 100 years. These features largely disappear and the different algorithms’ results converge at detection thesholds of higher FAR. For example, the detection efficiencies are worse for cWB for low-Q morphologies and high-Q morphologies because these injections are classified as C1 events. As shown in Fig. a, the C1 background extends to higher significances than in the other bins, meaning these high-Q and low-Q events must have large values of to meet the FAR threshold of 1 in 100 years. The oLIB detection efficiencies, while non-negligible across all morphologies, never quite reach 50% for some non-sine-Gaussian morphologies because the template mismatch residuals grow linearly with . Finally, the detection efficiencies of BW suffers for high-Q events since its prior range only extends to . However, almost every morphology can be detected efficiently by at least one of the algorithms.
Another way to interpret the search sensitivities is to map them into the minimum amount of energy that needs to be emitted through GWs for at least half of the sources to be detected within a given search volume. Assuming a fixed amount of energy is radiated isotropically away from the source in GWs of a fixed frequency , this distance can be converted into a value of via the relationship :
Here, we use the from Table 1, the central frequency of each morphology, and a fixed fiducial radius to calculate this energy via Eq. 1. Figure 2 shows this energy as a function of characteristic frequency assuming a galactic source at a distance of 10 kpc. When taking into account the results of all three algorithms, this emission energy is not strongly dependent on the type of waveform (with exceptions on an algorithm-by-algorithm basis, as described above). Fig. 2 can easily be converted to other distances by applying the scaling relation suggested by Eq. 1. Previous studies  have published similar emission-energy-versus-frequency plots at a detection threshold of 1 in 8 years. We note that the current results, when evaluated at this higher-FAR threshold, are roughly an order-of-magnitude more sensitive than these previous results, due mainly to the improvement in detector sensitivites.
iv.2 Binary black holes mergers
We also consider a set of astrophysical waveforms using models of merging of binary black hole systems.
Specifically, we choose the SEOBNRv2 model as implemented in the LAL software library ; .
The waveforms are generated with an initial frequency of Hz.
The simulated binary systems are isotropically located in the sky and isotropically oriented.
The total redshifted mass of the system in the detector frame
In Fig. 3 we compare the sensitive luminosity radius  as a function of the total redshifted mass in the detector frame. While systems inside this distance may be missed and systems outside of it may be detected depending on their sky position and orientation, this sensitive radius provides a “rule-of-thumb” determination on whether or not this burst search will detect a system’s GW transients. We can see that for systems like GW150914 ( ) and GW151226 ( ), the search ranges at the FAR of 1/100 years are approximately 500-700 Mpc and 100-200 Mpc, respectively. These ranges demonstrates why this search detects GW150914 ( Mpc ) but not GW151226 ( Mpc ). Even though the two sources are at a similar luminosity distance, this burst search is less efficient at detecting low-mass BBH systems. This behavior is true for two reasons: lower-mass systems emit less energy into GWs than higher-mass systems, and this energy is distributed over a longer duration of time. These two features make it more difficult for non-templated algorithms to extract the GW signal from the detector noise as compared to searches based on templates.
The most significant event and only detection established in this search is GW150914 ,
which is independently confirmed by all three algorithms. Specifically, it is found by cWB in the C3 class of
the low-frequency analysis with an estimated FAR of less than 1 in 350 years, by oLIB in the
“high-Q” class with an estimated FAR of less than 1 in 230 years, and by BayesWave with an estimated FAR
of less than 1 in 1000 years.
All other events generated by the analyses are consistent with the accidental noise coincidence rates. To be specific, there are no other events found above the SNR thresholds in either the “low-Q” class of oLIB or the entire BayesWave analysis bin. The rate of other events in the oLIB “high-Q” bin are consistent with the accidental noise coincidence rates within 1 sigma. The event in the cWB analysis with the second-lowest FAR belongs to the high frequency search, with a false-alarm probability of about 0.2.
These results set constraints on the population of transient GW sources within the volume of the Universe that the detectors were sensitive to during O1. Again, all of the results in this section refer to a FAR detection threshold of 1 in 100 years.
We estimate the limits on the rate density of generic non-BBH-like GW-burst sources in Fig. 4 by removing the known BBH detections GW150914 and GW151226 from our analysis. We use the sine-Gaussian injection set as a representative morphology, and present our cWB rate-density estimates as a function of their characteristic frequencies. The bands represent the 90% confidence intervals on rate density , calculated using the Feldman-Cousins formalism for 0 background events . The frequency-dependent variation among the upper-limits is due to the sine-Gaussians falling into different cWB search classes as a result of their specific value of . For a given value of , the results follow a smoother frequency dependence. These results are not directly comparable with those from previous runs  because of the different FAR detection thresholds. However, we note that at the previously-used FAR detection threshold of 1 in 8 years, our search lowers these upper-limits by about an order of magnitude across all frequencies. The sensitivity improvements of the detectors and pipelines allow us to make these stricter rate statements even though we analyzed less livetime compared to  (less than 50 days compared to 1.7 yr). Fig. 4 assumes of gravitational wave energy has been emitted from the source, but this can be scaled to any emission energy by using Eq. 1. Note that the rate density scales as .
This paper reports the results for the search for short duration GW in the first Advanced LIGO observing run, with minimal assumptions on the signal waveform, direction or arrival time. The two LIGO detectors, Livingston and Hanford, were operating from mid-September 2015 to mid-January 2016, with a greater sensitivity to GWs than any previous LIGO-Virgo run. This search has been performed considering two end-to-end algorithms and a follow-up algorithm.
The only detection established in this search is the GW150914 event, a binary system consisting of two black holes merging to form a single one . The other known black hole detection  falls below the sensitivity of this search, and all other events in the search result are consistent with accidental noise coincidences between the detectors.
We report the minimum GW emission energy needed to detect at least half of the transient events emitted within some fiducial distance. These energies depend primarily on the signal frequency and are approximately constant over the different models of GW emission morphology. We also estimate rate-density limits on non-BBH transient sources as a function of their frequency and their gravitational wave emission energy.
The interferometric detectors LIGO and Virgo are currently being upgraded for the next scientific run. LIGO should improve its sensitivity over the next few years, Virgo should soon come online, and the implementation of KAGRA and LIGO India is also in progress. All of these improvements will allow this type of un-modeled search to achieve a better sensitivity in the future .
Acknowledgments — The authors gratefully acknowledge the support of the United States National Science Foundation (NSF) for the construction and operation of the LIGO Laboratory and Advanced LIGO as well as the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Foundation for Fundamental Research on Matter supported by the Netherlands Organisation for Scientific Research, for the construction and operation of the Virgo detector and the creation and support of the EGO consortium. The authors also gratefully acknowledge research support from these agencies as well as by the Council of Scientific and Industrial Research of India, Department of Science and Technology, India, Science & Engineering Research Board (SERB), India, Ministry of Human Resource Development, India, the Spanish Ministerio de Economía y Competitividad, the Conselleria d’Economia i Competitivitat and Conselleria d’Educació, Cultura i Universitats of the Govern de les Illes Balears, the National Science Centre of Poland, the European Commission, the Royal Society, the Scottish Funding Council, the Scottish Universities Physics Alliance, the Hungarian Scientific Research Fund (OTKA), the Lyon Institute of Origins (LIO), the National Research Foundation of Korea, Industry Canada and the Province of Ontario through the Ministry of Economic Development and Innovation, the Natural Science and Engineering Research Council Canada, Canadian Institute for Advanced Research, the Brazilian Ministry of Science, Technology, and Innovation, Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Russian Foundation for Basic Research, the Leverhulme Trust, the Research Corporation, Ministry of Science and Technology (MOST), Taiwan and the Kavli Foundation. The authors gratefully acknowledge the support of the NSF, STFC, MPS, INFN, CNRS and the State of Niedersachsen/Germany for provision of computational resources.
This article has been assigned the document number P1600129.
- The cWB algorithm requires at least 600 s of continuous data to perform its analysis.
- Given the luminosity distance of the system, one can assume a cosmology and calculate its redshift . The system’s total mass in the source frame can then be obtained by dividing the total redshifted mass in the detector frame by .
- Because GW150914 was louder than any of the background events in this search, we can only provide the relatively un-precise upper-limits on FAR listed above.
- J. Aasi et al. Advanced LIGO. Class. Quant. Grav., 32:074001, 2015.
- F. Acernese et al. Advanced Virgo: a second-generation interferometric gravitational wave detector. Class. Quant. Grav., 32(2):024001, 2015.
- B. Abbott et al. Observation of Gravitational Waves from a Binary Black Hole Merger. Phys.Rev.Lett., 116:061102, 2016. LIGO-P150914.
- H. Lück et al. The upgrade of GEO600. J. Phys. Conf. Ser., 228:012012, 2010.
- Y. Aso, Y. Michimura, K. Somiya, M. Ando, O. Miyakawa, T. Sekiguchi, D. Tatsumi, and H. Yamamoto. Interferometer design of the KAGRA gravitational wave detector. Phys. Rev., D88(4):043007, 2013.
- J. Aasi et al. Prospects for Observing and Localizing Gravitational-Wave Transients with Advanced LIGO and Advanced Virgo. Living Rev. Relativity, 19:1, 2016.
- B. Abbott et al. Observing gravitational-wave transient GW150914 with minimal assumptions. Phys.Rev. D, 93:122004, 2016.
- J. Aasi et al. Search for gravitational radiation from intermediate mass black hole binaries in data from the second LIGO-Virgo joint science run. Phys.Rev., D89:122003, 2014.
- S. Mohapatra, L. Cadonati, S. Caudill, J. Clark, C. Hanna, S. Klimenko, C. Pankow, R. Vaulin, G. Vedovato, and S. Vitale. Sensitivity Comparison of Searches for Binary Black Hole Coalescences with Ground-based Gravitational-Wave Detectors. Phys. Rev., D90(2):022001, 2014.
- C. L. Fryer and K. C. B. New. Gravitational waves from gravitational collapse. Living Rev. Rel., 14:1, 2011.
- Thibault Damour and Alexander Vilenkin. Gravitational radiation from cosmic (super)strings: Bursts, stochastic background, and observational windows. Phys. Rev. D, 71:063510, Mar 2005.
- J. Abadie, B. P. Abbott, R. Abbott, T. D. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Affeldt, and et al. All-sky search for gravitational-wave bursts in the second joint LIGO-Virgo run. Phys. Rev. D, 85(12):122007, June 2012.
- J. Aasi et al. Characterization of the LIGO detectors during their sixth science run. Class. Quant. Grav., 32(11):115012, 2015.
- J.B. Kanner, T.B. Littenberg, N. Cornish, M. Millhouse, E. Xhakaj, F. Salemi, M. Drago, G. Vedovato, and S. Klimenko. Leveraging waveform complexity for confident detection of gravitational waves. Phys. Rev., D93(2):022002, 2016.
- T.B. Littenberg, J.B. Kanner, N.J. Cornish, and M. Millhouse. Enabling high confidence detections of gravitational-wave bursts. Phys. Rev., D94:044050, 2016.
- R. Lynch, S. Vitale, R. Essick, E. Katsavounidis, and F. Robinet. An information-theoretic approach to the gravitational-wave burst detection problem. arXiv:1511.05955, 2015.
- S. Klimenko, G. Vedovato, M. Drago, F. Salemi, V. Tiwari, G. A. Prodi, C. Lazzaro, S. Tiwari, F. Da Silva, and G. Mitselmakher. Method for detection and reconstruction of gravitational wave transients with networks of advanced detectors. Phys. Rev. D 93, 042004, 2015.
- B. Abbott et al. GW151226: Observation of Gravitational Waves from a 22-Solar-Mass Binary Black Hole Coalescence . Phys.Rev.Lett., 116:241103, 2016. LIGO-P151226.
- B. Abbott et al. GW150914: First results from the search for binary black hole coalescence with Advanced LIGO. arXiv:1602.03839, 2016.
- B. Abbott et al. Binary Black Hole Mergers in the first Advanced LIGO Observing Run. arXiv:1606.04856, 2016.
- B.P. Abbott et al. An all-sky search for long-duration gravitational wave transients with LIGO. Phys.Rev. D, 93:042005, 2015.
- J. Abadie et al. All-sky search for gravitational-wave bursts in the second joint LIGO-Virgo run. Phys. Rev., D85:122007, 2012.
- J. Abadie et al. All-sky search for gravitational-wave bursts in the first joint LIGO-GEO-Virgo run. Phys. Rev., D81:102001, 2010.
- B. Abbott et al. Characterization of transient noise in Advanced LIGO relevant to gravitational wave signal GW150914. Class. Quant. Grav., 33(13), 2016.
- B. Abbott et al. Calibration of the Advanced LIGO detectors for the discovery of the binary black-hole merger GW150914. arXiv:1602.03845, 2016.
- C. Callihane and J.Kissel. Private communication.
- S. Klimenko, I. Yakushin, A. Mercer, and Guenakh Mitselmakher. Coherent method for detection of gravitational wave bursts. Class. Quant. Grav., 25:114029, 2008.
- N.J. Cornish and T.B. Littenberg. BayesWave: Bayesian Inference for Gravitational Wave Bursts and Instrument Glitches. Class. Quant. Grav., 32(13):135012, 2015.
- T.B. Littenberg and N.J. Cornish. Bayesian inference for spectral estimation of gravitational wave detector noise. Phys. Rev., D91(8):084034, 2015.
- B. Abbott et al. Localization and broadband follow-up of the gravitational-wave transient GW150914. The Astrophysical Journal Letters, 826:L13, 2016.
- F. Robinet. Omicron: an algorithm to detect and characterize transient noise in gravitational-wave detectors. https://tds.ego-gw.it/ql/?c=10651, 2015.
- A. Taracchini et al. Effective-one-body model for black-hole binaries with generic mass ratios and spins. Phys. Rev., D89(6):061502, 2014.
- P. Kumar, K. Barkett, S. Bhagwat, N. Afshari, D.A. Brown, G. Lovelace, M.A. Scheel, and B. Szilágyi. Accuracy and precision of gravitational-wave models of inspiraling neutron star-black hole binaries with spin: Comparison with matter-free numerical relativity in the low-frequency regime. Phys. Rev., D92(10):102001, 2015.
- the LIGO Scientific Collaboration, the Virgo Collaboration: J. Abadie, B. P. Abbott, R. Abbott, T. D. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, and et al. Search for Gravitational Waves from Intermediate Mass Binary Black Holes. Phys. Rev. D, D85:102004, January 2012.
- B. Abbott et al. Properties of the Binary Black Hole Merger GW150914. Phys. Rev. Lett., 116(24):241102, 2016.
- G. J. Feldman and R. D. Cousins. A Unified Approach to the Classical Statistical Analysis of Small Signals. Phys.Rev., D57:3873, 1998.