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.

def make_phase_space_list():
    """
    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.
    """
    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]