EmpiricalFormula.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: Chris Bielow, Ahmed Khalil $
// $Authors: Andreas Bertsch, Chris Bielow $
// --------------------------------------------------------------------------
//
#pragma once
 
#include <iosfwd>
#include <map>
#include <set>
#include <string>
 
#include <OpenMS/CONCEPT/Types.h>
 
namespace OpenMS
{
  class String;
  class Element;
  class ElementDB;
  class IsotopeDistribution;
  class IsotopePatternGenerator;
  class CoarseIsotopePatternGenerator;
 
  class OPENMS_DLLAPI EmpiricalFormula
  {
 
protected:
    typedef std::map<const Element*, SignedSize> MapType_;
 
public:
    typedef MapType_::const_iterator ConstIterator;
    typedef MapType_::const_iterator const_iterator;
    typedef MapType_::iterator Iterator;
    typedef MapType_::iterator iterator;
 
    EmpiricalFormula();
 
    EmpiricalFormula(const EmpiricalFormula&) = default;
 
    EmpiricalFormula(EmpiricalFormula&&) = default;
 
    explicit EmpiricalFormula(const String& rhs);
 
    EmpiricalFormula(SignedSize number, const Element* element, SignedSize charge = 0);
 
    virtual ~EmpiricalFormula();
 
    static EmpiricalFormula fromString(const String& rhs)
    {
      EmpiricalFormula ef(rhs);
      return ef;
    }
 
    double getMonoWeight() const;
 
    double getLightestIsotopeWeight() const;
 
    double getAverageWeight() const;
 
    double calculateTheoreticalIsotopesNumber() const;
 
    bool estimateFromWeightAndComp(double average_weight, double C, double H, double N, double O, double S, double P);
 
    bool estimateFromMonoWeightAndComp(double mono_weight, double C, double H, double N, double O, double S, double P);
 
    bool estimateFromWeightAndCompAndS(double average_weight, UInt S, double C, double H, double N, double O, double P);
 
 
    IsotopeDistribution getIsotopeDistribution(const IsotopePatternGenerator& method) const;
 
    IsotopeDistribution getConditionalFragmentIsotopeDist(const EmpiricalFormula& precursor,
                                                          const std::set<UInt>& precursor_isotopes,
                                                          const CoarseIsotopePatternGenerator& method) const;
 
    SignedSize getNumberOf(const Element* element) const;
 
    SignedSize getNumberOfAtoms() const;
 
    Int getCharge() const;
 
    void setCharge(Int charge);
 
    String toString() const;
 
    std::map<std::string, int> toMap() const;
 
 
    EmpiricalFormula& operator=(const EmpiricalFormula&) = default;
 
    EmpiricalFormula& operator=(EmpiricalFormula&&) & = default;
 
    EmpiricalFormula& operator+=(const EmpiricalFormula& rhs);
 
    EmpiricalFormula operator*(const SignedSize& times) const;
 
    EmpiricalFormula operator+(const EmpiricalFormula& rhs) const;
 
    EmpiricalFormula& operator-=(const EmpiricalFormula& rhs);
 
    EmpiricalFormula operator-(const EmpiricalFormula& rhs) const;
 
 
    bool isEmpty() const;
 
    bool isCharged() const;
 
    bool hasElement(const Element* element) const;
 
    bool contains(const EmpiricalFormula& ef) const;
 
    bool operator==(const EmpiricalFormula& rhs) const;
 
    bool operator!=(const EmpiricalFormula& rhs) const;
 
    bool operator<(const EmpiricalFormula& rhs) const;
 
 
    friend OPENMS_DLLAPI std::ostream& operator<<(std::ostream& os, const EmpiricalFormula& formula);
 
    inline ConstIterator begin() const { return formula_.begin(); }
 
    inline ConstIterator end() const { return formula_.end(); }
 
    inline Iterator begin() { return formula_.begin(); }
 
    inline Iterator end() { return formula_.end(); }
 
    // @TODO: make these static member variables instead?
    static EmpiricalFormula hydrogen(int n_atoms = 1);
 
    static EmpiricalFormula water(int n_molecules = 1);
 
protected:
 
    void removeZeroedElements_();
 
    MapType_ formula_;
 
    Int charge_;
 
    Int parseFormula_(std::map<const Element*, SignedSize>& ef, const String& formula) const;
 
  };
 
  OPENMS_DLLAPI std::ostream& operator<<(std::ostream& os, const EmpiricalFormula& formula);
 
} // namespace OpenMS