| 0 | 1 <html> | 
|  | 2 <head> | 
|  | 3 <title>MayaChemTools:Code:InfoPDBFiles.pl</title> | 
|  | 4 <meta http-equiv="content-type" content="text/html;charset=utf-8"> | 
|  | 5 <link rel="stylesheet" type="text/css" href="../../../css/MayaChemToolsCode.css"> | 
|  | 6 </head> | 
|  | 7 <body leftmargin="20" rightmargin="20" topmargin="10" bottommargin="10"> | 
|  | 8 <br/> | 
|  | 9 <center> | 
|  | 10 <a href="http://www.mayachemtools.org" title="MayaChemTools Home"><img src="../../../images/MayaChemToolsLogo.gif" border="0" alt="MayaChemTools"></a> | 
|  | 11 </center> | 
|  | 12 <br/> | 
|  | 13 <pre> | 
|  | 14    1 #!/usr/bin/perl -w | 
|  | 15    2 <span class="c">#</span> | 
|  | 16    3 <span class="c"># $RCSfile: InfoPDBFiles.pl,v $</span> | 
|  | 17    4 <span class="c"># $Date: 2015/02/28 20:46:20 $</span> | 
|  | 18    5 <span class="c"># $Revision: 1.36 $</span> | 
|  | 19    6 <span class="c">#</span> | 
|  | 20    7 <span class="c"># Author: Manish Sud <msud@san.rr.com></span> | 
|  | 21    8 <span class="c">#</span> | 
|  | 22    9 <span class="c"># Copyright (C) 2015 Manish Sud. All rights reserved.</span> | 
|  | 23   10 <span class="c">#</span> | 
|  | 24   11 <span class="c"># This file is part of MayaChemTools.</span> | 
|  | 25   12 <span class="c">#</span> | 
|  | 26   13 <span class="c"># MayaChemTools is free software; you can redistribute it and/or modify it under</span> | 
|  | 27   14 <span class="c"># the terms of the GNU Lesser General Public License as published by the Free</span> | 
|  | 28   15 <span class="c"># Software Foundation; either version 3 of the License, or (at your option) any</span> | 
|  | 29   16 <span class="c"># later version.</span> | 
|  | 30   17 <span class="c">#</span> | 
|  | 31   18 <span class="c"># MayaChemTools is distributed in the hope that it will be useful, but without</span> | 
|  | 32   19 <span class="c"># any warranty; without even the implied warranty of merchantability of fitness</span> | 
|  | 33   20 <span class="c"># for a particular purpose.  See the GNU Lesser General Public License for more</span> | 
|  | 34   21 <span class="c"># details.</span> | 
|  | 35   22 <span class="c">#</span> | 
|  | 36   23 <span class="c"># You should have received a copy of the GNU Lesser General Public License</span> | 
|  | 37   24 <span class="c"># along with MayaChemTools; if not, see <http://www.gnu.org/licenses/> or</span> | 
|  | 38   25 <span class="c"># write to the Free Software Foundation Inc., 59 Temple Place, Suite 330,</span> | 
|  | 39   26 <span class="c"># Boston, MA, 02111-1307, USA.</span> | 
|  | 40   27 <span class="c">#</span> | 
|  | 41   28 | 
|  | 42   29 <span class="k">use</span> <span class="w">strict</span><span class="sc">;</span> | 
|  | 43   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">"$FindBin::Bin/../lib"</span><span class="sc">;</span> | 
|  | 44   31 <span class="k">use</span> <span class="w">Getopt::Long</span><span class="sc">;</span> | 
|  | 45   32 <span class="k">use</span> <span class="w">File::Basename</span><span class="sc">;</span> | 
|  | 46   33 <span class="k">use</span> <span class="w">Text::ParseWords</span><span class="sc">;</span> | 
|  | 47   34 <span class="k">use</span> <span class="w">Benchmark</span><span class="sc">;</span> | 
|  | 48   35 <span class="k">use</span> <span class="w">FileUtil</span><span class="sc">;</span> | 
|  | 49   36 <span class="k">use</span> <span class="w">TextUtil</span><span class="sc">;</span> | 
|  | 50   37 <span class="k">use</span> <span class="w">PDBFileUtil</span><span class="sc">;</span> | 
|  | 51   38 | 
|  | 52   39 <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> | 
|  | 53   40 | 
|  | 54   41 <span class="c"># Autoflush STDOUT</span> | 
|  | 55   42 <span class="i">$|</span> = <span class="n">1</span><span class="sc">;</span> | 
|  | 56   43 | 
|  | 57   44 <span class="c"># Starting message...</span> | 
|  | 58   45 <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> | 
|  | 59   46 <span class="k">print</span> <span class="q">"\n$ScriptName: Starting...\n\n"</span><span class="sc">;</span> | 
|  | 60   47 <span class="i">$StartTime</span> = <span class="w">new</span> <span class="w">Benchmark</span><span class="sc">;</span> | 
|  | 61   48 | 
|  | 62   49 <span class="c"># Get the options and setup script...</span> | 
|  | 63   50 <span class="i">SetupScriptUsage</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 64   51 <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">help</span>} || <span class="i">@ARGV</span> < <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 65   52   <span class="k">die</span> <span class="i">GetUsageFromPod</span><span class="s">(</span><span class="q">"$FindBin::Bin/$ScriptName"</span><span class="s">)</span><span class="sc">;</span> | 
|  | 66   53 <span class="s">}</span> | 
|  | 67   54 | 
|  | 68   55 <span class="k">my</span><span class="s">(</span><span class="i">@PDBFilesList</span><span class="s">)</span><span class="sc">;</span> | 
|  | 69   56 <span class="i">@PDBFilesList</span> = <span class="i">ExpandFileNames</span><span class="s">(</span>\<span class="i">@ARGV</span><span class="cm">,</span> <span class="q">"pdb"</span><span class="s">)</span><span class="sc">;</span> | 
|  | 70   57 | 
|  | 71   58 <span class="c"># Process options...</span> | 
|  | 72   59 <span class="k">print</span> <span class="q">"Processing options...\n"</span><span class="sc">;</span> | 
|  | 73   60 <span class="k">my</span><span class="s">(</span><span class="i">%OptionsInfo</span><span class="s">)</span><span class="sc">;</span> | 
|  | 74   61 <span class="i">ProcessOptions</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 75   62 | 
|  | 76   63 <span class="c"># Setup information about input files...</span> | 
|  | 77   64 <span class="k">my</span><span class="s">(</span><span class="i">%PDBFilesInfo</span><span class="s">)</span><span class="sc">;</span> | 
|  | 78   65 <span class="k">print</span> <span class="q">"Checking input PDB file(s)...\n"</span><span class="sc">;</span> | 
|  | 79   66 <span class="i">RetrievePDBFilesInfo</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 80   67 | 
|  | 81   68 <span class="c"># Process input files..</span> | 
|  | 82   69 <span class="k">my</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span> | 
|  | 83   70 <span class="k">if</span> <span class="s">(</span><span class="i">@PDBFilesList</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 84   71   <span class="k">print</span> <span class="q">"\nProcessing PDB files...\n"</span><span class="sc">;</span> | 
|  | 85   72 <span class="s">}</span> | 
|  | 86   73 <span class="k">for</span> <span class="i">$FileIndex</span> <span class="s">(</span><span class="n">0</span> .. <span class="i">$#PDBFilesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 87   74   <span class="k">if</span> <span class="s">(</span><span class="i">$PDBFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$FileIndex</span>]<span class="s">)</span> <span class="s">{</span> | 
|  | 88   75     <span class="k">print</span> <span class="q">"\nProcessing file $PDBFilesList[$FileIndex]...\n"</span><span class="sc">;</span> | 
|  | 89   76     <span class="i">ListPDBFileInfo</span><span class="s">(</span><span class="i">$FileIndex</span><span class="s">)</span><span class="sc">;</span> | 
|  | 90   77   <span class="s">}</span> | 
|  | 91   78 <span class="s">}</span> | 
|  | 92   79 <span class="i">ListTotalSizeOfFiles</span><span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 93   80 | 
|  | 94   81 <span class="k">print</span> <span class="q">"\n$ScriptName:Done...\n\n"</span><span class="sc">;</span> | 
|  | 95   82 | 
|  | 96   83 <span class="i">$EndTime</span> = <span class="w">new</span> <span class="w">Benchmark</span><span class="sc">;</span> | 
|  | 97   84 <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> | 
|  | 98   85 <span class="k">print</span> <span class="q">"Total time: "</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">"\n"</span><span class="sc">;</span> | 
|  | 99   86 | 
|  | 100   87 <span class="c">###############################################################################</span> | 
|  | 101   88 | 
|  | 102   89 <span class="c"># List appropriate information...</span> | 
|  | 103 <a name="ListPDBFileInfo-"></a>  90 <span class="k">sub </span><span class="m">ListPDBFileInfo</span> <span class="s">{</span> | 
|  | 104   91   <span class="k">my</span><span class="s">(</span><span class="i">$Index</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 105   92   <span class="k">my</span><span class="s">(</span><span class="i">$PDBFile</span><span class="cm">,</span> <span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 106   93 | 
|  | 107   94   <span class="i">$PDBFile</span> = <span class="i">$PDBFilesList</span>[<span class="i">$Index</span>]<span class="sc">;</span> | 
|  | 108   95   <span class="i">$PDBRecordLinesRef</span> = <span class="i">ReadPDBFile</span><span class="s">(</span><span class="i">$PDBFile</span><span class="s">)</span><span class="sc">;</span> | 
|  | 109   96 | 
|  | 110   97   <span class="c"># Header informaton...</span> | 
|  | 111   98   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">ListHeaderInfo</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 112   99     <span class="i">ListHeaderInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 113  100   <span class="s">}</span> | 
|  | 114  101 | 
|  | 115  102   <span class="c"># Experiment informaton...</span> | 
|  | 116  103   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">ListExperimentalTechniqueInfo</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 117  104     <span class="i">ListExperimentalTechniqueInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 118  105   <span class="s">}</span> | 
|  | 119  106 | 
|  | 120  107   <span class="c"># Total number of records...</span> | 
|  | 121  108   <span class="k">my</span><span class="s">(</span><span class="i">$TotalRecordsCount</span><span class="s">)</span> = <span class="k">scalar</span> <span class="i">@</span>{<span class="i">$PDBRecordLinesRef</span>}<span class="sc">;</span> | 
|  | 122  109   <span class="k">print</span> <span class="q">"\nTotal number of records: $TotalRecordsCount\n"</span><span class="sc">;</span> | 
|  | 123  110 | 
|  | 124  111   <span class="c"># List record type count information...</span> | 
|  | 125  112   <span class="i">ListRecordTypesInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 126  113 | 
|  | 127  114   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountChains</span>} || <span class="i">$OptionsInfo</span>{<span class="w">CountResiduesInChains</span>} || <span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyInChains</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 128  115     <span class="i">ListChainsAndResiduesInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 129  116   <span class="s">}</span> | 
|  | 130  117   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountResiduesAll</span>} || <span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyAll</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 131  118     <span class="i">ListAllResiduesInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 132  119   <span class="s">}</span> | 
|  | 133  120   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">ResidueNumbersInfo</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 134  121     <span class="i">ListResidueNumbersInfo</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 135  122   <span class="s">}</span> | 
|  | 136  123   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CalculateBoundingBox</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 137  124     <span class="i">ListBoundingBox</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 138  125   <span class="s">}</span> | 
|  | 139  126 | 
|  | 140  127   <span class="c"># File size and modification information...</span> | 
|  | 141  128   <span class="k">print</span> <span class="q">"\nFile size: "</span><span class="cm">,</span> <span class="i">FormatFileSize</span><span class="s">(</span><span class="i">$PDBFilesInfo</span>{<span class="w">FileSize</span>}[<span class="i">$Index</span>]<span class="s">)</span><span class="cm">,</span> <span class="q">" \n"</span><span class="sc">;</span> | 
|  | 142  129   <span class="k">print</span> <span class="q">"Last modified: "</span><span class="cm">,</span> <span class="i">$PDBFilesInfo</span>{<span class="w">FileLastModified</span>}[<span class="i">$Index</span>]<span class="cm">,</span> <span class="q">" \n"</span><span class="sc">;</span> | 
|  | 143  130 <span class="s">}</span> | 
|  | 144  131 | 
|  | 145 <a name="ListHeaderInfo-"></a> 132 <span class="k">sub </span><span class="m">ListHeaderInfo</span> <span class="s">{</span> | 
|  | 146  133   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 147  134   <span class="k">my</span><span class="s">(</span><span class="i">$HeaderRecordLine</span><span class="cm">,</span> <span class="i">$Classification</span><span class="cm">,</span> <span class="i">$DepositionDate</span><span class="cm">,</span> <span class="i">$IDCode</span><span class="s">)</span><span class="sc">;</span> | 
|  | 148  135 | 
|  | 149  136   <span class="s">(</span><span class="i">$Classification</span><span class="cm">,</span> <span class="i">$DepositionDate</span><span class="cm">,</span> <span class="i">$IDCode</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> | 
|  | 150  137   <span class="i">$HeaderRecordLine</span> = <span class="i">$PDBRecordLinesRef</span>->[<span class="n">0</span>]<span class="sc">;</span> | 
|  | 151  138   <span class="k">if</span> <span class="s">(</span><span class="i">IsHeaderRecordType</span><span class="s">(</span><span class="i">$HeaderRecordLine</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 152  139     <span class="s">(</span><span class="i">$Classification</span><span class="cm">,</span> <span class="i">$DepositionDate</span><span class="cm">,</span> <span class="i">$IDCode</span><span class="s">)</span> = <span class="i">ParseHeaderRecordLine</span><span class="s">(</span><span class="i">$HeaderRecordLine</span><span class="s">)</span><span class="sc">;</span> | 
|  | 153  140   <span class="s">}</span> | 
|  | 154  141 | 
|  | 155  142   <span class="i">$Classification</span> = <span class="i">IsEmpty</span><span class="s">(</span><span class="i">$Classification</span><span class="s">)</span> ? <span class="q">'Not available'</span> <span class="co">:</span> <span class="i">$Classification</span><span class="sc">;</span> | 
|  | 156  143   <span class="i">$DepositionDate</span> = <span class="i">IsEmpty</span><span class="s">(</span><span class="i">$DepositionDate</span><span class="s">)</span> ? <span class="q">'Not available'</span> <span class="co">:</span> <span class="i">$DepositionDate</span><span class="sc">;</span> | 
|  | 157  144   <span class="i">$IDCode</span> = <span class="i">IsEmpty</span><span class="s">(</span><span class="i">$IDCode</span><span class="s">)</span> ? <span class="q">'Not available'</span> <span class="co">:</span> <span class="i">$IDCode</span><span class="sc">;</span> | 
|  | 158  145 | 
|  | 159  146   <span class="k">print</span> <span class="q">"\nClassification: $Classification\nID: $IDCode\nDeposition date: $DepositionDate\n"</span><span class="sc">;</span> | 
|  | 160  147 <span class="s">}</span> | 
|  | 161  148 | 
|  | 162  149 <span class="c"># List experimental technique information info...</span> | 
|  | 163 <a name="ListExperimentalTechniqueInfo-"></a> 150 <span class="k">sub </span><span class="m">ListExperimentalTechniqueInfo</span> <span class="s">{</span> | 
|  | 164  151   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 165  152   <span class="k">my</span><span class="s">(</span><span class="i">$ExperimentalTechnique</span><span class="cm">,</span> <span class="i">$Resolution</span><span class="cm">,</span> <span class="i">$ResolutionUnits</span><span class="s">)</span><span class="sc">;</span> | 
|  | 166  153 | 
|  | 167  154   <span class="i">$ExperimentalTechnique</span> = <span class="i">GetExperimentalTechnique</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 168  155   <span class="k">print</span> <span class="q">"\nExperimental technique: "</span> . <span class="s">(</span><span class="i">$ExperimentalTechnique</span> ? <span class="i">$ExperimentalTechnique</span> <span class="co">:</span> <span class="q">"Not available"</span><span class="s">)</span> . <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 169  156 | 
|  | 170  157   <span class="s">(</span><span class="i">$Resolution</span><span class="cm">,</span> <span class="i">$ResolutionUnits</span><span class="s">)</span> = <span class="i">GetExperimentalTechniqueResolution</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 171  158   <span class="k">print</span> <span class="q">"Resolution: "</span> . <span class="s">(</span><span class="i">$Resolution</span> ? <span class="q">"$Resolution $ResolutionUnits"</span> <span class="co">:</span> <span class="q">"Not available"</span><span class="s">)</span> . <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 172  159 | 
|  | 173  160 <span class="s">}</span> | 
|  | 174  161 | 
|  | 175  162 <span class="c"># List record type info...</span> | 
|  | 176 <a name="ListRecordTypesInfo-"></a> 163 <span class="k">sub </span><span class="m">ListRecordTypesInfo</span> <span class="s">{</span> | 
|  | 177  164   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 178  165   <span class="k">my</span><span class="s">(</span><span class="i">$RecordType</span><span class="cm">,</span> <span class="i">$RecordCount</span><span class="cm">,</span> <span class="i">$RecordTypesCountRef</span><span class="cm">,</span> <span class="i">@RecordTypeCountInfo</span><span class="s">)</span><span class="sc">;</span> | 
|  | 179  166 | 
|  | 180  167   <span class="i">$RecordTypesCountRef</span> = <span class="i">GetRecordTypesCount</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 181  168 | 
|  | 182  169   <span class="i">@RecordTypeCountInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 183  170   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountRecordType</span>} =~ <span class="q">/^All$/i</span><span class="s">)</span> <span class="s">{</span> | 
|  | 184  171     <span class="k">for</span> <span class="i">$RecordType</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$RecordTypesCountRef</span>->{<span class="w">RecordTypes</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 185  172       <span class="i">$RecordCount</span> = <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordType</span>}<span class="sc">;</span> | 
|  | 186  173       <span class="k">push</span> <span class="i">@RecordTypeCountInfo</span><span class="cm">,</span> <span class="q">"$RecordType - $RecordCount"</span><span class="sc">;</span> | 
|  | 187  174     <span class="s">}</span> | 
|  | 188  175   <span class="s">}</span> | 
|  | 189  176   <span class="k">else</span> <span class="s">{</span> | 
|  | 190  177     <span class="k">for</span> <span class="i">$RecordType</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedRecordTypes</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 191  178       <span class="i">$RecordCount</span> = <span class="s">(</span><span class="k">exists</span> <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordType</span>}<span class="s">)</span> ? <span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordType</span>}<span class="s">)</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 192  179       <span class="k">push</span> <span class="i">@RecordTypeCountInfo</span><span class="cm">,</span> <span class="q">"$RecordType - $RecordCount"</span><span class="sc">;</span> | 
|  | 193  180     <span class="s">}</span> | 
|  | 194  181   <span class="s">}</span> | 
|  | 195  182   <span class="k">print</span> <span class="q">"Number of individual records: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@RecordTypeCountInfo</span><span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 196  183 | 
|  | 197  184   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CheckMasterRecord</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 198  185     <span class="i">CheckMasterRecord</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span><span class="cm">,</span> <span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 199  186   <span class="s">}</span> | 
|  | 200  187 <span class="s">}</span> | 
|  | 201  188 | 
|  | 202  189 <span class="c"># List information about residues and chains...</span> | 
|  | 203 <a name="ListChainsAndResiduesInfo-"></a> 190 <span class="k">sub </span><span class="m">ListChainsAndResiduesInfo</span> <span class="s">{</span> | 
|  | 204  191   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 205  192   <span class="k">my</span><span class="s">(</span><span class="i">$ResidueName</span><span class="cm">,</span> <span class="i">$ResidueCount</span><span class="cm">,</span> <span class="i">$ChainCount</span><span class="cm">,</span> <span class="i">$ChainID</span><span class="cm">,</span> <span class="i">$CollectChainResiduesBeyondTER</span><span class="cm">,</span> <span class="i">$ChainsAndResiduesInfoRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 206  193 | 
|  | 207  194   <span class="i">$CollectChainResiduesBeyondTER</span> = <span class="n">1</span><span class="sc">;</span> | 
|  | 208  195   <span class="i">$ChainsAndResiduesInfoRef</span> = <span class="i">GetChainsAndResidues</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="cm">,</span> <span class="q">'AtomAndHetatm'</span><span class="cm">,</span> <span class="i">$CollectChainResiduesBeyondTER</span><span class="s">)</span><span class="sc">;</span> | 
|  | 209  196   <span class="i">$ChainCount</span> = <span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="sc">;</span> | 
|  | 210  197   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountChains</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 211  198     <span class="k">print</span> <span class="q">"\nNumber of chains: $ChainCount \n"</span><span class="sc">;</span> | 
|  | 212  199     <span class="k">my</span><span class="s">(</span><span class="i">$ChainID</span><span class="cm">,</span> <span class="i">@ChainIDsList</span><span class="s">)</span><span class="sc">;</span> | 
|  | 213  200     <span class="i">@ChainIDsList</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 214  201     <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 215  202       <span class="k">push</span> <span class="i">@ChainIDsList</span><span class="cm">,</span> <span class="i">CleanupChainID</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span><span class="sc">;</span> | 
|  | 216  203     <span class="s">}</span> | 
|  | 217  204     <span class="k">print</span> <span class="q">"Chain IDs: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@ChainIDsList</span><span class="cm">,</span> <span class="q">', '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span><span class="q">"\n"</span><span class="sc">;</span> | 
|  | 218  205   <span class="s">}</span> | 
|  | 219  206 | 
|  | 220  207   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountResiduesInChains</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 221  208     <span class="k">my</span><span class="s">(</span><span class="i">$TotalResiduesCount</span><span class="cm">,</span> <span class="i">$ResidueCountInfo</span><span class="cm">,</span> <span class="i">@ResiduesCountInfo</span><span class="s">)</span><span class="sc">;</span> | 
|  | 222  209     <span class="i">@ResiduesCountInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 223  210     <span class="i">$TotalResiduesCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 224  211     <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 225  212       <span class="i">$ResidueCount</span> = <span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">Residues</span>}{<span class="i">$ChainID</span>}}<span class="sc">;</span> | 
|  | 226  213       <span class="i">$TotalResiduesCount</span> += <span class="i">$ResidueCount</span><span class="sc">;</span> | 
|  | 227  214       <span class="i">$ResidueCountInfo</span> =  <span class="q">"Chain ${ChainID} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 228  215       <span class="k">push</span> <span class="i">@ResiduesCountInfo</span><span class="cm">,</span> <span class="i">$ResidueCountInfo</span><span class="sc">;</span> | 
|  | 229  216     <span class="s">}</span> | 
|  | 230  217     <span class="k">print</span> <span class="q">"\nNumber of residues in chain(s): "</span><span class="sc">;</span> | 
|  | 231  218     <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 232  219       <span class="k">print</span> <span class="q">"Total - $TotalResiduesCount; "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@ResiduesCountInfo</span><span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span><span class="q">"\n"</span><span class="sc">;</span> | 
|  | 233  220     <span class="s">}</span> | 
|  | 234  221     <span class="k">else</span> <span class="s">{</span> | 
|  | 235  222       <span class="k">print</span> <span class="q">"$TotalResiduesCount\n"</span><span class="sc">;</span> | 
|  | 236  223     <span class="s">}</span> | 
|  | 237  224 | 
|  | 238  225     <span class="c"># List of residues in each chain...</span> | 
|  | 239  226     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">3</span><span class="s">)</span> <span class="s">{</span> | 
|  | 240  227       <span class="k">print</span> <span class="q">"List of residues in chain(s): \n"</span><span class="sc">;</span> | 
|  | 241  228       <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 242  229         <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 243  230           <span class="k">print</span> <span class="q">"Chain "</span><span class="cm">,</span> <span class="i">CleanupChainID</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span><span class="cm">,</span> <span class="q">": "</span><span class="sc">;</span> | 
|  | 244  231         <span class="s">}</span> | 
|  | 245  232         <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">Residues</span>}{<span class="i">$ChainID</span>}}<span class="cm">,</span> <span class="q">', '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span><span class="q">"\n"</span><span class="sc">;</span> | 
|  | 246  233       <span class="s">}</span> | 
|  | 247  234     <span class="s">}</span> | 
|  | 248  235   <span class="s">}</span> | 
|  | 249  236   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyInChains</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 250  237     <span class="c"># Setup a hash using residue count as key for sorting the values...</span> | 
|  | 251  238     <span class="k">my</span><span class="s">(</span><span class="i">%ResiduesCountToNameMap</span><span class="s">)</span><span class="sc">;</span> | 
|  | 252  239     <span class="i">%ResiduesCountToNameMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 253  240     <span class="i">@</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ChainIDs</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 254  241     <span class="i">%</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 255  242 | 
|  | 256  243     <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 257  244       <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ChainIDs</span>}}<span class="cm">,</span> <span class="i">$ChainID</span><span class="sc">;</span> | 
|  | 258  245       <span class="i">%</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 259  246 | 
|  | 260  247       <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">sort</span> <span class="k">keys</span> <span class="i">%</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}{<span class="i">$ChainID</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 261  248         <span class="i">$ResidueCount</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 262  249         <span class="c"># Setup count value for each chain...</span> | 
|  | 263  250         <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueCount</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 264  251           <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueCount</span>} .= <span class="q">"~${ResidueName}"</span><span class="sc">;</span> | 
|  | 265  252         <span class="s">}</span> | 
|  | 266  253         <span class="k">else</span> <span class="s">{</span> | 
|  | 267  254           <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueCount</span>} = <span class="i">$ResidueName</span><span class="sc">;</span> | 
|  | 268  255         <span class="s">}</span> | 
|  | 269  256       <span class="s">}</span> | 
|  | 270  257     <span class="s">}</span> | 
|  | 271  258     <span class="c"># Collect data for all the residues in all the chains...</span> | 
|  | 272  259     <span class="k">my</span><span class="s">(</span><span class="i">%AllResiduesNameToCountMap</span><span class="cm">,</span> <span class="i">%AllResiduesCountToNameMap</span><span class="s">)</span><span class="sc">;</span> | 
|  | 273  260     <span class="i">%AllResiduesNameToCountMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 274  261     <span class="i">%AllResiduesCountToNameMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 275  262     <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 276  263       <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 277  264         <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}{<span class="i">$ChainID</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 278  265           <span class="i">$ResidueCount</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 279  266           <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$AllResiduesNameToCountMap</span>{<span class="i">$ResidueName</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 280  267             <span class="i">$AllResiduesNameToCountMap</span>{<span class="i">$ResidueName</span>} += <span class="i">$ResidueCount</span><span class="sc">;</span> | 
|  | 281  268           <span class="s">}</span> | 
|  | 282  269           <span class="k">else</span> <span class="s">{</span> | 
|  | 283  270             <span class="i">$AllResiduesNameToCountMap</span>{<span class="i">$ResidueName</span>} = <span class="i">$ResidueCount</span><span class="sc">;</span> | 
|  | 284  271           <span class="s">}</span> | 
|  | 285  272         <span class="s">}</span> | 
|  | 286  273       <span class="s">}</span> | 
|  | 287  274       <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%AllResiduesNameToCountMap</span><span class="s">)</span> <span class="s">{</span> | 
|  | 288  275         <span class="i">$ResidueCount</span> = <span class="i">$AllResiduesNameToCountMap</span>{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 289  276         <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$AllResiduesCountToNameMap</span>{<span class="i">$ResidueCount</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 290  277           <span class="i">$AllResiduesCountToNameMap</span>{<span class="i">$ResidueCount</span>} .= <span class="q">"~${ResidueName}"</span><span class="sc">;</span> | 
|  | 291  278         <span class="s">}</span> | 
|  | 292  279         <span class="k">else</span> <span class="s">{</span> | 
|  | 293  280           <span class="i">$AllResiduesCountToNameMap</span>{<span class="i">$ResidueCount</span>} = <span class="i">$ResidueName</span><span class="sc">;</span> | 
|  | 294  281         <span class="s">}</span> | 
|  | 295  282       <span class="s">}</span> | 
|  | 296  283     <span class="s">}</span> | 
|  | 297  284 | 
|  | 298  285     <span class="c"># Setup distribution data for individual chains and the grand total as well...</span> | 
|  | 299  286     <span class="k">my</span><span class="s">(</span><span class="i">$ChainResidueCount</span><span class="cm">,</span> <span class="i">$PercentResidueCount</span><span class="cm">,</span> <span class="i">$TotalResidueCount</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="cm">,</span> <span class="i">@ResidueNamesList</span><span class="cm">,</span> <span class="i">%ResiduesFrequencyInfoMap</span><span class="s">)</span><span class="sc">;</span> | 
|  | 300  287     <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">ChainIDs</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 301  288     <span class="i">%</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">Frequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 302  289     <span class="i">%</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">PercentFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 303  290 | 
|  | 304  291     <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 305  292     <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsPercentFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 306  293 | 
|  | 307  294     <span class="i">$TotalResidueCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 308  295 | 
|  | 309  296     <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 310  297       <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">ChainIDs</span>}}<span class="cm">,</span> <span class="i">$ChainID</span><span class="sc">;</span> | 
|  | 311  298       <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">Frequency</span>}{<span class="i">$ChainID</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 312  299       <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">PercentFrequency</span>}{<span class="i">$ChainID</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 313  300 | 
|  | 314  301       <span class="i">$ChainResidueCount</span> = <span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">Residues</span>}{<span class="i">$ChainID</span>}}<span class="sc">;</span> | 
|  | 315  302       <span class="i">$TotalResidueCount</span> += <span class="i">$ChainResidueCount</span><span class="sc">;</span> | 
|  | 316  303 | 
|  | 317  304       <span class="k">for</span> <span class="i">$ResidueCount</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span><span class="i">$b</span> <=> <span class="i">$a</span><span class="s">}</span> <span class="k">keys</span> <span class="i">%</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 318  305         <span class="i">$ResidueNames</span> = <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ChainID</span>}{<span class="i">$ResidueCount</span>}<span class="sc">;</span> | 
|  | 319  306         <span class="i">@ResidueNamesList</span> = <span class="k">split</span> <span class="q">/~/</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="sc">;</span> | 
|  | 320  307         <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="i">@ResidueNamesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 321  308           <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">Frequency</span>}{<span class="i">$ChainID</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 322  309           <span class="i">$PercentResidueCount</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.1f"</span><span class="cm">,</span> <span class="s">(</span><span class="s">(</span><span class="i">$ResidueCount</span>/<span class="i">$ChainResidueCount</span><span class="s">)</span>*<span class="n">100</span><span class="s">)</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 323  310           <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">PercentFrequency</span>}{<span class="i">$ChainID</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${PercentResidueCount}%"</span><span class="sc">;</span> | 
|  | 324  311         <span class="s">}</span> | 
|  | 325  312       <span class="s">}</span> | 
|  | 326  313     <span class="s">}</span> | 
|  | 327  314     <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 328  315       <span class="k">for</span> <span class="i">$ResidueCount</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span><span class="i">$b</span> <=> <span class="i">$a</span><span class="s">}</span> <span class="k">keys</span> <span class="i">%AllResiduesCountToNameMap</span><span class="s">)</span> <span class="s">{</span> | 
|  | 329  316         <span class="i">$ResidueNames</span> = <span class="i">$AllResiduesCountToNameMap</span>{<span class="i">$ResidueCount</span>}<span class="sc">;</span> | 
|  | 330  317         <span class="i">@ResidueNamesList</span> = <span class="k">split</span> <span class="q">/~/</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="sc">;</span> | 
|  | 331  318         <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="i">@ResidueNamesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 332  319           <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsFrequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 333  320           <span class="i">$PercentResidueCount</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.1f"</span><span class="cm">,</span> <span class="s">(</span><span class="s">(</span><span class="i">$ResidueCount</span>/<span class="i">$TotalResidueCount</span><span class="s">)</span>*<span class="n">100</span><span class="s">)</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 334  321           <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsPercentFrequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${PercentResidueCount}%"</span><span class="sc">;</span> | 
|  | 335  322         <span class="s">}</span> | 
|  | 336  323       <span class="s">}</span> | 
|  | 337  324     <span class="s">}</span> | 
|  | 338  325 | 
|  | 339  326     <span class="c"># List distribution of residues</span> | 
|  | 340  327     <span class="k">print</span> <span class="q">"\nDistribution of residues in chain(s): \n"</span><span class="sc">;</span> | 
|  | 341  328     <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 342  329       <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 343  330         <span class="k">print</span> <span class="q">"Chain "</span><span class="cm">,</span> <span class="i">CleanupChainID</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span><span class="cm">,</span> <span class="q">": "</span><span class="sc">;</span> | 
|  | 344  331       <span class="s">}</span> | 
|  | 345  332       <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">Frequency</span>}{<span class="i">$ChainID</span>}}<span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 346  333     <span class="s">}</span> | 
|  | 347  334     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">2</span><span class="s">)</span> <span class="s">{</span> | 
|  | 348  335       <span class="k">print</span> <span class="q">"\nPercent distribution of residues in chain(s): \n"</span><span class="sc">;</span> | 
|  | 349  336       <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 350  337         <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 351  338           <span class="k">print</span> <span class="q">"Chain "</span><span class="cm">,</span> <span class="i">CleanupChainID</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span><span class="cm">,</span> <span class="q">": "</span><span class="sc">;</span> | 
|  | 352  339         <span class="s">}</span> | 
|  | 353  340         <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">PercentFrequency</span>}{<span class="i">$ChainID</span>}}<span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 354  341       <span class="s">}</span> | 
|  | 355  342     <span class="s">}</span> | 
|  | 356  343     <span class="k">if</span> <span class="s">(</span><span class="i">$ChainCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 357  344       <span class="k">print</span> <span class="q">"\nDistribution of residues across all chains: \n"</span><span class="sc">;</span> | 
|  | 358  345       <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsFrequency</span>}}<span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 359  346 | 
|  | 360  347       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">2</span><span class="s">)</span> <span class="s">{</span> | 
|  | 361  348         <span class="k">print</span> <span class="q">"\nPercent distribution of residues across all chains: \n"</span><span class="sc">;</span> | 
|  | 362  349         <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesFrequencyInfoMap</span>{<span class="w">AllChainsPercentFrequency</span>}}<span class="cm">,</span> <span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 363  350       <span class="s">}</span> | 
|  | 364  351     <span class="s">}</span> | 
|  | 365  352   <span class="s">}</span> | 
|  | 366  353 <span class="s">}</span> | 
|  | 367  354 | 
|  | 368  355 <span class="c"># List information about all the residues...</span> | 
|  | 369 <a name="ListAllResiduesInfo-"></a> 356 <span class="k">sub </span><span class="m">ListAllResiduesInfo</span> <span class="s">{</span> | 
|  | 370  357   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 371  358   <span class="k">my</span><span class="s">(</span><span class="i">$TotalResidueCount</span><span class="cm">,</span> <span class="i">$AtomResiduesCount</span><span class="cm">,</span> <span class="i">$HetatmResiduesCount</span><span class="cm">,</span> <span class="i">$ResiduesInfoRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 372  359 | 
|  | 373  360   <span class="i">$ResiduesInfoRef</span> = <span class="i">GetAllResidues</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 374  361   <span class="i">$TotalResidueCount</span> = <span class="i">@</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">ResidueNames</span>}}<span class="sc">;</span> | 
|  | 375  362   <span class="i">$AtomResiduesCount</span> = <span class="i">@</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">AtomResidueNames</span>}}<span class="sc">;</span> | 
|  | 376  363   <span class="i">$HetatmResiduesCount</span> = <span class="i">@</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">HetatmResidueNames</span>}}<span class="sc">;</span> | 
|  | 377  364 | 
|  | 378  365   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountResiduesAll</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 379  366     <span class="k">print</span> <span class="q">"\nTotal number of residues: Total - $TotalResidueCount; ATOM residues - $AtomResiduesCount; HETATM residues - $HetatmResiduesCount\n"</span><span class="sc">;</span> | 
|  | 380  367 | 
|  | 381  368     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">3</span><span class="s">)</span> <span class="s">{</span> | 
|  | 382  369       <span class="k">print</span> <span class="q">"List of residues: \n"</span><span class="sc">;</span> | 
|  | 383  370       <span class="k">if</span> <span class="s">(</span><span class="i">$AtomResiduesCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 384  371         <span class="k">print</span> <span class="q">"ATOM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">AtomResidueNames</span>}}<span class="cm">,</span> <span class="q">', '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 385  372       <span class="s">}</span> | 
|  | 386  373       <span class="k">if</span> <span class="s">(</span><span class="i">$HetatmResiduesCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 387  374         <span class="k">print</span> <span class="q">"HETATM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">HetatmResidueNames</span>}}<span class="cm">,</span> <span class="q">', '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 388  375       <span class="s">}</span> | 
|  | 389  376     <span class="s">}</span> | 
|  | 390  377   <span class="s">}</span> | 
|  | 391  378 | 
|  | 392  379   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyAll</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 393  380     <span class="k">my</span><span class="s">(</span><span class="i">$ResidueName</span><span class="cm">,</span> <span class="i">$ResidueCount</span><span class="s">)</span><span class="sc">;</span> | 
|  | 394  381 | 
|  | 395  382     <span class="c"># Setup a hash using residue count as key for sorting the values...</span> | 
|  | 396  383     <span class="k">my</span><span class="s">(</span><span class="i">%ResiduesCountToNameMap</span><span class="cm">,</span> <span class="i">%AtomResiduesCountToNameMap</span><span class="cm">,</span> <span class="i">%HetatmResiduesCountToNameMap</span><span class="s">)</span><span class="sc">;</span> | 
|  | 397  384     <span class="i">%ResiduesCountToNameMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 398  385     <span class="i">%</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 399  386 | 
|  | 400  387     <span class="i">%AtomResiduesCountToNameMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 401  388     <span class="i">%</span>{<span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 402  389 | 
|  | 403  390     <span class="i">%HetatmResiduesCountToNameMap</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 404  391     <span class="i">%</span>{<span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 405  392 | 
|  | 406  393     <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 407  394       <span class="i">$ResidueCount</span> = <span class="i">$ResiduesInfoRef</span>->{<span class="w">ResidueCount</span>}{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 408  395       <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 409  396         <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} .= <span class="q">"~${ResidueName}"</span><span class="sc">;</span> | 
|  | 410  397       <span class="s">}</span> | 
|  | 411  398       <span class="k">else</span> <span class="s">{</span> | 
|  | 412  399         <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} = <span class="i">$ResidueName</span><span class="sc">;</span> | 
|  | 413  400       <span class="s">}</span> | 
|  | 414  401     <span class="s">}</span> | 
|  | 415  402 | 
|  | 416  403     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 417  404       <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">AtomResidueCount</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 418  405         <span class="i">$ResidueCount</span> = <span class="i">$ResiduesInfoRef</span>->{<span class="w">AtomResidueCount</span>}{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 419  406         <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 420  407           <span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} .= <span class="q">"~${ResidueName}"</span><span class="sc">;</span> | 
|  | 421  408         <span class="s">}</span> | 
|  | 422  409         <span class="k">else</span> <span class="s">{</span> | 
|  | 423  410           <span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} = <span class="i">$ResidueName</span><span class="sc">;</span> | 
|  | 424  411         <span class="s">}</span> | 
|  | 425  412       <span class="s">}</span> | 
|  | 426  413       <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="k">keys</span> <span class="i">%</span>{<span class="i">$ResiduesInfoRef</span>->{<span class="w">HetatmResidueCount</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 427  414         <span class="i">$ResidueCount</span> = <span class="i">$ResiduesInfoRef</span>->{<span class="w">HetatmResidueCount</span>}{<span class="i">$ResidueName</span>}<span class="sc">;</span> | 
|  | 428  415         <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 429  416           <span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} .= <span class="q">"~${ResidueName}"</span><span class="sc">;</span> | 
|  | 430  417         <span class="s">}</span> | 
|  | 431  418         <span class="k">else</span> <span class="s">{</span> | 
|  | 432  419           <span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>} = <span class="i">$ResidueName</span><span class="sc">;</span> | 
|  | 433  420         <span class="s">}</span> | 
|  | 434  421       <span class="s">}</span> | 
|  | 435  422     <span class="s">}</span> | 
|  | 436  423 | 
|  | 437  424     <span class="c"># Setup distribution of residues info...</span> | 
|  | 438  425     <span class="k">my</span><span class="s">(</span><span class="i">$ResidueNames</span><span class="cm">,</span> <span class="i">$PercentResidueCount</span><span class="cm">,</span> <span class="i">@ResidueNamesList</span><span class="cm">,</span> <span class="i">%ResiduesCountInfoMap</span><span class="cm">,</span> <span class="i">%AtomResiduesCountInfoMap</span><span class="cm">,</span> <span class="i">%HetatmResiduesCountInfoMap</span><span class="s">)</span><span class="sc">;</span> | 
|  | 439  426 | 
|  | 440  427     <span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">Frequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 441  428     <span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 442  429     <span class="k">for</span> <span class="i">$ResidueCount</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span><span class="i">$b</span> <=> <span class="i">$a</span><span class="s">}</span> <span class="k">keys</span> <span class="i">%</span>{<span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 443  430       <span class="i">$PercentResidueCount</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.1f"</span><span class="cm">,</span> <span class="s">(</span><span class="s">(</span><span class="i">$ResidueCount</span>/<span class="i">$TotalResidueCount</span><span class="s">)</span>*<span class="n">100</span><span class="s">)</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 444  431       <span class="i">$ResidueNames</span> = <span class="i">$ResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="sc">;</span> | 
|  | 445  432       <span class="i">@ResidueNamesList</span> = <span class="k">split</span> <span class="q">/~/</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="sc">;</span> | 
|  | 446  433       <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="i">@ResidueNamesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 447  434         <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 448  435         <span class="k">push</span> <span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${PercentResidueCount}"</span><span class="sc">;</span> | 
|  | 449  436       <span class="s">}</span> | 
|  | 450  437     <span class="s">}</span> | 
|  | 451  438     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 452  439       <span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">Frequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 453  440       <span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 454  441       <span class="k">for</span> <span class="i">$ResidueCount</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span><span class="i">$b</span> <=> <span class="i">$a</span><span class="s">}</span> <span class="k">keys</span> <span class="i">%</span>{<span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 455  442         <span class="i">$PercentResidueCount</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.1f"</span><span class="cm">,</span> <span class="s">(</span><span class="s">(</span><span class="i">$ResidueCount</span>/<span class="i">$TotalResidueCount</span><span class="s">)</span>*<span class="n">100</span><span class="s">)</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 456  443         <span class="i">$ResidueNames</span> = <span class="i">$AtomResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="sc">;</span> | 
|  | 457  444         <span class="i">@ResidueNamesList</span> = <span class="k">split</span> <span class="q">/~/</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="sc">;</span> | 
|  | 458  445         <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="i">@ResidueNamesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 459  446           <span class="k">push</span> <span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 460  447           <span class="k">push</span> <span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${PercentResidueCount}"</span><span class="sc">;</span> | 
|  | 461  448         <span class="s">}</span> | 
|  | 462  449       <span class="s">}</span> | 
|  | 463  450       <span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">Frequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 464  451       <span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 465  452       <span class="k">for</span> <span class="i">$ResidueCount</span> <span class="s">(</span><span class="k">sort</span> <span class="s">{</span><span class="i">$b</span> <=> <span class="i">$a</span><span class="s">}</span> <span class="k">keys</span> <span class="i">%</span>{<span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 466  453         <span class="i">$PercentResidueCount</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.1f"</span><span class="cm">,</span> <span class="s">(</span><span class="s">(</span><span class="i">$ResidueCount</span>/<span class="i">$TotalResidueCount</span><span class="s">)</span>*<span class="n">100</span><span class="s">)</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 467  454         <span class="i">$ResidueNames</span> = <span class="i">$HetatmResiduesCountToNameMap</span>{<span class="w">ResidueNames</span>}{<span class="i">$ResidueCount</span>}<span class="sc">;</span> | 
|  | 468  455         <span class="i">@ResidueNamesList</span> = <span class="k">split</span> <span class="q">/~/</span><span class="cm">,</span> <span class="i">$ResidueNames</span><span class="sc">;</span> | 
|  | 469  456         <span class="k">for</span> <span class="i">$ResidueName</span> <span class="s">(</span><span class="i">@ResidueNamesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 470  457           <span class="k">push</span> <span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${ResidueCount}"</span><span class="sc">;</span> | 
|  | 471  458           <span class="k">push</span> <span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span> <span class="q">"${ResidueName} - ${PercentResidueCount}"</span><span class="sc">;</span> | 
|  | 472  459         <span class="s">}</span> | 
|  | 473  460       <span class="s">}</span> | 
|  | 474  461     <span class="s">}</span> | 
|  | 475  462 | 
|  | 476  463     <span class="c"># List distribution of residues</span> | 
|  | 477  464     <span class="k">print</span> <span class="q">"\nDistribution of residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 478  465     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">2</span><span class="s">)</span> <span class="s">{</span> | 
|  | 479  466       <span class="k">print</span> <span class="q">"\nPercent distribution of residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$ResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 480  467     <span class="s">}</span> | 
|  | 481  468 | 
|  | 482  469     <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 483  470       <span class="k">print</span> <span class="q">"\nDistribution of ATOM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 484  471       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">2</span><span class="s">)</span> <span class="s">{</span> | 
|  | 485  472         <span class="k">print</span> <span class="q">"\nPercent distribution of ATOM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$AtomResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 486  473       <span class="s">}</span> | 
|  | 487  474 | 
|  | 488  475       <span class="k">print</span> <span class="q">"\nDistribution of HETATM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">Frequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 489  476       <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} >= <span class="n">2</span><span class="s">)</span> <span class="s">{</span> | 
|  | 490  477         <span class="k">print</span> <span class="q">"\nPercent distribution of HETATM residues: "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@</span>{<span class="i">$HetatmResiduesCountInfoMap</span>{<span class="w">PercentFrequency</span>}}<span class="cm">,</span><span class="q">'; '</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 491  478       <span class="s">}</span> | 
|  | 492  479     <span class="s">}</span> | 
|  | 493  480   <span class="s">}</span> | 
|  | 494  481 <span class="s">}</span> | 
|  | 495  482 | 
|  | 496  483 <span class="c"># List information about residue numbers for each chain...</span> | 
|  | 497 <a name="ListResidueNumbersInfo-"></a> 484 <span class="k">sub </span><span class="m">ListResidueNumbersInfo</span> <span class="s">{</span> | 
|  | 498  485   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 499  486   <span class="k">my</span><span class="s">(</span><span class="i">$Index</span><span class="cm">,</span> <span class="i">$ResidueCount</span><span class="cm">,</span> <span class="i">$StartResidueNum</span><span class="cm">,</span> <span class="i">$EndResidueNum</span><span class="cm">,</span> <span class="i">$ChainID</span><span class="cm">,</span> <span class="i">$CollectChainResiduesBeyondTER</span><span class="cm">,</span> <span class="i">$ChainsAndResiduesInfoRef</span><span class="cm">,</span> <span class="i">$ResidueNum</span><span class="cm">,</span> <span class="i">$PreviousResidueNum</span><span class="cm">,</span> <span class="i">$ResidueName</span><span class="cm">,</span> <span class="i">$PreviousResidueName</span><span class="cm">,</span> <span class="i">$GapResiduePairsCount</span><span class="cm">,</span> <span class="i">$GapLength</span><span class="cm">,</span> <span class="i">$DescendingOrderResiduePairsCount</span><span class="cm">,</span> <span class="i">@DescendingOrderResiduePairs</span><span class="cm">,</span> <span class="i">@GapResiduePairs</span><span class="s">)</span><span class="sc">;</span> | 
|  | 500  487 | 
|  | 501  488   <span class="i">$CollectChainResiduesBeyondTER</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 502  489   <span class="i">$ChainsAndResiduesInfoRef</span> = <span class="i">GetChainsAndResidues</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="cm">,</span> <span class="q">'AtomAndHetatm'</span><span class="cm">,</span> <span class="i">$CollectChainResiduesBeyondTER</span><span class="s">)</span><span class="sc">;</span> | 
|  | 503  490 | 
|  | 504  491   <span class="k">print</span> <span class="q">"\nATOM/HETATM residue numbers information for chains:\n"</span><span class="sc">;</span> | 
|  | 505  492 | 
|  | 506  493   <span class="k">for</span> <span class="i">$ChainID</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ChainIDs</span>}}<span class="s">)</span> <span class="s">{</span> | 
|  | 507  494     <span class="k">print</span> <span class="q">"\nChain ID - "</span><span class="cm">,</span>  <span class="i">CleanupChainID</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span><span class="cm">,</span> <span class="q">""</span><span class="sc">;</span> | 
|  | 508  495 | 
|  | 509  496     <span class="i">$ResidueCount</span> = <span class="i">@</span>{<span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueNumbers</span>}{<span class="i">$ChainID</span>}}<span class="sc">;</span> | 
|  | 510  497 | 
|  | 511  498     <span class="c"># Start and end residue numbers...</span> | 
|  | 512  499     <span class="i">$StartResidueNum</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueNumbers</span>}{<span class="i">$ChainID</span>}[<span class="n">0</span>]<span class="sc">;</span> | 
|  | 513  500     <span class="i">$EndResidueNum</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueNumbers</span>}{<span class="i">$ChainID</span>}[<span class="i">$ResidueCount</span> - <span class="n">1</span>]<span class="sc">;</span> | 
|  | 514  501     <span class="k">print</span> <span class="q">"; Number of residues: $ResidueCount; Start residue number - $StartResidueNum; End residue number - $EndResidueNum\n"</span><span class="sc">;</span> | 
|  | 515  502 | 
|  | 516  503     <span class="c"># Identify any gaps in residue numbers or non-ascending order residue numbers...</span> | 
|  | 517  504     <span class="i">$GapResiduePairsCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 518  505     <span class="i">$DescendingOrderResiduePairsCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 519  506 | 
|  | 520  507     <span class="i">@DescendingOrderResiduePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 521  508     <span class="i">@GapResiduePairs</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 522  509 | 
|  | 523  510     <span class="j">RESIDUE:</span> <span class="k">for</span> <span class="i">$Index</span> <span class="s">(</span><span class="n">1</span> .. <span class="s">(</span><span class="i">$ResidueCount</span> - <span class="n">1</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 524  511       <span class="i">$ResidueNum</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueNumbers</span>}{<span class="i">$ChainID</span>}[<span class="i">$Index</span>]<span class="sc">;</span> | 
|  | 525  512       <span class="i">$PreviousResidueNum</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">ResidueNumbers</span>}{<span class="i">$ChainID</span>}[<span class="i">$Index</span> - <span class="n">1</span>]<span class="sc">;</span> | 
|  | 526  513 | 
|  | 527  514       <span class="i">$ResidueName</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">Residues</span>}{<span class="i">$ChainID</span>}[<span class="i">$Index</span>]<span class="sc">;</span> | 
|  | 528  515       <span class="i">$PreviousResidueName</span> = <span class="i">$ChainsAndResiduesInfoRef</span>->{<span class="w">Residues</span>}{<span class="i">$ChainID</span>}[<span class="i">$Index</span> - <span class="n">1</span>]<span class="sc">;</span> | 
|  | 529  516 | 
|  | 530  517       <span class="k">if</span> <span class="s">(</span><span class="i">$ResidueNum</span> == <span class="s">(</span><span class="i">$PreviousResidueNum</span> + <span class="n">1</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 531  518         <span class="c"># All is good...</span> | 
|  | 532  519         <span class="k">next</span> <span class="j">RESIDUE</span><span class="sc">;</span> | 
|  | 533  520       <span class="s">}</span> | 
|  | 534  521 | 
|  | 535  522       <span class="c"># Are residue in descending order?</span> | 
|  | 536  523       <span class="k">if</span> <span class="s">(</span><span class="i">$ResidueNum</span> < <span class="i">$PreviousResidueNum</span><span class="s">)</span> <span class="s">{</span> | 
|  | 537  524         <span class="i">$DescendingOrderResiduePairsCount</span>++<span class="sc">;</span> | 
|  | 538  525         <span class="k">push</span> <span class="i">@DescendingOrderResiduePairs</span><span class="cm">,</span> <span class="q">"<${PreviousResidueName}${PreviousResidueNum} - ${ResidueName}${ResidueNum}>"</span><span class="sc">;</span> | 
|  | 539  526       <span class="s">}</span> | 
|  | 540  527 | 
|  | 541  528       <span class="c"># Track gaps in residue pairs...</span> | 
|  | 542  529       <span class="i">$GapResiduePairsCount</span>++<span class="sc">;</span> | 
|  | 543  530       <span class="i">$GapLength</span> = <span class="k">abs</span><span class="s">(</span><span class="i">$ResidueNum</span> - <span class="i">$PreviousResidueNum</span><span class="s">)</span> - <span class="n">1</span><span class="sc">;</span> | 
|  | 544  531 | 
|  | 545  532       <span class="k">push</span> <span class="i">@GapResiduePairs</span><span class="cm">,</span> <span class="q">"<${PreviousResidueName}${PreviousResidueNum} - ${ResidueName}${ResidueNum}; $GapLength>"</span><span class="sc">;</span> | 
|  | 546  533     <span class="s">}</span> | 
|  | 547  534 | 
|  | 548  535     <span class="c"># Print gaps information...</span> | 
|  | 549  536     <span class="k">print</span> <span class="q">"Gaps in residue numbers: "</span><span class="cm">,</span> <span class="i">$GapResiduePairsCount</span> ? <span class="q">"Yes"</span> <span class="co">:</span> <span class="q">"None"</span><span class="sc">;</span> | 
|  | 550  537     <span class="k">if</span> <span class="s">(</span><span class="i">$GapResiduePairsCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 551  538       <span class="k">print</span> <span class="q">"; Number of gap residue number pairs: $GapResiduePairsCount; Gap residue pairs: <StartRes-EndRes; GapLength> - "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@GapResiduePairs</span><span class="cm">,</span> <span class="q">"; "</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="sc">;</span> | 
|  | 552  539     <span class="s">}</span> | 
|  | 553  540     <span class="k">print</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 554  541 | 
|  | 555  542     <span class="c"># Print descending residue order information...</span> | 
|  | 556  543     <span class="k">print</span> <span class="q">"Residue numbers in descending order: "</span><span class="cm">,</span> <span class="i">$DescendingOrderResiduePairsCount</span> ? <span class="q">"Yes"</span> <span class="co">:</span> <span class="q">"None"</span><span class="sc">;</span> | 
|  | 557  544     <span class="k">if</span> <span class="s">(</span><span class="i">$DescendingOrderResiduePairsCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 558  545       <span class="k">print</span> <span class="q">"; Number of descending residue number pairs: $DescendingOrderResiduePairsCount; Descending residue number pairs: <StartRes-EndRes> "</span><span class="cm">,</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@DescendingOrderResiduePairs</span><span class="cm">,</span> <span class="q">"; "</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="sc">;</span> | 
|  | 559  546     <span class="s">}</span> | 
|  | 560  547     <span class="k">print</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 561  548   <span class="s">}</span> | 
|  | 562  549 <span class="s">}</span> | 
|  | 563  550 | 
|  | 564  551 <span class="c"># List min/max XYZ coordinates for ATOM/HETATM records...</span> | 
|  | 565 <a name="ListBoundingBox-"></a> 552 <span class="k">sub </span><span class="m">ListBoundingBox</span> <span class="s">{</span> | 
|  | 566  553   <span class="k">my</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 567  554   <span class="k">my</span><span class="s">(</span><span class="i">$XMin</span><span class="cm">,</span> <span class="i">$YMin</span><span class="cm">,</span> <span class="i">$ZMin</span><span class="cm">,</span> <span class="i">$XMax</span><span class="cm">,</span> <span class="i">$YMax</span><span class="cm">,</span> <span class="i">$ZMax</span><span class="cm">,</span> <span class="i">$XSize</span><span class="cm">,</span> <span class="i">$YSize</span><span class="cm">,</span> <span class="i">$ZSize</span><span class="s">)</span><span class="sc">;</span> | 
|  | 568  555 | 
|  | 569  556   <span class="s">(</span><span class="i">$XMin</span><span class="cm">,</span> <span class="i">$YMin</span><span class="cm">,</span> <span class="i">$ZMin</span><span class="cm">,</span> <span class="i">$XMax</span><span class="cm">,</span> <span class="i">$YMax</span><span class="cm">,</span> <span class="i">$ZMax</span><span class="s">)</span> = <span class="i">GetMinMaxCoords</span><span class="s">(</span><span class="i">$PDBRecordLinesRef</span><span class="s">)</span><span class="sc">;</span> | 
|  | 570  557   <span class="i">$XSize</span> = <span class="k">abs</span><span class="s">(</span><span class="i">$XMax</span> - <span class="i">$XMin</span><span class="s">)</span><span class="sc">;</span> | 
|  | 571  558   <span class="i">$YSize</span> = <span class="k">abs</span><span class="s">(</span><span class="i">$YMax</span> - <span class="i">$YMin</span><span class="s">)</span><span class="sc">;</span> | 
|  | 572  559   <span class="i">$ZSize</span> = <span class="k">abs</span><span class="s">(</span><span class="i">$ZMax</span> - <span class="i">$ZMin</span><span class="s">)</span><span class="sc">;</span> | 
|  | 573  560 | 
|  | 574  561   <span class="i">$XMin</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$XMin</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> <span class="i">$XMax</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$XMax</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 575  562   <span class="i">$YMin</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$YMin</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> <span class="i">$YMax</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$YMax</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 576  563   <span class="i">$ZMin</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$ZMin</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> <span class="i">$ZMax</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$ZMax</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 577  564 | 
|  | 578  565   <span class="i">$XSize</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$XSize</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 579  566   <span class="i">$YSize</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$YSize</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 580  567   <span class="i">$ZSize</span> = <span class="k">sprintf</span><span class="s">(</span><span class="q">"%.3f"</span><span class="cm">,</span> <span class="i">$ZSize</span><span class="s">)</span> + <span class="n">0</span><span class="sc">;</span> | 
|  | 581  568 | 
|  | 582  569   <span class="k">print</span> <span class="q">"\nBounding box coordinates: <XMin, XMax> - <$XMin, $XMax>; <YMin, YMax> - <$YMin, $YMax>; <ZMin, ZMax> - <$ZMin, $ZMax>;\n"</span><span class="sc">;</span> | 
|  | 583  570   <span class="k">print</span> <span class="q">"Bounding box size in angstroms: XSize - $XSize; YSize - $YSize; ZSize - $ZSize\n"</span><span class="sc">;</span> | 
|  | 584  571 | 
|  | 585  572 <span class="s">}</span> | 
|  | 586  573 | 
|  | 587  574 <span class="c"># Check master record counts against actual record counts...</span> | 
|  | 588 <a name="CheckMasterRecord-"></a> 575 <span class="k">sub </span><span class="m">CheckMasterRecord</span> <span class="s">{</span> | 
|  | 589  576   <span class="k">my</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span><span class="cm">,</span> <span class="i">$PDBRecordLinesRef</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 590  577 | 
|  | 591  578   <span class="c"># Get master record information...</span> | 
|  | 592  579   <span class="k">my</span><span class="s">(</span><span class="i">$NumOfRemarkRecords</span><span class="cm">,</span> <span class="i">$NumOfHetRecords</span><span class="cm">,</span> <span class="i">$NumOfHelixRecords</span><span class="cm">,</span> <span class="i">$NumOfSheetRecords</span><span class="cm">,</span> <span class="i">$NumOfTurnRecords</span><span class="cm">,</span> <span class="i">$NumOfSiteRecords</span><span class="cm">,</span> <span class="i">$NumOfTransformationsRecords</span><span class="cm">,</span> <span class="i">$NumOfAtomAndHetatmRecords</span><span class="cm">,</span> <span class="i">$NumOfTerRecords</span><span class="cm">,</span> <span class="i">$NumOfConectRecords</span><span class="cm">,</span> <span class="i">$NumOfSeqresRecords</span><span class="s">)</span> = <span class="s">(</span><span class="k">undef</span><span class="s">)</span> x <span class="n">11</span><span class="sc">;</span> | 
|  | 593  580   <span class="k">my</span><span class="s">(</span><span class="i">$RecordLine</span><span class="cm">,</span> <span class="i">$MasterRecordFound</span><span class="s">)</span><span class="sc">;</span> | 
|  | 594  581   <span class="i">$MasterRecordFound</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 595  582 | 
|  | 596  583   <span class="j">LINE:</span> <span class="k">for</span> <span class="i">$RecordLine</span> <span class="s">(</span><span class="i">@</span>{<span class="i">$PDBRecordLinesRef</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 597  584       <span class="k">if</span> <span class="s">(</span><span class="i">IsMasterRecordType</span><span class="s">(</span><span class="i">$RecordLine</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 598  585         <span class="s">(</span><span class="i">$NumOfRemarkRecords</span><span class="cm">,</span> <span class="i">$NumOfHetRecords</span><span class="cm">,</span> <span class="i">$NumOfHelixRecords</span><span class="cm">,</span> <span class="i">$NumOfSheetRecords</span><span class="cm">,</span> <span class="i">$NumOfTurnRecords</span><span class="cm">,</span> <span class="i">$NumOfSiteRecords</span><span class="cm">,</span> <span class="i">$NumOfTransformationsRecords</span><span class="cm">,</span> <span class="i">$NumOfAtomAndHetatmRecords</span><span class="cm">,</span> <span class="i">$NumOfTerRecords</span><span class="cm">,</span> <span class="i">$NumOfConectRecords</span><span class="cm">,</span> <span class="i">$NumOfSeqresRecords</span><span class="s">)</span> = <span class="i">ParseMasterRecordLine</span><span class="s">(</span><span class="i">$RecordLine</span><span class="s">)</span><span class="sc">;</span> | 
|  | 599  586         <span class="i">$MasterRecordFound</span> = <span class="n">1</span><span class="sc">;</span> | 
|  | 600  587         <span class="k">last</span> <span class="j">LINE</span><span class="sc">;</span> | 
|  | 601  588       <span class="s">}</span> | 
|  | 602  589   <span class="s">}</span> | 
|  | 603  590   <span class="k">if</span> <span class="s">(</span>!<span class="i">$MasterRecordFound</span><span class="s">)</span> <span class="s">{</span> | 
|  | 604  591     <span class="k">print</span> <span class="q">"\nWarning: MASTER record is missing.\n"</span><span class="sc">;</span> | 
|  | 605  592     <span class="k">return</span><span class="sc">;</span> | 
|  | 606  593   <span class="s">}</span> | 
|  | 607  594   <span class="k">my</span><span class="s">(</span><span class="i">@MasterRecordValidationInfo</span><span class="s">)</span><span class="sc">;</span> | 
|  | 608  595   <span class="i">@MasterRecordValidationInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 609  596   <span class="i">$NumOfRemarkRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 610  597   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">REMARK</span>}<span class="s">)</span> && <span class="i">$NumOfRemarkRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">REMARK</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 611  598     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of REMARK records, $NumOfRemarkRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{REMARK}."</span><span class="sc">;</span> | 
|  | 612  599   <span class="s">}</span> | 
|  | 613  600   <span class="i">$NumOfHetRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 614  601   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">HET</span>}<span class="s">)</span> && <span class="i">$NumOfHetRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">HET</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 615  602     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of HET records, $NumOfHetRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{HET}."</span><span class="sc">;</span> | 
|  | 616  603   <span class="s">}</span> | 
|  | 617  604   <span class="i">$NumOfHelixRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 618  605   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">HELIX</span>}<span class="s">)</span> && <span class="i">$NumOfHelixRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">HELIX</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 619  606     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of HELIX records, $NumOfHelixRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{HELIX}."</span><span class="sc">;</span> | 
|  | 620  607   <span class="s">}</span> | 
|  | 621  608   <span class="i">$NumOfSheetRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 622  609   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SHEET</span>}<span class="s">)</span> && <span class="i">$NumOfSheetRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SHEET</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 623  610     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of SHEET records, $NumOfSheetRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{SHEET}."</span><span class="sc">;</span> | 
|  | 624  611   <span class="s">}</span> | 
|  | 625  612   <span class="i">$NumOfTurnRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 626  613   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">TURN</span>}<span class="s">)</span> && <span class="i">$NumOfTurnRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">TURN</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 627  614     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of TURN records, $NumOfTurnRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{REMARK}."</span><span class="sc">;</span> | 
|  | 628  615   <span class="s">}</span> | 
|  | 629  616   <span class="i">$NumOfSiteRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 630  617   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SITE</span>}<span class="s">)</span> && <span class="i">$NumOfSiteRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SITE</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 631  618     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of SITE records, $NumOfSiteRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{SITE}."</span><span class="sc">;</span> | 
|  | 632  619   <span class="s">}</span> | 
|  | 633  620 | 
|  | 634  621   <span class="i">$NumOfTransformationsRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 635  622   <span class="k">my</span><span class="s">(</span><span class="i">$RecordsCount</span><span class="cm">,</span> <span class="i">$ID</span><span class="cm">,</span> <span class="i">$RecordID</span><span class="cm">,</span> <span class="i">$RecordLabel</span><span class="s">)</span><span class="sc">;</span> | 
|  | 636  623   <span class="i">$RecordsCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 637  624   <span class="k">for</span> <span class="i">$RecordLabel</span> <span class="s">(</span><span class="q">'ORIGX'</span><span class="cm">,</span> <span class="q">'SCALE'</span><span class="cm">,</span> <span class="q">'MTRIX'</span><span class="s">)</span> <span class="s">{</span> | 
|  | 638  625     <span class="k">for</span> <span class="i">$ID</span> <span class="s">(</span><span class="n">1</span> .. <span class="n">3</span><span class="s">)</span> <span class="s">{</span> | 
|  | 639  626       <span class="i">$RecordID</span> = <span class="q">"${RecordLabel}${ID}"</span><span class="sc">;</span> | 
|  | 640  627       <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordID</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 641  628         <span class="i">$RecordsCount</span> += <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordID</span>}<span class="sc">;</span> | 
|  | 642  629       <span class="s">}</span> | 
|  | 643  630     <span class="s">}</span> | 
|  | 644  631   <span class="s">}</span> | 
|  | 645  632   <span class="k">if</span> <span class="s">(</span><span class="i">$NumOfTransformationsRecords</span> != <span class="i">$RecordsCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 646  633     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of transformation records (ORIGXn+SCALEn+MTRIXn), $NumOfTransformationsRecords, specified in MASTER record doen't match its explict count, $RecordsCount."</span><span class="sc">;</span> | 
|  | 647  634   <span class="s">}</span> | 
|  | 648  635 | 
|  | 649  636   <span class="i">$RecordsCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 650  637   <span class="k">for</span> <span class="i">$RecordLabel</span> <span class="s">(</span><span class="q">'ATOM'</span><span class="cm">,</span> <span class="q">'HETATM'</span><span class="s">)</span> <span class="s">{</span> | 
|  | 651  638       <span class="k">if</span> <span class="s">(</span><span class="k">exists</span> <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordLabel</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 652  639         <span class="i">$RecordsCount</span> += <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="i">$RecordLabel</span>}<span class="sc">;</span> | 
|  | 653  640       <span class="s">}</span> | 
|  | 654  641   <span class="s">}</span> | 
|  | 655  642   <span class="i">$NumOfAtomAndHetatmRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 656  643   <span class="k">if</span> <span class="s">(</span><span class="i">$NumOfAtomAndHetatmRecords</span> != <span class="i">$RecordsCount</span><span class="s">)</span> <span class="s">{</span> | 
|  | 657  644     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of ATOM + HETATM records, $NumOfAtomAndHetatmRecords, specified in MASTER record doen't match its explict count, $RecordsCount."</span><span class="sc">;</span> | 
|  | 658  645   <span class="s">}</span> | 
|  | 659  646   <span class="i">$NumOfTerRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 660  647   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">TER</span>}<span class="s">)</span> && <span class="i">$NumOfTerRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">TER</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 661  648     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of TER records, $NumOfTerRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{TER}."</span><span class="sc">;</span> | 
|  | 662  649   <span class="s">}</span> | 
|  | 663  650   <span class="i">$NumOfConectRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 664  651   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">CONECT</span>}<span class="s">)</span> && <span class="i">$NumOfConectRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">CONECT</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 665  652     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of CONECT records, $NumOfConectRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{CONECT}."</span><span class="sc">;</span> | 
|  | 666  653   <span class="s">}</span> | 
|  | 667  654   <span class="i">$NumOfSeqresRecords</span> += <span class="n">0</span><span class="sc">;</span> | 
|  | 668  655   <span class="k">if</span> <span class="s">(</span><span class="k">exists</span><span class="s">(</span><span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SEQRES</span>}<span class="s">)</span> && <span class="i">$NumOfSeqresRecords</span> != <span class="i">$RecordTypesCountRef</span>->{<span class="w">Count</span>}{<span class="w">SEQRES</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 669  656     <span class="k">push</span> <span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"Number of SITE records, $NumOfSeqresRecords, specified in MASTER record doen't match its explict count, $RecordTypesCountRef->{Count}{SEQRES}."</span><span class="sc">;</span> | 
|  | 670  657   <span class="s">}</span> | 
|  | 671  658 | 
|  | 672  659   <span class="k">if</span> <span class="s">(</span><span class="i">@MasterRecordValidationInfo</span><span class="s">)</span> <span class="s">{</span> | 
|  | 673  660     <span class="k">print</span> <span class="q">"\nMASTER record validation: Count mismatches found:\n"</span><span class="sc">;</span> | 
|  | 674  661     <span class="k">print</span> <span class="i">JoinWords</span><span class="s">(</span>\<span class="i">@MasterRecordValidationInfo</span><span class="cm">,</span> <span class="q">"\n"</span><span class="cm">,</span> <span class="n">0</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 675  662   <span class="s">}</span> | 
|  | 676  663   <span class="k">else</span> <span class="s">{</span> | 
|  | 677  664     <span class="k">print</span> <span class="q">"\nMASTER record validation: Count values match with the explicit count of the corresponding records.\n"</span><span class="sc">;</span> | 
|  | 678  665   <span class="s">}</span> | 
|  | 679  666 <span class="s">}</span> | 
|  | 680  667 | 
|  | 681  668 <span class="c"># Total size of all the files...</span> | 
|  | 682 <a name="ListTotalSizeOfFiles-"></a> 669 <span class="k">sub </span><span class="m">ListTotalSizeOfFiles</span> <span class="s">{</span> | 
|  | 683  670   <span class="k">my</span><span class="s">(</span><span class="i">$FileOkayCount</span><span class="cm">,</span> <span class="i">$TotalSize</span><span class="cm">,</span> <span class="i">$Index</span><span class="s">)</span><span class="sc">;</span> | 
|  | 684  671 | 
|  | 685  672   <span class="i">$FileOkayCount</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 686  673   <span class="i">$TotalSize</span> = <span class="n">0</span><span class="sc">;</span> | 
|  | 687  674 | 
|  | 688  675   <span class="k">for</span> <span class="i">$Index</span> <span class="s">(</span><span class="n">0</span> .. <span class="i">$#PDBFilesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 689  676     <span class="k">if</span> <span class="s">(</span><span class="i">$PDBFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$Index</span>]<span class="s">)</span> <span class="s">{</span> | 
|  | 690  677       <span class="i">$FileOkayCount</span>++<span class="sc">;</span> | 
|  | 691  678       <span class="i">$TotalSize</span> += <span class="i">$PDBFilesInfo</span>{<span class="w">FileSize</span>}[<span class="i">$Index</span>]<span class="sc">;</span> | 
|  | 692  679     <span class="s">}</span> | 
|  | 693  680   <span class="s">}</span> | 
|  | 694  681   <span class="k">if</span> <span class="s">(</span><span class="i">$FileOkayCount</span> > <span class="n">1</span><span class="s">)</span> <span class="s">{</span> | 
|  | 695  682     <span class="k">print</span> <span class="q">"\nTotal size of $FileOkayCount files: "</span><span class="cm">,</span> <span class="i">FormatFileSize</span><span class="s">(</span><span class="i">$TotalSize</span><span class="s">)</span><span class="cm">,</span> <span class="q">"\n"</span><span class="sc">;</span> | 
|  | 696  683   <span class="s">}</span> | 
|  | 697  684 | 
|  | 698  685 <span class="s">}</span> | 
|  | 699  686 | 
|  | 700  687 <span class="c"># Empty chain IDs are replaced with "None[1-9]". But for displaying purposes, take out any</span> | 
|  | 701  688 <span class="c"># numbers from label...</span> | 
|  | 702 <a name="CleanupChainID-"></a> 689 <span class="k">sub </span><span class="m">CleanupChainID</span> <span class="s">{</span> | 
|  | 703  690   <span class="k">my</span><span class="s">(</span><span class="i">$ChainID</span><span class="s">)</span> = <span class="i">@_</span><span class="sc">;</span> | 
|  | 704  691 | 
|  | 705  692   <span class="k">if</span> <span class="s">(</span><span class="i">$ChainID</span> =~ <span class="q">/^None/i</span><span class="s">)</span> <span class="s">{</span> | 
|  | 706  693     <span class="k">return</span> <span class="q">'None'</span><span class="sc">;</span> | 
|  | 707  694   <span class="s">}</span> | 
|  | 708  695   <span class="k">return</span> <span class="i">$ChainID</span><span class="sc">;</span> | 
|  | 709  696 <span class="s">}</span> | 
|  | 710  697 | 
|  | 711  698 <span class="c"># Process option values...</span> | 
|  | 712 <a name="ProcessOptions-"></a> 699 <span class="k">sub </span><span class="m">ProcessOptions</span> <span class="s">{</span> | 
|  | 713  700   <span class="i">%OptionsInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 714  701 | 
|  | 715  702   <span class="c"># Setup record types to count...</span> | 
|  | 716  703   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">count</span>}<span class="s">)</span> <span class="s">{</span> | 
|  | 717  704     <span class="i">$OptionsInfo</span>{<span class="w">CountRecordType</span>} = <span class="i">$Options</span>{<span class="w">count</span>}<span class="sc">;</span> | 
|  | 718  705   <span class="s">}</span> | 
|  | 719  706   <span class="k">else</span> <span class="s">{</span> | 
|  | 720  707     <span class="i">$OptionsInfo</span>{<span class="w">CountRecordType</span>} = <span class="i">$Options</span>{<span class="w">all</span>} ? <span class="q">'All'</span> <span class="co">:</span> <span class="q">'ATOM,HETATM'</span><span class="sc">;</span> | 
|  | 721  708   <span class="s">}</span> | 
|  | 722  709   <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedRecordTypes</span>}} =<span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 723  710   <span class="k">if</span> <span class="s">(</span><span class="i">$OptionsInfo</span>{<span class="w">CountRecordType</span>} !~ <span class="q">/^All$/i</span><span class="s">)</span> <span class="s">{</span> | 
|  | 724  711     <span class="k">my</span><span class="s">(</span><span class="i">@RecordTypes</span><span class="s">)</span><span class="sc">;</span> | 
|  | 725  712     <span class="i">@RecordTypes</span> = <span class="k">split</span> <span class="q">/\,/</span><span class="cm">,</span> <span class="i">$OptionsInfo</span>{<span class="w">CountRecordType</span>}<span class="sc">;</span> | 
|  | 726  713     <span class="k">push</span> <span class="i">@</span>{<span class="i">$OptionsInfo</span>{<span class="w">SpecifiedRecordTypes</span>}}<span class="cm">,</span> <span class="i">@RecordTypes</span><span class="sc">;</span> | 
|  | 727  714   <span class="s">}</span> | 
|  | 728  715   <span class="i">$OptionsInfo</span>{<span class="w">CountChains</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">chains</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 729  716   <span class="i">$OptionsInfo</span>{<span class="w">CheckMasterRecord</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">mastercheck</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 730  717 | 
|  | 731  718   <span class="c"># Residue count is the default. So $Options{residues} is simply ignored.</span> | 
|  | 732  719   <span class="k">my</span><span class="s">(</span><span class="i">$CountResidues</span><span class="s">)</span> = <span class="n">1</span><span class="sc">;</span> | 
|  | 733  720   <span class="i">$OptionsInfo</span>{<span class="w">CountResiduesInChains</span>} = <span class="s">(</span><span class="s">(</span><span class="i">$CountResidues</span> || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> && <span class="i">$Options</span>{<span class="w">residuesmode</span>} =~ <span class="q">/^(InChains|Both)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 734  721   <span class="i">$OptionsInfo</span>{<span class="w">CountResiduesAll</span>} = <span class="s">(</span><span class="s">(</span><span class="i">$CountResidues</span> || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> && <span class="i">$Options</span>{<span class="w">residuesmode</span>} =~ <span class="q">/^(All|Both)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 735  722 | 
|  | 736  723   <span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyInChains</span>} = <span class="s">(</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">frequency</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> && <span class="i">$Options</span>{<span class="w">residuesmode</span>} =~ <span class="q">/^(InChains|Both)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 737  724   <span class="i">$OptionsInfo</span>{<span class="w">ResiduesFrequencyAll</span>} = <span class="s">(</span><span class="s">(</span><span class="i">$Options</span>{<span class="w">frequency</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> && <span class="i">$Options</span>{<span class="w">residuesmode</span>} =~ <span class="q">/^(All|Both)$/i</span><span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 738  725 | 
|  | 739  726   <span class="i">$OptionsInfo</span>{<span class="w">ResidueNumbersInfo</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">residuenumbers</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span>  ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 740  727 | 
|  | 741  728   <span class="i">$OptionsInfo</span>{<span class="w">CalculateBoundingBox</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">boundingbox</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 742  729 | 
|  | 743  730   <span class="i">$OptionsInfo</span>{<span class="w">ListHeaderInfo</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">header</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 744  731   <span class="i">$OptionsInfo</span>{<span class="w">DetailLevel</span>} = <span class="i">$Options</span>{<span class="w">detail</span>}<span class="sc">;</span> | 
|  | 745  732 | 
|  | 746  733   <span class="i">$OptionsInfo</span>{<span class="w">ListExperimentalTechniqueInfo</span>} = <span class="s">(</span><span class="i">$Options</span>{<span class="w">experiment</span>} || <span class="i">$Options</span>{<span class="w">all</span>}<span class="s">)</span> ? <span class="n">1</span> <span class="co">:</span> <span class="n">0</span><span class="sc">;</span> | 
|  | 747  734 | 
|  | 748  735 <span class="s">}</span> | 
|  | 749  736 | 
|  | 750  737 <span class="c"># Retrieve information about PDB files...</span> | 
|  | 751 <a name="RetrievePDBFilesInfo-"></a> 738 <span class="k">sub </span><span class="m">RetrievePDBFilesInfo</span> <span class="s">{</span> | 
|  | 752  739   <span class="k">my</span><span class="s">(</span><span class="i">$Index</span><span class="cm">,</span> <span class="i">$PDBFile</span><span class="cm">,</span> <span class="i">$ModifiedTimeString</span><span class="cm">,</span> <span class="i">$ModifiedDateString</span><span class="s">)</span><span class="sc">;</span> | 
|  | 753  740 | 
|  | 754  741   <span class="i">%PDBFilesInfo</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 755  742   <span class="i">@</span>{<span class="i">$PDBFilesInfo</span>{<span class="w">FileOkay</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 756  743   <span class="i">@</span>{<span class="i">$PDBFilesInfo</span>{<span class="w">FileSize</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 757  744   <span class="i">@</span>{<span class="i">$PDBFilesInfo</span>{<span class="w">FileLastModified</span>}} = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 758  745 | 
|  | 759  746   <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">$#PDBFilesList</span><span class="s">)</span> <span class="s">{</span> | 
|  | 760  747     <span class="i">$PDBFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$Index</span>] = <span class="n">0</span><span class="sc">;</span> | 
|  | 761  748     <span class="i">$PDBFilesInfo</span>{<span class="w">FileSize</span>}[<span class="i">$Index</span>] = <span class="n">0</span><span class="sc">;</span> | 
|  | 762  749     <span class="i">$PDBFilesInfo</span>{<span class="w">FileLastModified</span>}[<span class="i">$Index</span>] = <span class="q">''</span><span class="sc">;</span> | 
|  | 763  750 | 
|  | 764  751     <span class="i">$PDBFile</span> = <span class="i">$PDBFilesList</span>[<span class="i">$Index</span>]<span class="sc">;</span> | 
|  | 765  752     <span class="k">if</span> <span class="s">(</span>!<span class="s">(</span><span class="k">-e</span> <span class="i">$PDBFile</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 766  753       <span class="k">warn</span> <span class="q">"Warning: Ignoring file $PDBFile: It doesn't exist\n"</span><span class="sc">;</span> | 
|  | 767  754       <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span> | 
|  | 768  755     <span class="s">}</span> | 
|  | 769  756     <span class="k">if</span> <span class="s">(</span>!<span class="i">CheckFileType</span><span class="s">(</span><span class="i">$PDBFile</span><span class="cm">,</span> <span class="q">"pdb"</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 770  757       <span class="k">warn</span> <span class="q">"Warning: Ignoring file $PDBFile: It's not a PDB file\n"</span><span class="sc">;</span> | 
|  | 771  758       <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span> | 
|  | 772  759     <span class="s">}</span> | 
|  | 773  760     <span class="k">if</span> <span class="s">(</span>! <span class="k">open</span> <span class="w">PDBFILE</span><span class="cm">,</span> <span class="q">"$PDBFile"</span><span class="s">)</span> <span class="s">{</span> | 
|  | 774  761       <span class="k">warn</span> <span class="q">"Warning: Ignoring file $PDBFile: Couldn't open it: $! \n"</span><span class="sc">;</span> | 
|  | 775  762       <span class="k">next</span> <span class="j">FILELIST</span><span class="sc">;</span> | 
|  | 776  763     <span class="s">}</span> | 
|  | 777  764     <span class="k">close</span> <span class="w">PDBFILE</span><span class="sc">;</span> | 
|  | 778  765 | 
|  | 779  766     <span class="i">$PDBFilesInfo</span>{<span class="w">FileOkay</span>}[<span class="i">$Index</span>] = <span class="n">1</span><span class="sc">;</span> | 
|  | 780  767     <span class="i">$PDBFilesInfo</span>{<span class="w">FileSize</span>}[<span class="i">$Index</span>] = <span class="i">FileSize</span><span class="s">(</span><span class="i">$PDBFile</span><span class="s">)</span><span class="sc">;</span> | 
|  | 781  768     <span class="s">(</span><span class="i">$ModifiedTimeString</span><span class="cm">,</span> <span class="i">$ModifiedDateString</span><span class="s">)</span> = <span class="i">FormattedFileModificationTimeAndDate</span><span class="s">(</span><span class="i">$PDBFile</span><span class="s">)</span><span class="sc">;</span> | 
|  | 782  769     <span class="i">$PDBFilesInfo</span>{<span class="w">FileLastModified</span>}[<span class="i">$Index</span>] = <span class="q">"$ModifiedTimeString; $ModifiedDateString"</span><span class="sc">;</span> | 
|  | 783  770   <span class="s">}</span> | 
|  | 784  771 <span class="s">}</span> | 
|  | 785  772 | 
|  | 786  773 | 
|  | 787  774 <span class="c"># Setup script usage  and retrieve command line arguments specified using various options...</span> | 
|  | 788 <a name="SetupScriptUsage-"></a> 775 <span class="k">sub </span><span class="m">SetupScriptUsage</span> <span class="s">{</span> | 
|  | 789  776 | 
|  | 790  777   <span class="c"># Retrieve all the options...</span> | 
|  | 791  778   <span class="i">%Options</span> = <span class="s">(</span><span class="s">)</span><span class="sc">;</span> | 
|  | 792  779   <span class="i">$Options</span>{<span class="w">count</span>} = <span class="q">''</span><span class="sc">;</span> | 
|  | 793  780   <span class="i">$Options</span>{<span class="w">detail</span>} = <span class="n">1</span><span class="sc">;</span> | 
|  | 794  781   <span class="i">$Options</span>{<span class="w">residuesmode</span>} = <span class="q">'Both'</span><span class="sc">;</span> | 
|  | 795  782 | 
|  | 796  783   <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">"all|a"</span><span class="cm">,</span> <span class="q">"boundingbox|b"</span><span class="cm">,</span> <span class="q">"count|c=s"</span><span class="cm">,</span> <span class="q">"chains"</span><span class="cm">,</span> <span class="q">"detail|d=i"</span><span class="cm">,</span> <span class="q">"experiment|e"</span><span class="cm">,</span> <span class="q">"frequency|f"</span><span class="cm">,</span> <span class="q">"mastercheck|m"</span><span class="cm">,</span> <span class="q">"header"</span><span class="cm">,</span> <span class="q">"help|h"</span><span class="cm">,</span> <span class="q">"residues"</span><span class="cm">,</span> <span class="q">"residuesmode=s"</span><span class="cm">,</span> <span class="q">"residuenumbers"</span><span class="cm">,</span> <span class="q">"workingdir|w=s"</span><span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 797  784     <span class="k">die</span> <span class="q">"\nTo get a list of valid options and their values, use \"$ScriptName -h\" or\n\"perl -S $ScriptName -h\" command and try again...\n"</span><span class="sc">;</span> | 
|  | 798  785   <span class="s">}</span> | 
|  | 799  786   <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> | 
|  | 800  787     <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> | 
|  | 801  788       <span class="k">die</span> <span class="q">"Error: The value specified, $Options{workingdir}, for option \"-w --workingdir\" is not a directory name.\n"</span><span class="sc">;</span> | 
|  | 802  789     <span class="s">}</span> | 
|  | 803  790     <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">"Error: Couldn't chdir $Options{workingdir}: $! \n"</span><span class="sc">;</span> | 
|  | 804  791   <span class="s">}</span> | 
|  | 805  792   <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">detail</span>}<span class="s">)</span><span class="s">)</span> <span class="s">{</span> | 
|  | 806  793     <span class="k">die</span> <span class="q">"Error: The value specified, $Options{detail}, for option \"-d --detail\" is not valid. Allowed values: > 0\n"</span><span class="sc">;</span> | 
|  | 807  794   <span class="s">}</span> | 
|  | 808  795   <span class="k">if</span> <span class="s">(</span><span class="i">$Options</span>{<span class="w">residuesmode</span>} !~ <span class="q">/^(InChains|All|Both)$/i</span><span class="s">)</span> <span class="s">{</span> | 
|  | 809  796     <span class="k">die</span> <span class="q">"Error: The value specified, $Options{residuesmode}, for option \"--ResiduesMode\" is not valid. Allowed values: InChains, All, or Both\n"</span><span class="sc">;</span> | 
|  | 810  797   <span class="s">}</span> | 
|  | 811  798 <span class="s">}</span> | 
|  | 812  799 | 
|  | 813 <a name="EOF-"></a></pre> | 
|  | 814 <p> </p> | 
|  | 815 <br /> | 
|  | 816 <center> | 
|  | 817 <img src="../../../images/h2o2.png"> | 
|  | 818 </center> | 
|  | 819 </body> | 
|  | 820 </html> |