File: test_tims_dda.cpp

package info (click to toggle)
libpappsomspp 0.11.10-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 18,844 kB
  • sloc: cpp: 85,012; xml: 44,164; python: 668; sql: 186; sh: 33; makefile: 32
file content (82 lines) | stat: -rw-r--r-- 2,771 bytes parent folder | download
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
//
// File: test_tims_dda.cpp
// Created by: Olivier Langella
// Created on: 6/9/2024
//
/*******************************************************************************
 * Copyright (c) 2024 Olivier Langella
 *<Olivier.Langella@universite-paris-saclay.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/>.
 *
 ******************************************************************************/


// ./tests/pappsotree/catch2-only-pappsotree [dda] -s


#include <catch2/catch_test_macros.hpp>
#include <catch2/matchers/catch_matchers_vector.hpp>
#include "../common.h"
#include "tests-config.h"


#include <QDebug>
#include <QCoreApplication>
#include <pappsomspp/core/exception/exceptionnotfound.h>
#include <pappsomspp/core/exception/exceptionoutofrange.h>


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

  int argc(0);
  char *argv[]          = {};
  QCoreApplication *app = new QCoreApplication(argc, argv);

  SECTION("..:: test DIA file access ::..", "[dda]")
  {

    pappso::MsFileAccessor dda_accessor(
      "/gorgone/pappso/jouy/raw/2024_timsTOF_Pro/20240826_30_Zuzana/8-30-2024_Gal_3_pool_9501.d",
      "a1");


    dda_accessor.setPreferredFileReaderType(Enums::MsDataFormat::brukerTims,
                                            Enums::FileReaderType::tims_ms2);
    pappso::MsRunReaderSPtr p_tims_reader =
      dda_accessor.msRunReaderSPtr(dda_accessor.getMsRunIds().front());

    REQUIRE(
      p_tims_reader.get()->getOboPsiModTermInstrumentModelName().getAccession().toStdString() ==
      "MS:1003005");
    REQUIRE(p_tims_reader.get()->getOboPsiModTermInstrumentModelName().isA("MS:1003123"));

    REQUIRE(p_tims_reader.get() != nullptr);

    pappso::MsRunSimpleStatistics stats;
    pappso::MsRunReadConfig config;
    config.setNeedPeakList(false);
    p_tims_reader.get()->readSpectrumCollection2(config, stats);

    REQUIRE(stats.getMsLevelCount(1) == 98511);
    REQUIRE(stats.getMsLevelCount(2) == 98511);
  }

  delete app;
}