!============================================================
!Tracking file
&change_particle
  name = electron
&end

&divide_elements name=*, type=CSBEND, divisions=5 &end

!&divide_elements name=*, type=KQUAD, divisions=5 &end
!&divide_elements name=*, type=KSEXT, divisions=5 &end

&run_setup
    lattice       = SR_cooler_v4.lte,
    use_beamline  = complete_ring,
!    p_central_mev = 54.997626,  ! 149.99912960, !54.997626
    p_central     = 303.3258147, !107.6276801, 
!    default_order = 3,
!    concat_order  = 3,
!    centroid      = "%s.cen",
!    output        = "beam.out",
!    sigma         = "%s.sig",
!    always_change_p0 = 1,
&end


!uncomment if sr effect is introduced while tracking
&twiss_output
   filename             = "%s.twi",
   concat_order         = 2,	
   matched              = 1,
!  cavities_are_drifts_if_matched=1,
   output_at_each_step  = 1, 
!  radiation_integrals  = 1,
   higher_order_chromaticity = 1,
   higher_order_chromaticity_range  = 3e-6,  
   higher_order_chromaticity_points = 7,
   compute_driving_terms     = 1, 
&end


&alter_elements name=*, type=CSBEND, item=N_KICKS, value=30, allow_missing_elements=1 &end
&alter_elements name=*, type=KSEXT, item=N_KICKS, value=25, allow_missing_elements=1 &end
&alter_elements name=*, type=KQUAD, item=N_KICKS, value=30, allow_missing_elements=1 &end

!&alter_elements name=B*, type=CSBEND, item=SYNCH_RAD, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=ISR, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=USE_RAD_DIST, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=ADD_OPENING_ANGLE, value=1 &end

!&alter_elements name=Q*, type=KQUAD, item=SYNCH_RAD, value=1 &end
!&alter_elements name=Q*, type=KQUAD, item=ISR, value=1 &end

!&alter_elements name=S*, type=*SEXT*, item=SYNCH_RAD, value=1 &end
!&alter_elements name="*", type=*SEXT*, item=K2, value=0, differential=<wsext>, &end

!set momentum offset
!&alter_elements name=MA, item=DP, value="(<offset> 100 /)" &end
!&alter_elements name=MA, item=DX, value=-2e-3 &end

&run_control
!    n_indices = 1
    n_steps = 1
    n_passes =2
&end

!&vary_element 
!  index_number = 0, 
!  index_limit = 1,      
!  name = ma, 
!  item = dp, 
!  initial = 0,
!  final = 0,
!  multiplicative = 1
!&

! &sdds_beam input=beam_in.sdds &end

&bunched_beam n_particles_per_bunch = 1 &end

!&bunched_beam n_particles_per_bunch = 1, centroid[0]=0.0,0.0,0.0,0.0,0.0,0.0 &end

&track &end
&stop &end






!=====================================================
!Twiss file

&change_particle
  name = electron
&end

&divide_elements name=*, type=CSBEND, divisions=5 &end

!&divide_elements name=*, type=KQUAD, divisions=5 &end
!&divide_elements name=*, type=KSEXT, divisions=5 &end

&run_setup
    lattice       = SR_cooler_v4.lte,
!    lattice      = "<fn>",
    use_beamline  = complete_ring_twiss,
!    p_central_mev = 54.997626,
    p_central     = 107.6276801,
    default_order = 3,
    concat_order  = 2,
!    centroid     = %s.out,
!    acceptance   = %s.acc,
!    final        = %s.fin,		
    parameters    = %s.par,
    magnets       = "%s.mag",
!    rootname=<rootname>
&end

!there is no need to set synch_rad=1 for twiss parameter output

&run_control
!        n_steps  = 1,
        n_passes = 128
&end

&twiss_output
   filename             = "%s.twi",
   concat_order         = 2,	
   matched              = 1,
   output_at_each_step  = 1, 
   radiation_integrals  = 1,
   higher_order_chromaticity = 1,
   higher_order_chromaticity_range  = 3e-6,  
   higher_order_chromaticity_points = 7,
   compute_driving_terms     = 1, 
&end

&bunched_beam n_particles_per_bunch = 1 &end
&track &end

&stop &end






















!===================================================
! with sr element in

&change_particle
  name = electron
&end

&divide_elements name=*, type=CSBEND, divisions=5 &end

!&divide_elements name=*, type=KQUAD, divisions=5 &end
!&divide_elements name=*, type=KSEXT, divisions=5 &end

&run_setup
    lattice       = SR_cooler_v4.lte,
    use_beamline  = complete_ring,
    p_central_mev = 54.997626, !154.999158
    default_order = 3,
    concat_order  = 2,
!  centroid      = "%s.cen"
!    sigma         = "%s.sig"
&end


&twiss_output
   filename             = "%s.twi",
   concat_order         = 2,	
   matched              = 1,
  output_at_each_step   = 1,
  radiation_integrals   = 1,
  higher_order_chromaticity = 1,
  higher_order_chromaticity_range  = 3e-6,  
  higher_order_chromaticity_points = 7,
  compute_driving_terms     = 1, 
&end

&alter_elements name=*, type=CSBEND, item=N_KICKS, value=25, allow_missing_elements=1 &end
!&alter_elements name=*, type=KSEXT, item=N_KICKS, value=2, allow_missing_elements=1 &end
&alter_elements name=*, type=KQUAD, item=N_KICKS, value=25, allow_missing_elements=1 &end

!&alter_elements name=B*, type=CSBEND, item=SYNCH_RAD, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=ISR, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=USE_RAD_DIST, value=1 &end
!&alter_elements name=B*, type=CSBEND, item=ADD_OPENING_ANGLE, value=1 &end

!&alter_elements name=Q*, type=KQUAD, item=SYNCH_RAD, value=1 &end
!&alter_elements name=Q*, type=KQUAD, item=ISR, value=1 &end

!&alter_elements name=S*, type=*SEXT*, item=SYNCH_RAD, value=1 &end
!&alter_elements name="*", type=*SEXT*, item=K2, value=0, differential=<wsext>, &end

! set momentum offset
!&alter_elements name=MA, item=DP, value="(<offset> 100 /)" &end
!&alter_elements name=MA, item=DX, value=-2e-3 &end

&run_control
    n_indices = 1
    n_steps = 1
    n_passes =500
&end

&vary_element 
  index_number = 0, 
  index_limit = 2,      
  name = ma, 
  item = dp, 
  initial = 0,
  final = 1,
  multiplicative = 1
&

&bunched_beam 
n_particles_per_bunch = 1 centroid[0]=0.0,0.0,0.0,0.0,0.0,0.0 
&end


&track &end

&stop &end

&matrix_output
  printout = "%s.mpr"
  printout_order = 1
&end