! Load the solution from the optimization.
! Track and compute twiss parameters.

&run_setup
	lattice = btsScat.lte,
	use_beamline = BTSSCAT,
	default_order = 2,
	magnets = %s.mag,
	element_divisions = 10,
	output = %s.out,
	sigma = %s.sig,
	centroid = %s.cen,
	print_statistics = 1,
	p_central = 5.8708e3   ! 3 GeV
&end

&load_parameters
	filename = btsScat04.paramOpt,
	change_defined_values = 1
&end

!&alter_elements name=SSWIN, type=MATTER,
!  		item=ELASTIC, value=0 &end
!&alter_elements name=SSWIN, type=MATTER,
!		item=ENERGY_STRAGGLE, value=1 &end
!&alter_elements name=KAPWIN, type=MATTER,
!  		item=ELASTIC, value=0 &end
!&alter_elements name=KAPWIN, type=MATTER,
!		item=ENERGY_STRAGGLE, value=1 &end
!&alter_elements name=AIRGAP1, type=MATTER,
!  		item=ELASTIC, value=0 &end
!&alter_elements name=AIRGAP1, type=MATTER,
!		item=ENERGY_STRAGGLE, value=1 &end
!&alter_elements name=AIRGAP2, type=MATTER,
!  		item=ELASTIC, value=0 &end
!&alter_elements name=AIRGAP2, type=MATTER,
!		item=ENERGY_STRAGGLE, value=1 &end
!&alter_elements name=HELIUM2, type=MATTER,
!  		item=ELASTIC, value=0 &end
!&alter_elements name=HELIUM2, type=MATTER,
!		item=ENERGY_STRAGGLE, value=1 &end


&twiss_output
	filename = %s.twi,
	matched = 0,
	beta_x=9.196, alpha_x=2.737,
	eta_x=-0.236, etap_x=.03,
	beta_y=2.832, alpha_y=-0.950	
&end

&run_control
	n_steps = 1
&end

! The booster should produce 0.553e-3 energy spread.  I increased this to 
! 1.0e-3 as a wild guess for the increase in energy spread from scattering.	
!	sigma_dp = 0.553e-3,
&bunched_beam
	n_particles_per_bunch = 10000,
	use_twiss_command_values = 1,
	emit_x = 0.170667e-6,
	emit_y = 0.085333e-6,
	sigma_dp = 1.0e-3,
	sigma_s = 0,
	distribution_cutoff[0] = 3*3,
	distribution_type[0] = 3*"gaussian",
&end

&track &end

&save_lattice filename = %s.new &end