The following example program shows one way of combining the functions detector_site(), noise_power(), overlap(), and optimal_filter() to calculate the spectrum of the optimal filter function for a given pair of detectors. Below we explictly calculate for the initial Hanford, WA and Livingston, LA LIGO detectors. (We also choose to normalize the magnitude of the spectrum to 1, for later convenience when making plots of the output data.) Noise power information for these two detectors is read from the input data file noise_init.dat. This file is specified by the information contained in detectors.dat. (See Sec. for more details.) The resulting optimal filter function data is stored as two columns of double precision numbers ( and ) in the file LIGO_filter.dat, where and . A plot of this data is shown in Fig. .
As usual, the user can modify the parameters in the #define statements listed at the beginning of the program to change the number of frequency points, the frequency spacing, etc. used when calculating . Also, by changing the site location identification numbers and the output file name, the user can calculate and save the spectrum of the optimal filter function for any pair of detectors. For example, Fig. is a plot of the optimal filter function for the advanced LIGO detectors.
Includes/optimal_filter.tex