&parameters 
  lat_file = '$DIST_BASE_DIR/lattice/cesr/bmad/cta_2085mev_20090516.lat'
  !lat_file = "a_lattice.bmad"
  bpm_rotation = 0e-3 ! typically 10.e-3 radians
  bpm_noise = 0e-6 ! typically 10.e-6 meters
  gain_sigma = 0e-2 ! typically 5% = 0.05
  timing_sigma = 0.e-12 ! typically 10.e-12 seconds
  n_turns = 4096
  skq02w_kick = 0.0
  track6x6 = .true.
  n_damping = 100000 ! # turns to damp before measurement

  auto_bookkeeper = .true. ! set to .false. to use intelligent bookkeeping

  init_vec(1) = 3.e-3
  init_vec(3) = 0.e-3

  output = 'cesrv'
  nIter = 1

  ! target tunes for q_tuning. comment out if you want to use lattice default tunes.
  ! set to (integer + fractional) tunes
  target_tunes(1) = 14.5701
  target_tunes(2) =  9.6281
  target_tunes(3) = -0.06588


! tune tracker parameters
n_TTs = 2 ! if n_TTs = 0, sim_cbar will turn off damping, even if 6x6 tracking is on

tt_params(1)%useSaveState = .false.
tt_params(1)%orientation = 'h'
tt_params(1)%kck_name = 'Q10W'
tt_params(1)%bpm_name = 'DET_05W'
tt_params(1)%LPinertia = 32768   ! LP filter time constant is LPinertia*fastPeriod.  2^15 = 32768
tt_params(1)%cyc_per_turn = 183  ! Number of times tune tracker cycles per CESR period.  Usually 183
tt_params(1)%Ki = 0.2 ! 0.2   ! integrator gain
tt_params(1)%Kp = 20.0 ! 20.0   ! proportional gain
tt_params(1)%Kvco = 10.0   ! VCO gain.  Redundant, because total gain looks like Kvco*(Ki*x + Kp*y + Kd*z)
tt_params(1)%modTfrac0 = 0.570  ! 0.567  ! guess of fractional tune.  used to compute w0 of modulator
tt_params(1)%kickAmplitude = 5.0E-8
tt_params(1)%mixmode = 'sin'
  !Parameters specific to D channel
tt_params(1)%use_D_chan = .false.  ! bad?
tt_params(1)%Kd = 0.1     ! differential gain
tt_params(1)%wls_N = 25   ! Number of data points for LS fit.  More DP makes smoother but increases delay
tt_params(1)%wls_order = 1  ! Order of fit polynomial.  Higher order is more accurate but noisier

tt_params(2)%useSaveState = .false.
tt_params(2)%orientation = 'v'
tt_params(2)%kck_name = 'Q10W'
tt_params(2)%bpm_name = 'DET_05W'
tt_params(2)%LPinertia = 32768   ! LP filter time constant is LPinertia*fastPeriod.  2^15 = 32768
tt_params(2)%cyc_per_turn = 183 ! Number of times tune tracker cycles per CESR period.  Usually 183
tt_params(2)%Ki = 0.2   ! integrator gain
tt_params(2)%Kp = 20.  !  20.0    ! proportional gain
tt_params(2)%Kvco = 10.0   ! VCO gain.  Redundant, because total gain looks like Kvco*(Ki*x + Kp*y + Kd*z)
tt_params(2)%modTfrac0 = 0.628 ! 0.628    ! guess of fractional tune.  used to compute w0 of modulator
tt_params(2)%kickAmplitude =  4.e-8 ! 5.0E-8
tt_params(2)%mixmode = 'sin'
  !Parameters specific to D channel
tt_params(2)%use_D_chan = .false.
tt_params(2)%Kd = 0.0010  ! differential gain
tt_params(2)%wls_N = 25   ! Number of data points for LS fit.  More DP makes smoother but increases delay
tt_params(2)%wls_order = 1  ! Order of fit polynomial.


! ma_params_struct parameters
   misalign_magnets = .true.
   seed = -1 ! seed < 1 = seed based on CPU time

! Keys:   sbend=2; quad=3; sextupole=5; wiggler=12
! Params: tilt=2; k1=4; roll=15; x_offset=36; y_offset=37;
!         x_pitch=19, y_pitch=20; s_offset=27;

! Quads
   ma(1)%key_name   = "quadrupole"
   ma(1)%mask       = "^Q"
   ma(1)%param_name = "x_offset"
   ma(1)%amp        = 150.e-6

   ma(2)%key_name   = "quadrupole"
   ma(2)%mask       = "^Q"
   ma(2)%param_name = "y_offset"
   ma(2)%amp        = 150.e-6

   ma(3)%key_name   = "quadrupole"
   ma(3)%mask       = "^Q"
   ma(3)%param_name = "tilt"
   ma(3)%amp        = 100.e-6

   ma(4)%key_name   = "quadrupole"
   ma(4)%mask       = "^Q"
   ma(4)%param_name = "k1"
   ma(4)%amp        = 4.e-4

! Sextupoles
   ma(5)%key_name   = "sextupole"
   ma(5)%mask       = ""
   ma(5)%param_name = "x_offset"
   ma(5)%amp        = 0 !300.e-6

   ma(6)%key_name   = "sextupole"
   ma(6)%mask       = ""
   ma(6)%param_name = "y_offset"
   ma(6)%amp        = 0 !300.e-6

   ma(7)%key_name   = "sextupole"
   ma(7)%mask       = ""
   ma(7)%param_name = "k2"
   ma(7)%amp        = 0 !100.e-6

! Bends
   ma(8)%key_name   = "sbend"
   ma(8)%mask       = ""
   ma(8)%param_name = "x_offset"
   ma(8)%amp        = 0 !150.e-6

   ma(9)%key_name   = "sbend"
   ma(9)%mask       = ""
   ma(9)%param_name = "x_offset"
   ma(9)%amp        = 0 !150.e-6

   ma(10)%key_name   = "sbend"
   ma(10)%mask       = ""
   ma(10)%param_name = "roll"
   ma(10)%amp        = 0 !100.e-6


/