! Optimization run to reduce the emittance of the PAR
! by adjusting the quadrupoles.

&run_setup
    lattice = "parTrack.lte",
    magnets = "%s.mag",
    p_central =1565.5577,
    use_beamline="Ring",
    final = %s.fin,
    losses = %s.los,
    default_order = 2
&end

&twiss_output 
        radiation_integrals=1,
        output_at_each_step=1
&end

&run_control
        n_passes=50,  
        n_steps = 10
&end

&optimization_setup
	equation = "Particles",
	mode = "maximize", method = "randomwalk",
	tolerance = 1e-10, n_passes = 2, n_evaluations = 15, log_file = "tty",
        n_restarts = 6, verbose = 0
&end

&optimization_variable
        name = q1, item=k1, 
        lower_limit=-8.0, upper_limit=8.0, step_size=0.001 
&end
&optimization_variable
        name = q2, item=k1, 
        lower_limit=-8.0, upper_limit=8.0, step_size=0.001 
&end
&optimization_variable
        name = q3, item=k1, 
        lower_limit=-8.0, upper_limit=8.0, step_size=0.001 
&end

&chromaticity
    sextupoles = "Sf1 Sd1"
    change_defined_values = 1,
    dnux_dp = 0.1,
    dnuy_dp = 0.1,
    tolerance = 0.1
&end

&bunched_beam
        n_particles_per_bunch = 1000,
        bunch = %s.bun,
        distribution_type[0] = "dynamic-aperture", "dynamic-aperture",
        distribution_cutoff[0] = 40, 20,
        beta_x = 1e6, emit_x = 50e-10,
        beta_y = 1e6, emit_y = 40e-10
&end

&optimize summarize_setup=1 &end


&track &end