srv_parser.cpp

//
// =================================================================
//
//    28.10.02   <--  Date of Last Modification.
//                   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//  -------------------------------------------------------------------
//
//  **** Module  :  SRV_Parser <implementation>
//       ~~~~~~~~~
//  **** Project :  SSM Server
//       ~~~~~~~~~
//  **** Functions : main  ( example parsing SSM XML output )
//       ~~~~~~~~~~~
//
//  E. Krissinel, 2002
//
// =================================================================
//


#ifndef  __SRV_PXML__
#include "srv_pxml.h"
#endif


void  PrintInstructions ( pstr argv0 )  {
  printf ( 
    "\n"
    "SSM XML parser example\n"
    "~~~~~~~~~~~~~~~~~~~~~~\n"
    "\n"
    "%s  ssmout.xml\n"
    "\n"
    "where 'ssmout.xml' is an output XML file.\n"
    "\n",argv0
         );
}


int main ( int argc, char ** argv, char ** env )  {
PCSSMXML SSMXML;
char     queryID[500],targetID[500];
char     queryCh[10] ,targetCh[10];
mat44    RT;
int      rc, nMatches,i;

  if (argc<2)  {
    PrintInstructions ( argv[0] );
    return 1;
  }

  //  1.  Make routine initializations, which must always be done
  //      before working with MMDB
  InitMatType();

  //  2.  Create an instance of XML object
  SSMXML = new CSSMXML();

  //  3.  Read XML file.
  rc = SSMXML->readSSMXML ( argv[1] );
  printf ( "  SSM XML file read rc=%i\n",rc );

  if (rc)  return rc;

  nMatches = SSMXML->getNofMatches();
  printf ( "  Total %i matches\n",nMatches );
  if (nMatches<=0)  return 2;

  rc = SSMXML->getQueryID ( queryID );
  if (rc==query_pdbentry)  printf ( "  Query: PDB code %s\n",queryID );
  else if (rc==query_file) printf ( "  Query: coordinate file %s\n",queryID );
                      else printf ( "  Error retrieving query ID #%i\n",rc );

  for (i=1;i<=nMatches;i++)  {

    printf ( "  MATCH #%i\n",i );

    rc = SSMXML->getTargetID ( targetID,i );
    switch (rc)  {
      case target_pdbentry   :
        printf ( "  Target: PDB code %s\n" ,targetID );   break;
      case target_file       :
        printf ( "  Target: file %s\n"     ,targetID );   break;
      case target_pdb        :
        printf ( "  Target: PDB entry %s\n",targetID );   break;
      case target_scop       :
        printf ( "  Target: SCOP entry %s\n",targetID );  break;
      case target_scopsubset :
        printf ( "  Target: SCOP subset entry %s\n",targetID );  break;
      default :
        printf ( "  Error retrueving target ID #%i\n",rc );
    }

    SSMXML->getQueryChain  ( queryCh ,i );
    SSMXML->getTargetChain ( targetCh,i );
    printf ( "          query  chain '%s'\n"
             "          target chain '%s'\n", queryCh,targetCh );

    printf ( "          P-score %11.4f, Z-score %11.4f, RPDBe %11.4f\n",
             SSMXML->getPScore(i),SSMXML->getZScore(i),SSMXML->getRMSD(i) );

    rc = SSMXML->getRTMatrix ( RT,i );
    if (rc)  {
      printf ( "          Error #%i at retrieving RT-matrix.\n",rc );
      break;
    }
    printf ( "          RT-matrix: %10.3f %10.3f %10.3f   %10.3f\n"
             "                     %10.3f %10.3f %10.3f   %10.3f\n"
             "                     %10.3f %10.3f %10.3f   %10.3f\n",
             RT[0][0], RT[0][1], RT[0][2], RT[0][3],
             RT[1][0], RT[1][1], RT[1][2], RT[1][3],
             RT[2][0], RT[2][1], RT[2][2], RT[2][3] );

  }

  return 0;

}

---