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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: accession_item.cpp,v $
  4.  * PRODUCTION Revision 1000.2  2004/06/01 19:43:44  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.5
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: accession_item.cpp,v 1000.2 2004/06/01 19:43:44 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. *   flat-file generator -- accession item implementation
  38. *
  39. */
  40. #include <ncbi_pch.hpp>
  41. #include <corelib/ncbistd.hpp>
  42. #include <objects/seq/Bioseq.hpp>
  43. #include <objects/seq/Seqdesc.hpp>
  44. #include <objects/seqblock/GB_block.hpp>
  45. #include <objects/seqblock/EMBL_block.hpp>
  46. #include <objects/seqloc/Textseq_id.hpp>
  47. #include <objmgr/util/sequence.hpp>
  48. #include <objmgr/seqdesc_ci.hpp>
  49. #include <objtools/format/formatter.hpp>
  50. #include <objtools/format/text_ostream.hpp>
  51. #include <objtools/format/items/accession_item.hpp>
  52. #include <objtools/format/context.hpp>
  53. #include "utils.hpp"
  54. BEGIN_NCBI_SCOPE
  55. BEGIN_SCOPE(objects)
  56. CAccessionItem::CAccessionItem(CBioseqContext& ctx) :
  57.     CFlatItem(&ctx), m_ExtraAccessions(0), m_IsSetRegion(false)
  58. {
  59.     x_GatherInfo(ctx);
  60. }
  61. void CAccessionItem::Format
  62. (IFormatter& formatter,
  63.  IFlatTextOStream& text_os) const
  64. {
  65.     formatter.FormatAccession(*this, text_os);
  66. }
  67. /***************************************************************************/
  68. /*                                  PRIVATE                                */
  69. /***************************************************************************/
  70. void CAccessionItem::x_GatherInfo(CBioseqContext& ctx)
  71. {
  72.     if ( ctx.GetPrimaryId() == 0 ) {
  73.         x_SetSkip();
  74.         return;
  75.     }
  76.     m_Accession = ctx.GetPrimaryId()->GetSeqIdString();
  77.     if ( ctx.IsWGS()  && ctx.GetLocation().IsWhole() ) {
  78.         size_t acclen = m_Accession.length();
  79.         m_WGSAccession = m_Accession;
  80.         if ( acclen == 12  &&  !NStr::EndsWith(m_WGSAccession, "000000") ) {
  81.             m_WGSAccession.replace(acclen - 6, 6, 6, '0');
  82.         } else if ( acclen == 13  &&  !NStr::EndsWith(m_WGSAccession, "0000000") ) {
  83.             m_WGSAccession.replace(acclen - 7, 7, 7, '0');
  84.         } else if ( acclen == 15  &&  !NStr::EndsWith(m_WGSAccession, "00000000") ) {
  85.             m_WGSAccession.replace(acclen - 8, 8, 8, '0');
  86.         } else {
  87.             m_WGSAccession.erase();
  88.         }
  89.     }
  90.     
  91.     const list<string>* xtra = 0;
  92.     CSeqdesc_CI gb_desc(ctx.GetHandle(), CSeqdesc::e_Genbank);
  93.     if ( gb_desc ) {
  94.         x_SetObject(*gb_desc);
  95.         xtra = &gb_desc->GetGenbank().GetExtra_accessions();
  96.     }
  97.     CSeqdesc_CI embl_desc(ctx.GetHandle(), CSeqdesc::e_Embl);
  98.     if ( embl_desc ) {
  99.         x_SetObject(*embl_desc);
  100.         xtra = &embl_desc->GetEmbl().GetExtra_acc();
  101.     }
  102.     if ( xtra != 0 ) {
  103.         ITERATE (list<string>, it, *xtra) {
  104.             if ( ValidateAccession(*it) ) {
  105.                 m_ExtraAccessions.push_back(*it);
  106.             }
  107.         }
  108.     }
  109.     if ( !ctx.GetLocation().IsWhole() ) {
  110.         // specific region is set
  111.         const CSeq_loc& loc = ctx.GetLocation();
  112.         m_Region.SetFrom(loc.GetStart(kInvalidSeqPos) + 1);
  113.         m_Region.SetTo(loc.GetEnd(kInvalidSeqPos) + 1);
  114.         m_IsSetRegion = true;
  115.     }
  116. }
  117. END_SCOPE(objects)
  118. END_NCBI_SCOPE
  119. /*
  120. * ===========================================================================
  121. *
  122. * $Log: accession_item.cpp,v $
  123. * Revision 1000.2  2004/06/01 19:43:44  gouriano
  124. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.5
  125. *
  126. * Revision 1.5  2004/05/21 21:42:54  gorelenk
  127. * Added PCH ncbi_pch.hpp
  128. *
  129. * Revision 1.4  2004/04/22 15:50:40  shomrat
  130. * Changes in context; + Region
  131. *
  132. * Revision 1.3  2004/04/13 16:44:59  shomrat
  133. * Added WGS accession
  134. *
  135. * Revision 1.2  2003/12/18 17:43:31  shomrat
  136. * context.hpp moved
  137. *
  138. * Revision 1.1  2003/12/17 20:17:45  shomrat
  139. * Initial Revision (adapted from flat lib)
  140. *
  141. *
  142. * ===========================================================================
  143. */