! Example of off-momentum DA calculation

&run_setup
	lattice =aps.lte,
	use_beamline = ringrf,
	p_central_mev = 3000,
        semaphore_file = %s.done,
          tracking_updates = 0,
&end

&twiss_output
        output_at_each_step = 1,
        filename = %s.twi,
&end

&alter_elements name=MALIN, item=DP, value=0.04 &end

&run_control
        n_passes = 1024
&end

&error_control error_log = %s.erl &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KQUAD, item = DX, bind = 0, amplitude = 40e-6 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KQUAD, item = DY, bind = 0, amplitude = 40e-6 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KSEXT, item = DX, bind = 0, amplitude = 40e-6 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KSEXT, item = DY, bind = 0, amplitude = 40e-6 &end
&error_element name = be1, type = gaussian, cutoff = 1.0, element_type=CSBEN, item = DX, bind = 0, amplitude = 40e-6 &end
&error_element name = be1, type = gaussian, cutoff = 1.0, element_type=CSBEN, item = DY, bind = 0, amplitude = 40e-6 &end
&error_element name = be2, type = gaussian, cutoff = 1.0, element_type=CSBEN, item = DX, bind = 0, amplitude = 40e-6 &end
&error_element name = be2, type = gaussian, cutoff = 1.0, element_type=CSBEN, item = DY, bind = 0, amplitude = 40e-6 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KQUAD, item = FSE, bind = 0, amplitude = 4.00e-4 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KSEXT, item = FSE, bind = 0, amplitude = 1.00e-3 &end
&error_element name = be1, type = gaussian, cutoff = 1.0, element_type=CSBEN,item = FSE, bind = 0, amplitude = 1e-4 &end
&error_element name = be2, type = gaussian, cutoff = 1.0, element_type=CSBEN,item = FSE, bind = 0, amplitude = 1e-4 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KQUAD, item = TILT, bind = 0, amplitude = 1.0e-4 &end
&error_element name = *, type = gaussian, cutoff = 1.0, element_type=KSEXT, item = TILT, bind = 0, amplitude = 1.0e-4 &end
&error_element name = be1, type = gaussian, cutoff = 1.0, element_type=CSBEN,item = ETILT, bind = 0, amplitude = 1.0e-4 &end
&error_element name = be2, type = gaussian, cutoff = 1.0, element_type=CSBEN,item = ETILT, bind = 0, amplitude = 1.0e-4 &end

&correct
    disable = 0,
    mode = "orbit", method = "global",
    verbose = 1, n_iterations = 4, n_xy_cycles = 100,
    correction_fraction[0] = 0.5, 0.5,
    bpm_noise[0] = 0e-6, 0e-6,
    correction_accuracy[0] = 1e-6, 1e-6,
    trajectory_output = "%s.orb",
    corrector_output = "%s.cor",
    statistics = "%s.scor",
    closed_orbit_iteration_fraction = 0.5
    closed_orbit_iterations = 500,
    closed_orbit_accuracy = 1e-9,
    prezero_correctors = 0,
    keep_largest_SVs[0] = 0, 0
&end


&closed_orbit
    output = %s.clo,
    closed_orbit_accuracy = 1e-12,
    closed_orbit_iterations = 500,
    iteration_fraction = 0.5 
&end

&find_aperture
        output = %s.aper
        mode = "n-line",
        n_lines = 21,
        xmax = 0.03
        ymax = 0.008
        nx = 17
        n_splits = 2
        split_fraction =0.1
&end