AminoAcidTest.cpp File Reference

#include "pwiz/utility/misc/unit.hpp"
#include "AminoAcid.hpp"
#include <cstring>
#include "pwiz/utility/misc/Std.hpp"
#include "boost/thread/thread.hpp"
#include "boost/thread/barrier.hpp"

Go to the source code of this file.

Classes
struct	TestAminoAcid

Functions
bool	hasLowerMass (const AminoAcid::Info::Record &a, const AminoAcid::Info::Record &b)
void	printRecord (ostream *os, const AminoAcid::Info::Record &record)
void	test ()
void	testThreadSafetyWorker (boost::barrier *testBarrier)
void	testThreadSafety (const int &testThreadCount)
int	main (int argc, char *argv[])

Variables
ostream *	os_ = 0
TestAminoAcid	testAminoAcids []

---

Function Documentation

bool hasLowerMass	(	const AminoAcid::Info::Record &	a,
		const AminoAcid::Info::Record &	b
	)

Definition at line 42 of file AminoAcidTest.cpp.

References pwiz::proteome::AminoAcid::Info::Record::formula, and pwiz::chemistry::Formula::monoisotopicMass().

Referenced by test().

{
    return a.formula.monoisotopicMass() < b.formula.monoisotopicMass();
}

void printRecord	(	ostream *	os,
		const AminoAcid::Info::Record &	record
	)

Definition at line 48 of file AminoAcidTest.cpp.

References pwiz::proteome::AminoAcid::Info::Record::abundance, pwiz::proteome::AminoAcid::Info::Record::formula, pwiz::chemistry::Formula::molecularWeight(), pwiz::chemistry::Formula::monoisotopicMass(), pwiz::proteome::AminoAcid::Info::Record::name, pwiz::proteome::AminoAcid::Info::Record::residueFormula, and pwiz::proteome::AminoAcid::Info::Record::symbol.

Referenced by test().

{
    if (!os) return;

    *os << record.symbol << ": " 
        << setw(14) << record.name << " " 
        << setw(11) << record.formula << " " 
        << setprecision(3)
        << setw(5) << record.abundance << " "
        << fixed << setprecision(6)
        << setw(10) << record.formula.monoisotopicMass() << " "
        << setw(10) << record.formula.molecularWeight() << " "
        << setw(10) << record.residueFormula.monoisotopicMass() << " "
        << setw(10) << record.residueFormula.molecularWeight() << endl;
}

void test

(

)

Definition at line 102 of file AminoAcidTest.cpp.

References pwiz::proteome::AminoAcid::Info::Record::abundance, Alanine, TestAminoAcid::avgMass, C, H, hasLowerMass(), pwiz::chemistry::Formula::molecularWeight(), pwiz::chemistry::Formula::monoisotopicMass(), TestAminoAcid::monoMass, N, O, os_, printRecord(), pwiz::proteome::AminoAcid::Info::record(), pwiz::proteome::AminoAcid::Info::Record::residueFormula, S, Se, Selenocysteine, sort(), TestAminoAcid::symbol, testAminoAcids, unit_assert, and unit_assert_equal.

