! Optimize the beamline by tracking.
! The beam is from btsScat03, which includes the effect of the
! scattering elements at the start of the BTS.

&run_setup
	lattice = btsScat.lte,
	use_beamline = BTSSCAT1,
	default_order = 1,
	magnets = %s.mag,
	final = %s.finOpt,
	parameters= %s.paramOpt,
	p_central = 5.8708e3   ! 3 GeV
&end

&run_control
	n_steps = 1
&end

&optimization_setup
	target = 1e-6,
	tolerance = 1e-12,
	log_file = /dev/tty,
	verbose = 1,
&end

&optimization_term
	term = "Sx Sxp / 4.0 - sqr" &end
&optimization_term
	term = "s12 Sx / Sxp / .139 - sqr" &end
&optimization_term
	term = "Sy Syp / 5.0 - sqr 0.25 *" &end
&optimization_term
	term = "s34 abs Sy / Syp / .287 - sqr 0.25 *" &end
&optimization_term
	term = "Transmission 1 - sqr 200 *" &end

&optimization_variable
	name = Q1F, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q2D, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q3F, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q4D, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q5D, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q6F, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q7D, item=K1, lower_limit=-1.9, upper_limit=1.9,
	step_size = 0.01 &end
&optimization_variable
	name = Q8, item=K1, lower_limit=-1.0, upper_limit=1.0,
	step_size = 0.01 &end
&optimization_variable
	name = Q9, item=K1, lower_limit=-1.0, upper_limit=1.0,
	step_size = 0.01 &end
	
&sdds_beam
	input = btsScat03.out,
	reuse_bunch = 1,
&end

&optimize &end