alg_arc = PyAlgos(windows=winds_arc, mask=mask_arc, pbits=2)
alg_arc.set_peak_selection_pars(npix_min=4, npix_max=500, amax_thr=0, atot_thr=600, son_min=5)

alg_equ = # all parameters are the same as for alg_arc

        #peaks_arc = alg_arc.peak_finder_v2(nda, thr=10, r0=6, dr=0.5)
        #peaks_equ = alg_equ.peak_finder_v3(nda, rank=5, r0=6, dr=0.5) ??????????? finds too many peaks
        peaks_arc = alg_arc.peak_finder_v4(nda, thr_low=10, thr_high=150, rank=5, r0=6, dr=0.5)