stage/include/ost/qa/torsion_potential.hh

Go to the documentation of this file.

//------------------------------------------------------------------------------
// This file is part of the OpenStructure project <www.openstructure.org>
//
// Copyright (C) 2008-2010 by the OpenStructure authors
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 3.0 of the License, or (at your option)
// any later version.
// This library is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
// details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this library; if not, write to the Free Software Foundation, Inc.,
// 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//------------------------------------------------------------------------------
#ifndef OST_QA_TORSION_POTENTIAL_HH
#define OST_QA_TORSION_POTENTIAL_HH
/*
  This implementation of the torsion angle potential combines 3
  torsion angle pairs with 1 residues in analogy to the combined
  potential described by:
  Pascal Benkert, Silvio C. E. Tosatto, Dietmar Schomburg
  QMEAN: A comprehensive scoring function for model quality assessment
  Proteins, 2008

  Author: Pascal, Benkert, Marco Biasini
*/
#include <ost/qa/torsion_statistics.hh>
#include <boost/shared_ptr.hpp>
#include <vector>
#include <boost/scoped_ptr.hpp>

namespace ost { namespace qa {

typedef std::vector<AminoAcidSet> AminoAcidAlphabet;

struct DLLEXPORT_OST_QA TorsionPotentialOpts {
  TorsionPotentialOpts();

  int                 prev_angular_bucket_size;
  int                 central_angular_bucket_size;
  int                 next_angular_bucket_size;
  AminoAcidAlphabet   alphabet;
  Real              sigma;
  template <typename DS>
  void Serialize(DS& ds);
};

class TorsionPotential;
typedef boost::shared_ptr<TorsionPotential> TorsionPotentialPtr;

class DLLEXPORT_OST_QA TorsionPotential {
public:
  static TorsionPotentialPtr Create(TorsionStatisticsPtr statistics,
                                    const TorsionPotentialOpts& opts,
                                    bool calculate_average_energy_flag);
  static TorsionPotentialPtr LoadFromFile(const String& path);

  Real GetTotalEnergy(mol::EntityHandle entity);
  
  Real GetTotalEnergy(mol::EntityView view);

  int GetEnergyCounts() const;

  void SaveToFile(const String& path);  
  
  template <typename DS>
  void Serialize(DS& ds);
  

public:
  typedef MultiClassifier<float, int, Real, Real,
                          Real, Real, Real, Real> TorsionEnergies;
private:
  void Fill(const TorsionStatisticsPtr& stat,
            bool calculate_average_energy_flag);

  TorsionPotentialOpts               options_;
  TorsionEnergies                    energies_;  
  int                                num_torsions_;
};


}}

#endif

---