Core/EventData/MultiTrajectoryTests.cpp

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/data/test_case.hpp>
0010 #include <boost/test/tools/context.hpp>
0011 #include <boost/test/unit_test.hpp>
0012
0013 #include "Acts/Definitions/TrackParametrization.hpp"
0014 #include "Acts/EventData/MultiTrajectory.hpp"
0015 #include "Acts/EventData/ProxyAccessor.hpp"
0016 #include "Acts/EventData/TrackStatePropMask.hpp"
0017 #include "Acts/EventData/VectorMultiTrajectory.hpp"
0018 #include "Acts/EventData/detail/MultiTrajectoryTestsCommon.hpp"
0019 #include "Acts/EventData/detail/TestTrackState.hpp"
0020 #include "Acts/Geometry/GeometryContext.hpp"
0021
0022 #include <random>
0023 #include <string>
0024 #include <type_traits>
0025 #include <utility>
0026
0027 namespace {
0028
0029 using namespace Acts;
0030 using namespace Acts::UnitLiterals;
0031 using namespace Acts::detail::Test;
0032 using namespace Acts::HashedStringLiteral;
0033 namespace bd = boost::unit_test::data;
0034
0035 using ParametersVector = BoundVector;
0036 using CovarianceMatrix = BoundMatrix;
0037 using Jacobian = BoundMatrix;
0038
0039 const auto gctx = GeometryContext::dangerouslyDefaultConstruct();
0040 // fixed seed for reproducible tests
0041 std::default_random_engine rng(31415);
0042
0043 struct Factory {
0044   using trajectory_t = VectorMultiTrajectory;
0045   using const_trajectory_t = ConstVectorMultiTrajectory;
0046
0047   VectorMultiTrajectory create() { return {}; }
0048   ConstVectorMultiTrajectory createConst() { return {}; }
0049 };
0050
0051 using CommonTests = MultiTrajectoryTestsCommon<Factory>;
0052
0053 }  // namespace
0054
0055 namespace ActsTests {
0056
0057 BOOST_AUTO_TEST_SUITE(EventDataSuite)
0058
0059 BOOST_AUTO_TEST_CASE(Build) {
0060   CommonTests ct;
0061   ct.testBuild();
0062 }
0063
0064 BOOST_AUTO_TEST_CASE(ConstCorrectness) {
0065   // make mutable
0066   VectorMultiTrajectory t;
0067   auto i0 = t.addTrackState();
0068
0069   BOOST_CHECK(!IsReadOnlyMultiTrajectory<decltype(t)>::value);
0070
0071   {
0072     VectorMultiTrajectory::TrackStateProxy tsp = t.getTrackState(i0);
0073     static_cast<void>(tsp);
0074     VectorMultiTrajectory::ConstTrackStateProxy ctsp{t.getTrackState(i0)};
0075     static_cast<void>(ctsp);
0076
0077     tsp.predicted().setRandom();
0078     // const auto& tsp_const = tsp;
0079     // tsp_const.predicted().setRandom();
0080     // ctsp.predicted().setRandom();
0081   }
0082
0083   // is this something we actually want?
0084   ConstVectorMultiTrajectory ct{t};
0085   BOOST_CHECK_EQUAL(ct.size(), t.size());
0086
0087   ConstVectorMultiTrajectory ctm{std::move(t)};
0088   BOOST_CHECK_EQUAL(ctm.size(), ct.size());
0089
0090   {
0091     static_assert(
0092         std::is_same_v<ConstVectorMultiTrajectory::ConstTrackStateProxy,
0093                        decltype(ctm.getTrackState(i0))>,
0094         "Got mutable track state proxy");
0095     ConstVectorMultiTrajectory::ConstTrackStateProxy ctsp =
0096         ctm.getTrackState(i0);
0097     static_cast<void>(ctsp);
0098
0099     // doesn't compile:
0100     // ctsp.predictedCovariance().setIdentity();
0101   }
0102
0103   // doesn't compile:
0104   // ct.clear();
0105   // ct.addTrackState();
0106 }
0107
0108 BOOST_AUTO_TEST_CASE(Clear) {
0109   CommonTests ct;
0110   ct.testClear();
0111 }
0112
0113 BOOST_AUTO_TEST_CASE(ApplyWithAbort) {
0114   CommonTests ct;
0115   ct.testApplyWithAbort();
0116 }
0117
0118 BOOST_AUTO_TEST_CASE(AddTrackStateWithBitMask) {
0119   CommonTests ct;
0120   ct.testAddTrackStateWithBitMask();
0121 }
0122
0123 BOOST_AUTO_TEST_CASE(AddTrackStateComponents) {
0124   CommonTests ct;
0125   ct.testAddTrackStateComponents();
0126 }
0127
0128 // assert expected "cross-talk" between trackstate proxies
0129 BOOST_AUTO_TEST_CASE(TrackStateProxyCrossTalk) {
0130   CommonTests ct;
0131   ct.testTrackStateProxyCrossTalk(rng);
0132 }
0133
0134 BOOST_AUTO_TEST_CASE(TrackStateReassignment) {
0135   CommonTests ct;
0136   ct.testTrackStateReassignment(rng);
0137 }
0138
0139 BOOST_DATA_TEST_CASE(TrackStateProxyStorage, bd::make({1u, 2u}),
0140                      nMeasurements) {
0141   CommonTests ct;
0142   ct.testTrackStateProxyStorage(rng, nMeasurements);
0143 }
0144
0145 BOOST_AUTO_TEST_CASE(TrackStateProxyAllocations) {
0146   CommonTests ct;
0147   ct.testTrackStateProxyAllocations(rng);
0148 }
0149
0150 BOOST_AUTO_TEST_CASE(TrackStateProxyGetMask) {
0151   CommonTests ct;
0152   ct.testTrackStateProxyGetMask();
0153 }
0154
0155 BOOST_AUTO_TEST_CASE(TrackStateProxyCopy) {
0156   CommonTests ct;
0157   ct.testTrackStateProxyCopy(rng);
0158 }
0159
0160 BOOST_AUTO_TEST_CASE(TrackStateCopyDynamicColumns) {
0161   CommonTests ct;
0162   ct.testTrackStateCopyDynamicColumns();
0163 }
0164
0165 BOOST_AUTO_TEST_CASE(TrackStateProxyCopyDiffMTJ) {
0166   CommonTests ct;
0167   ct.testTrackStateProxyCopyDiffMTJ();
0168 }
0169
0170 BOOST_AUTO_TEST_CASE(ProxyAssignment) {
0171   CommonTests ct;
0172   ct.testProxyAssignment();
0173 }
0174
0175 BOOST_AUTO_TEST_CASE(CopyFromConst) {
0176   CommonTests ct;
0177   ct.testCopyFromConst();
0178 }
0179
0180 BOOST_AUTO_TEST_CASE(TrackStateProxyShare) {
0181   CommonTests ct;
0182   ct.testTrackStateProxyShare(rng);
0183 }
0184
0185 BOOST_AUTO_TEST_CASE(MultiTrajectoryExtraColumns) {
0186   CommonTests ct;
0187   ct.testMultiTrajectoryExtraColumns();
0188 }
0189
0190 BOOST_AUTO_TEST_CASE(MultiTrajectoryExtraColumnsRuntime) {
0191   CommonTests ct;
0192   ct.testMultiTrajectoryExtraColumnsRuntime();
0193 }
0194
0195 BOOST_AUTO_TEST_CASE(MultiTrajectoryAllocateCalibratedInit) {
0196   CommonTests ct;
0197   ct.testMultiTrajectoryAllocateCalibratedInit(rng);
0198 }
0199
0200 BOOST_AUTO_TEST_CASE(MemoryStats) {
0201   using namespace boost::histogram;
0202   using cat = axis::category<std::string>;
0203
0204   VectorMultiTrajectory mt;
0205
0206   auto stats = mt.statistics();
0207
0208   std::stringstream ss;
0209   stats.toStream(ss);
0210   std::string out = ss.str();
0211   BOOST_CHECK(!out.empty());
0212   BOOST_CHECK_NE(out.find("total"), std::string::npos);
0213
0214   const auto& h = stats.hist;
0215
0216   auto column_axis = axis::get<cat>(h.axis(0));
0217   auto type_axis = axis::get<axis::category<>>(h.axis(1));
0218
0219   for (int t = 0; t < type_axis.size(); t++) {
0220     for (int c = 0; c < column_axis.size(); c++) {
0221       double v = h.at(c, t);
0222       BOOST_CHECK_EQUAL(v, 0.0);
0223     }
0224   }
0225
0226   TestTrackState pc(rng, 2u);
0227   auto ts = mt.makeTrackState();
0228   fillTrackState<VectorMultiTrajectory>(pc, TrackStatePropMask::All, ts);
0229
0230   stats = mt.statistics();
0231
0232   for (int t = 0; t < type_axis.size(); t++) {
0233     BOOST_TEST_CONTEXT((type_axis.bin(t) == 1 ? "meas" : "other"))
0234     for (int c = 0; c < column_axis.size(); c++) {
0235       std::string key = column_axis.bin(c);
0236       BOOST_TEST_CONTEXT("column: " << key) {
0237         double v = h.at(c, t);
0238         if (t == 0) {
0239           BOOST_CHECK_NE(v, 0.0);
0240         } else {
0241           BOOST_CHECK_EQUAL(v, 0.0);
0242         }
0243       }
0244     }
0245   }
0246 }
0247
0248 BOOST_AUTO_TEST_CASE(Accessors) {
0249   VectorMultiTrajectory mtj;
0250   mtj.addColumn<unsigned int>("ndof");
0251   mtj.addColumn<double>("super_chi2");
0252
0253   auto ts = mtj.makeTrackState();
0254
0255   ProxyAccessor<unsigned int> ndof("ndof");
0256   ConstProxyAccessor<unsigned int> ndofConst("ndof");
0257   ProxyAccessor<double> superChi2("super_chi2");
0258   ConstProxyAccessor<double> superChi2Const("super_chi2");
0259
0260   ndof(ts) = 65;
0261   BOOST_CHECK_EQUAL((ts.component<unsigned int, "ndof"_hash>()), 65);
0262   BOOST_CHECK_EQUAL(ndofConst(ts), 65);
0263
0264   // should not compile
0265   // ndofConst(ts) = 66;
0266
0267   superChi2(ts) = 123.45;
0268   BOOST_CHECK_EQUAL((ts.component<double, "super_chi2"_hash>()), 123.45);
0269   BOOST_CHECK_EQUAL(superChi2Const(ts), 123.45);
0270
0271   // should not compile
0272   // superChi2Const(ts) = 66.66;
0273 }
0274
0275 BOOST_AUTO_TEST_CASE(ChangeSourceLinkType) {
0276   VectorMultiTrajectory mtj;
0277   auto ts = mtj.makeTrackState();
0278
0279   int value = 5;
0280   ts.setUncalibratedSourceLink(SourceLink{value});
0281
0282   BOOST_CHECK_EQUAL(ts.getUncalibratedSourceLink().get<int>(), value);
0283   BOOST_CHECK_THROW(ts.getUncalibratedSourceLink().get<double>(),
0284                     std::bad_any_cast);
0285
0286   double otherValue = 42.42;
0287
0288   // this changes the stored type
0289   ts.setUncalibratedSourceLink(SourceLink{otherValue});
0290   BOOST_CHECK_EQUAL(ts.getUncalibratedSourceLink().get<double>(), otherValue);
0291   BOOST_CHECK_THROW(ts.getUncalibratedSourceLink().get<int>(),
0292                     std::bad_any_cast);
0293 }
0294
0295 BOOST_AUTO_TEST_SUITE_END()
0296
0297 }  // namespace ActsTests