g2_tracking example input files                  February 5, 2016

Example4
   g-2/examples/example4 - Start of tracking for a single particle is at the hole in the iron.
   Muon decay and spin tracking are off
   Muons are tracked until they hit an aperture (mostly defined in lattice file, but occasionally hard coded) 
   
The lattice file is bmad_short_kickers. Calorimeters can be included but are commented out in the lattice file.

! <lat_file_name> :character string -  BMAD format lattice file
  lat_file_name = 'bmad_short_kickers.'

! <nturns> :integer - Number of turns to track
  nturns   = 20           

! the number of muons is set to 1
! <nmuons> :integer - Number of muons to track
  nmuons = 1
! if the nmuons = 1, then there is no distribution, and parameters defining the distribution are ignored

The muon distribution is read from a file
! <muon_file> : If you choose to read an existing distribution, here is where the name is written. 
  muon_file = ' '       ! read the phase space distribution off the target from this file. If blank, use freshly created distribution

! since there is only a single muon, twiss parameters are irrelevant
! <twiss> :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.  


Tracking begins at the start of the injection channel (hole in the magnet iron) or in the ring the end of the inflector
Here tracking it will start at the hole. 
!<start_tracking_at_inflector_exit>: logical - if true then start tracking the distribution at the inflector exit, if false, start tracking at end of M5 line
  start_tracking_at_inflector_exit = F

Scattering and energy loss in the inflector end coils is turned off
! <inf_end_us> :logical - if true include scattering at  upstream  end of inflector (E821 coils)
  inf_end_us   = F                                    
! <inf_end_ds> :logical - if true include scattering at  downstream  end of inflector (E821 coils)
  inf_end_ds   = F                                     
! <enerloss> :logical - if true include energy loss with scattering in inflector(E821 coils)!
  enerloss = F

The coordinates of the single muon can be offset in each of the phase space coordinates
Since we are starting tracking through the injection channel, we need to choose 
just the right horizontal angle and offset for the trajectory to emerge from the center of the inflector
! <initial_offsets> :real - x, y, z, x-angle, y-angle, deltaE/E 
!   - offset of beam at entrance into the hole in the ring iron, with respect to the tangent through the center of the inflector 
    initial_offsets = -4.7671E-02  0. 0. 2.2134E-02 0. 0.  !these are the offsets you need if starting at the hole in the iron

! kicker fields 
!<kicker_params%kicker_field> : 3 reals - Bfield (T) for the three kicker mangets
    kicker_params%kicker_field      =   326.6e-4   304.8e-4  240.2e-4 ! Kick strength (Tesla) for 989 kick 1.27m kickers

Scattering in the quad plates is turned off
!<quad_plate> : logical - if true then include scattering in Q1 quad plates
  quad_plate    = F                                    ! include scatter in Q1 quad plates

The strengths of the three kickers are chosen to steer the muons onto the magic radius.
!<kicker_params%kicker_field> : 3 reals - Bfield (T) for the three kicker mangets
    kicker_params%kicker_field      =   326.6e-4   304.8e-4  240.2e-4 ! Kick strength (Tesla) for 989 kick 1.27m kickers

Muon spin tracking and decay are off

!<spin_tracking_on>: logical - if true then include spin tracking of muons. Default = true
 spin_tracking_on = F
!<muon_decay>: logical - if true then include muon decay. Default = true
 muon_decay_on = F

Some output is controlled by the following logicals.
!<write_phase_space>: logical - if true write phase space coordinates of each muon to a file at each element on first turn 
write_phase_space_file = F

The average and rms of each of the phase space coordinates computed on each turn is written to 'Beam_moments_tbt.dat'
If you want to see how the distribution changes going around the ring, set <save_element_by_element_info = T>
Then the average and rms, computed at the end of each element, are written to 'Beam_moments_by_element.dat' 
!<save_element_by_element_info> - if true allocate array muon_ele(:,:) to save distribution at each element. To reduce memory requirements - false
save_element_by_element_info = T


For all the details see the input.dat file

 
Files:
  input.dat
  bmad_short_kickers_trackers.

Also
  /g-2/files/pulse_127cm.dat
  /g-2/magneticfield/ring_inject_bfield3d_wuzheng_extended.dat 
  /g-2/magneticfield/inf_field_alone.dat 
  /g-2/files/KICKER_E989_NoRolledEdges.dat