OpenMS
ElutionPeakDetection.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: Timo Sachsenberg$
32 // $Authors: Erhan Kenar, Holger Franken $
33 // --------------------------------------------------------------------------
34 
35 
36 #pragma once
37 
42 
43 namespace OpenMS
44 {
76  class OPENMS_DLLAPI ElutionPeakDetection :
77  public DefaultParamHandler, public ProgressLogger
78  {
79 public:
82 
85 
95  void detectPeaks(MassTrace& mt, std::vector<MassTrace>& single_mtraces);
96 
106  void detectPeaks(std::vector<MassTrace>& mt_vec, std::vector<MassTrace>& single_mtraces);
107 
109  void filterByPeakWidth(std::vector<MassTrace>&, std::vector<MassTrace>&);
110 
113 
116 
118  double computeApexSNR(const MassTrace&);
119 
144  void findLocalExtrema(const MassTrace& tr, const Size& num_neighboring_peaks,
145  std::vector<Size>& chrom_maxes, std::vector<Size>& chrom_mins) const;
146 
148  void smoothData(MassTrace& mt, int win_size) const;
149 
150 protected:
151  void updateMembers_() override;
152 
153 private:
154  double chrom_fwhm_;
156 
157  double min_fwhm_;
158  double max_fwhm_;
159 
162 
165 
167  void detectElutionPeaks_(MassTrace&, std::vector<MassTrace>&);
168  };
169 
170 } // namespace OpenMS
171 
A base class for all classes handling default parameters.
Definition: DefaultParamHandler.h:92
Extracts chromatographic peaks from a mass trace.
Definition: ElutionPeakDetection.h:78
double computeApexSNR(const MassTrace &)
Compute the signal to noise ratio at the apex (estimated by computeMassTraceNoise)
bool mt_snr_filtering_
Whether to apply S/N filtering.
Definition: ElutionPeakDetection.h:164
void detectElutionPeaks_(MassTrace &, std::vector< MassTrace > &)
Main function to do the work.
double chrom_peak_snr_
Definition: ElutionPeakDetection.h:155
double max_fwhm_
Definition: ElutionPeakDetection.h:158
void findLocalExtrema(const MassTrace &tr, const Size &num_neighboring_peaks, std::vector< Size > &chrom_maxes, std::vector< Size > &chrom_mins) const
Computes local extrema on a mass trace.
void detectPeaks(std::vector< MassTrace > &mt_vec, std::vector< MassTrace > &single_mtraces)
Extracts chromatographic peaks from multiple MassTraces and stores the resulting split traces in a ve...
String pw_filtering_
Type of width filtering.
Definition: ElutionPeakDetection.h:161
void detectPeaks(MassTrace &mt, std::vector< MassTrace > &single_mtraces)
Extracts chromatographic peaks from a single MassTrace and stores the resulting split traces in a vec...
void smoothData(MassTrace &mt, int win_size) const
adds smoothed_intensities to internal data of mt
double min_fwhm_
Definition: ElutionPeakDetection.h:157
void updateMembers_() override
This method is used to update extra member variables at the end of the setParameters() method.
double computeMassTraceSNR(const MassTrace &)
Compute the signal to noise ratio (estimated by computeMassTraceNoise)
ElutionPeakDetection()
Default Constructor.
double computeMassTraceNoise(const MassTrace &)
Compute noise level (as RMSE of the actual signal and the smoothed signal)
double chrom_fwhm_
Definition: ElutionPeakDetection.h:154
~ElutionPeakDetection() override
Destructor.
void filterByPeakWidth(std::vector< MassTrace > &, std::vector< MassTrace > &)
Filter out mass traces below lower 5 % quartile and above upper 95 % quartile.
A container type that gathers peaks similar in m/z and moving along retention time.
Definition: MassTrace.h:62
Base class for all classes that want to report their progress.
Definition: ProgressLogger.h:53
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