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

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: flat_context.cpp,v $
  4.  * PRODUCTION Revision 1000.2  2004/06/01 19:42:58  gouriano
  5.  * PRODUCTION PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.7
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: flat_context.cpp,v 1000.2 2004/06/01 19:42:58 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:  Aaron Ucko, NCBI
  35. *
  36. * File Description:
  37. *   new (early 2003) flat-file generator -- context needed when (pre)formatting
  38. *
  39. * ===========================================================================
  40. */
  41. #include <ncbi_pch.hpp>
  42. #include <objtools/flat/flat_formatter.hpp>
  43. #include <corelib/ncbiutil.hpp>
  44. #include <objects/seq/Bioseq.hpp>
  45. #include <objects/seqset/Bioseq_set.hpp>
  46. #include <objects/seqset/Seq_entry.hpp>
  47. #include <objmgr/scope.hpp>
  48. BEGIN_NCBI_SCOPE
  49. BEGIN_SCOPE(objects)
  50. void CFlatContext::SetFlags(const CSeq_entry& entry, bool do_parents)
  51. {
  52.     if (do_parents  &&  entry.GetParentEntry()) {
  53.         SetFlags(*entry.GetParentEntry(), true);
  54.     }
  55.     if (entry.IsSet()) {
  56.         const CBioseq_set& bss = entry.GetSet();
  57.         m_SegmentCount = 0;
  58.         switch (bss.GetClass()) {
  59.         case CBioseq_set::eClass_segset:
  60.         case CBioseq_set::eClass_conset:
  61.             ITERATE (CBioseq_set::TSeq_set, it, bss.GetSeq_set()) {
  62.                 if ((*it)->IsSeq()) {
  63.                     m_SegMaster = &(*it)->GetSeq();
  64.                     break;
  65.                 }
  66.             }
  67.             break;
  68.         case CBioseq_set::eClass_parts:
  69.             m_SegmentNum   = 1;
  70.             m_SegmentCount = bss.GetSeq_set().size();
  71.             break;
  72.             // ...
  73.         default:
  74.             break;
  75.         }
  76.     } else {
  77.         if (m_SegmentCount > 0  &&  m_SegmentNum == 1) {
  78.             // Properly number individual segments
  79.             const CBioseq_set& bss = entry.GetParentEntry()->GetSet();
  80.             ITERATE (CBioseq_set::TSeq_set, it, bss.GetSeq_set()) {
  81.                 if (*it == &entry) {
  82.                     break;
  83.                 }
  84.                 ++m_SegmentNum;
  85.             }
  86.             _ASSERT(m_SegmentNum <= m_SegmentCount);
  87.         }
  88.     }
  89. }
  90. const CSeq_id& CFlatContext::GetPreferredSynonym(const CSeq_id& id) const
  91. {
  92.     if (id.IsGi()  &&  id.GetGi() == m_GI) {
  93.         return *m_PrimaryID;
  94.     }
  95.     CBioseq_Handle h = m_Handle.GetScope().GetBioseqHandle(id);
  96.     if (h == m_Handle) {
  97.         return *m_PrimaryID;
  98.     } else if (h  &&  h.GetSeqId().NotEmpty()) {
  99.         return *FindBestChoice(h.GetBioseqCore()->GetId(), CSeq_id::Score);
  100.     } else {
  101.         return id;
  102.     }
  103. }
  104. const char* CFlatContext::GetUnits(bool abbrev) const
  105. {
  106.     if (m_IsWGSMaster) {
  107.         return abbrev ? "rc" : "bases"; // XXX - not "records"?
  108.     } else if (m_IsProt) {
  109.         return abbrev ? "aa" : "residues";
  110.     } else {
  111.         return abbrev ? "bp" : "bases";
  112.     }
  113.     
  114. }
  115. END_SCOPE(objects)
  116. END_NCBI_SCOPE
  117. /*
  118. * ===========================================================================
  119. *
  120. * $Log: flat_context.cpp,v $
  121. * Revision 1000.2  2004/06/01 19:42:58  gouriano
  122. * PRODUCTION: UPGRADED [GCC34_MSVC7] Dev-tree R1.7
  123. *
  124. * Revision 1.7  2004/05/21 21:42:53  gorelenk
  125. * Added PCH ncbi_pch.hpp
  126. *
  127. * Revision 1.6  2003/11/11 01:31:31  ucko
  128. * Fix for MSVC now that CBioseq_Handle::GetSeqId returns a CConstRef.
  129. *
  130. * Revision 1.5  2003/10/08 21:10:07  ucko
  131. * For segmented sequences, find and note the "master" sequence rather
  132. * than just setting a flag.
  133. *
  134. * Revision 1.4  2003/06/02 16:06:42  dicuccio
  135. * Rearranged src/objects/ subtree.  This includes the following shifts:
  136. *     - src/objects/asn2asn --> arc/app/asn2asn
  137. *     - src/objects/testmedline --> src/objects/ncbimime/test
  138. *     - src/objects/objmgr --> src/objmgr
  139. *     - src/objects/util --> src/objmgr/util
  140. *     - src/objects/alnmgr --> src/objtools/alnmgr
  141. *     - src/objects/flat --> src/objtools/flat
  142. *     - src/objects/validator --> src/objtools/validator
  143. *     - src/objects/cddalignview --> src/objtools/cddalignview
  144. * In addition, libseq now includes six of the objects/seq... libs, and libmmdb
  145. * replaces the three libmmdb? libs.
  146. *
  147. * Revision 1.3  2003/03/21 18:49:17  ucko
  148. * Turn most structs into (accessor-requiring) classes; replace some
  149. * formerly copied fields with pointers to the original data.
  150. *
  151. * Revision 1.2  2003/03/11 15:37:51  kuznets
  152. * iterate -> ITERATE
  153. *
  154. * Revision 1.1  2003/03/10 16:39:09  ucko
  155. * Initial check-in of new flat-file generator
  156. *
  157. *
  158. * ===========================================================================
  159. */