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
38using namespace pappso;
39
40MassSpectrumWidget::MassSpectrumWidget(QWidget *parent)
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);
53 qDebug();
54}
56{
57}
58
59bool
60MassSpectrumWidget::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
74void
75MassSpectrumWidget::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}
86void
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}
118void
120{
121 qDebug();
123 qDebug();
124 _custom_plot->clearItems();
126 qDebug();
127}
128
129void
130MassSpectrumWidget::setPeptideCharge(unsigned int parent_ion_charge)
131{
132 _peptide_charge = parent_ion_charge;
133}
134void
135MassSpectrumWidget::setIonList(const std::list<PeptideIon> &ion_list)
136{
137 _ion_list = ion_list;
138}
139
140void
141MassSpectrumWidget::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 {
154 }
155
156 //_precision._precision = precision._precision;
157}
158void
160{
161 qDebug() << "precision->toString()=" << precision->toString();
162 _p_ms1_precision = precision;
163 //_precision._precision = precision._precision;
164}
165void
167{
168 qDebug() << "precision->toString()=" << precision->toString();
169 _p_ms2_precision = precision;
170 //_precision._precision = precision._precision;
171}
172
173void
174MassSpectrumWidget::setMaximumIsotopeNumber(unsigned int max_isotope_number)
175{
176 _max_isotope_number = max_isotope_number;
177}
178
179void
180MassSpectrumWidget::setMaximumIsotopeRank(unsigned int max_isotope_rank)
181{
182 _max_isotope_rank = max_isotope_rank;
183}
184void
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()),
201 _ion_list,
204
206 }
207 else
208 {
209 }
211 }
212 qDebug();
213}
214
215void
217{
218 qDebug();
219 _peptide_sp = peptide_sp;
220
221 // clearData();
222 qDebug();
223}
224
225void
227{
228 qDebug();
229 _spectrum_sp = spectrum;
230
231 clearData();
232 qDebug();
233}
234
235void
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
251void
253 const QualifiedMassSpectrum &spectrum)
254{
255 qDebug() << "spectrum.getPrecursorCharge()=" << spectrum.getPrecursorCharge();
256
257 setMsLevel(spectrum.getMsLevel());
259
260 qDebug();
261}
262
263void
265{
266
267 qDebug();
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}
364void
366{
367 emit mzChanged(mz);
368}
369
370void
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;
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
410void
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 =
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
453void
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)
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:18
pappso_double x
Definition: datapoint.h:23
pappso_double y
Definition: datapoint.h:24
DataPointCstSPtr makeDataPointCstSPtr() const
Definition: datapoint.cpp:70