56 const int MS2_INDEX_0 = 1;
57 const int MS2_INDEX_1 = 3;
62 SpectrumPtr s21 = centroidedPtr->spectrum(MS2_INDEX_1,
true);
65 s21->binaryDataArrayPtrs.clear();
70 newMzs = vector<double>({ 0.0, 2.0, 2.000001, 3.999999, 4.0, 4.000001, 6.0, 8.0, 10.0, 12.0, 14.0, 16.0, 18.0 });
71 for (
size_t mz = 0; mz < newMzs.
size(); ++mz)
76 vector<double> s21ExpectedIntensities = { 1.0, 2.0, 3.0 };
77 while (s21ExpectedIntensities.size() < s20->getIntensityArray()->data.size())
79 s21ExpectedIntensities.push_back(1.0);
92 signal.reset(
new MatrixType(numSpectra, mzsToDemux->data.size()));
95 peakExtractor(s20, *signal, 0);
96 peakExtractor(s21, *signal, 1);
101 for (
size_t i = 0; i < baseIntensities->data.size(); ++i)
107 for (
size_t i = 0; i < s21ExpectedIntensities.size(); ++i)
117 signal.reset(
new MatrixType(numSpectra, s21->getMZArray()->data.size()));
118 binExamplePeakExtractor(s21, *signal, 0);
119 binExamplePeakExtractor(s20, *signal, 1);
122 for (
size_t i = 0; i < s21->getIntensityArray()->data.size(); ++i)
124 unit_assert_equal(signal->row(0)[i], s21->getIntensityArray()->data.at(i), 0.0001);