OpenMS
PeakPickerMRM.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: Hannes Roest $
32 // $Authors: Hannes Roest $
33 // --------------------------------------------------------------------------
34 
35 #pragma once
36 
42 
46 
48 
49 #ifdef WITH_CRAWDAD
50 #include <CrawdadWrapper.h>
51 #endif
52 
53 namespace OpenMS
54 {
55 
68  class OPENMS_DLLAPI PeakPickerMRM :
69  public DefaultParamHandler
70  {
71 
72 public:
73 
75  PeakPickerMRM();
77 
79  ~PeakPickerMRM() override {}
81 
83  enum FLOATINDICES { IDX_FWHM = 0, IDX_ABUNDANCE = 1, IDX_LEFTBORDER = 2, IDX_RIGHTBORDER = 3, SIZE_OF_FLOATINDICES };
84 
92  void pickChromatogram(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
93 
94 
102  void pickChromatogram(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom, MSChromatogram& smoothed_chrom);
103 
104 protected:
105 
106  void pickChromatogramCrawdad_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
107 
108  void pickChromatogram_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
109 
113  void integratePeaks_(const MSChromatogram& chromatogram);
114 
123  Size findClosestPeak_(const MSChromatogram& chromatogram, double target_rt, Size current_peak = 0);
124 
129  void removeOverlappingPeaks_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
130 
131 
133  void updateMembers_() override;
134 
137 
138  // Members
144  double gauss_width_;
149 
151  double peak_width_;
154 
156  double sn_win_len_;
163 
165  std::vector<double> integrated_intensities_;
167  std::vector<int> left_width_;
169  std::vector<int> right_width_;
170 
175  };
176 }
177 
178 
A base class for all classes handling default parameters.
Definition: DefaultParamHandler.h:92
This class represents a Gaussian lowpass-filter which works on uniform as well as on non-uniform prof...
Definition: GaussFilter.h:73
The representation of a chromatogram.
Definition: MSChromatogram.h:57
This class implements a fast peak-picking algorithm best suited for high resolution MS data (FT-ICR-M...
Definition: PeakPickerHiRes.h:84
The PeakPickerMRM finds peaks a single chromatogram.
Definition: PeakPickerMRM.h:70
PeakPickerMRM & operator=(const PeakPickerMRM &rhs)
Assignment operator is protected for algorithm.
std::vector< int > right_width_
Temporary vector to hold the peak right widths.
Definition: PeakPickerMRM.h:169
GaussFilter gauss_
Definition: PeakPickerMRM.h:173
~PeakPickerMRM() override
Destructor.
Definition: PeakPickerMRM.h:79
double sn_win_len_
Signal to noise window length.
Definition: PeakPickerMRM.h:156
void removeOverlappingPeaks_(const MSChromatogram &chromatogram, MSChromatogram &picked_chrom)
Helper function to remove overlapping peaks in a single Chromatogram.
void pickChromatogram(const MSChromatogram &chromatogram, MSChromatogram &picked_chrom)
Finds peaks in a single chromatogram and annotates left/right borders.
UInt sgolay_frame_length_
Frame length for the SGolay smoothing.
Definition: PeakPickerMRM.h:140
String method_
Peak picker method.
Definition: PeakPickerMRM.h:162
double signal_to_noise_
Signal to noise threshold.
Definition: PeakPickerMRM.h:153
double gauss_width_
Width of the Gaussian smoothing.
Definition: PeakPickerMRM.h:144
void pickChromatogram(const MSChromatogram &chromatogram, MSChromatogram &picked_chrom, MSChromatogram &smoothed_chrom)
Finds peaks in a single chromatogram and annotates left/right borders.
bool use_gauss_
Whether to use Gaussian smoothing.
Definition: PeakPickerMRM.h:146
std::vector< double > integrated_intensities_
Temporary vector to hold the integrated intensities.
Definition: PeakPickerMRM.h:165
Size findClosestPeak_(const MSChromatogram &chromatogram, double target_rt, Size current_peak=0)
Helper function to find the closest peak in a chromatogram to "target_rt".
bool write_sn_log_messages_
Whether to write out log messages of the SN estimator.
Definition: PeakPickerMRM.h:160
void integratePeaks_(const MSChromatogram &chromatogram)
Compute peak area (peak integration)
void pickChromatogram_(const MSChromatogram &chromatogram, MSChromatogram &picked_chrom)
void pickChromatogramCrawdad_(const MSChromatogram &chromatogram, MSChromatogram &picked_chrom)
double peak_width_
Forced peak with.
Definition: PeakPickerMRM.h:151
bool remove_overlapping_
Whether to resolve overlapping peaks.
Definition: PeakPickerMRM.h:148
SignalToNoiseEstimatorMedian< MSChromatogram > snt_
Definition: PeakPickerMRM.h:174
FLOATINDICES
indices into FloatDataArrays of resulting picked chromatograms
Definition: PeakPickerMRM.h:83
void updateMembers_() override
Synchronize members with param class.
UInt sgolay_polynomial_order_
Polynomial order for the SGolay smoothing.
Definition: PeakPickerMRM.h:142
UInt sn_bin_count_
Signal to noise bin count.
Definition: PeakPickerMRM.h:158
std::vector< int > left_width_
Temporary vector to hold the peak left widths.
Definition: PeakPickerMRM.h:167
PeakPickerHiRes pp_
Definition: PeakPickerMRM.h:171
SavitzkyGolayFilter sgolay_
Definition: PeakPickerMRM.h:172
Computes the Savitzky-Golay filter coefficients using QR decomposition.
Definition: SavitzkyGolayFilter.h:104
Estimates the signal/noise (S/N) ratio of each data point in a scan by using the median (histogram ba...
Definition: SignalToNoiseEstimatorMedian.h:83
A more convenient string class.
Definition: String.h:60
unsigned int UInt
Unsigned integer type.
Definition: Types.h:94
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