| 115 | } |
| 116 | |
| 117 | void ProteinResolver::resolveConsensus(ConsensusMap & consensus) |
| 118 | { |
| 119 | vector<ProteinEntry> * protein_nodes = new vector<ProteinEntry>; |
| 120 | vector<PeptideEntry> * peptide_nodes = new vector<PeptideEntry>; |
| 121 | vector<ISDGroup> * isd_groups = new vector<ISDGroup>; |
| 122 | vector<MSDGroup> * msd_groups = new vector<MSDGroup>; |
| 123 | vector<Size> * reindexed_proteins = new vector<Size>; |
| 124 | vector<Size> * reindexed_peptides = new vector<Size>; |
| 125 | |
| 126 | // building ISD Groups |
| 127 | buildingISDGroups_(*protein_nodes, *peptide_nodes, *isd_groups); |
| 128 | |
| 129 | // Including all MSMS derived peptides into the graph |
| 130 | includeMSMSPeptides_(consensus, *peptide_nodes); |
| 131 | |
| 132 | // building MSD Groups |
| 133 | buildingMSDGroups_(*msd_groups, *isd_groups); |
| 134 | |
| 135 | // calculations + reindexing |
| 136 | reindexingNodes_(*msd_groups, *reindexed_proteins, *reindexed_peptides); |
| 137 | |
| 138 | // compute intensity of a msd group |
| 139 | computeIntensityOfMSD_(*msd_groups); |
| 140 | |
| 141 | primaryProteins_(*peptide_nodes, *reindexed_peptides); |
| 142 | //TODO indistinguishableProteins(msd_groups); |
| 143 | |
| 144 | // |
| 145 | countTargetDecoy(*msd_groups, consensus); |
| 146 | |
| 147 | ResolverResult result; |
| 148 | //result.identifier = file_identifier; |
| 149 | result.isds = isd_groups; |
| 150 | result.msds = msd_groups; |
| 151 | result.peptide_entries = peptide_nodes; |
| 152 | result.protein_entries = protein_nodes; |
| 153 | result.reindexed_peptides = reindexed_peptides; |
| 154 | result.reindexed_proteins = reindexed_proteins; |
| 155 | result.input_type = ResolverResult::Consensus; |
| 156 | result.consensus_map = &consensus; |
| 157 | |
| 158 | resolver_result_.push_back(result); |
| 159 | } |
| 160 | |
| 161 | void ProteinResolver::computeIntensityOfMSD_(vector<MSDGroup> & msd_groups) |
| 162 | { |
no test coverage detected