LoadANSTOHelper.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
 
//---------------------------------------------------
// Includes
//---------------------------------------------------
 
#include "MantidAPI/IFileLoader.h"
#include "MantidAPI/LogManager.h"
#include "MantidDataObjects/EventWorkspace.h"
#include "MantidGeometry/Instrument.h"
#include "MantidNexus/NexusClasses_fwd.h"
#include <algorithm>
#include <regex>
 
#define TarTypeFlag_NormalFile '0'
#define TarTypeFlag_HardLink '1'
#define TarTypeFlag_SymbolicLink '2'
#define TarTypeFlag_CharacterSpecial '3'
#define TarTypeFlag_BlockSpecial '4'
#define TarTypeFlag_Directory '5'
#define TarTypeFlag_FIFO '6'
#define TarTypeFlag_ContiguousFile '7'
 
namespace Mantid {
namespace DataHandling {
namespace ANSTO {
 
std::vector<std::string> filterDatasets(const Nexus::NXEntry &entry, const std::string &groupAddress,
                                        const std::string &regexFilter);
 
using EventVector_pt = std::vector<Types::Event::TofEvent> *;
 
class ProgressTracker {
private:
  // fields
  const std::string m_msg;
  size_t m_count;
  int64_t m_step;
  int64_t m_next;
  // matntid
  API::Progress &m_progBar;
 
public:
  // construction
  ProgressTracker(API::Progress &progBar, const char *msg, int64_t target, size_t count);
  ~ProgressTracker();
 
  // methods
  void update(int64_t position);
  void complete();
  void setTarget(int64_t target);
};
 
class EventProcessor {
protected:
  // fields
  const std::vector<bool> &m_roi;
  const size_t m_stride;
  // number of frames
  size_t m_frames;
  size_t m_framesValid;
  int64_t m_startTime;
  // tof correction
  const double m_period;
  const double m_phase;
  // boundaries
  const double m_tofMinBoundary;
  const double m_tofMaxBoundary;
  const double m_timeMinBoundary;
  const double m_timeMaxBoundary;
 
  // methods
  bool validFrame() const;
  virtual void addEventImpl(size_t id, int64_t pulse, double tof) = 0;
 
public:
  // construction
  EventProcessor(const std::vector<bool> &roi, size_t stride, const double period, const double phase,
                 const int64_t startTime, const double tofMinBoundary, const double tofMaxBoundary,
                 const double timeMinBoundary, const double timeMaxBoundary);
 
  // methods
  void newFrame();
  void addEvent(size_t x, size_t y, double tof);
};
 
class EventCounter : public EventProcessor {
protected:
  // fields
  std::vector<size_t> &m_eventCounts;
  // tof
  double m_tofMin;
  double m_tofMax;
 
  // methods
  void addEventImpl(size_t id, int64_t pulse, double tof) override;
 
public:
  // construction
  EventCounter(const std::vector<bool> &roi, const size_t stride, const double period, const double phase,
               const int64_t startTime, const double tofMinBoundary, const double tofMaxBoundary,
               const double timeMinBoundary, const double timeMaxBoundary, std::vector<size_t> &eventCounts);
 
  // properties
  size_t numFrames() const;
  double tofMin() const;
  double tofMax() const;
};
 
class EventAssigner : public EventProcessor {
protected:
  // fields
  std::vector<EventVector_pt> &m_eventVectors;
 
  // methods
  void addEventImpl(size_t id, int64_t pulse, double tof) override;
 
public:
  // construction
  EventAssigner(const std::vector<bool> &roi, const size_t stride, const double period, const double phase,
                int64_t startTime, const double tofMinBoundary, const double tofMaxBoundary,
                const double timeMinBoundary, const double timeMaxBoundary, std::vector<EventVector_pt> &eventVectors);
};
 
class EventAssignerFixedWavelength : public EventAssigner {
protected:
  // fields
  double m_wavelength;
 
  // methods
  void addEventImpl(size_t id, int64_t pulse, double tof) override;
 
public:
  // construction
  EventAssignerFixedWavelength(const std::vector<bool> &roi, const size_t stride, const double wavelength,
                               const double period, const double phase, const int64_t startTime,
                               const double tofMinBoundary, const double tofMaxBoundary, const double timeMinBoundary,
                               const double timeMaxBoundary, std::vector<EventVector_pt> &eventVectors);
};
 
class FastReadOnlyFile {
private:
#ifdef _WIN32
  HANDLE m_handle;
#else
  FILE *m_handle;
#endif
public:
  // construction
  FastReadOnlyFile(const char *filename);
  ~FastReadOnlyFile();
 
