Page History
...
2017-02-27-cxif5315-r0169-bad-events.txt
Comparison of analysis script with psocake
Code Block | ||||
---|---|---|---|---|
| ||||
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) |
Note |
---|
Conclusion: within mask-permitted region analysis script finds the same peaks as psocake |
Overview
Content Tools