libpappsomspp
Library for mass spectrometry
Loading...
Searching...
No Matches
pappso::PeptideRawFragmentMasses Class Reference

#include <peptiderawfragmentmasses.h>

Public Member Functions

 PeptideRawFragmentMasses (const Peptide &peptide, RawFragmentationMode mode)
 ~PeptideRawFragmentMasses ()
void pushBackIonMasses (std::vector< pappso_double > &mass_list, Enums::PeptideIon ion_type) const
void pushBackIonMz (std::vector< pappso_double > &mass_list, Enums::PeptideIon ion_type, unsigned int charge) const
void pushBackMatchSpectrum (std::vector< SimplePeakIonMatch > &peak_match_list, const MassSpectrum &spectrum, PrecisionPtr precision, Enums::PeptideIon ion_type, unsigned int charge) const

Static Public Member Functions

static pappso_double getDeltaMass (Enums::PeptideIon ion_type)

Private Types

using ionDeltatMzMassMap = pappso_double[20]

Private Attributes

std::vector< pappso_doublem_cumulativeNterMasses
 cumulative Nter masses (without internal Nter modification)
std::vector< pappso_doublem_cumulativeCterMasses
 cumulative Cter masses (without internal Cter modification)

Static Private Attributes

static ionDeltatMzMassMap m_ionDeltaMz

Detailed Description

Definition at line 59 of file peptiderawfragmentmasses.h.

Member Typedef Documentation

◆ ionDeltatMzMassMap

Constructor & Destructor Documentation

◆ PeptideRawFragmentMasses()

pappso::PeptideRawFragmentMasses::PeptideRawFragmentMasses ( const Peptide & peptide,
RawFragmentationMode mode )

Definition at line 61 of file peptiderawfragmentmasses.cpp.

63{
64
65 qDebug();
66 std::vector<Aa>::const_iterator it(peptide.begin());
67 std::vector<Aa>::const_iterator end(peptide.end());
68 if(it == end)
69 return;
70
71 qDebug();
72 pappso_double cumulative_mass = it->getMass();
73 if(peptide.getNterModification() != nullptr)
74 {
75 cumulative_mass += peptide.getNterModification()->getMass();
76 }
77 qDebug();
78 m_cumulativeNterMasses.push_back(cumulative_mass);
79 it++;
80
81 if(it != end)
82 {
83 end--;
84 qDebug();
85 while(it != end)
86 {
87
88 qDebug() << " cumulative_mass=" << cumulative_mass;
89 cumulative_mass += it->getMass();
90 if((mode == RawFragmentationMode::proline_effect) && (it->getLetter() == 'P'))
91 {
92 }
93 else
94 {
95 m_cumulativeNterMasses.push_back(cumulative_mass);
96 }
97 it++;
98 }
99 }
100
101 qDebug();
102 std::vector<Aa>::const_reverse_iterator rit(peptide.rbegin());
103 std::vector<Aa>::const_reverse_iterator ritf(peptide.rbegin());
104 std::vector<Aa>::const_reverse_iterator rend(peptide.rend());
105 ritf++;
106 cumulative_mass = rit->getMass();
107 if(peptide.getCterModification() != nullptr)
108 {
109 cumulative_mass += peptide.getCterModification()->getMass();
110 }
111 if((mode == RawFragmentationMode::proline_effect) && (ritf != rend) && (ritf->getLetter() == 'P'))
112 {
113 }
114 else
115 {
116 m_cumulativeCterMasses.push_back(cumulative_mass);
117 }
118 qDebug();
119
120 rit++;
121 ritf++;
122 if(rit != rend)
123 {
124 rend--;
125 while(rit != rend)
126 {
127 cumulative_mass += rit->getMass();
128 if((mode == RawFragmentationMode::proline_effect) && (ritf != peptide.rend()) &&
129 (ritf->getLetter() == 'P'))
130 {
131 }
132 else
133 {
134 m_cumulativeCterMasses.push_back(cumulative_mass);
135 }
136 rit++;
137 ritf++;
138 }
139 }
140
141 qDebug();
142}
std::vector< pappso_double > m_cumulativeCterMasses
cumulative Cter masses (without internal Cter modification)
std::vector< pappso_double > m_cumulativeNterMasses
cumulative Nter masses (without internal Nter modification)
double pappso_double
A type definition for doubles.
Definition types.h:60

References pappso::Peptide::begin(), pappso::Peptide::end(), pappso::Peptide::getCterModification(), pappso::AaModification::getMass(), pappso::Peptide::getNterModification(), m_cumulativeCterMasses, m_cumulativeNterMasses, pappso::proline_effect, pappso::Peptide::rbegin(), and pappso::Peptide::rend().

◆ ~PeptideRawFragmentMasses()

pappso::PeptideRawFragmentMasses::~PeptideRawFragmentMasses ( )

Definition at line 222 of file peptiderawfragmentmasses.cpp.

223{
224}

Member Function Documentation

◆ getDeltaMass()

pappso_double pappso::PeptideRawFragmentMasses::getDeltaMass ( Enums::PeptideIon ion_type)
static

◆ pushBackIonMasses()

void pappso::PeptideRawFragmentMasses::pushBackIonMasses ( std::vector< pappso_double > & mass_list,
Enums::PeptideIon ion_type ) const

Definition at line 145 of file peptiderawfragmentmasses.cpp.

147{
148 const std::vector<pappso_double> *p_mass_list = &m_cumulativeCterMasses;
149
150 if(peptideIonIsNter(ion_type))
151 {
152 // nter
153 p_mass_list = &m_cumulativeNterMasses;
154 }
155
156 for(pappso_double mass : *p_mass_list)
157 {
158 mass_list.push_back(mass + m_ionDeltaMz[(std::int8_t)ion_type]);
159 }
160}
bool peptideIonIsNter(Enums::PeptideIon ion_type)
tells if an ion is Nter
Definition peptide.cpp:87

References m_cumulativeCterMasses, m_cumulativeNterMasses, m_ionDeltaMz, and pappso::peptideIonIsNter().

◆ pushBackIonMz()

void pappso::PeptideRawFragmentMasses::pushBackIonMz ( std::vector< pappso_double > & mass_list,
Enums::PeptideIon ion_type,
unsigned int charge ) const

Definition at line 163 of file peptiderawfragmentmasses.cpp.

166{
167 const std::vector<pappso_double> *p_mass_list = &m_cumulativeCterMasses;
168
169 if(peptideIonIsNter(ion_type))
170 {
171 // nter
172 p_mass_list = &m_cumulativeNterMasses;
173 }
174
175 for(pappso_double mass : *p_mass_list)
176 {
177 mass_list.push_back((mass + m_ionDeltaMz[(std::int8_t)ion_type] + (MHPLUS * charge)) /
178 charge);
179 }
180}
const pappso_double MHPLUS(1.007276466879)

References m_cumulativeCterMasses, m_cumulativeNterMasses, m_ionDeltaMz, pappso::MHPLUS(), and pappso::peptideIonIsNter().

Referenced by pappso::MorpheusScore::MorpheusScore(), and pushBackMatchSpectrum().

◆ pushBackMatchSpectrum()

void pappso::PeptideRawFragmentMasses::pushBackMatchSpectrum ( std::vector< SimplePeakIonMatch > & peak_match_list,
const MassSpectrum & spectrum,
PrecisionPtr precision,
Enums::PeptideIon ion_type,
unsigned int charge ) const

Definition at line 183 of file peptiderawfragmentmasses.cpp.

188{
189 std::vector<pappso_double> mass_list;
190 pushBackIonMz(mass_list, ion_type, charge);
191
192 // no need to sort
193 // std::sort(mass_list.begin(), mass_list.end());
194
195 std::vector<pappso_double>::iterator it_mz = mass_list.begin();
196 std::vector<pappso_double>::iterator it_mz_end = mass_list.end();
197
198 // scan products over each peak in spectrum :
199 std::vector<DataPoint>::const_iterator it_peak = spectrum.begin();
200 std::vector<DataPoint>::const_iterator it_peak_end = spectrum.end();
201 unsigned int ion_size = 1;
202 while((it_peak != it_peak_end) && (it_mz != it_mz_end))
203 {
204 MzRange massrange(it_peak->x, precision);
205 if((*it_mz) > massrange.upper())
206 {
207 it_peak++;
208 continue;
209 }
210 if((*it_mz) < massrange.lower())
211 {
212 it_mz++;
213 ion_size++;
214 continue;
215 }
216 peak_match_list.push_back({(*it_peak), ion_type, ion_size, charge, (*it_mz)});
217 it_mz++;
218 ion_size++;
219 }
220}
void pushBackIonMz(std::vector< pappso_double > &mass_list, Enums::PeptideIon ion_type, unsigned int charge) const

References pappso::MzRange::lower(), pushBackIonMz(), and pappso::MzRange::upper().

Referenced by pappso::XtandemHyperscoreBis::computeXtandemHyperscore().

Member Data Documentation

◆ m_cumulativeCterMasses

std::vector<pappso_double> pappso::PeptideRawFragmentMasses::m_cumulativeCterMasses
private

cumulative Cter masses (without internal Cter modification)

Definition at line 72 of file peptiderawfragmentmasses.h.

Referenced by PeptideRawFragmentMasses(), pushBackIonMasses(), and pushBackIonMz().

◆ m_cumulativeNterMasses

std::vector<pappso_double> pappso::PeptideRawFragmentMasses::m_cumulativeNterMasses
private

cumulative Nter masses (without internal Nter modification)

Definition at line 68 of file peptiderawfragmentmasses.h.

Referenced by PeptideRawFragmentMasses(), pushBackIonMasses(), and pushBackIonMz().

◆ m_ionDeltaMz

PeptideRawFragmentMasses::ionDeltatMzMassMap pappso::PeptideRawFragmentMasses::m_ionDeltaMz
staticprivate
Initial value:
= {
[(std::int8_t)Enums::PeptideIon::b] = 0,
[(std::int8_t)Enums::PeptideIon::bstar] = -MASSNH3,
[(std::int8_t)Enums::PeptideIon::bo] = -MASSH2O,
[(std::int8_t)Enums::PeptideIon::a] = -MASSCO,
[(std::int8_t)Enums::PeptideIon::ao] = -MASSCO - MASSH2O,
[(std::int8_t)Enums::PeptideIon::bp] = -1,
[(std::int8_t)Enums::PeptideIon::c] = MASSNH3,
[(std::int8_t)Enums::PeptideIon::y] = MASSH2O,
[(std::int8_t)Enums::PeptideIon::yo] = 0,
[(std::int8_t)Enums::PeptideIon::yp] = -1,
@ a
Nter aldimine ions.
Definition types.h:290
@ y
Cter amino ions.
Definition types.h:295
@ c
Nter amino ions.
Definition types.h:294
@ astar
Nter aldimine ions + NH3 loss.
Definition types.h:291
@ ystar
Cter amino ions + NH3 loss.
Definition types.h:296
@ yo
Cter amino ions + H2O loss.
Definition types.h:297
@ bstar
Nter acylium ions + NH3 loss.
Definition types.h:288
@ b
Nter acylium ions.
Definition types.h:287
@ x
Cter acylium ions.
Definition types.h:300
@ bo
Nter acylium ions + H2O loss.
Definition types.h:289
@ ao
Nter aldimine ions + H2O loss.
Definition types.h:292
@ z
Cter carbocations.
Definition types.h:298
const pappso_double MASSCO(MASSCARBON+MASSOXYGEN)
const pappso_double MPROTIUM(1.007825032241)
const pappso_double MASSH2O((MPROTIUM *2)+MASSOXYGEN)
const pappso_double MASSNH3((MPROTIUM *3)+MASSNITROGEN)
const pappso_double MASSNITROGEN(14.0030740048)
const pappso_double MASSOXYGEN(15.99491461956)

Definition at line 64 of file peptiderawfragmentmasses.h.

Referenced by getDeltaMass(), pushBackIonMasses(), and pushBackIonMz().


The documentation for this class was generated from the following files: