! Orbit correction
! Quadrupoles: displacement(0.2mm); tilt(0.5mrad); strength(0.1%)
! Bend: displacement(0.2mm); tilt(0.5mrad); strength(0.1%)

&change_particle
		name="proton"
&end

&run_setup
		lattice = "beamline.lte",
		use_beamline = "PTF_F1",
		default_order = 2,
		print_statistics= 1,
		random_number_seed = 0,
		parameters = %s.param,
    	p_central = 0.7771,
		tracking_updates = 1,
		echo_lattice = 1,
		centroid = %s.cen,
&end

&twiss_output
		filename = "%s.twi",
		matched = 0,
		output_at_each_step= 1,
		final_values_only = 0,
		statistics=1,
		radiation_integrals=0,
		concat_order = 3,
		beta_x = 1.25,
		alpha_x = 0,
		eta_x = 0,
		etap_x = 0,
		beta_y = 0.2,
		alpha_y = 0,
		eta_y = 0,
		etap_y = 0,
&end

&run_control
        n_steps = 1000,
&end

&error_control
		clear_error_settings = 1,
		error_log = "%s.erl",
&end

&error_element
		name = *,
		element_type = KQUAD,
		item = "dx",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = KQUAD
		item = "dy",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = KQUAD
		item = "dz",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = KQUAD
		item = "FSE",
		type = "uniform",
		amplitude = 1e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = KQUAD
		item = "tilt",
		type = "uniform",
		amplitude = 0.5e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = CSBEND,
		item = "dx",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = CSBEND
		item = "dy",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = CSBEND,
		item = "dz",
		type = "uniform",
		amplitude = 0.2e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = CSBEND
		item = "FSE",
		type = "uniform",
		amplitude = 1e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = CSBEND
		item = "etilt",
		type = "uniform",
		amplitude = 0.5e-3,
		post_correction = 0,
		fractional = 0,
&end

&error_element
		name = *,
		element_type = MONI,
		item = "dx",
		type = "uniform",
		amplitude = 0.3e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = MONI,
		item = "dy",
		type = "uniform",
		amplitude = 0.3e-3,
		post_correction = 0,
		fractional = 0,
&end
&error_element
		name = *,
		element_type = MONI,
		item = "tilt",
		type = "uniform",
		amplitude = 5e-3,
		post_correction = 0,
		fractional = 0,
&end

&error_element
		name = *,
		element_type = HKICK,
		item = "tilt",
		type = "uniform",
		amplitude = 2e-3,
		post_correction = 0,
		fractional = 0,
&end

&error_element
		name = *,
		element_type = VKICK,
		item = "tilt",
		type = "uniform",
		amplitude = 2e-3,
		post_correction = 0,
		fractional = 0,
&end

&error_control
		summarize_error_settings = 1,
&end

&global_settings
        mpi_io_write_buffer_size = 1000000
&end

&correct
		mode = "trajectory",
		method = "global",
		trajectory_output = "%s.traj",
		corrector_output = "%s.cor",
		statistics = "%s.scor",
!		corrector_tweek =
		corrector_limit[0] = 5.0E-3, 5.0E-3
		correction_fraction[0] = 0.5, 0.5
		correction_accuracy[0] = 1e-6, 1e-6
		bpm_noise[0] = 0e-6, 0e-6,
		verbose = 1,
		n_iterations = 5,
		n_xy_cycles = 2,
		prezero_correctors = 1,
		track_before_and_after = 1,
&end

&bunched_beam
		bunch = "%s.bun",
		emit_x  = 1.6e-8;
		beta_x  = 1.25;
		alpha_x = 0.0;
		eta_x   = 0.0;
		etap_x  = 0.0;
		emit_y  = 1.6e-8;
		beta_y  = 0.2;
		alpha_y = 0.0;
		eta_y   = 0.0;
		etap_y  = 0.0;
		Po = 0.7771,
		n_particles_per_bunch = 1,
		distribution_cutoff[0] = 3*3,
        enforce_rms_values[0] = 3*1
&end

&track &end