Feature_dataFetcherTest.cpp File Reference

#include "Feature_dataFetcher.hpp"
#include "pwiz/utility/misc/unit.hpp"
#include "pwiz/utility/misc/Std.hpp"

Go to the source code of this file.

Classes
struct	IsFS

Functions
FeaturePtr	makeFeature (double mz, double retentionTime)
FeatureSequencedPtr	makeFeatureSequenced (const FeaturePtr feature, string _ms1_5="", string _ms2="")
void	test ()
void	testMerge ()
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0

---

Function Documentation

FeaturePtr makeFeature	(	double	mz,
		double	retentionTime
	)

Definition at line 34 of file Feature_dataFetcherTest.cpp.

References pwiz::chemistry::Ion::mz().

{
    FeaturePtr feature(new Feature());
    feature->mz = mz;
    feature->retentionTime = retentionTime;

    return feature;

}

FeatureSequencedPtr makeFeatureSequenced	(	const FeaturePtr	feature,
		string	_ms1_5 = "",
		string	_ms2 = ""
	)

Definition at line 44 of file Feature_dataFetcherTest.cpp.

Referenced by testMerge().

{
    FeatureSequencedPtr fs(new FeatureSequenced());
    fs->feature = feature;
    fs->ms1_5 = _ms1_5;
    fs->ms2 = _ms2;

    return fs;

}

void test

(

)

Definition at line 63 of file Feature_dataFetcherTest.cpp.

References pwiz::data::peakdata::FeatureFile::features, pwiz::eharmony::Feature_dataFetcher::getFeatures(), makeFeature(), os_, unit_assert, and pwiz::data::peakdata::FeatureFile::write().

{
    
    if (os_) *os_ << "\ntest() ... \n\n";

    // make a vector of features
    FeaturePtr a = makeFeature(1,2);
    FeaturePtr b = makeFeature(3,4);
    FeaturePtr c = makeFeature(5,6);

    // make FeatureSequenced objects
    FeatureSequenced fs_a(a);
    FeatureSequenced fs_b(b);
    FeatureSequenced fs_c(c);

    vector<FeaturePtr> features;
    features.push_back(a);
    features.push_back(b);
    features.push_back(c);

    // test vector<FeaturePtr> constructor
    if (os_) *os_ << "constructing Feature_dataFetcher ... " << endl;
    Feature_dataFetcher fdf(features);
    if (os_) *os_ << "constructed. " << endl;

    if (os_) *os_ << "testing getFeatures ... " << endl;
    vector<boost::shared_ptr<FeatureSequenced> > test_a = fdf.getFeatures(1,2);
    vector<boost::shared_ptr<FeatureSequenced> > test_b = fdf.getFeatures(3,4);
    vector<boost::shared_ptr<FeatureSequenced> > test_c = fdf.getFeatures(5,6);

    unit_assert(find_if(test_a.begin(), test_a.end(), IsFS(fs_a)) != test_a.end());
    unit_assert(find_if(test_b.begin(), test_b.end(), IsFS(fs_b)) != test_b.end());
    unit_assert(find_if(test_c.begin(), test_c.end(), IsFS(fs_c)) != test_c.end());

    if (os_) 
        {
            *os_ << "testing vector<Feature> constructor ... \n";
            ostringstream oss;
            XMLWriter writer(oss);
            vector<boost::shared_ptr<FeatureSequenced> >::iterator a_it = test_a.begin();
            for(; a_it != test_a.end(); ++a_it) (*a_it)->feature->write(writer);
            *os_ << oss.str() << endl;

        }


    // write a FeatureFile
    ostringstream oss;
    XMLWriter writer(oss);
    FeatureFile ff;
    ff.features = features;
    ff.write(writer);
    
    //read it to an istream
    istringstream iss(oss.str());

    // test istream constructor
    Feature_dataFetcher fdf_is(iss);

    vector<boost::shared_ptr<FeatureSequenced> > test_is_a = fdf_is.getFeatures(1,2);
    vector<boost::shared_ptr<FeatureSequenced> > test_is_b = fdf_is.getFeatures(3,4);
    vector<boost::shared_ptr<FeatureSequenced> > test_is_c = fdf_is.getFeatures(5,6);

    unit_assert(find_if(test_is_a.begin(), test_is_a.end(), IsFS(fs_a)) != test_is_a.end());
    unit_assert(find_if(test_is_b.begin(), test_is_b.end(), IsFS(fs_b)) != test_is_b.end());
    unit_assert(find_if(test_is_c.begin(), test_is_c.end(), IsFS(fs_c)) != test_is_c.end());


    if (os_)
        {
            *os_ << "\ntesting istream constructor ... \n";
      
        }


    return;

}

void testMerge

(

)

Definition at line 142 of file Feature_dataFetcherTest.cpp.

References pwiz::eharmony::Feature_dataFetcher::getFeatures(), makeFeature(), makeFeatureSequenced(), pwiz::eharmony::Feature_dataFetcher::merge(), os_, and unit_assert.

Referenced by main().

{
    if(os_) *os_ << "\ntestMerge()...\n" << endl;

    FeaturePtr a = makeFeature(5,1);
    FeaturePtr b = makeFeature(20,9);

    FeatureSequencedPtr fs_b = makeFeatureSequenced(b);

    vector<FeaturePtr> v_a;
    vector<FeaturePtr> v_b;

    v_a.push_back(a);
    v_b.push_back(b);

    Feature_dataFetcher fdf_a(v_a);
    Feature_dataFetcher fdf_b(v_b);
    
    fdf_a.merge(fdf_b);
    vector<boost::shared_ptr<FeatureSequenced> > binContents = fdf_a.getFeatures(b->mz, b->retentionTime);
    
    unit_assert(binContents.size() > 0);

    if (os_)
        {
            *os_ << "Merged FeatureSequenced:\n " << endl;

            XMLWriter writer(*os_);
            (*binContents.begin())->feature->write(writer);
            *os_ << "ms1_5: " <<  (*binContents.begin())->ms1_5 << endl;
            *os_ << "ms2: " << (*binContents.begin())->ms2 << endl;

            *os_ << "\nOriginal FeatureSequenced:\n " << endl;   
             
            fs_b->feature->write(writer);
            *os_ << "ms1_5: " << fs_b->ms1_5 << endl;
            *os_ << "ms2: " << fs_b->ms2 << endl;

        }  

    unit_assert((*binContents.back()) == *fs_b);

}

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

Definition at line 186 of file Feature_dataFetcherTest.cpp.

References e(), os_, test(), and testMerge().

{

    try
        {
            if (argc>1 && !strcmp(argv[1],"-v")) os_ = &cout;
            if (os_) *os_ << "Feature_dataFetcherTest ... \n";

            test();
            testMerge();

            return 0;
        }
    catch (exception& e)
        {
            cerr << e.what() << endl;
        }
    catch (...)
        {
            cerr << "Caught unknown exception.\n";
        }

    return 1;
}

---

Variable Documentation

ostream* os_ = 0

Definition at line 32 of file Feature_dataFetcherTest.cpp.