&input ! lat_file_name = 'bmad_short_wide_inf_ap_cor_tilt.' lat_file_name = 'bmad_short_inf_ap_collimator_backleg2.' !nmuons = 100000 nmuons = 10000 nturns = 500 ! number of turns (integer) create_new_distribution = .true. ! create a new phase space distribution off of the target / or not seed = 0 ! random number generator seed. If seed=0, the seed is set by the system clock. If seed>0, random numbers will be from a list with that starting point new_file = 'muons_off_target.dat' ! write the new distribution to this file, if the name is blank no file will be written ! muon_file = 'muons_off_target.dat' ! read the phase space distribution off the target from this file. If blank, use freshly created distribution ! muon_file = 'VDstop_DS_436_12000.dat' ! read the phase space distribution off the target from this file. If blank, use freshly created distribution tdistr = "e989" ! time/longitudinal distribution: "flat", "gaus", "e821", "e989" tlength = 120.e-9 ! base length of injected bunch (doesn't apply to tdistr="e821") tsigma = 0.e-9 ! sigma for tdistr="gaus" (doesn't apply to tdistr={"uniform","e821","e989"}) pzdistr = "gaus" ! momentum distribution: "flat", "gaus" pz = 0.06 ! cutoff of longitudinal momentum distribution pzsigma = 0.0112 ! sigma for pzdistr="gaus" epsdistr = "gaus" ! emittance distribution at target: "delta", "flat", "gaus" epsx = 14.e-7 ! emittance at inflector (m-rad) epsy = 14.e-7 ! emittance at inflector (m-rad) twiss = 1.5 10.0 0. 0.0 0.0 0. 0. 0.7 0.7 0. 0. ! betax, betay, alphax, alphay, etax, etapx, etay, etapy, phix, phiy, gammax, gammay at inflector midpoint twiss_ref = 'end' inf_aperture = "e821" ! inflector aperture: "none", "rect", "e821", "e989" inflector_width = 0.009 inf_end_us = .false. ! include scattering at upstream end of inflector (E821 coils) inf_end_ds = .false. ! include scattering at downstream end of inflector (E821 coils) enerloss = .true. ! include energy loss with scattering in inflector(E821 coils)! inflector_angle = -0.002 inflector_field = 1.002 ! multiplier of Bmagic. 1.02747 * Bmagic = B_wuzeng inflector_field = 1.012747 ! multiplier of Bmagic. 1.02747 * Bmagic = B_wuzeng inflector_field = 1.02 ! multiplier of Bmagic. 1.02747 * Bmagic = B_wuzeng ! initial_offsets = -5.9190E-02 0. 0. 2.5909E-02 0.0 0.0 !18mm inflector field = Wuzeng map, inflector_field = 1.002. If inflector field < 1.002 then trajector outside map ! initial_offsets = -0.048837 0. 0. 2.4804E-02 0.0 0.0 !36mm inflector field = Bmagic ! initial_offsets = -4.6985E-02 0.0 0.0 2.4231E-02 0.0 0.0 !36mm inflector field = Bmagic ! initial_offsets = -2.9035E-02 0.0 0.0 1.7898E-02 0. 0. !36mm inflector field same as Wuzeng map ! initial_offsets = -0.045937 0. 0. 0.021538 0.0 0.0 !18mm inflector field from Wuzeng map initial_offsets = -4.8627E-02 0. 0. 2.2473E-02 0. 0. !18mm,Wuzeng, -3mrad initial_offsets = -4.8452E-02 0. 0. 2.2402E-02 0. 0. !18mm, Wuzeng, -3mrad initial_offsets = -4.7671E-02 0. 0. 2.2134E-02 0. 0. ! initial_offsets = -4.5671E-02 0. 0. 2.2134E-02 0. 0. ! initial_offsets = -5.1985E-02 0.0 0.0 2.4231E-02 0.0 0.0 !36mm inflector field = Bmagic initial_offsets = -5.2011E-02 0.0 0.0 2.3681E-02 0.0 0.0 initial_offsets = -4.3579E-02 0.0 0.0 2.0900E-02 0.0 0.0 opt_incident = .true. inj_matrix_tracking = .true. ! if true, propagate trajectories through injection channel using matrices, otherwise step through field map ring_theta = 0. ! Rotation of ring elements with respect to inflector exit kickerPlates = 989 ! Controls shape of B-field: 1=uniform, 821, 989 kickerCircuit = 989 ! Controls shape of pulse: 1=square kick, 821=RLC, 989=Blumlein kickerFieldType = 2 ! Controls B-field calculation: 1=analytic, 2=mapped kicker_params%kicker_field = 220.e-4 220.e-4 220.e-4 ! Kick strength (Tesla) for 989 kick kicker_params%dtRise = 20.e-9 20.e-9 20.e-9 ! Kick rise (ns) kicker_params%dtFall = 20.e-9 20.e-9 20.e-9 ! Kick fall (ns) ! kicker_params%kicker_field = 215.e-4 215.e-4 215.e-4 ! Kick strength (Tesla) ! kicker_params%kicker_field = 114.5e-4 114.5e-4 114.5e-4 ! Kick strength (Tesla) kicker_params%kick_width = 120.e-9 120.e-9 120.e-9 ! Kick width (seconds) -- doesn't apply to kickerCircuit=821 kicker_tStart = -15.0e-9 -15.0e-9 -15.0e-9 ! If kicker_tstart(i) = -1., automatically calculate kicker start time quad_plate = .false. ! include scatter in Q1 quad plates quadPlates = 1 ! 821, 989 -- controls shape of E-field for Q1_Long quadCircuit = 1 ! 0=off at injection, ramped up with tauRC=5us; 1=storage potential at injection; todo:2=scraping at inj. quadFieldType = 1 ! 1=analytic, 2=mapped quad_params = 0.185 0.185 0.185 0.185 ! field_index(4) (For the time being quads 2,3 and 4 all have the value assigne to 2 twiss_file = .true. ! true => write twiss parameters every cm to "twiss_file.dat" vparam_id = 1 ! 0-nothing, 1-betax, 2-betay, 3-eta, 4-kick, 5-kick_width, 6-initial position offset, 7-initial angle offset is the parameter to vary vparam_min = 5 vparam_max = 25 ! maximum of parameter to vary delta_vparam = 0.5 ! parameter step loop = .false. ! fringe_file = 'injec_fld.dat' 0.5 1 0 !field map for fringe field, grid spacing, type, flag fringe_file = 'ring_inject_bfield3d_wuzheng_extended.dat' 0.5 2 0 !field map for fringe field, grid spacing, type = type of read format, flag ! inflector_file = 'inf_field_alone.dat' 0.5 2 0 !field map for inflector field, grid spacing, type, flag inflector_file = 'uniform' 0.5 4 0 !field map for inflector field, grid spacing, type, flag /