Skip to content

AssertionError in gstlal_compute_far_from_snr_chisq_histograms

This is an error coming from the offline tutorial test. I forget, what did you think a solution would be?

gstlal_compute_far_from_snr_chisq_histograms \
	--non-injection-db H1L1-ALL_LLOID-1000000000-2048.sqlite \
	--background-bins-file marginalized_likelihood.xml.gz \
	--tmp-space /tmp \
	--verbose
reading 'marginalized_likelihood.xml.gz' ...
beginning count of above-threshold events
1/1: H1L1-ALL_LLOID-1000000000-2048.sqlite
/usr/lib64/python2.7/site-packages/gstlal/far.py:747: RuntimeWarning: invalid value encountered in divide
  survival_probability = zl_counts / noise_counts
/usr/lib64/python2.7/site-packages/gstlal/far.py:878: RuntimeWarning: invalid value encountered in divide
  ccdf /= ccdf[0]
Traceback (most recent call last):
  File "/usr/bin/gstlal_compute_far_from_snr_chisq_histograms", line 164, in <module>
    fapfar = far.FAPFAR(rankingstatpdf.new_with_extinction())
  File "/usr/lib64/python2.7/site-packages/gstlal/far.py", line 887, in __init__
    ccdf = gstlalstats.poisson_p_not_0(ccdf)
  File "/usr/lib64/python2.7/site-packages/gstlal/stats/__init__.py", line 117, in g
    p = f(*args, **kwargs)
  File "/usr/lib64/python2.7/site-packages/numpy/lib/function_base.py", line 1873, in __call__
    return self._vectorize_call(func=func, args=vargs)
  File "/usr/lib64/python2.7/site-packages/numpy/lib/function_base.py", line 1933, in _vectorize_call
    ufunc, otypes = self._get_ufunc_and_otypes(func=func, args=args)
  File "/usr/lib64/python2.7/site-packages/numpy/lib/function_base.py", line 1897, in _get_ufunc_and_otypes
    outputs = func(*inputs)
  File "/usr/lib64/python2.7/site-packages/gstlal/stats/__init__.py", line 156, in poisson_p_not_0
    assert l >= 0.
AssertionError