BinnedSpectrum.h

Go to the documentation of this file.

// Copyright (c) 2002-2023, The OpenMS Team -- EKU Tuebingen, ETH Zurich, and FU Berlin
// SPDX-License-Identifier: BSD-3-Clause
//
// --------------------------------------------------------------------------
// $Maintainer: Timo Sachsenberg $
// $Authors: Timo Sachsenberg, Mathias Walzer $
// --------------------------------------------------------------------------
//
#pragma once
 
#include <OpenMS/KERNEL/StandardTypes.h>
#include <OpenMS/KERNEL/MSSpectrum.h>
#include <OpenMS/CONCEPT/Exception.h>
 
#include <cmath>
 
// forward decl
namespace Eigen
{
    template<typename _Scalar, int _Flags, typename _StorageIndex>
    class SparseVector;
}
 
namespace OpenMS
{
 
  class OPENMS_DLLAPI BinnedSpectrum
  {
    // smallest possible m/z value (needs to be >= 1)
    static constexpr const float MIN_MZ_ = 1.0;
 
public:
    // default bin width for low-resolution data (adapted from doi:10.1007/s13361-015-1179-x)
    static constexpr const float DEFAULT_BIN_WIDTH_LOWRES = 1.0005f;
 
    // default bin width for high-resolution data (adapted from doi:10.1007/s13361-015-1179-x)
    static constexpr const float DEFAULT_BIN_WIDTH_HIRES = 0.02f;
 
    static constexpr const float DEFAULT_BIN_OFFSET_HIRES = 0.0f;
 
    static constexpr const float DEFAULT_BIN_OFFSET_LOWRES = 0.4f;
 
    using SparseVectorType = Eigen::SparseVector<float, 0, int>;
 
    // static const SparseVectorType EmptySparseVector;
 
    // BinnedSpectrum() = delete;
    BinnedSpectrum() {};
 
    BinnedSpectrum(const PeakSpectrum& ps, float size, bool unit_ppm, UInt spread, float offset);
 
    BinnedSpectrum(const BinnedSpectrum&);
 
    virtual ~BinnedSpectrum();
 
    BinnedSpectrum& operator=(const BinnedSpectrum&);
 
    bool operator==(const BinnedSpectrum& rhs) const;
 
    bool operator!=(const BinnedSpectrum& rhs) const;
 
    float getBinIntensity(double mz);
 
    size_t getBinIndex(float mz) const;
 
    inline float getBinLowerMZ(size_t i) const
    {
      if (unit_ppm_)
      {
        // mz = MIN_MZ_ * (1.0 + bin_size_)^index for index
        return float(MIN_MZ_ * pow(1.0 + bin_size_ * 1e-6, i));
      }
      else 
      { 
        return ((static_cast<float>(i) - offset_) * bin_size_);
      }
    }
 
    inline float getBinSize() const { return bin_size_; }
 
    inline size_t getBinSpread() const { return bin_spread_; }
 
    const SparseVectorType* getBins() const;
 
    SparseVectorType* getBins();
 
    inline float getOffset() const { return offset_; }
 
    const std::vector<Precursor>& getPrecursors() const;
 
    std::vector<Precursor>& getPrecursors();
 
    //  returns true if bin size, unit and offset are equal, otherwise false
    static bool isCompatible(const BinnedSpectrum& a, const BinnedSpectrum& b);
 
private:
    UInt bin_spread_ {0};
 
    float bin_size_ {0};
 
    bool unit_ppm_ {false};
 
    float offset_ {0};
 
    SparseVectorType* bins_ {nullptr};
 
    void binSpectrum_(const PeakSpectrum& ps);
 
    std::vector<Precursor> precursors_;
  };
 
}