&tao_design_lattice
 n_universes = 11
 design_lattice(1:11)%file = 
  'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell'
  'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell' 'cell.bmad@cell'
  'cell.bmad@cell'
/

&tao_template_plot
 plot%name = 'all_orbits'
 plot%n_graph = 1
 plot%x_axis_type = 's'
 plot%autoscale_y = t
/

&tao_template_graph
 graph_index=1
 graph%n_curve=11
 curve(1:11)%data_type = 'orbit.x' 'orbit.x' 'orbit.x' 'orbit.x' 'orbit.x' 
  'orbit.x' 'orbit.x' 'orbit.x' 'orbit.x' 'orbit.x' 'orbit.x'
 curve(1:11)%ix_universe = 1 2 3 4 5 6 7 8 9 10 11
 curve(1:11)%legend_text = '10 MeV' '11 MeV' '12 MeV' '13 MeV' '14 MeV' 
  '15 MeV' '16 MeV' '17 MeV' '18 MeV' '19 MeV' '20 MeV'
/

&tao_template_plot
 plot%name = 'tune'
 plot%n_graph = 1
 plot%x_axis_type = 'lat'
 plot%n_graph = 1
 plot%x%min = -0.33381762648506844
 plot%x%max = 0.3336721723516679
/

&tao_template_graph
 graph_index=1
 graph%n_curve=2
 curve(1:2)%data_type_x='beam_start[pz]' 'beam_start[pz]'
 curve(1:2)%data_type='tune.a[0]/twopi' 'tune.b[0]/twopi'
/

&tao_var
 v1_var%name = 'b0'
 ix_min_var = 1
 ix_max_var = 2
 var(1:2)%ele_name = 'qf' 'qd'
 var(1:2)%attribute = 'b_field_err' 'b_field_err'
 var(1:2)%step = 1e-4 1e-4
/

&tao_var
 v1_var%name = 'b1'
 ix_min_var = 1
 ix_max_var = 2
 var(1:2)%ele_name = 'qf' 'qd'
 var(1:2)%attribute = 'b1_gradient' 'b1_gradient'
 var(1:2)%step = 1e-3 1e-3
/

&tao_d2_data
 d2_data%name = 'fit'
 universe = "1,11"
 n_d1_data = 2
/

&tao_d1_data
 ix_d1_data = 1
 d1_data%name = 'tune'
 default_weight = 1.0
 ix_min_data = 1 
 ix_max_data = 2
 datum(1)%data_type = 'expression: lat::r.11[beginning&end]+r.22[beginning&end]'
 datum(2)%data_type = 'expression: lat::r.33[beginning&end]+r.44[beginning&end]'
/

&tao_d1_data
 ix_d1_data = 2
 d1_data%name = 'orbit'
 default_weight = 1.0
 ix_min_data = 1 
 ix_max_data = 2
 datum(1:2)%data_type = 'orbit.x' 'orbit.x'
 datum(1:2)%ele_start_name = 'beginning' 'beginning'
 datum(1:2)%ele_name = 'end' 'end'
 datum(1:2)%meas = 0.0 0.0
 datum(1:2)%merit_type = 'min' 'max'
/

&tao_d2_data
 d2_data%name = 'fit2'
 universe = "11"
 n_d1_data = 1
/

&tao_d1_data
 ix_d1_data = 1
 d1_data%name = 'center'
 default_weight = 1.0
 ix_min_data = 1
 ix_max_data = 1
 datum(1)%data_type = 'expression: 0.5*(1@data::fit.orbit[1]+11@data::fit.orbit[2])'
 datum(1)%meas = 0.0
 datum(1)%good_user = f
/