/* $RCSfile$ * $Author$ * $Date$ * $Revision$ * * Copyright (C) 2003-2007 The Chemistry Development Kit (CDK) project * * Contact: cdk-devel@lists.sourceforge.net * * This program is free software; you can redistribute it and/or * modify it under the terms of the GNU Lesser General Public License * as published by the Free Software Foundation; either version 2.1 * of the License, or (at your option) any later version. * All we ask is that proper credit is given for our work, which includes * - but is not limited to - adding the above copyright notice to the beginning * of your source code files, and to any copyright notice that you may distribute * with programs based on this work. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. * */ package org.openscience.cdk.tools.manipulator; import java.util.ArrayList; import java.util.List; import org.openscience.cdk.annotations.TestClass; import org.openscience.cdk.annotations.TestMethod; import org.openscience.cdk.interfaces.IAtom; import org.openscience.cdk.interfaces.IAtomContainer; import org.openscience.cdk.interfaces.IBond; import org.openscience.cdk.interfaces.IChemObject; import org.openscience.cdk.interfaces.IElectronContainer; import org.openscience.cdk.interfaces.IMolecule; import org.openscience.cdk.interfaces.IMoleculeSet; import org.openscience.cdk.interfaces.IReaction; import org.openscience.cdk.interfaces.IReactionSet; /** * @cdk.module standard * @cdk.githash * * @see ChemModelManipulator */ @TestClass("org.openscience.cdk.tools.manipulator.ReactionSetManipulatorTest") public class ReactionSetManipulator { @TestMethod("testGetAtomCount_IReactionSet") public static int getAtomCount(IReactionSet set) { int count = 0; for (IReaction iReaction : set.reactions()) { count += ReactionManipulator.getAtomCount(iReaction); } return count; } @TestMethod("testGetBondCount_IReactionSet") public static int getBondCount(IReactionSet set) { int count = 0; for (IReaction iReaction : set.reactions()) { count += ReactionManipulator.getBondCount(iReaction); } return count; } @TestMethod("testRemoveAtomAndConnectedElectronContainers_IReactionSet_IAtom") public static void removeAtomAndConnectedElectronContainers(IReactionSet set, IAtom atom) { for (IReaction reaction : set.reactions()) { ReactionManipulator.removeAtomAndConnectedElectronContainers(reaction, atom); } } @TestMethod("testRemoveElectronContainer_IReactionSet_IElectronContainer") public static void removeElectronContainer(IReactionSet set, IElectronContainer electrons) { for (IReaction reaction : set.reactions()) { ReactionManipulator.removeElectronContainer(reaction, electrons); } } /** * get all Molecules object from a set of Reactions. * * @param set The set of reaction to inspect * @return The IMoleculeSet */ @TestMethod("testGetAllMolecules_IReactionSet") public static IMoleculeSet getAllMolecules(IReactionSet set) { IMoleculeSet moleculeSet = set.getBuilder().newMoleculeSet(); for (IReaction reaction : set.reactions()) { IMoleculeSet molecules = ReactionManipulator.getAllMolecules(reaction); for (IAtomContainer ac : molecules.molecules()) { boolean contain = false; for (IAtomContainer atomContainer : moleculeSet.molecules()) { if (atomContainer.equals(ac)) { contain = true; break; } } if (!contain) moleculeSet.addMolecule((IMolecule) (ac)); } } return moleculeSet; } @TestMethod("testGetAllIDs_IReactionSet") public static List<String> getAllIDs(IReactionSet set) { List<String> IDlist = new ArrayList<String>(); if (set.getID() != null) IDlist.add(set.getID()); for (IReaction reaction : set.reactions()) { IDlist.addAll(ReactionManipulator.getAllIDs(reaction)); } return IDlist; } /** * Returns all the AtomContainer's of a Reaction. * @param set the reaction set to get the molecules from * @return a List containing the IAtomContainer objects in the IReactionSet */ @TestMethod("testGetAllAtomContainers_IReactionSet") public static List<IAtomContainer> getAllAtomContainers(IReactionSet set) { return MoleculeSetManipulator.getAllAtomContainers( getAllMolecules(set) ); } @TestMethod("testGetRelevantReaction_IReactionSet_IAtom") public static IReaction getRelevantReaction(IReactionSet set, IAtom atom) { for (IReaction reaction : set.reactions()) { IAtomContainer container = ReactionManipulator.getRelevantAtomContainer(reaction, atom); if (container != null) { // a match! return reaction; } } return null; } @TestMethod("testGetRelevantReaction_IReactionSet_IBond") public static IReaction getRelevantReaction(IReactionSet set, IBond bond) { for (IReaction reaction : set.reactions()) { IAtomContainer container = ReactionManipulator.getRelevantAtomContainer(reaction, bond); if (container != null) { // a match! return reaction; } } return null; } /** * Get all Reactions object containing a Molecule from a set of Reactions. * * @param reactSet The set of reaction to inspect * @param molecule The molecule to find * @return The IReactionSet */ @TestMethod("testGetRelevantReactions_IReactionSet_IMolecule") public static IReactionSet getRelevantReactions(IReactionSet reactSet, IAtomContainer molecule) { IReactionSet newReactSet = reactSet.getBuilder().newReactionSet(); IReactionSet reactSetProd = getRelevantReactionsAsProduct(reactSet, molecule); for (IReaction reaction : reactSetProd.reactions()) newReactSet.addReaction(reaction); IReactionSet reactSetReact = getRelevantReactionsAsReactant(reactSet, molecule); for (IReaction reaction : reactSetReact.reactions()) newReactSet.addReaction(reaction); return newReactSet; } /** * Get all Reactions object containing a Molecule as a Reactant from a set * of Reactions. * * @param reactSet The set of reaction to inspect * @param molecule The molecule to find as a reactant * @return The IReactionSet */ @TestMethod("testGetRelevantReactionsAsReactant_IReactionSet_IMolecule") public static IReactionSet getRelevantReactionsAsReactant( IReactionSet reactSet, IAtomContainer molecule) { IReactionSet newReactSet = reactSet.getBuilder().newReactionSet(); for (IReaction reaction : reactSet.reactions()) { for(IAtomContainer atomContainer : reaction.getReactants().molecules()) if(atomContainer.equals(molecule)) newReactSet.addReaction(reaction); } return newReactSet; } /** * Get all Reactions object containing a Molecule as a Product from a set of * Reactions. * * @param reactSet The set of reaction to inspect * @param molecule The molecule to find as a product * @return The IReactionSet */ @TestMethod("testGetRelevantReactionsAsProduct_IReactionSet_IMolecule") public static IReactionSet getRelevantReactionsAsProduct( IReactionSet reactSet, IAtomContainer molecule) { IReactionSet newReactSet = reactSet.getBuilder().newReactionSet(); for (IReaction reaction : reactSet.reactions()) { for(IAtomContainer atomContainer : reaction.getProducts().molecules()) if(atomContainer.equals(molecule)) newReactSet.addReaction(reaction); } return newReactSet; } @TestMethod("testGetRelevantAtomContainer_IReactionSet_IAtom") public static IAtomContainer getRelevantAtomContainer(IReactionSet set, IAtom atom) { for (IReaction reaction : set.reactions()) { IAtomContainer container = ReactionManipulator.getRelevantAtomContainer(reaction, atom); if (container != null) { // a match! return container; } } return null; } @TestMethod("testGetRelevantAtomContainer_IReactionSet_IBond") public static IAtomContainer getRelevantAtomContainer(IReactionSet set, IBond bond) { for (IReaction reaction : set.reactions()) { IAtomContainer container = ReactionManipulator.getRelevantAtomContainer(reaction, bond); if (container != null) { // a match! return container; } } return null; } @TestMethod("testSetAtomProperties_IReactionSet_Object_Object") public static void setAtomProperties(IReactionSet set, Object propKey, Object propVal) { for (IReaction reaction : set.reactions()) { ReactionManipulator.setAtomProperties(reaction, propKey, propVal); } } @TestMethod("testGetAllChemObjects_IReactionSet") public static List<IChemObject> getAllChemObjects(IReactionSet set) { ArrayList<IChemObject> list = new ArrayList<IChemObject>(); list.add(set); for (IReaction reaction : set.reactions()) { list.addAll(ReactionManipulator.getAllChemObjects(reaction)); } return list; } /** * Gets a reaction from a ReactionSet by ID of any product or reactant. If several exist, * only the first one will be returned. * * @param reactionSet The reactionSet to search in * @param id The id to search for. * @return The Reaction or null; */ @TestMethod("testGetReactionByAtomContainerID_IReactionSet_String") public static IReaction getReactionByAtomContainerID(IReactionSet reactionSet, String id) { for (IReaction reaction : reactionSet.reactions()) { if(AtomContainerSetManipulator.containsByID(reaction.getProducts(),id)) return reaction; } for (IReaction reaction : reactionSet.reactions()) { if(AtomContainerSetManipulator.containsByID(reaction.getReactants(),id)) return reaction; } return null; } /** * Gets a reaction from a ReactionSet by ID. If several exist, * only the first one will be returned. * * @param reactionSet The reactionSet to search in * @param id The id to search for. * @return The Reaction or null; */ @TestMethod("testGetReactionByReactionID_IReactionSet_String") public static IReaction getReactionByReactionID(IReactionSet reactionSet, String id) { Iterable<IReaction> reactionIter = reactionSet.reactions(); for (IReaction reaction : reactionIter) { if (reaction.getID()!=null && reaction.getID().equals(id)) { return reaction; } } return null; } }