SpectrumList_SorterTest.cpp File Reference

#include "SpectrumList_Sorter.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"
#include <cstring>
#include "pwiz/data/msdata/examples.hpp"
#include "pwiz/data/msdata/TextWriter.hpp"
#include "boost/logic/tribool.hpp"

Go to the source code of this file.

Classes
struct	DefaultArrayLengthSorter
struct	MSLevelSorter

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

Variables
ostream *	os_ = 0

---

Function Documentation

void test

(

)

Definition at line 79 of file SpectrumList_SorterTest.cpp.

References pwiz::identdata::examples::initializeTiny(), os_, pwiz::msdata::MSData::run, pwiz::msdata::Run::spectrumListPtr, unit_assert, unit_assert_operator_equal, and write().

{
    MSData msd;
    examples::initializeTiny(msd);

    SpectrumListPtr originalList = msd.run.spectrumListPtr;

    SpectrumListPtr defaultArrayLengthSortedList(
        new SpectrumList_Sorter(originalList, DefaultArrayLengthSorter()));

    SpectrumListPtr msLevelUnstableSortedList(
        new SpectrumList_Sorter(originalList, MSLevelSorter()));

    SpectrumListPtr msLevelStableSortedList(
        new SpectrumList_Sorter(originalList, MSLevelSorter(), true));

    SpectrumListPtr sillySortedList(
        new SpectrumList_Sorter(msLevelStableSortedList, DefaultArrayLengthSorter()));

    if (os_)
    {
        *os_ << "Original spectrum list (" << originalList->size() << "):\n";
        TextWriter write(*os_);
        write(*originalList);
        *os_ << endl;
    }

    if (os_)
    {
        *os_ << "Default array length sorted spectrum list (" << defaultArrayLengthSortedList->size() << "):\n";
        TextWriter write(*os_);
        write(*defaultArrayLengthSortedList);
        *os_ << endl;
    }

    if (os_)
    {
        *os_ << "MS level unstable sorted spectrum list (" << msLevelUnstableSortedList->size() << "):\n";
        TextWriter write(*os_);
        write(*msLevelUnstableSortedList);
        *os_ << endl;
    }

    if (os_)
    {
        *os_ << "MS level stable sorted spectrum list (" << msLevelStableSortedList->size() << "):\n";
        TextWriter write(*os_);
        write(*msLevelStableSortedList);
        *os_ << endl;
    }

    if (os_)
    {
        *os_ << "Silly (nested) sorted spectrum list (" << sillySortedList->size() << "):\n";
        TextWriter write(*os_);
        write(*sillySortedList);
        *os_ << endl;
    }

    unit_assert_operator_equal(originalList->size(), defaultArrayLengthSortedList->size());
    unit_assert_operator_equal(originalList->size(), msLevelUnstableSortedList->size());
    unit_assert_operator_equal(originalList->size(), msLevelStableSortedList->size());
    unit_assert_operator_equal(originalList->size(), sillySortedList->size());

    SpectrumPtr s;

    // assert that the original list is unmodified
    unit_assert_operator_equal("scan=19", originalList->spectrumIdentity(0).id);
    unit_assert_operator_equal(0, originalList->spectrumIdentity(0).index);
    unit_assert_operator_equal("scan=20", originalList->spectrumIdentity(1).id);
    unit_assert_operator_equal(1, originalList->spectrumIdentity(1).index);
    unit_assert_operator_equal("scan=21", originalList->spectrumIdentity(2).id);
    unit_assert_operator_equal(2, originalList->spectrumIdentity(2).index);
    unit_assert_operator_equal("scan=22", originalList->spectrumIdentity(3).id);
    unit_assert_operator_equal(3, originalList->spectrumIdentity(3).index);
    s = originalList->spectrum(0);
    unit_assert_operator_equal("scan=19", s->id);
    unit_assert_operator_equal(0, s->index);

    // validate the default array length sorted list (ascending order, scan=19 and scan=22 are interchangeable)
    unit_assert_operator_equal("scan=21", defaultArrayLengthSortedList->spectrumIdentity(0).id);
    unit_assert_operator_equal(0, defaultArrayLengthSortedList->spectrumIdentity(0).index);
    unit_assert_operator_equal("scan=20", defaultArrayLengthSortedList->spectrumIdentity(1).id);
    unit_assert_operator_equal(1, defaultArrayLengthSortedList->spectrumIdentity(1).index);
    unit_assert_operator_equal(2, defaultArrayLengthSortedList->spectrumIdentity(2).index);
    unit_assert_operator_equal(3, defaultArrayLengthSortedList->spectrumIdentity(3).index);
    s = defaultArrayLengthSortedList->spectrum(0);
    unit_assert_operator_equal("scan=21", s->id);
    unit_assert_operator_equal(0, s->index);
    s = defaultArrayLengthSortedList->spectrum(1);
    unit_assert_operator_equal("scan=20", s->id);
    unit_assert_operator_equal(1, s->index);
    s = defaultArrayLengthSortedList->spectrum(2);
    unit_assert_operator_equal(2, s->index);
    s = defaultArrayLengthSortedList->spectrum(3);
    unit_assert_operator_equal(3, s->index);
    for (size_t i=1, end=defaultArrayLengthSortedList->size(); i < end; ++i)
        unit_assert(defaultArrayLengthSortedList->spectrum(i)->defaultArrayLength >=
                    defaultArrayLengthSortedList->spectrum(i-1)->defaultArrayLength);

    // validate the MS level unstable sorted list (scan=19, scan=21, and scan=22 are interchangeable)
    unit_assert_operator_equal(0, msLevelUnstableSortedList->spectrumIdentity(0).index);
    unit_assert_operator_equal(1, msLevelUnstableSortedList->spectrumIdentity(1).index);
    unit_assert_operator_equal(2, msLevelUnstableSortedList->spectrumIdentity(2).index);
    unit_assert_operator_equal("scan=20", msLevelUnstableSortedList->spectrumIdentity(3).id);
    unit_assert_operator_equal(3, msLevelUnstableSortedList->spectrumIdentity(3).index);
    s = msLevelUnstableSortedList->spectrum(0);
    unit_assert_operator_equal(0, s->index);
    s = msLevelUnstableSortedList->spectrum(1);
    unit_assert_operator_equal(1, s->index);
    s = msLevelUnstableSortedList->spectrum(2);
    unit_assert_operator_equal(2, s->index);
    s = msLevelUnstableSortedList->spectrum(3);
    unit_assert_operator_equal("scan=20", s->id);
    unit_assert_operator_equal(3, s->index);

    // validate the MS level stable sorted list (scan=19, scan=21, and scan=22 should stay in order)
    unit_assert_operator_equal("scan=19", msLevelStableSortedList->spectrumIdentity(0).id);
    unit_assert_operator_equal(0, msLevelStableSortedList->spectrumIdentity(0).index);
    unit_assert_operator_equal("scan=21", msLevelStableSortedList->spectrumIdentity(1).id);
    unit_assert_operator_equal(1, msLevelStableSortedList->spectrumIdentity(1).index);
    unit_assert_operator_equal("sample=1 period=1 cycle=23 experiment=1", msLevelStableSortedList->spectrumIdentity(2).id);
    unit_assert_operator_equal(2, msLevelStableSortedList->spectrumIdentity(2).index);
    unit_assert_operator_equal("scan=20", msLevelStableSortedList->spectrumIdentity(3).id);
    unit_assert_operator_equal(3, msLevelStableSortedList->spectrumIdentity(3).index);
    s = msLevelStableSortedList->spectrum(0);
    unit_assert_operator_equal("scan=19", s->id);
    unit_assert_operator_equal(0, s->index);
    s = msLevelStableSortedList->spectrum(1);
    unit_assert_operator_equal("scan=21", s->id);
    unit_assert_operator_equal(1, s->index);
    s = msLevelStableSortedList->spectrum(2);
    unit_assert_operator_equal("sample=1 period=1 cycle=23 experiment=1", s->id);
    unit_assert_operator_equal(2, s->index);
    s = msLevelStableSortedList->spectrum(3);
    unit_assert_operator_equal("scan=20", s->id);
    unit_assert_operator_equal(3, s->index);

    // validate the silly (nested) sorted list
    unit_assert_operator_equal("scan=21", sillySortedList->spectrumIdentity(0).id);
    unit_assert_operator_equal(0, sillySortedList->spectrumIdentity(0).index);
    unit_assert_operator_equal("scan=20", sillySortedList->spectrumIdentity(1).id);
    unit_assert_operator_equal(1, sillySortedList->spectrumIdentity(1).index);
    unit_assert_operator_equal(2, sillySortedList->spectrumIdentity(2).index);
    unit_assert_operator_equal(3, sillySortedList->spectrumIdentity(3).index);
    s = sillySortedList->spectrum(0);
    unit_assert_operator_equal("scan=21", s->id);
    unit_assert_operator_equal(0, s->index);
    s = sillySortedList->spectrum(1);
    unit_assert_operator_equal("scan=20", s->id);
    unit_assert_operator_equal(1, s->index);
    s = sillySortedList->spectrum(2);
    unit_assert_operator_equal(2, s->index);
    s = sillySortedList->spectrum(3);
    unit_assert_operator_equal(3, s->index);
    for (size_t i=1, end=sillySortedList->size(); i < end; ++i)
        unit_assert(sillySortedList->spectrum(i)->defaultArrayLength >=
                    sillySortedList->spectrum(i-1)->defaultArrayLength);
}

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

Definition at line 240 of file SpectrumList_SorterTest.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 39 of file SpectrumList_SorterTest.cpp.