FindDeadDetectors.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 "MantidAlgorithms/FindDeadDetectors.h"
#include "MantidAPI/MatrixWorkspace.h"
#include "MantidKernel/BoundedValidator.h"
 
#include <fstream>
 
namespace Mantid::Algorithms {
 
// Register the class into the algorithm factory
DECLARE_ALGORITHM(FindDeadDetectors)
 
using namespace Kernel;
using namespace API;
 
void FindDeadDetectors::init() {
  declareProperty(std::make_unique<WorkspaceProperty<>>("InputWorkspace", "", Direction::Input),
                  "Name of the input workspace");
  declareProperty(std::make_unique<WorkspaceProperty<>>("OutputWorkspace", "", Direction::Output),
                  "Each histogram from the input workspace maps to a histogram in this\n"
                  "workspace with one value that indicates if there was a dead detector");
 
  auto mustBePositive = std::make_shared<BoundedValidator<double>>();
  mustBePositive->setLower(0);
  declareProperty("DeadThreshold", 0.0, mustBePositive,
                  "The threshold against which to judge if a spectrum belongs to a dead\n"
                  "detector");
  // As the property takes ownership of the validator pointer, have to take care
  // to pass in a unique
  // pointer to each property.
  declareProperty("LiveValue", 0.0, mustBePositive,
                  "The value to assign to an integrated spectrum flagged as 'live'\n"
                  "(default 0.0)");
  declareProperty("DeadValue", 100.0, mustBePositive,
                  "The value to assign to an integrated spectrum flagged as 'dead'\n"
                  "(default 100.0)");
  // EMPTY_DBL() is a tag that tells us that no value has been set and we want
  // to use the default
  declareProperty("RangeLower", EMPTY_DBL(),
                  "No bin with a boundary at an x value less than this will be used\n"
                  "in the summation that decides if a detector is 'dead' (default: the\n"
                  "start of each histogram)");
  declareProperty("RangeUpper", EMPTY_DBL(),
                  "No bin with a boundary at an x value higher than this value will\n"
                  "be used in the summation that decides if a detector is 'dead'\n"
                  "(default: the end of each histogram)");
  declareProperty("OutputFile", "", "A filename to which to write the list of dead detector UDETs");
  // This output property will contain the list of UDETs for the dead detectors
  declareProperty("FoundDead", std::vector<detid_t>(), Direction::Output);
}
 
void FindDeadDetectors::exec() {
  double deadThreshold = getProperty("DeadThreshold");
  double liveValue = getProperty("LiveValue");
  double deadValue = getProperty("DeadValue");
 
  // Try and open the output file, if specified, and write a header
  std::ofstream file(getPropertyValue("OutputFile").c_str());
  file << "Index Spectrum UDET(S)\n";
 
  // Get the integrated input workspace
  MatrixWorkspace_sptr integratedWorkspace = integrateWorkspace();
 
  std::vector<detid_t> deadDets;
  int countSpec = 0, countDets = 0;
 
  // iterate over the data values setting the live and dead values
  g_log.information() << "Marking dead detectors\n";
  const int64_t numSpec = integratedWorkspace->getNumberHistograms();
  const auto numSpec_d = static_cast<double>(numSpec);
  int64_t iprogress_step = numSpec / 100;
  if (iprogress_step == 0)
    iprogress_step = 1;
  for (int64_t i = 0; i < int64_t(numSpec); ++i) {
    // Spectrum in the integratedWorkspace
    double &y = integratedWorkspace->mutableY(i)[0];
    if (y > deadThreshold) {
      y = liveValue;
    } else {
      ++countSpec;
      y = deadValue;
      const auto &spec = integratedWorkspace->getSpectrum(i);
      // Write the spectrum number to file
      file << i << " " << spec.getSpectrumNo();
      for (const auto &id : spec.getDetectorIDs()) {
        // Write the detector ID to file, log & the FoundDead output property
        file << " " << id;
        // we could write dead detectors to the log but if they are viewing the
        // log in the MantidPlot viewer it will crash MantidPlot
        deadDets.emplace_back(id);
        ++countDets;
      }
      file << '\n';
    }
    if (i % iprogress_step == 0) {
      progress(static_cast<double>(i) / numSpec_d);
      interruption_point();
    }
  }
 
  g_log.notice() << "Found a total of " << countDets << " 'dead' detectors within " << countSpec
                 << " 'dead' spectra.\n";
 
  // Assign it to the output workspace property
  setProperty("OutputWorkspace", integratedWorkspace);
  setProperty("FoundDead", deadDets);
 
  // Close the output file
  file.close();
}
 
MatrixWorkspace_sptr FindDeadDetectors::integrateWorkspace() {
  g_log.information() << "Integrating input workspace\n";
 
  auto childAlg = createChildAlgorithm("Integration");
  // Now execute integration.
  // pass inputed values straight to Integration, checking must be done there
  childAlg->setProperty<MatrixWorkspace_sptr>("InputWorkspace", getProperty("InputWorkspace"));
  childAlg->setProperty<MatrixWorkspace_sptr>("OutputWorkspace", getProperty("OutputWorkspace"));
  childAlg->setProperty<double>("RangeLower", getProperty("RangeLower"));
  childAlg->setProperty<double>("RangeUpper", getProperty("RangeUpper"));
  childAlg->executeAsChildAlg();
 
  return childAlg->getProperty("OutputWorkspace");
}
 
} // namespace Mantid::Algorithms