{
    // get a copy of all the records

    vector<AminoAcid::Info::Record> records;

    for (char symbol='A'; symbol<='Z'; symbol++)
    {
        try 
        {
            const AminoAcid::Info::Record& record = AminoAcid::Info::record(symbol);
            records.push_back(record);
        }
        catch (exception&)
        {}
    }

    for (vector<AminoAcid::Info::Record>::iterator it=records.begin(); it!=records.end(); ++it)
        printRecord(os_, *it);

    unit_assert(AminoAcid::Info::record(Alanine).residueFormula[C] == 3);
    unit_assert(AminoAcid::Info::record(Alanine).residueFormula[H] == 5);
    unit_assert(AminoAcid::Info::record(Alanine).residueFormula[N] == 1);
    unit_assert(AminoAcid::Info::record(Alanine).residueFormula[O] == 1);
    unit_assert(AminoAcid::Info::record(Alanine).residueFormula[S] == 0);

    unit_assert(AminoAcid::Info::record(Alanine).formula[C] == 3);
    unit_assert(AminoAcid::Info::record(Alanine).formula[H] == 7);
    unit_assert(AminoAcid::Info::record(Alanine).formula[N] == 1);
    unit_assert(AminoAcid::Info::record(Alanine).formula[O] == 2);
    unit_assert(AminoAcid::Info::record(Alanine).formula[S] == 0);

    unit_assert(AminoAcid::Info::record(Selenocysteine).formula[Se] == 1);

    // test single amino acids
    for (int i=0; i < 22; ++i) // skip X for now
    {
        TestAminoAcid& aa = testAminoAcids[i];
        Formula residueFormula = AminoAcid::Info::record(aa.symbol).residueFormula;
        unit_assert_equal(residueFormula.monoisotopicMass(), aa.monoMass, 0.00001);
        unit_assert_equal(residueFormula.molecularWeight(), aa.avgMass, 0.0001);
        //set<char> mmNames = mm2n.getNames(aa.monoMass, EPSILON);
        //set<char> amNames = am2n.getNames(aa.avgMass, EPSILON);
        //unit_assert(mmNames.count(aa.symbol) > 0);
        //unit_assert(amNames.count(aa.symbol) > 0);
    }


    // compute some averages

    double averageMonoisotopicMass = 0;
    double averageC = 0;
    double averageH = 0;
    double averageN = 0;
    double averageO = 0;
    double averageS = 0;

    for (vector<AminoAcid::Info::Record>::iterator it=records.begin(); it!=records.end(); ++it)
    {
        const AminoAcid::Info::Record& record = *it;

        Formula residueFormula = record.residueFormula;
        averageMonoisotopicMass += residueFormula.monoisotopicMass() * record.abundance; 
        averageC += residueFormula[C] * record.abundance;
        averageH += residueFormula[H] * record.abundance;
        averageN += residueFormula[N] * record.abundance;
        averageO += residueFormula[O] * record.abundance;
        averageS += residueFormula[S] * record.abundance;
    }

    if (os_) *os_ << setprecision(8) << endl;
    if (os_) *os_ << "average residue C: " << averageC << endl;    
    if (os_) *os_ << "average residue H: " << averageH << endl;    
    if (os_) *os_ << "average residue N: " << averageN << endl;    
    if (os_) *os_ << "average residue O: " << averageO << endl;    
    if (os_) *os_ << "average residue S: " << averageS << endl;    
    if (os_) *os_ << endl;

    if (os_) *os_ << "average monoisotopic mass: " << averageMonoisotopicMass << endl;
    double averageResidueMass = averageMonoisotopicMass;
    if (os_) *os_ << "average residue mass: " << averageResidueMass << endl << endl;

    // sort by monoisotopic mass and print again
    sort(records.begin(), records.end(), hasLowerMass);
    for (vector<AminoAcid::Info::Record>::iterator it=records.begin(); it!=records.end(); ++it)
        printRecord(os_, *it); 
}

void testThreadSafetyWorker

(

boost::barrier *

testBarrier

)

Definition at line 191 of file AminoAcidTest.cpp.

References e(), and test().

{
    testBarrier->wait(); // wait until all threads have started

    try
    {
        test();
    }
    catch (exception& e)
    {
        cerr << "Exception in worker thread: " << e.what() << endl;
    }
    catch (...)
    {
        cerr << "Unhandled exception in worker thread." << endl;
    }
}

void testThreadSafety

(

const int &

testThreadCount

)

Definition at line 209 of file AminoAcidTest.cpp.

References testThreadSafetyWorker().

{
    boost::barrier testBarrier(testThreadCount);
    boost::thread_group testThreadGroup;
    for (int i=0; i < testThreadCount; ++i)
        testThreadGroup.add_thread(new boost::thread(&testThreadSafetyWorker, &testBarrier));
    testThreadGroup.join_all();
}

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

Definition at line 219 of file AminoAcidTest.cpp.

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

{
    TEST_PROLOG(argc, argv)

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

    try
    {
        //testThreadSafety(1); // does not test thread-safety of singleton initialization
        testThreadSafety(2);
        testThreadSafety(4);
        testThreadSafety(8);
        testThreadSafety(16);
    }
    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 AminoAcidTest.cpp.

TestAminoAcid testAminoAcids[]

Initial value:

{
    { 71.0371140, 71.07790, 'A' },    
    { 156.101111, 156.1857, 'R' },    
    { 114.042927, 114.1026, 'N' },    
    { 115.026943, 115.0874, 'D' },    
    { 103.009185, 103.1429, 'C' },    
    { 129.042593, 129.1140, 'E' },    
    { 128.058578, 128.1292, 'Q' },    
    { 57.0214640, 57.05130, 'G' },    
    { 137.058912, 137.1393, 'H' },    
    { 113.084064, 113.1576, 'I' },    
    { 113.084064, 113.1576, 'L' },    
    { 128.094963, 128.1723, 'K' },    
    { 131.040485, 131.1961, 'M' },    
    { 147.068414, 147.1739, 'F' },    
    { 97.0527640, 97.11520, 'P' },    
    { 87.0320280, 87.07730, 'S' },    
    { 150.953636, 150.0379, 'U' },    
    { 101.047679, 101.1039, 'T' },    
    { 186.079313, 186.2099, 'W' },    
    { 163.063329, 163.1733, 'Y' },    
    { 99.0684140, 99.13110, 'V' },    
    { 114.042927, 114.1026, 'B' },    
    { 128.058578, 128.1292, 'Z' },    
    { 114.091900, 114.1674, 'X' }     
}

Definition at line 73 of file AminoAcidTest.cpp.

Referenced by test().