• 'trajectory_ele_by_ele.dat'
  nothing there
• 'trajectory.dat'
  If single particle,
• 'beam_trajectory.dat' br> N particles > 1,
• 'phasespace.dat'
  if(phase_space =T) phase space of all particles on each turn, <# turn = i, and then particle(:)%vec
• 'extracted_beam.dat'
  Nothing there
• 'target_beam.dat' phase space for particles that reach target on each turn