docs: def filter_muons(blob):
      """Write all muons from McTracks to Muons."""
      tracks = blob['McTracks']
      muons = tracks[tracks.type == -13]    # PDG particle code
      blob["Muons"] = Table(muons)
      return blob