-
Notifications
You must be signed in to change notification settings - Fork 494
Expand file tree
/
Copy pathMCTrack.cxx
More file actions
95 lines (85 loc) · 2.73 KB
/
MCTrack.cxx
File metadata and controls
95 lines (85 loc) · 2.73 KB
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
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
#define BOOST_TEST_MODULE Test MCTrack class
#define BOOST_TEST_MAIN
#define BOOST_TEST_DYN_LINK
#include "SimulationDataFormat/MCTrack.h"
#include <boost/test/unit_test.hpp>
#include <iomanip>
#include <ios>
#include <iostream>
#include "DetectorsCommonDataFormats/DetID.h"
#include "SimulationDataFormat/Stack.h"
#include "SimulationDataFormat/PrimaryChunk.h"
#include "TFile.h"
#include "TParticle.h"
#include "TMCProcess.h"
using namespace o2;
BOOST_AUTO_TEST_CASE(MCTrack_test)
{
MCTrack track;
// check properties on default constructed object
BOOST_CHECK(track.getStore() == false);
for (auto i = o2::detectors::DetID::First; i < o2::detectors::DetID::nDetectors; ++i) {
BOOST_CHECK(track.leftTrace(i) == false);
}
BOOST_CHECK(track.getNumDet() == 0);
BOOST_CHECK(track.hasHits() == false);
// check storing
track.setStore(true);
BOOST_CHECK(track.getStore() == true);
track.setStore(false);
BOOST_CHECK(track.getStore() == false);
track.setStore(true);
BOOST_CHECK(track.getStore() == true);
// set hit for first detector
BOOST_CHECK(track.leftTrace(1) == false);
track.setHit(1);
BOOST_CHECK(track.hasHits() == true);
BOOST_CHECK(track.leftTrace(1) == true);
BOOST_CHECK(track.getNumDet() == 1);
{
// serialize it
TFile f("MCTrackOut.root", "RECREATE");
f.WriteObject(&track, "MCTrack");
f.Close();
}
{
MCTrack* intrack = nullptr;
TFile f("MCTrackOut.root", "OPEN");
f.GetObject("MCTrack", intrack);
BOOST_CHECK(intrack->getStore() == true);
BOOST_CHECK(intrack->hasHits() == true);
}
}
// unit tests on stack
BOOST_AUTO_TEST_CASE(Stack_test)
{
o2::data::Stack st;
int a;
TMCProcess proc{kPPrimary};
// add a 2 primary particles
st.PushTrack(1, -1, 0, 0, 0., 0., 10., 5., 5., 5., 0.1, 0., 0., 0., proc, a, 1., 1);
st.PushTrack(1, -1, 0, 0, 0., 0., 10., 5., 5., 5., 0.1, 0., 0., 0., proc, a, 1., 1);
BOOST_CHECK(st.getPrimaries().size() == 2);
{
// serialize it
TFile f("StackOut.root", "RECREATE");
f.WriteObject(&st, "Stack");
f.Close();
}
{
o2::data::Stack* inst = nullptr;
TFile f("StackOut.root", "OPEN");
f.GetObject("Stack", inst);
BOOST_CHECK(inst->getPrimaries().size() == 2);
}
}