Add statistics function to Pythia interface.
Create opposite dir/equal energy particle gun setup.
Add argument parsing and refactor main.cc