function rho = mf_snr(filename, Sh)
% function rho = mf_snr(filename, Sh)
% Sh = N x 2 array of frequencies and noise power spectral density
% filename = points to a text file with columns of t, hp, hx

% load injection file
h = load(filename);

% parse injection data
t = h(:,1);
hp = h(:,2);
hx = h(:,3);

% get detector noise
sampleRate = 1 / (t(2)-t(1));
duration = t(end);
h = gaussian_noise(Sh, sampleRate, duration);



% deal with antenna factors
