MzIdentMLDOMHandler.h

Go to the documentation of this file.

// Copyright (c) 2002-present, OpenMS Inc. -- EKU Tuebingen, ETH Zurich, and FU Berlin
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Mathias Walzer $
// $Authors: Mathias Walzer$
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/CHEMISTRY/AASequence.h>
#include <OpenMS/CONCEPT/UniqueIdGenerator.h>
#include <OpenMS/FORMAT/ControlledVocabulary.h>
#include <OpenMS/FORMAT/HANDLERS/XMLHandler.h>
#include <OpenMS/KERNEL/StandardTypes.h>
#include <OpenMS/METADATA/CVTermList.h>
#include <OpenMS/METADATA/PeptideIdentification.h>
#include <OpenMS/METADATA/PeptideIdentificationList.h>
#include <OpenMS/METADATA/ProteinHit.h>
#include <OpenMS/METADATA/ProteinIdentification.h>
 
#include <xercesc/dom/DOM.hpp>
#include <xercesc/dom/DOMDocument.hpp>
#include <xercesc/dom/DOMDocumentType.hpp>
#include <xercesc/dom/DOMElement.hpp>
#include <xercesc/dom/DOMImplementation.hpp>
#include <xercesc/dom/DOMImplementationLS.hpp>
#include <xercesc/dom/DOMNodeIterator.hpp>
#include <xercesc/dom/DOMNodeList.hpp>
#include <xercesc/dom/DOMText.hpp>
#include <xercesc/framework/LocalFileFormatTarget.hpp>
#include <xercesc/framework/psvi/XSValue.hpp>
#include <xercesc/parsers/XercesDOMParser.hpp>
#include <xercesc/util/OutOfMemoryException.hpp>
#include <xercesc/util/PlatformUtils.hpp>
#include <xercesc/util/XMLString.hpp>
#include <xercesc/util/XMLUni.hpp>
 
#include <list>
#include <map>
#include <set>
#include <stdexcept>
#include <string>
#include <vector>
 
// Error codes
//enum {
//   ERROR_ARGS = 1,
//   ERROR_XERCES_INIT,
//   ERROR_PARSE,
//   ERROR_EMPTY_DOCUMENT
//};
 
namespace OpenMS
{
  class ProgressLogger;
 
  namespace Internal
  {
    class OPENMS_DLLAPI MzIdentMLDOMHandler
    {
public:
      MzIdentMLDOMHandler(const std::vector<ProteinIdentification>& pro_id, const PeptideIdentificationList& pep_id, const std::string& version, const ProgressLogger& logger);
 
      MzIdentMLDOMHandler(std::vector<ProteinIdentification>& pro_id, PeptideIdentificationList& pep_id, const std::string& version, const ProgressLogger& logger);
 
      virtual ~MzIdentMLDOMHandler();
 
      void readMzIdentMLFile(const std::string& mzid_file);
      void writeMzIdentMLFile(const std::string& mzid_file);
 
protected:
      const ProgressLogger& logger_;
 
      ControlledVocabulary cv_;
      ControlledVocabulary unimod_;
 
      std::vector<ProteinIdentification>* pro_id_ = nullptr;
      PeptideIdentificationList* pep_id_ = nullptr;
 
      const std::vector<ProteinIdentification>* cpro_id_ = nullptr;
      const PeptideIdentificationList* cpep_id_ = nullptr;
 
      const std::string schema_version_;
 
      ControlledVocabulary::CVTerm getChildWithName_(const std::string& parent_accession, const std::string& name) const;
 
      void initScoreTermCaches_();
 
 
      std::pair<CVTermList, std::map<std::string, DataValue> > parseParamGroup_(xercesc::DOMNodeList* paramGroup);
      CVTerm parseCvParam_(xercesc::DOMElement* param);
      std::pair<std::string, DataValue> parseUserParam_(xercesc::DOMElement* param);
      void parseAnalysisSoftwareList_(xercesc::DOMNodeList* analysisSoftwareElements);
      void parseDBSequenceElements_(xercesc::DOMNodeList* dbSequenceElements);
      void parsePeptideElements_(xercesc::DOMNodeList* peptideElements);
      //AASequence parsePeptideSiblings_(xercesc::DOMNodeList* peptideSiblings);
      AASequence parsePeptideSiblings_(xercesc::DOMElement* peptide);
      SignedSize inferModificationLocation_(xercesc::DOMElement* modification_element, Size peptide_length) const;
      void parsePeptideEvidenceElements_(xercesc::DOMNodeList* peptideEvidenceElements);
      void parseSpectrumIdentificationElements_(xercesc::DOMNodeList* spectrumIdentificationElements);
      void parseSpectrumIdentificationProtocolElements_(xercesc::DOMNodeList* spectrumIdentificationProtocolElements);
      void parseInputElements_(xercesc::DOMNodeList* inputElements);
      void parseSpectrumIdentificationListElements_(xercesc::DOMNodeList* spectrumIdentificationListElements);
      void parseSpectrumIdentificationItemSetXLMS(std::set<std::string>::const_iterator set_it, std::multimap<std::string, int> xl_val_map, xercesc::DOMElement* element_res, const std::string& spectrumID);
      void parseSpectrumIdentificationItemElement_(xercesc::DOMElement* spectrumIdentificationItemElement, PeptideIdentification& spectrum_identification, std::string& spectrumIdentificationList_ref);
      void parseProteinDetectionHypothesisElement_(xercesc::DOMElement* proteinDetectionHypothesisElement, ProteinIdentification& protein_identification);
      void parseProteinAmbiguityGroupElement_(xercesc::DOMElement* proteinAmbiguityGroupElement, ProteinIdentification& protein_identification);
      void parseProteinDetectionListElements_(xercesc::DOMNodeList* proteinDetectionListElements);
      static ProteinIdentification::SearchParameters findSearchParameters_(std::pair<CVTermList, std::map<std::string, DataValue> > as_params);
 
