IsobaricIsotopeCorrector.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: Timo Sachsenberg $
// $Authors: Stephan Aiche, Chris Bielow $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/DATASTRUCTURES/Matrix.h>
#include <memory>
#include <vector>
 
namespace OpenMS
{
  class IsobaricQuantitationMethod;
  class IsobaricQuantifierStatistics;
  class ConsensusMap;
  class ConsensusFeature;
 
  class OPENMS_DLLAPI IsobaricIsotopeCorrector
  {
public:
    static IsobaricQuantifierStatistics correctIsotopicImpurities(const ConsensusMap& consensus_map_in,
                                                                  ConsensusMap& consensus_map_out,
                                                                  const IsobaricQuantitationMethod* quant_method);
 
    static void
    correctIsotopicImpurities(std::vector<double>& intensities,
                              const IsobaricQuantitationMethod* quant_method);
 
    // No instance methods as this is a purely static class
 
private:
    static void fillInputVector_(std::vector<double>& b,
                                 Matrix<double>& m_b,
                                 const ConsensusFeature& cf,
                                 const ConsensusMap& cm);
 
    static std::vector<double> getIntensities_(const IsobaricQuantitationMethod* quant_method,
                                               const ConsensusFeature& cf,
                                               const ConsensusMap& cm);
 
    static void solveNNLS_(const Matrix<double>& correction_matrix,
                           const Matrix<double>& m_b, Matrix<double>& m_x);
 
    static void solveNNLS_(Matrix<double>& correction_matrix,
                           std::vector<double>& b,
                           std::vector<double>& x);
 
    static void computeStats_(const std::vector<double>& m_x,
                              const std::vector<double>& x_naive,
                              const float cf_intensity,
                              const IsobaricQuantitationMethod* quant_method,
                              IsobaricQuantifierStatistics& stats);
 
    static void computeStats_(const Matrix<double>& m_x,
                              const std::vector<double>& x_naive,
                              const float cf_intensity,
                              const IsobaricQuantitationMethod* quant_method,
                              IsobaricQuantifierStatistics& stats);
 
    static float updateOutputMap_(const ConsensusMap& consensus_map_in,
                                 ConsensusMap& consensus_map_out,
                                 Size current_cf,
                                 const std::vector<double>& m_x);
 
    static float updateOutputMap_(const ConsensusMap& consensus_map_in,
                                 ConsensusMap& consensus_map_out,
                                 Size current_cf,
                                 const Matrix<double>& m_x);
  };
} // namespace