OpenMS  2.4.0
IsotopeWavelet.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-2018.
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: Timo Sachsenberg$
32 // $Authors: $
33 // --------------------------------------------------------------------------
34 
35 #pragma once
36 
40 
41 namespace OpenMS
42 {
50  class OPENMS_DLLAPI IsotopeWavelet
51  {
52 public:
53 
55  static IsotopeWavelet * init(const double max_m, const UInt max_charge);
56 
59  {
60  return me_;
61  }
62 
64  static void destroy();
65 
66 
76  static double getValueByMass(const double t, const double m, const UInt z, const Int mode = +1)
77  {
78  return getValueByLambda(getLambdaL(m * z - z * mode * Constants::IW_PROTON_MASS), t * z + 1);
79  }
80 
89  static double getValueByLambda(const double lambda, const double tz1);
90 
103  static double getValueByLambdaExtrapol(const double lambda, const double tz1);
104 
105  static double getValueByLambdaExact(const double lambda, const double tz1);
106 
107 
110  {
111  return max_charge_;
112  }
113 
115  static void setMaxCharge(const UInt max_charge)
116  {
117  max_charge_ = max_charge;
118  }
119 
124  static double getTableSteps()
125  {
126  return table_steps_;
127  }
128 
133  static double getInvTableSteps()
134  {
135  return inv_table_steps_;
136  }
137 
139  static void setTableSteps(const double table_steps)
140  {
141  inv_table_steps_ = 1. / table_steps;
142  table_steps_ = table_steps;
143  }
144 
146  static double getLambdaL(const double m);
147 
148 
153  static const IsotopeDistribution::ContainerType & getAveragine(const double m, UInt * size = nullptr);
154 
155 
158  {
159  return gamma_table_max_index_;
160  }
161 
164  {
165  return exp_table_max_index_;
166  }
167 
170  static float myPow(float a, float b);
171 
172  static UInt getMzPeakCutOffAtMonoPos(const double mass, const UInt z);
173 
174  static UInt getNumPeakCutOff(const double mass, const UInt z);
175 
176  static UInt getNumPeakCutOff(const double mz);
177 
178 
179 protected:
180 
182  static IsotopeWavelet * me_;
183 
185  IsotopeWavelet();
186 
190  IsotopeWavelet(const double max_m, const UInt max_charge);
191 
192 
194  virtual ~IsotopeWavelet();
195 
196 
205  static void preComputeExpensiveFunctions_(const double max_m);
206 
207 
210  static void computeIsotopeDistributionSize_(const double max_m);
211 
212 
217  static float myPow2_(float i);
218 
223  static float myLog2_(float i);
224 
226  union fi_
227  {
229  float f;
230  };
231 
234 
236  static double table_steps_;
237  static double inv_table_steps_;
238 
240  static std::vector<double> gamma_table_;
241  static std::vector<double> gamma_table_new_;
242 
244  static std::vector<double> exp_table_;
245 
247  static std::vector<double> sine_table_;
248 
252 
255 
256  };
257 
258 } //namespace
259 
float f
Definition: IsotopeWavelet.h:229
Isotope pattern generator for coarse isotope distributions.
Definition: CoarseIsotopePatternGenerator.h:88
static std::vector< double > gamma_table_
Definition: IsotopeWavelet.h:240
static Size gamma_table_max_index_
Definition: IsotopeWavelet.h:253
const double IW_PROTON_MASS
Definition: IsotopeWaveletConstants.h:68
static std::vector< double > exp_table_
Definition: IsotopeWavelet.h:244
static void setMaxCharge(const UInt max_charge)
Sets the max_charge parameter.
Definition: IsotopeWavelet.h:115
unsigned int UInt
Unsigned integer type.
Definition: Types.h:94
Definition: IsotopeDistribution.h:72
Int i
Definition: IsotopeWavelet.h:228
std::vector< MassAbundance > ContainerType
Definition: IsotopeDistribution.h:80
static std::vector< double > gamma_table_new_
Definition: IsotopeWavelet.h:241
Main OpenMS namespace.
Definition: FeatureDeconvolution.h:46
static double inv_table_steps_
Definition: IsotopeWavelet.h:237
Implements the isotope wavelet function.
Definition: IsotopeWavelet.h:50
static IsotopeWavelet * me_
Definition: IsotopeWavelet.h:182
static CoarseIsotopePatternGenerator solver_
Definition: IsotopeWavelet.h:250
static Size getExpTableMaxIndex()
Returns the largest possible index for the pre-sampled exp table.
Definition: IsotopeWavelet.h:163
static UInt getMaxCharge()
Returns the largest charge state we will consider.
Definition: IsotopeWavelet.h:109
static IsotopeWavelet * getInstance()
Definition: IsotopeWavelet.h:58
static void setTableSteps(const double table_steps)
Sets the table_steps parameter.
Definition: IsotopeWavelet.h:139
static double table_steps_
Definition: IsotopeWavelet.h:236
static Size exp_table_max_index_
Definition: IsotopeWavelet.h:254
size_t Size
Size type e.g. used as variable which can hold result of size()
Definition: Types.h:127
static Size getGammaTableMaxIndex()
Returns the largest possible index for the pre-sampled gamma table.
Definition: IsotopeWavelet.h:157
Internal union for fast computation of the power function.
Definition: IsotopeWavelet.h:226
static UInt max_charge_
Definition: IsotopeWavelet.h:233
static double getInvTableSteps()
Returns the inv_table_steps_ parameter.
Definition: IsotopeWavelet.h:133
static double getTableSteps()
Returns the table_steps_ parameter.
Definition: IsotopeWavelet.h:124
int Int
Signed integer type.
Definition: Types.h:102
static double getValueByMass(const double t, const double m, const UInt z, const Int mode=+1)
Returns the value of the isotope wavelet at position t. Usually, you do not need to call this functio...
Definition: IsotopeWavelet.h:76
static IsotopeDistribution averagine_
Definition: IsotopeWavelet.h:251
static std::vector< double > sine_table_
Definition: IsotopeWavelet.h:247