&divide_elements name=*, type=CSBEND, divisions=50 &end
&run_setup
    	lattice = "STORAGE.lat",
    	magnets = %s.mag,
    	parameters = %s.param,
    	final = %s.fin,
	losses = %s.los,
	centroid=%s.cen,
	output = %s.out,
        sigma = %s.sig
    	!random_number_seed = 1985,
   	p_central_mev = 2500,
    	use_beamline="MS",
    	default_order = 3,
    	concat_order = 2,
        !concat_order = 3,
	element_divisions = 10
        rootname = track,
        !!! Set reference momentum for the beam used in optimization
        !expand_for = optimizedQ.bun
&end

&twiss_output 
	filename = %s.twi,
	matched = 1,
	statistics=1,
	!radiation_integrals=1,
	beta_x = 13.5982e+00;  
        alpha_x = 7.8005e-05;  
    	eta_x =0.53;  
    	beta_y = 1.651e+00;  
    	alpha_y = -1.9316e-07;  
        eta_y = 0.0;  
        output_at_each_step = 1
&end

!!! Enable the WATCH elements
&alter_elements name=*, type=WATCH, item=DISABLE, value=0 &end

&run_control	 
        n_steps = 1
&end

&optimization_setup
	tolerance = 1e-26,
	target = 40000,
	method = "simplex",
	mode = "minimize",
	n_passes = 3,
	n_restarts = 10,
	n_evaluations =1500,
	verbose = 1
        log_file = "quads.sdds"
	output_sparsing_factor = 50
&end

&optimization_variable  name=QD1, item=K1, lower_limit=-3.0, upper_limit=0, step_size=0.001 &end
&optimization_variable  name=QD2, item=K1, lower_limit=-3.0, upper_limit=0, step_size=0.001 &end
&optimization_variable  name=QD3, item=K1, lower_limit=-3.0, upper_limit=0, step_size=0.001 &end


&optimization_variable  name=QF1, item=K1, lower_limit=0, upper_limit=3.0, step_size=0.001 &end
&optimization_variable  name=QF2, item=K1, lower_limit=0, upper_limit=3.0, step_size=0.001 &end
&optimization_variable  name=QF3, item=K1, lower_limit=0, upper_limit=3.0, step_size=0.001 &end

&optimization_term term = "IN_WIGG#1.alphax 0 1e-6 sene" &end
&optimization_term term = "IN_WIGG#1.alphay 0 1e-6 sene" &end
&optimization_term term = "OUT_WIGG#1.alphax 0 1e-6 sene" &end
&optimization_term term = "OUT_WIGG#1.alphay 0 1e-6 sene" &end
!&optimization_term term = "IN_WIGG#1.etapx 0 1e-6 sene" &end
!&optimization_term term = "OUT_WIGG#1.etapx 0 1e-6 sene" &end
!&optimization_term term = "max.betax 15 1e-3 segt" &end
!&optimization_term term = "max.etax 0.6 1e-3 segt" &end
!&optimization_term term = "max.etay 0.6 1e-3 segt" &end

&optimization_term
 term = "M2#1.nux M1#1.nux 7.230000 1e-5 sene", 
 !term = "M2#1.nux M1#1.nux - abs 0 1e-6 sene ", weight=1.0
&end

&optimization_term
 term = "M2#1.nuy M1#1.nuy  6.1900000 1e-5 sene", 
 !term = "M2#1.nuy M1#1.nuy - abs 0 1e-6 sene ", weight=1.0
&end

&bunched_beam
        bunch= %s.bun,
        n_particles_per_bunch = 1000,
        distribution_type[0]     = "gaussian",
        distribution_type[1]     = "gaussian",
        distribution_type[2]     = "gaussian",
        distribution_cutoff[0]   = 3.0,
        distribution_cutoff[1]   = 3.0,
        distribution_cutoff[2]   = 3.0,
        emit_x = 25e-09,
        emit_y =0.25e-09,
        sigma_dp = 1e-04
        sigma_s = 11e-4,
        !one_random_bunch=0,
        beta_x = 13.5982e+00, alpha_x = 0,
        beta_y = 1.651e+00, alpha_y = 0,
        !momentum_chirp = 5,
        symmetrize = 1,
        enforce_rms_values[0] = 1,1,1,
     
&end


&optimize
	summarize_setup = 1
&end

&save_lattice
	filename="%s.new"
&end