Page History

...

python psana-cxif5315-r0169-cspad-ds2-NDArrDropletFinder.py

Example: get table of peaks for different regions

Image with peaks for separate and both regions can be plotted:

...

The 1st line in this file is a (commented) header with colon-titles.

Example: optimization of peak parameters

Peak selection parameters in module ImgAlgos.NDArrDropletFinder:

[ImgAlgos.NDArrDropletFinder:Arc]
...
threshold_low  = 10
threshold_high = 150
peak_radius    = 3
...

Peak list selection parameters:

nda_peaks_arc = peaks_filter(nda_droplets_arc, atot_thr=2000, npix_thr=20, npeaks_min=1, npeaks_max=10)  

def peakIsSelected() :
    """Apply peak selection criteria to each peak from file
    """
    amax, atot, npix, x, y, r, phi = sp.amax, sp.atot, sp.npix, sp.x, sp.y, sp.r, sp.phi
    if amax<150  : return False
    if atot<2500 : return False
    if npix<30   : return False
    if r<434     : return False
    if r>444     : return False
    if phi<150 and phi>-150 : return False
    return True

in module NDArrDropletFinder:Arc

threshold_low

threshold_high

peak_radius

atot_thr

npix_thr

npeaks_min

npeaks_max

peakIsSelected()

amax

atot

2500

2500

2500

npix

rmin

phi

Sel.1 (rpeak=5):

Image AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage Added

Sel.2 (rpeak=7):

Image AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage Added

Sel.3 (rpeak=3):

Remarks

• Peak finder parameters - all parameters in the peak finder are set without optimization. They need to be tuned in two places;
  
  1.  in the configuration file for two modules ImgAlgos.NDArrDropletFinder(:Arc/:Equ), and
  2.  in the python script in calls of peaks_filter(...) method.
• ROI for peakfinder is defined by combination of two mechanisms;
  
  1. mask - ndarray with values 1/0. Pixels masked by "0", are ignored.
  2. a set of rectangular windows on sensors, defined by the parameter windows in module ImgAlgos::NDArrDropletFinder.

...