Improve trajectory loading with replay ensemble

[OndrejMarsalek]

The replay ensemble should be able to take a set of trajectory files, one for each bead, rather than a single file with all the beads. This way, normal output from a run can be used.

We should also document better what the expected input is and add one example that uses it.