Maintenance will be performed on git.ligo.org, chat.ligo.org, containers.ligo.org, and docs.ligo.org on Tuesday 7th July 2020 starting at approximately 10am PDT and lasting for around 15 minutes. There will be a short period of downtime towards the end of the maintenance window. Please direct any comments, questions, or concerns to uwm-help@cgca.uwm.edu.

Commit 46e3d41f authored by Chad Hanna's avatar Chad Hanna

inspiral_lr, inspiral_extrinsics: more snr chisq tuning

parent 39bed653
Pipeline #75582 passed with stages
in 29 minutes and 27 seconds
......@@ -945,7 +945,7 @@ class NumeratorSNRCHIPDF(rate.BinnedLnPDF):
snr2 = snr**2.
ncparam_per_pf = snr2
# takes into account the mean depending on noncentrality parameter
snrchi2 = numpy.outer(snr2 * df * (1.0 + numpy.mean(pfs)), rcoss)
snrchi2 = numpy.outer(snr2 * df * (1.0 + max(pfs)), rcoss)
arr = numpy.zeros_like(lnpdf.array)
for pf in pfs:
......
......@@ -423,7 +423,7 @@ class LnSignalDensity(LnLRDensity):
vtdict = self.horizon_history.functional_integral_dict(window.shift(float(gps)), lambda D: D**3.)
return dict((instrument, (vt / t)**(1./3.)) for instrument, vt in vtdict.items())
def add_signal_model(self, prefactors_range = (0.001, 0.30), df = 150, inv_snr_pow = 4.):
def add_signal_model(self, prefactors_range = (0.001, 0.30), df = 100, inv_snr_pow = 4.):
# normalize to 10 *mi*llion signals. this count makes the
# density estimation code choose a suitable kernel size
inspiral_extrinsics.NumeratorSNRCHIPDF.add_signal_model(self.densities["snr_chi"], 1e12, prefactors_range, df, inv_snr_pow = inv_snr_pow, snr_min = self.snr_min)
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment