ConstantBFieldTests.cpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file ConstantBFieldTests.cpp

// This file is part of the Acts project.
//
// Copyright (C) 2017-2018 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.


#include <boost/test/data/test_case.hpp>
#include <boost/test/unit_test.hpp>

#include "Acts/Definitions/Algebra.hpp"
#include "Acts/Definitions/Units.hpp"
#include "Acts/MagneticField/ConstantBField.hpp"
#include "Acts/MagneticField/MagneticFieldContext.hpp"
#include "Acts/Utilities/Result.hpp"

#include <utility>

namespace bdata = boost::unit_test::data;
namespace tt = boost::test_tools;
using namespace Acts::UnitLiterals;

namespace Acts {
namespace Test {

// Create a test context
MagneticFieldContext mfContext = MagneticFieldContext();

BOOST_DATA_TEST_CASE(ConstantBField_components,
                     bdata::random(-2_T, 2_T) ^ bdata::random(-1_T, 4_T) ^
                         bdata::random(0_T, 10_T) ^ bdata::random(-10_m, 10_m) ^
                         bdata::random(-10_m, 10_m) ^
                         bdata::random(-10_m, 10_m) ^ bdata::xrange(10),
                     x, y, z, bx, by, bz, index) {
  (void)index;
  const Vector3 Btrue(bx, by, bz);
  const Vector3 pos(x, y, z);
  const ConstantBField BField(Btrue);

  auto bCache = BField.makeCache(mfContext);

  BOOST_CHECK_EQUAL(Btrue, BField.getField());

  BOOST_CHECK_EQUAL(Btrue, BField.getField(pos, bCache).value());
  BOOST_CHECK_EQUAL(Btrue, BField.getField(Vector3(0, 0, 0), bCache).value());
  BOOST_CHECK_EQUAL(Btrue, BField.getField(-2 * pos, bCache).value());
}

BOOST_DATA_TEST_CASE(ConstantBField_update,
                     bdata::random(-2_T, 2_T) ^ bdata::random(-1_T, 4_T) ^
                         bdata::random(0_T, 10_T) ^ bdata::random(-10_m, 10_m) ^
                         bdata::random(-10_m, 10_m) ^
                         bdata::random(-10_m, 10_m) ^ bdata::xrange(10),
                     x, y, z, bx, by, bz, index) {
  (void)index;

  ConstantBField BField{Vector3{0, 0, 0}};
  const Vector3 Btrue(bx, by, bz);
  const Vector3 pos(x, y, z);
  BField.setField(Vector3{bx, by, bz});

  auto bCache = BField.makeCache(mfContext);

  BOOST_CHECK_EQUAL(Btrue, BField.getField());

  BOOST_CHECK_EQUAL(Btrue, BField.getField(pos, bCache).value());
  BOOST_CHECK_EQUAL(Btrue, BField.getField(Vector3(0, 0, 0), bCache).value());
  BOOST_CHECK_EQUAL(Btrue, BField.getField(-2 * pos, bCache).value());
}

}  // namespace Test
}  // namespace Acts