OpenMS
Ms2IdentificationRate.h
Go to the documentation of this file.
1 // --------------------------------------------------------------------------
2 // OpenMS -- Open-Source Mass Spectrometry
3 // --------------------------------------------------------------------------
4 // Copyright The OpenMS Team -- Eberhard Karls University Tuebingen,
5 // ETH Zurich, and Freie Universitaet Berlin 2002-2023.
6 //
7 // This software is released under a three-clause BSD license:
8 // * Redistributions of source code must retain the above copyright
9 // notice, this list of conditions and the following disclaimer.
10 // * Redistributions in binary form must reproduce the above copyright
11 // notice, this list of conditions and the following disclaimer in the
12 // documentation and/or other materials provided with the distribution.
13 // * Neither the name of any author or any participating institution
14 // may be used to endorse or promote products derived from this software
15 // without specific prior written permission.
16 // For a full list of authors, refer to the file AUTHORS.
17 // --------------------------------------------------------------------------
18 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
19 // AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
20 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
21 // ARE DISCLAIMED. IN NO EVENT SHALL ANY OF THE AUTHORS OR THE CONTRIBUTING
22 // INSTITUTIONS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
23 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
24 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
25 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
26 // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
27 // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
28 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
29 //
30 // --------------------------------------------------------------------------
31 // $Maintainer: Chris Bielow$
32 // $Authors: Patricia Scheil, Swenja Wagner$
33 // --------------------------------------------------------------------------
34 
35 #pragma once
36 
37 #include "OpenMS/QC/QCBase.h"
38 
39 #include <OpenMS/CONCEPT/Types.h>
40 #include <string>
41 #include <vector>
42 
43 namespace OpenMS
44 {
45  class FeatureMap;
46  class MSExperiment;
47  class MzTabMetaData;
48  class PeptideIdentification;
49 
57  class OPENMS_DLLAPI Ms2IdentificationRate : public QCBase
58  {
59  public:
62  Size num_peptide_identification = 0;
63  Size num_ms2_spectra = 0;
64  double identification_rate = 0.;
65  };
66 
67  private:
69  const String name_ = "Ms2IdentificationRate";
70 
72  std::vector<IdentificationRateData> rate_result_;
73 
76 
77  /*
78  * @brief Checks pepID for target/decoy
79  *
80  * Only checks the first (!) hit, all other hits are ignored
81  * Is static so that it can be used with MapUtilities::applyFunctionOnPeptideIDs() without creating a new object for each ID
82  *
83  * @param id pepID to be checked
84  * @param all_targets always returns true (if the hits aren't empty)
85  * @return true/false
86  * @throws MissingInformation if target/decoy annotation is missing
87  */
88  static bool isTargetPeptide_(const PeptideIdentification& id, bool all_targets);
89 
90  /*
91  * @brief Calculates id-rate and writes the result into a IdentificationRateData object which is appended to rate_result_
92  *
93  * @param ms2_spectra_count number of found ms2 spectra
94  * @param pep_ids_count number of found (target) peptide identifications
95  */
96  void writeResults_(Size pep_ids_count, Size ms2_spectra_count);
97 
98  public:
101 
103  virtual ~Ms2IdentificationRate() = default;
104 
118  void compute(const FeatureMap& feature_map, const MSExperiment& exp, bool assume_all_target = false);
119 
133  void compute(const std::vector<PeptideIdentification>& pep_ids, const MSExperiment& exp, bool assume_all_target = false);
134 
136  const String& getName() const override;
137 
139  const std::vector<IdentificationRateData>& getResults() const;
140 
145  QCBase::Status requirements() const override;
146 
148  };
149 
150 } // namespace OpenMS
A container for features.
Definition: FeatureMap.h:106
Stores and handles combinations of enum values, e.g. a set of flags as bits flipped in an UInt64.
Definition: FlagSet.h:54
In-Memory representation of a mass spectrometry run.
Definition: MSExperiment.h:72
This class is a metric for the QualityControl-ToppTool.
Definition: Ms2IdentificationRate.h:58
Ms2IdentificationRate()=default
Default constructor.
std::vector< IdentificationRateData > rate_result_
container that stores results
Definition: Ms2IdentificationRate.h:72
void compute(const FeatureMap &feature_map, const MSExperiment &exp, bool assume_all_target=false)
computes Ms2 Identification Rate with FeatureMap
void compute(const std::vector< PeptideIdentification > &pep_ids, const MSExperiment &exp, bool assume_all_target=false)
computes Ms2 Identification Rate with PeptideIdentifications
void addMetaDataMetricsToMzTab(MzTabMetaData &meta) const
const String & getName() const override
returns the name of the metric
QCBase::Status requirements() const override
Returns the input data requirements of the compute(...) function.
void writeResults_(Size pep_ids_count, Size ms2_spectra_count)
const std::vector< IdentificationRateData > & getResults() const
returns results
virtual ~Ms2IdentificationRate()=default
Destructor.
static bool isTargetPeptide_(const PeptideIdentification &id, bool all_targets)
Size getMS2Count_(const MSExperiment &exp)
returns number of all ms2 spectra in an MSExperiment
Structure for storing results.
Definition: Ms2IdentificationRate.h:61
all meta data of a mzTab file. Please refer to specification for documentation.
Definition: MzTab.h:141
Represents the peptide hits for a spectrum.
Definition: PeptideIdentification.h:65
This class serves as an abstract base class for all QC classes.
Definition: QCBase.h:55
A more convenient string class.
Definition: String.h:60
size_t Size
Size type e.g. used as variable which can hold result of size()
Definition: Types.h:127
Main OpenMS namespace.
Definition: FeatureDeconvolution.h:48