...

ImgAlgos/examples/ex_peakfinder_cspad.py

Results

Arc region

Peak data processing script is under development. Results for ARC region:

...

Single-peak events:

Two-peak events: 

New peakfinder more than doubles the number of signal events in the distance distribution.

Equ region

Peak selection in equatorial region:

def peakIsSelectedEqu() :
    """Apply peak selection criteria to each peak from file
    """
    if sp.sonc<9    : return False
    if sp.atot<1800 : return False
    if sp.r<100     : return False
    if sp.r>450     : return False
    return True

Event selection:

def eventIsSelected() :
    """Apply selection criteria to entire event, based on list of peaks 
    """
    sp.event_is_selected = False
    if sp.count_arc_pks_sel  > 2 : return False
    if sp.count_equ_pks_sel  > 5 : return False
    #if sp.count_equ_pks_sel  < 0 : return False
    sp.event_is_selected = True
    return True

Plots:

Image AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage AddedImage Added

References

• class ImgAlgos.PyAlgos - collection of n-d array processing algorithms
• class Detector.PyDetector - direct access to detector data and calibration n-d arrays