  // Prevent copying of a file handle
  FastReadOnlyFile(const FastReadOnlyFile &) = delete;
  FastReadOnlyFile &operator=(FastReadOnlyFile) = delete;
 
  // properties
  void *handle() const;
 
  // methods
  void close();
  bool read(void *buffer, uint32_t size);
  bool seek(int64_t offset, int whence, int64_t *newPosition = nullptr);
};
 
namespace Tar {
 
struct EntryHeader {
  char FileName[100];
  char FileMode[8];
  char OwnerUserID[8];
  char OwnerGroupID[8];
  char FileSize[12];         // in bytes (octal base)
  char LastModification[12]; // time in numeric Unix time format (octal)
  char Checksum[8];
  char TypeFlag;
  char LinkedFileName[100];
  char UStar[8];
  char OwnerUserName[32];
  char OwnerGroupName[32];
  char DeviceMajorNumber[8];
  char DeviceMinorNumber[8];
  char FilenamePrefix[155];
 
  // methods
  void writeChecksum();
  void writeFileSize(int64_t value);
  int64_t readFileSize();
};
 
class File {
 
  static const auto BUFFER_SIZE = 4096;
 
  struct FileInfo {
    int64_t Offset;
    int64_t Size;
  };
 
private:
  // fields
  bool m_good;
  FastReadOnlyFile m_file;
  std::vector<std::string> m_fileNames;
  std::vector<FileInfo> m_fileInfos;
  // selected file
  size_t m_selected; // index
  int64_t m_position;
  int64_t m_size;
  // buffer
  uint8_t m_buffer[BUFFER_SIZE];
  size_t m_bufferPosition;
  size_t m_bufferAvailable;
 
  // not supported
  File(const File &) = delete;
  File &operator=(const File &) = delete;
 
public:
  // construction
  explicit File(const std::string &path);
  void close();
 
  // properties
  bool good() const;
  const std::vector<std::string> &files() const;
  // from selected file
  const std::string &selected_name() const;
  int64_t selected_position() const;
  int64_t selected_size() const;
 
  // methods
  bool select(const char *file);
  bool skip(uint64_t offset);
  size_t read(void *dst, size_t size);
  int read_byte();
 
  // helpers
  static bool append(const std::string &path, const std::string &name, const void *buffer, size_t size);
};
 
} // namespace Tar
 
namespace Anxs {
// options to capture timeseries data
enum class ScanLog { Start, End, Mean };
 
std::string extractWorkspaceTitle(const std::string &nxsFile);
 
int64_t epochRelDateTimeBase(int64_t epochInNanoSeconds);
 
template <typename T> bool loadNXDataSet(const Nexus::NXEntry &entry, const std::string &path, T &value, int index);
bool loadNXString(const Nexus::NXEntry &entry, const std::string &path, std::string &value);
 
bool isTimedDataSet(const Nexus::NXEntry &entry, const std::string &path);
std::pair<uint64_t, uint64_t> getTimeScanLimits(const Nexus::NXEntry &entry, int datasetIx);
std::pair<uint64_t, uint64_t> getHMScanLimits(const Nexus::NXEntry &entry, int datasetIx);
 
template <typename T>
uint64_t extractTimedDataSet(const Nexus::NXEntry &entry, const std::string &path, uint64_t startTime, uint64_t endTime,
                             std::vector<uint64_t> &times, std::vector<T> &events, std::string &units);
template <typename T>
bool extractTimedDataSet(const Nexus::NXEntry &entry, const std::string &path, uint64_t startTime, uint64_t endTime,
                         ScanLog valueOption, uint64_t &eventTime, T &eventValue, std::string &units);
 
void ReadEventData(ProgressTracker &prog, const Nexus::NXEntry &entry, EventProcessor *handler, uint64_t start_nsec,
                   uint64_t end_nsec, const std::string &neutron_path, int tube_resolution = 1024);
 
} // namespace Anxs
 
} // namespace ANSTO
} // namespace DataHandling
} // namespace Mantid