      void buildCvList_(xercesc::DOMElement* cvElements);
      void buildAnalysisSoftwareList_(xercesc::DOMElement* analysisSoftwareElements);
      void buildSequenceCollection_(xercesc::DOMElement* sequenceCollectionElements);
      void buildAnalysisCollection_(xercesc::DOMElement* analysisCollectionElements);
      void buildAnalysisProtocolCollection_(xercesc::DOMElement* protocolElements);
      void buildInputDataCollection_(xercesc::DOMElement* inputElements);
      void buildEnclosedCV_(xercesc::DOMElement* parentElement, const std::string& encel, const std::string& acc, const std::string& name, const std::string& cvref);
      void buildAnalysisDataCollection_(xercesc::DOMElement* analysisElements);
 
 
private:
      MzIdentMLDOMHandler();
      MzIdentMLDOMHandler(const MzIdentMLDOMHandler& rhs);
      MzIdentMLDOMHandler& operator=(const MzIdentMLDOMHandler& rhs);
 
      struct AnalysisSoftware
      {
        std::string name;
        std::string version;
      };
      struct PeptideEvidence
      {
        int start;
        int stop;
        char pre;
        char post;
        bool idec;
      };
      struct DBSequence
      {
        std::string sequence;
        std::string database_ref;
        std::string accession;
        CVTermList cvs;
      };
      struct SpectrumIdentification
      {
        std::string spectra_data_ref;
        std::string search_database_ref;
        std::string spectrum_identification_protocol_ref;
        std::string spectrum_identification_list_ref;
      };
      struct ModificationParam
      {
        std::string fixed_mod;
        long double mass_delta;
        std::string residues;
        CVTermList modification_param_cvs;
        CVTermList specificities;
      };
      struct SpectrumIdentificationProtocol
      {
        CVTerm searchtype;
        std::string enzyme;
        CVTermList parameter_cvs;
        std::map<std::string, DataValue> parameter_ups;
//        std::vector<ModificationParam> modification_parameter;
        CVTermList modification_parameter;
        long double precursor_tolerance;
        long double fragment_tolerance;
        CVTermList threshold_cvs;
        std::map<std::string, DataValue> threshold_ups;
      };
      struct DatabaseInput
      {
        std::string name;
        std::string location;
        std::string version;
        DateTime date;
      };
 
      XMLCh* xml_root_tag_ptr_;
      XMLCh* xml_cvparam_tag_ptr_;
      XMLCh* xml_name_attr_ptr_;
 
      xercesc::XercesDOMParser mzid_parser_;
 
      std::unique_ptr<XMLHandler> xml_handler_ = nullptr;
 
      //from AnalysisSoftware
      std::string search_engine_;
      std::string search_engine_version_;
      //mapping from AnalysisSoftware
      std::map<std::string, AnalysisSoftware> as_map_; 
 
      //mapping from DataCollection Inputs
      std::map<std::string, std::string> sr_map_; 
      std::map<std::string, std::string> sd_map_; 
      std::map<std::string, DatabaseInput> db_map_; 
 
      //mapping from SpectrumIdentification - SpectrumIdentification will be the new IdentificationRuns
      std::map<std::string, SpectrumIdentification> si_map_; 
      std::map<std::string, size_t> si_pro_map_; 
 
      //mapping from SpectrumIdentificationProtocol
      std::map<std::string, SpectrumIdentificationProtocol> sp_map_; 
 
      //mapping from SequenceCollection
      std::map<std::string, AASequence> pep_map_; 
      std::map<std::string, PeptideEvidence> pe_ev_map_; 
      std::map<std::string, std::string> pv_db_map_; 
      std::multimap<std::string, std::string> p_pv_map_; 
      std::map<std::string, DBSequence> db_sq_map_; 
 
      std::list<std::list<std::string> > hit_pev_; 
 
      bool xl_ms_search_; 
      std::map<std::string, std::string> xl_id_donor_map_; 
      //std::map<std::string, std::string> xl_id_acceptor_map_; ///< mapping Peptide id -> crosslink acceptor value
      std::map<std::string, std::string> xl_id_acceptor_map_; 
      std::map<std::string, SignedSize> xl_donor_pos_map_; 
      std::map<std::string, SignedSize> xl_acceptor_pos_map_; 
      std::map<std::string, double> xl_mass_map_; 
      std::map<std::string, std::string> xl_mod_map_; 
 
      std::set<std::string> q_score_child_terms_;
      std::set<std::string> e_score_child_terms_;
      std::set<std::string> specific_score_child_terms_;
    };
  } // namespace Internal
} // namespace OpenMS