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

#include <psmspecglobscan.h>

Inheritance diagram for pappso::cbor::psm::PsmSpecGlobScan:
pappso::cbor::psm::CborScanMapBase

Public Member Functions

 PsmSpecGlobScan (const PsmSpecGlob &psm_specgloc, pappso::PrecisionPtr fragment_tolerance)
 
virtual ~PsmSpecGlobScan ()
 
- Public Member Functions inherited from pappso::cbor::psm::CborScanMapBase
 CborScanMapBase (const PsmFileScanProcess &psm_file_scan_process)
 
virtual ~CborScanMapBase ()
 
std::vector< double >::iterator addPsmEvalVectorDouble (const QString &eval_name, const QString &eval_value_key, std::vector< double >::iterator begin, std::vector< double >::const_iterator end)
 add a new eval key and double values (from a vector) to each PSM
 
QCborMap getCborScanId () const
 
QCborMap getCborScanPrecursor () const
 
QCborArray getCborPsmList () const
 

Protected Member Functions

void process () override
 
- Protected Member Functions inherited from pappso::cbor::psm::CborScanMapBase
pappso::QualifiedMassSpectrumSPtr getCurrentQualifiedMassSpectrumSPtr () const
 

Private Attributes

const PsmSpecGlobmp_psmSpecGlob
 
pappso::PrecisionPtr m_fragmentTolerance
 
bool m_checkMutations = false
 

Additional Inherited Members

- Protected Attributes inherited from pappso::cbor::psm::CborScanMapBase
const PsmFileScanProcessm_PsmFileScanProcess
 

Detailed Description

Todo
write docs

Definition at line 37 of file psmspecglobscan.h.

Constructor & Destructor Documentation

◆ PsmSpecGlobScan()

pappso::cbor::psm::PsmSpecGlobScan::PsmSpecGlobScan ( const PsmSpecGlob & psm_specgloc,
pappso::PrecisionPtr fragment_tolerance )

Default constructor

Definition at line 35 of file psmspecglobscan.cpp.

37 : CborScanMapBase(psm_specgloc)
38{
39 m_fragmentTolerance = fragment_tolerance;
40 mp_psmSpecGlob = &psm_specgloc;
41}
CborScanMapBase(const PsmFileScanProcess &psm_file_scan_process)
pappso::PrecisionPtr m_fragmentTolerance

References pappso::cbor::psm::CborScanMapBase::CborScanMapBase(), m_fragmentTolerance, and mp_psmSpecGlob.

◆ ~PsmSpecGlobScan()

pappso::cbor::psm::PsmSpecGlobScan::~PsmSpecGlobScan ( )
virtual

Destructor

Definition at line 43 of file psmspecglobscan.cpp.

44{
45}

Member Function Documentation

◆ process()

void pappso::cbor::psm::PsmSpecGlobScan::process ( )
overrideprotectedvirtual

Reimplemented from pappso::cbor::psm::CborScanMapBase.

Definition at line 48 of file psmspecglobscan.cpp.

