&run_setup
   lattice       = "BubbleChamber.lte",
   use_beamline  = "BUBBLE_CHAMBER",
   p_central_mev = 9.0,
   centroid      = %s.cen,
   sigma         = %s.sig,
   final         = %s.fin,
   parameters    = %s.par,
   magnets       = %s.mgn,
   default_order = 2,
&end

&run_control 
   n_steps=1
&end

&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 8.67, alpha_x = -2.61,
   beta_y = 6.41, alpha_y = -1.74,
&end

&optimization_setup
  balance_terms = 1,
  tolerance = 1E-60,
  n_restarts = 10,
  n_passes = 3,
  n_evaluations = 1000,
  target = 1E-16,
  output_sparsing_factor = 30
&end

&optimization_variable
  name = MQJ0L02,  item=K1,  lower_limit=-15, upper_limit=15, step_size=1E-03
&end

&optimization_variable
  name = MQJ0L02A, item=K1,  lower_limit=-15, upper_limit=15, step_size=1E-03
&end

&optimization_variable
  name = MQD5D00,  item=K1,  lower_limit=-15, upper_limit=15, step_size=1E-03
&end

&optimization_variable
  name = MQD5D01,  item=K1,  lower_limit=-15, upper_limit=15, step_size=1E-03
&end

&optimization_term
 weight=1,
 term="RADIATOR#1.Sx 5E-04 - sqr",
 verbose = 1,
&end

&optimization_term
 weight=1,
 term="RADIATOR#1.Sy 5E-04 - sqr",
 verbose = 1,
&end

&bunched_beam
   n_particles_per_bunch = 5000,
   one_random_bunch = 1,
   emit_nx = 0.75e-06,
   emit_ny = 1.25e-06,
   use_twiss_command_values = 1,
   !
   sigma_dp = 0.000546,    
   sigma_s =  0.000920,
   !
   distribution_type[0]   = 3*"gaussian"
   distribution_cutoff[0] = 3*3,
   enforce_rms_values[0]  = 3*1,
&end

!&sdds_beam
    input="6.3MeVc_20fC.sdds",
!&end

&optimize
  summarize_setup=1
&end

&save_lattice
  filename=NEW.lte,
&end