PeakFitterTest.cpp File Reference

#include "PeakFitter.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	testParabola ()
void	testMultiplePeaks ()
void	test ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0

---

Function Documentation

void testParabola

(

)

Definition at line 39 of file PeakFitterTest.cpp.

References pwiz::analysis::Peak::area, e(), epsilon, pwiz::analysis::PeakFitter_Parabola::fitPeak(), os_, unit_assert, unit_assert_equal, and pwiz::analysis::PeakFitter_Parabola::Config::windowRadius.

Referenced by test().

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

    const double center = 2.1;
    const double height = 5;

    vector<OrderedPair> pairs;
    for (double i=0; i<5; i++)
        pairs.push_back(OrderedPair(i, height-(i-center)*(i-center))); // sampled parabola

    PeakFitter_Parabola::Config config;
    config.windowRadius = 1;

    PeakFitter_Parabola fitter(config);
    Peak peak;

    fitter.fitPeak(pairs, 2, peak);

    if (os_) 
    {
        *os_ << peak;
        copy(peak.data.begin(), peak.data.end(), ostream_iterator<OrderedPair>(*os_, " "));
        *os_ << endl;
    }

    const double epsilon = 1e-6;
    unit_assert_equal(peak.mz, center, epsilon);
    unit_assert_equal(peak.intensity, 5, epsilon);
    unit_assert_equal(peak.area, 12.97, epsilon);
    unit_assert_equal(peak.error, 0, epsilon);
    unit_assert_equal(peak.intensity, 5, epsilon);
    unit_assert(!peak.data.empty());
}

void testMultiplePeaks

(

)

Definition at line 75 of file PeakFitterTest.cpp.

References e(), epsilon, pwiz::analysis::PeakFitter::fitPeaks(), pwiz::chemistry::Ion::mz(), os_, unit_assert, and unit_assert_equal.

Referenced by test().

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

    const double center = 2.1;
    const double height = 5;

    vector<OrderedPair> pairs;
    for (double i=0; i<5; i++)
        pairs.push_back(OrderedPair(i, height-(i-center)*(i-center))); // sampled parabola
    for (double i=0; i<5; i++)
        pairs.push_back(OrderedPair(i+5, height-(i-center)*(i-center))); // sampled parabola
    for (double i=0; i<5; i++)
        pairs.push_back(OrderedPair(i+10, height-(i-center)*(i-center))); // sampled parabola

    vector<size_t> indices;
    indices.push_back(2);
    indices.push_back(7);
    indices.push_back(12);

    PeakFitter_Parabola fitter;
    vector<Peak> peaks;

    fitter.fitPeaks(pairs, indices, peaks);

    if (os_) 
    {
        for (vector<Peak>::const_iterator it=peaks.begin(); it!=peaks.end(); ++it)
        {
            *os_ << *it;
            copy(it->data.begin(), it->data.end(), ostream_iterator<OrderedPair>(*os_, " "));
            *os_ << endl;
        }
    }

    const double epsilon = 1e-6;
    unit_assert(peaks.size() == 3);
    unit_assert_equal(peaks[0].mz, center, epsilon);
    unit_assert_equal(peaks[1].mz, center+5, epsilon);
    unit_assert_equal(peaks[2].mz, center+10, epsilon);
}

void test

(

)

Definition at line 118 of file PeakFitterTest.cpp.

References testMultiplePeaks(), and testParabola().

{
    testParabola();
    testMultiplePeaks();
}

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

Definition at line 125 of file PeakFitterTest.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 PeakFitterTest.cpp.