&input !configuration - muons ! :character string - BMAD format lattice file lat_file_name = 'bmad_esquad_grid.' !curvilinear_grid ! :integer - Number of turns to track nturns = 100 ! The next several lines refer to the intial muon phase space distribution ! :integer - Number of muons to track nmuons = 10000 ! :logical - Create a new distribution create_new_distribution = T ! : integer - ! Random number seed. If seed =0 then random number set to system clock. If seed > 0 then it will be seed =0 ! :character string - If you choose to create a new distribution it will be saved in new_file. If the name is new_file = 'empty' ! : If you choose to read an existing distribution, here is where the name is written. ! muon_file = 'VDstop_DS_436_12000.dat' ! read the phase space distribution off the target from this file. If bla ! muon_file = 'particles_endm4m5_100.txt' ! muon_file = 'particles_M4M5End_400_mod.txt' ! : Time distribution of the muons, choices are "flat", "gaus", "e821", "e989" tdistr = "e989-P8" ! tdistr = "gaus" ! : Width in time of distribution : (doesn't apply to tdistr="e821") tlength = 120.e-9 ! : standard deviation (width) of time distribution if it is gaussian, (sigma for tdistr="gaus", doesn't apply tsigma = 22.e-9 !appproximates E821 pulse ! : momentum distribution: "flat", "gaus" pzdistr = "gaus" ! : cutoff of longitudinal momentum distribution, (fractional energy offset with respect to reference) pz = 0.01 ! : sigma for pzdistr="gaus" pzsigma = 0.013 ! : emittance distribution at production target: "delta", "flat", "gaus" epsdistr = "gaus" ! : horizontal emittance at inflector assuming no losses in injection channel (m-rad) epsx = 9.3e-6 ! : vertical emittance at inflector assuming no losses in injection channel (m-rad) epsy = 9.3e-6 ! :12-real - betax, betay, alphax, alphay, etax, etapx, etay, etapy, phix, phiy, gammax, gammay at inflector , twiss = 2.0, 10.0, 0., 0.0, 0.0, 0., 0., 0., 0.7, 0.7, 0., 0. ! : Reference point within inflector for twiss parameters specified in previous line. 'end' = at inflector e twiss_ref = 'brynn' !: logical - if true then start tracking the distribution at the inflector exit, if fa start_tracking_at_inflector_exit = F !: logical - if true, write twiss parameters every cm to "twiss_file.dat", and write transfer matrices for a twiss_file = F ! true => write twiss parameters every cm to "twiss_file.dat" !: character - if = 'open' then use input twiss parameters as starting point for propagation ar ! if = 'closed' then use closed ring twiss parameters as starting point for propag ring_twiss = 'closed' ! :character string - inflector aperture: "none", "rect" (rectangular), "e821"('D' shaped aperture of 821 inf_aperture = "e821" ! :real - defines offset of axis of inflector from magic radius. If inflector_width = 0.009, central a inflector_width = 0.009 ! :logical - if true include scattering at upstream end of inflector (E821 coils) inf_end_us = F ! :logical - if true include scattering at downstream end of inflector (E821 coils) inf_end_ds = F ! :logical - if true include energy loss with scattering in inflector(E821 coils)! enerloss = F ! :real - angle of inflector axis with respect to tangent line. The rotation axis is the downstream en ! If angle < 0 the inflector is rotated clockwise about its downstream end inflector_angle = -0.002 ! :6 reals - multiplier of Bmagic. If inflector_field = 1. then the inflector field = -Bmagic. Note th inflector_field = 1.0 ! multiplier of Bmagic. 1.012747 * Bmagic = B_wuzeng initial_offsets_ref = -2.9343E-02 0. 0. 1.6292E-02 0. 0. !Inf = -2mrad, 0mm, -2mrad initial_offsets = -2.9343E-02 0. 0. 1.6292E-02 0. 0. !Inf = -2mrad, 0mm, -2mrad initial_offsets_ref = -4.5563E-02 0. 0. 2.1401E-02 0. 0. !Inf = 0mrad, 0mm initial_offsets = -4.5563E-02 0. 0. 2.1401E-02 0. 0. !Inf = 0mrad, 0mm !for lattice that starts 10cm upstream from IBMS1 initial_offsets_ref= -6.1076E-02 0. 0. 2.5068E-02 0. 0. ! -2mm, -3.4mm, Inf=0.976 initial_offsets = -6.1076E-02 0. 0. 2.5068E-02 0. 0. ! -2mm, -3.4mm, Inf=0.976 inflector_field=0.976 initial_offsets_ref = -7.8630E-02 0. 0. 3.1185E-02 0. 0. ! 0mm,0mrad,inf=0.976 initial_offsets = -7.8630E-02 0. 0. 3.1185E-02 0. 0. ! 0mm,0mrad,inf=0.976 ! for nominal lattice starting 3mm upstream IBMS1 !Run 3 - lattice that starts 3mm upstream of IBMS1 inflector_field = 0.976 ! multiplier of Bmagic. 1.012747 * Bmagic = B_wuzeng initial_offsets_ref = -5.8643E-02 0. 0. 2.5068E-02 0. 0. ! -2mm,-3.4mrad,Inf=0.976, 4mrad, -2.3mm, -3.6mrad initial_offsets = -5.8643E-02 0. 0. 2.5068E-02 0. 0.0 ! -2mm,-3.4mrad,Inf=0.976, 4mrad, -2.3mm, -3.6mrad, crosses zero at s=0.5, 60% at inflector ! initial_offsets = 0. 0. 0. 0. 0. 0. inflector_end_target = 0.0 0. 0. 0.0 0. 0. ! :logical - if true, the initial_offsets so that the trajectory emerges from the inflector on its axis a opt_incident = F ! :logical - true: use matrices to propagate through injeciton channel, false: track using field ma inj_matrix_tracking = F ! if true, propagate trajectories through injection channel usei ! :real - Rotation of ring elements with respect to inflector exit ring_theta = 0. ! :integer - 989-curved plates, 821-flat plates, 1-uniform vertical field kickerPlates = 989 ! Controls shape of B-field: 1=uniform, 821, 989 ! :integer - Controls shape of pulse: 1=square kick, 821=RLC, 989=Blumlein kickerCircuit = 20150328 ! Controls shape of pulse: 1=square kick, 821=RLC, 9 kickerCircuit = 20171113 !20150328 ! 20171113 - measured 11/17 kickerPulseFile = 'v0/traces-Ir.txt' kickerPulseFile = 'kicker_pulse_run3b4_4col.dat' ! kickerPulseFile = 'e821_kicker_pulse.dat' ! :integer - Controls B-field calculation: 1=analytic, 2=mapped kickerFieldType = 2 ! Controls B-field calculation: 1=analytic, 2=mapped ! : 3 reals - Bfield (T) for the three kicker mangets ! kicker_params%kicker_field = 175.2e-4 175.2e-4 175.2e-4 ! Kick strength (Tesla) for 989 kick 60% (292) - kicker_params%kicker_field = 224.2e-4 224.2e-4 224.2e-4 ! Kick strength (Tesla) for 989 kick 60% (292) - kicker_params%kicker_field = 224.2e-4 224.2e-4 224.2e-4 ! kicker_params%kicker_field = 204.0e-4 204.0e-4 204.0e-4 ! 3.6mm with kick time 240ns kicker_params%kicker_field = 184.0e-4 184.0e-4 184.0e-4 ! 5.8mm with kick time 240ns kicker_params%kicker_field = 194.0e-4 194.0e-4 194.0e-4 ! 5.5mm with kick time 240ns and eta=7.97 kicker_params%kicker_field = 224.0e-4 224.0e-4 224.0e-4 !peak of scan with 180ns ! kicker_params%kicker_field = 0. 0. 0. ! Kick strength (Tesla) for 989 kick 60% (292) - ! kicker_params%kicker_field = 200.e-4 170.e-4 185.e-4 ! Renee ! : 3 reals - Rise time (seconds) of the flat top of the 3 kicker pulses kicker_params%dtRise = 20.e-9 20.e-9 20.e-9 ! Kick rise (s) ! : 3 reals - Fall time (seconds) of the flat top of the 3 kicker pulses kicker_params%dtFall = 20.e-9 20.e-9 20.e-9 ! Kick fall (s) ! : 3 reals - width of flat top of 3 kicker pulses kicker_params%kick_width = 50.e-9 50.e-9 50.e-9 ! Kick width (seconds) -- doesn't apply to kickerCircuit=821 ! : 3 reals - kicker start times, if kicker_tStart = -1 then automatically calculate kicker start time kicker_tStart = 2.0e-9 2.0e-9 2.0e-9 ! good choice for 127 cm kicker ! kicker_tStart = 10.0e-9 10.0e-9 10.0e-9 ! good choice for 60 cm kicker kicker_tStart = 120.0e-9 120.0e-9 120.0e-9 ! good choice for 60 cm kicker kicker_tStart = 200.0e-9 200.0e-9 200.0e-9 ! optimized for tdistr = "e821 kicker pulse" and kick 254.2G kicker_tStart = 240.0e-9 240.0e-9 240.0e-9 ! optimized for tdistr and 194G kick, "e821 kicker pulse" kicker_tStart = 40.0e-9 40.0e-9 40.0e-9 ! best for tdistr = "e989-P1" kicker_tStart = 70.0e-9 70.0e-9 70.0e-9 ! best for tdistr = "e989-P1", reoptimized energy_vs_time_vs_turn/config_025/20221112_081500 kicker_tStart = 120.0e-9 120.0e-9 120.0e-9 ! best for tdistr = "e989-P1", optimized for kicker_pulse_run3b kicker_tStart = 180.0e-9 180.0e-9 180.0e-9 ! best for tdistr = "e989-P1", centered by eye ! : logical - if true then include scattering in Q1 quad plates quad_plate = F ! include scatter in Q1 quad plates ! : integer - 821 or 989 (displaced Q1 outer plate) -- controls shape of E-field for Q1_Long, 1 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 p !: integer, 1=analytic, 2=mapped quadFieldType = 1 !: integer - Determines parameter to be varied. ! 0-nothing, 1-betax, 2-betay, 3-eta, 4-kick, 5-kick_width, 6-initial position offset, 7-initia ! 8-kicker rise time, 9-kicker fall time, 10-inflector field (1.0=-B_magic), 11-initial energy o ! 14-kicker_tStart vparam_id = 4 ! 0-nothing, 1-betax, 2-betay, 3-eta, 4-kick, 5-kick_width, 6-initial position offset, 7-initial ! : real - minimum value of parameter to be varied vparam_min = 0.6 ! : real - maximum value of parameter to be varied vparam_max =1.5 ! maximum of parameter to vary !: real - change in value of parameter in each step delta_vparam = 0.1 !20.e-9 ! parameter step !: integer - peak to peak CBO amplitude is determined over turns nturns - vparam_cbo_turns . default = vparam_cbo_turns = 10 !: logical - if true vary parameter (vparam_id) accoring to - vparam_min,vparam_max,delta_vparam. If false then d loop = F !: field_file_struct ! character_string - name of file with magnetic field map of fringe field in injection channel. ! -4 real - grid_spacing ! integer - type, type of read format = 1,2 (3 D map), 3 (extend 2D map to 3 dimensions) or 4 (uniform fie ! integer - flag (miscellaneous) ! possibilities include - 'injec_fld.dat' 0.5 0. 1 0, or 'ring_inject_bfield3d_wuzheng_extended.dat' 0.5 0. 2 0 fringe_file = 'ring_inject_bfield3d_wuzheng_extended.dat' 0.5 0. 2 0 !field map for fringe field, !: field_file_struct ! character_string - name of file with magnetic field map of inflector field. ! real - grid_spacing ! integer - type, type of read format = 1,2 (3 D map), 3 (extend 2D map to 3 dimensions) or 4 (uniform fie ! integer - flag (miscellaneous) ! possibilities include - 'inf_field_alone.dat' 0.5 0. 1 0, 'uniform' 0.5 0. 4 0 inflector_file = 'inf_field_alone.dat' 0.5 0. 2 0 !field map for inflector field, grid spacing, gradient, !: logical - if true then include spin tracking of muons. Default = true spin_tracking_on = F !: logical - if true then include muon decay. Default = true muon_decay_on = F !: logical - if true save trajectory and spin every 5cm and layout of quads and kicker. Create png files for make_movie = F !: logical - if true then write gnu_input.dat and stop. Default = false !: logical - if true use twiss parameters defined in lattice file to propagate through injection line use_lattice_twiss = T !: logical - if true write phase space coordinates of each muon to a file at each element on first tu write_phase_space_file = T ! - if true allocate array muon_ele(:,:) to save distribution at each element. To reduce m save_element_by_element_info = F !: real - if non-zero, save moments vs time at end of each turn in bins of width = time_bin time_bin_width = 10.e-9 scraping_on = T !init_quad_focus = -0.145315 ! if = -0.145315 and if init_quad_steer = -0.145315 then same as 32kV, 22.7 kV config !init_quad_steer = -0.145315 init_quad_focus = -0.18421 ! if = -0.18421 and if init_quad_steer = -0.18421 then same as 19kV, 12 kV config if index= init_quad_steer = -0.18421 init_quad_focus = -0.14634 ! if = -0.14634 and if init_quad_steer = -0.14634 then same as 20.5kV, 14.5 kV config if in init_quad_steer = -0.14634 init_quad_focus = -0.14208 ! if = -0.14208 and if init_quad_steer = -0.14208 then same as 18.3kV, 13.1 kV config (18 init_quad_steer = -0.14208 quad_ramp_start_time = 15.e-6 quad_ramp_end_time = 35.e-6 error_fields = F azimuthal_exp_z ='BzFourier20170721.dat' !'BzFourier20170628_LogID983.dat' azimuthal_exp_r= 'BrFourier2016.dat' FiberScatter0 = F int_betaCrossE = T betaCrossE_verbose = F no_energy_change = F initial_offset_from_co = F quad_fringe_energy_change=F input_call = 'quad_input.dat' BetaCrossE_start_time = 35.e-6 write_phase_space_every_n_turns = 28 !4 microseconds !write_phase_space_every_n_turns = 1 !0.149 microseconds /