Skip to content
Snippets Groups Projects

O4a online dtdphi plotter

Merged Leo Tsukada requested to merge o4a-online-dtdphi-plotter into o4a-online
1 file
+ 2
1
Compare changes
  • Side-by-side
  • Inline
@@ -63,8 +63,9 @@ def dtdphi2prob(dt, dphi, ifo1, ifo2, refsnr, refhorizon):
phi = {ifo1: 0, ifo2: dphi}
snr = {ifo1: refsnr[ifo1], ifo2: refsnr[ifo2]}
horizon = {ifo1: refhorizon[ifo1], ifo2: refhorizon[ifo2]}
snr_net = sum(x**2 for x in snr.values())**.5
# signature is (time, phase, snr, horizon)
p = TPDPDF.time_phase_snr(t, phi, snr, horizon) * (sum(s**2 for s in snr.values())**.5)**4
p = TPDPDF.time_phase_snr(t, phi, snr, horizon) * (snr_net)**4 / (snr_net / numpy.prod([x for x in snr.values()]))
return float(p)
Loading