contig_item.cpp
上传用户:yhdzpy8989
上传日期:2007-06-13
资源大小:13604k
文件大小:5k
源码类别:

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: contig_item.cpp,v $
  4.  * PRODUCTION Revision 1000.1  2004/06/01 19:43:55  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.5
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: contig_item.cpp,v 1000.1 2004/06/01 19:43:55 gouriano Exp $
  10. * ===========================================================================
  11. *
  12. *                            PUBLIC DOMAIN NOTICE
  13. *               National Center for Biotechnology Information
  14. *
  15. *  This software/database is a "United States Government Work" under the
  16. *  terms of the United States Copyright Act.  It was written as part of
  17. *  the author's official duties as a United States Government employee and
  18. *  thus cannot be copyrighted.  This software/database is freely available
  19. *  to the public for use. The National Library of Medicine and the U.S.
  20. *  Government have not placed any restriction on its use or reproduction.
  21. *
  22. *  Although all reasonable efforts have been taken to ensure the accuracy
  23. *  and reliability of the software and data, the NLM and the U.S.
  24. *  Government do not and cannot warrant the performance or results that
  25. *  may be obtained by using this software or data. The NLM and the U.S.
  26. *  Government disclaim all warranties, express or implied, including
  27. *  warranties of performance, merchantability or fitness for any particular
  28. *  purpose.
  29. *
  30. *  Please cite the author in any work or product based on this material.
  31. *
  32. * ===========================================================================
  33. *
  34. * Author:  Mati Shomrat, NCBI
  35. *
  36. * File Description:
  37. *   Contig item for flat-file
  38. *
  39. */
  40. #include <ncbi_pch.hpp>
  41. #include <corelib/ncbistd.hpp>
  42. #include <objects/seq/Delta_ext.hpp>
  43. #include <objects/seq/Delta_seq.hpp>
  44. #include <objects/seq/Seq_literal.hpp>
  45. #include <objects/seq/Seq_ext.hpp>
  46. #include <objects/seq/Seg_ext.hpp>
  47. #include <objtools/format/formatter.hpp>
  48. #include <objtools/format/text_ostream.hpp>
  49. #include <objtools/format/items/contig_item.hpp>
  50. #include <objtools/format/items/flat_seqloc.hpp>
  51. #include <objtools/format/context.hpp>
  52. BEGIN_NCBI_SCOPE
  53. BEGIN_SCOPE(objects)
  54. CContigItem::CContigItem(CBioseqContext& ctx) :
  55.     CFlatItem(&ctx), m_Loc(new CSeq_loc)
  56. {
  57.     x_GatherInfo(ctx);
  58. }
  59. void CContigItem::Format
  60. (IFormatter& formatter,
  61.  IFlatTextOStream& text_os) const
  62. {
  63.     formatter.FormatContig(*this, text_os);
  64. }
  65. void CContigItem::x_GatherInfo(CBioseqContext& ctx)
  66. {
  67.     typedef CRef<CSeq_loc>  TLoc;
  68.     if ( !ctx.GetHandle().IsSetInst_Ext() ) {
  69.         return;
  70.     }
  71.     CSeq_loc_mix::Tdata& data = m_Loc->SetMix().Set();
  72.     //const CSeq_inst& inst = ctx.GetHandle().GetSeq_inst();
  73.     //if ( !inst.CanGetExt() ) {
  74.     //    return;
  75.     //}
  76.     const CSeq_ext& ext = ctx.GetHandle().GetInst_Ext();
  77.     if ( ctx.IsSegmented()  &&  ctx.HasParts() ) {
  78.         ITERATE (CSeg_ext::Tdata, it, ext.GetSeg().Get()) {
  79.             data.push_back(*it);
  80.         }
  81.     } else if ( ctx.IsDelta() ) {
  82.         ITERATE (CDelta_ext::Tdata, it, ext.GetDelta().Get()) {
  83.             if ( (*it)->IsLoc() ) {
  84.                 CSeq_loc& l = const_cast<CSeq_loc&>((*it)->GetLoc());
  85.                 CRef<CSeq_loc> lr(&l);
  86.                 data.push_back(lr);
  87.             } else {  // literal
  88.                 const CSeq_literal& lit = (*it)->GetLiteral();
  89.                 TSeqPos len = lit.CanGetLength() ? lit.GetLength() : 0;
  90.                 if ( lit.CanGetSeq_data() ) {
  91.                     // data with 0 length => gap
  92.                     if ( len == 0 ) {
  93.                         data.push_back(TLoc(new CFlatGapLoc(0)));
  94.                     } else {
  95.                         // !!! don't know what to do here
  96.                     }
  97.                 } else {
  98.                     // no data => gap
  99.                     data.push_back(TLoc(new CFlatGapLoc(len)));
  100.                 }
  101.             }
  102.         }
  103.     }
  104. }
  105. END_SCOPE(objects)
  106. END_NCBI_SCOPE
  107. /*
  108. * ===========================================================================
  109. *
  110. * $Log: contig_item.cpp,v $
  111. * Revision 1000.1  2004/06/01 19:43:55  gouriano
  112. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.5
  113. *
  114. * Revision 1.5  2004/05/21 21:42:54  gorelenk
  115. * Added PCH ncbi_pch.hpp
  116. *
  117. * Revision 1.4  2004/04/22 15:51:43  shomrat
  118. * Changes in context
  119. *
  120. * Revision 1.3  2004/02/19 18:04:01  shomrat
  121. * Implemented Contig item
  122. *
  123. * Revision 1.2  2003/12/18 17:43:32  shomrat
  124. * context.hpp moved
  125. *
  126. * Revision 1.1  2003/12/17 20:19:04  shomrat
  127. * Initial Revision (adapted from flat lib)
  128. *
  129. *
  130. * ===========================================================================
  131. */