MDNormDirectSC.h

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2014 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 "MantidMDAlgorithms/SlicingAlgorithm.h"
 
namespace Mantid {
namespace DataObjects {
class EventWorkspace;
}
namespace MDAlgorithms {
 
class MANTID_MDALGORITHMS_DLL MDNormDirectSC : public SlicingAlgorithm {
public:
  MDNormDirectSC();
 
  const std::string name() const override;
  int version() const override;
  const std::vector<std::string> seeAlso() const override { return {"MDNormSCD", "MDNormSCDPreprocessIncoherent"}; }
  const std::string category() const override;
  const std::string summary() const override;
  static constexpr double CHARGEBINSIZE = 1.0;    // Proton charge bin size in uA.hr for useLogTimes normalization
  static constexpr double GONIOBINSTEP = 0.25;    // Bin size in angle in degrees for useLogTimes normalization
  static constexpr double MINPROTONCHARGE = 0.1;  // Proton charge below which an angle will be ignored as no data
  static constexpr double STATIONARYANGLIM = 0.1; // Angle size (degrees) below which a gonio is considered stationary
 
private:
  void init() override;
  void exec() override;
 
  void cacheInputs();
  std::string inputEnergyMode() const;
 
  DataObjects::MDHistoWorkspace_sptr binInputWS();
  void createNormalizationWS(const DataObjects::MDHistoWorkspace &dataWS);
  std::vector<coord_t> getValuesFromOtherDimensions(bool &skipNormalization, uint16_t expInfoIndex = 0) const;
  Kernel::Matrix<coord_t> findIntergratedDimensions(const std::vector<coord_t> &otherDimValues,
                                                    bool &skipNormalization);
  void initializeLogTimes();
  void cacheDimensionXValues();
  void calculateNormalization(const std::vector<coord_t> &otherValues, const Kernel::Matrix<coord_t> &affineTrans,
                              uint16_t expInfoIndex);
  void calculateNormContinuous(const std::vector<coord_t> &otherValues, const Kernel::Matrix<coord_t> &affineTrans,
                               uint16_t expInfoIndex);
  void calculateNormInner(const API::SpectrumInfo &spectrumInfo, const double protonCharge,
                          const std::vector<coord_t> &otherValues, const Kernel::Matrix<coord_t> &affineTrans);
 
  void calculateIntersections(std::vector<std::array<double, 4>> &intersections, const double theta, const double phi);
 
  DataObjects::MDHistoWorkspace_sptr m_normWS;
  API::IMDEventWorkspace_sptr m_inputWS;
  coord_t m_hmin, m_hmax, m_kmin, m_kmax, m_lmin, m_lmax, m_dEmin, m_dEmax;
  double m_Ei, m_ki, m_kfmin, m_kfmax;
  bool m_hIntegrated, m_kIntegrated, m_lIntegrated, m_dEIntegrated;
  Mantid::Kernel::DblMatrix m_rubw;
  size_t m_hIdx, m_kIdx, m_lIdx, m_eIdx;
  std::vector<double> m_hX, m_kX, m_lX, m_eX;
  Kernel::V3D m_samplePos;
  Kernel::V3D m_beamDir;
  std::string convention;
  std::vector<std::atomic<signal_t>> m_signalArray;
  uint16_t m_numExptInfos;
  std::unique_ptr<API::Progress> m_progress;
};
 
} // namespace MDAlgorithms
} // namespace Mantid