Read events from testIOGenEvent.input Select events containing a photon of pT > 25 GeV Add arbitrary PDF information to one of the good events Write the selected events and read them back in using an istream