add transform to compute the matched-filter (network) SNR based on the true/optimal SNR
Just have this draw random shifts from the appropriate distributions (some combination of non-central chi and/or non-central chi-squared distributions). Then each drawn injection will essentially represent a single noise realization. We can go from there...
With this in place, we can then threshold on the "detected network SNR" rather than the "optimal network SNR"