PeakExtractorTest.cpp File Reference

#include "PeakExtractor.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"
#include <cstring>

Go to the source code of this file.

Functions
void	test ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0
const double	testData_ []
const size_t	testDataSize_ = sizeof(testData_)/sizeof(double)

---

Function Documentation

void test

(

)

Definition at line 170 of file PeakExtractorTest.cpp.

References epsilon, pwiz::analysis::PeakExtractor::extractPeaks(), pwiz::chemistry::Ion::mz(), os_, testData_, testDataSize_, unit_assert, unit_assert_equal, pwiz::analysis::PeakFinder_SNR::Config::windowRadius, pwiz::analysis::PeakFitter_Parabola::Config::windowRadius, and pwiz::analysis::PeakFinder_SNR::Config::zValueThreshold.

{
    if (os_) *os_ << "test()\n";

    shared_ptr<NoiseCalculator> noiseCalculator(new NoiseCalculator_2Pass);

    PeakFinder_SNR::Config pfsnrConfig;
    pfsnrConfig.windowRadius = 2;
    pfsnrConfig.zValueThreshold = 2;

    shared_ptr<PeakFinder> peakFinder(new PeakFinder_SNR(noiseCalculator, pfsnrConfig));

    PeakFitter_Parabola::Config pfpConfig;
    pfpConfig.windowRadius = 1; // (windowRadius != 1) is not good for real data
    shared_ptr<PeakFitter> peakFitter(new PeakFitter_Parabola(pfpConfig));

    PeakExtractor peakExtractor(peakFinder, peakFitter);

    OrderedPairContainerRef data(testData_, testData_+testDataSize_);
    vector<Peak> peaks;
    peakExtractor.extractPeaks(data, peaks);

    if (os_)
    {
        *os_ << "peaks: " << peaks.size() << endl; 
        copy(peaks.begin(), peaks.end(), ostream_iterator<Peak>(*os_, "\n"));
    }
    
    const double epsilon = .01;
    unit_assert(peaks.size() == 3);
    unit_assert_equal(peaks[0].mz, 810.41, epsilon);
    unit_assert_equal(peaks[1].mz, 810.91, epsilon);
    unit_assert_equal(peaks[2].mz, 811.41, epsilon);
}

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 206 of file PeakExtractorTest.cpp.

References e(), os_, test(), TEST_EPILOG, TEST_FAILED, and TEST_PROLOG.

{
    TEST_PROLOG(argc, argv)

    try
    {
        if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
        test();
    }
    catch (exception& e)
    {
        TEST_FAILED(e.what())
    }
    catch (...)
    {
        TEST_FAILED("Caught unknown exception.")
    }

    TEST_EPILOG
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 36 of file PeakExtractorTest.cpp.

const double testData_[]

Definition at line 39 of file PeakExtractorTest.cpp.

const size_t testDataSize_ = sizeof(testData_)/sizeof(double)

Definition at line 167 of file PeakExtractorTest.cpp.