BFieldXYZUtils.cxx

#include "Tracking/Sim/BFieldXYZUtils.h"
 
Acts::Vector3 default_transformPos(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 default_transformBField(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 localToGlobalBin_xyz(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;
}