LayerDataBase.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: Marc Sturm $
// --------------------------------------------------------------------------
 
#pragma once
 
#include <OpenMS/DATASTRUCTURES/Param.h>
#include <OpenMS/VISUAL/OpenMS_GUIConfig.h>
 
#include <OpenMS/DATASTRUCTURES/StringUtils.h>
#include <OpenMS/PROCESSING/MISC/DataFilters.h>
#include <OpenMS/KERNEL/ConsensusMap.h>
#include <OpenMS/KERNEL/FeatureMap.h>
#include <OpenMS/KERNEL/MSExperiment.h>
#include <OpenMS/KERNEL/OnDiscMSExperiment.h>
#include <OpenMS/KERNEL/StandardTypes.h>
#include <OpenMS/METADATA/PeptideIdentification.h>
#include <OpenMS/METADATA/ProteinIdentification.h>
#include <OpenMS/METADATA/AnnotatedMSRun.h>
#include <OpenMS/VISUAL/LogWindow.h>
#include <OpenMS/VISUAL/MISC/CommonDefs.h>
#include <OpenMS/VISUAL/MultiGradient.h>
#include <OpenMS/VISUAL/OpenMS_GUIConfig.h>
 
#include <memory>
 
#include <bitset>
#include <vector>
 
class QWidget;
 
namespace OpenMS
{
  class LayerData1DBase;
  class LayerStoreData;
  class LayerStatistics;
  class OSWData;
  class Painter1DBase;
  class Painter2DBase;
 
  template <int N_DIM> class DimMapper;
 
  struct LayerDataDefs
  {
 
    struct ProjectionData
    {
      ProjectionData();
      ProjectionData(ProjectionData&&);
      ~ProjectionData(); // needs to be implemented in cpp, since inline would require a full definition of LayerData1DBase;
 
      std::unique_ptr<LayerData1DBase> projection_ontoX;
      std::unique_ptr<LayerData1DBase> projection_ontoY;
      struct Summary {
        UInt number_of_datapoints {0};
        Peak1D::IntensityType max_intensity {0};
        double sum_intensity {0}; // double since sum could get large
      } stats;
    };
 
    enum DataType
    {
      DT_PEAK,         
      DT_CHROMATOGRAM, 
      DT_FEATURE,      
      DT_CONSENSUS,    
      DT_IDENT,        
      DT_UNKNOWN       
    };
 
    enum Flags
    {
      F_HULL,        
      F_HULLS,       
      F_UNASSIGNED,  
      P_PRECURSORS,  
      P_PROJECTIONS, 
      C_ELEMENTS,    
      I_PEPTIDEMZ,   
      I_LABELS,      
      SIZE_OF_FLAGS
    };
 
    enum LabelType
    {
      L_NONE,       
      L_INDEX,      
      L_META_LABEL, 
      L_ID,         
      L_ID_ALL,     
      SIZE_OF_LABEL_TYPE
    };
 
    static const std::string NamesOfLabelType[SIZE_OF_LABEL_TYPE];
 
    typedef FeatureMap FeatureMapType;
 
    typedef std::shared_ptr<FeatureMap> FeatureMapSharedPtrType;
 
    typedef ConsensusMap ConsensusMapType;
 
    typedef std::shared_ptr<ConsensusMap> ConsensusMapSharedPtrType;
 
    typedef AnnotatedMSRun ExperimentType;
 
    typedef std::shared_ptr<ExperimentType> ExperimentSharedPtrType;
 
    typedef std::shared_ptr<const ExperimentType> ConstExperimentSharedPtrType;
 
    typedef std::shared_ptr<OnDiscMSExperiment> ODExperimentSharedPtrType;
 
    typedef std::shared_ptr<OSWData> OSWDataSharedPtrType;
  };
 
#ifdef _MSC_VER
  #pragma warning(disable : 4250) // 'class1' : inherits 'class2::member' via dominance
#endif
  class OPENMS_GUI_DLLAPI LayerDataBase : public LayerDataDefs
  {
  public:
    std::bitset<SIZE_OF_FLAGS> flags;
 
 
    LayerDataBase() = delete;     // <-- this is the problem. Call assignment op in 1DPeak???
    explicit LayerDataBase(const DataType type) : type(type) {}
    LayerDataBase(const LayerDataBase& ld) = default;
    LayerDataBase& operator=(const LayerDataBase& ld) = delete;
    LayerDataBase(LayerDataBase&& ld) = delete;
    LayerDataBase& operator=(LayerDataBase&& ld) = delete;
    virtual ~LayerDataBase() = default;
 
    virtual std::unique_ptr<Painter2DBase> getPainter2D() const = 0;
 
 
    virtual std::unique_ptr <LayerData1DBase> to1DLayer() const = 0;
 
    virtual std::unique_ptr<LayerStoreData> storeVisibleData(const RangeAllType& /*visible_range*/, const DataFilters& /*layer_filters*/) const
    {
      throw Exception::NotImplemented(__FILE__, __LINE__, OPENMS_PRETTY_FUNCTION);
    }
    virtual std::unique_ptr<LayerStoreData> storeFullData() const
    {
      throw Exception::NotImplemented(__FILE__, __LINE__, OPENMS_PRETTY_FUNCTION);
    }
 
