Page History
...
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:
Code Block |
---|
[ImgAlgos.NDArrDropletFinder:Arc]
...
threshold_low = 10
threshold_high = 150
peak_radius = 3
... |
Peak list selection parameters:
Remarks
- Peak finder parameters - all parameters in the peak finder are set without optimization. They need to be tuned in two places;
- in the configuration file for two modules
ImgAlgos.NDArrDropletFinder(:Arc/:Equ)
, and - in the python script in calls of
peaks_filter(...)
method.
- in the configuration file for two modules
- ROI for peakfinder is defined by combination of two mechanisms;
- mask - ndarray with values 1/0. Pixels masked by "0", are ignored.
- a set of rectangular windows on sensors, defined by the parameter windows in module ImgAlgos::NDArrDropletFinder.
- Peak finder parameters need to be optimized!
Example: optimization of peak parameters
Selection parameters
Selection of peaks and events is defined in few stages, as listed below.
1. Peak selection parameters in module ImgAlgos.NDArrDropletFinder:
Code Block |
---|
[ImgAlgos.NDArrDropletFinder:Arc]
...
threshold_low = 10
threshold_high = 150
peak_radius = 3
... |
2. Peak list selection parameters:
Code Block |
---|
nda_peaks_arc = peaks_filter(nda_droplets_arc, |
Code Block |
nda_peaks_arc = peaks_filter(nda_droplets_arc, atot_thr=2000, npix_thr=20, npeaks_min=1, npeaks_max=10) |
3. Peak/event selection parameters at peak-list file readout
Code Block |
---|
def peakIsSelected |
Code Block |
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 |
...
Parameter | Sel.1 (rpeak5) | Sel.2 (rpeak7) | Sel.3 (rpeak3) |
---|---|---|---|
in module NDArrDropletFinder:Arc | |||
threshold_low | 10 | 10 | 10 |
threshold_high | 100 | 100 | 150 |
peak_radius | 5 | 7 | 3 |
in method peaks_filter(...) | |||
atot_thr | 2000 | 2000 | 2000 |
npix_thr | 20 | 4030 | 20 |
npeaks_min | 1 | 1 | 1 |
npeaks_max | 10 | 4020 | 10 |
Save selected peaks in file | |||
Process peaks from file | |||
amax | 150 | 150 | 150 |
atot | 2500 | 2500 | 2500 |
npix | 30 | 40 | 20 |
rmin | 434 | 434 | 434 |
rmax | 444 | 444 | 444 |
phi | <-150 & >150 | <-150 & >150 | <-150 & >150 |
Sel.1 (rpeak=5):
Sel.2 (rpeak=7):
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;
- in the configuration file for two modules
ImgAlgos.NDArrDropletFinder(:Arc/:Equ)
, and - in the python script in calls of
peaks_filter(...)
method.
- in the configuration file for two modules
- ROI for peakfinder is defined by combination of two mechanisms;
- mask - ndarray with values 1/0. Pixels masked by "0", are ignored.
- a set of rectangular windows on sensors, defined by the parameter windows in module ImgAlgos::NDArrDropletFinder.
...
References
...
Overview
Content Tools