Page History
...
This script can be run in real-time while data is being taken, and will typically complete a few minutes after the run ends. NOTE: when running in parallel, the standard python "break" statement can cause hangs. Use the "break_after" command here to terminate data processing early.
| Code Block |
|---|
from psana import *
dsource = MPIDataSource('exp=xpptut15:run=54:smd')
cspaddet = Detector('cspad')
smldata = dsource.small_data('run54.h5',gather_interval=100)
dsource.break_after(3) # stop iteration after 3 events (break statements do not work reliably with MPIDataSource).
partial_run_sum = None
for nevt,evt in enumerate(dsource.events()):
calib = cspaddet.calib(evt)
if calib is None: continue
cspad_sum = calib.sum() # number
cspad_roi = calib[0][0][3:5] # array
if partial_run_sum is None:
partial_run_sum = cspad_roi
else:
partial_run_sum += cspad_roi
# save per-event data
smldata.event(cspad_sum=cspad_sum,cspad_roi=cspad_roi)
# get (optional) "summary" data
run_sum = smldata.sum(partial_run_sum)
# save HDF5 file, including summary data
smldata.save(run_sum=run_sum) |
...
Overview
Content Tools