49{
50 // qDebug() << keys();
51 if(!keys().contains("id"))
52 {
53 throw pappso::PappsoException(QObject::tr("missing scan id"));
54 }
55 if(keys().contains("psm_list"))
56 {
58
59
60 mp_psmSpecGlob->filterMassSpectrum(
61 *(qualified_mass_spectrum.get()->getMassSpectrumSPtr().get()));
62
63 QCborArray new_psm_arr;
64 for(QCborValue cbor_psm : value("psm_list").toArray())
65 {
66 QCborMap cbor_psm_map = cbor_psm.toMap();
67
68 if(!cbor_psm_map.keys().contains("proforma"))
69 {
70 throw pappso::PappsoException(
71 QObject::tr("missing proforma in psm %1").arg(cbor_psm_map.keys().size()));
72 }
73 QCborMap cbor_psm_eval_specglob;
74 pappso::PeptideSp peptide_sp =
75 pappso::PeptideProFormaParser::parseString(cbor_psm_map.value("proforma").toString());
76
78 std::make_shared<pappso::specglob::PeptideSpectrum>(peptide_sp);
79
80
81 pappso::specglob::ExperimentalSpectrumCsp experimental_spectrum =
82 std::make_shared<pappso::specglob::ExperimentalSpectrum>(*qualified_mass_spectrum.get(),
84 pappso::specglob::SpectralAlignment spectral_alignment(pappso::specglob::ScoreValues(),
86
87 spectral_alignment.align(peptide_spectra, experimental_spectrum);
88
89
90 cbor_psm_eval_specglob.insert(QString("max_score"), spectral_alignment.getMaxScore());
91
92
93 if(spectral_alignment.getMaxScore() > 0)
94 {
95 pappso::specglob::PeptideModel peptide_model = spectral_alignment.buildPeptideModel();
96
97 // peptide_model.eliminateComplementaryDelta(m_fragmentTolerance);
98 // peptide_model.matchExperimentalPeaks(m_fragmentTolerance);
99
100
101 pappso::specglob::PostTreatment post_treatment(m_fragmentTolerance, peptide_model);
103 post_treatment.findReplaceMutations();
104 peptide_model = post_treatment.getBetterPeptideModel();
105
106 cbor_psm_eval_specglob.insert(QString("count"),
107 (qint64)peptide_model.getCountSharedPeaks());
108
109 cbor_psm_eval_specglob.insert(QString("proforma"), peptide_model.toProForma());
110 cbor_psm_eval_specglob.insert(QString("bracket"), peptide_model.toString());
111
112 // peptide_model.checkForMutations()
113 }
114
115
116 QCborMap psm_eval = cbor_psm_map.value("eval").toMap();
117 psm_eval.remove(QString("specglob"));
118 psm_eval.insert(QString("specglob"), cbor_psm_eval_specglob);
119 cbor_psm_map.remove(QString("eval"));
120 cbor_psm_map.insert(QString("eval"), psm_eval);
121
122 new_psm_arr.push_back(cbor_psm_map);
123 }
124
125 // insert(QString("psm_list"), new_psm_arr);
126 remove(QString("psm_list"));
127 insert(QString("psm_list"), new_psm_arr);
128 }
129}
static PeptideSp parseString(const QString &pepstr)
pappso::QualifiedMassSpectrumSPtr getCurrentQualifiedMassSpectrumSPtr() const
std::size_t getCountSharedPeaks() const
QString toProForma() const
get the peptide model in ProForma notation https://github.com/HUPO-PSI/ProForma/blob/master/README....
std::shared_ptr< const PeptideSpectrum > PeptideSpectraCsp
std::shared_ptr< const ExperimentalSpectrum > ExperimentalSpectrumCsp
std::shared_ptr< QualifiedMassSpectrum > QualifiedMassSpectrumSPtr
std::shared_ptr< const Peptide > PeptideSp

References pappso::specglob::SpectralAlignment::align(), pappso::specglob::SpectralAlignment::buildPeptideModel(), pappso::specglob::PostTreatment::findReplaceMutations(), pappso::specglob::PostTreatment::getBetterPeptideModel(), pappso::specglob::PeptideModel::getCountSharedPeaks(), pappso::cbor::psm::CborScanMapBase::getCurrentQualifiedMassSpectrumSPtr(), pappso::specglob::SpectralAlignment::getMaxScore(), m_checkMutations, m_fragmentTolerance, mp_psmSpecGlob, pappso::PeptideProFormaParser::parseString(), pappso::specglob::PeptideModel::toProForma(), and pappso::specglob::PeptideModel::toString().

Member Data Documentation

◆ m_checkMutations

bool pappso::cbor::psm::PsmSpecGlobScan::m_checkMutations = false
private

Definition at line 57 of file psmspecglobscan.h.

Referenced by process().

◆ m_fragmentTolerance

pappso::PrecisionPtr pappso::cbor::psm::PsmSpecGlobScan::m_fragmentTolerance
private

Definition at line 56 of file psmspecglobscan.h.

Referenced by PsmSpecGlobScan(), and process().

◆ mp_psmSpecGlob

const PsmSpecGlob* pappso::cbor::psm::PsmSpecGlobScan::mp_psmSpecGlob
private

Definition at line 55 of file psmspecglobscan.h.

Referenced by PsmSpecGlobScan(), and process().


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