Would be useful to include functionality so that the code can calculate half and and third energy spectra in addition to the full component, a given beam file. Would require a slight rewrite of "read_beam.pro" as it currently only uses the first index (ie. full beam component) and adding an option for the "batch.xml" run file