Back to home page

EIC code displayed by LXR

 
 

    


File indexing completed on 2025-10-29 07:55:39

0001 // This file is part of the ACTS project.
0002 //
0003 // Copyright (C) 2016 CERN for the benefit of the ACTS project
0004 //
0005 // This Source Code Form is subject to the terms of the Mozilla Public
0006 // License, v. 2.0. If a copy of the MPL was not distributed with this
0007 // file, You can obtain one at https://mozilla.org/MPL/2.0/.
0008 
0009 #include <boost/test/unit_test.hpp>
0010 
0011 #include "Acts/Definitions/Algebra.hpp"
0012 #include "Acts/Geometry/CuboidVolumeBounds.hpp"
0013 #include "Acts/Geometry/CutoutCylinderVolumeBounds.hpp"
0014 #include "Acts/Geometry/CylinderVolumeBounds.hpp"
0015 #include "Acts/Utilities/BinAdjustmentVolume.hpp"
0016 #include "Acts/Utilities/BinUtility.hpp"
0017 #include "Acts/Utilities/BinningType.hpp"
0018 
0019 #include <cmath>
0020 #include <memory>
0021 #include <numbers>
0022 #include <vector>
0023 
0024 using namespace Acts;
0025 
0026 namespace ActsTests {
0027 
0028 BOOST_AUTO_TEST_SUITE(UtilitiesSuite)
0029 
0030 // Test Cylinder
0031 BOOST_AUTO_TEST_CASE(BinAdjustmentVolume_Cylinder) {
0032   CylinderVolumeBounds bound(10, 50, 150, std::numbers::pi / 2., 0);
0033   BinUtility bu;
0034   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisR);
0035   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisPhi);
0036   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisZ);
0037 
0038   BinUtility buAdjust = adjustBinUtility(bu, bound, Transform3::Identity());
0039 
0040   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].min, 10);
0041   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].max, 50);
0042   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].min,
0043                     -static_cast<float>(std::numbers::pi / 2.));
0044   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].max,
0045                     static_cast<float>(std::numbers::pi / 2.));
0046   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].min, -150);
0047   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].max, 150);
0048 }
0049 
0050 // Test Cutout Cylinder
0051 BOOST_AUTO_TEST_CASE(BinAdjustmentVolume_CutoutCylinder) {
0052   CutoutCylinderVolumeBounds bound(10, 20, 50, 100, 15);
0053   BinUtility bu;
0054   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisR);
0055   bu += BinUtility(1, 0, 1, closed, AxisDirection::AxisPhi);
0056   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisZ);
0057 
0058   BinUtility buAdjust = adjustBinUtility(bu, bound, Transform3::Identity());
0059 
0060   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].min, 10);
0061   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].max, 50);
0062   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].min, -std::numbers::pi_v<float>);
0063   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].max, std::numbers::pi_v<float>);
0064   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].min, -100);
0065   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].max, 100);
0066 }
0067 
0068 // Test Cuboid
0069 BOOST_AUTO_TEST_CASE(BinAdjustmentVolume_Cuboid) {
0070   CuboidVolumeBounds bound(13, 23, 42);
0071   BinUtility bu;
0072   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisX);
0073   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisY);
0074   bu += BinUtility(1, 0, 1, open, AxisDirection::AxisZ);
0075 
0076   BinUtility buAdjust = adjustBinUtility(bu, bound, Transform3::Identity());
0077 
0078   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].min, -13);
0079   BOOST_CHECK_EQUAL(buAdjust.binningData()[0].max, 13);
0080   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].min, -23);
0081   BOOST_CHECK_EQUAL(buAdjust.binningData()[1].max, 23);
0082   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].min, -42);
0083   BOOST_CHECK_EQUAL(buAdjust.binningData()[2].max, 42);
0084 }
0085 
0086 BOOST_AUTO_TEST_SUITE_END()
0087 
0088 }  // namespace ActsTests