Split up long function to select initial state partons
Also rename it (filter_partons -> allowed_incoming)