& init_params_namelist

   ! lattice file path relative to bmad directory (max. 200 characters)

   init_params%latfile        = 'src/lattice/bds-noCC-noSol.bmad'

!    init_params%latfile        = 'src/lattice/bds-noCC-ILD-simple.bmad'

!    init_params%latfile        = 'src/lattice/bds-noCC-SiD-simple.bmad'



   ! initial beam file path relative to bmad directory (max. 200 characters)

   ! if beamfile = '', init_beam_distribution is used

   init_params%beamfile       = ''

!    init_params%beamfile       = '/tmp/mani36c.beam'



   init_params%elestart       = -1        ! first element to be tracked through (-1 = from marker "IP")

   init_params%elemax         = -2       ! last element to be tracked through (-1 = track till the bitter end; -2 = till marker "MEXFOC" (incl.); -3 = till marker "IP" (excl.))

   init_params%human_readable = .false.  ! switch output between human/ROOT-script readable

   init_params%verbiage       = 1        ! 0: pol,IP only   1: leave out markers etc.

                                        ! 2: (out of order)   3: all elements



   init_params%sx_offset     = 2.e-6      ! random (Gaussian) magnet misalignment [m]

   init_params%sy_offset     = 2.e-6      ! random (Gaussian) magnet misalignment [m]

   init_params%ss_offset     = 0.         ! random (Gaussian) magnet misalignment [m]

   init_params%sx_pitch      = 0.       ! random (Gaussian) magnet misalignment [m]

   init_params%sy_pitch      = 0.       ! random (Gaussian) magnet misalignment [m]

   init_params%stilt         = 0.       ! random (Gaussian) magnet misalignment [m]

   init_params%sroll         = 0.       ! random (Gaussian) magnet misalignment [m]



   init_params%FFstart        = 1050     ! start of final focus region (element number): misalignments multiplied by 0.001

   init_params%FFend          = 1092     ! end of final focus region (element number)     set to -1 to be ignored



   init_params%randomseed     = 7593589       ! random seed (<0: generate from system time)

   init_params%charge         = -1.      ! particle charge [e]  (i.e. electron charge = -1)



   ! spin initialization only, if beamfile = ''

   init_params%polarization   = 0.8      ! beam polarization

   init_params%spin_model     = 3        ! 0: all same spin   1: Gaussian distr.   2: mix from P=+-1 (exact)   3: mix from P=+-1

!    init_params%sigma0         = 0.005    ! Width of polarization distribution (spin_model == 1)

/



&initial_values

   ! parameters from ILC RDR 3, 2.6.2

   ! init_beam_distribution is used only, if beamfile = ''

   beam_init%n_bunch = 1                     ! Number of bunches.

   beam_init%n_particle = 10000              ! Number of simulated particles per bunch.

   

   beam_init%a_norm_emitt = 9.4e-6           ! a-mode emittance [m?]

   beam_init%b_norm_emitt = 34.e-9           ! b-mode emittance [m?]

   beam_init%sig_z = 0.3e-3                  ! Z sigma [m]

   beam_init%sig_e = 0.14e-2                 ! e_sigma in dE/E.

   beam_init%bunch_charge = -2.E10           ! charge in a bunch. [e]    (conversion to Coulomb internally)

   beam_init%dt_bunch = 369.2E-9             ! Time between bunches.

/















!   beam_init%dPz_dz = 0                      ! Correlation of Pz with long position.

!   beam_init%center(6) = 0                   ! Bunch center offset relative to reference.

!   beam_init%dt_bunch                        ! Time between bunches.

!   beam_init%center_jitter(6) = 0.0          ! Bunch center rms jitter

!   beam_init%emitt_jitter(2)  = 0.0          ! a and b bunch emittance rms jitter normalized to emittance

!   beam_init%sig_z_jitter     = 0.0          ! bunch length RMS jitter

!   beam_init%sig_e_jitter     = 0.0          ! energy spread RMS jitter

!   beam_init%init_spin = .false.             ! initialize beam spinors

      !   type(beam_spin_struct)  spin        ! Initialize the spin



! Note: Make sure: |beam_init%dpz_dz| < mode%sigE_E / mode%sig_z

      

      ! type beam_init_struct

      !   real(rp) a_norm_emitt     ! a-mode emittance

      !   real(rp) b_norm_emitt     ! b-mode emittance

      !   real(rp) :: dPz_dz = 0    ! Correlation of Pz with long position.

      !   real(rp) :: center(6) = 0 ! Bench center offset relative to reference.

      !   real(rp) dt_bunch         ! Time between bunches.

      !   real(rp) sig_z            ! Z sigma in m.

      !   real(rp) sig_e            ! e_sigma in dE/E.

      !   real(rp) bunch_charge     ! charge in a bunch.

      !   real(rp) :: center_jitter(6) = 0.0 ! Bunch center rms jitter

      !   real(rp) :: emitt_jitter(2)  = 0.0 ! a and b bunch emittance rms jitter normalized to emittance

      !   real(rp) :: sig_z_jitter     = 0.0 ! bunch length RMS jitter

      !   real(rp) :: sig_e_jitter     = 0.0 ! energy spread RMS jitter

      !   type(beam_spin_struct)  spin       ! Initialize the spin

      !   integer :: n_particle = 0          ! Number of simulated particles per bunch.

      !   integer :: n_bunch = 1             ! Number of bunches.

      !   logical :: renorm_center = .true.  ! Renormalize centroid?

      !   logical :: renorm_sigma = .true.   ! Renormalize sigma?

      !   logical :: init_spin     = .false. ! initialize beam spinors

      !   logical :: is_random = .true.      ! Random distribution?  Or a tail-weighted distribution?

      !   character(16) :: random_engine = 'pseudo' ! Or 'quasi'. Random number engine to use.

      !   character(16) :: random_gauss_converter = 'exact'

      !                                             ! Or 'limited'. Uniform to gauss conversion method.

      !   real(rp) :: random_sigma_cutoff = 4.0     ! Used with 'limited' converter. Cut-off in sigmas.

      ! 

      !   type (tail_weighted_beam_init_struct) tw_beam_init ! Parameters for a tail-weighted beam distribution

      ! end type

      ! 

      ! type beam_spin_struct

      !   real(rp) :: polarization = 1.0 ! i.e. 80% polarized

      !   real(rp) :: theta = 0.0  ! polar coordinates

      !   real(rp) :: phi = 0.0    ! polar coordinates

      ! end type



      ! sqrt(0.5) = 0.7071067812