&run_setup
        lattice = "lattice.lte",
        default_order = 2,
        use_beamline = "L1OtoL3FS5",
        final = %s.fin,
        parameters = %s.param,
        always_change_p0 = 1,
        p_central = 7.809589e+01
&end

&run_control
n_steps = 1,
reset_rf_for_each_step = 1
&end

&optimization_setup
	mode = "minimize"
	method = "simplex"
	n_evaluations = 1500
	n_passes = 3
	n_restarts = 2
	log_file = /dev/tty
	verbose = 0
&end

&optimization_term term = "St 1e15 *" &end
&optimization_term term = "150 mev / pCentral 1 sene" &end

&optimization_variable name=L2CELL1, item=PHASE,  lower_limit=-180, upper_limit=180, step_size=0.1 &end
&optimization_covariable name=L2CELL12, item=PHASE, equation=L2CELL1.PHASE &end
&optimization_covariable name=L2CELL14, item=PHASE, equation=L2CELL1.PHASE &end
&optimization_covariable name=L2CELL16, item=PHASE, equation=L2CELL1.PHASE &end
&optimization_covariable name=L2CELL2, item=PHASE, equation=L2CELL1.PHASE &end
&optimization_covariable name=L2CELL4, item=PHASE, equation=L2CELL1.PHASE &end
&optimization_covariable name=L2CELL8, item=PHASE, equation=L2CELL1.PHASE &end

&optimization_variable name=L2CELL1, item=VOLT,  lower_limit=0, upper_limit=3e6, step_size=1e3 &end
&optimization_covariable name=L2CELL12, item=VOLT, equation=L2CELL1.VOLT &end
&optimization_covariable name=L2CELL14, item=VOLT, equation=L2CELL1.VOLT &end
&optimization_covariable name=L2CELL16, item=VOLT, equation=L2CELL1.VOLT &end
&optimization_covariable name=L2CELL2, item=VOLT, equation=L2CELL1.VOLT &end
&optimization_covariable name=L2CELL4, item=VOLT, equation=L2CELL1.VOLT &end
&optimization_covariable name=L2CELL8, item=VOLT, equation=L2CELL1.VOLT &end

&bunched_beam
n_particles_per_bunch = 1000,
emit_nx = 4.601407e-7,
emit_ny = 4.625478e-7,
beta_x = 1.614258e1,
alpha_x = -2.660602,
beta_y = 1.627933e1,
alpha_y = -2.684884,
sigma_dp = 4.357766e-3,
sigma_s = 1.263005e-3,
distribution_type[0] = "gaussian",
distribution_type[1] = "gaussian",
distribution_type[2] = "gaussian",
distribution_cutoff[0] = 2.0,
distribution_cutoff[1] = 2.0,
distribution_cutoff[2] = 2.0,
enforce_rms_values[0] = 1,
enforce_rms_values[1] = 1,
enforce_rms_values[2] = 1,
one_random_bunch = 1,
bunch = %s.bun
&end
 
&optimize &end

&save_lattice filename = %s.new &end

&run_setup
        lattice = optimize.new,
        default_order = 2,
        use_beamline = "L1OtoL3FS5",
        output = %s.out,
	sigma = %s.sig,
        magnets = %s.mag,
        always_change_p0 = 1,
        p_central = 7.809589e+01
&end

&run_control
&end

&sdds_beam
	input = optimize.bun,
&end

&track &end