Detector.cpp

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2018 ISIS Rutherford Appleton Laboratory UKRI,
//   NScD Oak Ridge National Laboratory, European Spallation Source,
//   Institut Laue - Langevin & CSNS, Institute of High Energy Physics, CAS
// SPDX - License - Identifier: GPL - 3.0 +
#include "MantidGeometry/Instrument/Detector.h"
#include "MantidBeamline/DetectorInfo.h"
#include "MantidGeometry/Instrument.h"
#include "MantidGeometry/Instrument/ComponentVisitor.h"
#include "MantidGeometry/Instrument/DetectorInfo.h"
#include "MantidGeometry/Instrument/ParameterMap.h"
#include "MantidKernel/EigenConversionHelpers.h"
#include "MantidKernel/Logger.h"
 
namespace Mantid::Geometry {
namespace {
// static logger object
Kernel::Logger g_log("Detector");
} // namespace
 
using Kernel::V3D;
 
Detector::Detector(const Detector *base, const ParameterMap *map) : ObjComponent(base, map), m_id(base->m_id) {}
 
Detector::Detector(const std::string &name, int id, IComponent *parent)
    : IDetector(), ObjComponent(name, parent), m_id(id) {}
 
Detector::Detector(const std::string &name, int id, const std::shared_ptr<IObject> &shape, IComponent *parent)
    : IDetector(), ObjComponent(name, shape, parent), m_id(id) {}
 
detid_t Detector::getID() const {
  if (m_map) {
    const auto *d = dynamic_cast<const Detector *>(m_base);
    if (d) {
      return d->getID();
    }
  }
 
  return m_id;
}
 
double Detector::getDistance(const IComponent &comp) const { return ObjComponent::getDistance(comp); }
 
double Detector::getTwoTheta(const V3D &observer, const V3D &axis) const {
  const V3D sampleDetVec = this->getPos() - observer;
  return sampleDetVec.angle(axis);
}
 
double Detector::getSignedTwoTheta(const V3D &observer, const V3D &axis, const V3D &instrumentUp) const {
  const V3D sampleDetVec = this->getPos() - observer;
  double angle = sampleDetVec.angle(axis);
 
  V3D cross = axis.cross_prod(sampleDetVec);
  V3D normToSurface = axis.cross_prod(instrumentUp);
  if (normToSurface.scalar_prod(cross) < 0) {
    angle *= -1;
  }
  return angle;
}
 
double Detector::getPhi() const {
  const Kernel::V3D pos = this->getPos();
  return std::atan2(pos[1], pos[0]);
}
 
double Detector::getPhiOffset(const double &offset) const {
  double avgPos = getPhi();
  return avgPos < 0 ? -(offset + avgPos) : offset - avgPos;
}
 
det_topology Detector::getTopology(V3D &center) const {
  center = this->getPos();
  return rect;
}
 
const ParameterMap &Detector::parameterMap() const { return *m_map; }
 
size_t Detector::index() const { return m_map->detectorIndex(m_id); }
 
size_t Detector::registerContents(ComponentVisitor &componentVisitor) const {
  return componentVisitor.registerDetector(*this);
}
 
bool Detector::hasDetectorInfo() const {
  const IComponent *root = m_base;
  while (auto parent = root->getBareParent())
    root = parent;
  auto instrument = dynamic_cast<const Instrument *>(root);
  return m_map->hasDetectorInfo(instrument);
}
 
} // namespace Mantid::Geometry