HKLFilterWavelength.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/Crystal/HKLFilterWavelength.h"
#include <sstream>
#include <stdexcept>
#include <utility>
 
namespace Mantid::Geometry {
 
using namespace Kernel;
 
HKLFilterWavelength::HKLFilterWavelength(Kernel::DblMatrix ub, double lambdaMin, double lambdaMax)
    : m_ub(std::move(ub)), m_lambdaMin(lambdaMin), m_lambdaMax(lambdaMax) {
  checkProperLambdaRangeValues();
}
 
std::string HKLFilterWavelength::getDescription() const noexcept {
  std::ostringstream strm;
  strm << "(" << m_lambdaMin << " <= lambda <= " << m_lambdaMax << ")";
 
  return strm.str();
}
 
bool HKLFilterWavelength::isAllowed(const Kernel::V3D &hkl) const noexcept {
  V3D q = m_ub * hkl;
  double lambda = (2.0 * q.Z()) / (q.norm2());
 
  return lambda >= m_lambdaMin && lambda <= m_lambdaMax;
}
 
void HKLFilterWavelength::checkProperLambdaRangeValues() const {
  if (m_lambdaMin <= 0.0) {
    throw std::range_error("LambdaMin cannot be <= 0.");
  }
 
  if (m_lambdaMax <= m_lambdaMin) {
    throw std::range_error("LambdaMax cannot be smaller than LambdaMin.");
  }
}
 
} // namespace Mantid::Geometry