...

 r0=7  dr=2 (108 pixels)   
 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0
 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0 0
 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0
 0 0 0 1 1 1 0 0 0 0 0 0 0 1 1 1 0 0 0
 0 0 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 0 0
 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 1 1 1 0 0 0 0 0 0 + 0 0 0 0 0 0 1 1 1
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0
 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0
 0 0 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 0 0
 0 0 0 1 1 1 0 0 0 0 0 0 0 1 1 1 0 0 0
 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0
 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0 0 0 0 0
 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0

Test of peak finders

• Test of Peak Finders
• Test of Peak Finders - V2

Photon counting

Photon conversion in pixel detectors quite often happens in the region between pixels. Ionization produced by photon splits between two or rarely more pixels. Photon counting algorithm described here is designed to account for this effect and return array with correct number of photons per pixel. Pythonic API for this algorithm is pretty simple:

# Import
from ImgAlgos.PyAlgos import photons

# Merges photons split among pixels and returns n-d array with integer number of photons per pixel.
nphotons_nda = photons(fphotons, mask)

Sphinx doc

Method photons receives (float) n-d numpy array fphotons representing image intensity in terms of (float) number of photons and associated mask of bad pixels. Both arrays should have the same shape. Two lowest dimensions represent pixel rows and columns in 2-d pixel matrix arrays. Algorithm works with good pixels defined by the mask array (1/0 = good/bad pixel).  Array fphotons is split for two - integer (floor) and leftover fractional number of photon arrays of the same shape. Then Chucks' algorithm is applied to join adjacent pixel intensities in fractional array and combine them in integer number of photons. At the end algorithm returns a sum of integer arrays with corrected photon splitting effect.

References

• ImgAlgos.PyAlgos - code example in Sphinx documentation
• Peak Finding - short announcement about peak finders
• Hit and Peak Finders - examples in Chris' tutorial
• Peak Finding Module - (depricated) psana module, it demonstaration examples and results
• Psana Module Catalog - (depricated) peak finding  psana modules
• Psana Module Examples - (depricated) peak finding examples in psana modules
• GUI for tuning peak finding - Chun's page in development

• pyimgalgos.PeakStore - class helping to save peak parameter records in the text file

• pyimgalgos.TDFileContainer - class helping to retrieve peak parameter records from the text file

• Test of Peak Finders - example of exploitation of peak finders

• Test of Peak Finders - V2 - example of exploitation of peak finders after revision 1 (uniformization)

...

• photons - sphinx doc