Example #1
0
void ReactionCmlLoader::loadReaction (Reaction &rxn)
{
   rxn.clear();

   QS_DEF(Array<char>, buf);
   _scanner.readAll(buf);
   buf.push(0);

   TiXmlDocument xml;

   xml.Parse(buf.ptr());

   if (xml.Error())
      throw Error("XML parsing error: %s", xml.ErrorDesc());

   TiXmlHandle hxml(&xml);
   TiXmlElement *elem;
   TiXmlHandle hroot(0);

   elem = hxml.FirstChild("reaction").Element();
   if (elem == 0)
      elem = hxml.FirstChild("cml").FirstChild("reaction").Element();
   if (elem == 0)
      throw Error("no <reaction>?");
   hroot = TiXmlHandle(elem);

   const char *title = elem->Attribute("title");

   if (title != 0)
      rxn.name.readString(title, true);

   QS_DEF(Molecule, mol);

   elem = hroot.FirstChild("reactantList").FirstChild().Element();
   for (; elem; elem = elem->NextSiblingElement())
   {
      if (strcasecmp(elem->Value(), "molecule") != 0)
         continue;
      TiXmlHandle handle(elem);
      MoleculeCmlLoader loader(handle);
      loader.stereochemistry_options = stereochemistry_options;
      loader.loadMolecule(mol);
      rxn.addReactantCopy(mol, 0, 0);
   }

   elem = hroot.FirstChild("productList").FirstChild().Element();
   for (; elem; elem = elem->NextSiblingElement())
   {
      if (strcasecmp(elem->Value(), "molecule") != 0)
         continue;
      TiXmlHandle handle(elem);
      MoleculeCmlLoader loader(handle);
      loader.stereochemistry_options = stereochemistry_options;
      loader.loadMolecule(mol);
      rxn.addProductCopy(mol, 0, 0);
   }

   elem = hroot.FirstChild("spectatorList").FirstChild().Element();
   for (; elem; elem = elem->NextSiblingElement())
   {
      if (strcasecmp(elem->Value(), "molecule") != 0)
         continue;
      TiXmlHandle handle(elem);
      MoleculeCmlLoader loader(handle);
      loader.stereochemistry_options = stereochemistry_options;
      loader.loadMolecule(mol);
      rxn.addCatalystCopy(mol, 0, 0);
   }
}