BFieldXYZUtils.cxx

#include "Tracking/Sim/BFieldXYZUtils.h"
 
#include "Acts/Utilities/Result.hpp"
 
Acts::Vector3 defaultTransformPos(const Acts::Vector3& pos_) {
  Acts::Vector3 rot_pos;
  rot_pos(0) = pos_(1);
  rot_pos(1) = pos_(2);
  rot_pos(2) = pos_(0) + DIPOLE_OFFSET;
 
  return rot_pos;
}
 
Acts::Vector3 defaultTransformBField(const Acts::Vector3& field,
                                     const Acts::Vector3& /*pos_*/) {
  Acts::Vector3 rot_field;
 
  rot_field(0) = field(2);
  rot_field(1) = field(0);
  rot_field(2) = field(1);
 
  // std::cout<<"PF::DEFAULT TRANSFORM"<<std::endl;
  // std::cout<<"PF::Check:: transforming"<<std::endl;
  // std::cout<<field<<std::endl;
  // std::cout<<"TO"<<std::endl;
  // std::cout<<rot_field<<std::endl;
 
  return rot_field;
}
 
size_t localToGlobalBinXyz(std::array<size_t, 3> bins,
                           std::array<size_t, 3> sizes) {
  return (bins[0] * (sizes[1] * sizes[2]) + bins[1] * sizes[2] +
          bins[2]);  // xyz - field space
  // return (bins[1] * (sizes[2] * sizes[0]) + bins[2] * sizes[0] + bins[0]);
  // //zxy
}
 
void testField(const std::shared_ptr<Acts::MagneticFieldProvider> bfield,
               const Acts::Vector3& eval_pos,
               const Acts::MagneticFieldContext bctx) {
  Acts::MagneticFieldProvider::Cache cache = bfield->makeCache(bctx);
  std::cout << "Pos::\n" << eval_pos << std::endl;
  std::cout << " BField::\n"
            << bfield->getField(eval_pos, cache).value() /
                   Acts::UnitConstants::T
            << std::endl;
}