...

Comparison of the pfv2 with Cheetah list of peaks

Code for comparison: cxif5315/peak-list-comp-cheetah.py

Peak lists for comparison:

• peakfinder: pfv2-cxif5315-r0169-2016-03-28T15:02:47.txt
• cheetah: r0169-cheetah-peaks.txt - with unknown origin...
  

Comparison algorithm and conditions

Use pyimgalgos.TDFileContainer for both peakfinder and cheetah

    fc_ch = TDFileContainer(fncheet, indhdr='frameNumber', objtype=TDCheetahPeakRecord) # , pbits=256)
    fc_pf = TDFileContainer(fnpeaks, indhdr='Evnum', objtype=TDPeakRecord) #, pbits=256)

• Loop over all groups (frames/events) in cheetah
  • Find associated peak-finder group: gr_pf = fc_pf.group(evnum-1)  # accounts for -1 offset in event numeration in Cheetah and pf.
  • count cases when peak-finder group does not exist (peaks not found)
  • discard frames if fidelity do not match
  • discard frames with more than 5 peaks
  • loop over cheetah peaks
    
    • discard cheetah peaks outside ROI mask
    • loop over pf-peaks and check that cheetah peak can be matched with one of them by location (s, r, c)
    • if not-matched - print peak main features
      
      
      

Statistics for number of peaksNumber of peak:

• total in Cheetah: 14005.0  14002 
• matched: 10593
  
• not-matched: 3412 3409
  
• fraction of matched: 0.756.757

Plots for number of peaks in peak-finder, Cheetah, and number of non-matched:matched Cheetah peaks.

Peak-finder number of peaks accounts for arc region.

References

• Radial Background Subtraction Algorithm