Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Code Block
titlePeak reconstruction/selection parameters
peaks_equ = alg_equ.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=4, r0=6, dr=2)
alg_equ.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)

 

Image Removed

Note

For further histograms discard peaks with pk.r<105

Image Added

Note

For further histograms discard peaks with pk.atot<600

...

Tentative peak selection criteria

Code Block
alg_arc.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)
peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=8, r0=10, dr=2)

Image Removed

Note

Further discard peaks with pk.atot<1000

 

 

Image RemovedImage RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage RemovedImage Removed

Image RemovedImage RemovedImage Removed

 

 

Image RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage Removed

 

 

 

 

...

        if pk.atot<600  : return False
        if pk.npix>60   : return False
        if pk.rms<10    : return False
        if pk.rms>25    : return False
        if math.fabs(pk.bkgd)>15 : return False
pk.npix ?
pk.son ?
sigma_hov ?

ARC

Code Block
alg_arc.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)
peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=8, r0=10, dr=2)

Image Added

Note

Further discard peaks with pk.atot<1000

Image AddedImage AddedImage Added

Image AddedImage Added

Image AddedImage AddedImage Added

Image AddedImage AddedImage Added

Image AddedImage Added

Image AddedImage Added

Image AddedImage Added

 

Tentative peak selection criteria

 

Code Block
titleTentative peak selection criteria
        if pk.atot<1000 : return False
        if pk.npix>200  : return False
        if pk.r<435     : return False
        if pk.r>443     : return False
        if pk.rms<10    : return False
        if pk.rms>30    : return False
        if math.fabs(pk.bkgd)>15: return False

pk.npix ?
pk.son ?
sigma_hov ?