sequence_impl.hh

Go to the documentation of this file.

//------------------------------------------------------------------------------
// This file is part of the OpenStructure project <www.openstructure.org>
//
// Copyright (C) 2008-2011 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_SEQUENCE_HH
#define OST_SEQUENCE_HH

/*
  Author: Marco Biasini
 */
#include <ost/message.hh> 
#include <list>

#include <boost/shared_ptr.hpp>
#include <ost/generic_property.hh>
#include <ost/config.hh>
#if(OST_INFO_ENABLED)
#include <ost/info/info_fw.hh>
#endif
#include <ost/mol/residue_prop.hh>
#include <ost/mol/entity_view.hh>
#include <ost/mol/residue_view.hh>

#include <ost/seq/module_config.hh>
#include <ost/seq/invalid_sequence.hh>
namespace ost { namespace seq { namespace impl {

class SequenceImpl;

typedef boost::shared_ptr<SequenceImpl> SequenceImplPtr;

class DLLEXPORT_OST_SEQ SequenceImpl : public GenericPropContainerImpl {
public:
  static SequenceImplPtr FromString(const String& seq_name,
                                    const String& sequence_string,
                                    const String& role="UNKNOWN");

  int GetResidueIndex(int pos) const;

  int GetPos(int residue_index) const;

  int GetFirstNonGap() const;
  
  int GetLastNonGap() const;  
  
  const String& GetName() const;

  void SetName(const String& name);
  
  void Cut(int start, int n);

  void ShiftRegion(int start, int end, int amount);
  const String& GetString() const {
    return seq_string_;
  }
  void SetString(const String& seq);

  void Normalise();

  void Replace(const String& str,int start, int end);

  String GetGaplessString() const;

  int GetOffset() const;
  
  int GetIndex(const String& substr) const;
  void SetOffset(int offset);

  int GetLength() const;
  SequenceImpl(const String& seq_name, const String& sequence_string, 
               const String& role);

  char GetOneLetterCode(int position) const;

  void SetOneLetterCode(int position, char new_char);
  mol::ResidueView GetResidue(int position) const;

  mol::EntityView GetAttachedView() const;

  void AttachView(const mol::EntityView& view);

  SequenceImplPtr Copy() const;
  void AttachView(const mol::EntityView& view, const String& chain_name);

  bool HasAttachedView() const;
  
  void Append(char olc);
  
  char& operator[](size_t index)
  {
    return seq_string_[index];
  }
  char operator[](size_t index) const
  {
    return seq_string_[index];
  }
  const String& GetRole() const
  {
    return seq_role_;
  }
  
  void SetRole(const String& role)
  {
    seq_role_=role;
  }
private:

  void ShiftsFromSequence();


  int GetPosNoBounds(int index) const;

  static bool IsSequenceStringSane(const String& seq_string);

  typedef struct {
    int start;
    int shift;
  } Shift;
  String              seq_name_;
  String              seq_string_;
  String              seq_role_;
  std::list<Shift>    shifts_;
  bool                editing_;
  int                 offset_;
  mol::EntityView     attached_view_;
};

typedef std::vector<SequenceImplPtr>  SequenceList;

#if(OST_INFO_ENABLED)
SequenceImplPtr DLLEXPORT_OST_SEQ SequenceImplFromInfo(const info::InfoGroup& group);

void DLLEXPORT_OST_SEQ SequenceImplToInfo(const SequenceImplPtr& sequence,
                                     info::InfoGroup& group);
#endif


}}} //ns
#endif