libpappsomspp
Library for mass spectrometry
massspectrumwidget.cpp
Go to the documentation of this file.
1 /**
2  * \file pappsomspp/widget/spectrumwidget/massspectrumwidget.cpp
3  * \date 22/12/2017
4  * \author Olivier Langella
5  * \brief plot a sectrum and annotate with peptide
6  */
7 
8 
9 /*******************************************************************************
10  * Copyright (c) 2017 Olivier Langella <Olivier.Langella@u-psud.fr>.
11  *
12  * This file is part of the PAPPSOms++ library.
13  *
14  * PAPPSOms++ is free software: you can redistribute it and/or modify
15  * it under the terms of the GNU General Public License as published by
16  * the Free Software Foundation, either version 3 of the License, or
17  * (at your option) any later version.
18  *
19  * PAPPSOms++ is distributed in the hope that it will be useful,
20  * but WITHOUT ANY WARRANTY; without even the implied warranty of
21  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
22  * GNU General Public License for more details.
23  *
24  * You should have received a copy of the GNU General Public License
25  * along with PAPPSOms++. If not, see <http://www.gnu.org/licenses/>.
26  *
27  * Contributors:
28  * Olivier Langella <Olivier.Langella@u-psud.fr> - initial API and
29  *implementation
30  ******************************************************************************/
31 #include "massspectrumwidget.h"
32 #include "../../pappsoexception.h"
33 #include "../../peptide/peptidenaturalisotopelist.h"
34 #include <QDebug>
35 #include <QWidget>
36 
37 
38 using namespace pappso;
39 
41  : GraphicDeviceWidget(parent)
42 {
43  qDebug();
44 
45  _ms_level = 1;
47  _custom_plot = nullptr;
48 
49  this->setLayout(new QHBoxLayout(this));
50 
51  this->layout()->setContentsMargins(0,0,0,0);
52  setVisibleMassDelta(false);
53  qDebug();
54 }
56 {
57 }
58 
59 bool
60 MassSpectrumWidget::savePdf(const QString &fileName, int width, int height)
61 {
62 
63  if(_custom_plot != nullptr)
64  {
65  return _custom_plot->savePdf(fileName, width, height);
66  }
67  else
68  {
69  return false;
70  }
71 }
72 
73 
74 void
75 MassSpectrumWidget::toQPaintDevice(QPaintDevice *device, const QSize &size)
76 {
77 
78  if(_custom_plot != nullptr)
79  {
80  QCPPainter painter;
81  painter.begin(device);
82  _custom_plot->toPainter(&painter, size.width(), size.height());
83  painter.end();
84  }
85 }
86 void
88 {
89  qDebug();
90  if(_custom_plot != nullptr)
91  {
92  if(visible == _is_visible_mass_delta)
93  return;
94  delete _custom_plot;
95  }
96  _is_visible_mass_delta = visible;
97  while(auto item = this->layout()->takeAt(0))
98  {
99  delete item->widget();
100  }
101  qDebug();
102  _custom_plot = new QCPSpectrum(this, visible);
103  this->layout()->addWidget(_custom_plot);
104  qDebug();
105  _custom_plot->xAxis->setLabel("m/z");
106  _custom_plot->yAxis->setLabel("intensity");
107  qDebug();
108  _custom_plot->setInteractions(QCP::iRangeDrag | QCP::iRangeZoom);
109  _custom_plot->axisRect()->setRangeDrag(Qt::Horizontal);
110  _custom_plot->axisRect()->setRangeZoom(Qt::Horizontal);
111  qDebug();
112  // legend->setVisible(false);
113  // legend->setFont(QFont("Helvetica",9));
114  // set locale to english, so we get english decimal separator:
115  // setLocale(QLocale(QLocale::English, QLocale::UnitedKingdom));
116  qDebug();
117 }
118 void
120 {
121  qDebug();
123  qDebug();
124  _custom_plot->clearItems();
126  qDebug();
127 }
128 
129 void
130 MassSpectrumWidget::setPeptideCharge(unsigned int parent_ion_charge)
131 {
132  _peptide_charge = parent_ion_charge;
133 }
134 void
135 MassSpectrumWidget::setIonList(const std::list<PeptideIon> &ion_list)
136 {
137  _ion_list = ion_list;
138 }
139 
140 void
141 MassSpectrumWidget::setMsLevel(unsigned int ms_level)
142 {
143  qDebug() << "ms_level=" << ms_level;
144  _ms_level = ms_level;
145 
146 
147  if(_ms_level == 1)
148  {
149  setVisibleMassDelta(false);
150  }
151  else
152  {
153  setVisibleMassDelta(true);
154  }
155 
156  //_precision._precision = precision._precision;
157 }
158 void
160 {
161  qDebug() << "precision->toString()=" << precision->toString();
162  _p_ms1_precision = precision;
163  //_precision._precision = precision._precision;
164 }
165 void
167 {
168  qDebug() << "precision->toString()=" << precision->toString();
169  _p_ms2_precision = precision;
170  //_precision._precision = precision._precision;
171 }
172 
173 void
174 MassSpectrumWidget::setMaximumIsotopeNumber(unsigned int max_isotope_number)
175 {
176  _max_isotope_number = max_isotope_number;
177 }
178 
179 void
180 MassSpectrumWidget::setMaximumIsotopeRank(unsigned int max_isotope_rank)
181 {
182  _max_isotope_rank = max_isotope_rank;
183 }
184 void
186 {
187  qDebug() << " _max_isotope_number=" << _max_isotope_number;
188  clearData();
190  if((_spectrum_sp == nullptr) || (_peptide_sp == nullptr))
191  {
192  }
193  else
194  {
195  if(_ms_level > 1)
196  {
197  PeptideIsotopeSpectrumMatch psm_match(*(_spectrum_sp.get()),
198  _peptide_sp,
201  _ion_list,
204 
206  }
207  else
208  {
209  }
211  }
212  qDebug();
213 }
214 
215 void
217 {
218  qDebug();
219  _peptide_sp = peptide_sp;
220 
221  // clearData();
222  qDebug();
223 }
224 
225 void
227 {
228  qDebug();
229  _spectrum_sp = spectrum;
230 
231  clearData();
232  qDebug();
233 }
234 
235 void
237 {
238  qDebug();
239 
241 
242  /*
243  if (_p_delta_axis_rect != nullptr) {
244  _p_delta_axis_rect->axis(QCPAxis::AxisType::atLeft)->rescale();
245  }
246  */
247  _custom_plot->replot();
248  qDebug();
249 }
250 
251 void
253  const QualifiedMassSpectrum &spectrum)
254 {
255  qDebug() << "spectrum.getPrecursorCharge()=" << spectrum.getPrecursorCharge();
256 
257  setMsLevel(spectrum.getMsLevel());
259 
260  qDebug();
261 }
262 
263 void
265 {
266 
267  qDebug();
268  peptideAnnotate();
269  if(_ms_level == 1)
270  {
271  if(_spectrum_sp != nullptr)
272  {
273  if(_isotope_mass_list.size() > 0)
274  {
275 
276  qDebug() << "_isotope_mass_list.size()="
277  << _isotope_mass_list.size();
278  std::sort(_isotope_mass_list.begin(),
279  _isotope_mass_list.end(),
282  return a.get()->getMz() < b.get()->getMz();
283  });
284 
285  if(_isotope_mass_list.size() > 0)
286  {
287  PeptideNaturalIsotopeAverageSp precursor_peptide =
288  _isotope_mass_list.at(0);
289  qDebug() << "precursor_peptide.get()->getMz()="
290  << precursor_peptide.get()->getMz();
291  MzRange precursor_mass(precursor_peptide.get()->getMz(),
293  DataPoint monoisotope_peak;
294  monoisotope_peak.y = 0;
295 
296  for(const DataPoint &peak : *(_spectrum_sp.get()))
297  {
298  if(precursor_mass.contains(peak.x))
299  {
300  if(peak.y > monoisotope_peak.y)
301  {
302  qDebug() << "SpectrumWidget::plot "
303  "(peak.intensity > "
304  "monoisotope_peak.intensity) ";
305  monoisotope_peak = peak;
306  }
307  }
308  }
309  if(monoisotope_peak.y > 0)
310  {
311  qDebug() << "addMs1IsotopePattern";
313  monoisotope_peak.y);
314  }
315  }
316  }
317  }
318  }
319  else
320  {
322  [](const PeakIonIsotopeMatch &a, const PeakIonIsotopeMatch &b) {
323  return a.getPeak().y > b.getPeak().y;
324  });
325  unsigned int i = 0;
326  for(const PeakIonIsotopeMatch &peak_ion_match :
328  {
329  _custom_plot->addPeakIonIsotopeMatch(peak_ion_match);
330 
331  _custom_plot->addMassDelta(peak_ion_match);
332  //_p_delta_graph->addData(peak_ion_match.getPeak().x,
333  // peak_ion_match.getPeak().y);
334  if(i < _tag_nmost_intense)
335  {
336  QCPItemText *text_label = new QCPItemText(_custom_plot);
337  text_label->setVisible(true);
338  //_custom_plot->addItem(text_label);
339  text_label->setPositionAlignment(Qt::AlignBottom |
340  Qt::AlignHCenter);
341  text_label->position->setType(QCPItemPosition::ptPlotCoords);
342  text_label->position->setCoords(
343  peak_ion_match.getPeak().x,
344  peak_ion_match.getPeak()
345  .y); // place position at center/top of axis rect
346  text_label->setFont(QFont(font().family(), 8));
347  text_label->setText(
348  peak_ion_match.getPeptideFragmentIonSp()
349  .get()
350  ->getCompletePeptideIonName(peak_ion_match.getCharge()));
351  // text_label->setPen(QPen(PeptideFragmentIon::getPeptideIonColor(peak_ion_match.getPeptideIonType()),
352  // 1)); // show black border around text
353  text_label->setColor(
355  peak_ion_match.getPeptideIonType())));
356  }
357  i++;
358  }
359  }
360 
361  _custom_plot->replot();
362  qDebug();
363 }
364 void
366 {
367  emit mzChanged(mz);
368 }
369 
370 void
372 {
373  qDebug() << "p_peak_match=" << p_peak_match;
374  if(_p_mouse_peak != p_peak_match)
375  {
376  _p_mouse_peak = p_peak_match;
377  DataPointCstSPtr peak_shp;
378  // emit peakChanged(peak_shp);
379  if(_p_mouse_peak != nullptr)
380  {
381  qDebug() << "_p_mouse_peak->x=" << _p_mouse_peak->x;
382  peak_shp = _p_mouse_peak->makeDataPointCstSPtr();
383  emit peakChanged(peak_shp);
384  // try to find matched ion (if it exists)
385  bool found = false;
386  for(const PeakIonIsotopeMatch &peak_ion_match :
388  {
389  if(peak_ion_match.getPeak().x == _p_mouse_peak->x)
390  {
391  // found
392  emit ionChanged(std::make_shared<const PeakIonIsotopeMatch>(
393  PeakIonIsotopeMatch(peak_ion_match)));
394  found = true;
395  }
396  }
397  if(!found)
398  {
399  emit ionChanged(std::shared_ptr<const PeakIonIsotopeMatch>());
400  }
401  }
402  else
403  {
404  qDebug() << "no peak";
405  emit peakChanged(peak_shp);
406  }
407  }
408 }
409 
410 void
412 {
413  qDebug() << "_p_ms1_precision->toString()=" << _p_ms1_precision->toString();
414  _isotope_mass_list.resize(0);
415  // compute isotope masses :
416  if(_peptide_sp != nullptr)
417  {
419  std::map<unsigned int, pappso::pappso_double> map_isotope_number =
420  isotope_list.getIntensityRatioPerIsotopeNumber();
421 
422  for(unsigned int i = 0; i < map_isotope_number.size(); i++)
423  {
424 
425  unsigned int asked_rank = 0;
426  unsigned int given_rank = 0;
427  bool more_rank = true;
428  while(more_rank)
429  {
430  asked_rank++;
431  pappso::PeptideNaturalIsotopeAverage isotopeAverageMono(
432  isotope_list, asked_rank, i, _peptide_charge, _p_ms1_precision);
433  given_rank = isotopeAverageMono.getIsotopeRank();
434  if(given_rank < asked_rank)
435  {
436  more_rank = false;
437  }
438  else if(isotopeAverageMono.getIntensityRatio() == 0)
439  {
440  more_rank = false;
441  }
442  else
443  {
444  // isotopeAverageMono.makePeptideNaturalIsotopeAverageSp();
445  _isotope_mass_list.push_back(
446  isotopeAverageMono.makePeptideNaturalIsotopeAverageSp());
447  }
448  }
449  }
450  }
451 }
452 
453 void
455 {
456 
457  if(_ms_level > 1)
458  {
459  pappso_double precursor_mz_1 = _peptide_sp->getMz(1);
460  _custom_plot->highlightPrecursorPeaks(
461  precursor_mz_1, 1, _p_ms2_precision);
462  pappso_double precursor_mz_charge = _peptide_sp->getMz(_peptide_charge);
463  _custom_plot->highlightPrecursorPeaks(
464  precursor_mz_charge, _peptide_charge, _p_ms2_precision);
465  }
466 }
void toQPaintDevice(QPaintDevice *device, const QSize &size) override
bool savePdf(const QString &fileName, int width=0, int height=0)
std::vector< pappso::PeptideNaturalIsotopeAverageSp > _isotope_mass_list
list of isotope precursors
void setQualifiedMassSpectrum(const QualifiedMassSpectrum &spectrum)
void peakChanged(pappso::DataPointCstSPtr peak_match) const
void setMs1Precision(PrecisionPtr precision)
void setMassSpectrumCstSPtr(const MassSpectrumCstSPtr &spectrum)
void setIonList(const std::list< PeptideIon > &ion_list)
void setMs2Precision(PrecisionPtr precision)
std::list< PeakIonIsotopeMatch > _peak_ion_isotope_match_list
void setPeptideCharge(unsigned int parent_ion_charge)
void setMaximumIsotopeRank(unsigned int max_isotope_rank)
MassSpectrumCstSPtr _spectrum_sp
void mzChanged(double mz) const
void setMsLevel(unsigned int ms_level)
MassSpectrumWidget(QWidget *parent=0)
void setPeptideSp(const PeptideSp &peptide_sp)
void mzChangeEvent(pappso_double mz) const
void ionChanged(pappso::PeakIonIsotopeMatchCstSPtr ion) const
void setVisibleMassDelta(bool visible)
void peakChangeEvent(const DataPoint *p_peak_match)
void setMaximumIsotopeNumber(unsigned int max_isotope_number)
std::list< PeptideIon > _ion_list
bool contains(pappso_double) const
Definition: mzrange.cpp:120
static std::list< PeptideIon > getCIDionList()
static const QColor getPeptideIonColor(PeptideIon ion_type)
const std::list< PeakIonIsotopeMatch > & getPeakIonIsotopeMatchList() const
PeptideNaturalIsotopeAverageSp makePeptideNaturalIsotopeAverageSp() const
const std::map< unsigned int, pappso_double > getIntensityRatioPerIsotopeNumber() const
virtual QString toString() const =0
void addPeakIonIsotopeMatch(const PeakIonIsotopeMatch &peak_ion_match)
void addMassDelta(const PeakIonIsotopeMatch &peak_ion_match)
void setSpectrumP(const MassSpectrum *spectrum)
void addMs1IsotopePattern(const std::vector< pappso::PeptideNaturalIsotopeAverageSp > &isotope_mass_list, pappso_double intensity)
Class representing a fully specified mass spectrum.
uint getMsLevel() const
Get the mass spectrum level.
MassSpectrumCstSPtr getMassSpectrumCstSPtr() const
Get the MassSpectrumCstSPtr.
uint getPrecursorCharge(bool *ok=nullptr) const
Get the precursor charge.
tries to keep as much as possible monoisotopes, removing any possible C13 peaks and changes multichar...
Definition: aa.cpp:39
std::shared_ptr< const Peptide > PeptideSp
double pappso_double
A type definition for doubles.
Definition: types.h:49
std::shared_ptr< const PeptideNaturalIsotopeAverage > PeptideNaturalIsotopeAverageSp
std::shared_ptr< const MassSpectrum > MassSpectrumCstSPtr
Definition: massspectrum.h:55
std::shared_ptr< const DataPoint > DataPointCstSPtr
Definition: datapoint.h:17
pappso_double x
Definition: datapoint.h:23
pappso_double y
Definition: datapoint.h:24
DataPointCstSPtr makeDataPointCstSPtr() const
Definition: datapoint.cpp:70