Glitch average task
This task seeks to answer the question: is there a high probability that data quality issues are present based on the long-term behavior of the detectors?
glitchaverage reports a \(p\)-value for the DQR based on the rate of Omicron triggers during the GW candidate time and comparing that to the rate of Omicron triggers seen in the past week of observing time. Should the rate of triggers be exceptional compared to the last week of data so that the \(p\)-value is \(0.05\) or less, this task raises alerts event validators.
Requirements
No additional dependencies are required. To view the help output python 3.9 is required.
Description
This task queries the last week of Omicron triggers to determine if the rate of triggers above a given signal to noise ratio near the GW candidate is exceptional. The Omicron paper describes in detail how these triggers are generated. Segments of observing time are broken down into smaller segments and an average glitch rate is calculated on each. The distribution of values of this average compared to the average rate of triggers during the GW time is presented as the main deliverable for the user. However, a bar chart of the background distribution with the rate of Omicron triggers in the foreground is given.
In addition, it may be useful to compare the rate of low-, medium- and high-frequency triggers for the past week of data versus the GW candidate time. The moving average of these triggers as a function of time is plotted along with the rate of low-, medium- and high- frequency triggers at the candidate time, as well.
Example command-line
This is the help message:
usage: dqr-glitch-average [-h] [-v] [-o OUTPUT_DIR] [--channel CHANNEL] [--analysis-ready-flag ANALYSIS_READY_FLAG] [--snr SNR] [--stride STRIDE] [--avg-len AVG_LEN] --gps GPS [--timespan TIMESPAN]
Main script function, glitchaverage
optional arguments:
-h, --help show this help message and exit
-v, --verbose increase verbose output
-o OUTPUT_DIR, --output-dir OUTPUT_DIR Directory for all output
--channel CHANNEL
--analysis-ready-flag ANALYSIS_READY_FLAG
--snr SNR SNR range. Default is 6.
--stride STRIDE Default is 60.
--avg-len AVG_LEN Default is 30.
--gps GPS
--timespan TIMESPAN