SequestOutfile.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: Timo Sachsenberg $
// $Authors: Martin Langwisch $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/CONCEPT/Exception.h>
#include <OpenMS/DATASTRUCTURES/StringUtils.h>
#include <OpenMS/METADATA/PeptideIdentification.h>
#include <OpenMS/METADATA/PeptideIdentificationList.h>
 
#include <map>
#include <vector>
#include <cmath>
 
namespace OpenMS
{
  class ProteinIdentification;
  class DateTime;
 
  class OPENMS_DLLAPI SequestOutfile
  {
public:
    SequestOutfile();
 
    SequestOutfile(const SequestOutfile & sequest_outfile);
 
    virtual ~SequestOutfile();
 
    SequestOutfile & operator=(const SequestOutfile & sequest_outfile);
 
    bool operator==(const SequestOutfile & sequest_outfile) const;
 
    void load(const std::string & result_filename, PeptideIdentificationList & peptide_identifications, ProteinIdentification & protein_identification, const double p_value_threshold, std::vector<double> & pvalues, const std::string & database = "", const bool ignore_proteins_per_peptide = false);
 
// /// retrieve the p-values from the out files
//          void getPValuesFromOutFiles(vector< pair < String, vector< double > > >& filenames_and_pvalues) throw (Exception::FileNotFound&, Exception::ParseError);
 
    bool getColumns(const std::string & line, std::vector<std::string> & substrings, Size number_of_columns, Size reference_column);
 
    void getSequences(const std::string & database_filename, const std::map<std::string, Size> & ac_position_map, std::vector<std::string> & sequences, std::vector<std::pair<std::string, Size> > & found, std::map<std::string, Size> & not_found);
 
    void getACAndACType(std::string line, std::string &accession, std::string &accession_type);
 
    void readOutHeader(const std::string & result_filename, DateTime & datetime, double & precursor_mz_value, Int & charge, Size & precursor_mass_type, Size & ion_mass_type, Size & displayed_peptides, std::string &sequest, std::string &sequest_version, std::string &database_type, Int & number_column, Int & rank_sp_column, Int & id_column, Int & mh_column, Int & delta_cn_column, Int & xcorr_column, Int & sp_column, Int & sf_column, Int & ions_column, Int & reference_column, Int & peptide_column, Int & score_column, Size & number_of_columns);
 
private:
 
    static double const_weights_[];
    static double xcorr_weights_[];
    static double delta_cn_weights_[];
    static double rank_sp_weights_[];
    static double delta_mass_weights_[];
    static Size max_pep_lens_[];
    static Size num_frags_[];
  };
 
} //namespace OpenMS