/**
* This method calculates the ionization potential of an atom.
*
* @param atom The IAtom to ionize.
* @param container Parameter is the IAtomContainer.
* @return The ionization potential. Not possible the ionization.
*/
@Override
public DescriptorValue calculate(IAtom atom, IAtomContainer container) {
double value = 0;
// FIXME: for now I'll cache a few modified atomic properties, and restore them at the end of
// this method
String originalAtomtypeName = atom.getAtomTypeName();
Integer originalNeighborCount = atom.getFormalNeighbourCount();
Integer originalValency = atom.getValency();
IAtomType.Hybridization originalHybrid = atom.getHybridization();
Double originalBondOrderSum = atom.getBondOrderSum();
Order originalMaxBondOrder = atom.getMaxBondOrder();
if (!isCachedAtomContainer(container)) {
try {
AtomContainerManipulator.percieveAtomTypesAndConfigureAtoms(container);
LonePairElectronChecker lpcheck = new LonePairElectronChecker();
lpcheck.saturate(container);
} catch (CDKException e) {
return new DescriptorValue(
getSpecification(),
getParameterNames(),
getParameters(),
new DoubleResult(Double.NaN),
getDescriptorNames(),
e);
}
}
try {
value = IonizationPotentialTool.predictIP(container, atom);
} catch (CDKException e) {
return new DescriptorValue(
getSpecification(),
getParameterNames(),
getParameters(),
new DoubleResult(Double.NaN),
getDescriptorNames(),
e);
}
// restore original props
atom.setAtomTypeName(originalAtomtypeName);
atom.setFormalNeighbourCount(originalNeighborCount);
atom.setValency(originalValency);
atom.setHybridization(originalHybrid);
atom.setMaxBondOrder(originalMaxBondOrder);
atom.setBondOrderSum(originalBondOrderSum);
return new DescriptorValue(
getSpecification(),
getParameterNames(),
getParameters(),
new DoubleResult(value),
getDescriptorNames());
}
@Test
public void testConfigure_IAtom() throws Exception {
IAtomType atomType;
IAtom atom = new org.openscience.cdk.Atom("X");
atom.setAtomTypeName("C.ar");
AtomTypeFactory factory =
AtomTypeFactory.getInstance(
"org/openscience/cdk/config/data/mol2_atomtypes.xml", new ChemObject().getBuilder());
atomType = factory.configure(atom);
Assert.assertNotNull(atomType);
Assert.assertEquals("C", atom.getSymbol());
}