    // todo: put this into a LayerData2DBase class, since a LayerData1DPeak should not implement this.
    virtual ProjectionData getProjection(const DIM_UNIT unit_x, const DIM_UNIT unit_y, const RangeAllType& area) const = 0;
 
    virtual PeakIndex findClosestDataPoint(const RangeAllType& area) const
    {
      (void)area; // allow doxygen to document the param
      throw Exception::NotImplemented(__FILE__, __LINE__, OPENMS_PRETTY_FUNCTION);
    }
 
    virtual PeakIndex findHighestDataPoint(const RangeAllType& area) const
    {
      (void)area; // allow doxygen to document the param
      throw Exception::NotImplemented(__FILE__, __LINE__, OPENMS_PRETTY_FUNCTION);
    }
 
 
    virtual PointXYType peakIndexToXY(const PeakIndex& peak, const DimMapper<2>& mapper) const = 0;
 
    virtual std::string getDataArrayDescription(const PeakIndex& peak_index)
    {
      (void)peak_index; // allow doxygen to document the param
      throw Exception::NotImplemented(__FILE__, __LINE__, OPENMS_PRETTY_FUNCTION);
    }
 
    virtual bool annotate(const PeptideIdentificationList& identifications,
                          const std::vector<ProteinIdentification>& protein_identifications);
 
 
    virtual void updateRanges() = 0;
 
    float getMinIntensity() const;
 
    float getMaxIntensity() const;
 
    using RangeAllType = RangeManager<RangeRT, RangeMZ, RangeIntensity, RangeMobility>;
 
    virtual RangeAllType getRange() const = 0;
 
    virtual std::unique_ptr<LayerStatistics> getStats() const = 0;
 
    const std::string& getName() const
    {
      return name_;
    }
    void setName(const std::string& new_name)
    {
      name_ = new_name;
    }
 
    const std::string& getNameSuffix() const
    {
      return name_suffix_;
    }
    void setNameSuffix(const std::string& decorator)
    {
      name_suffix_ = decorator;
    }
 
    virtual std::string getDecoratedName() const;
 
    bool visible = true;
 
    DataType type = DT_UNKNOWN;
 
    std::string filename;
 
    Param param;
 
    MultiGradient gradient;
 
    DataFilters filters;
 
    bool modifiable = false;
 
    bool modified = false;
 
    LabelType label = L_NONE;
 
    int peptide_id_index = -1;
    int peptide_hit_index = -1;
 
  private:
    std::string name_;
    std::string name_suffix_;
  };
 
  class LayerAnnotatorBase
  {
  public:
    LayerAnnotatorBase(const FileTypeList& supported_types, const std::string& file_dialog_text, QWidget* gui_lock);
    
    virtual ~LayerAnnotatorBase() = default;
 
    bool annotateWithFileDialog(LayerDataBase& layer, LogWindow& log, const std::string& current_path) const;
 
    bool annotateWithFilename(LayerDataBase& layer, LogWindow& log, const std::string& filename) const;
 
    static std::unique_ptr<LayerAnnotatorBase> getAnnotatorWhichSupports(const FileTypes::Type& type);
 
    static std::unique_ptr<LayerAnnotatorBase> getAnnotatorWhichSupports(const std::string& filename);
 
  protected:
    virtual bool annotateWorker_(LayerDataBase& layer, const std::string& filename, LogWindow& log) const = 0;
 
    const FileTypeList supported_types_;
    const std::string file_dialog_text_;
    QWidget* gui_lock_ = nullptr;
  };
 
  class LayerAnnotatorPeptideID : public LayerAnnotatorBase
  {
  public:
    LayerAnnotatorPeptideID(QWidget* gui_lock) :
        LayerAnnotatorBase(std::vector<FileTypes::Type>{FileTypes::IDXML, FileTypes::MZIDENTML},
                           "Select peptide identification data", gui_lock)
    {
    }
 
  protected:
    virtual bool annotateWorker_(LayerDataBase& layer, const std::string& filename, LogWindow& log) const;
  };
 
  class LayerAnnotatorAMS : public LayerAnnotatorBase
  {
  public:
    LayerAnnotatorAMS(QWidget* gui_lock) :
        LayerAnnotatorBase(std::vector<FileTypes::Type>{FileTypes::FEATUREXML},
                           "Select AccurateMassSearch's featureXML file", gui_lock)
    {
    }
 
  protected:
    virtual bool annotateWorker_(LayerDataBase& layer, const std::string& filename, LogWindow& log) const;
  };
 
  class LayerAnnotatorOSW : public LayerAnnotatorBase
  {
  public:
    LayerAnnotatorOSW(QWidget* gui_lock) :
        LayerAnnotatorBase(std::vector<FileTypes::Type>{FileTypes::OSW},
                           "Select OpenSwath/pyProphet output file", gui_lock)
    {
    }
 
  protected:
    virtual bool annotateWorker_(LayerDataBase& layer, const std::string& filename, LogWindow& log) const;
  };
 
  OPENMS_GUI_DLLAPI std::ostream& operator<<(std::ostream& os, const LayerDataBase& rhs);
 
}// namespace OpenMS