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

WIP: using msrunretention object

parent cd0259e3
No related branches found
No related tags found
No related merge requests found
......@@ -309,18 +309,46 @@ MsRun::getMsRunXicExtractorInterfaceSp()
}
void
MsRun::initializeMsRunRetentionTime(
MsRun::buildMsRunRetentionTime(
const PeptideEvidenceStore &peptide_evidence_store)
{
pappso::MsRunReaderSPtr reader = getMsRunReaderSPtr();
mpa_msrunRetentionTime =
new pappso::MsRunRetentionTime<const pappso::Peptide *>(reader);
if(mpa_msrunRetentionTime == nullptr)
{
mpa_msrunRetentionTime =
new pappso::MsRunRetentionTime<const pappso::Peptide *>(reader);
}
for(auto &peptide_evidence : peptide_evidence_store.getPeptideEvidenceList())
{
mpa_msrunRetentionTime->addPeptideAsSeamark(
peptide_evidence.get()->getPeptideXtpSp().get()->getNativePeptideP(),
reader.get()->scanNumber2SpectrumIndex(
peptide_evidence.get()->getScanNumber()));
if(peptide_evidence.get()->getMsRunP() == this)
{
mpa_msrunRetentionTime->addPeptideAsSeamark(
peptide_evidence.get()
->getPeptideXtpSp()
.get()
->getNativePeptideP(),
reader.get()->scanNumber2SpectrumIndex(
peptide_evidence.get()->getScanNumber()));
}
}
// mpa_msrunRetentionTime->computePeptideRetentionTimes();
}
void
MsRun::computeMsRunRetentionTime()
{
if(mpa_msrunRetentionTime == nullptr)
{
mpa_msrunRetentionTime->computePeptideRetentionTimes();
}
}
void
MsRun::clearMsRunRetentionTime()
{
if(mpa_msrunRetentionTime == nullptr)
{
delete mpa_msrunRetentionTime;
mpa_msrunRetentionTime = nullptr;
}
}
......@@ -81,8 +81,10 @@ class MsRun : public pappso::MsRunId
pappso::MsRunXicExtractorInterfaceSp getMsRunXicExtractorInterfaceSp();
void initializeMsRunRetentionTime(
const PeptideEvidenceStore &peptide_evidence_store);
void
buildMsRunRetentionTime(const PeptideEvidenceStore &peptide_evidence_store);
void computeMsRunRetentionTime();
void clearMsRunRetentionTime();
private:
std::map<MsRunStatistics, QVariant> _param_stats;
......
......@@ -277,22 +277,22 @@ ExportMasschroqDialog::doFindBestMsrunForAlignment()
// p_identification_group =
// mp_project->getIdentificationGroupList()[0]->get
}
std::vector<MsRunSp> msrun_list;
bool is_ok = true;
for(IdentificationDataSourceSp p_ident_data_source :
mp_project->getIdentificationDataSourceStore()
.getIdentificationDataSourceList())
{
MsRunSp msrun_sp = p_ident_data_source->getMsRunSp();
msrun_list.push_back(msrun_sp);
if(msrun_sp->findMsRunFile())
{
// msrun_sp->checkMsRunStatistics();
msrun_sp->initializeMsRunRetentionTime(
msrun_sp->buildMsRunRetentionTime(
p_ident_data_source->getPeptideEvidenceStore());
}
else
{
is_ok = false;
ui->warningDirectoryLabel->setText(
tr("%1 not found : Please choose the directory to look for this MS "
"run file.")
......@@ -301,4 +301,19 @@ ExportMasschroqDialog::doFindBestMsrunForAlignment()
break;
}
}
if(is_ok)
{
for(auto &msrun_sp : mp_project->getMsRunStore().getMsRunList())
{
msrun_sp->computeMsRunRetentionTime();
}
//find the best reference
}
else
{
for(auto &msrun_sp : mp_project->getMsRunStore().getMsRunList())
{
msrun_sp->clearMsRunRetentionTime();
}
}
}
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