|
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
|
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexSeed.h"
#include "catch.hpp"
static constexpr auto s_tag = "[PFDisplacedVertexSeed]";
TEST_CASE("Check adding elements", s_tag) {
reco::PFDisplacedVertexSeed seed;
REQUIRE(seed.elements().empty());
SECTION("updateSeedPoint") {
//empty tracks are fine
std::vector<reco::Track> tracks(5);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 1)));
REQUIRE(seed.elements().size() == 2);
REQUIRE(seed.nTracks() == 2);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 1)));
REQUIRE(seed.elements().size() == 2);
REQUIRE(seed.nTracks() == 2);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 2)));
REQUIRE(seed.elements().size() == 3);
REQUIRE(seed.nTracks() == 3);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 3)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 4)));
REQUIRE(seed.elements().size() == 5);
REQUIRE(seed.nTracks() == 5);
}
SECTION("addElement") {
//empty tracks are fine
std::vector<reco::Track> tracks(3);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 1)));
REQUIRE(seed.elements().size() == 2);
REQUIRE(seed.nTracks() == 2);
seed.addElement(reco::TrackBaseRef(reco::TrackRef(&tracks, 0)));
REQUIRE(seed.elements().size() == 2);
REQUIRE(seed.nTracks() == 2);
seed.addElement(reco::TrackBaseRef(reco::TrackRef(&tracks, 2)));
REQUIRE(seed.elements().size() == 3);
REQUIRE(seed.nTracks() == 3);
}
SECTION("mergeWith") {
std::vector<reco::Track> tracks(5);
seed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 1)));
SECTION("completely overlapping seeds") {
reco::PFDisplacedVertexSeed otherSeed;
otherSeed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 1)));
seed.mergeWith(otherSeed);
REQUIRE(seed.elements().size() == 2);
}
SECTION("partially overlapping seeds") {
reco::PFDisplacedVertexSeed otherSeed;
otherSeed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 0)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 2)));
seed.mergeWith(otherSeed);
REQUIRE(seed.elements().size() == 3);
}
SECTION("non overlapping seeds") {
REQUIRE(seed.elements().size() == 2);
reco::PFDisplacedVertexSeed otherSeed;
otherSeed.updateSeedPoint(GlobalPoint(0.01, 0.01, 0.01),
reco::TrackBaseRef(reco::TrackRef(&tracks, 3)),
reco::TrackBaseRef(reco::TrackRef(&tracks, 2)));
REQUIRE(otherSeed.elements().size() == 2);
seed.mergeWith(otherSeed);
REQUIRE(seed.elements().size() == 4);
}
}
}
|