MultiplexResolverAlgorithm.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: Lars Nilse, Timo Sachsenberg $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/DATASTRUCTURES/DefaultParamHandler.h>
#include <OpenMS/FEATUREFINDER/MultiplexDeltaMasses.h>
#include <OpenMS/KERNEL/ConsensusMap.h>
#include <OpenMS/KERNEL/MSExperiment.h>
 
#include <map>
#include <string>
#include <vector>
 
namespace OpenMS
{
  class AASequence;
  class MultiplexDeltaMassesGenerator;
 
  class OPENMS_DLLAPI MultiplexResolverAlgorithm : public DefaultParamHandler
  {
  public:
    MultiplexResolverAlgorithm();
 
    void resolve(const ConsensusMap& map_in, ConsensusMap& map_out, ConsensusMap& map_conflicts,
                 const MSExperiment& blacklist = MSExperiment()) const;
 
    std::vector<bool> checkIdentifications(const ConsensusFeature& consensus) const;
 
    std::vector<double> expectedChannelMZs(const AASequence& sequence, int charge, Size& position) const;
 
  protected:
    MultiplexDeltaMassesGenerator createGenerator_() const;
 
    static double deltaMassFromMapIndex_(const ConsensusFeature::HandleSetType& feature_handles, unsigned idx);
 
    static double matchLabelSet_(const std::vector<MultiplexDeltaMasses::DeltaMass>& pattern,
                                 const MultiplexDeltaMasses::LabelSet& label_set, int& index_label_set);
 
    bool matchDeltaMasses_(const ConsensusFeature& consensus, double detected_delta_mass_at_label_set,
                           const std::vector<MultiplexDeltaMasses::DeltaMass>& pattern,
                           double theoretical_delta_mass_at_label_set, std::vector<bool>& delta_mass_matched) const;
 
    int findMatchingPattern_(const ConsensusFeature& consensus, double detected_delta_mass_at_label_set,
                             const AASequence& sequence, const MultiplexDeltaMasses::LabelSet& label_set,
                             const std::vector<MultiplexDeltaMasses>& theoretical_patterns,
                             std::vector<bool>& delta_mass_matched, int& index_label_set) const;
 
    static bool matchesResidues_(const std::vector<MultiplexDeltaMasses::DeltaMass>& pattern, const AASequence& sequence);
 
    static double findNewMZ_(double mz, int charge, const std::vector<MultiplexDeltaMasses::DeltaMass>& pattern,
                             const std::vector<bool>& delta_mass_matched);
 
    bool isBlacklisted_(const MSExperiment& blacklist, double rt, double mz, size_t charge, const ConsensusFeature& consensus) const;
 
    ConsensusFeature completeConsensus_(const ConsensusFeature& consensus, const std::vector<MultiplexDeltaMasses::DeltaMass>& pattern,
                                        const std::vector<bool>& delta_mass_matched, int index_label_set,
                                        const MSExperiment& blacklist) const;
 
    double mass_tolerance_ = 0.1;
    double mz_tolerance_ = 10.0;
    double rt_tolerance_ = 5.0;
 
    void updateMembers_() override;
  };
} // namespace OpenMS