OpenMS  2.7.0
MassExplainer.h
Go to the documentation of this file.
1 // --------------------------------------------------------------------------
2 // OpenMS -- Open-Source Mass Spectrometry
3 // --------------------------------------------------------------------------
4 // Copyright The OpenMS Team -- Eberhard Karls University Tuebingen,
5 // ETH Zurich, and Freie Universitaet Berlin 2002-2021.
6 //
7 // This software is released under a three-clause BSD license:
8 // * Redistributions of source code must retain the above copyright
9 // notice, this list of conditions and the following disclaimer.
10 // * Redistributions in binary form must reproduce the above copyright
11 // notice, this list of conditions and the following disclaimer in the
12 // documentation and/or other materials provided with the distribution.
13 // * Neither the name of any author or any participating institution
14 // may be used to endorse or promote products derived from this software
15 // without specific prior written permission.
16 // For a full list of authors, refer to the file AUTHORS.
17 // --------------------------------------------------------------------------
18 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
19 // AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
20 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
21 // ARE DISCLAIMED. IN NO EVENT SHALL ANY OF THE AUTHORS OR THE CONTRIBUTING
22 // INSTITUTIONS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
23 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
24 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
25 // OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
26 // WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
27 // OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
28 // ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
29 //
30 // --------------------------------------------------------------------------
31 // $Maintainer: Chris Bielow $
32 // $Authors: Chris Bielow $
33 // --------------------------------------------------------------------------
34 
35 #pragma once
36 
38 
39 #include <OpenMS/CONCEPT/Types.h>
41 
43 #include <OpenMS/OpenMSConfig.h>
44 
45 namespace OpenMS
46 {
47  class Compomer;
48 
55  class OPENMS_DLLAPI MassExplainer
56  {
57 
58 public:
59 
60  typedef Adduct::AdductsType AdductsType; //vector<Adduct>
61  typedef std::vector<Compomer>::const_iterator CompomerIterator;
62 
64 
65  MassExplainer();
67 
69  MassExplainer(AdductsType adduct_base);
70 
72  MassExplainer(Int q_min, Int q_max, Int max_span, double thresh_logp);
73 
75  MassExplainer(AdductsType adduct_base, Int q_min, Int q_max, Int max_span, double thresh_logp, Size max_neutrals);
76 
77 
78 private:
81  void init_(bool init_thresh_p);
82 public:
85 
87  virtual ~MassExplainer();
89 
90 
92  void compute();
93 
94 
95  //@name Accessors
97 
99  void setAdductBase(AdductsType adduct_base);
102 
104  const Compomer& getCompomerById(Size id) const;
106 
107 
115  SignedSize query(const Int net_charge,
116  const float mass_to_explain,
117  const float mass_delta,
118  const float thresh_log_p,
119  std::vector<Compomer>::const_iterator& firstExplanation,
120  std::vector<Compomer>::const_iterator& lastExplanation) const;
121 protected:
122 
124  bool compomerValid_(const Compomer& cmp);
125 
127  Adduct createAdduct_(const String& formula, const Int charge, const double p) const;
128 
130  std::vector<Compomer> explanations_;
140  double thresh_p_;
143 
144  };
145 
146 
147 } // namespace OpenMS
148 
Definition: Adduct.h:45
std::vector< Adduct > AdductsType
Definition: Adduct.h:48
Holds information on an edge connecting two features from a (putative) charge ladder.
Definition: Compomer.h:59
computes empirical formulas for given mass differences using a set of allowed elements
Definition: MassExplainer.h:56
MassExplainer(AdductsType adduct_base)
Constructor.
Adduct createAdduct_(const String &formula, const Int charge, const double p) const
create a proper adduct from formula and charge and probability
Int q_max_
maximal expected charge
Definition: MassExplainer.h:136
MassExplainer & operator=(const MassExplainer &rhs)
Assignment operator.
void compute()
fill map with possible mass-differences along with their explanation
Int q_min_
minimal expected charge
Definition: MassExplainer.h:134
SignedSize query(const Int net_charge, const float mass_to_explain, const float mass_delta, const float thresh_log_p, std::vector< Compomer >::const_iterator &firstExplanation, std::vector< Compomer >::const_iterator &lastExplanation) const
const Compomer & getCompomerById(Size id) const
return a compomer by its Id (useful after a query() ).
AdductsType getAdductBase() const
Returns the set of adducts.
AdductsType adduct_base_
all allowed adducts, whose combination explains the mass difference
Definition: MassExplainer.h:132
void init_(bool init_thresh_p)
bool compomerValid_(const Compomer &cmp)
check if the generated compomer is valid judged by its probability, charges etc
MassExplainer(Int q_min, Int q_max, Int max_span, double thresh_logp)
Constructor.
MassExplainer(AdductsType adduct_base, Int q_min, Int q_max, Int max_span, double thresh_logp, Size max_neutrals)
Constructor.
void setAdductBase(AdductsType adduct_base)
Sets the set of possible adducts.
Size max_neutrals_
Maximum number of neutral(q=0) adducts.
Definition: MassExplainer.h:142
double thresh_p_
minimum required probability of a compound (all other compounds are discarded)
Definition: MassExplainer.h:140
Int max_span_
maximal span (in terms of charge) for co-features, e.g. a cluster with q={3,6} has span=4
Definition: MassExplainer.h:138
virtual ~MassExplainer()
Destructor.
std::vector< Compomer > explanations_
store possible explanations (as formula) for a certain ChargeDifference and MassDifference
Definition: MassExplainer.h:130
Adduct::AdductsType AdductsType
Definition: MassExplainer.h:60
std::vector< Compomer >::const_iterator CompomerIterator
Definition: MassExplainer.h:61
A more convenient string class.
Definition: String.h:61
int Int
Signed integer type.
Definition: Types.h:102
ptrdiff_t SignedSize
Signed Size type e.g. used as pointer difference.
Definition: Types.h:134
size_t Size
Size type e.g. used as variable which can hold result of size()
Definition: Types.h:127
Main OpenMS namespace.
Definition: FeatureDeconvolution.h:47