File: test_big_hyperscore.cpp

package info (click to toggle)
libpappsomspp 0.11.9-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 18,580 kB
  • sloc: cpp: 82,463; xml: 44,164; python: 668; sql: 186; sh: 33; makefile: 32
file content (111 lines) | stat: -rw-r--r-- 3,490 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
//
// File: test_big_hyperscore.cpp
// Created by: Olivier Langella
// Created on: 19/3/2015
//
/*******************************************************************************
 * Copyright (c) 2015 Olivier Langella <Olivier.Langella@moulon.inra.fr>.
 *
 * This file is part of the PAPPSOms++ library.
 *
 *     PAPPSOms++ is free software: you can redistribute it and/or modify
 *     it under the terms of the GNU General Public License as published by
 *     the Free Software Foundation, either version 3 of the License, or
 *     (at your option) any later version.
 *
 *     PAPPSOms++ is distributed in the hope that it will be useful,
 *     but WITHOUT ANY WARRANTY; without even the implied warranty of
 *     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *     GNU General Public License for more details.
 *
 *     You should have received a copy of the GNU General Public License
 *     along with PAPPSOms++.  If not, see <http://www.gnu.org/licenses/>.
 *
 * Contributors:
 *     Olivier Langella <Olivier.Langella@moulon.inra.fr> - initial API and
 *implementation
 ******************************************************************************/


// ./tests/catch2-only-tests [bighyperscore] -s

#ifdef CATCH2_MAJOR_VERSION_2
#include <catch2/catch.hpp>
using namespace Catch;
#elif CATCH2_MAJOR_VERSION_3
#include <catch2/catch_all.hpp>
using namespace Catch;
using namespace Catch::Matchers;
#endif


#include <pappsomspp/core/mzrange.h>
#include <pappsomspp/core/peptide/peptidenaturalisotopelist.h>
#include <pappsomspp/core/peptide/peptidefragmentionlistbase.h>
#include <pappsomspp/core/massspectrum/massspectrum.h>
#include <pappsomspp/core/psm/xtandem/xtandemhyperscore.h>
#include <pappsomspp/core/pappsoexception.h>
#include <iostream>
#include <iomanip>
#include <cmath>
#include <set>
#include <QDebug>
#include <QString>
#include <QFileInfo>
#include "config.h"
#include "saxparsers/xtandemresultshandler.h"

using namespace pappso;
using namespace std;
// using namespace pwiz::msdata;


TEST_CASE("bighyperscore test suite.", "[bighyperscore]")
{
  // Set the debugging message formatting pattern.
  qSetMessagePattern(QString("%{file}@%{line}, %{function}(): %{message}"));

  SECTION("..:: test bighyperscore ::..", "[bighyperscore]")
  {

#if USEPAPPSOTREE == 1
    // return 1;
    std::cout
      << std::endl
      << "..:: test all identified peptides in an XML tandem result file ::.."
      << std::endl;
    PrecisionPtr precision = PrecisionFactory::getDaltonInstance(0.02);

    // bool refine_spectrum_synthesis = false;
    QFileInfo fileinfo(
      "/gorgone/pappso/formation/Janvier2014/TD/xml_tandem/"
      "20120906_balliau_extract_1_A01_urnb-1.xml");

    XtandemResultsHandler *parser = new XtandemResultsHandler(precision);

    QXmlSimpleReader simplereader;
    simplereader.setContentHandler(parser);
    simplereader.setErrorHandler(parser);

    qDebug() << "Read tandem XML result file '" << fileinfo.filePath() << "'";

    QFile qfile(fileinfo.absoluteFilePath());
    QXmlInputSource xmlInputSource(&qfile);

    if(simplereader.parse(xmlInputSource))
      {
      }
    else
      {
        qDebug() << parser->errorString();
        // throw PappsoException(
        //    QObject::tr("error reading tandem XML result file :\n").append(
        //         parser->errorString()));
      }
    qfile.close();

#endif
    // big hyperscore test :
    // check every X!Tandem match in an XML result file
  }
}