code: def make_phase_space_list():
      
      with open('dynac.short') as f:
          data_str = ''.join(line for line in f.readlines())
          data_str_array = data_str.split('beam (emit card)')[1:]
          data_str_matrix = [[j.strip().split() for j in i] for i in[chunk.split('\n')[1:8] for chunk in data_str_array]]
  
          return [PhaseSpace(data) for data in data_str_matrix]
docs: Extract all the phase space information (due to ``EMIT`` commands in the input
      file), and create a list of PhaseSpace objects.  The primary purpose of this
      is for interactive explorations of the data produced during Pynac simulations.