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

&floor_coordinates
   filename="%s.flr",
   X0 =  80.6,
   Z0 = -249.14700, ! exit of quarter cryounit
   Y0 =  100.0,
   magnet_centers = 1
&end

!&floor_coordinates
   filename="%s.flr",
   X0 =  0,
   Z0 =  0, ! exit of quarter cryounit
   Y0 =  0,
   magnet_centers = 1
!&end

&run_control 
   n_steps=1
&end

!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 11.75, alpha_x = 25.00,
   beta_y = 44.78, alpha_y = -4.03,
!&end

! 8.3 MeV/c
!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 2.29, alpha_x =-0.092,
   beta_y = 2.04, alpha_y =-0.019,!
!&end

! 6.3 MeV/c
!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 3.12, alpha_x =-0.57,
   beta_y = 2.83, alpha_y =-0.40,
!&end

! 6.3 MeV/c FROM FIT
&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 22.96, alpha_x = 10.14,
   beta_y = 29.69, alpha_y = -4.56,
&end

! 5.5 MeV/c
!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 4.91, alpha_x =-1.22,
   beta_y = 4.62, alpha_y =-1.00,
!&end

! 4.2 MeV/c
!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 9.26, alpha_x =-2.85,
   beta_y = 9.06, alpha_y =-2.64,
!&end

! 3.2 MeV/c
!&twiss_output
   filename = "%s.twi",
   matched=0,radiation_integrals=1,
   beta_x = 8.12, alpha_x =-3.01,
   beta_y = 7.86, alpha_y =-2.77,
!&end

&matrix_output
    !SDDS_output="%s.mpr",
    !SDDS_output_order=2,
    printout="%s.mpr",
    printout_order=2,
&end

&bunched_beam
   n_particles_per_bunch = 25000,
   one_random_bunch = 1,
   emit_nx = 1.00e-07,
   emit_ny = 1.00e-07,
   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

&track
&end