Skip to content
Snippets Groups Projects
Commit 24d1ad41 authored by Leo Tsukada's avatar Leo Tsukada
Browse files

gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs : exclude...

gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs : exclude jacobina determinat from norm calculation
parent e6630474
No related branches found
No related tags found
1 merge request!457gstlal-inspiral/bin/gstlal_inspiral_add_dt_dphi_snr_ratio_pdfs : exclude...
Pipeline #526859 passed
......@@ -55,7 +55,8 @@ for i in range(len(time["H1"])):
p2 = extract_dict(p, ifos)
s2 = extract_dict(s, ifos)
h2 = extract_dict(h, ifos)
norm[frozenset(ifos)] += TPS(t2,p2,s2,h2) * (sum(x**2 for x in s2.values())**.5)**4 / len(time["H1"]) * 1. / (16. * numpy.pi**2)
snr_net = sum(x**2 for x in s2.values())**.5
norm[frozenset(ifos)] += TPS(t2,p2,s2,h2) * snr_net**4 / len(time["H1"]) * 1. / (16. * numpy.pi**2) / (snr_net / numpy.prod([snr for snr in s2.values()]))
TPS.norm = norm
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment