libpappsomspp
Library for mass spectrometry
Loading...
Searching...
No Matches
pappso::cbor::psm::PsmCborUtils Class Reference

utilities to help managing data in PSM CBOR data file More...

#include <psmcborutils.h>

Classes

struct  PsmProteinRef

Static Public Member Functions

static void prepareCborScanWithSpectrum (QCborMap &cbor_scan, const pappso::QualifiedMassSpectrum &ms2_qualified_mass_spectrum)
static void prepareCborScanWithSpectrumAndPeakList (QCborMap &cbor_scan, const pappso::QualifiedMassSpectrum &ms2_qualified_mass_spectrum)
static std::vector< PsmProteinRefgetPsmProteinRefList (const QCborMap &cbor_psm)
static void setPsmProteinRefList (QCborMap &cbor_psm, const std::vector< PsmProteinRef > &protein_ref_list)
static void mergePsmProteinRefList (QCborMap &cbor_psm_destination, const QCborMap &cbor_psm_source)

Detailed Description

utilities to help managing data in PSM CBOR data file

Definition at line 44 of file psmcborutils.h.

Member Function Documentation

◆ getPsmProteinRefList()

std::vector< PsmCborUtils::PsmProteinRef > pappso::cbor::psm::PsmCborUtils::getPsmProteinRefList ( const QCborMap & cbor_psm)
static

Definition at line 92 of file psmcborutils.cpp.

93{
94 std::vector<PsmCborUtils::PsmProteinRef> protein_ref_list;
95
96 if(cbor_psm.contains(QString("protein_list")))
97 {
98 for(auto it : cbor_psm.value("protein_list").toArray())
99 {
100 QCborMap cbor_protein_ref = it.toMap();
101 PsmCborUtils::PsmProteinRef protein_ref;
102 protein_ref.accession = cbor_protein_ref.value("accession").toString();
103
104 for(auto ref_position : cbor_protein_ref.value("positions").toArray())
105 {
106 protein_ref.positions.push_back(ref_position.toInteger());
107 }
108 protein_ref_list.push_back(protein_ref);
109 }
110 }
111
112 return protein_ref_list;
113}

References pappso::cbor::psm::PsmCborUtils::PsmProteinRef::accession, and pappso::cbor::psm::PsmCborUtils::PsmProteinRef::positions.

Referenced by mergePsmProteinRefList().

◆ mergePsmProteinRefList()

void pappso::cbor::psm::PsmCborUtils::mergePsmProteinRefList ( QCborMap & cbor_psm_destination,
const QCborMap & cbor_psm_source )
static

Definition at line 141 of file psmcborutils.cpp.

143{
144 std::vector<PsmCborUtils::PsmProteinRef> protein_ref_list =
145 getPsmProteinRefList(cbor_psm_destination);
146 std::vector<PsmCborUtils::PsmProteinRef> protein_ref_list_source =
147 getPsmProteinRefList(cbor_psm_source);
148
149 protein_ref_list.insert(
150 protein_ref_list.end(), protein_ref_list_source.begin(), protein_ref_list_source.end());
151
152
153 std::sort(protein_ref_list.begin(),
154 protein_ref_list.end(),
155 [](PsmCborUtils::PsmProteinRef &a, PsmCborUtils::PsmProteinRef &b) {
156 return a.accession > b.accession;
157 });
158
159 std::vector<PsmCborUtils::PsmProteinRef> unique_protein_ref_list;
160
161 for(auto it = protein_ref_list.begin(); it != protein_ref_list.end(); it++)
162 {
163 // qDebug() << it->proforma;
164 if(unique_protein_ref_list.size() > 0)
165 {
166 if(unique_protein_ref_list.back().accession == it->accession)
167 {
168 // merge positions
169 unique_protein_ref_list.back().positions.insert(
170 unique_protein_ref_list.back().positions.end(),
171 it->positions.begin(),
172 it->positions.end());
173
174 std::sort(unique_protein_ref_list.back().positions.begin(),
175 unique_protein_ref_list.back().positions.end());
176
177 auto last = std::unique(unique_protein_ref_list.back().positions.begin(),
178 unique_protein_ref_list.back().positions.end());
179 // v now holds {1 2 3 4 5 x x}, where 'x' is indeterminate
180 unique_protein_ref_list.back().positions.erase(
181 last, unique_protein_ref_list.back().positions.end());
182 }
183 else
184 {
185 unique_protein_ref_list.push_back(*it);
186 }
187 }
188 else
189 {
190 unique_protein_ref_list.push_back(*it);
191 }
192
193 qDebug();
194 }
195
196 setPsmProteinRefList(cbor_psm_destination, unique_protein_ref_list);
197}
static void setPsmProteinRefList(QCborMap &cbor_psm, const std::vector< PsmProteinRef > &protein_ref_list)
static std::vector< PsmProteinRef > getPsmProteinRefList(const QCborMap &cbor_psm)

References pappso::a, pappso::b, getPsmProteinRefList(), pappso::last, and setPsmProteinRefList().

Referenced by pappso::cbor::psm::CborScanMapBase::filterPsmListUniqueUniqueProforma().

◆ prepareCborScanWithSpectrum()

void pappso::cbor::psm::PsmCborUtils::prepareCborScanWithSpectrum ( QCborMap & cbor_scan,
const pappso::QualifiedMassSpectrum & ms2_qualified_mass_spectrum )
static

Definition at line 38 of file psmcborutils.cpp.

40{
41 // id
42 QCborMap cbor_scan_id;
43 cbor_scan_id.insert(QString("index"),
44 (qint64)ms2_qualified_mass_spectrum.getMassSpectrumId().getSpectrumIndex());
45 cbor_scan_id.insert(QString("native_id"),
46 ms2_qualified_mass_spectrum.getMassSpectrumId().getNativeId());
47 cbor_scan.insert(QString("id"), cbor_scan_id.toCborValue());
48
49
50 // precursor
51 QCborMap cbor_scan_precursor;
52 cbor_scan_precursor.insert(QString("z"), ms2_qualified_mass_spectrum.getPrecursorCharge());
53 cbor_scan_precursor.insert(QString("mz"), ms2_qualified_mass_spectrum.getPrecursorMz());
54 cbor_scan_precursor.insert(QString("mh"), ms2_qualified_mass_spectrum.getPrecursorMass());
55 cbor_scan_precursor.insert(QString("intensity"),
56 ms2_qualified_mass_spectrum.getPrecursorIntensity());
57 cbor_scan.insert(QString("precursor"), cbor_scan_precursor.toCborValue());
58
59
60 // ms2
61 QCborMap cbor_scan_ms2;
62 cbor_scan_ms2.insert(QString("rt"), ms2_qualified_mass_spectrum.getRtInSeconds());
63 cbor_scan.insert(QString("ms2"), cbor_scan_ms2.toCborValue());
64}
std::size_t getSpectrumIndex() const
const QString & getNativeId() const
uint getPrecursorCharge(bool *ok=nullptr) const
get precursor charge
pappso_double getPrecursorIntensity(bool *ok=nullptr) const
get precursor intensity
double getPrecursorMass(bool *ok_p=nullptr) const
get precursor mass given the charge stats and precursor mz
const MassSpectrumId & getMassSpectrumId() const
Get the MassSpectrumId.
pappso_double getPrecursorMz(bool *ok=nullptr) const
get precursor mz
pappso_double getRtInSeconds() const
Get the retention time in seconds.

References pappso::QualifiedMassSpectrum::getMassSpectrumId(), pappso::MassSpectrumId::getNativeId(), pappso::QualifiedMassSpectrum::getPrecursorCharge(), pappso::QualifiedMassSpectrum::getPrecursorIntensity(), pappso::QualifiedMassSpectrum::getPrecursorMass(), pappso::QualifiedMassSpectrum::getPrecursorMz(), pappso::QualifiedMassSpectrum::getRtInSeconds(), and pappso::MassSpectrumId::getSpectrumIndex().

Referenced by prepareCborScanWithSpectrumAndPeakList().

◆ prepareCborScanWithSpectrumAndPeakList()

void pappso::cbor::psm::PsmCborUtils::prepareCborScanWithSpectrumAndPeakList ( QCborMap & cbor_scan,
const pappso::QualifiedMassSpectrum & ms2_qualified_mass_spectrum )
static

Definition at line 67 of file psmcborutils.cpp.

69{
70 prepareCborScanWithSpectrum(cbor_scan, ms2_qualified_mass_spectrum);
71 QCborMap spectrum_cbor;
72
73 QCborArray mz_cbor;
74 QCborArray intensity_cbor;
75 for(const pappso::DataPoint &data_point :
76 *(ms2_qualified_mass_spectrum.getMassSpectrumCstSPtr().get()))
77 {
78 mz_cbor.append(data_point.x);
79 intensity_cbor.append(data_point.y);
80 }
81 spectrum_cbor.insert(QString("mz"), mz_cbor);
82 spectrum_cbor.insert(QString("intensity"), intensity_cbor);
83
84 QCborMap new_ms2_map = cbor_scan.value("ms2").toMap();
85 new_ms2_map.insert(QString("spectrum"), spectrum_cbor.toCborValue());
86
87 cbor_scan.insert(QString("ms2"), new_ms2_map);
88}
MassSpectrumCstSPtr getMassSpectrumCstSPtr() const
Get the MassSpectrumCstSPtr.
static void prepareCborScanWithSpectrum(QCborMap &cbor_scan, const pappso::QualifiedMassSpectrum &ms2_qualified_mass_spectrum)

References pappso::QualifiedMassSpectrum::getMassSpectrumCstSPtr(), and prepareCborScanWithSpectrum().

◆ setPsmProteinRefList()

void pappso::cbor::psm::PsmCborUtils::setPsmProteinRefList ( QCborMap & cbor_psm,
const std::vector< PsmProteinRef > & protein_ref_list )
static

Definition at line 117 of file psmcborutils.cpp.

119{
120 QCborArray protein_list;
121
122 for(auto it : protein_ref_list)
123 {
124 QCborMap protein_ref;
125 protein_ref.insert(QString("accession"), it.accession);
126 QCborArray positions_arr;
127 for(auto position : it.positions)
128 {
129 positions_arr.append((qint64)position);
130 }
131 protein_ref.insert(QString("positions"), positions_arr);
132 protein_list.append(protein_ref);
133 }
134
135 cbor_psm.remove(QString("protein_list"));
136 cbor_psm.insert(QString("protein_list"), protein_list);
137}

Referenced by mergePsmProteinRefList().


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