RF-Track 2.7.2__cp313-cp313-win_amd64.whl

RF_Track.py

@@ -0,0 +1,4509 @@
+# This file was automatically generated by SWIG (https://www.swig.org).
+# Version 4.4.1
+#
+# Do not make changes to this file unless you know what you are doing - modify
+# the SWIG interface file instead.
+
+from sys import version_info as _swig_python_version_info
+# Import the low-level C/C++ module
+if getattr(globals().get("__spec__"), "parent", None) or __package__ or "." in __name__:
+    from . import _RF_Track
+else:
+    import _RF_Track
+
+try:
+    import builtins as __builtin__
+except ImportError:
+    import __builtin__
+
+_swig_new_instance_method = _RF_Track.SWIG_PyInstanceMethod_New
+_swig_new_static_method = _RF_Track.SWIG_PyStaticMethod_New
+
+def _swig_repr(self):
+    try:
+        strthis = "proxy of " + self.this.__repr__()
+    except __builtin__.Exception:
+        strthis = ""
+    return "<%s.%s; %s >" % (self.__class__.__module__, self.__class__.__name__, strthis,)
+
+
+def _swig_setattr_nondynamic_instance_variable(set):
+    def set_instance_attr(self, name, value):
+        if name == "this":
+            set(self, name, value)
+        elif name == "thisown":
+            self.this.own(value)
+        elif hasattr(self, name) and isinstance(getattr(type(self), name), property):
+            set(self, name, value)
+        else:
+            raise AttributeError("You cannot add instance attributes to %s" % self)
+    return set_instance_attr
+
+
+def _swig_setattr_nondynamic_class_variable(set):
+    def set_class_attr(cls, name, value):
+        if hasattr(cls, name) and not isinstance(getattr(cls, name), property):
+            set(cls, name, value)
+        else:
+            raise AttributeError("You cannot add class attributes to %s" % cls)
+    return set_class_attr
+
+
+def _swig_add_metaclass(metaclass):
+    """Class decorator for adding a metaclass to a SWIG wrapped class - a slimmed down version of six.add_metaclass"""
+    def wrapper(cls):
+        return metaclass(cls.__name__, cls.__bases__, cls.__dict__.copy())
+    return wrapper
+
+
+class _SwigNonDynamicMeta(type):
+    """Meta class to enforce nondynamic attributes (no new attributes) for a class"""
+    __setattr__ = _swig_setattr_nondynamic_class_variable(type.__setattr__)
+
+
+import weakref
+
+SHARED_PTR_DISOWN = _RF_Track.SHARED_PTR_DISOWN
+
+RF_TRACK_VERSION = _RF_Track.RF_TRACK_VERSION
+
+relativistic_velocity_addition = _RF_Track.relativistic_velocity_addition
+lorentz_boost = _RF_Track.lorentz_boost
+lorentz_boost_matrix = _RF_Track.lorentz_boost_matrix
+rng_set = _RF_Track.rng_set
+rng_set_seed = _RF_Track.rng_set_seed
+rng_get = _RF_Track.rng_get
+class Frame(object):
+    r"""
+
+    A rigid coordinate frame with an origin and rotated axes. Within RF-Track, origins are
+    in mm and angles in rad. Multiplication composes frames or transforms a local point into
+    its parent frame; division applies the inverse transform.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Frame self, double x, double y, double z, double phi_x, double theta_z, double phi_z) -> Frame
+        __init__(Frame self, double x, double y, double z, Rotation const & _axes) -> Frame
+        __init__(Frame self, Vector3d const & _origin, double phi_x, double theta_z, double phi_z) -> Frame
+        __init__(Frame self, Vector3d const & _origin, Rotation const & _axes) -> Frame
+        __init__(Frame self) -> Frame
+        __init__(Frame self, Vector3d const & _origin) -> Frame
+        __init__(Frame self, Rotation const & _axes) -> Frame
+
+        A rigid coordinate frame with an origin and rotated axes. Within RF-Track, origins are
+        in mm and angles in rad. Multiplication composes frames or transforms a local point into
+        its parent frame; division applies the inverse transform.
+
+
+        """
+        _RF_Track.Frame_swiginit(self, _RF_Track.new_Frame(*args))
+    __swig_destroy__ = _RF_Track.delete_Frame
+    get_origin = _swig_new_instance_method(_RF_Track.Frame_get_origin)
+    get_nx = _swig_new_instance_method(_RF_Track.Frame_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.Frame_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.Frame_get_nz)
+    __imul__ = _swig_new_instance_method(_RF_Track.Frame___imul__)
+
+    def __itruediv__(self, *args):
+        return _RF_Track.Frame___itruediv__(self, *args)
+    __idiv__ = __itruediv__
+
+
+    __eq__ = _swig_new_instance_method(_RF_Track.Frame___eq__)
+    __ne__ = _swig_new_instance_method(_RF_Track.Frame___ne__)
+    __mul__ = _swig_new_instance_method(_RF_Track.Frame___mul__)
+    __truediv__ = _swig_new_instance_method(_RF_Track.Frame___truediv__)
+    __rtruediv__ = _swig_new_instance_method(_RF_Track.Frame___rtruediv__)
+
+# Register Frame in _RF_Track:
+_RF_Track.Frame_swigregister(Frame)
+cvar = _RF_Track.cvar
+version = cvar.version
+max_number_of_threads = cvar.max_number_of_threads
+clight = cvar.clight
+muonmass = cvar.muonmass
+protonmass = cvar.protonmass
+electronmass = cvar.electronmass
+muonlifetime = cvar.muonlifetime
+s = cvar.s
+ms = cvar.ms
+us = cvar.us
+ns = cvar.ns
+ps = cvar.ps
+fs = cvar.fs
+C = cvar.C
+mC = cvar.mC
+uC = cvar.uC
+nC = cvar.nC
+pC = cvar.pC
+electron_anomalous_magnetic_moment = cvar.electron_anomalous_magnetic_moment
+proton_anomalous_magnetic_moment = cvar.proton_anomalous_magnetic_moment
+muon_anomalous_magnetic_moment = cvar.muon_anomalous_magnetic_moment
+
+class Parallel_ODE_Solver(object):
+    r"""
+
+    ODE integration controls inherited by numerically tracked field elements. Select an
+    algorithm with set_odeint_algorithm and GSL tolerances with set_odeint_epsabs and
+    set_odeint_epsrel. Volume tracking uses its TrackingOptions controls.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_Parallel_ODE_Solver
+    set_odeint_epsabs = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_set_odeint_epsabs)
+    set_odeint_epsrel = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_set_odeint_epsrel)
+    get_odeint_epsabs = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_get_odeint_epsabs)
+    get_odeint_epsrel = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_get_odeint_epsrel)
+    get_odeint_algorithm = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_get_odeint_algorithm)
+    set_odeint_algorithm = _swig_new_instance_method(_RF_Track.Parallel_ODE_Solver_set_odeint_algorithm)
+
+# Register Parallel_ODE_Solver in _RF_Track:
+_RF_Track.Parallel_ODE_Solver_swigregister(Parallel_ODE_Solver)
+class TrackingOptions(object):
+    r"""
+
+    Time-integration settings for Volume.track and Volume.btrack. Times and step intervals
+    ending in _mm are in mm/c. Configure a TrackingOptions object and pass it as the second
+    tracking argument. Defaults are given by the public tracking_options.hh declaration.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    odeint_algorithm = property(_RF_Track.TrackingOptions_odeint_algorithm_get, _RF_Track.TrackingOptions_odeint_algorithm_set, doc=r"""
+    odeint_algorithm : std::string
+
+    ODE integration algorithm; default rk2.
+
+
+    """)
+    odeint_epsabs = property(_RF_Track.TrackingOptions_odeint_epsabs_get, _RF_Track.TrackingOptions_odeint_epsabs_set, doc=r"""
+    odeint_epsabs : double
+
+    Absolute error tolerance for GSL integration algorithms; default 0.001.
+
+
+    """)
+    odeint_epsrel = property(_RF_Track.TrackingOptions_odeint_epsrel_get, _RF_Track.TrackingOptions_odeint_epsrel_set, doc=r"""
+    odeint_epsrel : double
+
+    Relative error tolerance for GSL integration algorithms; default 0.
+
+
+    """)
+    backtrack_at_entrance = property(_RF_Track.TrackingOptions_backtrack_at_entrance_get, _RF_Track.TrackingOptions_backtrack_at_entrance_set, doc=r"""
+    backtrack_at_entrance : bool
+
+    At tracking start, backtrack the bunch so its foremost particle is at the entrance; default false.
+
+
+    """)
+    verbosity = property(_RF_Track.TrackingOptions_verbosity_get, _RF_Track.TrackingOptions_verbosity_set, doc=r"""
+    verbosity : int
+
+    Diagnostic output level; 0 is silent.
+
+
+    """)
+    dt_mm = property(_RF_Track.TrackingOptions_dt_mm_get, _RF_Track.TrackingOptions_dt_mm_set, doc=r"""
+    dt_mm : double
+
+    Integration time step in mm/c; default infinity.
+
+
+    """)
+    t_max_mm = property(_RF_Track.TrackingOptions_t_max_mm_get, _RF_Track.TrackingOptions_t_max_mm_set, doc=r"""
+    t_max_mm : double
+
+    Upper tracking-time limit in mm/c; default infinity.
+
+
+    """)
+    t_min_mm = property(_RF_Track.TrackingOptions_t_min_mm_get, _RF_Track.TrackingOptions_t_min_mm_set, doc=r"""
+    t_min_mm : double
+
+    Lower tracking-time limit in mm/c for backward tracking; default minus infinity.
+
+
+    """)
+    tt_dt_mm = property(_RF_Track.TrackingOptions_tt_dt_mm_get, _RF_Track.TrackingOptions_tt_dt_mm_set, doc=r"""
+    tt_dt_mm : double
+
+    Time interval between transport-table samples in mm/c; default 0.
+
+
+    """)
+    tt_select = property(_RF_Track.TrackingOptions_tt_select_get, _RF_Track.TrackingOptions_tt_select_set, doc=r"""
+    tt_select : std::string
+
+    Particles used for transport statistics: all, all_in_volume, active, or
+    active_in_volume. Active particles already exist. Default all.
+
+
+    """)
+    sc_dt_mm = property(_RF_Track.TrackingOptions_sc_dt_mm_get, _RF_Track.TrackingOptions_sc_dt_mm_set, doc=r"""
+    sc_dt_mm : double
+
+    Time interval between space-charge kicks in mm/c; default 0.
+
+
+    """)
+    cfx_dt_mm = property(_RF_Track.TrackingOptions_cfx_dt_mm_get, _RF_Track.TrackingOptions_cfx_dt_mm_set, doc=r"""
+    cfx_dt_mm : double
+
+    Time interval between other collective-effect kicks in mm/c; default 0.
+
+
+    """)
+    emission_nsteps = property(_RF_Track.TrackingOptions_emission_nsteps_get, _RF_Track.TrackingOptions_emission_nsteps_set, doc=r"""
+    emission_nsteps : size_t
+
+    Number of emission steps; default 10.
+
+
+    """)
+    emission_range = property(_RF_Track.TrackingOptions_emission_range_get, _RF_Track.TrackingOptions_emission_range_set, doc=r"""
+    emission_range : double
+
+    Range over which reduced emission space-charge steps apply; default 10.
+
+
+    """)
+    wp_dt_mm = property(_RF_Track.TrackingOptions_wp_dt_mm_get, _RF_Track.TrackingOptions_wp_dt_mm_set, doc=r"""
+    wp_dt_mm : double
+
+    Time interval between saved beam watchpoints in mm/c; default 0.
+
+
+    """)
+    wp_basename = property(_RF_Track.TrackingOptions_wp_basename_get, _RF_Track.TrackingOptions_wp_basename_set, doc=r"""
+    wp_basename : std::string
+
+    Base filename for watchpoint distributions; default watch_beam.
+
+
+    """)
+    wp_gzip = property(_RF_Track.TrackingOptions_wp_gzip_get, _RF_Track.TrackingOptions_wp_gzip_set, doc=r"""
+    wp_gzip : bool
+
+    Compress watchpoint files with gzip; default false.
+
+
+    """)
+
+    def __init__(self):
+        r"""
+        __init__(TrackingOptions self) -> TrackingOptions
+
+        Time-integration settings for Volume.track and Volume.btrack. Times and step intervals
+        ending in _mm are in mm/c. Configure a TrackingOptions object and pass it as the second
+        tracking argument. Defaults are given by the public tracking_options.hh declaration.
+
+
+        """
+        _RF_Track.TrackingOptions_swiginit(self, _RF_Track.new_TrackingOptions())
+    __swig_destroy__ = _RF_Track.delete_TrackingOptions
+
+# Register TrackingOptions in _RF_Track:
+_RF_Track.TrackingOptions_swigregister(TrackingOptions)
+class Element(object):
+    r"""
+
+    Base interface for beamline elements. Provides naming, placement, apertures,
+    tracking-step settings, electromagnetic-field queries and collective effects. Geometry
+    lengths are in metres; field-query coordinates are in millimetres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined - class is abstract")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_Element
+    clone = _swig_new_instance_method(_RF_Track.Element_clone)
+    is_point_inside_aperture = _swig_new_instance_method(_RF_Track.Element_is_point_inside_aperture)
+    is_particle_inside_aperture = _swig_new_instance_method(_RF_Track.Element_is_particle_inside_aperture)
+    get_aperture_area = _swig_new_instance_method(_RF_Track.Element_get_aperture_area)
+    get_aperture_x = _swig_new_instance_method(_RF_Track.Element_get_aperture_x)
+    get_aperture_y = _swig_new_instance_method(_RF_Track.Element_get_aperture_y)
+    get_aperture_shape = _swig_new_instance_method(_RF_Track.Element_get_aperture_shape)
+    set_aperture_x = _swig_new_instance_method(_RF_Track.Element_set_aperture_x)
+    set_aperture_y = _swig_new_instance_method(_RF_Track.Element_set_aperture_y)
+    set_aperture = _swig_new_instance_method(_RF_Track.Element_set_aperture)
+    set_aperture_shape = _swig_new_instance_method(_RF_Track.Element_set_aperture_shape)
+    set_width = _swig_new_instance_method(_RF_Track.Element_set_width)
+    set_height = _swig_new_instance_method(_RF_Track.Element_set_height)
+    get_width = _swig_new_instance_method(_RF_Track.Element_get_width)
+    get_height = _swig_new_instance_method(_RF_Track.Element_get_height)
+    __assign__ = _swig_new_instance_method(_RF_Track.Element___assign__)
+    set_offsets = _swig_new_instance_method(_RF_Track.Element_set_offsets)
+    scatter_element = _swig_new_instance_method(_RF_Track.Element_scatter_element)
+    get_offsets = _swig_new_instance_method(_RF_Track.Element_get_offsets)
+    get_name = _swig_new_instance_method(_RF_Track.Element_get_name)
+    get_nsteps = _swig_new_instance_method(_RF_Track.Element_get_nsteps)
+    get_tt_nsteps = _swig_new_instance_method(_RF_Track.Element_get_tt_nsteps)
+    get_sc_nsteps = _swig_new_instance_method(_RF_Track.Element_get_sc_nsteps)
+    get_cfx_nsteps = _swig_new_instance_method(_RF_Track.Element_get_cfx_nsteps)
+    get_length = _swig_new_instance_method(_RF_Track.Element_get_length)
+    set_name = _swig_new_instance_method(_RF_Track.Element_set_name)
+    set_tt_nsteps = _swig_new_instance_method(_RF_Track.Element_set_tt_nsteps)
+    set_sc_nsteps = _swig_new_instance_method(_RF_Track.Element_set_sc_nsteps)
+    set_cfx_nsteps = _swig_new_instance_method(_RF_Track.Element_set_cfx_nsteps)
+    set_nsteps = _swig_new_instance_method(_RF_Track.Element_set_nsteps)
+    get_field = _swig_new_instance_method(_RF_Track.Element_get_field)
+    set_sc_engine = _swig_new_instance_method(_RF_Track.Element_set_sc_engine)
+    add_collective_effect = _swig_new_instance_method(_RF_Track.Element_add_collective_effect)
+    clear_collective_effects = _swig_new_instance_method(_RF_Track.Element_clear_collective_effects)
+    get_collective_effects = _swig_new_instance_method(_RF_Track.Element_get_collective_effects)
+    get_S = _swig_new_instance_method(_RF_Track.Element_get_S)
+    remove = _swig_new_instance_method(_RF_Track.Element_remove)
+    replace_with = _swig_new_instance_method(_RF_Track.Element_replace_with)
+    insert = _swig_new_instance_method(_RF_Track.Element_insert)
+    insert_ref = _swig_new_instance_method(_RF_Track.Element_insert_ref)
+    insert_after = _swig_new_instance_method(_RF_Track.Element_insert_after)
+    insert_ref_after = _swig_new_instance_method(_RF_Track.Element_insert_ref_after)
+
+# Register Element in _RF_Track:
+_RF_Track.Element_swigregister(Element)
+class Lattice(Element):
+    r"""
+
+    An ordered beamline for tracking Bunch6d or Beam in longitudinal position S.
+
+    Construct an empty lattice, copy a lattice, or import a MAD-X Twiss file.
+    append copies an element; append_ref shares it. Appending a lattice keeps
+    it as a nested element; insert flattens its elements into this lattice.
+    Python indexing starts at L[0]; Octave uses L{1}. String indexing selects
+    names with wildcards. track returns the transported beam or bunch.
+
+    Python example:
+        L = RF_Track.Lattice()
+        L.append(RF_Track.Drift(1.0))
+        B1 = L.track(B0)
+    Octave example:
+        L = Lattice();
+        L.append(Drift(1.0));
+        B1 = L.track(B0);
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    verbosity = property(_RF_Track.Lattice_verbosity_get, _RF_Track.Lattice_verbosity_set, doc=r"""verbosity : int""")
+    __swig_destroy__ = _RF_Track.delete_Lattice
+
+    def __init__(self, *args):
+        r"""
+        __init__(Lattice self) -> Lattice
+        __init__(Lattice self, Lattice lattice) -> Lattice
+        __init__(Lattice self, std::string const & madx_twiss_file) -> Lattice
+
+        An ordered beamline for tracking Bunch6d or Beam in longitudinal position S.
+
+        Construct an empty lattice, copy a lattice, or import a MAD-X Twiss file.
+        append copies an element; append_ref shares it. Appending a lattice keeps
+        it as a nested element; insert flattens its elements into this lattice.
+        Python indexing starts at L[0]; Octave uses L{1}. String indexing selects
+        names with wildcards. track returns the transported beam or bunch.
+
+        Python example:
+            L = RF_Track.Lattice()
+            L.append(RF_Track.Drift(1.0))
+            B1 = L.track(B0)
+        Octave example:
+            L = Lattice();
+            L.append(Drift(1.0));
+            B1 = L.track(B0);
+
+
+        """
+        _RF_Track.Lattice_swiginit(self, _RF_Track.new_Lattice(*args))
+    __assign__ = _swig_new_instance_method(_RF_Track.Lattice___assign__)
+    append = _swig_new_instance_method(_RF_Track.Lattice_append)
+    append_ref = _swig_new_instance_method(_RF_Track.Lattice_append_ref)
+    insert = _swig_new_instance_method(_RF_Track.Lattice_insert)
+    insert_ref = _swig_new_instance_method(_RF_Track.Lattice_insert_ref)
+    accept = _swig_new_instance_method(_RF_Track.Lattice_accept)
+    get_transport_table = _swig_new_instance_method(_RF_Track.Lattice_get_transport_table)
+    get_lost_particles = _swig_new_instance_method(_RF_Track.Lattice_get_lost_particles)
+    size = _swig_new_instance_method(_RF_Track.Lattice_size)
+    track = _swig_new_instance_method(_RF_Track.Lattice_track)
+    btrack = _swig_new_instance_method(_RF_Track.Lattice_btrack)
+    set_t0 = _swig_new_instance_method(_RF_Track.Lattice_set_t0)
+    unset_t0 = _swig_new_instance_method(_RF_Track.Lattice_unset_t0)
+    distribute_sc_nsteps = _swig_new_instance_method(_RF_Track.Lattice_distribute_sc_nsteps)
+    distribute_cfx_nsteps = _swig_new_instance_method(_RF_Track.Lattice_distribute_cfx_nsteps)
+    autophase = _swig_new_instance_method(_RF_Track.Lattice_autophase)
+    get_ncorrectors = _swig_new_instance_method(_RF_Track.Lattice_get_ncorrectors)
+    get_nbpms = _swig_new_instance_method(_RF_Track.Lattice_get_nbpms)
+    get_rf_elements = _swig_new_instance_method(_RF_Track.Lattice_get_rf_elements)
+    get_quadrupoles = _swig_new_instance_method(_RF_Track.Lattice_get_quadrupoles)
+    get_correctors = _swig_new_instance_method(_RF_Track.Lattice_get_correctors)
+    get_multipoles = _swig_new_instance_method(_RF_Track.Lattice_get_multipoles)
+    get_sextupoles = _swig_new_instance_method(_RF_Track.Lattice_get_sextupoles)
+    get_solenoids = _swig_new_instance_method(_RF_Track.Lattice_get_solenoids)
+    get_lattices = _swig_new_instance_method(_RF_Track.Lattice_get_lattices)
+    get_volumes = _swig_new_instance_method(_RF_Track.Lattice_get_volumes)
+    get_screens = _swig_new_instance_method(_RF_Track.Lattice_get_screens)
+    get_sbends = _swig_new_instance_method(_RF_Track.Lattice_get_sbends)
+    get_drifts = _swig_new_instance_method(_RF_Track.Lattice_get_drifts)
+    get_coils = _swig_new_instance_method(_RF_Track.Lattice_get_coils)
+    get_bpms = _swig_new_instance_method(_RF_Track.Lattice_get_bpms)
+    get_elements_by_name = _swig_new_instance_method(_RF_Track.Lattice_get_elements_by_name)
+    get_beam_at_screens = _swig_new_instance_method(_RF_Track.Lattice_get_beam_at_screens)
+    get_bunch_at_screens = _swig_new_instance_method(_RF_Track.Lattice_get_bunch_at_screens)
+    set_bpm_resolution = _swig_new_instance_method(_RF_Track.Lattice_set_bpm_resolution)
+    get_bpm_resolution = _swig_new_instance_method(_RF_Track.Lattice_get_bpm_resolution)
+    get_bpm_readings = _swig_new_instance_method(_RF_Track.Lattice_get_bpm_readings)
+    get_correctors_strengths = _swig_new_instance_method(_RF_Track.Lattice_get_correctors_strengths)
+    set_correctors_strengths = _swig_new_instance_method(_RF_Track.Lattice_set_correctors_strengths)
+    vary_correctors_strengths = _swig_new_instance_method(_RF_Track.Lattice_vary_correctors_strengths)
+    reset_correctors = _swig_new_instance_method(_RF_Track.Lattice_reset_correctors)
+    get_response_matrix = _swig_new_instance_method(_RF_Track.Lattice_get_response_matrix)
+    scatter_elements = _swig_new_instance_method(_RF_Track.Lattice_scatter_elements)
+    align_elements = _swig_new_instance_method(_RF_Track.Lattice_align_elements)
+    get_elements_offsets = _swig_new_instance_method(_RF_Track.Lattice_get_elements_offsets)
+    set_elements_offsets = _swig_new_instance_method(_RF_Track.Lattice_set_elements_offsets)
+    orbit_correction = _swig_new_instance_method(_RF_Track.Lattice_orbit_correction)
+    __getitem__ = _swig_new_instance_method(_RF_Track.Lattice___getitem__)
+
+# Register Lattice in _RF_Track:
+_RF_Track.Lattice_swigregister(Lattice)
+class Lattice_view(object):
+    r"""
+
+    A view of a lattice between selected beginning and ending elements. The view refers to
+    the underlying lattice rather than duplicating its elements.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, lattice_ptr, begin, end):
+        r"""
+        __init__(Lattice_view self, std::shared_ptr< Lattice > lattice_ptr, std::shared_ptr< Element > begin, std::shared_ptr< Element > end) -> Lattice_view
+
+        A view of a lattice between selected beginning and ending elements. The view refers to
+        the underlying lattice rather than duplicating its elements.
+
+
+        """
+        _RF_Track.Lattice_view_swiginit(self, _RF_Track.new_Lattice_view(lattice_ptr, begin, end))
+    track = _swig_new_instance_method(_RF_Track.Lattice_view_track)
+    autophase = _swig_new_instance_method(_RF_Track.Lattice_view_autophase)
+    get_transport_table = _swig_new_instance_method(_RF_Track.Lattice_view_get_transport_table)
+    __swig_destroy__ = _RF_Track.delete_Lattice_view
+
+# Register Lattice_view in _RF_Track:
+_RF_Track.Lattice_view_swigregister(Lattice_view)
+class Volume(Element, TrackingOptions):
+    r"""
+
+    A three-dimensional arrangement of electromagnetic elements for time tracking.
+
+    Overlapping element fields are superposed. add copies elements and add_ref
+    shares them; positions are in metres and placement angles in radians.
+    A Lattice or Volume added to a Volume is expanded into individual elements.
+    Use Bunch6dT or BeamT to retain a simultaneous spatial distribution.
+    Bunch6d/Beam tracking overloads support passage between boundary planes.
+    Set boundary planes with set_s0 and set_s1, and integration controls with
+    TrackingOptions. track returns a transported copy.
+
+    Python example:
+        V = RF_Track.Volume()
+        V.add(RF_Track.Drift(1.0), 0., 0., 0.)
+        V.set_s0(0.); V.set_s1(1.)
+        B1 = V.track(B0)
+    Octave example:
+        V = Volume();
+        V.add(Drift(1.0), 0, 0, 0);
+        V.set_s0(0); V.set_s1(1);
+        B1 = V.track(B0);
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_Volume
+
+    def __init__(self, *args):
+        r"""
+        __init__(Volume self) -> Volume
+        __init__(Volume self, Volume v) -> Volume
+
+        A three-dimensional arrangement of electromagnetic elements for time tracking.
+
+        Overlapping element fields are superposed. add copies elements and add_ref
+        shares them; positions are in metres and placement angles in radians.
+        A Lattice or Volume added to a Volume is expanded into individual elements.
+        Use Bunch6dT or BeamT to retain a simultaneous spatial distribution.
+        Bunch6d/Beam tracking overloads support passage between boundary planes.
+        Set boundary planes with set_s0 and set_s1, and integration controls with
+        TrackingOptions. track returns a transported copy.
+
+        Python example:
+            V = RF_Track.Volume()
+            V.add(RF_Track.Drift(1.0), 0., 0., 0.)
+            V.set_s0(0.); V.set_s1(1.)
+            B1 = V.track(B0)
+        Octave example:
+            V = Volume();
+            V.add(Drift(1.0), 0, 0, 0);
+            V.set_s0(0); V.set_s1(1);
+            B1 = V.track(B0);
+
+
+        """
+        _RF_Track.Volume_swiginit(self, _RF_Track.new_Volume(*args))
+    __assign__ = _swig_new_instance_method(_RF_Track.Volume___assign__)
+    add = _swig_new_instance_method(_RF_Track.Volume_add)
+    add_ref = _swig_new_instance_method(_RF_Track.Volume_add_ref)
+    size = _swig_new_instance_method(_RF_Track.Volume_size)
+    clear = _swig_new_instance_method(_RF_Track.Volume_clear)
+    btrack = _swig_new_instance_method(_RF_Track.Volume_btrack)
+    track = _swig_new_instance_method(_RF_Track.Volume_track)
+    accept = _swig_new_instance_method(_RF_Track.Volume_accept)
+    get_transport_table = _swig_new_instance_method(_RF_Track.Volume_get_transport_table)
+    get_transport_table_from_screens = _swig_new_instance_method(_RF_Track.Volume_get_transport_table_from_screens)
+    get_lost_particles = _swig_new_instance_method(_RF_Track.Volume_get_lost_particles)
+    get_bunch_at_s0 = _swig_new_instance_method(_RF_Track.Volume_get_bunch_at_s0)
+    get_bunch_at_s1 = _swig_new_instance_method(_RF_Track.Volume_get_bunch_at_s1)
+    get_beam_at_s0 = _swig_new_instance_method(_RF_Track.Volume_get_beam_at_s0)
+    get_beam_at_s1 = _swig_new_instance_method(_RF_Track.Volume_get_beam_at_s1)
+    distribute_screens_along_path = _swig_new_instance_method(_RF_Track.Volume_distribute_screens_along_path)
+    get_beam_at_screens = _swig_new_instance_method(_RF_Track.Volume_get_beam_at_screens)
+    get_bunch_at_screens = _swig_new_instance_method(_RF_Track.Volume_get_bunch_at_screens)
+    get_beam_at_tt_screens = _swig_new_instance_method(_RF_Track.Volume_get_beam_at_tt_screens)
+    get_bunch_at_tt_screens = _swig_new_instance_method(_RF_Track.Volume_get_bunch_at_tt_screens)
+    get_path_length = _swig_new_instance_method(_RF_Track.Volume_get_path_length)
+    get_rf_elements = _swig_new_instance_method(_RF_Track.Volume_get_rf_elements)
+    get_quadrupoles = _swig_new_instance_method(_RF_Track.Volume_get_quadrupoles)
+    get_multipoles = _swig_new_instance_method(_RF_Track.Volume_get_multipoles)
+    get_sextupoles = _swig_new_instance_method(_RF_Track.Volume_get_sextupoles)
+    get_correctors = _swig_new_instance_method(_RF_Track.Volume_get_correctors)
+    get_solenoids = _swig_new_instance_method(_RF_Track.Volume_get_solenoids)
+    get_screens = _swig_new_instance_method(_RF_Track.Volume_get_screens)
+    get_sbends = _swig_new_instance_method(_RF_Track.Volume_get_sbends)
+    get_drifts = _swig_new_instance_method(_RF_Track.Volume_get_drifts)
+    get_coils = _swig_new_instance_method(_RF_Track.Volume_get_coils)
+    get_bpms = _swig_new_instance_method(_RF_Track.Volume_get_bpms)
+    get_ncorrectors = _swig_new_instance_method(_RF_Track.Volume_get_ncorrectors)
+    get_nbpms = _swig_new_instance_method(_RF_Track.Volume_get_nbpms)
+    get_elements_by_name = _swig_new_instance_method(_RF_Track.Volume_get_elements_by_name)
+    scatter_elements = _swig_new_instance_method(_RF_Track.Volume_scatter_elements)
+    align_elements = _swig_new_instance_method(_RF_Track.Volume_align_elements)
+    set_s0 = _swig_new_instance_method(_RF_Track.Volume_set_s0)
+    set_s1 = _swig_new_instance_method(_RF_Track.Volume_set_s1)
+    set_s0_from_s1 = _swig_new_instance_method(_RF_Track.Volume_set_s0_from_s1)
+    set_s1_from_s0 = _swig_new_instance_method(_RF_Track.Volume_set_s1_from_s0)
+    get_s0 = _swig_new_instance_method(_RF_Track.Volume_get_s0)
+    get_s1 = _swig_new_instance_method(_RF_Track.Volume_get_s1)
+    set_length = _swig_new_instance_method(_RF_Track.Volume_set_length)
+    get_static_Efield = _swig_new_instance_method(_RF_Track.Volume_get_static_Efield)
+    get_static_Bfield = _swig_new_instance_method(_RF_Track.Volume_get_static_Bfield)
+    set_static_Efield = _swig_new_instance_method(_RF_Track.Volume_set_static_Efield)
+    set_static_Bfield = _swig_new_instance_method(_RF_Track.Volume_set_static_Bfield)
+    set_t0 = _swig_new_instance_method(_RF_Track.Volume_set_t0)
+    unset_t0 = _swig_new_instance_method(_RF_Track.Volume_unset_t0)
+    autophase = _swig_new_instance_method(_RF_Track.Volume_autophase)
+    __getitem__ = _swig_new_instance_method(_RF_Track.Volume___getitem__)
+
+# Register Volume in _RF_Track:
+_RF_Track.Volume_swigregister(Volume)
+class Plasma(object):
+    r"""
+
+    Plasma properties used by the electron-cooling model.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    null_cell = property(_RF_Track.Plasma_null_cell_get, _RF_Track.Plasma_null_cell_set, doc=r"""null_cell : Plasma::Cell""")
+
+    def __init__(self, length=0.0, rx=0.0, ry=0.0, ne=0.0, Vz=0.0, D_ele_r=0.0, D_ele_l=0.0):
+        r"""
+        __init__(Plasma self, double length=0.0, double rx=0.0, double ry=0.0, double ne=0.0, double Vz=0.0, double D_ele_r=0.0, double D_ele_l=0.0) -> Plasma
+
+        Plasma properties used by the electron-cooling model.
+
+
+        """
+        _RF_Track.Plasma_swiginit(self, _RF_Track.new_Plasma(length, rx, ry, ne, Vz, D_ele_r, D_ele_l))
+    __swig_destroy__ = _RF_Track.delete_Plasma
+    get_self_fields = _swig_new_instance_method(_RF_Track.Plasma_get_self_fields)
+    get_rx = _swig_new_instance_method(_RF_Track.Plasma_get_rx)
+    get_ry = _swig_new_instance_method(_RF_Track.Plasma_get_ry)
+    get_area = _swig_new_instance_method(_RF_Track.Plasma_get_area)
+    set_radii = _swig_new_instance_method(_RF_Track.Plasma_set_radii)
+    is_point_inside_area = _swig_new_instance_method(_RF_Track.Plasma_is_point_inside_area)
+    set_Q = _swig_new_instance_method(_RF_Track.Plasma_set_Q)
+    set_mass = _swig_new_instance_method(_RF_Track.Plasma_set_mass)
+    set_length = _swig_new_instance_method(_RF_Track.Plasma_set_length)
+    set_temperature = _swig_new_instance_method(_RF_Track.Plasma_set_temperature)
+    set_plasma_mesh = _swig_new_instance_method(_RF_Track.Plasma_set_plasma_mesh)
+    set_nsteps = _swig_new_instance_method(_RF_Track.Plasma_set_nsteps)
+    get_Q = _swig_new_instance_method(_RF_Track.Plasma_get_Q)
+    get_mass = _swig_new_instance_method(_RF_Track.Plasma_get_mass)
+    get_length = _swig_new_instance_method(_RF_Track.Plasma_get_length)
+    get_density = _swig_new_instance_method(_RF_Track.Plasma_get_density)
+    get_D_ele_r = _swig_new_instance_method(_RF_Track.Plasma_get_D_ele_r)
+    get_D_ele_l = _swig_new_instance_method(_RF_Track.Plasma_get_D_ele_l)
+    get_temperature = _swig_new_instance_method(_RF_Track.Plasma_get_temperature)
+    get_debye_length = _swig_new_instance_method(_RF_Track.Plasma_get_debye_length)
+    get_plasma_parameter = _swig_new_instance_method(_RF_Track.Plasma_get_plasma_parameter)
+    get_velocity = _swig_new_instance_method(_RF_Track.Plasma_get_velocity)
+    get_current_density = _swig_new_instance_method(_RF_Track.Plasma_get_current_density)
+    get_density_mesh = _swig_new_instance_method(_RF_Track.Plasma_get_density_mesh)
+    get_static_Bfield = _swig_new_instance_method(_RF_Track.Plasma_get_static_Bfield)
+    set_static_Bfield = _swig_new_instance_method(_RF_Track.Plasma_set_static_Bfield)
+    project = _swig_new_instance_method(_RF_Track.Plasma_project)
+    advect = _swig_new_instance_method(_RF_Track.Plasma_advect)
+    apply_momentum_through_dt = _swig_new_instance_method(_RF_Track.Plasma_apply_momentum_through_dt)
+    apply_momentum_through_dS = _swig_new_instance_method(_RF_Track.Plasma_apply_momentum_through_dS)
+
+# Register Plasma in _RF_Track:
+_RF_Track.Plasma_swigregister(Plasma)
+class ElectronCooler(Element, Plasma):
+    r"""
+
+    Electron-cooling element. Length and transverse radii are in m, electron density in
+    1/m^3 and longitudinal velocity in units of c.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length=0.0, rx=0.0, ry=0.0, ne=0.0, Vz=0.0):
+        r"""
+        __init__(ElectronCooler self, double length=0.0, double rx=0.0, double ry=0.0, double ne=0.0, double Vz=0.0) -> ElectronCooler
+
+        Electron-cooling element. Length and transverse radii are in m, electron density in
+        1/m^3 and longitudinal velocity in units of c.
+
+
+        """
+        _RF_Track.ElectronCooler_swiginit(self, _RF_Track.new_ElectronCooler(length, rx, ry, ne, Vz))
+    get_current = _swig_new_instance_method(_RF_Track.ElectronCooler_get_current)
+    set_electron_mesh = _swig_new_instance_method(_RF_Track.ElectronCooler_set_electron_mesh)
+    compute_cooling_force = _swig_new_instance_method(_RF_Track.ElectronCooler_compute_cooling_force)
+    __swig_destroy__ = _RF_Track.delete_ElectronCooler
+
+# Register ElectronCooler in _RF_Track:
+_RF_Track.ElectronCooler_swigregister(ElectronCooler)
+class SBend(Element, ):
+    r"""
+
+    A sector bending dipole. L is in m, angle and pole-face angles E1/E2 in rad, reference
+    P_over_Q in MV/c, and K1 in 1/m^2. Supports edge focusing and a quadrupole component.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SBend self, double L, double angle, double ref_P_over_Q, double E1_=0.0, double E2_=0.0, double K1_=0.0) -> SBend
+        __init__(SBend self, double L=0.0) -> SBend
+
+        A sector bending dipole. L is in m, angle and pole-face angles E1/E2 in rad, reference
+        P_over_Q in MV/c, and K1 in 1/m^2. Supports edge focusing and a quadrupole component.
+
+
+        """
+        _RF_Track.SBend_swiginit(self, _RF_Track.new_SBend(*args))
+    get_E1 = _swig_new_instance_method(_RF_Track.SBend_get_E1)
+    get_E2 = _swig_new_instance_method(_RF_Track.SBend_get_E2)
+    get_E1d = _swig_new_instance_method(_RF_Track.SBend_get_E1d)
+    get_E2d = _swig_new_instance_method(_RF_Track.SBend_get_E2d)
+    get_h = _swig_new_instance_method(_RF_Track.SBend_get_h)
+    get_K0 = _swig_new_instance_method(_RF_Track.SBend_get_K0)
+    get_K1 = _swig_new_instance_method(_RF_Track.SBend_get_K1)
+    get_K1L = _swig_new_instance_method(_RF_Track.SBend_get_K1L)
+    get_gradient = _swig_new_instance_method(_RF_Track.SBend_get_gradient)
+    get_hgap = _swig_new_instance_method(_RF_Track.SBend_get_hgap)
+    get_rho = _swig_new_instance_method(_RF_Track.SBend_get_rho)
+    get_angle = _swig_new_instance_method(_RF_Track.SBend_get_angle)
+    get_angled = _swig_new_instance_method(_RF_Track.SBend_get_angled)
+    get_Brho = _swig_new_instance_method(_RF_Track.SBend_get_Brho)
+    get_Bfield = _swig_new_instance_method(_RF_Track.SBend_get_Bfield)
+    get_P_over_Q = _swig_new_instance_method(_RF_Track.SBend_get_P_over_Q)
+    set_E1 = _swig_new_instance_method(_RF_Track.SBend_set_E1)
+    set_E2 = _swig_new_instance_method(_RF_Track.SBend_set_E2)
+    set_E1d = _swig_new_instance_method(_RF_Track.SBend_set_E1d)
+    set_E2d = _swig_new_instance_method(_RF_Track.SBend_set_E2d)
+    set_h = _swig_new_instance_method(_RF_Track.SBend_set_h)
+    set_hgap = _swig_new_instance_method(_RF_Track.SBend_set_hgap)
+    set_fint = _swig_new_instance_method(_RF_Track.SBend_set_fint)
+    set_K0 = _swig_new_instance_method(_RF_Track.SBend_set_K0)
+    set_K1 = _swig_new_instance_method(_RF_Track.SBend_set_K1)
+    set_K1L = _swig_new_instance_method(_RF_Track.SBend_set_K1L)
+    set_gradient = _swig_new_instance_method(_RF_Track.SBend_set_gradient)
+    set_angle = _swig_new_instance_method(_RF_Track.SBend_set_angle)
+    set_angled = _swig_new_instance_method(_RF_Track.SBend_set_angled)
+    set_length = _swig_new_instance_method(_RF_Track.SBend_set_length)
+    set_Brho = _swig_new_instance_method(_RF_Track.SBend_set_Brho)
+    set_Bfield = _swig_new_instance_method(_RF_Track.SBend_set_Bfield)
+    set_P_over_Q = _swig_new_instance_method(_RF_Track.SBend_set_P_over_Q)
+    get_s0 = _swig_new_instance_method(_RF_Track.SBend_get_s0)
+    get_s1 = _swig_new_instance_method(_RF_Track.SBend_get_s1)
+    __swig_destroy__ = _RF_Track.delete_SBend
+
+# Register SBend in _RF_Track:
+_RF_Track.SBend_swigregister(SBend)
+RBend = _RF_Track.RBend
+class TransferLine(Element):
+    r"""
+
+    A transfer-map element constructed from Twiss data or a MAD-X Twiss file. Reference
+    momentum is in MeV/c. Supports chromaticity and momentum-compaction parameters.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(TransferLine self, MatrixNd const & twiss_matrix_, double Pref, double DQx_, double DQy_, double momentum_compaction_) -> TransferLine
+        __init__(TransferLine self, MatrixNd const & twiss_matrix_, double Pref) -> TransferLine
+        __init__(TransferLine self, std::string const & twiss_file, double Pref=-1.0) -> TransferLine
+        __init__(TransferLine self, TransferLine t) -> TransferLine
+        __init__(TransferLine self) -> TransferLine
+
+        A transfer-map element constructed from Twiss data or a MAD-X Twiss file. Reference
+        momentum is in MeV/c. Supports chromaticity and momentum-compaction parameters.
+
+
+        """
+        _RF_Track.TransferLine_swiginit(self, _RF_Track.new_TransferLine(*args))
+    __swig_destroy__ = _RF_Track.delete_TransferLine
+    __assign__ = _swig_new_instance_method(_RF_Track.TransferLine___assign__)
+
+# Register TransferLine in _RF_Track:
+_RF_Track.TransferLine_swigregister(TransferLine)
+class Quadrupole(Element):
+    r"""
+
+    A matrix-based quadrupole. Constructors take (L [m], integrated strength [MV/c/m]) or (L
+    [m], P_over_Q [MV/c], k1 [1/m^2]). Strength equals P_over_Q*k1*L. Use set_gradient for a
+    gradient in T/m.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Quadrupole self, double L=0.0, double strength_=0.0) -> Quadrupole
+        __init__(Quadrupole self, double L, double P_over_Q, double k1_) -> Quadrupole
+
+        A matrix-based quadrupole. Constructors take (L [m], integrated strength [MV/c/m]) or (L
+        [m], P_over_Q [MV/c], k1 [1/m^2]). Strength equals P_over_Q*k1*L. Use set_gradient for a
+        gradient in T/m.
+
+
+        """
+        _RF_Track.Quadrupole_swiginit(self, _RF_Track.new_Quadrupole(*args))
+    get_K1 = _swig_new_instance_method(_RF_Track.Quadrupole_get_K1)
+    get_K1L = _swig_new_instance_method(_RF_Track.Quadrupole_get_K1L)
+    get_strength = _swig_new_instance_method(_RF_Track.Quadrupole_get_strength)
+    get_gradient = _swig_new_instance_method(_RF_Track.Quadrupole_get_gradient)
+    set_K1 = _swig_new_instance_method(_RF_Track.Quadrupole_set_K1)
+    set_K1L = _swig_new_instance_method(_RF_Track.Quadrupole_set_K1L)
+    set_P_over_Q = _swig_new_instance_method(_RF_Track.Quadrupole_set_P_over_Q)
+    set_length = _swig_new_instance_method(_RF_Track.Quadrupole_set_length)
+    set_strength = _swig_new_instance_method(_RF_Track.Quadrupole_set_strength)
+    set_gradient = _swig_new_instance_method(_RF_Track.Quadrupole_set_gradient)
+    __swig_destroy__ = _RF_Track.delete_Quadrupole
+
+# Register Quadrupole in _RF_Track:
+_RF_Track.Quadrupole_swigregister(Quadrupole)
+class GenericField(Element, Parallel_ODE_Solver):
+    r"""
+
+    Base element for numerical tracking through an electromagnetic field. Length is in
+    metres; inherited ODE-solver controls select the integration algorithm.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined")
+    __repr__ = _swig_repr
+    set_length = _swig_new_instance_method(_RF_Track.GenericField_set_length)
+    set_scaling_factor = _swig_new_instance_method(_RF_Track.GenericField_set_scaling_factor)
+    get_scaling_factor = _swig_new_instance_method(_RF_Track.GenericField_get_scaling_factor)
+    get_static_Efield = _swig_new_instance_method(_RF_Track.GenericField_get_static_Efield)
+    get_static_Bfield = _swig_new_instance_method(_RF_Track.GenericField_get_static_Bfield)
+    set_static_Efield = _swig_new_instance_method(_RF_Track.GenericField_set_static_Efield)
+    set_static_Bfield = _swig_new_instance_method(_RF_Track.GenericField_set_static_Bfield)
+    get_field = _swig_new_instance_method(_RF_Track.GenericField_get_field)
+    enable_end_fields = _swig_new_instance_method(_RF_Track.GenericField_enable_end_fields)
+    disable_end_fields = _swig_new_instance_method(_RF_Track.GenericField_disable_end_fields)
+    __swig_destroy__ = _RF_Track.delete_GenericField
+
+# Register GenericField in _RF_Track:
+_RF_Track.GenericField_swigregister(GenericField)
+class Solenoid(Element, ):
+    r"""
+
+    A solenoid magnetic field. Length and radii are in metres; B0 is in T. A
+    finite-thickness winding can be represented by concentric current sheets.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Solenoid self) -> Solenoid
+        __init__(Solenoid self, double L, double B0_=0.0, double R_=0.0) -> Solenoid
+        __init__(Solenoid self, double L, double B0_, double inner_R, double outer_R, int n_sheets) -> Solenoid
+
+        A solenoid magnetic field. Length and radii are in metres; B0 is in T. A
+        finite-thickness winding can be represented by concentric current sheets.
+
+
+        """
+        _RF_Track.Solenoid_swiginit(self, _RF_Track.new_Solenoid(*args))
+    get_KS = _swig_new_instance_method(_RF_Track.Solenoid_get_KS)
+    get_inner_R = _swig_new_instance_method(_RF_Track.Solenoid_get_inner_R)
+    get_outer_R = _swig_new_instance_method(_RF_Track.Solenoid_get_outer_R)
+    get_max_current = _swig_new_instance_method(_RF_Track.Solenoid_get_max_current)
+    get_current_density = _swig_new_instance_method(_RF_Track.Solenoid_get_current_density)
+    set_KS = _swig_new_instance_method(_RF_Track.Solenoid_set_KS)
+    set_R = _swig_new_instance_method(_RF_Track.Solenoid_set_R)
+    set_length = _swig_new_instance_method(_RF_Track.Solenoid_set_length)
+    set_max_current = _swig_new_instance_method(_RF_Track.Solenoid_set_max_current)
+    set_current_density = _swig_new_instance_method(_RF_Track.Solenoid_set_current_density)
+    __swig_destroy__ = _RF_Track.delete_Solenoid
+
+# Register Solenoid in _RF_Track:
+_RF_Track.Solenoid_swigregister(Solenoid)
+class Coil(GenericField, ):
+    r"""
+
+    Magnetic field of a circular current coil. Geometry parameters are in metres and the
+    magnetic-field amplitude is in T.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Coil self) -> Coil
+        __init__(Coil self, double length, double B0_, double R_) -> Coil
+
+        Magnetic field of a circular current coil. Geometry parameters are in metres and the
+        magnetic-field amplitude is in T.
+
+
+        """
+        _RF_Track.Coil_swiginit(self, _RF_Track.new_Coil(*args))
+    get_R = _swig_new_instance_method(_RF_Track.Coil_get_R)
+    get_B0 = _swig_new_instance_method(_RF_Track.Coil_get_B0)
+    set_R = _swig_new_instance_method(_RF_Track.Coil_set_R)
+    set_B0 = _swig_new_instance_method(_RF_Track.Coil_set_B0)
+    set_R_and_current = _swig_new_instance_method(_RF_Track.Coil_set_R_and_current)
+    __swig_destroy__ = _RF_Track.delete_Coil
+
+# Register Coil in _RF_Track:
+_RF_Track.Coil_swigregister(Coil)
+class SpaceCharge_Field(GenericField, ):
+    r"""
+
+    An electromagnetic field reconstructed from a Bunch6dT distribution on a mesh with Nx,
+    Ny, Nz cells. Vz_slices controls velocity slicing.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SpaceCharge_Field self) -> SpaceCharge_Field
+        __init__(SpaceCharge_Field self, Bunch6dT bunch, size_t Nx, size_t Ny, size_t Nz, size_t Vz_slices=1) -> SpaceCharge_Field
+
+        An electromagnetic field reconstructed from a Bunch6dT distribution on a mesh with Nx,
+        Ny, Nz cells. Vz_slices controls velocity slicing.
+
+
+        """
+        _RF_Track.SpaceCharge_Field_swiginit(self, _RF_Track.new_SpaceCharge_Field(*args))
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_Field
+    get_Q = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_get_Q)
+    get_V = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_get_V)
+
+# Register SpaceCharge_Field in _RF_Track:
+_RF_Track.SpaceCharge_Field_swigregister(SpaceCharge_Field)
+class SpaceCharge_Field_LongCylinder(GenericField, ):
+    r"""
+
+    Space-charge field reconstructed inside a long conducting cylinder. The explicit
+    aperture argument is in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SpaceCharge_Field_LongCylinder self) -> SpaceCharge_Field_LongCylinder
+        __init__(SpaceCharge_Field_LongCylinder self, Bunch6dT bunch, size_t Nx, size_t Ny, size_t Nz, size_t Vz_slices=1) -> SpaceCharge_Field_LongCylinder
+        __init__(SpaceCharge_Field_LongCylinder self, Bunch6dT bunch, size_t Nx, size_t Ny, size_t Nz, size_t Vz_slices, double aperture) -> SpaceCharge_Field_LongCylinder
+
+        Space-charge field reconstructed inside a long conducting cylinder. The explicit
+        aperture argument is in metres.
+
+
+        """
+        _RF_Track.SpaceCharge_Field_LongCylinder_swiginit(self, _RF_Track.new_SpaceCharge_Field_LongCylinder(*args))
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_Field_LongCylinder
+    get_Q = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_LongCylinder_get_Q)
+    get_V = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_LongCylinder_get_V)
+    get_aperture = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_LongCylinder_get_aperture)
+
+# Register SpaceCharge_Field_LongCylinder in _RF_Track:
+_RF_Track.SpaceCharge_Field_LongCylinder_swigregister(SpaceCharge_Field_LongCylinder)
+class SpaceCharge_Field_HorizontalPlates(GenericField, ):
+    r"""
+
+    Space-charge field reconstructed between conducting horizontal plates. The explicit
+    half-gap argument is in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SpaceCharge_Field_HorizontalPlates self) -> SpaceCharge_Field_HorizontalPlates
+        __init__(SpaceCharge_Field_HorizontalPlates self, Bunch6dT bunch, size_t Nx, size_t Ny, size_t Nz, size_t Vz_slices=1) -> SpaceCharge_Field_HorizontalPlates
+        __init__(SpaceCharge_Field_HorizontalPlates self, Bunch6dT bunch, size_t Nx, size_t Ny, size_t Nz, size_t Vz_slices, double hgap) -> SpaceCharge_Field_HorizontalPlates
+
+        Space-charge field reconstructed between conducting horizontal plates. The explicit
+        half-gap argument is in metres.
+
+
+        """
+        _RF_Track.SpaceCharge_Field_HorizontalPlates_swiginit(self, _RF_Track.new_SpaceCharge_Field_HorizontalPlates(*args))
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_Field_HorizontalPlates
+    get_Q = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_HorizontalPlates_get_Q)
+    get_V = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_HorizontalPlates_get_V)
+    get_hgap = _swig_new_instance_method(_RF_Track.SpaceCharge_Field_HorizontalPlates_get_hgap)
+
+# Register SpaceCharge_Field_HorizontalPlates in _RF_Track:
+_RF_Track.SpaceCharge_Field_HorizontalPlates_swigregister(SpaceCharge_Field_HorizontalPlates)
+class ToroidalHarmonics(GenericField, ):
+    r"""
+
+    Magnetic-field element represented through toroidal harmonics.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    M00 = property(_RF_Track.ToroidalHarmonics_M00_get, _RF_Track.ToroidalHarmonics_M00_set, doc=r"""M00 : double""")
+    Mcc = property(_RF_Track.ToroidalHarmonics_Mcc_get, _RF_Track.ToroidalHarmonics_Mcc_set, doc=r"""Mcc : MatrixNd""")
+    Mcs = property(_RF_Track.ToroidalHarmonics_Mcs_get, _RF_Track.ToroidalHarmonics_Mcs_set, doc=r"""Mcs : MatrixNd""")
+    Msc = property(_RF_Track.ToroidalHarmonics_Msc_get, _RF_Track.ToroidalHarmonics_Msc_set, doc=r"""Msc : MatrixNd""")
+    Mss = property(_RF_Track.ToroidalHarmonics_Mss_get, _RF_Track.ToroidalHarmonics_Mss_set, doc=r"""Mss : MatrixNd""")
+
+    def __init__(self, Rmajor, Rminor):
+        r"""
+        __init__(ToroidalHarmonics self, double Rmajor, double Rminor) -> ToroidalHarmonics
+
+        Magnetic-field element represented through toroidal harmonics.
+
+
+        """
+        _RF_Track.ToroidalHarmonics_swiginit(self, _RF_Track.new_ToroidalHarmonics(Rmajor, Rminor))
+    get_a = _swig_new_instance_method(_RF_Track.ToroidalHarmonics_get_a)
+    set_Mcc = _swig_new_instance_method(_RF_Track.ToroidalHarmonics_set_Mcc)
+    set_Mcs = _swig_new_instance_method(_RF_Track.ToroidalHarmonics_set_Mcs)
+    set_Msc = _swig_new_instance_method(_RF_Track.ToroidalHarmonics_set_Msc)
+    set_Mss = _swig_new_instance_method(_RF_Track.ToroidalHarmonics_set_Mss)
+    __swig_destroy__ = _RF_Track.delete_ToroidalHarmonics
+
+# Register ToroidalHarmonics in _RF_Track:
+_RF_Track.ToroidalHarmonics_swigregister(ToroidalHarmonics)
+class ConstantField(GenericField):
+    r"""
+
+    A region of constant electric and magnetic fields. Length is in metres, electric field
+    in V/m and magnetic field in T.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length=0.0):
+        r"""
+        __init__(ConstantField self, double length=0.0) -> ConstantField
+
+        A region of constant electric and magnetic fields. Length is in metres, electric field
+        in V/m and magnetic field in T.
+
+
+        """
+        _RF_Track.ConstantField_swiginit(self, _RF_Track.new_ConstantField(length))
+    __swig_destroy__ = _RF_Track.delete_ConstantField
+
+# Register ConstantField in _RF_Track:
+_RF_Track.ConstantField_swigregister(ConstantField)
+class AdiabaticMatchingDevice(GenericField):
+    r"""
+
+    An axially varying solenoidal field for adiabatic beam matching. Geometry is specified in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length_=0.0, Bmax_=0.0, fMu_=0.0):
+        r"""
+        __init__(AdiabaticMatchingDevice self, double length_=0.0, double Bmax_=0.0, double fMu_=0.0) -> AdiabaticMatchingDevice
+
+        An axially varying solenoidal field for adiabatic beam matching. Geometry is specified in metres.
+
+
+        """
+        _RF_Track.AdiabaticMatchingDevice_swiginit(self, _RF_Track.new_AdiabaticMatchingDevice(length_, Bmax_, fMu_))
+    get_mu = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_get_mu)
+    get_B0 = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_get_B0)
+    set_mu = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_set_mu)
+    set_B0 = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_set_B0)
+    get_entrance_aperture = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_get_entrance_aperture)
+    get_exit_aperture = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_get_exit_aperture)
+    set_entrance_aperture = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_set_entrance_aperture)
+    set_exit_aperture = _swig_new_instance_method(_RF_Track.AdiabaticMatchingDevice_set_exit_aperture)
+    __swig_destroy__ = _RF_Track.delete_AdiabaticMatchingDevice
+
+# Register AdiabaticMatchingDevice in _RF_Track:
+_RF_Track.AdiabaticMatchingDevice_swigregister(AdiabaticMatchingDevice)
+class BroadbandDielectricStructure(GenericField):
+    r"""
+
+    A dielectric structure represented by a broadband electromagnetic field.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, x0, y0, hx, hy, hz, length, t0, t1, N, P_map=1.0, P_actual=1.0):
+        r"""
+        __init__(BroadbandDielectricStructure self, double x0, double y0, double hx, double hy, double hz, double length, double t0, double t1, size_t N, double P_map=1.0, double P_actual=1.0) -> BroadbandDielectricStructure
+
+        A dielectric structure represented by a broadband electromagnetic field.
+
+
+        """
+        _RF_Track.BroadbandDielectricStructure_swiginit(self, _RF_Track.new_BroadbandDielectricStructure(x0, y0, hx, hy, hz, length, t0, t1, N, P_map, P_actual))
+    resize = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_resize)
+    set_field_n = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_field_n)
+    get_x0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_x0)
+    get_y0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_y0)
+    get_z0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_z1)
+    get_hx = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_hz)
+    get_t0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_t0)
+    get_t1 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_t1)
+    get_dt = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_dt)
+    get_P_map = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_get_P_actual)
+    set_x0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_x0)
+    set_y0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_y0)
+    set_z0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_z1)
+    set_hx = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_hx)
+    set_hy = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_hy)
+    set_hz = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_hz)
+    set_t0 = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_t0)
+    set_dt = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_dt)
+    set_P_map = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_P_map)
+    set_P_actual = _swig_new_instance_method(_RF_Track.BroadbandDielectricStructure_set_P_actual)
+    __swig_destroy__ = _RF_Track.delete_BroadbandDielectricStructure
+
+# Register BroadbandDielectricStructure in _RF_Track:
+_RF_Track.BroadbandDielectricStructure_swigregister(BroadbandDielectricStructure)
+class Multipole(GenericField):
+    r"""
+
+    A magnetic multipole element with complex normal and skew coefficients. Length is in
+    metres. Supports field coefficients Bn or integrated normalized coefficients KnL at a
+    specified P_over_Q.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Multipole self, double length_=0.0, std::vector< std::complex< double > > const & Bn_=std::vector< std::complex< double > >(size_t(0))) -> Multipole
+        __init__(Multipole self, double length_, double P_over_Q, std::vector< std::complex< double > > const & KnL) -> Multipole
+
+        A magnetic multipole element with complex normal and skew coefficients. Length is in
+        metres. Supports field coefficients Bn or integrated normalized coefficients KnL at a
+        specified P_over_Q.
+
+
+        """
+        _RF_Track.Multipole_swiginit(self, _RF_Track.new_Multipole(*args))
+    get_Bn = _swig_new_instance_method(_RF_Track.Multipole_get_Bn)
+    get_KnL = _swig_new_instance_method(_RF_Track.Multipole_get_KnL)
+    get_strengths = _swig_new_instance_method(_RF_Track.Multipole_get_strengths)
+    set_Bn = _swig_new_instance_method(_RF_Track.Multipole_set_Bn)
+    set_KnL = _swig_new_instance_method(_RF_Track.Multipole_set_KnL)
+    set_P_over_Q = _swig_new_instance_method(_RF_Track.Multipole_set_P_over_Q)
+    set_strengths = _swig_new_instance_method(_RF_Track.Multipole_set_strengths)
+    __swig_destroy__ = _RF_Track.delete_Multipole
+
+# Register Multipole in _RF_Track:
+_RF_Track.Multipole_swigregister(Multipole)
+class TimeDependent_Field(GenericField):
+    r"""
+
+    Base interface for elements with time-dependent fields and an adjustable synchronization time t0 in mm/c.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length=0.0):
+        r"""
+        __init__(TimeDependent_Field self, double length=0.0) -> TimeDependent_Field
+
+        Base interface for elements with time-dependent fields and an adjustable synchronization time t0 in mm/c.
+
+
+        """
+        _RF_Track.TimeDependent_Field_swiginit(self, _RF_Track.new_TimeDependent_Field(length))
+    get_t0 = _swig_new_instance_method(_RF_Track.TimeDependent_Field_get_t0)
+    set_t0 = _swig_new_instance_method(_RF_Track.TimeDependent_Field_set_t0)
+    unset_t0 = _swig_new_instance_method(_RF_Track.TimeDependent_Field_unset_t0)
+    synchronize_with = _swig_new_instance_method(_RF_Track.TimeDependent_Field_synchronize_with)
+    t0_is_synchronized = _swig_new_instance_method(_RF_Track.TimeDependent_Field_t0_is_synchronized)
+    t0_is_set = _swig_new_instance_method(_RF_Track.TimeDependent_Field_t0_is_set)
+    t0_is_unset = _swig_new_instance_method(_RF_Track.TimeDependent_Field_t0_is_unset)
+    __swig_destroy__ = _RF_Track.delete_TimeDependent_Field
+
+# Register TimeDependent_Field in _RF_Track:
+_RF_Track.TimeDependent_Field_swigregister(TimeDependent_Field)
+class Sextupole(GenericField):
+    r"""
+
+    A sextupole magnet. Length is in m; the rigidity-based constructor takes P_over_Q in
+    MV/c and k2 in 1/m^3. Integrated strength equals P_over_Q*k2*L.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Sextupole self, double L=0.0, double strength_=0.0) -> Sextupole
+        __init__(Sextupole self, double L, double P_over_Q, double k2) -> Sextupole
+
+        A sextupole magnet. Length is in m; the rigidity-based constructor takes P_over_Q in
+        MV/c and k2 in 1/m^3. Integrated strength equals P_over_Q*k2*L.
+
+
+        """
+        _RF_Track.Sextupole_swiginit(self, _RF_Track.new_Sextupole(*args))
+    get_K2 = _swig_new_instance_method(_RF_Track.Sextupole_get_K2)
+    get_K2L = _swig_new_instance_method(_RF_Track.Sextupole_get_K2L)
+    get_strength = _swig_new_instance_method(_RF_Track.Sextupole_get_strength)
+    get_gradient = _swig_new_instance_method(_RF_Track.Sextupole_get_gradient)
+    set_K2 = _swig_new_instance_method(_RF_Track.Sextupole_set_K2)
+    set_K2L = _swig_new_instance_method(_RF_Track.Sextupole_set_K2L)
+    set_P_over_Q = _swig_new_instance_method(_RF_Track.Sextupole_set_P_over_Q)
+    set_strength = _swig_new_instance_method(_RF_Track.Sextupole_set_strength)
+    set_gradient = _swig_new_instance_method(_RF_Track.Sextupole_set_gradient)
+    __swig_destroy__ = _RF_Track.delete_Sextupole
+
+# Register Sextupole in _RF_Track:
+_RF_Track.Sextupole_swigregister(Sextupole)
+class Aperture_1d(GenericField):
+    r"""
+
+    A longitudinally varying aperture supplied on a one-dimensional mesh. Aperture sizes and
+    longitudinal mesh spacing are in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Aperture_1d self, Mesh1d const & aperture_1d, double hz, double length=-1) -> Aperture_1d
+        __init__(Aperture_1d self) -> Aperture_1d
+
+        A longitudinally varying aperture supplied on a one-dimensional mesh. Aperture sizes and
+        longitudinal mesh spacing are in metres.
+
+
+        """
+        _RF_Track.Aperture_1d_swiginit(self, _RF_Track.new_Aperture_1d(*args))
+    get_hz = _swig_new_instance_method(_RF_Track.Aperture_1d_get_hz)
+    get_nz = _swig_new_instance_method(_RF_Track.Aperture_1d_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.Aperture_1d_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Aperture_1d_get_z1)
+    set_hz = _swig_new_instance_method(_RF_Track.Aperture_1d_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.Aperture_1d_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Aperture_1d_set_z1)
+    __swig_destroy__ = _RF_Track.delete_Aperture_1d
+
+# Register Aperture_1d in _RF_Track:
+_RF_Track.Aperture_1d_swigregister(Aperture_1d)
+class Static_Electric_FieldMap(GenericField):
+    r"""
+
+    Static electric field on a regular three-dimensional mesh, provided as potential [V] or
+    electric-field components [V/m]. Geometry arguments are in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Electric_FieldMap self, Mesh3d const & Phi=DefaultMesh3d, double x0=0.0, double y0=0.0, double hx=1.0, double hy=1.0, double hz=1.0, double length=-1.0) -> Static_Electric_FieldMap
+        __init__(Static_Electric_FieldMap self, Mesh3d const & Ex, Mesh3d const & Ey, Mesh3d const & Ez, double x0, double y0, double hx, double hy, double hz, double length=-1.0, double quality=1.0) -> Static_Electric_FieldMap
+
+        Static electric field on a regular three-dimensional mesh, provided as potential [V] or
+        electric-field components [V/m]. Geometry arguments are in metres.
+
+
+        """
+        _RF_Track.Static_Electric_FieldMap_swiginit(self, _RF_Track.new_Static_Electric_FieldMap(*args))
+    __swig_destroy__ = _RF_Track.delete_Static_Electric_FieldMap
+    get_hx = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_hz)
+    get_nx = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_nz)
+    get_x0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_x0)
+    get_y0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_y0)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_z0)
+    get_x1 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_x1)
+    get_y1 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_y1)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_z1)
+    set_hx = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_hx)
+    set_hy = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_hy)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_hz)
+    set_x0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_x0)
+    set_y0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_y0)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_z1)
+    set_Ex_Ey_Ez = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_set_Ex_Ey_Ez)
+    get_Phi = _swig_new_instance_method(_RF_Track.Static_Electric_FieldMap_get_Phi)
+
+# Register Static_Electric_FieldMap in _RF_Track:
+_RF_Track.Static_Electric_FieldMap_swigregister(Static_Electric_FieldMap)
+class Static_Magnetic_FieldMap(GenericField):
+    r"""
+
+    Static magnetic field on a regular three-dimensional mesh, provided as vector potential
+    [T.m] or magnetic-field components [T]. Geometry arguments are in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Magnetic_FieldMap self, Mesh3d const & Ax=DefaultMesh3d, Mesh3d const & Ay=DefaultMesh3d, Mesh3d const & Az=DefaultMesh3d, Mesh3d const & PhiM=DefaultMesh3d, double x0=0.0, double y0=0.0, double hx=1.0, double hy=1.0, double hz=1.0, double length=-1.0) -> Static_Magnetic_FieldMap
+        __init__(Static_Magnetic_FieldMap self, Mesh3d const & Bx, Mesh3d const & By, Mesh3d const & Bz, double x0, double y0, double hx, double hy, double hz, double length=-1.0, double quality=1.0) -> Static_Magnetic_FieldMap
+
+        Static magnetic field on a regular three-dimensional mesh, provided as vector potential
+        [T.m] or magnetic-field components [T]. Geometry arguments are in metres.
+
+
+        """
+        _RF_Track.Static_Magnetic_FieldMap_swiginit(self, _RF_Track.new_Static_Magnetic_FieldMap(*args))
+    __swig_destroy__ = _RF_Track.delete_Static_Magnetic_FieldMap
+    get_hx = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_hz)
+    get_nx = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_nz)
+    get_x0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_x0)
+    get_y0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_y0)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_z0)
+    get_x1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_x1)
+    get_y1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_y1)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_z1)
+    set_hx = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_hx)
+    set_hy = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_hy)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_hz)
+    set_x0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_x0)
+    set_y0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_y0)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_z1)
+    set_Bx_By_Bz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_set_Bx_By_Bz)
+    get_Ax = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_Ax)
+    get_Ay = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_Ay)
+    get_Az = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_Az)
+    get_PhiM = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_PhiM)
+    get_divB = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_get_divB)
+
+# Register Static_Magnetic_FieldMap in _RF_Track:
+_RF_Track.Static_Magnetic_FieldMap_swigregister(Static_Magnetic_FieldMap)
+class Undulator(GenericField):
+    r"""Proxy of C++ Undulator class."""
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, lperiod, K, nperiods, kx2=0.0):
+        r"""__init__(Undulator self, double lperiod, double K, double nperiods, double kx2=0.0) -> Undulator"""
+        _RF_Track.Undulator_swiginit(self, _RF_Track.new_Undulator(lperiod, K, nperiods, kx2))
+    get_B = _swig_new_instance_method(_RF_Track.Undulator_get_B)
+    set_B = _swig_new_instance_method(_RF_Track.Undulator_set_B)
+    __swig_destroy__ = _RF_Track.delete_Undulator
+
+# Register Undulator in _RF_Track:
+_RF_Track.Undulator_swigregister(Undulator)
+class Absorber(ConstantField):
+    r"""
+
+    Material region with scattering, energy loss and straggling. Construct from length [m]
+    and material name, or radiation length [cm], atomic Z, A [g/mol], density [g/cm^3] and
+    mean excitation energy [eV].
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Absorber self, double length=0.0, double radiation_length=0.0, double Z=0.0, double A=0.0, double density=0.0, double mean_excitation_energy=-1.0) -> Absorber
+        __init__(Absorber self, double length, std::string const & material_name) -> Absorber
+
+        Material region with scattering, energy loss and straggling. Construct from length [m]
+        and material name, or radiation length [cm], atomic Z, A [g/mol], density [g/cm^3] and
+        mean excitation energy [eV].
+
+
+        """
+        _RF_Track.Absorber_swiginit(self, _RF_Track.new_Absorber(*args))
+    enable_log_term = _swig_new_instance_method(_RF_Track.Absorber_enable_log_term)
+    disable_log_term = _swig_new_instance_method(_RF_Track.Absorber_disable_log_term)
+    set_shape = _swig_new_instance_method(_RF_Track.Absorber_set_shape)
+    get_shape = _swig_new_instance_method(_RF_Track.Absorber_get_shape)
+    enable_multiple_coulomb_scattering = _swig_new_instance_method(_RF_Track.Absorber_enable_multiple_coulomb_scattering)
+    enable_energy_straggling = _swig_new_instance_method(_RF_Track.Absorber_enable_energy_straggling)
+    enable_stopping_power = _swig_new_instance_method(_RF_Track.Absorber_enable_stopping_power)
+    enable_fruehwirth_model = _swig_new_instance_method(_RF_Track.Absorber_enable_fruehwirth_model)
+    disable_fruehwirth_model = _swig_new_instance_method(_RF_Track.Absorber_disable_fruehwirth_model)
+    enable_wentzel_model = _swig_new_instance_method(_RF_Track.Absorber_enable_wentzel_model)
+    disable_wentzel_model = _swig_new_instance_method(_RF_Track.Absorber_disable_wentzel_model)
+    disable_multiple_coulomb_scattering = _swig_new_instance_method(_RF_Track.Absorber_disable_multiple_coulomb_scattering)
+    disable_energy_straggling = _swig_new_instance_method(_RF_Track.Absorber_disable_energy_straggling)
+    disable_stopping_power = _swig_new_instance_method(_RF_Track.Absorber_disable_stopping_power)
+    __swig_destroy__ = _RF_Track.delete_Absorber
+
+# Register Absorber in _RF_Track:
+_RF_Track.Absorber_swigregister(Absorber)
+class Bpm(ConstantField):
+    r"""
+
+    Beam position monitor. Constructor length is in m, spatial resolution in mm and time
+    resolution in mm/c. Readings are collected during tracking.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length=0.0, bpm_resolution=0.0, bpm_resolution_t=0.0):
+        r"""
+        __init__(Bpm self, double length=0.0, double const bpm_resolution=0.0, double const bpm_resolution_t=0.0) -> Bpm
+
+        Beam position monitor. Constructor length is in m, spatial resolution in mm and time
+        resolution in mm/c. Readings are collected during tracking.
+
+
+        """
+        _RF_Track.Bpm_swiginit(self, _RF_Track.new_Bpm(length, bpm_resolution, bpm_resolution_t))
+    get_reading = _swig_new_instance_method(_RF_Track.Bpm_get_reading)
+    get_reading_t = _swig_new_instance_method(_RF_Track.Bpm_get_reading_t)
+    set_resolution = _swig_new_instance_method(_RF_Track.Bpm_set_resolution)
+    set_resolution_t = _swig_new_instance_method(_RF_Track.Bpm_set_resolution_t)
+    set_scaling_factor = _swig_new_instance_method(_RF_Track.Bpm_set_scaling_factor)
+    get_resolution = _swig_new_instance_method(_RF_Track.Bpm_get_resolution)
+    get_resolution_t = _swig_new_instance_method(_RF_Track.Bpm_get_resolution_t)
+    get_total_charge = _swig_new_instance_method(_RF_Track.Bpm_get_total_charge)
+    __swig_destroy__ = _RF_Track.delete_Bpm
+
+# Register Bpm in _RF_Track:
+_RF_Track.Bpm_swigregister(Bpm)
+class Drift(ConstantField):
+    r"""
+
+    A drift space of length in metres, field-free by default. Inherited methods can add
+    static electric or magnetic fields and collective effects.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, length=0.0):
+        r"""
+        __init__(Drift self, double length=0.0) -> Drift
+
+        A drift space of length in metres, field-free by default. Inherited methods can add
+        static electric or magnetic fields and collective effects.
+
+
+        """
+        _RF_Track.Drift_swiginit(self, _RF_Track.new_Drift(length))
+    __swig_destroy__ = _RF_Track.delete_Drift
+
+# Register Drift in _RF_Track:
+_RF_Track.Drift_swigregister(Drift)
+class Corrector(Multipole):
+    r"""
+
+    A steering magnet. Length is in m and horizontal/vertical integrated strengths Sx, Sy in
+    T.mm. The rigidity-based overload takes P_over_Q in MV/c and kicks in mrad.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Corrector self, double length=0.0, double Sx=0.0, double Sy=0.0) -> Corrector
+        __init__(Corrector self, double length, double P_over_Q, double kxL, double kyL) -> Corrector
+
+        A steering magnet. Length is in m and horizontal/vertical integrated strengths Sx, Sy in
+        T.mm. The rigidity-based overload takes P_over_Q in MV/c and kicks in mrad.
+
+
+        """
+        _RF_Track.Corrector_swiginit(self, _RF_Track.new_Corrector(*args))
+    set_kick = _swig_new_instance_method(_RF_Track.Corrector_set_kick)
+    set_strength = _swig_new_instance_method(_RF_Track.Corrector_set_strength)
+    vary_strength = _swig_new_instance_method(_RF_Track.Corrector_vary_strength)
+    get_kick = _swig_new_instance_method(_RF_Track.Corrector_get_kick)
+    get_strength = _swig_new_instance_method(_RF_Track.Corrector_get_strength)
+    __swig_destroy__ = _RF_Track.delete_Corrector
+
+# Register Corrector in _RF_Track:
+_RF_Track.Corrector_swigregister(Corrector)
+class Screen(TimeDependent_Field):
+    r"""
+
+    A zero-length diagnostic plane that records beam distributions crossing it. Retrieve
+    recorded distributions from the enclosing Lattice or Volume.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self):
+        r"""
+        __init__(Screen self) -> Screen
+
+        A zero-length diagnostic plane that records beam distributions crossing it. Retrieve
+        recorded distributions from the enclosing Lattice or Volume.
+
+
+        """
+        _RF_Track.Screen_swiginit(self, _RF_Track.new_Screen())
+    get_bunch = _swig_new_instance_method(_RF_Track.Screen_get_bunch)
+    get_beam = _swig_new_instance_method(_RF_Track.Screen_get_beam)
+    set_time_window = _swig_new_instance_method(_RF_Track.Screen_set_time_window)
+    __swig_destroy__ = _RF_Track.delete_Screen
+
+# Register Screen in _RF_Track:
+_RF_Track.Screen_swigregister(Screen)
+class RF_Field(TimeDependent_Field):
+    r"""
+
+    Base interface for radiofrequency fields. Frequency is in Hz. set_phi uses radians and
+    set_phid degrees; the phase offset is relative to the synchronized accelerating crest.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined")
+    __repr__ = _swig_repr
+    set_P_map = _swig_new_instance_method(_RF_Track.RF_Field_set_P_map)
+    set_P_actual = _swig_new_instance_method(_RF_Track.RF_Field_set_P_actual)
+    get_phi = _swig_new_instance_method(_RF_Track.RF_Field_get_phi)
+    get_phid = _swig_new_instance_method(_RF_Track.RF_Field_get_phid)
+    get_period = _swig_new_instance_method(_RF_Track.RF_Field_get_period)
+    set_phi = _swig_new_instance_method(_RF_Track.RF_Field_set_phi)
+    set_phid = _swig_new_instance_method(_RF_Track.RF_Field_set_phid)
+    get_frequency = _swig_new_instance_method(_RF_Track.RF_Field_get_frequency)
+    set_frequency = _swig_new_instance_method(_RF_Track.RF_Field_set_frequency)
+    get_field_abs = _swig_new_instance_method(_RF_Track.RF_Field_get_field_abs)
+    get_field_complex = _swig_new_instance_method(_RF_Track.RF_Field_get_field_complex)
+    __swig_destroy__ = _RF_Track.delete_RF_Field
+
+# Register RF_Field in _RF_Track:
+_RF_Track.RF_Field_swigregister(RF_Field)
+class LaserBeam(TimeDependent_Field):
+    r"""
+
+    A laser pulse for inverse Compton scattering. Configure the pulse properties and
+    interaction geometry before tracking. See manual chapter 7 for photon collection and
+    polarization conventions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    tLaser_IP = property(_RF_Track.LaserBeam_tLaser_IP_get, _RF_Track.LaserBeam_tLaser_IP_set, doc=r"""tLaser_IP : double""")
+    length = property(_RF_Track.LaserBeam_length_get, _RF_Track.LaserBeam_length_set, doc=r"""length : double""")
+    pulse_energy = property(_RF_Track.LaserBeam_pulse_energy_get, _RF_Track.LaserBeam_pulse_energy_set, doc=r"""pulse_energy : double""")
+    pulse_length = property(_RF_Track.LaserBeam_pulse_length_get, _RF_Track.LaserBeam_pulse_length_set, doc=r"""pulse_length : double""")
+    wavelength = property(_RF_Track.LaserBeam_wavelength_get, _RF_Track.LaserBeam_wavelength_set, doc=r"""wavelength : double""")
+    position = property(_RF_Track.LaserBeam_position_get, _RF_Track.LaserBeam_position_set, doc=r"""position : StaticVector<(3)>""")
+    direction = property(_RF_Track.LaserBeam_direction_get, _RF_Track.LaserBeam_direction_set, doc=r"""direction : StaticVector<(3)>""")
+    Rx = property(_RF_Track.LaserBeam_Rx_get, _RF_Track.LaserBeam_Rx_set, doc=r"""Rx : double""")
+    Ry = property(_RF_Track.LaserBeam_Ry_get, _RF_Track.LaserBeam_Ry_set, doc=r"""Ry : double""")
+    R = property(_RF_Track.LaserBeam_R_get, _RF_Track.LaserBeam_R_set, doc=r"""R : double""")
+    M2x = property(_RF_Track.LaserBeam_M2x_get, _RF_Track.LaserBeam_M2x_set, doc=r"""M2x : double""")
+    M2y = property(_RF_Track.LaserBeam_M2y_get, _RF_Track.LaserBeam_M2y_set, doc=r"""M2y : double""")
+    M2 = property(_RF_Track.LaserBeam_M2_get, _RF_Track.LaserBeam_M2_set, doc=r"""M2 : double""")
+    P = property(_RF_Track.LaserBeam_P_get, _RF_Track.LaserBeam_P_set, doc=r"""P : double""")
+    lambda_e = property(_RF_Track.LaserBeam_lambda_e_get, _RF_Track.LaserBeam_lambda_e_set, doc=r"""lambda_e : double""")
+    rep_frequency = property(_RF_Track.LaserBeam_rep_frequency_get, _RF_Track.LaserBeam_rep_frequency_set, doc=r"""rep_frequency : double""")
+    number_of_pulses = property(_RF_Track.LaserBeam_number_of_pulses_get, _RF_Track.LaserBeam_number_of_pulses_set, doc=r"""number_of_pulses : size_t""")
+    min_number_of_gammas_per_slice = property(_RF_Track.LaserBeam_min_number_of_gammas_per_slice_get, _RF_Track.LaserBeam_min_number_of_gammas_per_slice_set, doc=r"""min_number_of_gammas_per_slice : size_t""")
+    Nmin = property(_RF_Track.LaserBeam_Nmin_get, _RF_Track.LaserBeam_Nmin_set, doc=r"""Nmin : double""")
+    check_R = _swig_new_instance_method(_RF_Track.LaserBeam_check_R)
+    check_M2 = _swig_new_instance_method(_RF_Track.LaserBeam_check_M2)
+
+    def __init__(self, length=0.0):
+        r"""
+        __init__(LaserBeam self, double length=0.0) -> LaserBeam
+
+        A laser pulse for inverse Compton scattering. Configure the pulse properties and
+        interaction geometry before tracking. See manual chapter 7 for photon collection and
+        polarization conventions.
+
+
+        """
+        _RF_Track.LaserBeam_swiginit(self, _RF_Track.new_LaserBeam(length))
+    get_M2 = _swig_new_instance_method(_RF_Track.LaserBeam_get_M2)
+    get_M2x = _swig_new_instance_method(_RF_Track.LaserBeam_get_M2x)
+    get_M2y = _swig_new_instance_method(_RF_Track.LaserBeam_get_M2y)
+    get_zR = _swig_new_instance_method(_RF_Track.LaserBeam_get_zR)
+    get_zRx = _swig_new_instance_method(_RF_Track.LaserBeam_get_zRx)
+    get_zRy = _swig_new_instance_method(_RF_Track.LaserBeam_get_zRy)
+    get_direction = _swig_new_instance_method(_RF_Track.LaserBeam_get_direction)
+    get_position = _swig_new_instance_method(_RF_Track.LaserBeam_get_position)
+    set_direction = _swig_new_instance_method(_RF_Track.LaserBeam_set_direction)
+    set_position = _swig_new_instance_method(_RF_Track.LaserBeam_set_position)
+    set_M2 = _swig_new_instance_method(_RF_Track.LaserBeam_set_M2)
+    set_M2x = _swig_new_instance_method(_RF_Track.LaserBeam_set_M2x)
+    set_M2y = _swig_new_instance_method(_RF_Track.LaserBeam_set_M2y)
+    set_zR = _swig_new_instance_method(_RF_Track.LaserBeam_set_zR)
+    set_zRx = _swig_new_instance_method(_RF_Track.LaserBeam_set_zRx)
+    set_zRy = _swig_new_instance_method(_RF_Track.LaserBeam_set_zRy)
+    __swig_destroy__ = _RF_Track.delete_LaserBeam
+
+# Register LaserBeam in _RF_Track:
+_RF_Track.LaserBeam_swigregister(LaserBeam)
+class Pillbox_Cavity(RF_Field):
+    r"""
+
+    Analytic pillbox RF cavity represented by mode coefficients. Frequency is in Hz and cell
+    length in metres. Inherits RF phase and timing controls.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Pillbox_Cavity self, MatrixNd const & coeffs_, double frequency, double cell_length, double number_of_cells) -> Pillbox_Cavity
+        __init__(Pillbox_Cavity self) -> Pillbox_Cavity
+
+        Analytic pillbox RF cavity represented by mode coefficients. Frequency is in Hz and cell
+        length in metres. Inherits RF phase and timing controls.
+
+
+        """
+        _RF_Track.Pillbox_Cavity_swiginit(self, _RF_Track.new_Pillbox_Cavity(*args))
+    get_cell_length = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_get_cell_length)
+    get_z0 = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_get_z1)
+    get_coefficients = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_get_coefficients)
+    set_coefficients = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_set_coefficients)
+    set_z0 = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_set_z1)
+    set_cell_length = _swig_new_instance_method(_RF_Track.Pillbox_Cavity_set_cell_length)
+    __swig_destroy__ = _RF_Track.delete_Pillbox_Cavity
+
+# Register Pillbox_Cavity in _RF_Track:
+_RF_Track.Pillbox_Cavity_swigregister(Pillbox_Cavity)
+class SW_Structure(RF_Field):
+    r"""
+
+    Analytic standing-wave accelerating structure. Fourier coefficients are in V/m,
+    frequency in Hz and cell_length in m; number_of_cells can be nonintegral. Phase and
+    timing are controlled by inherited RF_Field methods.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SW_Structure self, std::vector< double > const & coeffs_, double frequency, double cell_length, double number_of_cells) -> SW_Structure
+        __init__(SW_Structure self, double coeffs_, double frequency, double cell_length, double number_of_cells) -> SW_Structure
+        __init__(SW_Structure self) -> SW_Structure
+
+        Analytic standing-wave accelerating structure. Fourier coefficients are in V/m,
+        frequency in Hz and cell_length in m; number_of_cells can be nonintegral. Phase and
+        timing are controlled by inherited RF_Field methods.
+
+
+        """
+        _RF_Track.SW_Structure_swiginit(self, _RF_Track.new_SW_Structure(*args))
+    get_cell_length = _swig_new_instance_method(_RF_Track.SW_Structure_get_cell_length)
+    get_z0 = _swig_new_instance_method(_RF_Track.SW_Structure_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.SW_Structure_get_z1)
+    get_coefficients = _swig_new_instance_method(_RF_Track.SW_Structure_get_coefficients)
+    set_coefficients = _swig_new_instance_method(_RF_Track.SW_Structure_set_coefficients)
+    set_z0 = _swig_new_instance_method(_RF_Track.SW_Structure_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.SW_Structure_set_z1)
+    set_cell_length = _swig_new_instance_method(_RF_Track.SW_Structure_set_cell_length)
+    __swig_destroy__ = _RF_Track.delete_SW_Structure
+
+# Register SW_Structure in _RF_Track:
+_RF_Track.SW_Structure_swigregister(SW_Structure)
+class TW_Structure(RF_Field):
+    r"""
+
+    Analytic travelling-wave accelerating structure described by spatial Fourier
+    coefficients. Field coefficients are in V/m, frequency in Hz and cell length in metres.
+    Consult the constructor signature for phase advance and cell count.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(TW_Structure self, std::vector< double > const & coeffs_, int first_suffix, double frequency, double phase_advance, double number_of_cells) -> TW_Structure
+        __init__(TW_Structure self, double coeffs_, int first_suffix, double frequency, double phase_advance, double number_of_cells) -> TW_Structure
+        __init__(TW_Structure self) -> TW_Structure
+
+        Analytic travelling-wave accelerating structure described by spatial Fourier
+        coefficients. Field coefficients are in V/m, frequency in Hz and cell length in metres.
+        Consult the constructor signature for phase advance and cell count.
+
+
+        """
+        _RF_Track.TW_Structure_swiginit(self, _RF_Track.new_TW_Structure(*args))
+    get_cell_length = _swig_new_instance_method(_RF_Track.TW_Structure_get_cell_length)
+    get_z0 = _swig_new_instance_method(_RF_Track.TW_Structure_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.TW_Structure_get_z1)
+    get_coefficients = _swig_new_instance_method(_RF_Track.TW_Structure_get_coefficients)
+    set_coefficients = _swig_new_instance_method(_RF_Track.TW_Structure_set_coefficients)
+    set_z0 = _swig_new_instance_method(_RF_Track.TW_Structure_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.TW_Structure_set_z1)
+    __swig_destroy__ = _RF_Track.delete_TW_Structure
+
+# Register TW_Structure in _RF_Track:
+_RF_Track.TW_Structure_swigregister(TW_Structure)
+class TM_mnp(RF_Field):
+    r"""Proxy of C++ TM_mnp class."""
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    get_field_complex_analytic = _swig_new_instance_method(_RF_Track.TM_mnp_get_field_complex_analytic)
+
+    def __init__(self, a_mnp, first_m, first_p, frequency, phase_advance, number_of_cells, radius, cell_length=-1.0, nx=41, ny=41, nz=0):
+        r"""__init__(TM_mnp self, ComplexMesh3d const & a_mnp, int first_m, int first_p, double frequency, double phase_advance, double number_of_cells, double radius, double cell_length=-1.0, size_t nx=41, size_t ny=41, size_t nz=0) -> TM_mnp"""
+        _RF_Track.TM_mnp_swiginit(self, _RF_Track.new_TM_mnp(a_mnp, first_m, first_p, frequency, phase_advance, number_of_cells, radius, cell_length, nx, ny, nz))
+    get_coefficients = _swig_new_instance_method(_RF_Track.TM_mnp_get_coefficients)
+    get_first_m = _swig_new_instance_method(_RF_Track.TM_mnp_get_first_m)
+    get_first_p = _swig_new_instance_method(_RF_Track.TM_mnp_get_first_p)
+    get_cell_length = _swig_new_instance_method(_RF_Track.TM_mnp_get_cell_length)
+    get_phase_advance = _swig_new_instance_method(_RF_Track.TM_mnp_get_phase_advance)
+    get_number_of_cells = _swig_new_instance_method(_RF_Track.TM_mnp_get_number_of_cells)
+    get_radius = _swig_new_instance_method(_RF_Track.TM_mnp_get_radius)
+    get_nx = _swig_new_instance_method(_RF_Track.TM_mnp_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.TM_mnp_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.TM_mnp_get_nz)
+    get_hx = _swig_new_instance_method(_RF_Track.TM_mnp_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.TM_mnp_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.TM_mnp_get_hz)
+    set_coefficients = _swig_new_instance_method(_RF_Track.TM_mnp_set_coefficients)
+    set_cell_length = _swig_new_instance_method(_RF_Track.TM_mnp_set_cell_length)
+    set_phase_advance = _swig_new_instance_method(_RF_Track.TM_mnp_set_phase_advance)
+    set_number_of_cells = _swig_new_instance_method(_RF_Track.TM_mnp_set_number_of_cells)
+    set_radius = _swig_new_instance_method(_RF_Track.TM_mnp_set_radius)
+    set_mesh_size = _swig_new_instance_method(_RF_Track.TM_mnp_set_mesh_size)
+    __swig_destroy__ = _RF_Track.delete_TM_mnp
+
+# Register TM_mnp in _RF_Track:
+_RF_Track.TM_mnp_swigregister(TM_mnp)
+class TE_mnp(RF_Field):
+    r"""Proxy of C++ TE_mnp class."""
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    get_field_complex_analytic = _swig_new_instance_method(_RF_Track.TE_mnp_get_field_complex_analytic)
+
+    def __init__(self, a_mnp, first_m, first_p, frequency, phase_advance, number_of_cells, radius, cell_length=-1.0, nx=41, ny=41, nz=0):
+        r"""__init__(TE_mnp self, ComplexMesh3d const & a_mnp, int first_m, int first_p, double frequency, double phase_advance, double number_of_cells, double radius, double cell_length=-1.0, size_t nx=41, size_t ny=41, size_t nz=0) -> TE_mnp"""
+        _RF_Track.TE_mnp_swiginit(self, _RF_Track.new_TE_mnp(a_mnp, first_m, first_p, frequency, phase_advance, number_of_cells, radius, cell_length, nx, ny, nz))
+    get_coefficients = _swig_new_instance_method(_RF_Track.TE_mnp_get_coefficients)
+    get_first_m = _swig_new_instance_method(_RF_Track.TE_mnp_get_first_m)
+    get_first_p = _swig_new_instance_method(_RF_Track.TE_mnp_get_first_p)
+    get_cell_length = _swig_new_instance_method(_RF_Track.TE_mnp_get_cell_length)
+    get_phase_advance = _swig_new_instance_method(_RF_Track.TE_mnp_get_phase_advance)
+    get_number_of_cells = _swig_new_instance_method(_RF_Track.TE_mnp_get_number_of_cells)
+    get_radius = _swig_new_instance_method(_RF_Track.TE_mnp_get_radius)
+    get_nx = _swig_new_instance_method(_RF_Track.TE_mnp_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.TE_mnp_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.TE_mnp_get_nz)
+    get_hx = _swig_new_instance_method(_RF_Track.TE_mnp_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.TE_mnp_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.TE_mnp_get_hz)
+    set_coefficients = _swig_new_instance_method(_RF_Track.TE_mnp_set_coefficients)
+    set_cell_length = _swig_new_instance_method(_RF_Track.TE_mnp_set_cell_length)
+    set_phase_advance = _swig_new_instance_method(_RF_Track.TE_mnp_set_phase_advance)
+    set_number_of_cells = _swig_new_instance_method(_RF_Track.TE_mnp_set_number_of_cells)
+    set_radius = _swig_new_instance_method(_RF_Track.TE_mnp_set_radius)
+    set_mesh_size = _swig_new_instance_method(_RF_Track.TE_mnp_set_mesh_size)
+    __swig_destroy__ = _RF_Track.delete_TE_mnp
+
+# Register TE_mnp in _RF_Track:
+_RF_Track.TE_mnp_swigregister(TE_mnp)
+class RF_FieldMap_1d_CINT(RF_Field):
+    r"""
+
+    RF field reconstructed from a regular on-axis Ez map in V/m, using cylindrical symmetry
+    and Maxwell equations. hz and length are in m; length=-1 uses the map length. Frequency
+    is in Hz, direction is +1/-1 for travelling fields or 0 for static fields, and
+    P_map/P_actual are in W. The default exported class uses linear interpolation; the _CINT
+    variant uses cubic interpolation.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap_1d_CINT self, ComplexMesh1d const & Ez, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_1d_CINT
+        __init__(RF_FieldMap_1d_CINT self) -> RF_FieldMap_1d_CINT
+
+        RF field reconstructed from a regular on-axis Ez map in V/m, using cylindrical symmetry
+        and Maxwell equations. hz and length are in m; length=-1 uses the map length. Frequency
+        is in Hz, direction is +1/-1 for travelling fields or 0 for static fields, and
+        P_map/P_actual are in W. The default exported class uses linear interpolation; the _CINT
+        variant uses cubic interpolation.
+
+
+        """
+        _RF_Track.RF_FieldMap_1d_CINT_swiginit(self, _RF_Track.new_RF_FieldMap_1d_CINT(*args))
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_hz)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_z1)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_direction)
+    get_Ez = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_get_Ez)
+    set_smooth = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_set_smooth)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_set_z1)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_CINT_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap_1d_CINT
+
+# Register RF_FieldMap_1d_CINT in _RF_Track:
+_RF_Track.RF_FieldMap_1d_CINT_swigregister(RF_FieldMap_1d_CINT)
+class RF_FieldMap_1d(RF_Field):
+    r"""
+
+
+
+    RF field reconstructed from a regular on-axis Ez map in V/m, using cylindrical symmetry
+    and Maxwell equations. hz and length are in m; length=-1 uses the map length. Frequency
+    is in Hz, direction is +1/-1 for travelling fields or 0 for static fields, and
+    P_map/P_actual are in W. The default exported class uses linear interpolation; the _CINT
+    variant uses cubic interpolation.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap_1d self, ComplexMesh1d const & Ez, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_1d
+        __init__(RF_FieldMap_1d self) -> RF_FieldMap_1d
+
+        RF field reconstructed from a regular on-axis Ez map in V/m, using cylindrical symmetry
+        and Maxwell equations. hz and length are in m; length=-1 uses the map length. Frequency
+        is in Hz, direction is +1/-1 for travelling fields or 0 for static fields, and
+        P_map/P_actual are in W. The default exported class uses linear interpolation; the _CINT
+        variant uses cubic interpolation.
+
+
+        """
+        _RF_Track.RF_FieldMap_1d_swiginit(self, _RF_Track.new_RF_FieldMap_1d(*args))
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_hz)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_z1)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_direction)
+    get_Ez = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_get_Ez)
+    set_smooth = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_set_smooth)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_set_z1)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_1d_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap_1d
+
+# Register RF_FieldMap_1d in _RF_Track:
+_RF_Track.RF_FieldMap_1d_swigregister(RF_FieldMap_1d)
+class TW_Field(RF_FieldMap_1d_CINT):
+    r"""
+
+    Travelling-wave RF field constructed from input power and longitudinal profiles of
+    quality factor, shunt impedance over Q, and group velocity. Inherits RF field-map phase
+    and power controls.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(TW_Field self, double P_in, Mesh1d const & QQ, Mesh1d const & R_Q, Mesh1d const & VG, double frequency, double phase_advance, double number_of_cells, double z0_L=0.0) -> TW_Field
+        __init__(TW_Field self, Mesh1d const & Pinput, double dt_Pinput, double t_inj, Mesh1d const & QQ, Mesh1d const & R_Q, Mesh1d const & VG, double frequency, double phase_advance, double number_of_cells) -> TW_Field
+        __init__(TW_Field self) -> TW_Field
+
+        Travelling-wave RF field constructed from input power and longitudinal profiles of
+        quality factor, shunt impedance over Q, and group velocity. Inherits RF field-map phase
+        and power controls.
+
+
+        """
+        _RF_Track.TW_Field_swiginit(self, _RF_Track.new_TW_Field(*args))
+    __swig_destroy__ = _RF_Track.delete_TW_Field
+
+# Register TW_Field in _RF_Track:
+_RF_Track.TW_Field_swigregister(TW_Field)
+class Particle(object):
+    r"""
+
+    One macroparticle in the spatial tracking model: x,y [mm], xp,yp [mrad], t [mm/c], Pc
+    [MeV/c], mass [MeV/c^2], Q [e] and population N. Usually accessed through Bunch6d.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    mass = property(_RF_Track.Particle_mass_get, _RF_Track.Particle_mass_set, doc=r"""mass : double""")
+    Q = property(_RF_Track.Particle_Q_get, _RF_Track.Particle_Q_set, doc=r"""Q : double""")
+    N = property(_RF_Track.Particle_N_get, _RF_Track.Particle_N_set, doc=r"""N : double""")
+    x = property(_RF_Track.Particle_x_get, _RF_Track.Particle_x_set, doc=r"""x : double""")
+    xp = property(_RF_Track.Particle_xp_get, _RF_Track.Particle_xp_set, doc=r"""xp : double""")
+    y = property(_RF_Track.Particle_y_get, _RF_Track.Particle_y_set, doc=r"""y : double""")
+    yp = property(_RF_Track.Particle_yp_get, _RF_Track.Particle_yp_set, doc=r"""yp : double""")
+    t = property(_RF_Track.Particle_t_get, _RF_Track.Particle_t_set, doc=r"""t : double""")
+    Pc = property(_RF_Track.Particle_Pc_get, _RF_Track.Particle_Pc_set, doc=r"""Pc : double""")
+    S_lost = property(_RF_Track.Particle_S_lost_get, _RF_Track.Particle_S_lost_set, doc=r"""S_lost : double""")
+    lifetime = property(_RF_Track.Particle_lifetime_get, _RF_Track.Particle_lifetime_set, doc=r"""lifetime : double""")
+    Sx = property(_RF_Track.Particle_Sx_get, _RF_Track.Particle_Sx_set, doc=r"""Sx : double""")
+    Sy = property(_RF_Track.Particle_Sy_get, _RF_Track.Particle_Sy_set, doc=r"""Sy : double""")
+    Sz = property(_RF_Track.Particle_Sz_get, _RF_Track.Particle_Sz_set, doc=r"""Sz : double""")
+    anomalous_magnetic_moment = property(_RF_Track.Particle_anomalous_magnetic_moment_get, _RF_Track.Particle_anomalous_magnetic_moment_set, doc=r"""anomalous_magnetic_moment : double""")
+
+    def __init__(self, *args):
+        r"""
+        __init__(Particle self) -> Particle
+        __init__(Particle self, double x_, double xp_, double y_, double yp_, double t_, double Pc_) -> Particle
+
+        One macroparticle in the spatial tracking model: x,y [mm], xp,yp [mrad], t [mm/c], Pc
+        [MeV/c], mass [MeV/c^2], Q [e] and population N. Usually accessed through Bunch6d.
+
+
+        """
+        _RF_Track.Particle_swiginit(self, _RF_Track.new_Particle(*args))
+
+    def __nonzero__(self):
+        return _RF_Track.Particle_is_valid(self)
+    __bool__ = __nonzero__
+
+
+    lost_at = _swig_new_instance_method(_RF_Track.Particle_lost_at)
+    get_Pc_sqr = _swig_new_instance_method(_RF_Track.Particle_get_Pc_sqr)
+    get_Pc = _swig_new_instance_method(_RF_Track.Particle_get_Pc)
+    get_total_energy_sqr = _swig_new_instance_method(_RF_Track.Particle_get_total_energy_sqr)
+    get_total_energy = _swig_new_instance_method(_RF_Track.Particle_get_total_energy)
+    get_kinetic_energy = _swig_new_instance_method(_RF_Track.Particle_get_kinetic_energy)
+    get_delta_plus_1 = _swig_new_instance_method(_RF_Track.Particle_get_delta_plus_1)
+    get_delta = _swig_new_instance_method(_RF_Track.Particle_get_delta)
+    get_beta_sqr = _swig_new_instance_method(_RF_Track.Particle_get_beta_sqr)
+    get_beta = _swig_new_instance_method(_RF_Track.Particle_get_beta)
+    get_gamma_sqr = _swig_new_instance_method(_RF_Track.Particle_get_gamma_sqr)
+    get_gamma = _swig_new_instance_method(_RF_Track.Particle_get_gamma)
+    get_pt = _swig_new_instance_method(_RF_Track.Particle_get_pt)
+    get_xp_yp_zp = _swig_new_instance_method(_RF_Track.Particle_get_xp_yp_zp)
+    get_Vx_Vy_Vz = _swig_new_instance_method(_RF_Track.Particle_get_Vx_Vy_Vz)
+    get_Px_Py_Pz = _swig_new_instance_method(_RF_Track.Particle_get_Px_Py_Pz)
+    get_px_py_pz = _swig_new_instance_method(_RF_Track.Particle_get_px_py_pz)
+    get_Sx_Sy_Sz = _swig_new_instance_method(_RF_Track.Particle_get_Sx_Sy_Sz)
+    get_four_momentum = _swig_new_instance_method(_RF_Track.Particle_get_four_momentum)
+    get_four_velocity = _swig_new_instance_method(_RF_Track.Particle_get_four_velocity)
+    get_id = _swig_new_instance_method(_RF_Track.Particle_get_id)
+    set_id = _swig_new_instance_method(_RF_Track.Particle_set_id)
+    get_pdg_id = _swig_new_instance_method(_RF_Track.Particle_get_pdg_id)
+    set_pdg_id = _swig_new_instance_method(_RF_Track.Particle_set_pdg_id)
+    __swig_destroy__ = _RF_Track.delete_Particle
+
+# Register Particle in _RF_Track:
+_RF_Track.Particle_swigregister(Particle)
+class ParticleT(object):
+    r"""
+
+    One macroparticle in the time tracking model: X,Y,Z [mm], Px,Py,Pz [MeV/c], mass
+    [MeV/c^2], Q [e] and population N. Creation time is in mm/c. Usually accessed through
+    Bunch6dT.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    mass = property(_RF_Track.ParticleT_mass_get, _RF_Track.ParticleT_mass_set, doc=r"""mass : double""")
+    Q = property(_RF_Track.ParticleT_Q_get, _RF_Track.ParticleT_Q_set, doc=r"""Q : double""")
+    N = property(_RF_Track.ParticleT_N_get, _RF_Track.ParticleT_N_set, doc=r"""N : double""")
+    X = property(_RF_Track.ParticleT_X_get, _RF_Track.ParticleT_X_set, doc=r"""X : double""")
+    Px = property(_RF_Track.ParticleT_Px_get, _RF_Track.ParticleT_Px_set, doc=r"""Px : double""")
+    Y = property(_RF_Track.ParticleT_Y_get, _RF_Track.ParticleT_Y_set, doc=r"""Y : double""")
+    Py = property(_RF_Track.ParticleT_Py_get, _RF_Track.ParticleT_Py_set, doc=r"""Py : double""")
+    Z = property(_RF_Track.ParticleT_Z_get, _RF_Track.ParticleT_Z_set, doc=r"""Z : double""")
+    Pz = property(_RF_Track.ParticleT_Pz_get, _RF_Track.ParticleT_Pz_set, doc=r"""Pz : double""")
+    t0 = property(_RF_Track.ParticleT_t0_get, _RF_Track.ParticleT_t0_set, doc=r"""t0 : double""")
+    t1 = property(_RF_Track.ParticleT_t1_get, _RF_Track.ParticleT_t1_set, doc=r"""t1 : double""")
+    t_lost = property(_RF_Track.ParticleT_t_lost_get, _RF_Track.ParticleT_t_lost_set, doc=r"""t_lost : double""")
+    lifetime = property(_RF_Track.ParticleT_lifetime_get, _RF_Track.ParticleT_lifetime_set, doc=r"""lifetime : double""")
+    Sx = property(_RF_Track.ParticleT_Sx_get, _RF_Track.ParticleT_Sx_set, doc=r"""Sx : double""")
+    Sy = property(_RF_Track.ParticleT_Sy_get, _RF_Track.ParticleT_Sy_set, doc=r"""Sy : double""")
+    Sz = property(_RF_Track.ParticleT_Sz_get, _RF_Track.ParticleT_Sz_set, doc=r"""Sz : double""")
+    anomalous_magnetic_moment = property(_RF_Track.ParticleT_anomalous_magnetic_moment_get, _RF_Track.ParticleT_anomalous_magnetic_moment_set, doc=r"""anomalous_magnetic_moment : double""")
+
+    def __init__(self, *args):
+        r"""
+        __init__(ParticleT self) -> ParticleT
+        __init__(ParticleT self, double X_, double xp_, double Y_, double yp_, double Z_, double Pc_) -> ParticleT
+
+        One macroparticle in the time tracking model: X,Y,Z [mm], Px,Py,Pz [MeV/c], mass
+        [MeV/c^2], Q [e] and population N. Creation time is in mm/c. Usually accessed through
+        Bunch6dT.
+
+
+        """
+        _RF_Track.ParticleT_swiginit(self, _RF_Track.new_ParticleT(*args))
+
+    def __nonzero__(self):
+        return _RF_Track.ParticleT_is_valid(self)
+    __bool__ = __nonzero__
+
+
+    lost_at = _swig_new_instance_method(_RF_Track.ParticleT_lost_at)
+    get_Pc_sqr = _swig_new_instance_method(_RF_Track.ParticleT_get_Pc_sqr)
+    get_Pc = _swig_new_instance_method(_RF_Track.ParticleT_get_Pc)
+    get_total_energy_sqr = _swig_new_instance_method(_RF_Track.ParticleT_get_total_energy_sqr)
+    get_total_energy = _swig_new_instance_method(_RF_Track.ParticleT_get_total_energy)
+    get_kinetic_energy = _swig_new_instance_method(_RF_Track.ParticleT_get_kinetic_energy)
+    get_delta_plus_1 = _swig_new_instance_method(_RF_Track.ParticleT_get_delta_plus_1)
+    get_delta = _swig_new_instance_method(_RF_Track.ParticleT_get_delta)
+    get_beta_sqr = _swig_new_instance_method(_RF_Track.ParticleT_get_beta_sqr)
+    get_beta = _swig_new_instance_method(_RF_Track.ParticleT_get_beta)
+    get_gamma_sqr = _swig_new_instance_method(_RF_Track.ParticleT_get_gamma_sqr)
+    get_gamma = _swig_new_instance_method(_RF_Track.ParticleT_get_gamma)
+    get_pt = _swig_new_instance_method(_RF_Track.ParticleT_get_pt)
+    get_xp_yp_zp = _swig_new_instance_method(_RF_Track.ParticleT_get_xp_yp_zp)
+    get_Vx_Vy_Vz = _swig_new_instance_method(_RF_Track.ParticleT_get_Vx_Vy_Vz)
+    get_Px_Py_Pz = _swig_new_instance_method(_RF_Track.ParticleT_get_Px_Py_Pz)
+    get_px_py_pz = _swig_new_instance_method(_RF_Track.ParticleT_get_px_py_pz)
+    get_Sx_Sy_Sz = _swig_new_instance_method(_RF_Track.ParticleT_get_Sx_Sy_Sz)
+    get_four_momentum = _swig_new_instance_method(_RF_Track.ParticleT_get_four_momentum)
+    get_four_velocity = _swig_new_instance_method(_RF_Track.ParticleT_get_four_velocity)
+    get_id = _swig_new_instance_method(_RF_Track.ParticleT_get_id)
+    set_id = _swig_new_instance_method(_RF_Track.ParticleT_set_id)
+    get_pdg_id = _swig_new_instance_method(_RF_Track.ParticleT_get_pdg_id)
+    set_pdg_id = _swig_new_instance_method(_RF_Track.ParticleT_set_pdg_id)
+    __swig_destroy__ = _RF_Track.delete_ParticleT
+
+# Register ParticleT in _RF_Track:
+_RF_Track.ParticleT_swigregister(ParticleT)
+class Bunch6dT(object):
+    r"""
+
+    A bunch represented at a common time t, for Volume tracking.
+
+    Each phase-space matrix row is one macroparticle. The six standard columns
+    are [X, Px, Y, Py, Z, Pz]: positions in mm and momenta in MeV/c.
+    Extended rows add [mass, Q, N, T0] and optionally ID; T0 is the creation
+    time in mm/c. Mass is in MeV/c^2, Q in units of e, and N is the number of
+    real particles represented by that macroparticle. The common clock t is
+    in mm/c. One mm/c is approximately 3.336 ps.
+    Construct from (mass, population, charge, matrix), an extended matrix,
+    Twiss parameters, Bunch6dT_Generator, Bunch6d, or a saved RF-Track file.
+    The no-argument constructor creates an empty bunch. The conversion from
+    Bunch6d reconstructs a simultaneous distribution; its optional S is in m.
+    For Twiss generation, Pref is in MeV/c and sigma_cut=0 disables truncation.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    t = property(_RF_Track.Bunch6dT_t_get, _RF_Track.Bunch6dT_t_set, doc=r"""
+    t : double
+
+    Common bunch time in mm/c.
+
+
+    """)
+
+    def __init__(self, *args):
+        r"""
+        __init__(Bunch6dT self, size_t N=0) -> Bunch6dT
+        __init__(Bunch6dT self, double mass, double population, double charge, MatrixNd const & X) -> Bunch6dT
+        __init__(Bunch6dT self, double mass, double population, double charge, double Pref, Bunch6dT_twiss T, size_t N, double sigma_cut=0.0) -> Bunch6dT
+        __init__(Bunch6dT self, MatrixNd const & X) -> Bunch6dT
+        __init__(Bunch6dT self, Bunch6dT_Generator g, size_t Nparticles) -> Bunch6dT
+        __init__(Bunch6dT self, Bunch6d b, double S=std::numeric_limits< double >::quiet_NaN(), bool backtracking=False) -> Bunch6dT
+        __init__(Bunch6dT self, Bunch6d b, Frame f, bool backtrackng=False) -> Bunch6dT
+        __init__(Bunch6dT self, std::string const & rft_filename) -> Bunch6dT
+
+        A bunch represented at a common time t, for Volume tracking.
+
+        Each phase-space matrix row is one macroparticle. The six standard columns
+        are [X, Px, Y, Py, Z, Pz]: positions in mm and momenta in MeV/c.
+        Extended rows add [mass, Q, N, T0] and optionally ID; T0 is the creation
+        time in mm/c. Mass is in MeV/c^2, Q in units of e, and N is the number of
+        real particles represented by that macroparticle. The common clock t is
+        in mm/c. One mm/c is approximately 3.336 ps.
+        Construct from (mass, population, charge, matrix), an extended matrix,
+        Twiss parameters, Bunch6dT_Generator, Bunch6d, or a saved RF-Track file.
+        The no-argument constructor creates an empty bunch. The conversion from
+        Bunch6d reconstructs a simultaneous distribution; its optional S is in m.
+        For Twiss generation, Pref is in MeV/c and sigma_cut=0 disables truncation.
+
+
+        """
+        _RF_Track.Bunch6dT_swiginit(self, _RF_Track.new_Bunch6dT(*args))
+    resize = _swig_new_instance_method(_RF_Track.Bunch6dT_resize)
+    clear = _swig_new_instance_method(_RF_Track.Bunch6dT_clear)
+    size = _swig_new_instance_method(_RF_Track.Bunch6dT_size)
+    displaced = _swig_new_instance_method(_RF_Track.Bunch6dT_displaced)
+    append = _swig_new_instance_method(_RF_Track.Bunch6dT_append)
+    set_lifetime = _swig_new_instance_method(_RF_Track.Bunch6dT_set_lifetime)
+    set_polarization = _swig_new_instance_method(_RF_Track.Bunch6dT_set_polarization)
+    coasting = _swig_new_instance_method(_RF_Track.Bunch6dT_coasting)
+    set_coasting = _swig_new_instance_method(_RF_Track.Bunch6dT_set_coasting)
+    get_info = _swig_new_instance_method(_RF_Track.Bunch6dT_get_info)
+    set_phase_space = _swig_new_instance_method(_RF_Track.Bunch6dT_set_phase_space)
+    get_particle = _swig_new_instance_method(_RF_Track.Bunch6dT_get_particle)
+    get_S_min = _swig_new_instance_method(_RF_Track.Bunch6dT_get_S_min)
+    get_S_max = _swig_new_instance_method(_RF_Track.Bunch6dT_get_S_max)
+    get_S_mean = _swig_new_instance_method(_RF_Track.Bunch6dT_get_S_mean)
+    get_t = _swig_new_instance_method(_RF_Track.Bunch6dT_get_t)
+    get_average_particle = _swig_new_instance_method(_RF_Track.Bunch6dT_get_average_particle)
+    get_reference_particle = _swig_new_instance_method(_RF_Track.Bunch6dT_get_reference_particle)
+    get_bunch_temperature = _swig_new_instance_method(_RF_Track.Bunch6dT_get_bunch_temperature)
+    get_phase_space = _swig_new_instance_method(_RF_Track.Bunch6dT_get_phase_space)
+    get_population = _swig_new_instance_method(_RF_Track.Bunch6dT_get_population)
+    get_ngood = _swig_new_instance_method(_RF_Track.Bunch6dT_get_ngood)
+    get_nlost = _swig_new_instance_method(_RF_Track.Bunch6dT_get_nlost)
+    get_total_charge = _swig_new_instance_method(_RF_Track.Bunch6dT_get_total_charge)
+    apply_force = _swig_new_instance_method(_RF_Track.Bunch6dT_apply_force)
+    kick = _swig_new_instance_method(_RF_Track.Bunch6dT_kick)
+    drift = _swig_new_instance_method(_RF_Track.Bunch6dT_drift)
+    load = _swig_new_instance_method(_RF_Track.Bunch6dT_load)
+    save = _swig_new_instance_method(_RF_Track.Bunch6dT_save)
+    save_as_dst_file = _swig_new_instance_method(_RF_Track.Bunch6dT_save_as_dst_file)
+    save_as_sdds_file = _swig_new_instance_method(_RF_Track.Bunch6dT_save_as_sdds_file)
+    __getitem__ = _swig_new_instance_method(_RF_Track.Bunch6dT___getitem__)
+    __len__ = _swig_new_instance_method(_RF_Track.Bunch6dT___len__)
+    __iter__ = _swig_new_instance_method(_RF_Track.Bunch6dT___iter__)
+    __swig_destroy__ = _RF_Track.delete_Bunch6dT
+
+# Register Bunch6dT in _RF_Track:
+_RF_Track.Bunch6dT_swigregister(Bunch6dT)
+class Bunch6d(object):
+    r"""
+
+    A bunch represented at a common longitudinal position S, for Lattice tracking.
+
+    Construct from (mass, population, charge, X), an extended matrix X, Twiss
+    parameters, a generator, or a saved RF-Track file. With no arguments the bunch
+    is empty. Each matrix row is one macroparticle.
+    The six standard columns are [x, xp, y, yp, t, Pc]: positions in mm, slopes
+    in mrad, arrival time in mm/c, and total momentum in MeV/c.
+    Extended rows add [mass, Q, N] and optionally ID; mass is in MeV/c^2,
+    Q is single-particle charge in units of e, and N is macroparticle population.
+    In the separate-species constructor, population is the total number of real
+    particles, not the number of macroparticles. S is in metres.
+    Twiss generation takes Pref in MeV/c and N macroparticles; sigma_cut=0
+    means no truncation. See Bunch6d_twiss for the distribution parameters.
+
+    Python example:
+        B = RF_Track.Bunch6d(RF_Track.electronmass, 1e6, -1,
+                            numpy.array([[0., 0., 0., 0., 0., 100.]]))
+        X = B.get_phase_space("%x %xp %y %yp %t %Pc")
+    Octave example (after loading RF_Track):
+        B = Bunch6d(electronmass, 1e6, -1, [0 0 0 0 0 100]);
+        X = B.get_phase_space('%x %xp %y %yp %t %Pc');
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    S = property(_RF_Track.Bunch6d_S_get, _RF_Track.Bunch6d_S_set, doc=r"""
+    S : double
+
+    Longitudinal bunch coordinate in metres.
+
+
+    """)
+
+    def __init__(self, *args):
+        r"""
+        __init__(Bunch6d self, size_t N=0) -> Bunch6d
+        __init__(Bunch6d self, double mass, double population, double charge, MatrixNd const & X) -> Bunch6d
+        __init__(Bunch6d self, double mass, double population, double charge, double Pref, Bunch6d_twiss T, size_t N, double sigma_cut=0.0) -> Bunch6d
+        __init__(Bunch6d self, MatrixNd const & X, double S_position=0.0) -> Bunch6d
+        __init__(Bunch6d self, Bunch6dT_Generator g, size_t Nparticles) -> Bunch6d
+        __init__(Bunch6d self, std::string const & rft_filename) -> Bunch6d
+
+        A bunch represented at a common longitudinal position S, for Lattice tracking.
+
+        Construct from (mass, population, charge, X), an extended matrix X, Twiss
+        parameters, a generator, or a saved RF-Track file. With no arguments the bunch
+        is empty. Each matrix row is one macroparticle.
+        The six standard columns are [x, xp, y, yp, t, Pc]: positions in mm, slopes
+        in mrad, arrival time in mm/c, and total momentum in MeV/c.
+        Extended rows add [mass, Q, N] and optionally ID; mass is in MeV/c^2,
+        Q is single-particle charge in units of e, and N is macroparticle population.
+        In the separate-species constructor, population is the total number of real
+        particles, not the number of macroparticles. S is in metres.
+        Twiss generation takes Pref in MeV/c and N macroparticles; sigma_cut=0
+        means no truncation. See Bunch6d_twiss for the distribution parameters.
+
+        Python example:
+            B = RF_Track.Bunch6d(RF_Track.electronmass, 1e6, -1,
+                                numpy.array([[0., 0., 0., 0., 0., 100.]]))
+            X = B.get_phase_space("%x %xp %y %yp %t %Pc")
+        Octave example (after loading RF_Track):
+            B = Bunch6d(electronmass, 1e6, -1, [0 0 0 0 0 100]);
+            X = B.get_phase_space('%x %xp %y %yp %t %Pc');
+
+
+        """
+        _RF_Track.Bunch6d_swiginit(self, _RF_Track.new_Bunch6d(*args))
+    resize = _swig_new_instance_method(_RF_Track.Bunch6d_resize)
+    clear = _swig_new_instance_method(_RF_Track.Bunch6d_clear)
+    size = _swig_new_instance_method(_RF_Track.Bunch6d_size)
+    set_lifetime = _swig_new_instance_method(_RF_Track.Bunch6d_set_lifetime)
+    set_polarization = _swig_new_instance_method(_RF_Track.Bunch6d_set_polarization)
+    displaced = _swig_new_instance_method(_RF_Track.Bunch6d_displaced)
+    append = _swig_new_instance_method(_RF_Track.Bunch6d_append)
+    coasting = _swig_new_instance_method(_RF_Track.Bunch6d_coasting)
+    set_coasting = _swig_new_instance_method(_RF_Track.Bunch6d_set_coasting)
+    set_phase_space = _swig_new_instance_method(_RF_Track.Bunch6d_set_phase_space)
+    set_S = _swig_new_instance_method(_RF_Track.Bunch6d_set_S)
+    get_particle = _swig_new_instance_method(_RF_Track.Bunch6d_get_particle)
+    get_t_mean = _swig_new_instance_method(_RF_Track.Bunch6d_get_t_mean)
+    get_t_min = _swig_new_instance_method(_RF_Track.Bunch6d_get_t_min)
+    get_t_max = _swig_new_instance_method(_RF_Track.Bunch6d_get_t_max)
+    get_S = _swig_new_instance_method(_RF_Track.Bunch6d_get_S)
+    get_average_particle = _swig_new_instance_method(_RF_Track.Bunch6d_get_average_particle)
+    get_reference_particle = _swig_new_instance_method(_RF_Track.Bunch6d_get_reference_particle)
+    get_bunch_temperature = _swig_new_instance_method(_RF_Track.Bunch6d_get_bunch_temperature)
+    get_phase_space = _swig_new_instance_method(_RF_Track.Bunch6d_get_phase_space)
+    get_population = _swig_new_instance_method(_RF_Track.Bunch6d_get_population)
+    get_ngood = _swig_new_instance_method(_RF_Track.Bunch6d_get_ngood)
+    get_nlost = _swig_new_instance_method(_RF_Track.Bunch6d_get_nlost)
+    get_total_charge = _swig_new_instance_method(_RF_Track.Bunch6d_get_total_charge)
+    get_info = _swig_new_instance_method(_RF_Track.Bunch6d_get_info)
+    apply_force = _swig_new_instance_method(_RF_Track.Bunch6d_apply_force)
+    kick = _swig_new_instance_method(_RF_Track.Bunch6d_kick)
+    drift = _swig_new_instance_method(_RF_Track.Bunch6d_drift)
+    load = _swig_new_instance_method(_RF_Track.Bunch6d_load)
+    save = _swig_new_instance_method(_RF_Track.Bunch6d_save)
+    save_as_dst_file = _swig_new_instance_method(_RF_Track.Bunch6d_save_as_dst_file)
+    save_as_sdds_file = _swig_new_instance_method(_RF_Track.Bunch6d_save_as_sdds_file)
+    __getitem__ = _swig_new_instance_method(_RF_Track.Bunch6d___getitem__)
+    __len__ = _swig_new_instance_method(_RF_Track.Bunch6d___len__)
+    __iter__ = _swig_new_instance_method(_RF_Track.Bunch6d___iter__)
+    __swig_destroy__ = _RF_Track.delete_Bunch6d
+
+# Register Bunch6d in _RF_Track:
+_RF_Track.Bunch6d_swigregister(Bunch6d)
+Bunch6d_QR = _RF_Track.Bunch6d_QR
+class Bunch6d_info(object):
+    r"""
+
+    Bunch statistics returned by get_info(). Contains population, mean coordinates, RMS
+    sizes, normalized emittances, Twiss parameters and correlations. Coordinate units follow
+    the corresponding bunch model; normalized emittances are in mm.mrad.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    S = property(_RF_Track.Bunch6d_info_S_get, _RF_Track.Bunch6d_info_S_set, doc=r"""S : double""")
+    mean_x = property(_RF_Track.Bunch6d_info_mean_x_get, _RF_Track.Bunch6d_info_mean_x_set, doc=r"""mean_x : double""")
+    mean_y = property(_RF_Track.Bunch6d_info_mean_y_get, _RF_Track.Bunch6d_info_mean_y_set, doc=r"""mean_y : double""")
+    mean_t = property(_RF_Track.Bunch6d_info_mean_t_get, _RF_Track.Bunch6d_info_mean_t_set, doc=r"""mean_t : double""")
+    mean_xp = property(_RF_Track.Bunch6d_info_mean_xp_get, _RF_Track.Bunch6d_info_mean_xp_set, doc=r"""mean_xp : double""")
+    mean_yp = property(_RF_Track.Bunch6d_info_mean_yp_get, _RF_Track.Bunch6d_info_mean_yp_set, doc=r"""mean_yp : double""")
+    mean_Px = property(_RF_Track.Bunch6d_info_mean_Px_get, _RF_Track.Bunch6d_info_mean_Px_set, doc=r"""mean_Px : double""")
+    mean_Py = property(_RF_Track.Bunch6d_info_mean_Py_get, _RF_Track.Bunch6d_info_mean_Py_set, doc=r"""mean_Py : double""")
+    mean_Pz = property(_RF_Track.Bunch6d_info_mean_Pz_get, _RF_Track.Bunch6d_info_mean_Pz_set, doc=r"""mean_Pz : double""")
+    mean_P = property(_RF_Track.Bunch6d_info_mean_P_get, _RF_Track.Bunch6d_info_mean_P_set, doc=r"""mean_P : double""")
+    mean_K = property(_RF_Track.Bunch6d_info_mean_K_get, _RF_Track.Bunch6d_info_mean_K_set, doc=r"""mean_K : double""")
+    mean_E = property(_RF_Track.Bunch6d_info_mean_E_get, _RF_Track.Bunch6d_info_mean_E_set, doc=r"""mean_E : double""")
+    mean_Sx = property(_RF_Track.Bunch6d_info_mean_Sx_get, _RF_Track.Bunch6d_info_mean_Sx_set, doc=r"""mean_Sx : double""")
+    mean_Sy = property(_RF_Track.Bunch6d_info_mean_Sy_get, _RF_Track.Bunch6d_info_mean_Sy_set, doc=r"""mean_Sy : double""")
+    mean_Sz = property(_RF_Track.Bunch6d_info_mean_Sz_get, _RF_Track.Bunch6d_info_mean_Sz_set, doc=r"""mean_Sz : double""")
+    sigma_x = property(_RF_Track.Bunch6d_info_sigma_x_get, _RF_Track.Bunch6d_info_sigma_x_set, doc=r"""sigma_x : double""")
+    sigma_y = property(_RF_Track.Bunch6d_info_sigma_y_get, _RF_Track.Bunch6d_info_sigma_y_set, doc=r"""sigma_y : double""")
+    sigma_t = property(_RF_Track.Bunch6d_info_sigma_t_get, _RF_Track.Bunch6d_info_sigma_t_set, doc=r"""sigma_t : double""")
+    sigma_px = property(_RF_Track.Bunch6d_info_sigma_px_get, _RF_Track.Bunch6d_info_sigma_px_set, doc=r"""sigma_px : double""")
+    sigma_py = property(_RF_Track.Bunch6d_info_sigma_py_get, _RF_Track.Bunch6d_info_sigma_py_set, doc=r"""sigma_py : double""")
+    sigma_pt = property(_RF_Track.Bunch6d_info_sigma_pt_get, _RF_Track.Bunch6d_info_sigma_pt_set, doc=r"""sigma_pt : double""")
+    sigma_E = property(_RF_Track.Bunch6d_info_sigma_E_get, _RF_Track.Bunch6d_info_sigma_E_set, doc=r"""sigma_E : double""")
+    sigma_P = property(_RF_Track.Bunch6d_info_sigma_P_get, _RF_Track.Bunch6d_info_sigma_P_set, doc=r"""sigma_P : double""")
+    emitt_x = property(_RF_Track.Bunch6d_info_emitt_x_get, _RF_Track.Bunch6d_info_emitt_x_set, doc=r"""emitt_x : double""")
+    emitt_y = property(_RF_Track.Bunch6d_info_emitt_y_get, _RF_Track.Bunch6d_info_emitt_y_set, doc=r"""emitt_y : double""")
+    emitt_z = property(_RF_Track.Bunch6d_info_emitt_z_get, _RF_Track.Bunch6d_info_emitt_z_set, doc=r"""emitt_z : double""")
+    emitt_4d = property(_RF_Track.Bunch6d_info_emitt_4d_get, _RF_Track.Bunch6d_info_emitt_4d_set, doc=r"""emitt_4d : double""")
+    emitt_6d = property(_RF_Track.Bunch6d_info_emitt_6d_get, _RF_Track.Bunch6d_info_emitt_6d_set, doc=r"""emitt_6d : double""")
+    alpha_x = property(_RF_Track.Bunch6d_info_alpha_x_get, _RF_Track.Bunch6d_info_alpha_x_set, doc=r"""alpha_x : double""")
+    alpha_y = property(_RF_Track.Bunch6d_info_alpha_y_get, _RF_Track.Bunch6d_info_alpha_y_set, doc=r"""alpha_y : double""")
+    alpha_z = property(_RF_Track.Bunch6d_info_alpha_z_get, _RF_Track.Bunch6d_info_alpha_z_set, doc=r"""alpha_z : double""")
+    beta_x = property(_RF_Track.Bunch6d_info_beta_x_get, _RF_Track.Bunch6d_info_beta_x_set, doc=r"""beta_x : double""")
+    beta_y = property(_RF_Track.Bunch6d_info_beta_y_get, _RF_Track.Bunch6d_info_beta_y_set, doc=r"""beta_y : double""")
+    beta_z = property(_RF_Track.Bunch6d_info_beta_z_get, _RF_Track.Bunch6d_info_beta_z_set, doc=r"""beta_z : double""")
+    disp_x = property(_RF_Track.Bunch6d_info_disp_x_get, _RF_Track.Bunch6d_info_disp_x_set, doc=r"""disp_x : double""")
+    disp_y = property(_RF_Track.Bunch6d_info_disp_y_get, _RF_Track.Bunch6d_info_disp_y_set, doc=r"""disp_y : double""")
+    disp_z = property(_RF_Track.Bunch6d_info_disp_z_get, _RF_Track.Bunch6d_info_disp_z_set, doc=r"""disp_z : double""")
+    disp_px = property(_RF_Track.Bunch6d_info_disp_px_get, _RF_Track.Bunch6d_info_disp_px_set, doc=r"""disp_px : double""")
+    disp_py = property(_RF_Track.Bunch6d_info_disp_py_get, _RF_Track.Bunch6d_info_disp_py_set, doc=r"""disp_py : double""")
+    rmax = property(_RF_Track.Bunch6d_info_rmax_get, _RF_Track.Bunch6d_info_rmax_set, doc=r"""rmax : double""")
+    rmax90 = property(_RF_Track.Bunch6d_info_rmax90_get, _RF_Track.Bunch6d_info_rmax90_set, doc=r"""rmax90 : double""")
+    rmax99 = property(_RF_Track.Bunch6d_info_rmax99_get, _RF_Track.Bunch6d_info_rmax99_set, doc=r"""rmax99 : double""")
+    rmax99p9 = property(_RF_Track.Bunch6d_info_rmax99p9_get, _RF_Track.Bunch6d_info_rmax99p9_set, doc=r"""rmax99p9 : double""")
+    xmax = property(_RF_Track.Bunch6d_info_xmax_get, _RF_Track.Bunch6d_info_xmax_set, doc=r"""xmax : double""")
+    xmax90 = property(_RF_Track.Bunch6d_info_xmax90_get, _RF_Track.Bunch6d_info_xmax90_set, doc=r"""xmax90 : double""")
+    xmax99 = property(_RF_Track.Bunch6d_info_xmax99_get, _RF_Track.Bunch6d_info_xmax99_set, doc=r"""xmax99 : double""")
+    ymax = property(_RF_Track.Bunch6d_info_ymax_get, _RF_Track.Bunch6d_info_ymax_set, doc=r"""ymax : double""")
+    ymax90 = property(_RF_Track.Bunch6d_info_ymax90_get, _RF_Track.Bunch6d_info_ymax90_set, doc=r"""ymax90 : double""")
+    ymax99 = property(_RF_Track.Bunch6d_info_ymax99_get, _RF_Track.Bunch6d_info_ymax99_set, doc=r"""ymax99 : double""")
+    transmission = property(_RF_Track.Bunch6d_info_transmission_get, _RF_Track.Bunch6d_info_transmission_set, doc=r"""transmission : double""")
+    get_bunch_info = _swig_new_instance_method(_RF_Track.Bunch6d_info_get_bunch_info)
+
+    def __init__(self):
+        r"""
+        __init__(Bunch6d_info self) -> Bunch6d_info
+
+        Bunch statistics returned by get_info(). Contains population, mean coordinates, RMS
+        sizes, normalized emittances, Twiss parameters and correlations. Coordinate units follow
+        the corresponding bunch model; normalized emittances are in mm.mrad.
+
+
+        """
+        _RF_Track.Bunch6d_info_swiginit(self, _RF_Track.new_Bunch6d_info())
+    __swig_destroy__ = _RF_Track.delete_Bunch6d_info
+
+# Register Bunch6d_info in _RF_Track:
+_RF_Track.Bunch6d_info_swigregister(Bunch6d_info)
+class Bunch6dT_info(object):
+    r"""
+
+    Bunch statistics returned by get_info(). Contains population, mean coordinates, RMS
+    sizes, normalized emittances, Twiss parameters and correlations. Coordinate units follow
+    the corresponding bunch model; normalized emittances are in mm.mrad.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    t = property(_RF_Track.Bunch6dT_info_t_get, _RF_Track.Bunch6dT_info_t_set, doc=r"""t : double""")
+    mean_X = property(_RF_Track.Bunch6dT_info_mean_X_get, _RF_Track.Bunch6dT_info_mean_X_set, doc=r"""mean_X : double""")
+    mean_Y = property(_RF_Track.Bunch6dT_info_mean_Y_get, _RF_Track.Bunch6dT_info_mean_Y_set, doc=r"""mean_Y : double""")
+    mean_Z = property(_RF_Track.Bunch6dT_info_mean_Z_get, _RF_Track.Bunch6dT_info_mean_Z_set, doc=r"""mean_Z : double""")
+    mean_Px = property(_RF_Track.Bunch6dT_info_mean_Px_get, _RF_Track.Bunch6dT_info_mean_Px_set, doc=r"""mean_Px : double""")
+    mean_Py = property(_RF_Track.Bunch6dT_info_mean_Py_get, _RF_Track.Bunch6dT_info_mean_Py_set, doc=r"""mean_Py : double""")
+    mean_Pz = property(_RF_Track.Bunch6dT_info_mean_Pz_get, _RF_Track.Bunch6dT_info_mean_Pz_set, doc=r"""mean_Pz : double""")
+    mean_K = property(_RF_Track.Bunch6dT_info_mean_K_get, _RF_Track.Bunch6dT_info_mean_K_set, doc=r"""mean_K : double""")
+    mean_E = property(_RF_Track.Bunch6dT_info_mean_E_get, _RF_Track.Bunch6dT_info_mean_E_set, doc=r"""mean_E : double""")
+    mean_Sx = property(_RF_Track.Bunch6dT_info_mean_Sx_get, _RF_Track.Bunch6dT_info_mean_Sx_set, doc=r"""mean_Sx : double""")
+    mean_Sy = property(_RF_Track.Bunch6dT_info_mean_Sy_get, _RF_Track.Bunch6dT_info_mean_Sy_set, doc=r"""mean_Sy : double""")
+    mean_Sz = property(_RF_Track.Bunch6dT_info_mean_Sz_get, _RF_Track.Bunch6dT_info_mean_Sz_set, doc=r"""mean_Sz : double""")
+    sigma_X = property(_RF_Track.Bunch6dT_info_sigma_X_get, _RF_Track.Bunch6dT_info_sigma_X_set, doc=r"""sigma_X : double""")
+    sigma_Y = property(_RF_Track.Bunch6dT_info_sigma_Y_get, _RF_Track.Bunch6dT_info_sigma_Y_set, doc=r"""sigma_Y : double""")
+    sigma_Z = property(_RF_Track.Bunch6dT_info_sigma_Z_get, _RF_Track.Bunch6dT_info_sigma_Z_set, doc=r"""sigma_Z : double""")
+    sigma_Px = property(_RF_Track.Bunch6dT_info_sigma_Px_get, _RF_Track.Bunch6dT_info_sigma_Px_set, doc=r"""sigma_Px : double""")
+    sigma_Py = property(_RF_Track.Bunch6dT_info_sigma_Py_get, _RF_Track.Bunch6dT_info_sigma_Py_set, doc=r"""sigma_Py : double""")
+    sigma_Pz = property(_RF_Track.Bunch6dT_info_sigma_Pz_get, _RF_Track.Bunch6dT_info_sigma_Pz_set, doc=r"""sigma_Pz : double""")
+    sigma_E = property(_RF_Track.Bunch6dT_info_sigma_E_get, _RF_Track.Bunch6dT_info_sigma_E_set, doc=r"""sigma_E : double""")
+    emitt_x = property(_RF_Track.Bunch6dT_info_emitt_x_get, _RF_Track.Bunch6dT_info_emitt_x_set, doc=r"""emitt_x : double""")
+    emitt_y = property(_RF_Track.Bunch6dT_info_emitt_y_get, _RF_Track.Bunch6dT_info_emitt_y_set, doc=r"""emitt_y : double""")
+    emitt_z = property(_RF_Track.Bunch6dT_info_emitt_z_get, _RF_Track.Bunch6dT_info_emitt_z_set, doc=r"""emitt_z : double""")
+    emitt_4d = property(_RF_Track.Bunch6dT_info_emitt_4d_get, _RF_Track.Bunch6dT_info_emitt_4d_set, doc=r"""emitt_4d : double""")
+    emitt_6d = property(_RF_Track.Bunch6dT_info_emitt_6d_get, _RF_Track.Bunch6dT_info_emitt_6d_set, doc=r"""emitt_6d : double""")
+    alpha_x = property(_RF_Track.Bunch6dT_info_alpha_x_get, _RF_Track.Bunch6dT_info_alpha_x_set, doc=r"""alpha_x : double""")
+    alpha_y = property(_RF_Track.Bunch6dT_info_alpha_y_get, _RF_Track.Bunch6dT_info_alpha_y_set, doc=r"""alpha_y : double""")
+    alpha_z = property(_RF_Track.Bunch6dT_info_alpha_z_get, _RF_Track.Bunch6dT_info_alpha_z_set, doc=r"""alpha_z : double""")
+    beta_x = property(_RF_Track.Bunch6dT_info_beta_x_get, _RF_Track.Bunch6dT_info_beta_x_set, doc=r"""beta_x : double""")
+    beta_y = property(_RF_Track.Bunch6dT_info_beta_y_get, _RF_Track.Bunch6dT_info_beta_y_set, doc=r"""beta_y : double""")
+    beta_z = property(_RF_Track.Bunch6dT_info_beta_z_get, _RF_Track.Bunch6dT_info_beta_z_set, doc=r"""beta_z : double""")
+    disp_x = property(_RF_Track.Bunch6dT_info_disp_x_get, _RF_Track.Bunch6dT_info_disp_x_set, doc=r"""disp_x : double""")
+    disp_y = property(_RF_Track.Bunch6dT_info_disp_y_get, _RF_Track.Bunch6dT_info_disp_y_set, doc=r"""disp_y : double""")
+    disp_z = property(_RF_Track.Bunch6dT_info_disp_z_get, _RF_Track.Bunch6dT_info_disp_z_set, doc=r"""disp_z : double""")
+    disp_px = property(_RF_Track.Bunch6dT_info_disp_px_get, _RF_Track.Bunch6dT_info_disp_px_set, doc=r"""disp_px : double""")
+    disp_py = property(_RF_Track.Bunch6dT_info_disp_py_get, _RF_Track.Bunch6dT_info_disp_py_set, doc=r"""disp_py : double""")
+    rmax = property(_RF_Track.Bunch6dT_info_rmax_get, _RF_Track.Bunch6dT_info_rmax_set, doc=r"""rmax : double""")
+    rmax90 = property(_RF_Track.Bunch6dT_info_rmax90_get, _RF_Track.Bunch6dT_info_rmax90_set, doc=r"""rmax90 : double""")
+    rmax99 = property(_RF_Track.Bunch6dT_info_rmax99_get, _RF_Track.Bunch6dT_info_rmax99_set, doc=r"""rmax99 : double""")
+    rmax99p9 = property(_RF_Track.Bunch6dT_info_rmax99p9_get, _RF_Track.Bunch6dT_info_rmax99p9_set, doc=r"""rmax99p9 : double""")
+    xmax = property(_RF_Track.Bunch6dT_info_xmax_get, _RF_Track.Bunch6dT_info_xmax_set, doc=r"""xmax : double""")
+    xmax90 = property(_RF_Track.Bunch6dT_info_xmax90_get, _RF_Track.Bunch6dT_info_xmax90_set, doc=r"""xmax90 : double""")
+    xmax99 = property(_RF_Track.Bunch6dT_info_xmax99_get, _RF_Track.Bunch6dT_info_xmax99_set, doc=r"""xmax99 : double""")
+    ymax = property(_RF_Track.Bunch6dT_info_ymax_get, _RF_Track.Bunch6dT_info_ymax_set, doc=r"""ymax : double""")
+    ymax90 = property(_RF_Track.Bunch6dT_info_ymax90_get, _RF_Track.Bunch6dT_info_ymax90_set, doc=r"""ymax90 : double""")
+    ymax99 = property(_RF_Track.Bunch6dT_info_ymax99_get, _RF_Track.Bunch6dT_info_ymax99_set, doc=r"""ymax99 : double""")
+    transmission = property(_RF_Track.Bunch6dT_info_transmission_get, _RF_Track.Bunch6dT_info_transmission_set, doc=r"""transmission : double""")
+    get_bunch_info = _swig_new_instance_method(_RF_Track.Bunch6dT_info_get_bunch_info)
+
+    def __init__(self):
+        r"""
+        __init__(Bunch6dT_info self) -> Bunch6dT_info
+
+        Bunch statistics returned by get_info(). Contains population, mean coordinates, RMS
+        sizes, normalized emittances, Twiss parameters and correlations. Coordinate units follow
+        the corresponding bunch model; normalized emittances are in mm.mrad.
+
+
+        """
+        _RF_Track.Bunch6dT_info_swiginit(self, _RF_Track.new_Bunch6dT_info())
+    __swig_destroy__ = _RF_Track.delete_Bunch6dT_info
+
+# Register Bunch6dT_info in _RF_Track:
+_RF_Track.Bunch6dT_info_swigregister(Bunch6dT_info)
+class Bunch6d_twiss(object):
+    r"""
+
+    Twiss parameters for bunch generation. Normalized emittances are in mm.mrad; beta
+    functions and spatial dispersions are in metres. See the individual fields for
+    longitudinal conventions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    emitt_x = property(_RF_Track.Bunch6d_twiss_emitt_x_get, _RF_Track.Bunch6d_twiss_emitt_x_set, doc=r"""
+    emitt_x : double
+
+    Normalized x emittance in mm.mrad. Default: 0.
+
+
+    """)
+    emitt_y = property(_RF_Track.Bunch6d_twiss_emitt_y_get, _RF_Track.Bunch6d_twiss_emitt_y_set, doc=r"""
+    emitt_y : double
+
+    Normalized y emittance in mm.mrad. Default: 0.
+
+
+    """)
+    emitt_z = property(_RF_Track.Bunch6d_twiss_emitt_z_get, _RF_Track.Bunch6d_twiss_emitt_z_set, doc=r"""
+    emitt_z : double
+
+    Normalized z emittance in mm.mrad. Default: 0.
+
+
+    """)
+    alpha_x = property(_RF_Track.Bunch6d_twiss_alpha_x_get, _RF_Track.Bunch6d_twiss_alpha_x_set, doc=r"""
+    alpha_x : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    alpha_y = property(_RF_Track.Bunch6d_twiss_alpha_y_get, _RF_Track.Bunch6d_twiss_alpha_y_set, doc=r"""
+    alpha_y : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    alpha_z = property(_RF_Track.Bunch6d_twiss_alpha_z_get, _RF_Track.Bunch6d_twiss_alpha_z_set, doc=r"""
+    alpha_z : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    beta_x = property(_RF_Track.Bunch6d_twiss_beta_x_get, _RF_Track.Bunch6d_twiss_beta_x_set, doc=r"""
+    beta_x : double
+
+    Twiss beta_x in metres. Default: 0.
+
+
+    """)
+    beta_y = property(_RF_Track.Bunch6d_twiss_beta_y_get, _RF_Track.Bunch6d_twiss_beta_y_set, doc=r"""
+    beta_y : double
+
+    Twiss beta_y in metres. Default: 0.
+
+
+    """)
+    beta_z = property(_RF_Track.Bunch6d_twiss_beta_z_get, _RF_Track.Bunch6d_twiss_beta_z_set, doc=r"""
+    beta_z : double
+
+    Twiss beta_z in metres. Default: 0.
+
+
+    """)
+    sigma_t = property(_RF_Track.Bunch6d_twiss_sigma_t_get, _RF_Track.Bunch6d_twiss_sigma_t_set, doc=r"""
+    sigma_t : double
+
+    RMS arrival-time spread in mm/c. Default: 0.
+
+
+    """)
+    sigma_pt = property(_RF_Track.Bunch6d_twiss_sigma_pt_get, _RF_Track.Bunch6d_twiss_sigma_pt_set, doc=r"""
+    sigma_pt : double
+
+    RMS relative energy coordinate (E-E0)/(P0*c), expressed in per mille. Default: 0.
+
+
+    """)
+    disp_x = property(_RF_Track.Bunch6d_twiss_disp_x_get, _RF_Track.Bunch6d_twiss_disp_x_set, doc=r"""
+    disp_x : double
+
+    Spatial x dispersion in metres. Default: 0.
+
+
+    """)
+    disp_px = property(_RF_Track.Bunch6d_twiss_disp_px_get, _RF_Track.Bunch6d_twiss_disp_px_set, doc=r"""
+    disp_px : double
+
+    Angular x dispersion in radians. Default: 0.
+
+
+    """)
+    disp_y = property(_RF_Track.Bunch6d_twiss_disp_y_get, _RF_Track.Bunch6d_twiss_disp_y_set, doc=r"""
+    disp_y : double
+
+    Spatial y dispersion in metres. Default: 0.
+
+
+    """)
+    disp_py = property(_RF_Track.Bunch6d_twiss_disp_py_get, _RF_Track.Bunch6d_twiss_disp_py_set, doc=r"""
+    disp_py : double
+
+    Angular y dispersion in radians. Default: 0.
+
+
+    """)
+    disp_z = property(_RF_Track.Bunch6d_twiss_disp_z_get, _RF_Track.Bunch6d_twiss_disp_z_set, doc=r"""
+    disp_z : double
+
+    Spatial z dispersion in metres. Default: 0.
+
+
+    """)
+
+    def __init__(self):
+        r"""
+        __init__(Bunch6d_twiss self) -> Bunch6d_twiss
+
+        Twiss parameters for bunch generation. Normalized emittances are in mm.mrad; beta
+        functions and spatial dispersions are in metres. See the individual fields for
+        longitudinal conventions.
+
+
+        """
+        _RF_Track.Bunch6d_twiss_swiginit(self, _RF_Track.new_Bunch6d_twiss())
+    __swig_destroy__ = _RF_Track.delete_Bunch6d_twiss
+
+# Register Bunch6d_twiss in _RF_Track:
+_RF_Track.Bunch6d_twiss_swigregister(Bunch6d_twiss)
+class Bunch6dT_twiss(object):
+    r"""
+
+    Twiss parameters for bunch generation. Normalized emittances are in mm.mrad; beta
+    functions and spatial dispersions are in metres. See the individual fields for
+    longitudinal conventions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    emitt_x = property(_RF_Track.Bunch6dT_twiss_emitt_x_get, _RF_Track.Bunch6dT_twiss_emitt_x_set, doc=r"""
+    emitt_x : double
+
+    Normalized x emittance in mm.mrad. Default: 0.
+
+
+    """)
+    emitt_y = property(_RF_Track.Bunch6dT_twiss_emitt_y_get, _RF_Track.Bunch6dT_twiss_emitt_y_set, doc=r"""
+    emitt_y : double
+
+    Normalized y emittance in mm.mrad. Default: 0.
+
+
+    """)
+    emitt_z = property(_RF_Track.Bunch6dT_twiss_emitt_z_get, _RF_Track.Bunch6dT_twiss_emitt_z_set, doc=r"""
+    emitt_z : double
+
+    Normalized z emittance in mm.mrad. Default: 0.
+
+
+    """)
+    alpha_x = property(_RF_Track.Bunch6dT_twiss_alpha_x_get, _RF_Track.Bunch6dT_twiss_alpha_x_set, doc=r"""
+    alpha_x : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    alpha_y = property(_RF_Track.Bunch6dT_twiss_alpha_y_get, _RF_Track.Bunch6dT_twiss_alpha_y_set, doc=r"""
+    alpha_y : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    alpha_z = property(_RF_Track.Bunch6dT_twiss_alpha_z_get, _RF_Track.Bunch6dT_twiss_alpha_z_set, doc=r"""
+    alpha_z : double
+
+    Dimensionless Twiss alpha parameter. Default: 0.0.
+
+
+    """)
+    beta_x = property(_RF_Track.Bunch6dT_twiss_beta_x_get, _RF_Track.Bunch6dT_twiss_beta_x_set, doc=r"""
+    beta_x : double
+
+    Twiss beta_x in metres. Default: 0.
+
+
+    """)
+    beta_y = property(_RF_Track.Bunch6dT_twiss_beta_y_get, _RF_Track.Bunch6dT_twiss_beta_y_set, doc=r"""
+    beta_y : double
+
+    Twiss beta_y in metres. Default: 0.
+
+
+    """)
+    beta_z = property(_RF_Track.Bunch6dT_twiss_beta_z_get, _RF_Track.Bunch6dT_twiss_beta_z_set, doc=r"""
+    beta_z : double
+
+    Twiss beta_z in metres. Default: 0.
+
+
+    """)
+    sigma_z = property(_RF_Track.Bunch6dT_twiss_sigma_z_get, _RF_Track.Bunch6dT_twiss_sigma_z_set, doc=r"""
+    sigma_z : double
+
+    RMS bunch length in mm. Default: 0.
+
+
+    """)
+    sigma_pz = property(_RF_Track.Bunch6dT_twiss_sigma_pz_get, _RF_Track.Bunch6dT_twiss_sigma_pz_set, doc=r"""
+    sigma_pz : double
+
+    RMS relative longitudinal momentum spread (Pz-P0)/P0, expressed in per mille. Default:
+    0.
+
+
+    """)
+    disp_x = property(_RF_Track.Bunch6dT_twiss_disp_x_get, _RF_Track.Bunch6dT_twiss_disp_x_set, doc=r"""
+    disp_x : double
+
+    Spatial x dispersion in metres. Default: 0.
+
+
+    """)
+    disp_px = property(_RF_Track.Bunch6dT_twiss_disp_px_get, _RF_Track.Bunch6dT_twiss_disp_px_set, doc=r"""
+    disp_px : double
+
+    Angular x dispersion in radians. Default: 0.
+
+
+    """)
+    disp_y = property(_RF_Track.Bunch6dT_twiss_disp_y_get, _RF_Track.Bunch6dT_twiss_disp_y_set, doc=r"""
+    disp_y : double
+
+    Spatial y dispersion in metres. Default: 0.
+
+
+    """)
+    disp_py = property(_RF_Track.Bunch6dT_twiss_disp_py_get, _RF_Track.Bunch6dT_twiss_disp_py_set, doc=r"""
+    disp_py : double
+
+    Angular y dispersion in radians. Default: 0.
+
+
+    """)
+    disp_z = property(_RF_Track.Bunch6dT_twiss_disp_z_get, _RF_Track.Bunch6dT_twiss_disp_z_set, doc=r"""
+    disp_z : double
+
+    Spatial z dispersion in metres. Default: 0.
+
+
+    """)
+
+    def __init__(self):
+        r"""
+        __init__(Bunch6dT_twiss self) -> Bunch6dT_twiss
+
+        Twiss parameters for bunch generation. Normalized emittances are in mm.mrad; beta
+        functions and spatial dispersions are in metres. See the individual fields for
+        longitudinal conventions.
+
+
+        """
+        _RF_Track.Bunch6dT_twiss_swiginit(self, _RF_Track.new_Bunch6dT_twiss())
+    __swig_destroy__ = _RF_Track.delete_Bunch6dT_twiss
+
+# Register Bunch6dT_twiss in _RF_Track:
+_RF_Track.Bunch6dT_twiss_swigregister(Bunch6dT_twiss)
+class Bunch6dT_Generator(object):
+    r"""
+
+    Distribution and photocathode-emission parameters for bunch generation.
+    Pass this object and the requested macroparticle count to Bunch6dT.
+    The generator uses mixed units: q_total [nC], ref_zpos [m], ref_clock [ns],
+    ref_ekin [MeV], transverse sizes [mm], momentum spreads [keV/c] and emission
+    times [ns]. These differ from the bunch storage units. cathode enables
+    emission-time generation. See individual property help for units.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    species = property(_RF_Track.Bunch6dT_Generator_species_get, _RF_Track.Bunch6dT_Generator_species_set, doc=r"""
+    species : std::string
+
+    Type of particle. Default: "electron".
+
+
+    """)
+    noise_reduc = property(_RF_Track.Bunch6dT_Generator_noise_reduc_get, _RF_Track.Bunch6dT_Generator_noise_reduc_set, doc=r"""
+    noise_reduc : bool
+
+    Enable noise reduction. Default: true.
+
+
+    """)
+    rand_generator = property(_RF_Track.Bunch6dT_Generator_rand_generator_get, _RF_Track.Bunch6dT_Generator_rand_generator_set, doc=r"""
+    rand_generator : std::string
+
+    generator type. Default: "halton".
+
+
+    """)
+    cathode = property(_RF_Track.Bunch6dT_Generator_cathode_get, _RF_Track.Bunch6dT_Generator_cathode_set, doc=r"""
+    cathode : bool
+
+    Enable photocathode emission. Default: false.
+
+
+    """)
+    q_total = property(_RF_Track.Bunch6dT_Generator_q_total_get, _RF_Track.Bunch6dT_Generator_q_total_set, doc=r"""
+    q_total : double
+
+    nC bunch charge. Default: 0.
+
+
+    """)
+    ref_zpos = property(_RF_Track.Bunch6dT_Generator_ref_zpos_get, _RF_Track.Bunch6dT_Generator_ref_zpos_set, doc=r"""
+    ref_zpos : double
+
+    m position of reference particle. Default: 0.
+
+
+    """)
+    ref_clock = property(_RF_Track.Bunch6dT_Generator_ref_clock_get, _RF_Track.Bunch6dT_Generator_ref_clock_set, doc=r"""
+    ref_clock : double
+
+    ns clock of reference particle. Default: 0.
+
+
+    """)
+    ref_ekin = property(_RF_Track.Bunch6dT_Generator_ref_ekin_get, _RF_Track.Bunch6dT_Generator_ref_ekin_set, doc=r"""
+    ref_ekin : double
+
+    MeV energy of reference particle. Default: 0.
+
+
+    """)
+    dist_z = property(_RF_Track.Bunch6dT_Generator_dist_z_get, _RF_Track.Bunch6dT_Generator_dist_z_set, doc=r"""
+    dist_z : std::string
+
+    longitudinal particle distribution. Default: "gaussian".
+
+
+    """)
+    sig_z = property(_RF_Track.Bunch6dT_Generator_sig_z_get, _RF_Track.Bunch6dT_Generator_sig_z_set, doc=r"""
+    sig_z : double
+
+    mm rms value of the bunch length. Default: 0.
+
+
+    """)
+    shape_z = property(_RF_Track.Bunch6dT_Generator_shape_z_get, _RF_Track.Bunch6dT_Generator_shape_z_set, doc=r"""
+    shape_z : double
+
+    shape of the Generalized Gaussian distribution. Default: 2.
+
+
+    """)
+    c_sig_z = property(_RF_Track.Bunch6dT_Generator_c_sig_z_get, _RF_Track.Bunch6dT_Generator_c_sig_z_set, doc=r"""
+    c_sig_z : double
+
+    cuts off a Gaussian longitudinal distribution. Default: 4.
+
+
+    """)
+    lz = property(_RF_Track.Bunch6dT_Generator_lz_get, _RF_Track.Bunch6dT_Generator_lz_set, doc=r"""
+    lz : double
+
+    mm length of the bunch. Default: 0.
+
+
+    """)
+    rz = property(_RF_Track.Bunch6dT_Generator_rz_get, _RF_Track.Bunch6dT_Generator_rz_set, doc=r"""
+    rz : double
+
+    mm rising of the bunch distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    sig_t = property(_RF_Track.Bunch6dT_Generator_sig_t_get, _RF_Track.Bunch6dT_Generator_sig_t_set, doc=r"""
+    sig_t : double
+
+    ns rms value of the emission time,. Default: 0.
+
+
+    """)
+    shape_t = property(_RF_Track.Bunch6dT_Generator_shape_t_get, _RF_Track.Bunch6dT_Generator_shape_t_set, doc=r"""
+    shape_t : double
+
+    shape of the Generalized Gaussian distribution. Default: 2.
+
+
+    """)
+    c_sig_t = property(_RF_Track.Bunch6dT_Generator_c_sig_t_get, _RF_Track.Bunch6dT_Generator_c_sig_t_set, doc=r"""
+    c_sig_t : double
+
+    cuts off a Gaussian longitudinal distribution. Default: 4.
+
+
+    """)
+    lt = property(_RF_Track.Bunch6dT_Generator_lt_get, _RF_Track.Bunch6dT_Generator_lt_set, doc=r"""
+    lt : double
+
+    ns emission time. Default: 0.
+
+
+    """)
+    rt = property(_RF_Track.Bunch6dT_Generator_rt_get, _RF_Track.Bunch6dT_Generator_rt_set, doc=r"""
+    rt : double
+
+    ns rise time of the bunch; only for plateau distribution. Default: 0.
+
+
+    """)
+    dist_pz = property(_RF_Track.Bunch6dT_Generator_dist_pz_get, _RF_Track.Bunch6dT_Generator_dist_pz_set, doc=r"""
+    dist_pz : std::string
+
+    longitudinal energy and momentum distribution. Default: "gaussian".
+
+
+    """)
+    sig_ekin = property(_RF_Track.Bunch6dT_Generator_sig_ekin_get, _RF_Track.Bunch6dT_Generator_sig_ekin_set, doc=r"""
+    sig_ekin : double
+
+    keV rms value of the energy spread. Default: 0.
+
+
+    """)
+    c_sig_ekin = property(_RF_Track.Bunch6dT_Generator_c_sig_ekin_get, _RF_Track.Bunch6dT_Generator_c_sig_ekin_set, doc=r"""
+    c_sig_ekin : double
+
+    cuts off a Gaussian energy and momentum distribution. Default: 4.
+
+
+    """)
+    le = property(_RF_Track.Bunch6dT_Generator_le_get, _RF_Track.Bunch6dT_Generator_le_set, doc=r"""
+    le : double
+
+    keV width of the energy distribution. Default: 0.
+
+
+    """)
+    re = property(_RF_Track.Bunch6dT_Generator_re_get, _RF_Track.Bunch6dT_Generator_re_set, doc=r"""
+    re : double
+
+    keV rising of the energy distribution, only for plateau distribution. Default: 0.
+
+
+    """)
+    emit_z = property(_RF_Track.Bunch6dT_Generator_emit_z_get, _RF_Track.Bunch6dT_Generator_emit_z_set, doc=r"""
+    emit_z : double
+
+    keV mm longitudinal emittance. Default: 0.
+
+
+    """)
+    emit_t = property(_RF_Track.Bunch6dT_Generator_emit_t_get, _RF_Track.Bunch6dT_Generator_emit_t_set, doc=r"""
+    emit_t : double
+
+    keV ps longitudinal emittance. Default: 0.
+
+
+    """)
+    cor_ekin = property(_RF_Track.Bunch6dT_Generator_cor_ekin_get, _RF_Track.Bunch6dT_Generator_cor_ekin_set, doc=r"""
+    cor_ekin : double
+
+    keV correlated energy spread. Default: 0.
+
+
+    """)
+    e_photon = property(_RF_Track.Bunch6dT_Generator_e_photon_get, _RF_Track.Bunch6dT_Generator_e_photon_set, doc=r"""
+    e_photon : double
+
+    eV photon energy for Fermi-Dirac distribution. Default: 0.
+
+
+    """)
+    phi_eff = property(_RF_Track.Bunch6dT_Generator_phi_eff_get, _RF_Track.Bunch6dT_Generator_phi_eff_set, doc=r"""
+    phi_eff : double
+
+    eV effective work function for Fermi-Dirac distribution. Default: 0.
+
+
+    """)
+    dist_x = property(_RF_Track.Bunch6dT_Generator_dist_x_get, _RF_Track.Bunch6dT_Generator_dist_x_set, doc=r"""
+    dist_x : std::string
+
+    transverse particle distribution. Default: "gaussian".
+
+
+    """)
+    sig_x = property(_RF_Track.Bunch6dT_Generator_sig_x_get, _RF_Track.Bunch6dT_Generator_sig_x_set, doc=r"""
+    sig_x : double
+
+    mm rms bunch size in the horizontal direction. Default: 0.
+
+
+    """)
+    shape_x = property(_RF_Track.Bunch6dT_Generator_shape_x_get, _RF_Track.Bunch6dT_Generator_shape_x_set, doc=r"""
+    shape_x : double
+
+    shape of the Generalized Gaussian distribution  . Default: 2.
+
+
+    """)
+    c_sig_x = property(_RF_Track.Bunch6dT_Generator_c_sig_x_get, _RF_Track.Bunch6dT_Generator_c_sig_x_set, doc=r"""
+    c_sig_x : double
+
+    cuts off a Gaussian horizontal distribution. Default: 4.
+
+
+    """)
+    lx = property(_RF_Track.Bunch6dT_Generator_lx_get, _RF_Track.Bunch6dT_Generator_lx_set, doc=r"""
+    lx : double
+
+    mm width of the horizontal particle distribution. Default: 0.
+
+
+    """)
+    rx = property(_RF_Track.Bunch6dT_Generator_rx_get, _RF_Track.Bunch6dT_Generator_rx_set, doc=r"""
+    rx : double
+
+    mm rising of the horizontal particle distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    x_off = property(_RF_Track.Bunch6dT_Generator_x_off_get, _RF_Track.Bunch6dT_Generator_x_off_set, doc=r"""
+    x_off : double
+
+    mm horizontal offset of the particle distribution. Default: 0.
+
+
+    """)
+    disp_x = property(_RF_Track.Bunch6dT_Generator_disp_x_get, _RF_Track.Bunch6dT_Generator_disp_x_set, doc=r"""
+    disp_x : double
+
+    m horizontal dispersion;. Default: 0.
+
+
+    """)
+    dist_y = property(_RF_Track.Bunch6dT_Generator_dist_y_get, _RF_Track.Bunch6dT_Generator_dist_y_set, doc=r"""
+    dist_y : std::string
+
+    vertical particle distribution. Default: "gaussian".
+
+
+    """)
+    sig_y = property(_RF_Track.Bunch6dT_Generator_sig_y_get, _RF_Track.Bunch6dT_Generator_sig_y_set, doc=r"""
+    sig_y : double
+
+    mm rms bunch size in the vertical direction. Default: 0.
+
+
+    """)
+    shape_y = property(_RF_Track.Bunch6dT_Generator_shape_y_get, _RF_Track.Bunch6dT_Generator_shape_y_set, doc=r"""
+    shape_y : double
+
+    shape of the Generalized Gaussian distribution  . Default: 2.
+
+
+    """)
+    c_sig_y = property(_RF_Track.Bunch6dT_Generator_c_sig_y_get, _RF_Track.Bunch6dT_Generator_c_sig_y_set, doc=r"""
+    c_sig_y : double
+
+    cuts off a Gaussian vertical distribution. Default: 4.
+
+
+    """)
+    ly = property(_RF_Track.Bunch6dT_Generator_ly_get, _RF_Track.Bunch6dT_Generator_ly_set, doc=r"""
+    ly : double
+
+    mm width of the vertical particle distribution. Default: 0.
+
+
+    """)
+    ry = property(_RF_Track.Bunch6dT_Generator_ry_get, _RF_Track.Bunch6dT_Generator_ry_set, doc=r"""
+    ry : double
+
+    mm rising of the vertical particle distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    y_off = property(_RF_Track.Bunch6dT_Generator_y_off_get, _RF_Track.Bunch6dT_Generator_y_off_set, doc=r"""
+    y_off : double
+
+    mm vertical offset of the particle distribution. Default: 0.
+
+
+    """)
+    disp_y = property(_RF_Track.Bunch6dT_Generator_disp_y_get, _RF_Track.Bunch6dT_Generator_disp_y_set, doc=r"""
+    disp_y : double
+
+    m vertical dispersion;. Default: 0.
+
+
+    """)
+    dist_px = property(_RF_Track.Bunch6dT_Generator_dist_px_get, _RF_Track.Bunch6dT_Generator_dist_px_set, doc=r"""
+    dist_px : std::string
+
+    transverse momentum distribution in the horizontal direction. Default: "gaussian".
+
+
+    """)
+    nemit_x = property(_RF_Track.Bunch6dT_Generator_nemit_x_get, _RF_Track.Bunch6dT_Generator_nemit_x_set, doc=r"""
+    nemit_x : double
+
+    mrad mm normalized transverse emittance in the horizontal direction. Default: 0.
+
+
+    """)
+    sig_px = property(_RF_Track.Bunch6dT_Generator_sig_px_get, _RF_Track.Bunch6dT_Generator_sig_px_set, doc=r"""
+    sig_px : double
+
+    keV/c rms value of the horizontal momentum distribution. Default: 0.
+
+
+    """)
+    c_sig_px = property(_RF_Track.Bunch6dT_Generator_c_sig_px_get, _RF_Track.Bunch6dT_Generator_c_sig_px_set, doc=r"""
+    c_sig_px : double
+
+    cuts off the horizontal momentum distribution. Default: 4.
+
+
+    """)
+    lpx = property(_RF_Track.Bunch6dT_Generator_lpx_get, _RF_Track.Bunch6dT_Generator_lpx_set, doc=r"""
+    lpx : double
+
+    keV/c width of the horizontal momentum distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    rpx = property(_RF_Track.Bunch6dT_Generator_rpx_get, _RF_Track.Bunch6dT_Generator_rpx_set, doc=r"""
+    rpx : double
+
+    keV/c rising of the horizontal momentum distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    cor_px = property(_RF_Track.Bunch6dT_Generator_cor_px_get, _RF_Track.Bunch6dT_Generator_cor_px_set, doc=r"""
+    cor_px : double
+
+    mrad correlated beam divergence in the horizontal direction. Default: 0.
+
+
+    """)
+    dist_py = property(_RF_Track.Bunch6dT_Generator_dist_py_get, _RF_Track.Bunch6dT_Generator_dist_py_set, doc=r"""
+    dist_py : std::string
+
+    transverse momentum distribution in the vertical direction. Default: "gaussian".
+
+
+    """)
+    nemit_y = property(_RF_Track.Bunch6dT_Generator_nemit_y_get, _RF_Track.Bunch6dT_Generator_nemit_y_set, doc=r"""
+    nemit_y : double
+
+    mrad mm normalized transverse emittance in the vertical direction. Default: 0.
+
+
+    """)
+    sig_py = property(_RF_Track.Bunch6dT_Generator_sig_py_get, _RF_Track.Bunch6dT_Generator_sig_py_set, doc=r"""
+    sig_py : double
+
+    keV/c rms value of the vertical momentum distribution. Default: 0.
+
+
+    """)
+    c_sig_py = property(_RF_Track.Bunch6dT_Generator_c_sig_py_get, _RF_Track.Bunch6dT_Generator_c_sig_py_set, doc=r"""
+    c_sig_py : double
+
+    cuts off the vertical momentum distribution. Default: 4.
+
+
+    """)
+    lpy = property(_RF_Track.Bunch6dT_Generator_lpy_get, _RF_Track.Bunch6dT_Generator_lpy_set, doc=r"""
+    lpy : double
+
+    keV/c width of the vertical momentum distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    rpy = property(_RF_Track.Bunch6dT_Generator_rpy_get, _RF_Track.Bunch6dT_Generator_rpy_set, doc=r"""
+    rpy : double
+
+    keV/c rising of the vertical momentum distribution; only for plateau distribution. Default: 0.
+
+
+    """)
+    cor_py = property(_RF_Track.Bunch6dT_Generator_cor_py_get, _RF_Track.Bunch6dT_Generator_cor_py_set, doc=r"""
+    cor_py : double
+
+    mrad correlated beam divergence in the vertical direction. Default: 0.
+
+
+    """)
+    charge = property(_RF_Track.Bunch6dT_Generator_charge_get, _RF_Track.Bunch6dT_Generator_charge_set, doc=r"""charge : double""")
+    mass = property(_RF_Track.Bunch6dT_Generator_mass_get, _RF_Track.Bunch6dT_Generator_mass_set, doc=r"""mass : double""")
+    p_ref = property(_RF_Track.Bunch6dT_Generator_p_ref_get, _RF_Track.Bunch6dT_Generator_p_ref_set, doc=r"""p_ref : double""")
+    nmacro = property(_RF_Track.Bunch6dT_Generator_nmacro_get, _RF_Track.Bunch6dT_Generator_nmacro_set, doc=r"""nmacro : double""")
+
+    def __init__(self):
+        r"""
+        __init__(Bunch6dT_Generator self) -> Bunch6dT_Generator
+
+        Distribution and photocathode-emission parameters for bunch generation.
+        Pass this object and the requested macroparticle count to Bunch6dT.
+        The generator uses mixed units: q_total [nC], ref_zpos [m], ref_clock [ns],
+        ref_ekin [MeV], transverse sizes [mm], momentum spreads [keV/c] and emission
+        times [ns]. These differ from the bunch storage units. cathode enables
+        emission-time generation. See individual property help for units.
+
+
+        """
+        _RF_Track.Bunch6dT_Generator_swiginit(self, _RF_Track.new_Bunch6dT_Generator())
+    __swig_destroy__ = _RF_Track.delete_Bunch6dT_Generator
+
+# Register Bunch6dT_Generator in _RF_Track:
+_RF_Track.Bunch6dT_Generator_swigregister(Bunch6dT_Generator)
+randu = _RF_Track.randu
+qrand = _RF_Track.qrand
+qrandn = _RF_Track.qrandn
+class BeamT(object):
+    r"""
+
+    A multi-bunch beam containing Bunch6dT objects. Use this container to track bunch trains
+    and interactions between bunches. Bunch populations and macroparticle counts remain
+    distinct.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(BeamT self, Beam beam) -> BeamT
+        __init__(BeamT self, Bunch6dT bunch) -> BeamT
+        __init__(BeamT self, Beam b, Frame f, bool backtrackng=False) -> BeamT
+        __init__(BeamT self, size_t n, Bunch6dT bunch, double bunch_spacing) -> BeamT
+        __init__(BeamT self) -> BeamT
+
+        A multi-bunch beam containing Bunch6dT objects. Use this container to track bunch trains
+        and interactions between bunches. Bunch populations and macroparticle counts remain
+        distinct.
+
+
+        """
+        _RF_Track.BeamT_swiginit(self, _RF_Track.new_BeamT(*args))
+    __swig_destroy__ = _RF_Track.delete_BeamT
+    append = _swig_new_instance_method(_RF_Track.BeamT_append)
+    size = _swig_new_instance_method(_RF_Track.BeamT_size)
+    load = _swig_new_instance_method(_RF_Track.BeamT_load)
+    save = _swig_new_instance_method(_RF_Track.BeamT_save)
+    __getitem__ = _swig_new_instance_method(_RF_Track.BeamT___getitem__)
+    __len__ = _swig_new_instance_method(_RF_Track.BeamT___len__)
+    __iter__ = _swig_new_instance_method(_RF_Track.BeamT___iter__)
+
+# Register BeamT in _RF_Track:
+_RF_Track.BeamT_swigregister(BeamT)
+class Beam(object):
+    r"""
+
+    A multi-bunch beam containing Bunch6d objects. Use this container to track bunch trains
+    and interactions between bunches. Bunch populations and macroparticle counts remain
+    distinct.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Beam self, Bunch6d bunch) -> Beam
+        __init__(Beam self, size_t n, Bunch6d bunch, double bunch_spacing) -> Beam
+        __init__(Beam self) -> Beam
+
+        A multi-bunch beam containing Bunch6d objects. Use this container to track bunch trains
+        and interactions between bunches. Bunch populations and macroparticle counts remain
+        distinct.
+
+
+        """
+        _RF_Track.Beam_swiginit(self, _RF_Track.new_Beam(*args))
+    __swig_destroy__ = _RF_Track.delete_Beam
+    __assign__ = _swig_new_instance_method(_RF_Track.Beam___assign__)
+    append = _swig_new_instance_method(_RF_Track.Beam_append)
+    size = _swig_new_instance_method(_RF_Track.Beam_size)
+    load = _swig_new_instance_method(_RF_Track.Beam_load)
+    save = _swig_new_instance_method(_RF_Track.Beam_save)
+    __getitem__ = _swig_new_instance_method(_RF_Track.Beam___getitem__)
+    __len__ = _swig_new_instance_method(_RF_Track.Beam___len__)
+    __iter__ = _swig_new_instance_method(_RF_Track.Beam___iter__)
+
+# Register Beam in _RF_Track:
+_RF_Track.Beam_swigregister(Beam)
+class CollectiveEffect(object):
+    r"""
+
+    Base interface for forces acting collectively on a bunch or beam. Attach an effect with
+    add_collective_effect and configure the kick interval on the tracking environment.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined - class is abstract")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_CollectiveEffect
+    compute_force = _swig_new_instance_method(_RF_Track.CollectiveEffect_compute_force)
+
+# Register CollectiveEffect in _RF_Track:
+_RF_Track.CollectiveEffect_swigregister(CollectiveEffect)
+class SpaceCharge(object):
+    r"""
+
+    Base interface for space-charge solvers. Assign a concrete solver using set_sc_engine.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined - class is abstract")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge
+    compute_force = _swig_new_instance_method(_RF_Track.SpaceCharge_compute_force)
+
+# Register SpaceCharge in _RF_Track:
+_RF_Track.SpaceCharge_swigregister(SpaceCharge)
+class SpaceCharge_P2P(SpaceCharge):
+    r"""
+
+    Direct particle-to-particle space-charge solver. Computes pair interactions without a particle-in-cell mesh.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    set_Dmin = _swig_new_instance_method(_RF_Track.SpaceCharge_P2P_set_Dmin)
+    get_Dmin = _swig_new_instance_method(_RF_Track.SpaceCharge_P2P_get_Dmin)
+
+    def __init__(self):
+        r"""
+        __init__(SpaceCharge_P2P self) -> SpaceCharge_P2P
+
+        Direct particle-to-particle space-charge solver. Computes pair interactions without a particle-in-cell mesh.
+
+
+        """
+        _RF_Track.SpaceCharge_P2P_swiginit(self, _RF_Track.new_SpaceCharge_P2P())
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_P2P
+
+# Register SpaceCharge_P2P in _RF_Track:
+_RF_Track.SpaceCharge_P2P_swigregister(SpaceCharge_P2P)
+class SpaceCharge_PIC_FreeSpace(SpaceCharge):
+    r"""
+
+    Particle-in-cell space-charge solver. The exported specializations select free-space,
+    conducting-cylinder or horizontal-plate boundary conditions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(SpaceCharge_PIC_FreeSpace self, size_t Nx_=16, size_t Ny_=16, size_t Nz_=16) -> SpaceCharge_PIC_FreeSpace
+        __init__(SpaceCharge_PIC_FreeSpace self, SpaceCharge_PIC_FreeSpace sc) -> SpaceCharge_PIC_FreeSpace
+
+        Particle-in-cell space-charge solver. The exported specializations select free-space,
+        conducting-cylinder or horizontal-plate boundary conditions.
+
+
+        """
+        _RF_Track.SpaceCharge_PIC_FreeSpace_swiginit(self, _RF_Track.new_SpaceCharge_PIC_FreeSpace(*args))
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_PIC_FreeSpace
+    __assign__ = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace___assign__)
+    set_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_set_mirror)
+    set_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_set_smooth)
+    set_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_set_velocity_slices)
+    set_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_set_mirror_charge_tolerance)
+    get_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_get_mirror)
+    get_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_get_smooth)
+    get_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_get_velocity_slices)
+    get_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_FreeSpace_get_mirror_charge_tolerance)
+
+# Register SpaceCharge_PIC_FreeSpace in _RF_Track:
+_RF_Track.SpaceCharge_PIC_FreeSpace_swigregister(SpaceCharge_PIC_FreeSpace)
+class SpaceCharge_PIC_LongCylinder_T(SpaceCharge):
+    r"""
+
+    Particle-in-cell space-charge solver. The exported specializations select free-space,
+    conducting-cylinder or horizontal-plate boundary conditions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_PIC_LongCylinder_T
+    __assign__ = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T___assign__)
+    set_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_set_mirror)
+    set_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_set_smooth)
+    set_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_set_velocity_slices)
+    set_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_set_mirror_charge_tolerance)
+    get_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_get_mirror)
+    get_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_get_smooth)
+    get_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_get_velocity_slices)
+    get_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_T_get_mirror_charge_tolerance)
+
+# Register SpaceCharge_PIC_LongCylinder_T in _RF_Track:
+_RF_Track.SpaceCharge_PIC_LongCylinder_T_swigregister(SpaceCharge_PIC_LongCylinder_T)
+class SpaceCharge_PIC_HorizontalPlates_T(SpaceCharge):
+    r"""
+
+    Particle-in-cell space-charge solver. The exported specializations select free-space,
+    conducting-cylinder or horizontal-plate boundary conditions.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+
+    def __init__(self, *args, **kwargs):
+        raise AttributeError("No constructor defined")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_PIC_HorizontalPlates_T
+    __assign__ = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T___assign__)
+    set_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_set_mirror)
+    set_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_set_smooth)
+    set_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_set_velocity_slices)
+    set_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_set_mirror_charge_tolerance)
+    get_mirror = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_get_mirror)
+    get_smooth = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_get_smooth)
+    get_velocity_slices = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_get_velocity_slices)
+    get_mirror_charge_tolerance = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_get_mirror_charge_tolerance)
+
+# Register SpaceCharge_PIC_HorizontalPlates_T in _RF_Track:
+_RF_Track.SpaceCharge_PIC_HorizontalPlates_T_swigregister(SpaceCharge_PIC_HorizontalPlates_T)
+class MultipoleKick(CollectiveEffect):
+    r"""
+
+    Collective-effect representation of magnetic multipolar errors, applied as kicks.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(MultipoleKick self, std::vector< std::complex< double > > const & Bn_=std::vector< std::complex< double > >(size_t(0))) -> MultipoleKick
+
+        Collective-effect representation of magnetic multipolar errors, applied as kicks.
+
+
+        """
+        _RF_Track.MultipoleKick_swiginit(self, _RF_Track.new_MultipoleKick(*args))
+    set_Bn = _swig_new_instance_method(_RF_Track.MultipoleKick_set_Bn)
+    get_Bn = _swig_new_instance_method(_RF_Track.MultipoleKick_get_Bn)
+    compute_force = _swig_new_instance_method(_RF_Track.MultipoleKick_compute_force)
+    __swig_destroy__ = _RF_Track.delete_MultipoleKick
+
+# Register MultipoleKick in _RF_Track:
+_RF_Track.MultipoleKick_swigregister(MultipoleKick)
+class BeamLoading(CollectiveEffect):
+    r"""
+
+    Travelling-wave beam-loading model for relativistic beams. Constructor inputs are cell
+    count Ncells, frequency f0 [Hz], phaseadvance [rad], and longitudinal profiles QQ
+    (dimensionless), R_Q [ohm/m] and VG [c]. nbins controls the short-range wake binning.
+    The steady-state overload additionally takes single-particle charge q_single [e],
+    physical particles_bunch and bunch frequency fb [Hz]. Supports both spatial and
+    time-based bunch models.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(BeamLoading self, size_t Ncells, double f0, double phaseadvance, Mesh1d const & QQ, Mesh1d const & R_Q, Mesh1d const & VG, size_t nbins=16) -> BeamLoading
+        __init__(BeamLoading self, size_t Ncells, double f0, double phaseadvance, Mesh1d const & QQ, Mesh1d const & R_Q, Mesh1d const & VG, double q_single, double particles_bunch, double fb, size_t nbins=16) -> BeamLoading
+
+        Travelling-wave beam-loading model for relativistic beams. Constructor inputs are cell
+        count Ncells, frequency f0 [Hz], phaseadvance [rad], and longitudinal profiles QQ
+        (dimensionless), R_Q [ohm/m] and VG [c]. nbins controls the short-range wake binning.
+        The steady-state overload additionally takes single-particle charge q_single [e],
+        physical particles_bunch and bunch frequency fb [Hz]. Supports both spatial and
+        time-based bunch models.
+
+
+        """
+        _RF_Track.BeamLoading_swiginit(self, _RF_Track.new_BeamLoading(*args))
+    compute_force = _swig_new_instance_method(_RF_Track.BeamLoading_compute_force)
+    solve_pde_steady = _swig_new_instance_method(_RF_Track.BeamLoading_solve_pde_steady)
+    solve_pde_transient = _swig_new_instance_method(_RF_Track.BeamLoading_solve_pde_transient)
+    get_vg = _swig_new_instance_method(_RF_Track.BeamLoading_get_vg)
+    get_Q = _swig_new_instance_method(_RF_Track.BeamLoading_get_Q)
+    get_r_Q = _swig_new_instance_method(_RF_Track.BeamLoading_get_r_Q)
+    get_dr_Q = _swig_new_instance_method(_RF_Track.BeamLoading_get_dr_Q)
+    get_G = _swig_new_instance_method(_RF_Track.BeamLoading_get_G)
+    get_G_unloaded = _swig_new_instance_method(_RF_Track.BeamLoading_get_G_unloaded)
+    get_G_steady = _swig_new_instance_method(_RF_Track.BeamLoading_get_G_steady)
+    get_wake_function = _swig_new_instance_method(_RF_Track.BeamLoading_get_wake_function)
+    get_Lcell = _swig_new_instance_method(_RF_Track.BeamLoading_get_Lcell)
+    get_dt = _swig_new_instance_method(_RF_Track.BeamLoading_get_dt)
+    get_z0 = _swig_new_instance_method(_RF_Track.BeamLoading_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.BeamLoading_get_z1)
+    get_tfill = _swig_new_instance_method(_RF_Track.BeamLoading_get_tfill)
+    set_z0 = _swig_new_instance_method(_RF_Track.BeamLoading_set_z0)
+    set_unloaded_gradient = _swig_new_instance_method(_RF_Track.BeamLoading_set_unloaded_gradient)
+    set_early_injection = _swig_new_instance_method(_RF_Track.BeamLoading_set_early_injection)
+    disable_short_range = _swig_new_instance_method(_RF_Track.BeamLoading_disable_short_range)
+    enable_short_range = _swig_new_instance_method(_RF_Track.BeamLoading_enable_short_range)
+    __swig_destroy__ = _RF_Track.delete_BeamLoading
+
+# Register BeamLoading in _RF_Track:
+_RF_Track.BeamLoading_swigregister(BeamLoading)
+class BeamLoadingSW(CollectiveEffect):
+    r"""
+
+    Standing-wave beam-loading model for Bunch6dT or BeamT in a Volume. Initialize with an
+    RF_Field, quality factor QQ, RHO=R/Q per unit length [ohm/m], cell count Ncells, mass
+    [MeV/c^2], single-particle charge q [e], and tinj_frac. The injection time is tinj_frac
+    times the filling time. Frequency, QQ, RHO, Ncells, mass and tinj_frac must be positive;
+    q must be nonzero. The Bunch6d and Beam force overloads are not implemented.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, Ez, QQ, RHO, Ncells, mass, q, tinj_frac):
+        r"""
+        __init__(BeamLoadingSW self, RF_Field Ez, double QQ, double RHO, double Ncells, double mass, double q, double tinj_frac) -> BeamLoadingSW
+
+        Standing-wave beam-loading model for Bunch6dT or BeamT in a Volume. Initialize with an
+        RF_Field, quality factor QQ, RHO=R/Q per unit length [ohm/m], cell count Ncells, mass
+        [MeV/c^2], single-particle charge q [e], and tinj_frac. The injection time is tinj_frac
+        times the filling time. Frequency, QQ, RHO, Ncells, mass and tinj_frac must be positive;
+        q must be nonzero. The Bunch6d and Beam force overloads are not implemented.
+
+
+        """
+        _RF_Track.BeamLoadingSW_swiginit(self, _RF_Track.new_BeamLoadingSW(Ez, QQ, RHO, Ncells, mass, q, tinj_frac))
+    compute_force = _swig_new_instance_method(_RF_Track.BeamLoadingSW_compute_force)
+    get_TT1 = _swig_new_instance_method(_RF_Track.BeamLoadingSW_get_TT1)
+    get_TT2 = _swig_new_instance_method(_RF_Track.BeamLoadingSW_get_TT2)
+    get_Lcell = _swig_new_instance_method(_RF_Track.BeamLoadingSW_get_Lcell)
+    get_tfill = _swig_new_instance_method(_RF_Track.BeamLoadingSW_get_tfill)
+    get_tinj = _swig_new_instance_method(_RF_Track.BeamLoadingSW_get_tinj)
+    __swig_destroy__ = _RF_Track.delete_BeamLoadingSW
+
+# Register BeamLoadingSW in _RF_Track:
+_RF_Track.BeamLoadingSW_swigregister(BeamLoadingSW)
+class LongRangeWakefield(CollectiveEffect):
+    r"""
+
+    Long-range wakefield effect for interactions between separated bunches.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(LongRangeWakefield self, MatrixNd const & A_, MatrixNd const & f1_, MatrixNd const & Q__, MatrixNd const & angle_, double z0_center, double z1_center) -> LongRangeWakefield
+        __init__(LongRangeWakefield self, MatrixNd const & A_, MatrixNd const & f1_, MatrixNd const & Q__, double z0_center, double z1_center) -> LongRangeWakefield
+        __init__(LongRangeWakefield self, MatrixNd const & A_, MatrixNd const & f1_, MatrixNd const & Q__, MatrixNd const & angle_) -> LongRangeWakefield
+        __init__(LongRangeWakefield self, MatrixNd const & A_, MatrixNd const & f1_, MatrixNd const & Q__) -> LongRangeWakefield
+
+        Long-range wakefield effect for interactions between separated bunches.
+
+
+        """
+        _RF_Track.LongRangeWakefield_swiginit(self, _RF_Track.new_LongRangeWakefield(*args))
+    set_nbins = _swig_new_instance_method(_RF_Track.LongRangeWakefield_set_nbins)
+    compute_force = _swig_new_instance_method(_RF_Track.LongRangeWakefield_compute_force)
+    w_transv = _swig_new_instance_method(_RF_Track.LongRangeWakefield_w_transv)
+    __swig_destroy__ = _RF_Track.delete_LongRangeWakefield
+
+# Register LongRangeWakefield in _RF_Track:
+_RF_Track.LongRangeWakefield_swigregister(LongRangeWakefield)
+class ResistiveWallWake(CollectiveEffect):
+    r"""
+
+    Wakefield effect from the finite conductivity of the beam-pipe wall.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(ResistiveWallWake self, std::string material) -> ResistiveWallWake
+        __init__(ResistiveWallWake self, double sigma_c) -> ResistiveWallWake
+        __init__(ResistiveWallWake self, double b, double sigma_c) -> ResistiveWallWake
+        __init__(ResistiveWallWake self, double b, std::string material) -> ResistiveWallWake
+        __init__(ResistiveWallWake self) -> ResistiveWallWake
+
+        Wakefield effect from the finite conductivity of the beam-pipe wall.
+
+
+        """
+        _RF_Track.ResistiveWallWake_swiginit(self, _RF_Track.new_ResistiveWallWake(*args))
+    set_nbins = _swig_new_instance_method(_RF_Track.ResistiveWallWake_set_nbins)
+    w_transv = _swig_new_instance_method(_RF_Track.ResistiveWallWake_w_transv)
+    w_long = _swig_new_instance_method(_RF_Track.ResistiveWallWake_w_long)
+    __swig_destroy__ = _RF_Track.delete_ResistiveWallWake
+
+# Register ResistiveWallWake in _RF_Track:
+_RF_Track.ResistiveWallWake_swigregister(ResistiveWallWake)
+class ShortRangeWakefield(CollectiveEffect):
+    r"""
+
+    Short-range wakefield effect acting within a bunch.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(ShortRangeWakefield self, double a=1, double g=1, double l=1) -> ShortRangeWakefield
+        __init__(ShortRangeWakefield self, VectorNd const & a_, VectorNd const & g_, VectorNd const & l_, size_t ncells_) -> ShortRangeWakefield
+
+        Short-range wakefield effect acting within a bunch.
+
+
+        """
+        _RF_Track.ShortRangeWakefield_swiginit(self, _RF_Track.new_ShortRangeWakefield(*args))
+    set_nbins = _swig_new_instance_method(_RF_Track.ShortRangeWakefield_set_nbins)
+    scatter_cells = _swig_new_instance_method(_RF_Track.ShortRangeWakefield_scatter_cells)
+    compute_force = _swig_new_instance_method(_RF_Track.ShortRangeWakefield_compute_force)
+    w_long = _swig_new_instance_method(_RF_Track.ShortRangeWakefield_w_long)
+    w_transv = _swig_new_instance_method(_RF_Track.ShortRangeWakefield_w_transv)
+    __swig_destroy__ = _RF_Track.delete_ShortRangeWakefield
+
+# Register ShortRangeWakefield in _RF_Track:
+_RF_Track.ShortRangeWakefield_swigregister(ShortRangeWakefield)
+class MultipleCoulombScattering(CollectiveEffect):
+    r"""
+
+    Multiple Coulomb scattering in matter. Material parameters include radiation length in cm.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(MultipleCoulombScattering self, double radiation_length=0.0, double Z=0.0, double A=0.0, double density=0.0, double mean_excitation_energy=-1.0) -> MultipleCoulombScattering
+        __init__(MultipleCoulombScattering self, std::string const & material_name) -> MultipleCoulombScattering
+        __init__(MultipleCoulombScattering self, Material const & m) -> MultipleCoulombScattering
+
+        Multiple Coulomb scattering in matter. Material parameters include radiation length in cm.
+
+
+        """
+        _RF_Track.MultipleCoulombScattering_swiginit(self, _RF_Track.new_MultipleCoulombScattering(*args))
+    enable_log_term = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_enable_log_term)
+    disable_log_term = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_disable_log_term)
+    enable_fruehwirth_model = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_enable_fruehwirth_model)
+    disable_fruehwirth_model = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_disable_fruehwirth_model)
+    enable_wentzel_model = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_enable_wentzel_model)
+    disable_wentzel_model = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_disable_wentzel_model)
+    compute_force = _swig_new_instance_method(_RF_Track.MultipleCoulombScattering_compute_force)
+    __swig_destroy__ = _RF_Track.delete_MultipleCoulombScattering
+
+# Register MultipleCoulombScattering in _RF_Track:
+_RF_Track.MultipleCoulombScattering_swigregister(MultipleCoulombScattering)
+class IntraBeamScattering(CollectiveEffect):
+    r"""
+
+    Collective effect modelling intra-beam scattering. Attach with add_collective_effect and
+    configure collective-effect kick sampling.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, nx=16, ny=16, nz=16, ncol=-1):
+        r"""
+        __init__(IntraBeamScattering self, size_t nx=16, size_t ny=16, size_t nz=16, int ncol=-1) -> IntraBeamScattering
+
+        Collective effect modelling intra-beam scattering. Attach with add_collective_effect and
+        configure collective-effect kick sampling.
+
+
+        """
+        _RF_Track.IntraBeamScattering_swiginit(self, _RF_Track.new_IntraBeamScattering(nx, ny, nz, ncol))
+    __swig_destroy__ = _RF_Track.delete_IntraBeamScattering
+
+# Register IntraBeamScattering in _RF_Track:
+_RF_Track.IntraBeamScattering_swigregister(IntraBeamScattering)
+class EnergyStraggling(CollectiveEffect):
+    r"""
+
+    Fluctuations of energy loss through matter.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(EnergyStraggling self, double radiation_length=0.0, double Z=0.0, double A=0.0, double density=0.0, double mean_excitation_energy=-1.0) -> EnergyStraggling
+        __init__(EnergyStraggling self, std::string const & material_name) -> EnergyStraggling
+        __init__(EnergyStraggling self, Material const & m) -> EnergyStraggling
+
+        Fluctuations of energy loss through matter.
+
+
+        """
+        _RF_Track.EnergyStraggling_swiginit(self, _RF_Track.new_EnergyStraggling(*args))
+    compute_force = _swig_new_instance_method(_RF_Track.EnergyStraggling_compute_force)
+    __swig_destroy__ = _RF_Track.delete_EnergyStraggling
+
+# Register EnergyStraggling in _RF_Track:
+_RF_Track.EnergyStraggling_swigregister(EnergyStraggling)
+class StoppingPower(CollectiveEffect):
+    r"""
+
+    Mean energy loss through matter due to stopping power.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(StoppingPower self, double radiation_length=0.0, double Z=0.0, double A=0.0, double density=0.0, double mean_excitation_energy=-1.0) -> StoppingPower
+        __init__(StoppingPower self, std::string const & material_name) -> StoppingPower
+        __init__(StoppingPower self, StoppingPower m) -> StoppingPower
+        __init__(StoppingPower self, Material const & m) -> StoppingPower
+
+        Mean energy loss through matter due to stopping power.
+
+
+        """
+        _RF_Track.StoppingPower_swiginit(self, _RF_Track.new_StoppingPower(*args))
+    __swig_destroy__ = _RF_Track.delete_StoppingPower
+    __assign__ = _swig_new_instance_method(_RF_Track.StoppingPower___assign__)
+
+# Register StoppingPower in _RF_Track:
+_RF_Track.StoppingPower_swigregister(StoppingPower)
+class IncoherentSynchrotronRadiation(CollectiveEffect):
+    r"""
+
+    Radiative energy loss and quantum excitation from incoherent synchrotron radiation.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, quantum=False, create_photons=False):
+        r"""
+        __init__(IncoherentSynchrotronRadiation self, bool quantum=False, bool create_photons=False) -> IncoherentSynchrotronRadiation
+
+        Radiative energy loss and quantum excitation from incoherent synchrotron radiation.
+
+
+        """
+        _RF_Track.IncoherentSynchrotronRadiation_swiginit(self, _RF_Track.new_IncoherentSynchrotronRadiation(quantum, create_photons))
+    compute_force = _swig_new_instance_method(_RF_Track.IncoherentSynchrotronRadiation_compute_force)
+    __swig_destroy__ = _RF_Track.delete_IncoherentSynchrotronRadiation
+
+# Register IncoherentSynchrotronRadiation in _RF_Track:
+_RF_Track.IncoherentSynchrotronRadiation_swigregister(IncoherentSynchrotronRadiation)
+class SpaceCharge_PIC_LongCylinder(SpaceCharge_PIC_LongCylinder_T):
+    r"""
+
+    Particle-in-cell space-charge solver with a long conducting cylindrical boundary.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, Nx_=1, Ny_=1, Nz_=1, a=1.0):
+        r"""
+        __init__(SpaceCharge_PIC_LongCylinder self, size_t Nx_=1, size_t Ny_=1, size_t Nz_=1, double a=1.0) -> SpaceCharge_PIC_LongCylinder
+
+        Particle-in-cell space-charge solver with a long conducting cylindrical boundary.
+
+
+        """
+        _RF_Track.SpaceCharge_PIC_LongCylinder_swiginit(self, _RF_Track.new_SpaceCharge_PIC_LongCylinder(Nx_, Ny_, Nz_, a))
+    set_aperture = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_set_aperture)
+    get_aperture = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_LongCylinder_get_aperture)
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_PIC_LongCylinder
+
+# Register SpaceCharge_PIC_LongCylinder in _RF_Track:
+_RF_Track.SpaceCharge_PIC_LongCylinder_swigregister(SpaceCharge_PIC_LongCylinder)
+class SpaceCharge_PIC_HorizontalPlates(SpaceCharge_PIC_HorizontalPlates_T):
+    r"""
+
+    Particle-in-cell space-charge solver between conducting horizontal plates.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, Nx_=1, Ny_=1, Nz_=1, h=1.0):
+        r"""
+        __init__(SpaceCharge_PIC_HorizontalPlates self, size_t Nx_=1, size_t Ny_=1, size_t Nz_=1, double h=1.0) -> SpaceCharge_PIC_HorizontalPlates
+
+        Particle-in-cell space-charge solver between conducting horizontal plates.
+
+
+        """
+        _RF_Track.SpaceCharge_PIC_HorizontalPlates_swiginit(self, _RF_Track.new_SpaceCharge_PIC_HorizontalPlates(Nx_, Ny_, Nz_, h))
+    set_half_gap = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_set_half_gap)
+    get_half_gap = _swig_new_instance_method(_RF_Track.SpaceCharge_PIC_HorizontalPlates_get_half_gap)
+    __swig_destroy__ = _RF_Track.delete_SpaceCharge_PIC_HorizontalPlates
+
+# Register SpaceCharge_PIC_HorizontalPlates in _RF_Track:
+_RF_Track.SpaceCharge_PIC_HorizontalPlates_swigregister(SpaceCharge_PIC_HorizontalPlates)
+class Visitor(object):
+    r"""
+
+    Base interface for visiting beamline elements.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+    __swig_destroy__ = _RF_Track.delete_Visitor
+    visit = _swig_new_instance_method(_RF_Track.Visitor_visit)
+    should_stop = _swig_new_instance_method(_RF_Track.Visitor_should_stop)
+
+    def __init__(self):
+        r"""
+        __init__(Visitor self) -> Visitor
+
+        Base interface for visiting beamline elements.
+
+
+        """
+        _RF_Track.Visitor_swiginit(self, _RF_Track.new_Visitor())
+
+# Register Visitor in _RF_Track:
+_RF_Track.Visitor_swigregister(Visitor)
+class UserEffect(CollectiveEffect):
+    r"""
+
+    Subclass in Python or Octave to implement a collective force. Force outputs have one row
+    per macroparticle and three columns in MeV/m. Time and spatial tracking use different
+    compute_force overloads; consult the extension examples in manual chapter 9.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self):
+        r"""
+        __init__(UserEffect self) -> UserEffect
+
+        Subclass in Python or Octave to implement a collective force. Force outputs have one row
+        per macroparticle and three columns in MeV/m. Time and spatial tracking use different
+        compute_force overloads; consult the extension examples in manual chapter 9.
+
+
+        """
+        if self.__class__ == UserEffect:
+            _self = None
+        else:
+            _self = self
+        _RF_Track.UserEffect_swiginit(self, _RF_Track.new_UserEffect(_self, ))
+    __swig_destroy__ = _RF_Track.delete_UserEffect
+    compute_force = _swig_new_instance_method(_RF_Track.UserEffect_compute_force)
+    def __disown__(self):
+        self.this.disown()
+        _RF_Track.disown_UserEffect(self)
+        return weakref.proxy(self)
+
+# Register UserEffect in _RF_Track:
+_RF_Track.UserEffect_swigregister(UserEffect)
+class UserElement(TimeDependent_Field):
+    r"""
+
+    Subclass in Python or Octave to implement a custom lattice element. Initialize the
+    parent with the element length in metres. Override track(bunch) and modify its Bunch6d
+    argument in place; no return value is used. Override btrack for backward tracking. See
+    manual chapter 9 for director/subclass examples.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, L=0.0):
+        r"""
+        __init__(UserElement self, double L=0.0) -> UserElement
+
+        Subclass in Python or Octave to implement a custom lattice element. Initialize the
+        parent with the element length in metres. Override track(bunch) and modify its Bunch6d
+        argument in place; no return value is used. Override btrack for backward tracking. See
+        manual chapter 9 for director/subclass examples.
+
+
+        """
+        if self.__class__ == UserElement:
+            _self = None
+        else:
+            _self = self
+        _RF_Track.UserElement_swiginit(self, _RF_Track.new_UserElement(_self, L))
+    __swig_destroy__ = _RF_Track.delete_UserElement
+    get_length = _swig_new_instance_method(_RF_Track.UserElement_get_length)
+    track = _swig_new_instance_method(_RF_Track.UserElement_track)
+    btrack = _swig_new_instance_method(_RF_Track.UserElement_btrack)
+    def __disown__(self):
+        self.this.disown()
+        _RF_Track.disown_UserElement(self)
+        return weakref.proxy(self)
+
+# Register UserElement in _RF_Track:
+_RF_Track.UserElement_swigregister(UserElement)
+class UserField(TimeDependent_Field):
+    r"""
+
+    Subclass in Python or Octave to implement an electromagnetic field. Override
+    get_field(x,y,z,t), with positions in mm and time in mm/c, returning E [V/m] and B [T].
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, L=0):
+        r"""
+        __init__(UserField self, double L=0) -> UserField
+
+        Subclass in Python or Octave to implement an electromagnetic field. Override
+        get_field(x,y,z,t), with positions in mm and time in mm/c, returning E [V/m] and B [T].
+
+
+        """
+        if self.__class__ == UserField:
+            _self = None
+        else:
+            _self = self
+        _RF_Track.UserField_swiginit(self, _RF_Track.new_UserField(_self, L))
+    __swig_destroy__ = _RF_Track.delete_UserField
+    get_length = _swig_new_instance_method(_RF_Track.UserField_get_length)
+    get_field = _swig_new_instance_method(_RF_Track.UserField_get_field)
+    def __disown__(self):
+        self.this.disown()
+        _RF_Track.disown_UserField(self)
+        return weakref.proxy(self)
+
+# Register UserField in _RF_Track:
+_RF_Track.UserField_swigregister(UserField)
+class UserVisitor(Visitor):
+    r"""
+
+    Subclass to traverse and inspect or modify beamline elements through the visitor
+    callbacks. Pass the visitor to the beamline accept method.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self):
+        r"""
+        __init__(UserVisitor self) -> UserVisitor
+
+        Subclass to traverse and inspect or modify beamline elements through the visitor
+        callbacks. Pass the visitor to the beamline accept method.
+
+
+        """
+        if self.__class__ == UserVisitor:
+            _self = None
+        else:
+            _self = self
+        _RF_Track.UserVisitor_swiginit(self, _RF_Track.new_UserVisitor(_self, ))
+    __swig_destroy__ = _RF_Track.delete_UserVisitor
+    stop_at = _swig_new_instance_method(_RF_Track.UserVisitor_stop_at)
+    visit = _swig_new_instance_method(_RF_Track.UserVisitor_visit)
+    def __disown__(self):
+        self.this.disown()
+        _RF_Track.disown_UserVisitor(self)
+        return weakref.proxy(self)
+
+# Register UserVisitor in _RF_Track:
+_RF_Track.UserVisitor_swigregister(UserVisitor)
+class Wakefield_1d_CINT(CollectiveEffect):
+    r"""
+
+    Wakefield effect from one-dimensional sampled wake functions. Linear and cubic
+    interpolation variants are exposed by the interface.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Wakefield_1d_CINT self, Mesh1d const & W_transv_, Mesh1d const & W_long_, double hz_, size_t nbins_=63) -> Wakefield_1d_CINT
+        __init__(Wakefield_1d_CINT self, Mesh1d const & W_transv_, double hz_, size_t nbins_=63) -> Wakefield_1d_CINT
+        __init__(Wakefield_1d_CINT self, size_t nbins_=63) -> Wakefield_1d_CINT
+
+        Wakefield effect from one-dimensional sampled wake functions. Linear and cubic
+        interpolation variants are exposed by the interface.
+
+
+        """
+        _RF_Track.Wakefield_1d_CINT_swiginit(self, _RF_Track.new_Wakefield_1d_CINT(*args))
+    w_transv = _swig_new_instance_method(_RF_Track.Wakefield_1d_CINT_w_transv)
+    w_long = _swig_new_instance_method(_RF_Track.Wakefield_1d_CINT_w_long)
+    set_nbins = _swig_new_instance_method(_RF_Track.Wakefield_1d_CINT_set_nbins)
+    set_hz = _swig_new_instance_method(_RF_Track.Wakefield_1d_CINT_set_hz)
+    get_hz = _swig_new_instance_method(_RF_Track.Wakefield_1d_CINT_get_hz)
+    __swig_destroy__ = _RF_Track.delete_Wakefield_1d_CINT
+
+# Register Wakefield_1d_CINT in _RF_Track:
+_RF_Track.Wakefield_1d_CINT_swigregister(Wakefield_1d_CINT)
+class Wakefield_1d(CollectiveEffect):
+    r"""
+
+
+
+    Wakefield effect from one-dimensional sampled wake functions. Linear and cubic
+    interpolation variants are exposed by the interface.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Wakefield_1d self, Mesh1d const & W_transv_, Mesh1d const & W_long_, double hz_, size_t nbins_=63) -> Wakefield_1d
+        __init__(Wakefield_1d self, Mesh1d const & W_transv_, double hz_, size_t nbins_=63) -> Wakefield_1d
+        __init__(Wakefield_1d self, size_t nbins_=63) -> Wakefield_1d
+
+        Wakefield effect from one-dimensional sampled wake functions. Linear and cubic
+        interpolation variants are exposed by the interface.
+
+
+        """
+        _RF_Track.Wakefield_1d_swiginit(self, _RF_Track.new_Wakefield_1d(*args))
+    w_transv = _swig_new_instance_method(_RF_Track.Wakefield_1d_w_transv)
+    w_long = _swig_new_instance_method(_RF_Track.Wakefield_1d_w_long)
+    set_nbins = _swig_new_instance_method(_RF_Track.Wakefield_1d_set_nbins)
+    set_hz = _swig_new_instance_method(_RF_Track.Wakefield_1d_set_hz)
+    get_hz = _swig_new_instance_method(_RF_Track.Wakefield_1d_get_hz)
+    __swig_destroy__ = _RF_Track.delete_Wakefield_1d
+
+# Register Wakefield_1d in _RF_Track:
+_RF_Track.Wakefield_1d_swigregister(Wakefield_1d)
+class RF_FieldMap_CINT(RF_Field):
+    r"""
+
+    Three-dimensional RF field map on a regular Cartesian mesh. Electric components are in
+    V/m and magnetic components in T; origins, mesh spacings and element length are in
+    metres. Frequency is in Hz. P_map and P_actual are in W; field amplitude scales with
+    sqrt(P_actual/P_map). Direction selects forward (+1), backward (-1), or static (0)
+    fields. NaN field values mark material boundaries.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap_CINT self, ComplexMesh3d const & Ex, ComplexMesh3d const & Ey, ComplexMesh3d const & Ez, ComplexMesh3d const & Bx, ComplexMesh3d const & By, ComplexMesh3d const & Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_CINT
+        __init__(RF_FieldMap_CINT self, double Ex, double Ey, double Ez, ComplexMesh3d const & Bx, ComplexMesh3d const & By, ComplexMesh3d const & Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_CINT
+        __init__(RF_FieldMap_CINT self, ComplexMesh3d const & Ex, ComplexMesh3d const & Ey, ComplexMesh3d const & Ez, double Bx, double By, double Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_CINT
+        __init__(RF_FieldMap_CINT self) -> RF_FieldMap_CINT
+
+        Three-dimensional RF field map on a regular Cartesian mesh. Electric components are in
+        V/m and magnetic components in T; origins, mesh spacings and element length are in
+        metres. Frequency is in Hz. P_map and P_actual are in W; field amplitude scales with
+        sqrt(P_actual/P_map). Direction selects forward (+1), backward (-1), or static (0)
+        fields. NaN field values mark material boundaries.
+
+
+        """
+        _RF_Track.RF_FieldMap_CINT_swiginit(self, _RF_Track.new_RF_FieldMap_CINT(*args))
+    get_hx = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_hz)
+    get_nx = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_nz)
+    get_xa = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_xa)
+    get_ya = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_ya)
+    get_x0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_x0)
+    get_y0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_y0)
+    get_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_z0)
+    get_x1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_x1)
+    get_y1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_y1)
+    get_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_z1)
+    is_cylindrical = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_is_cylindrical)
+    set_cylindrical = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_cylindrical)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_get_direction)
+    set_hx = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_hx)
+    set_hy = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_hy)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_hz)
+    set_xa = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_xa)
+    set_ya = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_ya)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_z1)
+    set_smooth = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_smooth)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_CINT_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap_CINT
+
+# Register RF_FieldMap_CINT in _RF_Track:
+_RF_Track.RF_FieldMap_CINT_swigregister(RF_FieldMap_CINT)
+class RF_FieldMap(RF_Field):
+    r"""
+
+
+
+    Three-dimensional RF field map on a regular Cartesian mesh. Electric components are in
+    V/m and magnetic components in T; origins, mesh spacings and element length are in
+    metres. Frequency is in Hz. P_map and P_actual are in W; field amplitude scales with
+    sqrt(P_actual/P_map). Direction selects forward (+1), backward (-1), or static (0)
+    fields. NaN field values mark material boundaries.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap self, ComplexMesh3d const & Ex, ComplexMesh3d const & Ey, ComplexMesh3d const & Ez, ComplexMesh3d const & Bx, ComplexMesh3d const & By, ComplexMesh3d const & Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap
+        __init__(RF_FieldMap self, double Ex, double Ey, double Ez, ComplexMesh3d const & Bx, ComplexMesh3d const & By, ComplexMesh3d const & Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap
+        __init__(RF_FieldMap self, ComplexMesh3d const & Ex, ComplexMesh3d const & Ey, ComplexMesh3d const & Ez, double Bx, double By, double Bz, double xa, double ya, double hx, double hy, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap
+        __init__(RF_FieldMap self) -> RF_FieldMap
+
+        Three-dimensional RF field map on a regular Cartesian mesh. Electric components are in
+        V/m and magnetic components in T; origins, mesh spacings and element length are in
+        metres. Frequency is in Hz. P_map and P_actual are in W; field amplitude scales with
+        sqrt(P_actual/P_map). Direction selects forward (+1), backward (-1), or static (0)
+        fields. NaN field values mark material boundaries.
+
+
+        """
+        _RF_Track.RF_FieldMap_swiginit(self, _RF_Track.new_RF_FieldMap(*args))
+    get_hx = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_hx)
+    get_hy = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_hy)
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_hz)
+    get_nx = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_nx)
+    get_ny = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_ny)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_nz)
+    get_xa = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_xa)
+    get_ya = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_ya)
+    get_x0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_x0)
+    get_y0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_y0)
+    get_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_z0)
+    get_x1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_x1)
+    get_y1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_y1)
+    get_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_z1)
+    is_cylindrical = _swig_new_instance_method(_RF_Track.RF_FieldMap_is_cylindrical)
+    set_cylindrical = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_cylindrical)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_get_direction)
+    set_hx = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_hx)
+    set_hy = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_hy)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_hz)
+    set_xa = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_xa)
+    set_ya = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_ya)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_z1)
+    set_smooth = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_smooth)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap
+
+# Register RF_FieldMap in _RF_Track:
+_RF_Track.RF_FieldMap_swigregister(RF_FieldMap)
+class RF_FieldMap_2d_CINT(RF_Field):
+    r"""
+
+    Cylindrically symmetric RF field map on a regular radial/longitudinal mesh. Electric
+    fields are in V/m, magnetic fields in T, geometry in m, frequency in Hz and
+    map/operating powers in W.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap_2d_CINT self, ComplexMesh2d const & Er, ComplexMesh2d const & Ez, ComplexMesh2d const & Bt, ComplexMesh2d const & Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d_CINT
+        __init__(RF_FieldMap_2d_CINT self, double Er, double Ez, ComplexMesh2d const & Bt, ComplexMesh2d const & Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d_CINT
+        __init__(RF_FieldMap_2d_CINT self, ComplexMesh2d const & Er, ComplexMesh2d const & Ez, double Bt, double Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d_CINT
+        __init__(RF_FieldMap_2d_CINT self) -> RF_FieldMap_2d_CINT
+
+        Cylindrically symmetric RF field map on a regular radial/longitudinal mesh. Electric
+        fields are in V/m, magnetic fields in T, geometry in m, frequency in Hz and
+        map/operating powers in W.
+
+
+        """
+        _RF_Track.RF_FieldMap_2d_CINT_swiginit(self, _RF_Track.new_RF_FieldMap_2d_CINT(*args))
+    get_hr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_hr)
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_hz)
+    get_nr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_nr)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_nz)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_get_direction)
+    set_hr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_set_hr)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_set_z1)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_CINT_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap_2d_CINT
+
+# Register RF_FieldMap_2d_CINT in _RF_Track:
+_RF_Track.RF_FieldMap_2d_CINT_swigregister(RF_FieldMap_2d_CINT)
+class RF_FieldMap_2d(RF_Field):
+    r"""
+
+
+
+    Cylindrically symmetric RF field map on a regular radial/longitudinal mesh. Electric
+    fields are in V/m, magnetic fields in T, geometry in m, frequency in Hz and
+    map/operating powers in W.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(RF_FieldMap_2d self, ComplexMesh2d const & Er, ComplexMesh2d const & Ez, ComplexMesh2d const & Bt, ComplexMesh2d const & Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d
+        __init__(RF_FieldMap_2d self, double Er, double Ez, ComplexMesh2d const & Bt, ComplexMesh2d const & Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d
+        __init__(RF_FieldMap_2d self, ComplexMesh2d const & Er, ComplexMesh2d const & Ez, double Bt, double Bz, double hr, double hz, double length, double freq, double direction, double P_map=1.0, double P_actual=1.0) -> RF_FieldMap_2d
+        __init__(RF_FieldMap_2d self) -> RF_FieldMap_2d
+
+        Cylindrically symmetric RF field map on a regular radial/longitudinal mesh. Electric
+        fields are in V/m, magnetic fields in T, geometry in m, frequency in Hz and
+        map/operating powers in W.
+
+
+        """
+        _RF_Track.RF_FieldMap_2d_swiginit(self, _RF_Track.new_RF_FieldMap_2d(*args))
+    get_hr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_hr)
+    get_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_hz)
+    get_nr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_nr)
+    get_nz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_nz)
+    get_P_map = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_P_map)
+    get_P_actual = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_P_actual)
+    get_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_get_direction)
+    set_hr = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_set_hr)
+    set_hz = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_set_z1)
+    set_direction = _swig_new_instance_method(_RF_Track.RF_FieldMap_2d_set_direction)
+    __swig_destroy__ = _RF_Track.delete_RF_FieldMap_2d
+
+# Register RF_FieldMap_2d in _RF_Track:
+_RF_Track.RF_FieldMap_2d_swigregister(RF_FieldMap_2d)
+class Static_Magnetic_FieldMap_1d_CINT(GenericField):
+    r"""
+
+    Axisymmetric static magnetic field reconstructed from an on-axis Bz map in T.
+    Longitudinal spacing and length are in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Magnetic_FieldMap_1d_CINT self, Mesh1d const & static_Bz, double hz, double length=-1) -> Static_Magnetic_FieldMap_1d_CINT
+        __init__(Static_Magnetic_FieldMap_1d_CINT self) -> Static_Magnetic_FieldMap_1d_CINT
+
+        Axisymmetric static magnetic field reconstructed from an on-axis Bz map in T.
+        Longitudinal spacing and length are in metres.
+
+
+        """
+        _RF_Track.Static_Magnetic_FieldMap_1d_CINT_swiginit(self, _RF_Track.new_Static_Magnetic_FieldMap_1d_CINT(*args))
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_get_hz)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_get_z1)
+    set_smooth = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_set_smooth)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_CINT_set_z1)
+    __swig_destroy__ = _RF_Track.delete_Static_Magnetic_FieldMap_1d_CINT
+
+# Register Static_Magnetic_FieldMap_1d_CINT in _RF_Track:
+_RF_Track.Static_Magnetic_FieldMap_1d_CINT_swigregister(Static_Magnetic_FieldMap_1d_CINT)
+class Static_Magnetic_FieldMap_1d(GenericField):
+    r"""
+
+
+
+    Axisymmetric static magnetic field reconstructed from an on-axis Bz map in T.
+    Longitudinal spacing and length are in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Magnetic_FieldMap_1d self, Mesh1d const & static_Bz, double hz, double length=-1) -> Static_Magnetic_FieldMap_1d
+        __init__(Static_Magnetic_FieldMap_1d self) -> Static_Magnetic_FieldMap_1d
+
+        Axisymmetric static magnetic field reconstructed from an on-axis Bz map in T.
+        Longitudinal spacing and length are in metres.
+
+
+        """
+        _RF_Track.Static_Magnetic_FieldMap_1d_swiginit(self, _RF_Track.new_Static_Magnetic_FieldMap_1d(*args))
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_get_hz)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_get_z1)
+    set_smooth = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_set_smooth)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_1d_set_z1)
+    __swig_destroy__ = _RF_Track.delete_Static_Magnetic_FieldMap_1d
+
+# Register Static_Magnetic_FieldMap_1d in _RF_Track:
+_RF_Track.Static_Magnetic_FieldMap_1d_swigregister(Static_Magnetic_FieldMap_1d)
+class Static_Magnetic_FieldMap_2d_CINT(GenericField):
+    r"""
+
+    Axisymmetric static magnetic field supplied as Br and Bz on a regular
+    radial/longitudinal mesh. Fields are in T and geometry arguments in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Magnetic_FieldMap_2d_CINT self, Mesh2d const & static_Br, Mesh2d const & static_Bz, double hr, double hz, double length=-1) -> Static_Magnetic_FieldMap_2d_CINT
+        __init__(Static_Magnetic_FieldMap_2d_CINT self) -> Static_Magnetic_FieldMap_2d_CINT
+
+        Axisymmetric static magnetic field supplied as Br and Bz on a regular
+        radial/longitudinal mesh. Fields are in T and geometry arguments in metres.
+
+
+        """
+        _RF_Track.Static_Magnetic_FieldMap_2d_CINT_swiginit(self, _RF_Track.new_Static_Magnetic_FieldMap_2d_CINT(*args))
+    get_hr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_hr)
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_hz)
+    get_nr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_nr)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_get_z1)
+    set_hr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_set_hr)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_CINT_set_z1)
+    __swig_destroy__ = _RF_Track.delete_Static_Magnetic_FieldMap_2d_CINT
+
+# Register Static_Magnetic_FieldMap_2d_CINT in _RF_Track:
+_RF_Track.Static_Magnetic_FieldMap_2d_CINT_swigregister(Static_Magnetic_FieldMap_2d_CINT)
+class Static_Magnetic_FieldMap_2d(GenericField):
+    r"""
+
+
+
+    Axisymmetric static magnetic field supplied as Br and Bz on a regular
+    radial/longitudinal mesh. Fields are in T and geometry arguments in metres.
+
+
+    """
+
+    thisown = property(lambda x: x.this.own(), lambda x, v: x.this.own(v), doc="The membership flag")
+    __repr__ = _swig_repr
+
+    def __init__(self, *args):
+        r"""
+        __init__(Static_Magnetic_FieldMap_2d self, Mesh2d const & static_Br, Mesh2d const & static_Bz, double hr, double hz, double length=-1) -> Static_Magnetic_FieldMap_2d
+        __init__(Static_Magnetic_FieldMap_2d self) -> Static_Magnetic_FieldMap_2d
+
+        Axisymmetric static magnetic field supplied as Br and Bz on a regular
+        radial/longitudinal mesh. Fields are in T and geometry arguments in metres.
+
+
+        """
+        _RF_Track.Static_Magnetic_FieldMap_2d_swiginit(self, _RF_Track.new_Static_Magnetic_FieldMap_2d(*args))
+    get_hr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_hr)
+    get_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_hz)
+    get_nr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_nr)
+    get_nz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_nz)
+    get_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_z0)
+    get_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_get_z1)
+    set_hr = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_set_hr)
+    set_hz = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_set_hz)
+    set_z0 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_set_z0)
+    set_z1 = _swig_new_instance_method(_RF_Track.Static_Magnetic_FieldMap_2d_set_z1)
+    __swig_destroy__ = _RF_Track.delete_Static_Magnetic_FieldMap_2d
+
+# Register Static_Magnetic_FieldMap_2d in _RF_Track:
+_RF_Track.Static_Magnetic_FieldMap_2d_swigregister(Static_Magnetic_FieldMap_2d)
+
+Element.deepcopy = lambda self: self.clone()
+Particle.id  = property(Particle.get_id,  Particle.set_id)
+ParticleT.id = property(ParticleT.get_id, ParticleT.set_id)
+Particle.pdg_id  = property(Particle.get_pdg_id,  Particle.set_pdg_id)
+ParticleT.pdg_id = property(ParticleT.get_pdg_id, ParticleT.set_pdg_id)
+
+# SWIG exposes writable C++ globals only through ``cvar``.  Forward them at
+# module level as well, while retaining ``cvar`` for backwards compatibility.
+# A module subclass is needed so that assignments update the C++ variable
+# instead of merely replacing a Python module attribute.
+import sys as _sys
+import types as _types
+
+_rft_writable_globals = frozenset(("number_of_threads", "SC_engine"))
+
+class _RFTrackModule(_types.ModuleType):
+    def __getattribute__(self, name):
+        if name in _rft_writable_globals:
+            return getattr(cvar, name)
+        return _types.ModuleType.__getattribute__(self, name)
+
+    def __setattr__(self, name, value):
+        if name in _rft_writable_globals:
+            setattr(cvar, name, value)
+        else:
+            _types.ModuleType.__setattr__(self, name, value)
+
+    def __dir__(self):
+        return sorted(set(_types.ModuleType.__dir__(self)) |
+                      _rft_writable_globals)
+
+_sys.modules[__name__].__class__ = _RFTrackModule
+
+