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
|
# Copyright (c) 2019-2025, Saransh Chopra, Henry Schreiner, Eduardo Rodrigues, Jonas Eschle, and Jim Pivarski.
#
# Distributed under the 3-clause BSD license, see accompanying file LICENSE
# or https://github.com/scikit-hep/vector for details.
from __future__ import annotations
import pytest
import vector.backends.object
def test_spatial_object():
vec = vector.backends.object.VectorObject3D(
azimuthal=vector.backends.object.AzimuthalObjectXY(0.5, 0.6),
longitudinal=vector.backends.object.LongitudinalObjectZ(0.7),
)
out = vec.rotate_quaternion(0.1, 0.2, 0.3, 0.4)
assert isinstance(out, vector.backends.object.VectorObject3D)
assert isinstance(out.azimuthal, vector.backends.object.AzimuthalObjectXY)
assert isinstance(out.longitudinal, vector.backends.object.LongitudinalObjectZ)
assert out.x == pytest.approx(0.078)
assert out.y == pytest.approx(0.18)
assert out.z == pytest.approx(0.246)
for t in "xyz", "xytheta", "xyeta", "rhophiz", "rhophitheta", "rhophieta":
tvec = getattr(vec, "to_" + t)()
out = tvec.rotate_quaternion(0.1, 0.2, 0.3, 0.4)
assert isinstance(out, vector.backends.object.VectorObject3D)
assert isinstance(out.azimuthal, vector.backends.object.AzimuthalObjectXY)
assert isinstance(out.longitudinal, vector.backends.object.LongitudinalObjectZ)
assert out.x == pytest.approx(0.078)
assert out.y == pytest.approx(0.18)
assert out.z == pytest.approx(0.246)
|
