from __future__ import division
from libtbx.utils import Sorry, Usage
import libtbx.phil.command_line
import sys

master_phil = libtbx.phil.parse("""
resname = None
  .type = str
d_max = None
  .type = float
polymeric_type = *Any Free Polymeric
  .type = choice
xray_only = True
  .type = bool
data_only = False
  .type = bool
identity_cutoff = None
  .type = int
quiet = False
  .type = bool
""")

def run (args, out=sys.stdout) :
  if (len(args) == 0) or ("--help" in args) :
    raise Usage("""mmtbx.find_residue_in_pdb RESNAME [options]

Use the RCSB web services to retrieve a list of PDB structures containing the
specified chemical ID.

Full parameters:
%s
""" % master_phil.as_str(prefix="  "))
  sources = []
  def process_unknown (arg) :
    if (1 <= len(arg) <= 3) and (arg.isalnum()) :
      return libtbx.phil.parse("resname=%s" % arg)
  cai = libtbx.phil.command_line.argument_interpreter(master_phil=master_phil)
  working_phil = cai.process_and_fetch(args=args,
    custom_processor=process_unknown)
  params = working_phil.extract()
  if (params.resname is None) :
    raise Sorry("No residue ID specified.")
  from mmtbx.wwpdb import rcsb_web_services
  pdb_ids = rcsb_web_services.chemical_id_search(
    resname=params.resname,
    d_max=params.d_max,
    polymeric_type=params.polymeric_type,
    xray_only=params.xray_only,
    data_only=params.data_only,
    identity_cutoff=params.identity_cutoff)
  pdb_ids = [ id.lower() for id in pdb_ids ]
  if (len(pdb_ids) == 0) :
    raise Sorry("No structures found matching the specified criteria.")
  else :
    if (not params.quiet) :
      print >> out, "%d PDB IDs retrieved:" % len(pdb_ids)
      i = 0
      while (i < len(pdb_ids)) :
        print >> out, "  %s" % " ".join(pdb_ids[i:i+16])
        i += 16
    else :
      print >> out, "%d PDB IDs matching" % len(pdb_ids)

if (__name__ == "__main__") :
  run(sys.argv[1:])