Page History

    alg_arc = PyAlgos(windows=winds_arc, mask=mask_arc, pbits=2)
    alg_arc.set_peak_selection_pars(npix_min=20, npix_max=1000, amax_thr=0, atot_thr=2000, son_min=5)

    alg_equ = PyAlgos(windows=winds_equ, mask=mask_equ, pbits=0)
    alg_equ.set_peak_selection_pars(npix_min=5, npix_max=1000, amax_thr=0, atot_thr=600, son_min=5)

    for i, evt in enumerate(ds.events()) :
        nda_data = det.calib(evt, cmpars=(5,50))

        if nda_data is not None :
            nda =  np.array(nda_data, dtype=np.float32, copy=True)

            if BKG_MODE == BKG_MODE_NORMAL : 
                nda = subtract_bkgd(nda, nda_bkgd, mask=nda_smask, winds=winds_bkgd, pbits=0)

            peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=160, rank=8, r0=10, dr=2)
            peaks_equ = alg_equ.peak_finder_v4r2(nda, thr_low=60, thr_high=100, rank=4, r0=6, dr=2)

Conclusion: within mask-permitted region analysis script finds the same peaks as psocake