diff mayachemtools/docs/scripts/html/code/TopologicalAtomPairsFingerprints.html @ 0:73ae111cf86f draft

Uploaded
author deepakjadmin
date Wed, 20 Jan 2016 11:55:01 -0500
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mayachemtools/docs/scripts/html/code/TopologicalAtomPairsFingerprints.html	Wed Jan 20 11:55:01 2016 -0500
@@ -0,0 +1,764 @@
+<html>
+<head>
+<title>MayaChemTools:Code:TopologicalAtomPairsFingerprints.pl</title>
+<meta http-equiv="content-type" content="text/html;charset=utf-8">
+<link rel="stylesheet" type="text/css" href="../../../css/MayaChemToolsCode.css">
+</head>
+<body leftmargin="20" rightmargin="20" topmargin="10" bottommargin="10">
+<br/>
+<center>
+<a href="http://www.mayachemtools.org" title="MayaChemTools Home"><img src="../../../images/MayaChemToolsLogo.gif" border="0" alt="MayaChemTools"></a>
+</center>
+<br/>
+<pre>
+   1 #!/usr/bin/perl -w
+   2 <span class="c">#</span>
+   3 <span class="c"># $RCSfile: TopologicalAtomPairsFingerprints.pl,v $</span>
+   4 <span class="c"># $Date: 2015/02/28 20:46:22 $</span>
+   5 <span class="c"># $Revision: 1.34 $</span>
+   6 <span class="c">#</span>
+   7 <span class="c"># Author: Manish Sud &lt;msud@san.rr.com&gt;</span>
+   8 <span class="c">#</span>
+   9 <span class="c"># Copyright (C) 2015 Manish Sud. All rights reserved.</span>
+  10 <span class="c">#</span>
+  11 <span class="c"># This file is part of MayaChemTools.</span>
+  12 <span class="c">#</span>
+  13 <span class="c"># MayaChemTools is free software; you can redistribute it and/or modify it under</span>
+  14 <span class="c"># the terms of the GNU Lesser General Public License as published by the Free</span>
+  15 <span class="c"># Software Foundation; either version 3 of the License, or (at your option) any</span>
+  16 <span class="c"># later version.</span>
+  17 <span class="c">#</span>
+  18 <span class="c"># MayaChemTools is distributed in the hope that it will be useful, but without</span>
+  19 <span class="c"># any warranty; without even the implied warranty of merchantability of fitness</span>
+  20 <span class="c"># for a particular purpose.  See the GNU Lesser General Public License for more</span>
+  21 <span class="c"># details.</span>
+  22 <span class="c">#</span>
+  23 <span class="c"># You should have received a copy of the GNU Lesser General Public License</span>
+  24 <span class="c"># along with MayaChemTools; if not, see &lt;http://www.gnu.org/licenses/&gt; or</span>
+  25 <span class="c"># write to the Free Software Foundation Inc., 59 Temple Place, Suite 330,</span>
+  26 <span class="c"># Boston, MA, 02111-1307, USA.</span>
+  27 <span class="c">#</span>
+  28 
+  29 <span class="k">use</span> <span class="w">strict</span><span class="sc">;</span>
+  30 <span class="k">use</span> <span class="w">FindBin</span><span class="sc">;</span> <span class="k">use</span> <span class="w">lib</span> <span class="q">&quot;$FindBin::Bin/../lib&quot;</span><span class="sc">;</span>
+  31 <span class="k">use</span> <span class="w">Getopt::Long</span><span class="sc">;</span>
+  32 <span class="k">use</span> <span class="w">File::Basename</span><span class="sc">;</span>
+  33 <span class="k">use</span> <span class="w">Text::ParseWords</span><span class="sc">;</span>
+  34 <span class="k">use</span> <span class="w">Benchmark</span><span class="sc">;</span>
+  35 <span class="k">use</span> <span class="w">FileUtil</span><span class="sc">;</span>
+  36 <span class="k">use</span> <span class="w">TextUtil</span><span class="sc">;</span>
+  37 <span class="k">use</span> <span class="w">SDFileUtil</span><span class="sc">;</span>
+  38 <span class="k">use</span> <span class="w">MoleculeFileIO</span><span class="sc">;</span>
+  39 <span class="k">use</span> <span class="w">FileIO::FingerprintsSDFileIO</span><span class="sc">;</span>
+  40 <span class="k">use</span> <span class="w">FileIO::FingerprintsTextFileIO</span><span class="sc">;</span>
+  41 <span class="k">use</span> <span class="w">FileIO::FingerprintsFPFileIO</span><span class="sc">;</span>
+  42 <span class="k">use</span> <span class="w">AtomTypes::AtomicInvariantsAtomTypes</span><span class="sc">;</span>
+  43 <span class="k">use</span> <span class="w">AtomTypes::FunctionalClassAtomTypes</span><span class="sc">;</span>
+  44 <span class="k">use</span> <span class="w">Fingerprints::TopologicalAtomPairsFingerprints</span><span class="sc">;</span>
+  45 
+  46 <span class="k">my</span><span class="s">(</span><span class="i">$ScriptName</span><span class="cm">,</span> <span class="i">%Options</span><span class="cm">,</span> <span class="i">$StartTime</span><span class="cm">,</span> <span class="i">$EndTime</span><span class="cm">,</span> <span class="i">$TotalTime</span><span class="s">)</span><span class="sc">;</span>
+  47 
+  48 <span class="c"># Autoflush STDOUT</span>
+  49 <span class="i">$|</span> = <span class="n">1</span><span class="sc">;</span>
+  50 
+  51 <span class="c"># Starting message...</span>
+  52 <span class="i">$ScriptName</span> = <span class="i">basename</span><span class="s">(</span><span class="i">$0</span><span class="s">)</span><span class="sc">;</span>
+  53 <span class="k">print</span> <span class="q">&quot;\n$ScriptName: Starting...\n\n&quot;</span><span class="sc">;</span>
+  54 <span class="i">$StartTime</span> = <span class="w">new</span> <span class="w">Benchmark</span><span class="sc">;</span>
+  55 
+  56 <span class="c"># Get the options and setup script...</span>
+  57 <span class="i">SetupScriptUsage</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+  58 <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">help</span>} || <span class="i">@ARGV</span> &lt; <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
+  59   <span class="k">die</span> <span class="i">GetUsageFromPod</span><span class="s">(</span><span class="q">&quot;$FindBin::Bin/$ScriptName&quot;</span><span class="s">)</span><span class="sc">;</span>
+  60 <span class="s">}</span>
+  61 
+  62 <span class="k">my</span><span class="s">(</span><span class="i">@SDFilesList</span><span class="s">)</span><span class="sc">;</span>
+  63 <span class="i">@SDFilesList</span> = <span class="i">ExpandFileNames</span><span class="s">(</span>\<span class="i">@ARGV</span><span class="cm">,</span> <span class="q">&quot;sdf sd&quot;</span><span class="s">)</span><span class="sc">;</span>
+  64 
+  65 <span class="c"># Process options...</span>
+  66 <span class="k">print</span> <span class="q">&quot;Processing options...\n&quot;</span><span class="sc">;</span>
+  67 <span class="k">my</span><span class="s">(</span><span class="i">%OptionsInfo</span><span class="s">)</span><span class="sc">;</span>
+  68 <span class="i">ProcessOptions</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+  69 
+  70 <span class="c"># Setup information about input files...</span>
+  71 <span class="k">print</span> <span class="q">&quot;Checking input SD file(s)...\n&quot;</span><span class="sc">;</span>
+  72 <span class="k">my</span><span class="s">(</span><span class="i">%SDFilesInfo</span><span class="s">)</span><span class="sc">;</span>
+  73 <span class="i">RetrieveSDFilesInfo</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+  74 
+  75 <span class="c"># Process input files..</span>
+  76 <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span>
+  77 <span class="k">if</span> <span class="s">(</span><span class="i">@SDFilesList</span> &gt; <span class="n">1</span><span class="s">)</span> <span class="s">{</span>
+  78   <span class="k">print</span> <span class="q">&quot;\nProcessing SD files...\n&quot;</span><span class="sc">;</span>
+  79 <span class="s">}</span>
+  80 <span class="k">for</span> <span class="i">$FileIndex</span> <span class="s">(</span><span class="n">0</span> .. <span class="i">$#SDFilesList</span><span class="s">)</span> <span class="s">{</span>
+  81   <span class="k">if</span> <span class="s">(</span><span class="i">$SDFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$FileIndex</span>]<span class="s">)</span> <span class="s">{</span>
+  82     <span class="k">print</span> <span class="q">&quot;\nProcessing file $SDFilesList[$FileIndex]...\n&quot;</span><span class="sc">;</span>
+  83     <span class="i">GenerateTopologicalAtomPairsFingerprints</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span>
+  84   <span class="s">}</span>
+  85 <span class="s">}</span>
+  86 <span class="k">print</span> <span class="q">&quot;\n$ScriptName:Done...\n\n&quot;</span><span class="sc">;</span>
+  87 
+  88 <span class="i">$EndTime</span> = <span class="w">new</span> <span class="w">Benchmark</span><span class="sc">;</span>
+  89 <span class="i">$TotalTime</span> = <span class="w">timediff</span> <span class="s">(</span><span class="i">$EndTime</span><span class="cm">,</span> <span class="i">$StartTime</span><span class="s">)</span><span class="sc">;</span>
+  90 <span class="k">print</span> <span class="q">&quot;Total time: &quot;</span><span class="cm">,</span> <span class="i">timestr</span><span class="s">(</span><span class="i">$TotalTime</span><span class="s">)</span><span class="cm">,</span> <span class="q">&quot;\n&quot;</span><span class="sc">;</span>
+  91 
+  92 <span class="c">###############################################################################</span>
+  93 
+  94 <span class="c"># Generate fingerprints for a SD file...</span>
+  95 <span class="c">#</span>
+<a name="GenerateTopologicalAtomPairsFingerprints-"></a>  96 <span class="k">sub </span><span class="m">GenerateTopologicalAtomPairsFingerprints</span> <span class="s">{</span>
+  97   <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+  98   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$IgnoredCmpdCount</span><span class="cm">,</span> <span class="i">$SDFile</span><span class="cm">,</span> <span class="i">$MoleculeFileIO</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span><span class="sc">;</span>
+  99 
+ 100   <span class="i">$SDFile</span> = <span class="i">$SDFilesList</span>[<span class="i">$FileIndex</span>]<span class="sc">;</span>
+ 101 
+ 102   <span class="c"># Setup output files...</span>
+ 103   <span class="c">#</span>
+ 104   <span class="s">(</span><span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span> = <span class="i">SetupAndOpenOutputFiles</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span>
+ 105 
+ 106   <span class="i">$MoleculeFileIO</span> = <span class="w">new</span> <span class="i">MoleculeFileIO</span><span class="s">(</span><span class="q">&#39;Name&#39;</span> <span class="cm">=&gt;</span> <span class="i">$SDFile</span><span class="s">)</span><span class="sc">;</span>
+ 107   <span class="i">$MoleculeFileIO</span><span class="i">-&gt;Open</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 108 
+ 109   <span class="i">$CmpdCount</span> = <span class="n">0</span><span class="sc">;</span>
+ 110   <span class="i">$IgnoredCmpdCount</span> = <span class="n">0</span><span class="sc">;</span>
+ 111 
+ 112   <span class="j">COMPOUND:</span> <span class="k">while</span> <span class="s">(</span><span class="i">$Molecule</span> = <span class="i">$MoleculeFileIO</span><span class="i">-&gt;ReadMolecule</span><span class="s">(</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 113     <span class="i">$CmpdCount</span>++<span class="sc">;</span>
+ 114 
+ 115     <span class="c"># Filter compound data before calculating fingerprints...</span>
+ 116     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">Filter</span>}<span class="s">)</span> <span class="s">{</span>
+ 117       <span class="k">if</span> <span class="s">(</span><span class="i">CheckAndFilterCompound</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 118         <span class="i">$IgnoredCmpdCount</span>++<span class="sc">;</span>
+ 119         <span class="k">next</span> <span class="j">COMPOUND</span><span class="sc">;</span>
+ 120       <span class="s">}</span>
+ 121     <span class="s">}</span>
+ 122 
+ 123     <span class="i">$TopologicalAtomPairsFingerprints</span> = <span class="i">GenerateMoleculeFingerprints</span><span class="s">(</span><span class="i">$Molecule</span><span class="s">)</span><span class="sc">;</span>
+ 124     <span class="k">if</span> <span class="s">(</span>!<span class="i">$TopologicalAtomPairsFingerprints</span><span class="s">)</span> <span class="s">{</span>
+ 125       <span class="i">$IgnoredCmpdCount</span>++<span class="sc">;</span>
+ 126       <span class="i">ProcessIgnoredCompound</span><span class="s">(</span><span class="q">&#39;FingerprintsGenerationFailed&#39;</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span><span class="sc">;</span>
+ 127       <span class="k">next</span> <span class="j">COMPOUND</span><span class="sc">;</span>
+ 128     <span class="s">}</span>
+ 129 
+ 130     <span class="i">WriteDataToOutputFiles</span><span class="s">(</span><span class="i">$FileIndex</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span><span class="sc">;</span>
+ 131   <span class="s">}</span>
+ 132   <span class="i">$MoleculeFileIO</span><span class="i">-&gt;Close</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 133 
+ 134   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPSDFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 135     <span class="i">$NewFPSDFileIO</span><span class="i">-&gt;Close</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 136   <span class="s">}</span>
+ 137   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPTextFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 138     <span class="i">$NewFPTextFileIO</span><span class="i">-&gt;Close</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 139   <span class="s">}</span>
+ 140   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 141     <span class="i">$NewFPFileIO</span><span class="i">-&gt;Close</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 142   <span class="s">}</span>
+ 143 
+ 144   <span class="i">WriteFingerprintsGenerationSummaryStatistics</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$IgnoredCmpdCount</span><span class="s">)</span><span class="sc">;</span>
+ 145 <span class="s">}</span>
+ 146 
+ 147 <span class="c"># Process compound being ignored due to problems in fingerprints geneation...</span>
+ 148 <span class="c">#</span>
+<a name="ProcessIgnoredCompound-"></a> 149 <span class="k">sub </span><span class="m">ProcessIgnoredCompound</span> <span class="s">{</span>
+ 150   <span class="k">my</span><span class="s">(</span><span class="i">$Mode</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 151   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdID</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 152 
+ 153   <span class="i">$DataFieldLabelAndValuesRef</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetDataFieldLabelAndValues</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 154   <span class="i">$CmpdID</span> = <span class="i">SetupCmpdIDForOutputFiles</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 155 
+ 156   <span class="j">MODE:</span> <span class="s">{</span>
+ 157     <span class="k">if</span> <span class="s">(</span><span class="i">$Mode</span> =~ <span class="q">/^ContainsNonElementalData$/i</span><span class="s">)</span> <span class="s">{</span>
+ 158       <span class="k">warn</span> <span class="q">&quot;\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Compound contains atom data corresponding to non-elemental atom symbol(s)...\n\n&quot;</span><span class="sc">;</span>
+ 159       <span class="k">next</span> <span class="j">MODE</span><span class="sc">;</span>
+ 160     <span class="s">}</span>
+ 161 
+ 162     <span class="k">if</span> <span class="s">(</span><span class="i">$Mode</span> =~ <span class="q">/^ContainsNoElementalData$/i</span><span class="s">)</span> <span class="s">{</span>
+ 163       <span class="k">warn</span> <span class="q">&quot;\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Compound contains no atom data...\n\n&quot;</span><span class="sc">;</span>
+ 164       <span class="k">next</span> <span class="j">MODE</span><span class="sc">;</span>
+ 165     <span class="s">}</span>
+ 166 
+ 167     <span class="k">if</span> <span class="s">(</span><span class="i">$Mode</span> =~ <span class="q">/^FingerprintsGenerationFailed$/i</span><span class="s">)</span> <span class="s">{</span>
+ 168       <span class="k">warn</span> <span class="q">&quot;\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Fingerprints generation didn&#39;t succeed...\n\n&quot;</span><span class="sc">;</span>
+ 169       <span class="k">next</span> <span class="j">MODE</span><span class="sc">;</span>
+ 170     <span class="s">}</span>
+ 171     <span class="k">warn</span> <span class="q">&quot;\nWarning: Ignoring compound record number $CmpdCount with ID $CmpdID: Fingerprints generation didn&#39;t succeed...\n\n&quot;</span><span class="sc">;</span>
+ 172   <span class="s">}</span>
+ 173 <span class="s">}</span>
+ 174 
+ 175 <span class="c"># Check and filter compounds....</span>
+ 176 <span class="c">#</span>
+<a name="CheckAndFilterCompound-"></a> 177 <span class="k">sub </span><span class="m">CheckAndFilterCompound</span> <span class="s">{</span>
+ 178   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 179   <span class="k">my</span><span class="s">(</span><span class="i">$ElementCount</span><span class="cm">,</span> <span class="i">$NonElementCount</span><span class="s">)</span><span class="sc">;</span>
+ 180 
+ 181   <span class="s">(</span><span class="i">$ElementCount</span><span class="cm">,</span> <span class="i">$NonElementCount</span><span class="s">)</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetNumOfElementsAndNonElements</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 182 
+ 183   <span class="k">if</span> <span class="s">(</span><span class="i">$NonElementCount</span><span class="s">)</span> <span class="s">{</span>
+ 184     <span class="i">ProcessIgnoredCompound</span><span class="s">(</span><span class="q">&#39;ContainsNonElementalData&#39;</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span><span class="sc">;</span>
+ 185     <span class="k">return</span> <span class="n">1</span><span class="sc">;</span>
+ 186   <span class="s">}</span>
+ 187 
+ 188   <span class="k">if</span> <span class="s">(</span>!<span class="i">$ElementCount</span><span class="s">)</span> <span class="s">{</span>
+ 189     <span class="i">ProcessIgnoredCompound</span><span class="s">(</span><span class="q">&#39;ContainsNoElementalData&#39;</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="s">)</span><span class="sc">;</span>
+ 190     <span class="k">return</span> <span class="n">1</span><span class="sc">;</span>
+ 191   <span class="s">}</span>
+ 192 
+ 193   <span class="k">return</span> <span class="n">0</span><span class="sc">;</span>
+ 194 <span class="s">}</span>
+ 195 
+ 196 <span class="c"># Write out compounds fingerprints generation summary statistics...</span>
+ 197 <span class="c">#</span>
+<a name="WriteFingerprintsGenerationSummaryStatistics-"></a> 198 <span class="k">sub </span><span class="m">WriteFingerprintsGenerationSummaryStatistics</span> <span class="s">{</span>
+ 199   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$IgnoredCmpdCount</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 200   <span class="k">my</span><span class="s">(</span><span class="i">$ProcessedCmpdCount</span><span class="s">)</span><span class="sc">;</span>
+ 201 
+ 202   <span class="i">$ProcessedCmpdCount</span> = <span class="i">$CmpdCount</span> - <span class="i">$IgnoredCmpdCount</span><span class="sc">;</span>
+ 203 
+ 204   <span class="k">print</span> <span class="q">&quot;\nNumber of compounds: $CmpdCount\n&quot;</span><span class="sc">;</span>
+ 205   <span class="k">print</span> <span class="q">&quot;Number of compounds processed successfully during fingerprints generation: $ProcessedCmpdCount\n&quot;</span><span class="sc">;</span>
+ 206   <span class="k">print</span> <span class="q">&quot;Number of compounds ignored during fingerprints generation: $IgnoredCmpdCount\n&quot;</span><span class="sc">;</span>
+ 207 <span class="s">}</span>
+ 208 
+ 209 <span class="c"># Open output files...</span>
+ 210 <span class="c">#</span>
+<a name="SetupAndOpenOutputFiles-"></a> 211 <span class="k">sub </span><span class="m">SetupAndOpenOutputFiles</span> <span class="s">{</span>
+ 212   <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 213   <span class="k">my</span><span class="s">(</span><span class="i">$NewFPSDFile</span><span class="cm">,</span> <span class="i">$NewFPFile</span><span class="cm">,</span> <span class="i">$NewFPTextFile</span><span class="cm">,</span> <span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="cm">,</span> <span class="i">%FingerprintsFileIOParams</span><span class="s">)</span><span class="sc">;</span>
+ 214 
+ 215   <span class="s">(</span><span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span> = <span class="s">(</span><span class="k">undef</span><span class="s">)</span> x <span class="n">3</span><span class="sc">;</span>
+ 216 
+ 217   <span class="c"># Setup common parameters for fingerprints file IO objects...</span>
+ 218   <span class="c">#</span>
+ 219   <span class="i">%FingerprintsFileIOParams</span> = <span class="s">(</span><span class="q">&#39;Mode&#39;</span> <span class="cm">=&gt;</span> <span class="q">&#39;Write&#39;</span><span class="cm">,</span> <span class="q">&#39;Overwrite&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">OverwriteFiles</span>}<span class="cm">,</span> <span class="q">&#39;FingerprintsStringMode&#39;</span> <span class="cm">=&gt;</span> <span class="q">&#39;FingerprintsVectorString&#39;</span><span class="cm">,</span> <span class="q">&#39;VectorStringFormat&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">VectorStringFormat</span>}<span class="s">)</span><span class="sc">;</span>
+ 220 
+ 221   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">SDOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 222     <span class="i">$NewFPSDFile</span> = <span class="i">$SDFilesInfo</span>{<span class="w">SDOutFileNames</span>}[<span class="i">$FileIndex</span>]<span class="sc">;</span>
+ 223     <span class="k">print</span> <span class="q">&quot;Generating SD file $NewFPSDFile...\n&quot;</span><span class="sc">;</span>
+ 224     <span class="i">$NewFPSDFileIO</span> = <span class="w">new</span> <span class="i">FileIO::FingerprintsSDFileIO</span><span class="s">(</span><span class="q">&#39;Name&#39;</span> <span class="cm">=&gt;</span> <span class="i">$NewFPSDFile</span><span class="cm">,</span> <span class="i">%FingerprintsFileIOParams</span><span class="cm">,</span> <span class="q">&#39;FingerprintsFieldLabel&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">FingerprintsLabel</span>}<span class="s">)</span><span class="sc">;</span>
+ 225     <span class="i">$NewFPSDFileIO</span><span class="i">-&gt;Open</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 226   <span class="s">}</span>
+ 227 
+ 228   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">FPOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 229     <span class="i">$NewFPFile</span> = <span class="i">$SDFilesInfo</span>{<span class="w">FPOutFileNames</span>}[<span class="i">$FileIndex</span>]<span class="sc">;</span>
+ 230     <span class="k">print</span> <span class="q">&quot;Generating FP file $NewFPFile...\n&quot;</span><span class="sc">;</span>
+ 231     <span class="i">$NewFPFileIO</span> = <span class="w">new</span> <span class="i">FileIO::FingerprintsFPFileIO</span><span class="s">(</span><span class="q">&#39;Name&#39;</span> <span class="cm">=&gt;</span> <span class="i">$NewFPFile</span><span class="cm">,</span> <span class="i">%FingerprintsFileIOParams</span><span class="s">)</span><span class="sc">;</span>
+ 232     <span class="i">$NewFPFileIO</span><span class="i">-&gt;Open</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 233   <span class="s">}</span>
+ 234 
+ 235   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">TextOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 236     <span class="k">my</span><span class="s">(</span><span class="i">$ColLabelsRef</span><span class="s">)</span><span class="sc">;</span>
+ 237 
+ 238     <span class="i">$NewFPTextFile</span> = <span class="i">$SDFilesInfo</span>{<span class="w">TextOutFileNames</span>}[<span class="i">$FileIndex</span>]<span class="sc">;</span>
+ 239     <span class="i">$ColLabelsRef</span> = <span class="i">SetupFPTextFileCoulmnLabels</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span>
+ 240 
+ 241     <span class="k">print</span> <span class="q">&quot;Generating text file $NewFPTextFile...\n&quot;</span><span class="sc">;</span>
+ 242     <span class="i">$NewFPTextFileIO</span> = <span class="w">new</span> <span class="i">FileIO::FingerprintsTextFileIO</span><span class="s">(</span><span class="q">&#39;Name&#39;</span> <span class="cm">=&gt;</span> <span class="i">$NewFPTextFile</span><span class="cm">,</span> <span class="i">%FingerprintsFileIOParams</span><span class="cm">,</span> <span class="q">&#39;DataColLabels&#39;</span> <span class="cm">=&gt;</span> <span class="i">$ColLabelsRef</span><span class="cm">,</span> <span class="q">&#39;OutDelim&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">OutDelim</span>}<span class="cm">,</span> <span class="q">&#39;OutQuote&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">OutQuote</span>}<span class="s">)</span><span class="sc">;</span>
+ 243     <span class="i">$NewFPTextFileIO</span><span class="i">-&gt;Open</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 244   <span class="s">}</span>
+ 245 
+ 246   <span class="k">return</span> <span class="s">(</span><span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span><span class="sc">;</span>
+ 247 <span class="s">}</span>
+ 248 
+ 249 <span class="c"># Write fingerpritns and other data to appropriate output files...</span>
+ 250 <span class="c">#</span>
+<a name="WriteDataToOutputFiles-"></a> 251 <span class="k">sub </span><span class="m">WriteDataToOutputFiles</span> <span class="s">{</span>
+ 252   <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$NewFPSDFileIO</span><span class="cm">,</span> <span class="i">$NewFPTextFileIO</span><span class="cm">,</span> <span class="i">$NewFPFileIO</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 253   <span class="k">my</span><span class="s">(</span><span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 254 
+ 255   <span class="i">$DataFieldLabelAndValuesRef</span> = <span class="k">undef</span><span class="sc">;</span>
+ 256   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPTextFileIO</span> || <span class="i">$NewFPFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 257     <span class="i">$DataFieldLabelAndValuesRef</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetDataFieldLabelAndValues</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 258   <span class="s">}</span>
+ 259 
+ 260   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPSDFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 261     <span class="k">my</span><span class="s">(</span><span class="i">$CmpdString</span><span class="s">)</span><span class="sc">;</span>
+ 262 
+ 263     <span class="i">$CmpdString</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetInputMoleculeString</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 264     <span class="i">$NewFPSDFileIO</span><span class="i">-&gt;WriteFingerprints</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$CmpdString</span><span class="s">)</span><span class="sc">;</span>
+ 265   <span class="s">}</span>
+ 266 
+ 267   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPTextFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 268     <span class="k">my</span><span class="s">(</span><span class="i">$ColValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 269 
+ 270     <span class="i">$ColValuesRef</span> = <span class="i">SetupFPTextFileCoulmnValues</span><span class="s">(</span><span class="i">$FileIndex</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 271     <span class="i">$NewFPTextFileIO</span><span class="i">-&gt;WriteFingerprints</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$ColValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 272   <span class="s">}</span>
+ 273 
+ 274   <span class="k">if</span> <span class="s">(</span><span class="i">$NewFPFileIO</span><span class="s">)</span> <span class="s">{</span>
+ 275     <span class="k">my</span><span class="s">(</span><span class="i">$CompoundID</span><span class="s">)</span><span class="sc">;</span>
+ 276 
+ 277     <span class="i">$CompoundID</span> = <span class="i">SetupCmpdIDForOutputFiles</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 278     <span class="i">$NewFPFileIO</span><span class="i">-&gt;WriteFingerprints</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="cm">,</span> <span class="i">$CompoundID</span><span class="s">)</span><span class="sc">;</span>
+ 279   <span class="s">}</span>
+ 280 <span class="s">}</span>
+ 281 
+ 282 <span class="c"># Generate approriate column labels for FPText output file...</span>
+ 283 <span class="c">#</span>
+<a name="SetupFPTextFileCoulmnLabels-"></a> 284 <span class="k">sub </span><span class="m">SetupFPTextFileCoulmnLabels</span> <span class="s">{</span>
+ 285   <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 286   <span class="k">my</span><span class="s">(</span><span class="i">$Line</span><span class="cm">,</span> <span class="i">@ColLabels</span><span class="s">)</span><span class="sc">;</span>
+ 287 
+ 288   <span class="i">@ColLabels</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 289   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^All$/i</span><span class="s">)</span> <span class="s">{</span>
+ 290     <span class="k">push</span> <span class="i">@ColLabels</span><span class="cm">,</span> <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">AllDataFieldsRef</span>}[<span class="i">$FileIndex</span>]}<span class="sc">;</span>
+ 291   <span class="s">}</span>
+ 292   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^Common$/i</span><span class="s">)</span> <span class="s">{</span>
+ 293     <span class="k">push</span> <span class="i">@ColLabels</span><span class="cm">,</span> <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">CommonDataFieldsRef</span>}[<span class="i">$FileIndex</span>]}<span class="sc">;</span>
+ 294   <span class="s">}</span>
+ 295   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^Specify$/i</span><span class="s">)</span> <span class="s">{</span>
+ 296     <span class="k">push</span> <span class="i">@ColLabels</span><span class="cm">,</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedDataFields</span>}}<span class="sc">;</span>
+ 297   <span class="s">}</span>
+ 298   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^CompoundID$/i</span><span class="s">)</span> <span class="s">{</span>
+ 299     <span class="k">push</span> <span class="i">@ColLabels</span><span class="cm">,</span> <span class="i">$OptionsInfo</span>{<span class="w">CompoundIDLabel</span>}<span class="sc">;</span>
+ 300   <span class="s">}</span>
+ 301   <span class="c"># Add fingerprints label...</span>
+ 302   <span class="k">push</span> <span class="i">@ColLabels</span><span class="cm">,</span> <span class="i">$OptionsInfo</span>{<span class="w">FingerprintsLabel</span>}<span class="sc">;</span>
+ 303 
+ 304   <span class="k">return</span> \<span class="i">@ColLabels</span><span class="sc">;</span>
+ 305 <span class="s">}</span>
+ 306 
+ 307 <span class="c"># Generate column values FPText output file..</span>
+ 308 <span class="c">#</span>
+<a name="SetupFPTextFileCoulmnValues-"></a> 309 <span class="k">sub </span><span class="m">SetupFPTextFileCoulmnValues</span> <span class="s">{</span>
+ 310   <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="cm">,</span> <span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 311   <span class="k">my</span><span class="s">(</span><span class="i">@ColValues</span><span class="s">)</span><span class="sc">;</span>
+ 312 
+ 313   <span class="i">@ColValues</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 314   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^CompoundID$/i</span><span class="s">)</span> <span class="s">{</span>
+ 315     <span class="k">push</span> <span class="i">@ColValues</span><span class="cm">,</span> <span class="i">SetupCmpdIDForOutputFiles</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span><span class="sc">;</span>
+ 316   <span class="s">}</span>
+ 317   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^All$/i</span><span class="s">)</span> <span class="s">{</span>
+ 318     <span class="i">@ColValues</span> = <span class="k">map</span> <span class="s">{</span> <span class="k">exists</span> <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} ? <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} <span class="co">:</span> <span class="q">&#39;&#39;</span><span class="s">}</span> <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">AllDataFieldsRef</span>}[<span class="i">$FileIndex</span>]}<span class="sc">;</span>
+ 319   <span class="s">}</span>
+ 320   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^Common$/i</span><span class="s">)</span> <span class="s">{</span>
+ 321     <span class="i">@ColValues</span> = <span class="k">map</span> <span class="s">{</span> <span class="k">exists</span> <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} ? <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} <span class="co">:</span> <span class="q">&#39;&#39;</span><span class="s">}</span> <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">CommonDataFieldsRef</span>}[<span class="i">$FileIndex</span>]}<span class="sc">;</span>
+ 322   <span class="s">}</span>
+ 323   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^Specify$/i</span><span class="s">)</span> <span class="s">{</span>
+ 324     <span class="i">@ColValues</span> = <span class="k">map</span> <span class="s">{</span> <span class="k">exists</span> <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} ? <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$_</span>} <span class="co">:</span> <span class="q">&#39;&#39;</span><span class="s">}</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedDataFields</span>}}<span class="sc">;</span>
+ 325   <span class="s">}</span>
+ 326 
+ 327   <span class="k">return</span> \<span class="i">@ColValues</span><span class="sc">;</span>
+ 328 <span class="s">}</span>
+ 329 
+ 330 <span class="c"># Generate compound ID for FP and FPText output files..</span>
+ 331 <span class="c">#</span>
+<a name="SetupCmpdIDForOutputFiles-"></a> 332 <span class="k">sub </span><span class="m">SetupCmpdIDForOutputFiles</span> <span class="s">{</span>
+ 333   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="i">$DataFieldLabelAndValuesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 334   <span class="k">my</span><span class="s">(</span><span class="i">$CmpdID</span><span class="s">)</span><span class="sc">;</span>
+ 335 
+ 336   <span class="i">$CmpdID</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 337   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} =~ <span class="q">/^MolNameOrLabelPrefix$/i</span><span class="s">)</span> <span class="s">{</span>
+ 338     <span class="k">my</span><span class="s">(</span><span class="i">$MolName</span><span class="s">)</span><span class="sc">;</span>
+ 339     <span class="i">$MolName</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetName</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 340     <span class="i">$CmpdID</span> = <span class="i">$MolName</span> ? <span class="i">$MolName</span> <span class="co">:</span> <span class="q">&quot;$OptionsInfo{CompoundID}${CmpdCount}&quot;</span><span class="sc">;</span>
+ 341   <span class="s">}</span>
+ 342   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} =~ <span class="q">/^LabelPrefix$/i</span><span class="s">)</span> <span class="s">{</span>
+ 343     <span class="i">$CmpdID</span> = <span class="q">&quot;$OptionsInfo{CompoundID}${CmpdCount}&quot;</span><span class="sc">;</span>
+ 344   <span class="s">}</span>
+ 345   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} =~ <span class="q">/^DataField$/i</span><span class="s">)</span> <span class="s">{</span>
+ 346     <span class="k">my</span><span class="s">(</span><span class="i">$SpecifiedDataField</span><span class="s">)</span><span class="sc">;</span>
+ 347     <span class="i">$SpecifiedDataField</span> = <span class="i">$OptionsInfo</span>{<span class="w">CompoundID</span>}<span class="sc">;</span>
+ 348     <span class="i">$CmpdID</span> = <span class="k">exists</span> <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$SpecifiedDataField</span>} ? <span class="i">$DataFieldLabelAndValuesRef</span>-&gt;{<span class="i">$SpecifiedDataField</span>} <span class="co">:</span> <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 349   <span class="s">}</span>
+ 350   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} =~ <span class="q">/^MolName$/i</span><span class="s">)</span> <span class="s">{</span>
+ 351     <span class="i">$CmpdID</span> = <span class="i">$Molecule</span><span class="i">-&gt;GetName</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 352   <span class="s">}</span>
+ 353   <span class="k">return</span> <span class="i">$CmpdID</span><span class="sc">;</span>
+ 354 <span class="s">}</span>
+ 355 
+ 356 <span class="c"># Generate fingerprints for molecule...</span>
+ 357 <span class="c">#</span>
+<a name="GenerateMoleculeFingerprints-"></a> 358 <span class="k">sub </span><span class="m">GenerateMoleculeFingerprints</span> <span class="s">{</span>
+ 359   <span class="k">my</span><span class="s">(</span><span class="i">$Molecule</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 360   <span class="k">my</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="s">)</span><span class="sc">;</span>
+ 361 
+ 362   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">KeepLargestComponent</span>}<span class="s">)</span> <span class="s">{</span>
+ 363     <span class="i">$Molecule</span><span class="i">-&gt;KeepLargestComponent</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 364   <span class="s">}</span>
+ 365   <span class="k">if</span> <span class="s">(</span>!<span class="i">$Molecule</span><span class="i">-&gt;DetectRings</span><span class="s">(</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 366     <span class="k">return</span> <span class="k">undef</span><span class="sc">;</span>
+ 367   <span class="s">}</span>
+ 368   <span class="i">$Molecule</span><span class="i">-&gt;SetAromaticityModel</span><span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">AromaticityModel</span>}<span class="s">)</span><span class="sc">;</span>
+ 369   <span class="i">$Molecule</span><span class="i">-&gt;DetectAromaticity</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 370 
+ 371   <span class="i">$TopologicalAtomPairsFingerprints</span> = <span class="w">new</span> <span class="i">Fingerprints::TopologicalAtomPairsFingerprints</span><span class="s">(</span><span class="q">&#39;Molecule&#39;</span> <span class="cm">=&gt;</span> <span class="i">$Molecule</span><span class="cm">,</span> <span class="q">&#39;MinDistance&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">MinDistance</span>}<span class="cm">,</span>  <span class="q">&#39;MaxDistance&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">MaxDistance</span>}<span class="cm">,</span> <span class="q">&#39;AtomIdentifierType&#39;</span> <span class="cm">=&gt;</span> <span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>}<span class="s">)</span><span class="sc">;</span>
+ 372   <span class="i">SetAtomIdentifierTypeValuesToUse</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="s">)</span><span class="sc">;</span>
+ 373 
+ 374   <span class="c"># Generate fingerprints...</span>
+ 375   <span class="i">$TopologicalAtomPairsFingerprints</span><span class="i">-&gt;GenerateFingerprints</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 376 
+ 377   <span class="c"># Make sure fingerprints generation is successful...</span>
+ 378   <span class="k">if</span> <span class="s">(</span>!<span class="i">$TopologicalAtomPairsFingerprints</span><span class="i">-&gt;IsFingerprintsGenerationSuccessful</span><span class="s">(</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 379     <span class="k">return</span> <span class="k">undef</span><span class="sc">;</span>
+ 380   <span class="s">}</span>
+ 381 
+ 382   <span class="k">return</span> <span class="i">$TopologicalAtomPairsFingerprints</span><span class="sc">;</span>
+ 383 <span class="s">}</span>
+ 384 
+ 385 <span class="c"># Set atom identifier type to use for generating fingerprints...</span>
+ 386 <span class="c">#</span>
+<a name="SetAtomIdentifierTypeValuesToUse-"></a> 387 <span class="k">sub </span><span class="m">SetAtomIdentifierTypeValuesToUse</span> <span class="s">{</span>
+ 388   <span class="k">my</span><span class="s">(</span><span class="i">$TopologicalAtomPairsFingerprints</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span>
+ 389 
+ 390   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>} =~ <span class="q">/^AtomicInvariantsAtomTypes$/i</span><span class="s">)</span> <span class="s">{</span>
+ 391     <span class="i">$TopologicalAtomPairsFingerprints</span><span class="i">-&gt;SetAtomicInvariantsToUse</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">AtomicInvariantsToUse</span>}}<span class="s">)</span><span class="sc">;</span>
+ 392   <span class="s">}</span>
+ 393   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>} =~ <span class="q">/^FunctionalClassAtomTypes$/i</span><span class="s">)</span> <span class="s">{</span>
+ 394     <span class="i">$TopologicalAtomPairsFingerprints</span><span class="i">-&gt;SetFunctionalClassesToUse</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">FunctionalClassesToUse</span>}}<span class="s">)</span><span class="sc">;</span>
+ 395   <span class="s">}</span>
+ 396   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>} =~ <span class="q">/^(DREIDINGAtomTypes|EStateAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 397     <span class="c"># Nothing to do for now...</span>
+ 398   <span class="s">}</span>
+ 399   <span class="k">else</span> <span class="s">{</span>
+ 400     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{atomidentifiertype}, for option \&quot;-a, --AtomIdentifierType\&quot; is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n&quot;</span><span class="sc">;</span>
+ 401   <span class="s">}</span>
+ 402 <span class="s">}</span>
+ 403 
+ 404 <span class="c"># Retrieve information about SD files...</span>
+ 405 <span class="c">#</span>
+<a name="RetrieveSDFilesInfo-"></a> 406 <span class="k">sub </span><span class="m">RetrieveSDFilesInfo</span> <span class="s">{</span>
+ 407   <span class="k">my</span><span class="s">(</span><span class="i">$SDFile</span><span class="cm">,</span> <span class="i">$Index</span><span class="cm">,</span> <span class="i">$FileDir</span><span class="cm">,</span> <span class="i">$FileExt</span><span class="cm">,</span> <span class="i">$FileName</span><span class="cm">,</span> <span class="i">$OutFileRoot</span><span class="cm">,</span> <span class="i">$TextOutFileExt</span><span class="cm">,</span> <span class="i">$SDOutFileExt</span><span class="cm">,</span> <span class="i">$FPOutFileExt</span><span class="cm">,</span> <span class="i">$NewSDFileName</span><span class="cm">,</span> <span class="i">$NewFPFileName</span><span class="cm">,</span> <span class="i">$NewTextFileName</span><span class="cm">,</span> <span class="i">$CheckDataField</span><span class="cm">,</span> <span class="i">$CollectDataFields</span><span class="cm">,</span> <span class="i">$AllDataFieldsRef</span><span class="cm">,</span> <span class="i">$CommonDataFieldsRef</span><span class="s">)</span><span class="sc">;</span>
+ 408 
+ 409   <span class="i">%SDFilesInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 410   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">FileOkay</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 411   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">OutFileRoot</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 412   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">SDOutFileNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 413   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">FPOutFileNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 414   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">TextOutFileNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 415   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">AllDataFieldsRef</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 416   <span class="i">@</span>{<span class="i">$SDFilesInfo</span>{<span class="w">CommonDataFieldsRef</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 417 
+ 418   <span class="i">$CheckDataField</span> = <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">TextOutput</span>} &amp;&amp; <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^CompoundID$/i</span><span class="s">)</span> &amp;&amp; <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} =~ <span class="q">/^DataField$/i</span><span class="s">)</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 419   <span class="i">$CollectDataFields</span> = <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">TextOutput</span>} &amp;&amp; <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} =~ <span class="q">/^(All|Common)$/i</span><span class="s">)</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 420 
+ 421   <span class="j">FILELIST:</span> <span class="k">for</span> <span class="i">$Index</span> <span class="s">(</span><span class="n">0</span> .. <span class="i">$#SDFilesList</span><span class="s">)</span> <span class="s">{</span>
+ 422     <span class="i">$SDFile</span> = <span class="i">$SDFilesList</span>[<span class="i">$Index</span>]<span class="sc">;</span>
+ 423 
+ 424     <span class="i">$SDFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$Index</span>] = <span class="n">0</span><span class="sc">;</span>
+ 425     <span class="i">$SDFilesInfo</span>{<span class="w">OutFileRoot</span>}[<span class="i">$Index</span>] = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 426     <span class="i">$SDFilesInfo</span>{<span class="w">SDOutFileNames</span>}[<span class="i">$Index</span>] = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 427     <span class="i">$SDFilesInfo</span>{<span class="w">FPOutFileNames</span>}[<span class="i">$Index</span>] = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 428     <span class="i">$SDFilesInfo</span>{<span class="w">TextOutFileNames</span>}[<span class="i">$Index</span>] = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 429 
+ 430     <span class="i">$SDFile</span> = <span class="i">$SDFilesList</span>[<span class="i">$Index</span>]<span class="sc">;</span>
+ 431     <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">-e</span> <span class="i">$SDFile</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 432       <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: It doesn&#39;t exist\n&quot;</span><span class="sc">;</span>
+ 433       <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 434     <span class="s">}</span>
+ 435     <span class="k">if</span> <span class="s">(</span>!<span class="i">CheckFileType</span><span class="s">(</span><span class="i">$SDFile</span><span class="cm">,</span> <span class="q">&quot;sd sdf&quot;</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 436       <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: It&#39;s not a SD file\n&quot;</span><span class="sc">;</span>
+ 437       <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 438     <span class="s">}</span>
+ 439 
+ 440     <span class="k">if</span> <span class="s">(</span><span class="i">$CheckDataField</span><span class="s">)</span> <span class="s">{</span>
+ 441       <span class="c"># Make sure data field exists in SD file..</span>
+ 442       <span class="k">my</span><span class="s">(</span><span class="i">$CmpdString</span><span class="cm">,</span> <span class="i">$SpecifiedDataField</span><span class="cm">,</span> <span class="i">@CmpdLines</span><span class="cm">,</span> <span class="i">%DataFieldValues</span><span class="s">)</span><span class="sc">;</span>
+ 443 
+ 444       <span class="i">@CmpdLines</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 445       <span class="k">open</span> <span class="w">SDFILE</span><span class="cm">,</span> <span class="q">&quot;$SDFile&quot;</span> <span class="k">or</span> <span class="k">die</span> <span class="q">&quot;Error: Couldn&#39;t open $SDFile: $! \n&quot;</span><span class="sc">;</span>
+ 446       <span class="i">$CmpdString</span> = <span class="i">ReadCmpdString</span><span class="s">(</span>\<span class="i">*SDFILE</span><span class="s">)</span><span class="sc">;</span>
+ 447       <span class="k">close</span> <span class="w">SDFILE</span><span class="sc">;</span>
+ 448       <span class="i">@CmpdLines</span> = <span class="k">split</span> <span class="q">&quot;\n&quot;</span><span class="cm">,</span> <span class="i">$CmpdString</span><span class="sc">;</span>
+ 449       <span class="i">%DataFieldValues</span> = <span class="i">GetCmpdDataHeaderLabelsAndValues</span><span class="s">(</span>\<span class="i">@CmpdLines</span><span class="s">)</span><span class="sc">;</span>
+ 450       <span class="i">$SpecifiedDataField</span> = <span class="i">$OptionsInfo</span>{<span class="w">CompoundID</span>}<span class="sc">;</span>
+ 451       <span class="k">if</span> <span class="s">(</span>!<span class="k">exists</span> <span class="i">$DataFieldValues</span>{<span class="i">$SpecifiedDataField</span>}<span class="s">)</span> <span class="s">{</span>
+ 452         <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: Data field value, $SpecifiedDataField, using  \&quot;--CompoundID\&quot; option in \&quot;DataField\&quot; \&quot;--CompoundIDMode\&quot; doesn&#39;t exist\n&quot;</span><span class="sc">;</span>
+ 453         <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 454       <span class="s">}</span>
+ 455     <span class="s">}</span>
+ 456 
+ 457     <span class="i">$AllDataFieldsRef</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 458     <span class="i">$CommonDataFieldsRef</span> = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 459     <span class="k">if</span> <span class="s">(</span><span class="i">$CollectDataFields</span><span class="s">)</span> <span class="s">{</span>
+ 460       <span class="k">my</span><span class="s">(</span><span class="i">$CmpdCount</span><span class="s">)</span><span class="sc">;</span>
+ 461       <span class="k">open</span> <span class="w">SDFILE</span><span class="cm">,</span> <span class="q">&quot;$SDFile&quot;</span> <span class="k">or</span> <span class="k">die</span> <span class="q">&quot;Error: Couldn&#39;t open $SDFile: $! \n&quot;</span><span class="sc">;</span>
+ 462       <span class="s">(</span><span class="i">$CmpdCount</span><span class="cm">,</span> <span class="i">$AllDataFieldsRef</span><span class="cm">,</span> <span class="i">$CommonDataFieldsRef</span><span class="s">)</span> = <span class="i">GetAllAndCommonCmpdDataHeaderLabels</span><span class="s">(</span>\<span class="i">*SDFILE</span><span class="s">)</span><span class="sc">;</span>
+ 463       <span class="k">close</span> <span class="w">SDFILE</span><span class="sc">;</span>
+ 464     <span class="s">}</span>
+ 465 
+ 466     <span class="c"># Setup output file names...</span>
+ 467     <span class="i">$FileDir</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span> <span class="i">$FileName</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span> <span class="i">$FileExt</span> = <span class="q">&quot;&quot;</span><span class="sc">;</span>
+ 468     <span class="s">(</span><span class="i">$FileDir</span><span class="cm">,</span> <span class="i">$FileName</span><span class="cm">,</span> <span class="i">$FileExt</span><span class="s">)</span> = <span class="i">ParseFileName</span><span class="s">(</span><span class="i">$SDFile</span><span class="s">)</span><span class="sc">;</span>
+ 469 
+ 470     <span class="i">$TextOutFileExt</span> = <span class="q">&quot;csv&quot;</span><span class="sc">;</span>
+ 471     <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">outdelim</span>} =~ <span class="q">/^tab$/i</span><span class="s">)</span> <span class="s">{</span>
+ 472       <span class="i">$TextOutFileExt</span> = <span class="q">&quot;tsv&quot;</span><span class="sc">;</span>
+ 473     <span class="s">}</span>
+ 474     <span class="i">$SDOutFileExt</span> = <span class="i">$FileExt</span><span class="sc">;</span>
+ 475     <span class="i">$FPOutFileExt</span> = <span class="q">&quot;fpf&quot;</span><span class="sc">;</span>
+ 476 
+ 477     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">OutFileRoot</span>} &amp;&amp; <span class="s">(</span><span class="i">@SDFilesList</span> == <span class="n">1</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 478       <span class="k">my</span> <span class="s">(</span><span class="i">$RootFileDir</span><span class="cm">,</span> <span class="i">$RootFileName</span><span class="cm">,</span> <span class="i">$RootFileExt</span><span class="s">)</span> = <span class="i">ParseFileName</span><span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">OutFileRoot</span>}<span class="s">)</span><span class="sc">;</span>
+ 479       <span class="k">if</span> <span class="s">(</span><span class="i">$RootFileName</span> &amp;&amp; <span class="i">$RootFileExt</span><span class="s">)</span> <span class="s">{</span>
+ 480         <span class="i">$FileName</span> = <span class="i">$RootFileName</span><span class="sc">;</span>
+ 481       <span class="s">}</span>
+ 482       <span class="k">else</span> <span class="s">{</span>
+ 483         <span class="i">$FileName</span> = <span class="i">$OptionsInfo</span>{<span class="w">OutFileRoot</span>}<span class="sc">;</span>
+ 484       <span class="s">}</span>
+ 485       <span class="i">$OutFileRoot</span> = <span class="i">$FileName</span><span class="sc">;</span>
+ 486     <span class="s">}</span>
+ 487     <span class="k">else</span> <span class="s">{</span>
+ 488       <span class="i">$OutFileRoot</span> = <span class="q">&quot;${FileName}TopologicalAtomPairsFP&quot;</span><span class="sc">;</span>
+ 489     <span class="s">}</span>
+ 490 
+ 491     <span class="i">$NewSDFileName</span> = <span class="q">&quot;${OutFileRoot}.${SDOutFileExt}&quot;</span><span class="sc">;</span>
+ 492     <span class="i">$NewFPFileName</span> = <span class="q">&quot;${OutFileRoot}.${FPOutFileExt}&quot;</span><span class="sc">;</span>
+ 493     <span class="i">$NewTextFileName</span> = <span class="q">&quot;${OutFileRoot}.${TextOutFileExt}&quot;</span><span class="sc">;</span>
+ 494 
+ 495     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">SDOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 496       <span class="k">if</span> <span class="s">(</span><span class="i">$SDFile</span> =~ <span class="q">/$NewSDFileName/i</span><span class="s">)</span> <span class="s">{</span>
+ 497         <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring input file $SDFile: Same output, $NewSDFileName, and input file names.\n&quot;</span><span class="sc">;</span>
+ 498         <span class="k">print</span> <span class="q">&quot;Specify a different name using \&quot;-r --root\&quot; option or use default name.\n&quot;</span><span class="sc">;</span>
+ 499         <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 500       <span class="s">}</span>
+ 501     <span class="s">}</span>
+ 502 
+ 503     <span class="k">if</span> <span class="s">(</span>!<span class="i">$OptionsInfo</span>{<span class="w">OverwriteFiles</span>}<span class="s">)</span> <span class="s">{</span>
+ 504       <span class="c"># Check SD and text outout files...</span>
+ 505       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">SDOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 506         <span class="k">if</span> <span class="s">(</span><span class="k">-e</span> <span class="i">$NewSDFileName</span><span class="s">)</span> <span class="s">{</span>
+ 507           <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: The file $NewSDFileName already exists\n&quot;</span><span class="sc">;</span>
+ 508           <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 509         <span class="s">}</span>
+ 510       <span class="s">}</span>
+ 511       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">FPOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 512         <span class="k">if</span> <span class="s">(</span><span class="k">-e</span> <span class="i">$NewFPFileName</span><span class="s">)</span> <span class="s">{</span>
+ 513           <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: The file $NewFPFileName already exists\n&quot;</span><span class="sc">;</span>
+ 514           <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 515         <span class="s">}</span>
+ 516       <span class="s">}</span>
+ 517       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">TextOutput</span>}<span class="s">)</span> <span class="s">{</span>
+ 518         <span class="k">if</span> <span class="s">(</span><span class="k">-e</span> <span class="i">$NewTextFileName</span><span class="s">)</span> <span class="s">{</span>
+ 519           <span class="k">warn</span> <span class="q">&quot;Warning: Ignoring file $SDFile: The file $NewTextFileName already exists\n&quot;</span><span class="sc">;</span>
+ 520           <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span>
+ 521         <span class="s">}</span>
+ 522       <span class="s">}</span>
+ 523     <span class="s">}</span>
+ 524 
+ 525     <span class="i">$SDFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$Index</span>] = <span class="n">1</span><span class="sc">;</span>
+ 526 
+ 527     <span class="i">$SDFilesInfo</span>{<span class="w">OutFileRoot</span>}[<span class="i">$Index</span>] = <span class="i">$OutFileRoot</span><span class="sc">;</span>
+ 528     <span class="i">$SDFilesInfo</span>{<span class="w">SDOutFileNames</span>}[<span class="i">$Index</span>] = <span class="i">$NewSDFileName</span><span class="sc">;</span>
+ 529     <span class="i">$SDFilesInfo</span>{<span class="w">FPOutFileNames</span>}[<span class="i">$Index</span>] = <span class="i">$NewFPFileName</span><span class="sc">;</span>
+ 530     <span class="i">$SDFilesInfo</span>{<span class="w">TextOutFileNames</span>}[<span class="i">$Index</span>] = <span class="i">$NewTextFileName</span><span class="sc">;</span>
+ 531 
+ 532     <span class="i">$SDFilesInfo</span>{<span class="w">AllDataFieldsRef</span>}[<span class="i">$Index</span>] = <span class="i">$AllDataFieldsRef</span><span class="sc">;</span>
+ 533     <span class="i">$SDFilesInfo</span>{<span class="w">CommonDataFieldsRef</span>}[<span class="i">$Index</span>] = <span class="i">$CommonDataFieldsRef</span><span class="sc">;</span>
+ 534   <span class="s">}</span>
+ 535 <span class="s">}</span>
+ 536 
+ 537 <span class="c"># Process option values...</span>
+<a name="ProcessOptions-"></a> 538 <span class="k">sub </span><span class="m">ProcessOptions</span> <span class="s">{</span>
+ 539   <span class="i">%OptionsInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 540 
+ 541   <span class="i">ProcessAtomIdentifierTypeOptions</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 542 
+ 543   <span class="i">$OptionsInfo</span>{<span class="w">AromaticityModel</span>} = <span class="i">$Options</span>{<span class="w">aromaticitymodel</span>}<span class="sc">;</span>
+ 544 
+ 545   <span class="i">$OptionsInfo</span>{<span class="w">CompoundIDMode</span>} = <span class="i">$Options</span>{<span class="w">compoundidmode</span>}<span class="sc">;</span>
+ 546   <span class="i">$OptionsInfo</span>{<span class="w">CompoundIDLabel</span>} = <span class="i">$Options</span>{<span class="w">compoundidlabel</span>}<span class="sc">;</span>
+ 547   <span class="i">$OptionsInfo</span>{<span class="w">DataFieldsMode</span>} = <span class="i">$Options</span>{<span class="w">datafieldsmode</span>}<span class="sc">;</span>
+ 548 
+ 549   <span class="k">my</span><span class="s">(</span><span class="i">@SpecifiedDataFields</span><span class="s">)</span><span class="sc">;</span>
+ 550   <span class="i">@SpecifiedDataFields</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 551 
+ 552   <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedDataFields</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 553   <span class="i">$OptionsInfo</span>{<span class="w">CompoundID</span>} = <span class="q">&#39;&#39;</span><span class="sc">;</span>
+ 554 
+ 555   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">datafieldsmode</span>} =~ <span class="q">/^CompoundID$/i</span><span class="s">)</span> <span class="s">{</span>
+ 556     <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">compoundidmode</span>} =~ <span class="q">/^DataField$/i</span><span class="s">)</span> <span class="s">{</span>
+ 557       <span class="k">if</span> <span class="s">(</span>!<span class="i">$Options</span>{<span class="w">compoundid</span>}<span class="s">)</span> <span class="s">{</span>
+ 558         <span class="k">die</span> <span class="q">&quot;Error: You must specify a value for \&quot;--CompoundID\&quot; option in \&quot;DataField\&quot; \&quot;--CompoundIDMode\&quot;. \n&quot;</span><span class="sc">;</span>
+ 559       <span class="s">}</span>
+ 560       <span class="i">$OptionsInfo</span>{<span class="w">CompoundID</span>} = <span class="i">$Options</span>{<span class="w">compoundid</span>}<span class="sc">;</span>
+ 561     <span class="s">}</span>
+ 562     <span class="k">elsif</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">compoundidmode</span>} =~ <span class="q">/^(LabelPrefix|MolNameOrLabelPrefix)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 563       <span class="i">$OptionsInfo</span>{<span class="w">CompoundID</span>} = <span class="i">$Options</span>{<span class="w">compoundid</span>} ? <span class="i">$Options</span>{<span class="w">compoundid</span>} <span class="co">:</span> <span class="q">&#39;Cmpd&#39;</span><span class="sc">;</span>
+ 564     <span class="s">}</span>
+ 565   <span class="s">}</span>
+ 566   <span class="k">elsif</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">datafieldsmode</span>} =~ <span class="q">/^Specify$/i</span><span class="s">)</span> <span class="s">{</span>
+ 567     <span class="k">if</span> <span class="s">(</span>!<span class="i">$Options</span>{<span class="w">datafields</span>}<span class="s">)</span> <span class="s">{</span>
+ 568       <span class="k">die</span> <span class="q">&quot;Error: You must specify a value for \&quot;--DataFields\&quot; option in \&quot;Specify\&quot; \&quot;-d, --DataFieldsMode\&quot;. \n&quot;</span><span class="sc">;</span>
+ 569     <span class="s">}</span>
+ 570     <span class="i">@SpecifiedDataFields</span> = <span class="k">split</span> <span class="q">/\,/</span><span class="cm">,</span> <span class="i">$Options</span>{<span class="w">datafields</span>}<span class="sc">;</span>
+ 571     <span class="k">push</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedDataFields</span>}}<span class="cm">,</span> <span class="i">@SpecifiedDataFields</span><span class="sc">;</span>
+ 572   <span class="s">}</span>
+ 573 
+ 574   <span class="i">$OptionsInfo</span>{<span class="w">Filter</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">filter</span>} =~ <span class="q">/^Yes$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 575 
+ 576   <span class="i">$OptionsInfo</span>{<span class="w">FingerprintsLabel</span>} = <span class="i">$Options</span>{<span class="w">fingerprintslabel</span>} ? <span class="i">$Options</span>{<span class="w">fingerprintslabel</span>} <span class="co">:</span> <span class="q">&#39;TopologicalAtomPairsFingerprints&#39;</span><span class="sc">;</span>
+ 577 
+ 578   <span class="i">$OptionsInfo</span>{<span class="w">KeepLargestComponent</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">keeplargestcomponent</span>} =~ <span class="q">/^Yes$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 579 
+ 580   <span class="i">$OptionsInfo</span>{<span class="w">MinDistance</span>} = <span class="i">$Options</span>{<span class="w">mindistance</span>}<span class="sc">;</span>
+ 581   <span class="i">$OptionsInfo</span>{<span class="w">MaxDistance</span>} = <span class="i">$Options</span>{<span class="w">maxdistance</span>}<span class="sc">;</span>
+ 582 
+ 583   <span class="i">$OptionsInfo</span>{<span class="w">Output</span>} = <span class="i">$Options</span>{<span class="w">output</span>}<span class="sc">;</span>
+ 584   <span class="i">$OptionsInfo</span>{<span class="w">SDOutput</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">output</span>} =~ <span class="q">/^(SD|All)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 585   <span class="i">$OptionsInfo</span>{<span class="w">FPOutput</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">output</span>} =~ <span class="q">/^(FP|All)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 586   <span class="i">$OptionsInfo</span>{<span class="w">TextOutput</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">output</span>} =~ <span class="q">/^(Text|All)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 587 
+ 588   <span class="i">$OptionsInfo</span>{<span class="w">OutDelim</span>} = <span class="i">$Options</span>{<span class="w">outdelim</span>}<span class="sc">;</span>
+ 589   <span class="i">$OptionsInfo</span>{<span class="w">OutQuote</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">quote</span>} =~ <span class="q">/^Yes$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 590 
+ 591   <span class="i">$OptionsInfo</span>{<span class="w">OverwriteFiles</span>} = <span class="i">$Options</span>{<span class="w">overwrite</span>} ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 592   <span class="i">$OptionsInfo</span>{<span class="w">OutFileRoot</span>} = <span class="i">$Options</span>{<span class="w">root</span>} ? <span class="i">$Options</span>{<span class="w">root</span>} <span class="co">:</span> <span class="n">0</span><span class="sc">;</span>
+ 593 
+ 594   <span class="i">$OptionsInfo</span>{<span class="w">VectorStringFormat</span>} = <span class="i">$Options</span>{<span class="w">vectorstringformat</span>}<span class="sc">;</span>
+ 595 <span class="s">}</span>
+ 596 
+ 597 <span class="c"># Process atom identifier type and related options...</span>
+ 598 <span class="c">#</span>
+<a name="ProcessAtomIdentifierTypeOptions-"></a> 599 <span class="k">sub </span><span class="m">ProcessAtomIdentifierTypeOptions</span> <span class="s">{</span>
+ 600 
+ 601   <span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>} = <span class="i">$Options</span>{<span class="w">atomidentifiertype</span>}<span class="sc">;</span>
+ 602 
+ 603   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">atomidentifiertype</span>} =~ <span class="q">/^AtomicInvariantsAtomTypes$/i</span><span class="s">)</span> <span class="s">{</span>
+ 604     <span class="i">ProcessAtomicInvariantsToUseOption</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 605   <span class="s">}</span>
+ 606   <span class="k">elsif</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">atomidentifiertype</span>} =~ <span class="q">/^FunctionalClassAtomTypes$/i</span><span class="s">)</span> <span class="s">{</span>
+ 607     <span class="i">ProcessFunctionalClassesToUse</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 608   <span class="s">}</span>
+ 609   <span class="k">elsif</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">AtomIdentifierType</span>} =~ <span class="q">/^(DREIDINGAtomTypes|EStateAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 610     <span class="c"># Nothing to do for now...</span>
+ 611   <span class="s">}</span>
+ 612   <span class="k">else</span> <span class="s">{</span>
+ 613     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{atomidentifiertype}, for option \&quot;-a, --AtomIdentifierType\&quot; is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n&quot;</span><span class="sc">;</span>
+ 614   <span class="s">}</span>
+ 615 <span class="s">}</span>
+ 616 
+ 617 <span class="c"># Process specified atomic invariants to use...</span>
+ 618 <span class="c">#</span>
+<a name="ProcessAtomicInvariantsToUseOption-"></a> 619 <span class="k">sub </span><span class="m">ProcessAtomicInvariantsToUseOption</span> <span class="s">{</span>
+ 620   <span class="k">my</span><span class="s">(</span><span class="i">$AtomicInvariant</span><span class="cm">,</span> <span class="i">$AtomSymbolSpecified</span><span class="cm">,</span> <span class="i">@AtomicInvariantsWords</span><span class="s">)</span><span class="sc">;</span>
+ 621 
+ 622   <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">AtomicInvariantsToUse</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 623   <span class="k">if</span> <span class="s">(</span><span class="i">IsEmpty</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">atomicinvariantstouse</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 624     <span class="k">die</span> <span class="q">&quot;Error: Atomic invariants value specified using \&quot;--AtomicInvariantsToUse\&quot; option is empty\n&quot;</span><span class="sc">;</span>
+ 625   <span class="s">}</span>
+ 626   <span class="i">$AtomSymbolSpecified</span> = <span class="n">0</span><span class="sc">;</span>
+ 627   <span class="i">@AtomicInvariantsWords</span> = <span class="k">split</span> <span class="q">/\,/</span><span class="cm">,</span> <span class="i">$Options</span>{<span class="w">atomicinvariantstouse</span>}<span class="sc">;</span>
+ 628   <span class="k">for</span> <span class="i">$AtomicInvariant</span> <span class="s">(</span><span class="i">@AtomicInvariantsWords</span><span class="s">)</span> <span class="s">{</span>
+ 629     <span class="k">if</span> <span class="s">(</span>!<span class="i">AtomTypes::AtomicInvariantsAtomTypes::IsAtomicInvariantAvailable</span><span class="s">(</span><span class="i">$AtomicInvariant</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 630       <span class="k">die</span> <span class="q">&quot;Error: Atomic invariant specified, $AtomicInvariant, using \&quot;--AtomicInvariantsToUse\&quot; option is not valid...\n &quot;</span><span class="sc">;</span>
+ 631     <span class="s">}</span>
+ 632     <span class="k">if</span> <span class="s">(</span><span class="i">$AtomicInvariant</span> =~ <span class="q">/^(AS|AtomSymbol)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 633       <span class="i">$AtomSymbolSpecified</span> = <span class="n">1</span><span class="sc">;</span>
+ 634     <span class="s">}</span>
+ 635     <span class="k">push</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">AtomicInvariantsToUse</span>}}<span class="cm">,</span> <span class="i">$AtomicInvariant</span><span class="sc">;</span>
+ 636   <span class="s">}</span>
+ 637   <span class="k">if</span> <span class="s">(</span>!<span class="i">$AtomSymbolSpecified</span><span class="s">)</span> <span class="s">{</span>
+ 638     <span class="k">die</span> <span class="q">&quot;Error: Atomic invariant, AS or AtomSymbol, must be specified as using \&quot;--AtomicInvariantsToUse\&quot; option...\n &quot;</span><span class="sc">;</span>
+ 639   <span class="s">}</span>
+ 640 <span class="s">}</span>
+ 641 
+ 642 <span class="c"># Process specified functional classes invariants to use...</span>
+ 643 <span class="c">#</span>
+<a name="ProcessFunctionalClassesToUse-"></a> 644 <span class="k">sub </span><span class="m">ProcessFunctionalClassesToUse</span> <span class="s">{</span>
+ 645   <span class="k">my</span><span class="s">(</span><span class="i">$FunctionalClass</span><span class="cm">,</span> <span class="i">@FunctionalClassesToUseWords</span><span class="s">)</span><span class="sc">;</span>
+ 646 
+ 647   <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">FunctionalClassesToUse</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 648   <span class="k">if</span> <span class="s">(</span><span class="i">IsEmpty</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">functionalclassestouse</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 649     <span class="k">die</span> <span class="q">&quot;Error: Functional classes value specified using \&quot;--FunctionalClassesToUse\&quot; option is empty\n&quot;</span><span class="sc">;</span>
+ 650   <span class="s">}</span>
+ 651   <span class="i">@FunctionalClassesToUseWords</span> = <span class="k">split</span> <span class="q">/\,/</span><span class="cm">,</span> <span class="i">$Options</span>{<span class="w">functionalclassestouse</span>}<span class="sc">;</span>
+ 652   <span class="k">for</span> <span class="i">$FunctionalClass</span> <span class="s">(</span><span class="i">@FunctionalClassesToUseWords</span><span class="s">)</span> <span class="s">{</span>
+ 653     <span class="k">if</span> <span class="s">(</span>!<span class="i">AtomTypes::FunctionalClassAtomTypes::IsFunctionalClassAvailable</span><span class="s">(</span><span class="i">$FunctionalClass</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 654       <span class="k">die</span> <span class="q">&quot;Error: Functional class specified, $FunctionalClass, using \&quot;--FunctionalClassesToUse\&quot; option is not valid...\n &quot;</span><span class="sc">;</span>
+ 655     <span class="s">}</span>
+ 656     <span class="k">push</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">FunctionalClassesToUse</span>}}<span class="cm">,</span> <span class="i">$FunctionalClass</span><span class="sc">;</span>
+ 657   <span class="s">}</span>
+ 658 <span class="s">}</span>
+ 659 
+ 660 <span class="c"># Setup script usage  and retrieve command line arguments specified using various options...</span>
+<a name="SetupScriptUsage-"></a> 661 <span class="k">sub </span><span class="m">SetupScriptUsage</span> <span class="s">{</span>
+ 662 
+ 663   <span class="c"># Retrieve all the options...</span>
+ 664   <span class="i">%Options</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 665 
+ 666   <span class="i">$Options</span>{<span class="w">aromaticitymodel</span>} = <span class="q">&#39;MayaChemToolsAromaticityModel&#39;</span><span class="sc">;</span>
+ 667 
+ 668   <span class="i">$Options</span>{<span class="w">atomidentifiertype</span>} = <span class="q">&#39;AtomicInvariantsAtomTypes&#39;</span><span class="sc">;</span>
+ 669   <span class="i">$Options</span>{<span class="w">atomicinvariantstouse</span>} = <span class="q">&#39;AS,X,BO,H,FC&#39;</span><span class="sc">;</span>
+ 670 
+ 671   <span class="i">$Options</span>{<span class="w">functionalclassestouse</span>} = <span class="q">&#39;HBD,HBA,PI,NI,Ar,Hal&#39;</span><span class="sc">;</span>
+ 672 
+ 673   <span class="i">$Options</span>{<span class="w">compoundidmode</span>} = <span class="q">&#39;LabelPrefix&#39;</span><span class="sc">;</span>
+ 674   <span class="i">$Options</span>{<span class="w">compoundidlabel</span>} = <span class="q">&#39;CompoundID&#39;</span><span class="sc">;</span>
+ 675   <span class="i">$Options</span>{<span class="w">datafieldsmode</span>} = <span class="q">&#39;CompoundID&#39;</span><span class="sc">;</span>
+ 676 
+ 677   <span class="i">$Options</span>{<span class="w">filter</span>} = <span class="q">&#39;Yes&#39;</span><span class="sc">;</span>
+ 678 
+ 679   <span class="i">$Options</span>{<span class="w">keeplargestcomponent</span>} = <span class="q">&#39;Yes&#39;</span><span class="sc">;</span>
+ 680 
+ 681   <span class="i">$Options</span>{<span class="w">mindistance</span>} = <span class="n">1</span><span class="sc">;</span>
+ 682   <span class="i">$Options</span>{<span class="w">maxdistance</span>} = <span class="n">10</span><span class="sc">;</span>
+ 683 
+ 684   <span class="i">$Options</span>{<span class="w">output</span>} = <span class="q">&#39;text&#39;</span><span class="sc">;</span>
+ 685   <span class="i">$Options</span>{<span class="w">outdelim</span>} = <span class="q">&#39;comma&#39;</span><span class="sc">;</span>
+ 686   <span class="i">$Options</span>{<span class="w">quote</span>} = <span class="q">&#39;yes&#39;</span><span class="sc">;</span>
+ 687 
+ 688   <span class="i">$Options</span>{<span class="w">vectorstringformat</span>} = <span class="q">&#39;IDsAndValuesString&#39;</span><span class="sc">;</span>
+ 689 
+ 690   <span class="k">if</span> <span class="s">(</span>!<span class="i">GetOptions</span><span class="s">(</span>\<span class="i">%Options</span><span class="cm">,</span> <span class="q">&quot;aromaticitymodel=s&quot;</span><span class="cm">,</span> <span class="q">&quot;atomidentifiertype|a=s&quot;</span><span class="cm">,</span> <span class="q">&quot;atomicinvariantstouse=s&quot;</span><span class="cm">,</span> <span class="q">&quot;functionalclassestouse=s&quot;</span><span class="cm">,</span> <span class="q">&quot;compoundid=s&quot;</span><span class="cm">,</span> <span class="q">&quot;compoundidlabel=s&quot;</span><span class="cm">,</span> <span class="q">&quot;compoundidmode=s&quot;</span><span class="cm">,</span> <span class="q">&quot;datafields=s&quot;</span><span class="cm">,</span> <span class="q">&quot;datafieldsmode|d=s&quot;</span><span class="cm">,</span> <span class="q">&quot;filter|f=s&quot;</span><span class="cm">,</span> <span class="q">&quot;fingerprintslabel=s&quot;</span><span class="cm">,</span>  <span class="q">&quot;help|h&quot;</span><span class="cm">,</span> <span class="q">&quot;keeplargestcomponent|k=s&quot;</span><span class="cm">,</span>  <span class="q">&quot;mindistance=s&quot;</span><span class="cm">,</span> <span class="q">&quot;maxdistance=s&quot;</span><span class="cm">,</span> <span class="q">&quot;outdelim=s&quot;</span><span class="cm">,</span> <span class="q">&quot;output=s&quot;</span><span class="cm">,</span> <span class="q">&quot;overwrite|o&quot;</span><span class="cm">,</span> <span class="q">&quot;quote|q=s&quot;</span><span class="cm">,</span> <span class="q">&quot;root|r=s&quot;</span><span class="cm">,</span> <span class="q">&quot;vectorstringformat|v=s&quot;</span><span class="cm">,</span> <span class="q">&quot;workingdir|w=s&quot;</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 691     <span class="k">die</span> <span class="q">&quot;\nTo get a list of valid options and their values, use \&quot;$ScriptName -h\&quot; or\n\&quot;perl -S $ScriptName -h\&quot; command and try again...\n&quot;</span><span class="sc">;</span>
+ 692   <span class="s">}</span>
+ 693   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">workingdir</span>}<span class="s">)</span> <span class="s">{</span>
+ 694     <span class="k">if</span> <span class="s">(</span>! <span class="k">-d</span> <span class="i">$Options</span>{<span class="w">workingdir</span>}<span class="s">)</span> <span class="s">{</span>
+ 695       <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{workingdir}, for option \&quot;-w --workingdir\&quot; is not a directory name.\n&quot;</span><span class="sc">;</span>
+ 696     <span class="s">}</span>
+ 697     <span class="k">chdir</span> <span class="i">$Options</span>{<span class="w">workingdir</span>} <span class="k">or</span> <span class="k">die</span> <span class="q">&quot;Error: Couldn&#39;t chdir $Options{workingdir}: $! \n&quot;</span><span class="sc">;</span>
+ 698   <span class="s">}</span>
+ 699   <span class="k">if</span> <span class="s">(</span>!<span class="i">Molecule::IsSupportedAromaticityModel</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">aromaticitymodel</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 700     <span class="k">my</span><span class="s">(</span><span class="i">@SupportedModels</span><span class="s">)</span> = <span class="i">Molecule::GetSupportedAromaticityModels</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span>
+ 701     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{aromaticitymodel}, for option \&quot;--AromaticityModel\&quot; is not valid. Supported aromaticity models in current release of MayaChemTools: @SupportedModels\n&quot;</span><span class="sc">;</span>
+ 702   <span class="s">}</span>
+ 703   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">atomidentifiertype</span>} !~ <span class="q">/^(AtomicInvariantsAtomTypes|DREIDINGAtomTypes|EStateAtomTypes|FunctionalClassAtomTypes|MMFF94AtomTypes|SLogPAtomTypes|SYBYLAtomTypes|TPSAAtomTypes|UFFAtomTypes)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 704     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{atomidentifiertype}, for option \&quot;-a, --AtomIdentifierType\&quot; is not valid. Supported atom identifier types in current release of MayaChemTools: AtomicInvariantsAtomTypes, DREIDINGAtomTypes, EStateAtomTypes, FunctionalClassAtomTypes, MMFF94AtomTypes, SLogPAtomTypes, SYBYLAtomTypes, TPSAAtomTypes, UFFAtomTypes\n&quot;</span><span class="sc">;</span>
+ 705   <span class="s">}</span>
+ 706   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">compoundidmode</span>} !~ <span class="q">/^(DataField|MolName|LabelPrefix|MolNameOrLabelPrefix)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 707     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{compoundidmode}, for option \&quot;--CompoundIDMode\&quot; is not valid. Allowed values: DataField, MolName, LabelPrefix or MolNameOrLabelPrefix\n&quot;</span><span class="sc">;</span>
+ 708   <span class="s">}</span>
+ 709   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">datafieldsmode</span>} !~ <span class="q">/^(All|Common|Specify|CompoundID)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 710     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{datafieldsmode}, for option \&quot;-d, --DataFieldsMode\&quot; is not valid. Allowed values: All, Common, Specify or CompoundID\n&quot;</span><span class="sc">;</span>
+ 711   <span class="s">}</span>
+ 712   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">filter</span>} !~ <span class="q">/^(Yes|No)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 713     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{filter}, for option \&quot;-f, --Filter\&quot; is not valid. Allowed values: Yes or No\n&quot;</span><span class="sc">;</span>
+ 714   <span class="s">}</span>
+ 715   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">keeplargestcomponent</span>} !~ <span class="q">/^(Yes|No)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 716     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{keeplargestcomponent}, for option \&quot;-k, --KeepLargestComponent\&quot; is not valid. Allowed values: Yes or No\n&quot;</span><span class="sc">;</span>
+ 717   <span class="s">}</span>
+ 718   <span class="k">if</span> <span class="s">(</span>!<span class="i">IsPositiveInteger</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">mindistance</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 719     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{mindistance}, for option \&quot;--MinDistance\&quot; is not valid. Allowed values: &gt; 0 \n&quot;</span><span class="sc">;</span>
+ 720   <span class="s">}</span>
+ 721   <span class="k">if</span> <span class="s">(</span>!<span class="i">IsPositiveInteger</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">maxdistance</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span>
+ 722     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{maxdistance}, for option \&quot;--MaxDistance\&quot; is not valid. Allowed values: &gt; 0 \n&quot;</span><span class="sc">;</span>
+ 723   <span class="s">}</span>
+ 724   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">mindistance</span>} &gt; <span class="i">$Options</span>{<span class="w">maxdistance</span>}<span class="s">)</span> <span class="s">{</span>
+ 725     <span class="k">die</span> <span class="q">&quot;Error: The value specified, specified, $Options{mindistance}, for option \&quot;--MinDistance\&quot; must be less than the value specified, $Options{maxdistance}, for option \&quot;--MaxDistance\&quot; \n&quot;</span><span class="sc">;</span>
+ 726   <span class="s">}</span>
+ 727   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">output</span>} !~ <span class="q">/^(SD|FP|text|all)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 728     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{output}, for option \&quot;--output\&quot; is not valid. Allowed values: SD, FP, text, or all\n&quot;</span><span class="sc">;</span>
+ 729   <span class="s">}</span>
+ 730   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">outdelim</span>} !~ <span class="q">/^(comma|semicolon|tab)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 731     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{outdelim}, for option \&quot;--outdelim\&quot; is not valid. Allowed values: comma, tab, or semicolon\n&quot;</span><span class="sc">;</span>
+ 732   <span class="s">}</span>
+ 733   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">quote</span>} !~ <span class="q">/^(Yes|No)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 734     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{quote}, for option \&quot;-q --quote\&quot; is not valid. Allowed values: Yes or No\n&quot;</span><span class="sc">;</span>
+ 735   <span class="s">}</span>
+ 736   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">outdelim</span>} =~ <span class="q">/semicolon/i</span> &amp;&amp; <span class="i">$Options</span>{<span class="w">quote</span>} =~ <span class="q">/^No$/i</span><span class="s">)</span> <span class="s">{</span>
+ 737     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{quote}, for option \&quot;-q --quote\&quot; is not allowed with, semicolon value of \&quot;--outdelim\&quot; option: Fingerprints string use semicolon as delimiter for various data fields and must be quoted.\n&quot;</span><span class="sc">;</span>
+ 738   <span class="s">}</span>
+ 739   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">vectorstringformat</span>} !~ <span class="q">/^(IDsAndValuesString|IDsAndValuesPairsString|ValuesAndIDsString|ValuesAndIDsPairsString)$/i</span><span class="s">)</span> <span class="s">{</span>
+ 740     <span class="k">die</span> <span class="q">&quot;Error: The value specified, $Options{vectorstringformat}, for option \&quot;-v, --VectorStringFormat\&quot; is not valid. Allowed values: IDsAndValuesString, IDsAndValuesPairsString, ValuesAndIDsString or ValuesAndIDsPairsString\n&quot;</span><span class="sc">;</span>
+ 741   <span class="s">}</span>
+ 742 <span class="s">}</span>
+ 743 
+<a name="EOF-"></a></pre>
+<p>&nbsp;</p>
+<br />
+<center>
+<img src="../../../images/h2o2.png">
+</center>
+</body>
+</html>