def METHOD_NAME(options, outFile=None, errorFile=None):
  """The main track for handling PDB and CIF input formats, which will involve
     parsing the model hierarchy to get the dihedral angles for ourselves"""
  setOptions(options)
  import suiteninput  # must be AFTER setOptions
  if not outFile: outFile = sys.stdout
  if not errorFile: errorFile = sys.stderr
  inFile = options.infile
  model = loadModel(inFile)
  residues = getResidueDihedrals(model, options.altid,
                                 name=os.path.splitext(inFile)[0],
                                 errorFile=errorFile)
  ### to print mp_geo-like output:
  # for r in residues:
  #   print(residueString(r))
  # useful for seeing what suites were generated
  if residues is not None and len(residues) > 0:
    suiteList = suiteninput.buildSuites(residues)
    suiteList = suiteList[:-1]
    suiteList = compute(suiteList)
    finalStats()
    write(outFile, suiteList)
    clearStats()