PeakFamilyDetectorFTTest.cpp

Go to the documentation of this file.

//
// $Id: PeakFamilyDetectorFTTest.cpp 4129 2012-11-20 00:05:37Z chambm $
//
//
// Original author: Darren Kessner <darren@proteowizard.org>
//
// Copyright 2008 Spielberg Family Center for Applied Proteomics
//   Cedars-Sinai Medical Center, Los Angeles, California  90048
//
// Licensed under the Apache License, Version 2.0 (the "License"); 
// you may not use this file except in compliance with the License. 
// You may obtain a copy of the License at 
//
// http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software 
// distributed under the License is distributed on an "AS IS" BASIS, 
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 
// See the License for the specific language governing permissions and 
// limitations under the License.
//


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

using namespace pwiz::util;
using namespace pwiz::analysis;
using namespace pwiz::data;
using namespace pwiz::data::peakdata;
using namespace pwiz;
using namespace pwiz::msdata;


ostream* os_ = 0;


extern double peakFamilyDetectorFTTestData_[];
extern int peakFamilyDetectorFTTestDataSize_;


void test()
{
    // instantiate PeakFamilyDetectorFT

    PeakFamilyDetectorFT::Config config;
    config.log = os_;
    config.cp = CalibrationParameters::thermo_FT();
    PeakFamilyDetectorFT detector(config); 

    // detect 

    vector<PeakFamily> result;
    const MZIntensityPair* begin = 
        reinterpret_cast<const MZIntensityPair*>(&peakFamilyDetectorFTTestData_[0]);
    const MZIntensityPair* end = begin + peakFamilyDetectorFTTestDataSize_/2;

    detector.detect(begin, end, result);

    if (os_)
    {
        *os_ << setprecision(10) << "result: " << result.size() << endl;
        copy(result.begin(), result.end(), ostream_iterator<PeakFamily>(*os_, "\n"));
    }

    unit_assert(result.size() == 1);
    unit_assert_equal(result[0].mzMonoisotopic, 810.4148, .005);
}


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

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

    TEST_EPILOG
}