clear all;

% define the injection file
filename = 'adi_D.dat';

% load aLIGO strain noise
Sh = load('ALIGO-HighP_PSD.txt');

% calculate matched filter SNR
rho = mf_snr(filename, Sh, 1e-5);

