PeakPickerChromatogram.h

Go to the documentation of this file.

// Copyright (c) 2002-present, OpenMS Inc. -- EKU Tuebingen, ETH Zurich, and FU Berlin
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Hannes Roest $
// $Authors: Hannes Roest $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/DATASTRUCTURES/DefaultParamHandler.h>
#include <OpenMS/CONCEPT/ProgressLogger.h>
#include <OpenMS/KERNEL/MSSpectrum.h>
#include <OpenMS/KERNEL/MSChromatogram.h>
#include <OpenMS/KERNEL/ChromatogramPeak.h>
 
#include <OpenMS/PROCESSING/NOISEESTIMATION/SignalToNoiseEstimatorMedian.h>
#include <OpenMS/PROCESSING/SMOOTHING/SavitzkyGolayFilter.h>
#include <OpenMS/PROCESSING/SMOOTHING/GaussFilter.h>
 
#include <OpenMS/PROCESSING/CENTROIDING/PeakPickerHiRes.h>
 
#ifdef WITH_CRAWDAD
#include <CrawdadWrapper.h>
#endif
 
namespace OpenMS
{
 
  class OPENMS_DLLAPI PeakPickerChromatogram :
    public DefaultParamHandler
  {
 
public:
 
 
    PeakPickerChromatogram();
 
    ~PeakPickerChromatogram() override {}
 
  enum FLOATINDICES { IDX_FWHM = 0, IDX_ABUNDANCE = 1, IDX_LEFTBORDER = 2, IDX_RIGHTBORDER = 3, IDX_SN = 4, SIZE_OF_FLOATINDICES };
 
    void pickChromatogram(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
    
 
    void pickChromatogram(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom, MSChromatogram& smoothed_chrom);
 
protected:
 
    void pickChromatogramCrawdad_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
 
    void pickChromatogram_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
 
    void integratePeaks_(const MSChromatogram& chromatogram);
 
    Size findClosestPeak_(const MSChromatogram& chromatogram, double target_rt, Size current_peak = 0);
 
    void removeOverlappingPeaks_(const MSChromatogram& chromatogram, MSChromatogram& picked_chrom);
 
 
    void updateMembers_() override;
 
    PeakPickerChromatogram& operator=(const PeakPickerChromatogram& rhs);
 
    // Members
    UInt sgolay_frame_length_;
    UInt sgolay_polynomial_order_;
    double gauss_width_;
    bool use_gauss_;
    bool remove_overlapping_;
 
    double peak_width_;
    double signal_to_noise_;
 
    double sn_win_len_;
    UInt sn_bin_count_;
    bool write_sn_log_messages_;
    bool report_sn_;
    std::string method_;
 
    std::vector<double> integrated_intensities_;
    std::vector<int> left_width_;
    std::vector<int> right_width_;
    std::vector<double> apex_sn_;
 
    PeakPickerHiRes pp_;
    SavitzkyGolayFilter sgolay_;
    GaussFilter gauss_;
    SignalToNoiseEstimatorMedian<MSChromatogram > snt_;
  };
}