Ejemplos de ConsensusFeature::getQuality en C++ (Cpp)

Ejemplo n.º 1

Archivo: StablePairFinder_test.cpp Proyecto: BioinformaticsArchive/OpenMS

	TEST_EQUAL(result[0].getQuality(), 0.0);

	// intermediate cases:
	// basis: feat1 < feat2 < feat3
	input[1].clear();
	input[0].push_back(ConsensusFeature(0, feat1));
	input[1].push_back(ConsensusFeature(1, feat2));
	spf.run(input, result);
	ConsensusFeature cons1 = result[0];
	TEST_EQUAL(cons1.size(), 2);
	input[0] = result;
	input[1][0] = ConsensusFeature(2, feat3);
	spf.run(input, result);
	ConsensusFeature cons2 = result[0];
	TEST_EQUAL(cons2.size(), 3);
	TEST_EQUAL(cons1.getQuality() > 0.0, true);
	TEST_EQUAL(cons2.getQuality() > 0.0, true);
	TEST_EQUAL(cons1.getQuality() < 1.0, true);
	TEST_EQUAL(cons2.getQuality() < 1.0, true);
	// quality(feat1, feat2) > quality((feat1, feat2), feat3):
	TEST_EQUAL(cons1.getQuality() > cons2.getQuality(), true);
	input[0].clear();
	input[0].push_back(ConsensusFeature(1, feat2));
	spf.run(input, result);
	ConsensusFeature cons3 = result[0];
	// quality(feat2, feat3) > quality(feat1, feat2), feat3):
	TEST_EQUAL(cons3.getQuality() > cons2.getQuality(), true);
	input[0][0] = ConsensusFeature(0, feat1);
	spf.run(input, result);
	ConsensusFeature cons4 = result[0];
	// quality(feat1, feat3) < quality(feat1, feat2), feat3):

Ejemplo n.º 2

Archivo: ConsensusXMLFile_test.C Proyecto: aiche/open-ms-mirror

TEST_EQUAL(map.getUnassignedPeptideIdentifications()[1].getHits()[0].getSequence(), "G")
TEST_EQUAL(map.getUnassignedPeptideIdentifications()[1].getHits()[1].getSequence(), "H")

//features
TEST_EQUAL(map.size(), 6)
ConsensusFeature cons_feature = map[0];
TEST_REAL_SIMILAR(cons_feature.getRT(), 1273.27)
TEST_REAL_SIMILAR(cons_feature.getMZ(), 904.47)
TEST_REAL_SIMILAR(cons_feature.getIntensity(), 3.12539e+07)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().minPosition()[0], 1273.27)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().maxPosition()[0], 1273.27)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().minPosition()[1], 904.47)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().maxPosition()[1], 904.47)
TEST_REAL_SIMILAR(cons_feature.getIntensityRange().minPosition()[0], 3.12539e+07)
TEST_REAL_SIMILAR(cons_feature.getIntensityRange().maxPosition()[0], 3.12539e+07)
TEST_REAL_SIMILAR(cons_feature.getQuality(), 1.1)
TEST_EQUAL(cons_feature.getMetaValue("peptide_id") == DataValue("RefSeq:NC_1234"), true)
ConsensusFeature::HandleSetType::const_iterator it = cons_feature.begin();
TEST_REAL_SIMILAR(it->getIntensity(), 3.12539e+07)

cons_feature = map[5];
TEST_REAL_SIMILAR(cons_feature.getRT(), 1194.82)
TEST_REAL_SIMILAR(cons_feature.getMZ(), 777.101)
TEST_REAL_SIMILAR(cons_feature.getIntensity(), 1.78215e+07)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().minPosition()[0], 1194.82)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().maxPosition()[0], 1194.82)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().minPosition()[1], 777.101)
TEST_REAL_SIMILAR(cons_feature.getPositionRange().maxPosition()[1], 777.101)
TEST_REAL_SIMILAR(cons_feature.getIntensityRange().minPosition()[0], 1.78215e+07)
TEST_REAL_SIMILAR(cons_feature.getIntensityRange().maxPosition()[0], 1.78215e+07)
TEST_REAL_SIMILAR(cons_feature.getQuality(), 0.0)