IntegratePeaksCWSD.h

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 +
#pragma once
 
#include "MantidAPI/Algorithm.h"
#include "MantidAPI/CompositeFunction.h"
#include "MantidAPI/DeprecatedAlgorithm.h"
#include "MantidAPI/IMDEventWorkspace_fwd.h"
#include "MantidDataObjects/MDEventWorkspace.h"
#include "MantidDataObjects/MaskWorkspace.h"
#include "MantidDataObjects/PeaksWorkspace.h"
#include "MantidMDAlgorithms/DllConfig.h"
 
namespace Mantid {
namespace MDAlgorithms {
 
class MANTID_MDALGORITHMS_DLL IntegratePeaksCWSD final : public API::Algorithm, public API::DeprecatedAlgorithm {
public:
  IntegratePeaksCWSD();
 
  const std::string name() const override { return "IntegratePeaksCWSD"; }
 
  const std::string summary() const override {
    return "Integrate single-crystal peaks in reciprocal space, for "
           "MDEventWorkspaces from reactor-source single crystal "
           "diffractometer.";
  }
 
  int version() const override { return 1; }
  const std::vector<std::string> seeAlso() const override {
    return {"IntegratePeaksHybrid", "IntegratePeaksMDHKL", "IntegratePeaksMD", "IntegratePeaksUsingClusters"};
  }
 
  const std::string category() const override { return "MDAlgorithms\\Peaks;Crystal\\Integration"; }
 
private:
  void init() override;
  void exec() override;
 
  void processInputs();
 
  void simplePeakIntegration(const std::vector<detid_t> &vecMaskedDetID,
                             const std::map<int, signal_t> &run_monitor_map);
  template <typename MDE, size_t nd>
  void integrate(typename DataObjects::MDEventWorkspace<MDE, nd>::sptr ws,
                 const std::map<uint16_t, signal_t> &run_monitor_map);
 
  std::map<int, signal_t> getMonitorCounts();
 
  std::map<int, double> getMeasureTime();
 
  std::vector<detid_t> processMaskWorkspace(const DataObjects::MaskWorkspace_const_sptr &maskws);
 
  void getPeakInformation();
 
  DataObjects::PeaksWorkspace_sptr createOutputs();
 
  void mergePeaks();
 
  void normalizePeaksIntensities();
 
  DataObjects::PeaksWorkspace_sptr createPeakworkspace(Kernel::V3D peakCenter, const API::IMDEventWorkspace_sptr &mdws);
 
  Mantid::API::IMDEventWorkspace_sptr m_inputWS;
 
  Mantid::DataObjects::PeaksWorkspace_sptr m_peaksWS;
 
  bool m_haveMultipleRun;
  std::map<int, double> m_runNormMap;
 
  std::map<int, Kernel::V3D> m_runPeakCenterMap;
  bool m_useSinglePeakCenterFmUser;
  Kernel::V3D m_peakCenter;
  double m_peakRadius;
  bool m_doMergePeak;
  bool m_normalizeByMonitor;
  bool m_normalizeByTime; // NormalizeByTime
  double m_scaleFactor;   // ScaleFactor
 
  std::vector<DataObjects::Peak> m_vecPeaks;
  std::map<int, double> m_runPeakCountsMap;
 
  DataObjects::MaskWorkspace_sptr m_maskWS;
  std::vector<detid_t> vecMaskedDetID;
 
  bool m_haveInputPeakWS;
};
 
} // namespace MDAlgorithms
} // namespace Mantid