|
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
|
#include "CommonTools/BaseParticlePropagator/interface/makeMuon.h"
#include "CommonTools/BaseParticlePropagator/interface/RawParticle.h"
#define CATCH_CONFIG_MAIN
#include "catch.hpp"
static constexpr const double kMuonMass = 0.10566;
TEST_CASE("makeMuon tests", "[makeMuon]") {
SECTION("muon at rest") {
auto m1 = rawparticle::makeMuon(true, math::XYZTLorentzVector{}, math::XYZTLorentzVector{});
REQUIRE(m1.charge() == -1.);
REQUIRE(m1.mass() == kMuonMass);
REQUIRE(m1.px() == 0.);
REQUIRE(m1.py() == 0.);
REQUIRE(m1.pz() == 0.);
REQUIRE(m1.x() == 0.);
REQUIRE(m1.y() == 0.);
REQUIRE(m1.z() == 0.);
REQUIRE(m1.t() == 0.);
//NOTE: energy is incorrectly calculated!
//REQUIRE(m1.e() == kMuonMass);
REQUIRE(m1.e() == 0.);
}
SECTION("anti-muon at rest") {
auto m1 = rawparticle::makeMuon(false, math::XYZTLorentzVector{}, math::XYZTLorentzVector{});
REQUIRE(m1.charge() == 1.);
REQUIRE(m1.mass() == kMuonMass);
REQUIRE(m1.px() == 0.);
REQUIRE(m1.py() == 0.);
REQUIRE(m1.pz() == 0.);
REQUIRE(m1.x() == 0.);
REQUIRE(m1.y() == 0.);
REQUIRE(m1.z() == 0.);
REQUIRE(m1.t() == 0.);
//NOTE: energy is incorrectly calculated!
//REQUIRE(m1.e() == kMuonMass);
REQUIRE(m1.e() == 0.);
}
SECTION("muon at rest, transposed") {
auto m1 = rawparticle::makeMuon(true, math::XYZTLorentzVector{}, math::XYZTLorentzVector{1., 2., 3., 4.});
REQUIRE(m1.charge() == -1.);
REQUIRE(m1.mass() == kMuonMass);
REQUIRE(m1.px() == 0.);
REQUIRE(m1.py() == 0.);
REQUIRE(m1.pz() == 0.);
REQUIRE(m1.x() == 1.);
REQUIRE(m1.y() == 2.);
REQUIRE(m1.z() == 3.);
REQUIRE(m1.t() == 4.);
//NOTE: energy is incorrectly calculated!
//REQUIRE(m1.e() == kMuonMass);
REQUIRE(m1.e() == 0.);
}
SECTION("muon in motion") {
auto m1 = rawparticle::makeMuon(true, math::XYZTLorentzVector{1., 2., 3., 8.}, math::XYZTLorentzVector{});
REQUIRE(m1.charge() == -1.);
REQUIRE(m1.mass() == kMuonMass);
REQUIRE(m1.px() == 1.);
REQUIRE(m1.py() == 2.);
REQUIRE(m1.pz() == 3.);
REQUIRE(m1.x() == 0.);
REQUIRE(m1.y() == 0.);
REQUIRE(m1.z() == 0.);
REQUIRE(m1.t() == 0.);
//NOTE: energy is incorrectly calculated!
//REQUIRE(m1.e() == sqrt(kMuonMass*kMuonMass+m1.px()*m1.px()+m1.py()*m1.py()+m1.pz()*m1.pz()));
REQUIRE(m1.e() == 8.);
}
}
|