Angular Integration

from psana import *
from pypsalg.AngularIntegrationM import *
ds  = DataSource('exp=xpptut15:run=54:smd')
det = Detector('cspad', ds.env())
for nevent,evt in enumerate(ds.events()):
    img = det.image(evt)
    # create a mask of real pixels in the image                                 
    goodPixels = det.image(evt,np.ones_like(det.calib(evt)))
    ai = AngularIntegratorM()
    ai.setParameters(img.shape[0],img.shape[1],
                     mask=goodPixels)
    bins,intensity = ai.getRadialHistogramArrays(img)
    print intensity
    break