Merge all options/reconstruction.*py and use multiple PodioOutput writers with branch filters

That will result in more consistent treatments.