Page History
...
| Code Block |
|---|
alg_arc = PyAlgos(windows=winds_arc, mask=mask_arc, pbits=2)
alg_arc.set_peak_selection_pars(npix_min=0, npix_max=1e6, amax_thr=0, atot_thr=0, son_min=6)
#alg_arc.set_peak_selection_pars(npix_min=0, npix_max=1e6, amax_thr=0, atot_thr=500, son_min=6) # for v2r1
alg_equ = ... the same
# in the event loop
#peaks_arc = alg_arc.peak_finder_v2r1(nda, thr=30, r0=7, dr=2)
#peaks_arc = alg_arc.peak_finder_v3r1(nda, rank=5, r0=7, dr=2)
peaks_arc = alg_arc.peak_finder_v4r1(nda, thr_low=10, thr_high=150, rank=5, r0=7, dr=2) |
Results
ARC
pre-selection for histograms
| Code Block |
|---|
def procPeakDataArc(pk) :
""" Process peak for ARC region; accumulate peak statistics in histogram arrays.
"""
#===================
sp.lst_arc_atot.append(pk.atot)
if pk.atot<1500 : return
#===================
sp.lst_arc_amax.append(pk.amax)
sp.lst_arc_npix.append(pk.npix)
sp.lst_arc_r .append(pk.r)
sp.lst_arc_phi .append(pk.phi)
... |
plots for pfv2
...
/3/4r1
pfv3
pfv4
EQU
EQU
pre-selection for histograms
| Code Block |
|---|
def procPeakDataEqu(pk) :
""" Process peak for EQU region; accumulate peak data
"""
# discard small radius peaks from all histograms
#===================
sp.lst_equ_atot.append(pk.atot)
if pk.atot<1500 : return
if pk.r<100 : return
#===================
sp.lst_equ_amax.append(pk.amax)
sp.lst_equ_npix.append(pk.npix)
sp.lst_equ_r .append(pk.r)
sp.lst_equ_phi .append(pk.phi)
... |
plots for pfv2/3/4r1
Overview
Content Tools