OpenMS  2.5.0
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-2020.
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 
136  PeakPickerMRM& operator=(const PeakPickerMRM& rhs);
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 
DefaultParamHandler.h
OpenMS::PeakPickerMRM::write_sn_log_messages_
bool write_sn_log_messages_
Whether to write out log messages of the SN estimator.
Definition: PeakPickerMRM.h:160
OpenMS::PeakPickerMRM::sgolay_polynomial_order_
UInt sgolay_polynomial_order_
Polynomial order for the SGolay smoothing.
Definition: PeakPickerMRM.h:142
OpenMS::PeakPickerMRM::left_width_
std::vector< int > left_width_
Temporary vector to hold the peak left widths.
Definition: PeakPickerMRM.h:167
OpenMS::PeakPickerMRM
The PeakPickerMRM finds peaks a single chromatogram.
Definition: PeakPickerMRM.h:68
ChromatogramPeak.h
OpenMS::PeakPickerMRM::right_width_
std::vector< int > right_width_
Temporary vector to hold the peak right widths.
Definition: PeakPickerMRM.h:169
OpenMS::PeakPickerMRM::snt_
SignalToNoiseEstimatorMedian< MSChromatogram > snt_
Definition: PeakPickerMRM.h:174
OpenMS::PeakPickerMRM::method_
String method_
Peak picker method.
Definition: PeakPickerMRM.h:162
OpenMS::String
A more convenient string class.
Definition: String.h:58
OpenMS::PeakPickerMRM::~PeakPickerMRM
~PeakPickerMRM() override
Destructor.
Definition: PeakPickerMRM.h:79
OpenMS::Size
size_t Size
Size type e.g. used as variable which can hold result of size()
Definition: Types.h:127
OpenMS::PeakPickerMRM::use_gauss_
bool use_gauss_
Whether to use Gaussian smoothing.
Definition: PeakPickerMRM.h:146
GaussFilter.h
OpenMS::PeakPickerMRM::peak_width_
double peak_width_
Forced peak with.
Definition: PeakPickerMRM.h:151
OpenMS::DefaultParamHandler
A base class for all classes handling default parameters.
Definition: DefaultParamHandler.h:91
OpenMS
Main OpenMS namespace.
Definition: FeatureDeconvolution.h:46
SavitzkyGolayFilter.h
OpenMS::PeakPickerHiRes
This class implements a fast peak-picking algorithm best suited for high resolution MS data (FT-ICR-M...
Definition: PeakPickerHiRes.h:73
OpenMS::PeakPickerMRM::integrated_intensities_
std::vector< double > integrated_intensities_
Temporary vector to hold the integrated intensities.
Definition: PeakPickerMRM.h:165
ProgressLogger.h
OpenMS::PeakPickerMRM::pp_
PeakPickerHiRes pp_
Definition: PeakPickerMRM.h:171
OpenMS::SignalToNoiseEstimatorMedian
Estimates the signal/noise (S/N) ratio of each data point in a scan by using the median (histogram ba...
Definition: SignalToNoiseEstimatorMedian.h:80
OpenMS::PeakPickerMRM::gauss_width_
double gauss_width_
Width of the Gaussian smoothing.
Definition: PeakPickerMRM.h:144
OpenMS::PeakPickerMRM::sgolay_frame_length_
UInt sgolay_frame_length_
Frame length for the SGolay smoothing.
Definition: PeakPickerMRM.h:140
OpenMS::PeakPickerMRM::gauss_
GaussFilter gauss_
Definition: PeakPickerMRM.h:173
OpenMS::PeakPickerMRM::remove_overlapping_
bool remove_overlapping_
Whether to resolve overlapping peaks.
Definition: PeakPickerMRM.h:148
OpenMS::PeakPickerMRM::signal_to_noise_
double signal_to_noise_
Signal to noise threshold.
Definition: PeakPickerMRM.h:153
OpenMS::UInt
unsigned int UInt
Unsigned integer type.
Definition: Types.h:94
OpenMS::PeakPickerMRM::sn_win_len_
double sn_win_len_
Signal to noise window length.
Definition: PeakPickerMRM.h:156
MSChromatogram.h
OpenMS::MSChromatogram
The representation of a chromatogram.
Definition: MSChromatogram.h:54
OpenMS::SavitzkyGolayFilter
Computes the Savitzky-Golay filter coefficients using QR decomposition.
Definition: SavitzkyGolayFilter.h:101
OpenMS::PeakPickerMRM::sn_bin_count_
UInt sn_bin_count_
Signal to noise bin count.
Definition: PeakPickerMRM.h:158
PeakPickerHiRes.h
OpenMS::PeakPickerMRM::sgolay_
SavitzkyGolayFilter sgolay_
Definition: PeakPickerMRM.h:172
OpenMS::GaussFilter
This class represents a Gaussian lowpass-filter which works on uniform as well as on non-uniform prof...
Definition: GaussFilter.h:74
MSSpectrum.h
OpenMS::PeakPickerMRM::FLOATINDICES
FLOATINDICES
indices into FloatDataArrays of resulting picked chromatograms
Definition: PeakPickerMRM.h:83
SignalToNoiseEstimatorMedian.h