InteractionRange.h

Go to the documentation of this file.

 
#ifndef INTARNA_INTERACTIONRANGE_H_
#define INTARNA_INTERACTIONRANGE_H_
 
#include "IntaRNA/general.h"
#include "IntaRNA/IndexRange.h"
#include "IntaRNA/RnaSequence.h"
 
namespace IntaRNA {
 
// dummy declaration to avoid inclusion loop
class Interaction;
 
class InteractionRange {
public:
 
    const RnaSequence * s1;
 
    const RnaSequence * s2;
 
    IndexRange r1;
 
    IndexRange r2;
 
    E_type energy;
 
public:
 
    InteractionRange(
                const RnaSequence & rna1
                , const RnaSequence & rna2
                , const IndexRange & r1 = IndexRange(0,RnaSequence::lastPos)
                , const IndexRange & r2 = IndexRange(RnaSequence::lastPos,0)
                , const E_type energy = E_INF )
        : s1(&rna1)
        , s2(&rna2)
        , r1(r1.from, (r1.to > rna1.size()?rna1.size()-1:r1.to))
        , r2((r2.from > rna2.size()?rna2.size()-1:r2.from), r2.to)
        , energy(energy)
    {
#if INTARNA_IN_DEBUG_MODE
        if (!isSane())
            throw std::runtime_error("InteractionRange("+toString(*this)+") not sane!");
#endif
    }
 
    InteractionRange( const Interaction& interaction )
        : s1(NULL), s2(NULL), r1(), r2(), energy(E_INF)
    {
        // copy data
        this->operator =(interaction);
    }
 
    virtual ~InteractionRange()
    {}
 
 
    bool isSane() const
    {
        // range of first sequence ascending
        // range of second sequence descending
        return r1.isAscending() && r2.isDescending();
    }
 
    friend std::ostream& operator<<(std::ostream& out, const InteractionRange& r)
    {
        return (out <<"("<<r.r1<<"#"<<r.r2<<")");
    }
 
    const bool operator < ( const InteractionRange &r ) const
    {
        return ( r1 < r.r1 || (!(r.r1<r1) && r2<r.r2) );
    }
 
    InteractionRange & operator= ( const Interaction & interaction );
 
};
 
} // namespace
 
#endif /* INTERACTIONRANGE_H_ */