B dL@sddlZddlZddlZddlZddlmZmZmZm Z ddl m Z m Z m Z mZddlmZddlZdZddZddd Zeed d Zeed d ZddZdadadadddZddZdddZGdddeZGdddeZ dS)N)zerosreal_same_precision_as TimeSeriescomplex_same_precision_as)sigmasq_seriesmake_frequency_seriesmatched_filter_coreget_cutoff_indices)schemedzpycbc.vetoes.chisq_cCsL||d}td|||}tj||||dd}||7}t||S)aReturns bins of equal power for use with the chisq functions Parameters ---------- sigmasq_series: FrequencySeries A frequency series containing the cumulative power of a filter template preweighted by a psd. num_bins: int The number of chisq bins to calculate. kmin: int DOCUMENTME kmax: int DOCUMENTME Returns ------- bins: List of ints A list of the edges of the chisq bins is returned. rright)Zside)numpyZarangeZ searchsortedappend)rnum_binskminkmaxZsigmasqZedge_vecbinsr_/work/yifan.wang/ringdown/master-ringdown-env/lib/python3.7/site-packages/pycbc/vetoes/chisq.py$power_chisq_bins_from_sigmasq_series!s  rcCs@t||||}t|||jt|dd\}}t||||S)aReturns bins of equal power for use with the chisq functions Parameters ---------- htilde: FrequencySeries A frequency series containing the template waveform num_bins: int The number of chisq bins to calculate. psd: FrequencySeries A frequency series containing the psd. Its length must be commensurate with the template waveform. low_frequency_cutoff: {None, float}, optional The low frequency cutoff to apply high_frequency_cutoff: {None, float}, optional The high frequency cutoff to apply Returns ------- bins: List of ints A list of the edges of the chisq bins is returned. r )rr r delta_flenr)htilderpsdlow_frequency_cutoffhigh_frequency_cutoffZ sigma_vecrrrrrpower_chisq_bins>s rcCsd}|d7}t|dS)NzcCsft|}t|}t|||||} tt|dd|jd} t|||||| d\} } } t| | | | |dS)aCalculate the chisq timeseries Parameters ---------- template: FrequencySeries or TimeSeries A time or frequency series that contains the filter template. data: FrequencySeries or TimeSeries A time or frequency series that contains the data to filter. The length must be commensurate with the template. (EXPLAINME - does this mean 'the same as' or something else?) num_bins: int The number of frequency bins used for chisq. The number of statistical degrees of freedom ('dof') is 2*num_bins-2. psd: FrequencySeries The psd of the data. low_frequency_cutoff: {None, float}, optional The low frequency cutoff for the filter high_frequency_cutoff: {None, float}, optional The high frequency cutoff for the filter return_bins: {boolean, False}, optional Return a list of the individual chisq bins Returns ------- chisq: TimeSeries TimeSeries containing the chisq values for all times. r r)r*)Zcorr_out)r6)rrrrr*rr8)templatedatarrrrr6rZstilderZcorraZ total_snrr%Ztnormrrr power_chisqs rAc@s6eZdZdZd ddZeddZdd Zd d ZdS) SingleDetPowerChisqzClass that handles precomputation and memory management for efficiently running the power chisq in a single detector inspiral analysis. rNcCs:|dks*|dks*d|_d|_d|_||_nd|_||_dS)N0rTrZ chisq_dofF)doZ column_nameZtable_dof_namer snr_threshold)selfrrErrr__init__6szSingleDetPowerChisq.__init__cCs@ttd}||j|tj|tjjt|ddi|S)N)maxmin __builtins__)rHrIupdate__dict__mathr4Zpnutilseval)rowargZ safe_dictrrr parse_option@s    z SingleDetPowerChisq.parse_optionc Csddlm}t|}t|ds$i|_t|ds:|dd|_||jksTt|j|jkrd|jt|j<t|||j }t|dr|j |j krt|j |j }|j}t|j |j |||}nt||||j }||j|<|j|S) Nr)LimitedSizeDict_chisq_cached_key _bin_cache)Z size_limitT sigmasq_vec)Z pycbc.optrRidhasattrrSrTparamsr3rQr approximantrVf_lowerrend_idxrr) rFr?rrRkeyrrrrrrrcached_chisq_binsHs*        z%SingleDetPowerChisq.cached_chisq_binscCs|jrt|}|jrft|||jk}|}td|||} ||} tjt|tj d} d} n|} |} |dkr| ||} t| ddd} t || || | }|jr|dkr|| |<n|} | t | t|fSdSdS) a Calculate the chisq at points given by indices. Returns ------- chisq: Array Chisq values, one for each sample index, or zero for points below the specified SNR threshold chisq_dof: Array Number of statistical degrees of freedom for the chisq test in the given template, equal to 2 * num_bins - 2 z#%s above chisq activation threshold)r*irr r)NNN) rDrrEabssumlogginginfor rfloat32r^r)repeat)rFr%r=r'rr(r? num_aboveabove above_indices above_snrvZ chisq_outdofrZ_chisqrrrvaluesfs.    zSingleDetPowerChisq.values)rN) __name__ __module__ __qualname____doc__rG staticmethodrQr^rjrrrrrB2s  rBcs0eZdZdZfddZddZddZZS)SingleDetSkyMaxPowerChisqzClass that handles precomputation and memory management for efficiently running the power chisq in a single detector inspiral analysis when maximizing analytically over sky location. c s"tt|jf|d|_d|_dS)N)superrprG template_memcorr_mem)rFkwds) __class__rrrGsz"SingleDetSkyMaxPowerChisq.__init__cCsjt|||j}t|drV|j|jkrVt|j|j}|j}t |j|j|||}nt ||||j}|S)z: Obtain the chisq bins for this template and PSD. rV) r3rQrrXrZrVr[rr\rr)rFr?rrrrrrrrcalculate_chisq_binss  z.SingleDetSkyMaxPowerChisq.calculate_chisq_binsc "Csn|jrft|} |jrlt||jk} | } td| || }|| }|| }tjt|tj d}d}n|}|}| dkr2g}d}d}|j dkst|j t|kstt|t |d|_ |j dkst|j t|kstt|t |d|_ |j }|j }t|j j |j t|j j |j |j j |_|j j |_xt|D]\}}t|g}t||g}||}||||}|}|j|_|j|_||d}|d|| 7}| |d }|| | }||||}|}|||} t| ddd}t||||| |}!||!dqrAobjectrBrprrrrs*   !  R& )a