Skip to content
Snippets Groups Projects
Commit 7f7e2e71 authored by Langella Olivier's avatar Langella Olivier
Browse files

new column added to display experimental precursor m/z

parent ecebd14d
No related branches found
No related tags found
No related merge requests found
...@@ -48,9 +48,17 @@ public: ...@@ -48,9 +48,17 @@ public:
PtmSampleScan(const PtmSampleScan& other); PtmSampleScan(const PtmSampleScan& other);
~PtmSampleScan(); ~PtmSampleScan();
/** @brief get the best peptide Match (Evalue)
* peptide match are sorted in this sample/scan by Evalue
* The representative peptide match is the first in the sorted array
*/
const PeptideMatch & getRepresentativePeptideMatch() const; const PeptideMatch & getRepresentativePeptideMatch() const;
bool add(const PeptideMatch & peptide_match); bool add(const PeptideMatch & peptide_match);
std::vector<unsigned int> getBestPtmPositionList(const PtmGroupingExperiment * p_ptm_grouping_experiment) const; std::vector<unsigned int> getBestPtmPositionList(const PtmGroupingExperiment * p_ptm_grouping_experiment) const;
/** @brief get all observed PTM positions for this sample/scan
* positions are intended in a peptide (not protein), from 0 (first amino acid) to size-1 (last amino acid)
*/
std::vector<unsigned int> getObservedPtmPositionList(const PtmGroupingExperiment * p_ptm_grouping_experiment) const; std::vector<unsigned int> getObservedPtmPositionList(const PtmGroupingExperiment * p_ptm_grouping_experiment) const;
const QString getHtmlSequence(const PtmGroupingExperiment * p_ptm_grouping_experiment) const; const QString getHtmlSequence(const PtmGroupingExperiment * p_ptm_grouping_experiment) const;
std::vector<PeptideMatch> & getPeptideMatchList(); std::vector<PeptideMatch> & getPeptideMatchList();
......
...@@ -68,7 +68,7 @@ int PeptideTableModel::rowCount(const QModelIndex &parent ) const { ...@@ -68,7 +68,7 @@ int PeptideTableModel::rowCount(const QModelIndex &parent ) const {
return 0; return 0;
} }
int PeptideTableModel::columnCount(const QModelIndex &parent ) const { int PeptideTableModel::columnCount(const QModelIndex &parent ) const {
return 19; return 20;
} }
QVariant PeptideTableModel::headerData(int section, Qt::Orientation orientation, int role) const QVariant PeptideTableModel::headerData(int section, Qt::Orientation orientation, int role) const
{ {
...@@ -148,6 +148,9 @@ QVariant PeptideTableModel::data(const QModelIndex &index, int role ) const { ...@@ -148,6 +148,9 @@ QVariant PeptideTableModel::data(const QModelIndex &index, int role ) const {
case (std::int8_t)PeptideListColumn::charge: case (std::int8_t)PeptideListColumn::charge:
return QVariant ((quint32) _p_protein_match->getPeptideMatchList().at(row).getPeptideEvidence()->getCharge()); return QVariant ((quint32) _p_protein_match->getPeptideMatchList().at(row).getPeptideEvidence()->getCharge());
break; break;
case (std::int8_t)PeptideListColumn::experimental_mz:
return QVariant ((qreal) _p_protein_match->getPeptideMatchList().at(row).getPeptideEvidence()->getExperimentalMz());
break;
case (std::int8_t)PeptideListColumn::sequence: case (std::int8_t)PeptideListColumn::sequence:
return _p_protein_match->getPeptideMatchList().at(row).getPeptideEvidence()->getPeptideXtpSp().get()->getSequence(); return _p_protein_match->getPeptideMatchList().at(row).getPeptideEvidence()->getPeptideXtpSp().get()->getSequence();
break; break;
...@@ -251,6 +254,9 @@ const QString PeptideTableModel::getTitle(std::int8_t column) { ...@@ -251,6 +254,9 @@ const QString PeptideTableModel::getTitle(std::int8_t column) {
case (std::int8_t) PeptideListColumn::charge: case (std::int8_t) PeptideListColumn::charge:
return "charge"; return "charge";
break; break;
case (std::int8_t) PeptideListColumn::experimental_mz:
return "observed m/z";
break;
case (std::int8_t) PeptideListColumn::sequence: case (std::int8_t) PeptideListColumn::sequence:
return "sequence"; return "sequence";
break; break;
...@@ -317,6 +323,9 @@ const QString PeptideTableModel::getDescription(std::int8_t column) { ...@@ -317,6 +323,9 @@ const QString PeptideTableModel::getDescription(std::int8_t column) {
case (std::int8_t) PeptideListColumn::charge: case (std::int8_t) PeptideListColumn::charge:
return "peptide charge"; return "peptide charge";
break; break;
case (std::int8_t) PeptideListColumn::experimental_mz:
return "observed m/z (experimental m/z of precursor)";
break;
case (std::int8_t) PeptideListColumn::sequence: case (std::int8_t) PeptideListColumn::sequence:
return "peptide sequence"; return "peptide sequence";
break; break;
......
...@@ -39,18 +39,19 @@ enum class PeptideListColumn: std::int8_t { ...@@ -39,18 +39,19 @@ enum class PeptideListColumn: std::int8_t {
rtmin=4, rtmin=4,
rt=5, rt=5,
charge=6, charge=6,
sequence=7, experimental_mz=7,
modifs=8, sequence=8,
label=9, modifs=9,
start=10, label=10,
length=11, start=11,
used=12, length=12,
subgroups=13, used=13,
Evalue=14, subgroups=14,
theoretical_mhplus=15, Evalue=15,
delta_mhplus=16, theoretical_mhplus=16,
hyperscore=17, delta_mhplus=17,
mascot_score=18 hyperscore=18,
mascot_score=19
}; };
class PeptideListWindow; class PeptideListWindow;
......
...@@ -56,7 +56,8 @@ void PtmSpectraSheet::writeCellHeader(PeptideListColumn column) { ...@@ -56,7 +56,8 @@ void PtmSpectraSheet::writeCellHeader(PeptideListColumn column) {
void PtmSpectraSheet::writeHeaders(IdentificationGroup * p_ident) { void PtmSpectraSheet::writeHeaders(IdentificationGroup * p_ident) {
// Group ID Peptide ID Sample Scan Rt (minutes) Sequence (top) Modifs (top) Number of subgroups Sub-groups Ids Best E-value Best hyperscore m/z Obs Charge MH+ Obs MH+ theo DeltaMH+ Delta-ppm // Group ID Sub-group ID Phosphopeptide ID Sample Scan Rt Sequence (top) Modifs (top) Best position in peptide All observed positions in phosphopeptide Number of phosphoislands Phosphoislands Ids Best e-value Charge MH+ Obs MH+ theo DeltaMH+ Delta-ppm
//MS Sample : 20120906_balliau_extract_1_A01_urnb-1 //MS Sample : 20120906_balliau_extract_1_A01_urnb-1
const std::vector<MsRunSp> msrun_list = p_ident->getMsRunSpList(); const std::vector<MsRunSp> msrun_list = p_ident->getMsRunSpList();
...@@ -73,7 +74,7 @@ void PtmSpectraSheet::writeHeaders(IdentificationGroup * p_ident) { ...@@ -73,7 +74,7 @@ void PtmSpectraSheet::writeHeaders(IdentificationGroup * p_ident) {
_p_writer->writeCell("Peptide ID"); _p_writer->writeCell("Peptide ID");
_p_writer->setCellAnnotation("MS sample name (MS run)"); _p_writer->setCellAnnotation("MS sample name (MS run)");
_p_writer->writeCell("Sample"); _p_writer->writeCell("Sample");
_p_writer->writeCell("Scan"); writeCellHeader(PeptideListColumn::scan);
_p_writer->setCellAnnotation("retention time in seconds"); _p_writer->setCellAnnotation("retention time in seconds");
_p_writer->writeCell("Rt (seconds)"); _p_writer->writeCell("Rt (seconds)");
_p_writer->setCellAnnotation("experimental precursor m/z"); _p_writer->setCellAnnotation("experimental precursor m/z");
......
...@@ -71,13 +71,11 @@ void SpectraSheet::writeHeaders(IdentificationGroup * p_ident) { ...@@ -71,13 +71,11 @@ void SpectraSheet::writeHeaders(IdentificationGroup * p_ident) {
_p_writer->writeLine(); _p_writer->writeLine();
_p_writer->writeCell("Group ID"); _p_writer->writeCell("Group ID");
_p_writer->writeCell("Peptide ID"); _p_writer->writeCell("Peptide ID");
_p_writer->setCellAnnotation("MS sample name (MS run)");
_p_writer->writeCell("Sample"); writeCellHeader(PeptideListColumn::sample);
_p_writer->writeCell("Scan"); writeCellHeader(PeptideListColumn::scan);
_p_writer->setCellAnnotation("retention time in seconds"); writeCellHeader(PeptideListColumn::rt);
_p_writer->writeCell("Rt (seconds)"); writeCellHeader(PeptideListColumn::experimental_mz);
_p_writer->setCellAnnotation("experimental precursor m/z");
_p_writer->writeCell("Precursor m/z");
_p_writer->setCellAnnotation("peptide sequence (best match)"); _p_writer->setCellAnnotation("peptide sequence (best match)");
_p_writer->writeCell("Sequence (top)"); _p_writer->writeCell("Sequence (top)");
_p_writer->setCellAnnotation("peptide charge (best match)"); _p_writer->setCellAnnotation("peptide charge (best match)");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment