def runGUI(filename, camname, imgnum, verbose=False, loglevel=0):
# ======================================
# Set up logging
# ======================================
logger = mt.mylogger(filename='runGUI') # noqa
# ======================================
# Load data and run program
# ======================================
data = E200.E200_load_data(filename, local=False)
# Generate or retrieve qt app
app = mt.qt.get_app()
window = ButterflyGUI(
analyzefcn = analyze_matlab ,
savefcn = save_analysis ,
dataset = data ,
camname = camname ,
imgnum = imgnum ,
verbose = verbose
)
# img = np.random.randn(10, 10)
# window = ButterflyGUI(analyze_matlab, image=img)
# window.ui.imageview_mpl.img.image=img
window.show()
app.exec_()
import analyze_matlab as analyze_matlab
from save_analysis import save_analysis
sets = [
['20140625', '13438'],
['20140625', '13449'],
['20140625', '13450'],
['20140629', '13537']
]
# sets = [['20140625', '13438'],
sets = [['20140629', '13537']]
# sets = [['20140625', '13450']]
logger = mt.mylogger(filename='reprocess')
logger.setLevel(51)
logger.debug('Beginning reprocessing...')
for pair in sets:
setdate = pair[0]
setnumber = pair[1]
loadfile = 'nas/nas-li20-pm00/E200/2014/{}/E200_{}'.format(setdate, setnumber)
data = E200.E200_load_data(loadfile, local=True)
wf = data.write_file
data_wf = wf['data']
processed_wf = data_wf['processed']
arrays_wf = processed_wf['arrays']
vectors_wf = processed_wf['vectors']
#!/usr/bin/env python -m pdb
import E200 as E200
import mytools as mt
mt.mylogger('basictest')
try:
data_a.close() # noqa
data_b.close() # noqa
except:
pass
# ======================================
# Diagnose set A
# ======================================
set_a = ['20140625', '13438']
emit_n_str_a = 'ss_ELANEX_emit_n'
setdate = set_a[0]
setnumber = set_a[1]
loadfile = 'nas/nas-li20-pm00/E200/2014/{}/E200_{}'.format(setdate, setnumber)
data_a = E200.E200_load_data(loadfile)
wf_a = data_a.write_file
processed_a = wf_a['data']['processed']
vectors_a = processed_a['vectors']
arrays_a = processed_a['arrays']
scalars_a = processed_a['scalars']
emit_n_a_str = scalars_a[emit_n_str_a]
emit_n_a = E200.E200_api_getdat(emit_n_a_str)
#!/usr/bin/env python import ButterflyEmittancePython as bt import matplotlib.pyplot as plt import matplotlib.gridspec as gridspec import mytools as mt import numpy as np import os import ElegantPy import shlex import slactrac as sltr import subprocess logger = mt.mylogger(filename='worksheet') # Set up initial beam energy_offset = np.float(3) energy0 = np.float(20.35) gamma = sltr.GeV2gamma(energy0+energy_offset) emitx = np.float(100e-6) emity = emitx QS1_K1_default = 3.8743331090707228e-1 QS2_K1_default = -2.5439067538354171e-1 PEXT_Z = 1994.97 QS1_Z = 1998.71 AL_Z = 2015.16 BE_Z = 1996.34 ELANEX_Z = 2015.22 # IP2QS1_length = 5.4217 IP2QS1_length = QS1_Z-PEXT_Z beam_x = sltr.BeamParams(
# wf.close()
# except:
# pass
sets = [['20140625', '13438'],
# ['20140625', '13449'],
['20140629', '13537'],
['20140625', '13450']]
# sets = [['20140625', '13438'],
sets = [['20140629', '13537']]
# sets = [['20140625', '13450']]
# ======================================
# Set up logging
# ======================================
logger = mt.mylogger(filename='postproc')
logger.critical('Logging set up')
for pair in sets:
setdate = pair[0]
setnumber = pair[1]
loadfile = 'nas/nas-li20-pm00/E200/2014/{}/E200_{}'.format(setdate, setnumber)
data = E200.E200_load_data(loadfile)
wf = data.write_file
processed = wf['data']['processed']
vectors = processed['vectors']
arrays = processed['arrays']
scalars = processed['scalars']
#!/usr/bin/env python
# #!/usr/bin/env python #-m pdb
# import sys
# from IPython.core import ultratb
# sys.excepthook = ultratb.FormattedTB(mode='Verbose',
# color_scheme='Linux', call_pdb=1)
import E200
import matplotlib.pyplot as plt
import matplotlib.gridspec as gs
import mytools as mt
import numpy as np
import scipy.io
logger = mt.mylogger(filename='compare_SCorde')
filebase = 'E200_13450'
dataset = '13438'
filebase = 'E200_{}'.format(dataset)
filename = '{}.mat'.format(filebase)
f = scipy.io.loadmat(filename)
# g = h5.File(filename, 'r', driver='core', backing_store=False)
camname = f[filebase].dtype.names[0]
cam = f[filebase][camname][0, 0]
# ======================================
# Get Sebastien's emittance, UID
