Probe files

from opal import load_dataset

ds = load_dataset('./RingCyclotron/', fname='RRI2.h5')

print( ds )

        H5 dataset.

        Number of steps: 1

        Available step attributes (0):



        Available step datasets (10):

        bunchNumber
        id
        px
        py
        pz
        time
        turn
        x
        y
        z

Set the plotting style

from opal.visualization.styles import load_style

load_style('jupyter')

/home/opalci/pyOPALTools/opal/visualization/styles/default.py:186: MatplotlibDeprecationWarning:
The mathtext.fallback_to_cm rcparam was deprecated in Matplotlib 3.3 and will be removed two minor releases later.
  mpl.rcParams['mathtext.fallback_to_cm'] = True
/home/opalci/pyOPALTools/opal/visualization/styles/default.py:186: MatplotlibDeprecationWarning: Support for setting the 'mathtext.fallback_to_cm' rcParam is deprecated since 3.3 and will be removed two minor releases later; use 'mathtext.fallback : 'cm' instead.
  mpl.rcParams['mathtext.fallback_to_cm'] = True
/home/opalci/pyOPALTools/opal/visualization/styles/default.py:186: MatplotlibDeprecationWarning:
The validate_bool_maybe_none function was deprecated in Matplotlib 3.3 and will be removed two minor releases later.
  mpl.rcParams['mathtext.fallback_to_cm'] = True
/home/opalci/pyOPALTools/opal/visualization/styles/default.py:419: MatplotlibDeprecationWarning:
The savefig.jpeg_quality rcparam was deprecated in Matplotlib 3.3 and will be removed two minor releases later.
  mpl.rcParams['savefig.jpeg_quality']    = 95        # when a jpeg is saved, the default quality parameter.

Do some plotting …

plt = ds.plot_probe_histogram(bins=300)
plt.show()
plt.close()

from opal.analysis.cyclotron import eval_radius
import matplotlib.pyplot as plt

t = ds.getData('turn')

x = ds.getData('x')
y = ds.getData('y')

# second peak of above figure
x = x[t == 1]
y = y[t == 1]

r = eval_radius(x, y)

plt.hist(r, bins=400)
plt.show()
plt.close()