Clipper
|
00001 00004 //C Copyright (C) 2000-2006 Kevin Cowtan and University of York 00005 //L 00006 //L This library is free software and is distributed under the terms 00007 //L and conditions of version 2.1 of the GNU Lesser General Public 00008 //L Licence (LGPL) with the following additional clause: 00009 //L 00010 //L `You may also combine or link a "work that uses the Library" to 00011 //L produce a work containing portions of the Library, and distribute 00012 //L that work under terms of your choice, provided that you give 00013 //L prominent notice with each copy of the work that the specified 00014 //L version of the Library is used in it, and that you include or 00015 //L provide public access to the complete corresponding 00016 //L machine-readable source code for the Library including whatever 00017 //L changes were used in the work. (i.e. If you make changes to the 00018 //L Library you must distribute those, but you do not need to 00019 //L distribute source or object code to those portions of the work 00020 //L not covered by this licence.)' 00021 //L 00022 //L Note that this clause grants an additional right and does not impose 00023 //L any additional restriction, and so does not affect compatibility 00024 //L with the GNU General Public Licence (GPL). If you wish to negotiate 00025 //L other terms, please contact the maintainer. 00026 //L 00027 //L You can redistribute it and/or modify the library under the terms of 00028 //L the GNU Lesser General Public License as published by the Free Software 00029 //L Foundation; either version 2.1 of the License, or (at your option) any 00030 //L later version. 00031 //L 00032 //L This library is distributed in the hope that it will be useful, but 00033 //L WITHOUT ANY WARRANTY; without even the implied warranty of 00034 //L MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00035 //L Lesser General Public License for more details. 00036 //L 00037 //L You should have received a copy of the CCP4 licence and/or GNU 00038 //L Lesser General Public License along with this library; if not, write 00039 //L to the CCP4 Secretary, Daresbury Laboratory, Warrington WA4 4AD, UK. 00040 //L The GNU Lesser General Public can also be obtained by writing to the 00041 //L Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, 00042 //L MA 02111-1307 USA 00043 00044 00045 #ifndef CLIPPER_RAMACHANDRAN 00046 #define CLIPPER_RAMACHANDRAN 00047 00048 00049 #include "clipper_types.h" 00050 00051 00052 namespace clipper { 00053 00055 00057 class Prob_phi_2d 00058 { 00059 public: 00061 void init( const int& size ); 00063 void accumulate( const ftype32 table[] ); 00065 void accumulate( const ftype& phi1, const ftype& phi2, ftype wgt = 1.0 ); 00067 void normalise(); 00069 ftype probability( const ftype& phi1, const ftype& phi2 ) const; 00071 String format() const; 00073 const ftype& data( const int& i, const int& j ) const 00074 { return data_[n*i+j]; } 00076 ftype& data( const int& i, const int& j ) 00077 { return data_[n*i+j]; } 00078 private: 00079 int n; 00080 std::vector<ftype> data_; 00081 }; 00082 00083 00085 00093 class Ramachandran : private Prob_phi_2d 00094 { 00095 public: 00097 enum TYPE { Gly, Pro, NonGlyPro, NonGly, All, Gly5, Pro5, NonGlyPro5, NonGly5, All5 }; 00099 Ramachandran() {} 00101 Ramachandran( TYPE type ); 00103 void init( TYPE type ); 00105 void set_thresholds( ftype prob_favored = 0.01, 00106 ftype prob_allowed = 0.0005 ); 00108 ftype probability( const ftype& phi, const ftype& psi ) const 00109 { return Prob_phi_2d::probability( phi, psi ); } 00111 bool favored( const ftype& phi, const ftype& psi ) const 00112 { return ( probability( phi, psi ) > p_favored ); } 00114 bool allowed( const ftype& phi, const ftype& psi ) const 00115 { return ( probability( phi, psi ) > p_allowed ); } 00116 private: 00117 ftype p_favored, p_allowed; 00118 }; 00119 00120 00121 } // namespace clipper 00122 00123 00124 #endif