File: atomCount.cpp

package info (click to toggle)
massxpert 2.3.6-1squeeze1
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 20,736 kB
  • ctags: 3,541
  • sloc: cpp: 44,108; xml: 7,381; sh: 604; makefile: 108; ansic: 7
file content (220 lines) | stat: -rw-r--r-- 5,686 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
/* massXpert - the true massist's program.
   -------------------------------------- 
   Copyright(C) 2006,2007 Filippo Rusconi

   http://www.massxpert.org/massXpert

   This file is part of the massXpert project.
   
   The massxpert project is the successor to the "GNU polyxmass"
   project that is an official GNU project package(see
   www.gnu.org). The massXpert project is not endorsed by the GNU
   project, although it is released ---in its entirety--- under the
   GNU General Public License. A huge part of the code in massXpert
   is actually a C++ rewrite of code in GNU polyxmass. As such
   massXpert was started at the Centre National de la Recherche
   Scientifique(FRANCE), that granted me the formal authorization to
   publish it under this Free Software License.

   This software is free software; you can redistribute it and/or
   modify it under the terms of the GNU  General Public
   License version 3, as published by the Free Software Foundation.
   
   This software is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
   General Public License for more details.
   
   You should have received a copy of the GNU General Public License
   along with this software; if not, write to the 

   Free Software Foundation, Inc.,
   
   51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA.
*/


/////////////////////// Local includes
#include "atomCount.hpp"


namespace massXpert
{

  //! Constructs an atom-count. Initializes the count to 0.
  AtomCount::AtomCount()
  {
    m_count = 0;
  }


  //! Constructs a copy of \p other.
  /*!  \param other atom-count to be used as a mold.
   */
  AtomCount::AtomCount(const AtomCount &other)
    : Atom(other), m_count(other.m_count)
  {
  }


  //! Creates a new atom count initialized with \c this.
  /*! The initialization of the new atom count involves duplicating all
    the data of \c this, including all the isotopes of the isotope list.

    \return The newly created atom count, which should be deleted when no
    longer in use.
  */
  AtomCount *
  AtomCount::clone() const
  {
    AtomCount *other = new AtomCount(*this);

    return other;
  }


  //! Modifies \p other to be identical to \p this.
  /*!  \param other atom-count.
   */
  void
  AtomCount::clone(AtomCount *other) const
  {
    Q_ASSERT(other);
  
    if (other == this)
      return;

    other->m_count = m_count;

    Atom::clone(other);
  }


  //! Modifies \p this  to be identical to \p other.
  /*!  \param other atom-count to be used as a mold.
   */
  void
  AtomCount::mold(const AtomCount &other)
  {
    if (&other == this)
      return;

    m_count = other.m_count;
  
    Atom::mold(other);
  }


  //! Assigns other to \p this atom-count.
  /*! \param other atom-count used as the mold to set values to \p this
    instance.
  
    \return a reference to \p this atom-count.
  */
  AtomCount & 
  AtomCount::operator =(const AtomCount &other)
  {
    if (&other != this)
      mold(other);
  
    return *this;
  }


  //! Sets the count of atom.
  /*!  \param count new count.
   */
  void
  AtomCount::setCount(int count)
  {
    m_count = count;
  }


  //! Increments the count by \p times.
  /*!  \param times value by which the count is incremented. Defaults to
    1.
  */
  void
  AtomCount::incrementCount(int times)
  {
    m_count += times;
  }


  //! Decrements the count by \p times.
  /*!  \param times value by which the count is incremented. Defaults to
    1.
  */
  void 
  AtomCount::decrementCount(int times)
  {
    m_count -= times;
  }


  //! Returns the count.
  /*! 
    \return the count.
  */
  int
  AtomCount::count() const
  {
    return m_count;
  }


  //! Increments the masses in the arguments.
  /*! This function performs the same that the base class function does
    unless the computation involves the atom count member data. The
    masses for \p this atom-count are not found in \p this instance,
    actually. First a reference atom is searched in \p reflist using \p
    this atom-count symbol as search criterion. The monoisotopic and
    average masses of the found atom are used for the computation. The
    arguments \p mono and \p avg are updated by incrementation of their
    values using the masses(monoisotopic and average) compounded by the
    \p times argument and by the count.
  
    For example, if \p times is 2 ; \p *mono is 100 and \p *avg is 101 ;
    \p this atom's monoisotopic mass is 200 and average mass is 202 ;
    the atom occurence count is 3, then the computation leads to \p
    *mono = 100 + 2 * 3 * 200 and \p *avg = 101 + 2 * 3 * 202.

    \param refList list of reference atoms.

    \param mono monoisotopic mass to update. Defaults to 0, in which
    case this mass is not updated.

    \param avg average mass to update. Defaults to 0, in which case this
    mass is not updated.

    \param times times that the increment should be performed. Defaults
    to 1.

    \return true if successful, false otherwise(that is, if no
    reference atom could be found in \p refList).

    \sa Atom::accountMasses(const QList<Atom *> &refList, double
    *mono, double *avg, int times).
    */
  bool 
  AtomCount::accountMasses(const QList<Atom *> &refList,
			    double *mono, double *avg, int times) const
  {
    int idx = -1;
  
    idx = isSymbolKnown(refList);
  
    if (idx == -1)
      return false;
  
    if (mono)
      *mono += refList.at(idx)->mono() * m_count * times;
  
    if (avg)
      *avg += refList.at(idx)->avg() * m_count * times;
  
    return true;
  }

} // namespace massXpert