Scaling Measurements

import datetime as dt
t = dt.datetime.now()
print('Start at',t.strftime('%H:%M:%S'))

from psana import *
import sys
myrun=sys.argv[1]

t = dt.datetime.now()
print('Import done at',t.strftime('%H:%M:%S'))

dsource = MPIDataSource('exp=cxilu9218:run='+myrun+':smd:dir=/cds/data/drpsrcf/cxi/cxilu9218/xtc')
jungfraudet = Detector('jungfrau4M')

t = dt.datetime.now()
print('Detector done at',t.strftime('%H:%M:%S'))

for nevt,evt in enumerate(dsource.events()):
   calib = jungfraudet.calib(evt)

t = dt.datetime.now()
print('Job done at',t.strftime('%H:%M:%S'))

from psana import *

dsource = MPIDataSource('exp=cxilu9218:run=237:smd:dir=/cds/drpsrcf/cxi/cxilu921
8/xtc/')
jungfraudet = Detector('jungfrau4M')

smldata = dsource.small_data('/cds/drpsrcf/cxi/cxilu9218/scratch/run237.h5',gath
er_interval=100)

partial_run_sum = None
for nevt,evt in enumerate(dsource.events()):

   calib = jungfraudet.calib(evt)
   if calib is None: continue
   jungfrau_sum = calib.sum()      # number
   smldata.event(jungfrau_sum=jungfrau_sum)