AlgebraJsonConverter.cpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2021 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 "Acts/Plugins/Json/AlgebraJsonConverter.hpp"

#include <array>
#include <stdexcept>

void Acts::to_json(nlohmann::json& j, const Acts::Transform3& t) {
  auto translation = t.translation();
  if (translation != Acts::Vector3(0., 0., 0)) {
    std::array<Acts::ActsScalar, 3> tdata = {translation.x(), translation.y(),
                                             translation.z()};
    j["translation"] = tdata;
  } else {
    j["translation"] = nlohmann::json();
  }

  auto rotation = t.rotation();
  if (rotation != Acts::RotationMatrix3::Identity()) {
    std::array<Acts::ActsScalar, 9> rdata = {
        rotation(0, 0), rotation(0, 1), rotation(0, 2),
        rotation(1, 0), rotation(1, 1), rotation(1, 2),
        rotation(2, 0), rotation(2, 1), rotation(2, 2)};
    j["rotation"] = rdata;
  } else {
    j["rotation"] = nlohmann::json();
  }
}

void Acts::from_json(const nlohmann::json& j, Acts::Transform3& t) {
  t = Acts::Transform3::Identity();
  if (j.find("rotation") != j.end() and not j["rotation"].empty()) {
    std::array<Acts::ActsScalar, 9> rdata = j["rotation"];
    Acts::RotationMatrix3 rot;
    rot << rdata[0], rdata[1], rdata[2], rdata[3], rdata[4], rdata[5], rdata[6],
        rdata[7], rdata[8];
    t.prerotate(rot);
  }
  if (j.find("translation") != j.end() and not j["translation"].empty()) {
    std::array<Acts::ActsScalar, 3> tdata = j["translation"];
    t.pretranslate(Acts::Vector3(tdata[0], tdata[1], tdata[2]));
  }
}

nlohmann::json Acts::Transform3JsonConverter::toJson(const Transform3& t,
                                                     const Options& options) {
  nlohmann::json jTransform;
  // Write out the translation
  auto translation = t.translation();
  if (translation != Acts::Vector3(0., 0., 0) or options.writeIdentity) {
    std::array<Acts::ActsScalar, 3> tdata = {translation.x(), translation.y(),
                                             translation.z()};
    jTransform["translation"] = tdata;
  } else {
    jTransform["translation"] = nlohmann::json();
  }
  // Write out the rotation, could be transposed
  auto rotation = options.transpose ? t.rotation().transpose() : t.rotation();
  if (rotation != Acts::RotationMatrix3::Identity() or options.writeIdentity) {
    std::array<Acts::ActsScalar, 9> rdata = {
        rotation(0, 0), rotation(0, 1), rotation(0, 2),
        rotation(1, 0), rotation(1, 1), rotation(1, 2),
        rotation(2, 0), rotation(2, 1), rotation(2, 2)};
    jTransform["rotation"] = rdata;
  } else {
    jTransform["rotation"] = nlohmann::json();
  }
  // Return the converted transform
  return jTransform;
}

Acts::Transform3 Acts::Transform3JsonConverter::fromJson(
    const nlohmann::json& jTransform) {
  Transform3 t = Transform3::Identity();
  Acts::from_json(jTransform, t);
  return t;
}