MayersSampleCorrectionStrategy.h

Go to the documentation of this file.

// Mantid Repository : https://github.com/mantidproject/mantid
//
// Copyright © 2015 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 "MantidAlgorithms/DllConfig.h"
#include "MantidHistogramData/Histogram.h"
#include <memory>
#include <utility>
#include <vector>
 
namespace Mantid {
namespace Kernel {
class PseudoRandomNumberGenerator;
}
namespace Algorithms {
 
class MANTID_ALGORITHMS_DLL MayersSampleCorrectionStrategy {
public:
  struct Parameters {
    bool mscat;       
    double l1;        
    double l2;        
    double twoTheta;  
    double azimuth;   
    double rho;       
    double sigmaSc;   
    double sigmaAbs;  
    double cylRadius; 
    double cylHeight; 
    size_t msNEvents; 
    size_t msNRuns;   
  };
 
  MayersSampleCorrectionStrategy(MayersSampleCorrectionStrategy::Parameters params, HistogramData::Histogram inputHist);
  ~MayersSampleCorrectionStrategy();
 
  Mantid::HistogramData::Histogram getCorrectedHisto();
 
  double calculateSelfAttenuation(const double muR);
  std::pair<double, double> calculateMS(const size_t irp, const double muR, const double abs);
 
private:
  inline double muRmin() const { return m_muRrange.first; }
  inline double muRmax() const { return m_muRrange.second; }
 
  std::pair<double, double> calculateMuRange() const;
  double muR(const double flightPath, const double tof) const;
  double muR(const double sigt) const;
  double sigmaTotal(const double flightPath, const double tof) const;
  void seedRNG(const size_t seed);
 
  const Parameters m_pars;
  // Holds histogram to process
  const HistogramData::Histogram m_histogram;
  const HistogramData::Points m_tofVals;
  const size_t m_histoYSize;
  const std::pair<double, double> m_muRrange;
  std::unique_ptr<Kernel::PseudoRandomNumberGenerator> m_rng;
};
 
} // namespace Algorithms
} // namespace Mantid