Page History
...
Code Block |
---|
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 (vs. fake pixels between panels)
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 |
Overview
Content Tools