OpenMS
FragmentMassError.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 
38 #include <OpenMS/QC/QCBase.h>
39 #include <vector>
40 
41 namespace OpenMS
42 {
43  class FeatureMap;
44  class MSExperiment;
45  class PeptideIdentification;
46  class WindowMower;
47 
48  class OPENMS_DLLAPI FragmentMassError : public QCBase
49  {
50  public:
52  FragmentMassError() = default;
53 
55  virtual ~FragmentMassError() = default;
56 
60  struct Statistics {
61  double average_ppm = 0;
62  double variance_ppm = 0;
63  };
64 
87  void compute(FeatureMap& fmap, const MSExperiment& exp, const QCBase::SpectraMap& map_to_spectrum, ToleranceUnit tolerance_unit = ToleranceUnit::AUTO, double tolerance = 20);
88 
112  void compute(std::vector<PeptideIdentification>& pep_ids, const ProteinIdentification::SearchParameters& search_params, const MSExperiment& exp, const QCBase::SpectraMap& map_to_spectrum,
113  ToleranceUnit tolerance_unit = ToleranceUnit::AUTO, double tolerance = 20);
114 
116  const String& getName() const override;
117 
119  const std::vector<Statistics>& getResults() const;
120 
121 
126  QCBase::Status requirements() const override;
127 
128  private:
130  std::vector<Statistics> results_;
131 
132  static void calculateFME_(PeptideIdentification& pep_id, const MSExperiment& exp, const QCBase::SpectraMap& map_to_spectrum, bool& print_warning, double tolerance,
133  FragmentMassError::ToleranceUnit tolerance_unit, double& accumulator_ppm, UInt32& counter_ppm, WindowMower& window_mower_filter);
134 
135  static void calculateVariance_(FragmentMassError::Statistics& result, const PeptideIdentification& pep_id, const UInt num_ppm);
136  };
137 
138 } // 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
Definition: FragmentMassError.h:49
void compute(std::vector< PeptideIdentification > &pep_ids, const ProteinIdentification::SearchParameters &search_params, const MSExperiment &exp, const QCBase::SpectraMap &map_to_spectrum, ToleranceUnit tolerance_unit=ToleranceUnit::AUTO, double tolerance=20)
computes FragmentMassError (FME) in ppm and Dalton (only of the first PeptideHit of each PepID)
void compute(FeatureMap &fmap, const MSExperiment &exp, const QCBase::SpectraMap &map_to_spectrum, ToleranceUnit tolerance_unit=ToleranceUnit::AUTO, double tolerance=20)
computes FragmentMassError (FME) in ppm and Dalton (only of the first PeptideHit of each PepID)
const String & getName() const override
returns the name of the metric
static void calculateVariance_(FragmentMassError::Statistics &result, const PeptideIdentification &pep_id, const UInt num_ppm)
static void calculateFME_(PeptideIdentification &pep_id, const MSExperiment &exp, const QCBase::SpectraMap &map_to_spectrum, bool &print_warning, double tolerance, FragmentMassError::ToleranceUnit tolerance_unit, double &accumulator_ppm, UInt32 &counter_ppm, WindowMower &window_mower_filter)
QCBase::Status requirements() const override
Returns the input data requirements of the compute(...) function.
const std::vector< Statistics > & getResults() const
returns results
std::vector< Statistics > results_
container that stores results
Definition: FragmentMassError.h:130
FragmentMassError()=default
Default constructor.
virtual ~FragmentMassError()=default
Destructor.
Structure for storing results: average and variance of all FragmentMassErrors in ppm.
Definition: FragmentMassError.h:60
In-Memory representation of a mass spectrometry run.
Definition: MSExperiment.h:72
Represents the peptide hits for a spectrum.
Definition: PeptideIdentification.h:65
Map to find a spectrum via its NativeID.
Definition: QCBase.h:89
This class serves as an abstract base class for all QC classes.
Definition: QCBase.h:55
ToleranceUnit
Definition: QCBase.h:75
A more convenient string class.
Definition: String.h:60
WindowMower augments the highest peaks in a sliding or jumping window.
Definition: WindowMower.h:56
OPENMS_UINT32_TYPE UInt32
Unsigned integer type (32bit)
Definition: Types.h:63
unsigned int UInt
Unsigned integer type.
Definition: Types.h:94
Main OpenMS namespace.
Definition: FeatureDeconvolution.h:48
Search parameters of the DB search.
Definition: ProteinIdentification.h:273