OutlierMoni -- Quick reference (draft!)
OutlierMoni -- Quick reference
OutlierMoni
Summary
OutlierMoni monitors short transients of large amplitude in a set of channels. OutlierMoni computes a whitened version of each input channel and it looks for outliers (i.e., points larger than a given threshold) in the resulting stream.
The output of OutlierMoni can be accessed with dataDisplay by opening raw.ffl. It is also available by a web interface.
OutlierMoni's home page: http://wwwcascina.virgo.infn.it/MonitoringWeb/Noise/outlierMoni
Principle
Whiteningwhitening = make the spectrum flat. in the frequency domain, this amounts to dividing by 1/sqrt(PSD).
since the PSD is unknown a priori, PSD has to be estimated : using running Welch. averaging time can be tuned by the
computation and output frame wise (one per GPS second)
output homogeneous to a SNR.
Detailed description
Several parameters are required OutlierMoniThe OutlierMoni's configuration is defined by /virgoData/Noise/outlierMoni.cfg
We detail a typical configuration.
Vetoing
a condition can be put on OutlierMoni execution. for instance, OutlierMoni can be enabled only in Science Mode.
/* keyword veto_channel threshold */Whitening
#VETO_CHANNEL Qc_Moni..ScienceMode 1 1
since the PSD is unknown a priori, PSD has to be estimated : using running Welch. averaging time can be tuned by the
parameter "Average" (see Configuration below).
this is an adaptive procedure: follows changes in the PSD with typical time longer than averaging time.
/* nsamples for FFT nAverage min/max frequency (percent of Nyquist frequency) Nb Nt */
/* keyword Nf a_max fmin fmax Nb Nt */
OUTLIER 70000 10 0. 0.5 0 0
when OutlierMoni enabled, whitening stream not available immediately: require "warming cycles" (number depends on averaging time).
/* choose single whitening */if so configured, the whitening stream is sent also to the output.
OUTLIER_SINGLE 1
/* channel for which whitened signal will be put in the data */Thresholding
OUTLIER_WHITENED_CHANNEL Pr_B1_ACp
/* keyword channel threshold */all computation of OutlierMoni is based on the frame.
MONITOR Pr_B1_ACp 6
compute the max in the current frame, threshold the max.
1 flag per frame.
Output
detection flag, SNR (calibrated)
Configuration
The monitor load its parameters from a configuration file.
Software
The source code of OutlierMoni can be found in /virgoDev/outlierMoni.Related presentations
http://wwwcascina.virgo.infn.it/collmeetings/presentations/May2005/DataAnalysis/ChassandeMottin_transients_in_C5-M1.pdfhttp://wwwcascina.virgo.infn.it/collmeetings/presentations/2005/2005-10/DataAnalysis/Chassande_outlier_C6_and_C7.pdf