LoadSpiceAscii.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 <boost/algorithm/string.hpp>
#include <fstream>
 
#include "MantidAPI/FileLoaderRegistry.h"
#include "MantidAPI/FileProperty.h"
#include "MantidAPI/MatrixWorkspace.h"
#include "MantidAPI/Run.h"
#include "MantidAPI/TableRow.h"
#include "MantidAPI/WorkspaceFactory.h"
#include "MantidAPI/WorkspaceProperty.h"
#include "MantidDataHandling/LoadSpiceAscii.h"
#include "MantidDataObjects/TableWorkspace.h"
#include "MantidKernel/ArrayProperty.h"
 
#include <boost/algorithm/string/finder.hpp>
#include <boost/algorithm/string/iter_find.hpp>
 
using namespace boost::algorithm;
 
using namespace Mantid::API;
using namespace Mantid::Kernel;
using namespace Mantid::DataHandling;
using Mantid::Types::Core::DateAndTime;
 
namespace Mantid::DataHandling {
 
// DECLARE_FILELOADER_ALGORITHM(LoadSpiceAscii)
DECLARE_ALGORITHM(LoadSpiceAscii)
 
static bool endswith(const std::string &s, const std::string &subs) {
  // s is not long enough
  if (s.size() < subs.size())
    return false;
 
  // get a substring
  std::string tail = s.substr(s.size() - subs.size());
 
  return tail == subs;
}
 
static bool checkIntersection(std::vector<std::string> v1, std::vector<std::string> v2) {
  // Sort
  std::sort(v1.begin(), v1.end());
  std::sort(v2.begin(), v2.end());
 
  // Check intersectiom
  std::vector<std::string> intersectvec(v1.size() + v2.size());
  auto outiter = std::set_intersection(v1.begin(), v1.end(), v2.begin(), v2.end(), intersectvec.begin());
  return static_cast<int>(outiter - intersectvec.begin()) != 0;
}
 
//----------------------------------------------------------------------------------------------
const std::string LoadSpiceAscii::name() const { return "LoadSpiceAscii"; }
 
//----------------------------------------------------------------------------------------------
int LoadSpiceAscii::version() const { return 1; }
 
//----------------------------------------------------------------------------------------------
const std::string LoadSpiceAscii::category() const { return "DataHandling\\Text"; }
 
//----------------------------------------------------------------------------------------------
const std::string LoadSpiceAscii::summary() const { return "Load Spice data to workspaces in general."; }
 
//----------------------------------------------------------------------------------------------
void LoadSpiceAscii::init() {
  declareProperty(std::make_unique<FileProperty>("Filename", "", API::FileProperty::Load, ".dat"),
                  "Name of SPICE data file.");
 
  // Logs to be float type sample log
  auto floatspckeyprop = std::make_unique<ArrayProperty<std::string>>("FloatSampleLogNames", Direction::Input);
  declareProperty(std::move(floatspckeyprop), "List of log names that will be imported as float property.");
 
  // Logs to be integer type sample log
  auto intspckeyprop = std::make_unique<ArrayProperty<std::string>>("IntegerSampleLogNames", Direction::Input);
  declareProperty(std::move(intspckeyprop), "List of log names that will be imported as integer property.");
 
  // Logs to be string type sample log
  auto strspckeyprop = std::make_unique<ArrayProperty<std::string>>("StringSampleLogNames", Direction::Input);
  declareProperty(std::move(strspckeyprop), "List of log names that will be imported as string property.");
 
  declareProperty("IgnoreUnlistedLogs", false,
                  "If it is true, all log names are not listed in any of above "
                  "3 input lists will be ignored. "
                  "Otherwise, any log name is not listed will be treated as "
                  "string property.");
 
  // date: MM/DD/YYYY,time: HH:MM:SS AM is the standard SPICE date time log name
  // and format
  std::vector<std::string> defaultlogformat(4);
  defaultlogformat[0] = "date";
  defaultlogformat[1] = "MM/DD/YYYY";
  defaultlogformat[2] = "time";
  defaultlogformat[3] = "HH:MM:SS AM";
  declareProperty(std::make_unique<ArrayProperty<std::string>>("DateAndTimeLog", std::move(defaultlogformat)),
                  "Name and format for date and time");
 
  // Output
  declareProperty(std::make_unique<WorkspaceProperty<ITableWorkspace>>("OutputWorkspace", "", Direction::Output),
                  "Name of TableWorkspace containing experimental data.");
 
  declareProperty(std::make_unique<WorkspaceProperty<MatrixWorkspace>>("RunInfoWorkspace", "", Direction::Output),
                  "Name of TableWorkspace containing experimental information.");
}
 
//----------------------------------------------------------------------------------------------
void LoadSpiceAscii::exec() {
  // Input properties and validate
  std::string filename = getPropertyValue("Filename");
  std::vector<std::string> strlognames = getProperty("StringSampleLogNames");
  std::vector<std::string> intlognames = getProperty("IntegerSampleLogNames");
  std::vector<std::string> floatlognames = getProperty("FloatSampleLogNames");
  bool ignoreunlisted = getProperty("IgnoreUnlistedLogs");
  std::vector<std::string> datetimeprop = getProperty("DateAndTimeLog");
 
  bool valid = validateLogNamesType(floatlognames, intlognames, strlognames);
  if (!valid)
    throw std::runtime_error("At one log name appears in multiple log type lists");
 
  // Parse
  std::vector<std::vector<std::string>> datalist;
  std::vector<std::string> titles;
  std::map<std::string, std::string> runinfodict;
  parseSPICEAscii(filename, datalist, titles, runinfodict);
 
  // Build output workspaces
  API::ITableWorkspace_sptr outws = createDataWS(datalist, titles);
 
  // Build run information workspace
  API::MatrixWorkspace_sptr runinfows =
      createRunInfoWS(runinfodict, floatlognames, intlognames, strlognames, ignoreunlisted);
 
  // Process date and time for run start explicitly
  setupRunStartTime(runinfows, datetimeprop);
 
  // Set properties
  setProperty("OutputWorkspace", outws);
  setProperty("RunInfoWorkspace", runinfows);
}
 
//----------------------------------------------------------------------------------------------
bool LoadSpiceAscii::validateLogNamesType(const std::vector<std::string> &floatlognames,
                                          const std::vector<std::string> &intlognames,
                                          const std::vector<std::string> &strlognames) {
  std::vector<std::vector<std::string>> vec_lognamelist;
  vec_lognamelist.emplace_back(floatlognames);
  vec_lognamelist.emplace_back(intlognames);
  vec_lognamelist.emplace_back(strlognames);
 
  // Check whther there is any intersction among 3 sets
  bool hascommon = false;
  for (size_t i = 0; i < 3; ++i)
    for (size_t j = i + 1; j < 3; ++j) {
      hascommon = checkIntersection(vec_lognamelist[i], vec_lognamelist[j]);
      if (hascommon) {
        std::stringstream ess;
        ess << "logsets[" << i << "] and log sets[" << j << "] has intersection.";
        g_log.error(ess.str());
        break;
      }
    }
 
  return (!hascommon);
}
 
//----------------------------------------------------------------------------------------------
void LoadSpiceAscii::parseSPICEAscii(const std::string &filename, std::vector<std::vector<std::string>> &datalist,
                                     std::vector<std::string> &titles,
                                     std::map<std::string, std::string> &runinfodict) {
  // Import file
  std::ifstream spicefile(filename.c_str());
  if (!spicefile.is_open()) {
    std::stringstream ess;
    ess << "File  " << filename << " cannot be opened.";
    throw std::runtime_error(ess.str());
  }
 
  std::string line;
  while (std::getline(spicefile, line)) {
    // Parse one line
 
    // Strip
    boost::trim(line);
    // skip for empyt line
    if (line.empty())
      continue;
 
    // Comment line for run information
    if (line[0] == '#') {
      // remove comment flag # and trim space
      line.erase(0, 1);
      boost::trim(line);
 
      if (line.find('=') != std::string::npos) {
        // run information line
        std::vector<std::string> terms;
        boost::split(terms, line, boost::is_any_of("="));
        boost::trim(terms[0]);
        g_log.debug() << "Title = " << terms[0] << ", number of splitted terms = " << terms.size() << "\n";
        std::string infovalue;
        if (terms.size() == 2) {
          infovalue = terms[1];
          boost::trim(infovalue);
        } else if (terms.size() > 2) {
          // Content contains '='
          for (size_t j = 1; j < terms.size(); ++j) {
            if (j > 1)
              infovalue += "=";
            infovalue += terms[j];
          }
        } else {
          std::stringstream wss;
          wss << "Line '" << line << "' is hard to parse.  It has more than 1 '='.";
          g_log.warning(wss.str());
        }
        runinfodict.emplace(terms[0], infovalue);
      } else if (line.find("Pt.") != std::string::npos) {
        // Title line
        boost::split(titles, line, boost::is_any_of("\t\n "), boost::token_compress_on);
      } else if (endswith(line, "scan completed.")) {
        std::vector<std::string> terms;
        boost::iter_split(terms, line, boost::algorithm::first_finder("scan completed."));
        std::string time = terms.front();
        boost::trim(time);
        runinfodict.emplace("runend", time);
      } else {
        // Not supported
        std::stringstream wss;
        wss << "File " << filename << ": line \"" << line << "\" cannot be parsed. It is ignored then.";
        g_log.warning(wss.str());
      }
    } // If for run info
    else {
      // data line
      std::vector<std::string> terms;
      boost::split(terms, line, boost::is_any_of(" \t\n"), boost::token_compress_on);
      datalist.emplace_back(terms);
    }
  }
 
  g_log.debug() << "Run info dictionary has " << runinfodict.size() << " entries."
                << "\n";
}
 
//----------------------------------------------------------------------------------------------
API::ITableWorkspace_sptr LoadSpiceAscii::createDataWS(const std::vector<std::vector<std::string>> &datalist,
                                                       const std::vector<std::string> &titles) {
  // Create a table workspace with columns defined
  DataObjects::TableWorkspace_sptr outws = std::make_shared<DataObjects::TableWorkspace>();
  size_t ipt = -1;
  for (size_t i = 0; i < titles.size(); ++i) {
    if (titles[i] == "Pt.") {
      outws->addColumn("int", titles[i]);
      ipt = i;
    } else {
      outws->addColumn("double", titles[i]);
    }
  }
 
  // Add rows
  size_t numrows = datalist.size();
  size_t numcols = outws->columnCount();
  for (size_t irow = 0; irow < numrows; ++irow) {
    TableRow newrow = outws->appendRow();
    for (size_t icol = 0; icol < numcols; ++icol) {
      std::string item = datalist[irow][icol];
      if (icol == ipt)
        newrow << std::stoi(item);
      else
        newrow << std::stod(item);
    }
  }
 
  ITableWorkspace_sptr tablews = std::dynamic_pointer_cast<ITableWorkspace>(outws);
  return tablews;
}
 
//----------------------------------------------------------------------------------------------
API::MatrixWorkspace_sptr LoadSpiceAscii::createRunInfoWS(std::map<std::string, std::string> runinfodict,
                                                          std::vector<std::string> &floatlognamelist,
                                                          std::vector<std::string> &intlognamelist,
                                                          std::vector<std::string> &strlognamelist,
                                                          bool ignoreunlisted) {
  // Create an empty workspace
  API::MatrixWorkspace_sptr infows = WorkspaceFactory::Instance().create("Workspace2D", 1, 2, 1);
 
  // Sort
  std::sort(floatlognamelist.begin(), floatlognamelist.end());
  std::sort(intlognamelist.begin(), intlognamelist.end());
  std::sort(strlognamelist.begin(), strlognamelist.end());
 
  // Create sample log properties
  std::map<std::string, std::string>::iterator miter;
  for (miter = runinfodict.begin(); miter != runinfodict.end(); ++miter) {
    const std::string title = miter->first;
    const std::string strvalue = miter->second;
 
    g_log.debug() << "Trying to add property " << title << " with value " << strvalue << "\n";
 
    if (std::binary_search(floatlognamelist.begin(), floatlognamelist.end(), title)) {
      // Case as a double property
      bool adderrorvalue = false;
      double value, error;
 
      // Convert to float value and error (if exists)
      if (strvalue.find("+/-") != std::string::npos) {
        adderrorvalue = true;
 
        std::vector<std::string> terms;
        boost::iter_split(terms, strvalue, boost::algorithm::first_finder("+/-"));
        value = std::stod(terms[0]);
        error = std::stod(terms[1]);
      } else {
        value = std::stod(strvalue);
        error = 0;
      }
 
      // Add properties
      addProperty<double>(infows, title, value);
      if (adderrorvalue) {
        std::stringstream tss;
        tss << title << ".error";
        addProperty<double>(infows, tss.str(), error);
      }
    } else if (std::binary_search(intlognamelist.begin(), intlognamelist.end(), title)) {
      // It is an integer log
      addProperty<int>(infows, title, std::stoi(strvalue));
    } else if (!ignoreunlisted || std::binary_search(strlognamelist.begin(), strlognamelist.end(), title)) {
      // It is a string log or it is not defined but not ignored either
      addProperty<std::string>(infows, title, strvalue);
    }
  }
 
  return infows;
}
 
//----------------------------------------------------------------------------------------------
void LoadSpiceAscii::setupRunStartTime(const API::MatrixWorkspace_sptr &runinfows,
                                       const std::vector<std::string> &datetimeprop) {
  // Check if no need to process run start time
  if (datetimeprop.empty()) {
    g_log.information("User chooses not to set up run start date and time.");
    return;
  }
 
  // Parse property vector
  if (datetimeprop.size() != 4) {
    g_log.warning() << "Run start date and time property must contain 4 "
                       "strings.  User only specifies "
                    << datetimeprop.size() << ".  Set up failed."
                    << "\n";
    return;
  }
 
  // Parse
  std::string datelogname = datetimeprop[0];
  std::string timelogname = datetimeprop[2];
  if (!(runinfows->run().hasProperty(datelogname) && runinfows->run().hasProperty(timelogname))) {
    std::stringstream errss;
    errss << "Unable to locate user specified date and time sample logs " << datelogname << " and " << timelogname
          << "."
          << "run_start will not be set up.";
    g_log.error(errss.str());
    return;
  }
 
  const std::string &rawdatestring = runinfows->run().getProperty(datelogname)->value();
  const std::string &dateformat = datetimeprop[1];
  std::string mtddatestring = processDateString(rawdatestring, dateformat);
 
  const std::string &rawtimestring = runinfows->run().getProperty(timelogname)->value();
  const std::string &timeformat = datetimeprop[3];
  std::string mtdtimestring = processTimeString(rawtimestring, timeformat);
 
  std::string mtddatetimestr = mtddatestring + "T" + mtdtimestring;
 
  // Set up property
  DateAndTime runstart(mtddatetimestr);
  addProperty<std::string>(runinfows, "run_start", runstart.toISO8601String());
}
 
//----------------------------------------------------------------------------------------------
std::string LoadSpiceAscii::processDateString(const std::string &rawdate, const std::string &dateformat) {
  // Identify splitter
  std::string splitter;
  if (dateformat.find('/') != std::string::npos)
    splitter += "/";
  else if (dateformat.find('-') != std::string::npos)
    splitter += "-";
  else if (dateformat.find('.') != std::string::npos)
    splitter += ".";
  else
    throw std::runtime_error("Input date format does not contain any of / - "
                             "or '.'.  Format unsupported.");
 
  // Split
  std::vector<std::string> dateterms;
  std::vector<std::string> formatterms;
  boost::split(dateterms, rawdate, boost::is_any_of(splitter));
  boost::split(formatterms, dateformat, boost::is_any_of(splitter));
 
  if (dateterms.size() != formatterms.size() || dateterms.size() != 3)
    throw std::runtime_error("Unsupported date string and format");
  std::string year;
  std::string month;
  std::string day;
  for (size_t i = 0; i < 3; ++i) {
    if (formatterms[i].find('Y') != std::string::npos)
      year = dateterms[i];
    else if (formatterms[i].find('M') != std::string::npos) {
      month = dateterms[i];
      if (month.size() == 1)
        month.insert(0, 1, '0');
    } else {
      day = dateterms[i];
      if (day.size() == 1)
        day.insert(0, 1, '0');
    }
  }
 
  std::string formatdate = year + "-" + month + "-" + day;
 
  return formatdate;
}
 
//----------------------------------------------------------------------------------------------
std::string LoadSpiceAscii::processTimeString(const std::string &rawtime, const std::string &timeformat) {
  // Process time format to find out it is 12 hour or 24 hour format
  std::string timeformatcpy(timeformat);
  boost::trim(timeformatcpy);
 
  int format;
  if (timeformatcpy.find(' ') == std::string::npos)
    format = 24;
  else
    format = 12;
 
  // Process
  std::string mtdtime;
  if (format == 24)
    mtdtime = rawtime;
  else {
    std::vector<std::string> terms;
    boost::split(terms, rawtime, boost::is_any_of(" "));
    bool pm = false;
    if (terms[1] == "PM")
      pm = true;
 
    std::vector<std::string> terms2;
    boost::split(terms2, terms[0], boost::is_any_of(":"));
    int hour = std::stoi(terms[0]);
    if (hour < 12 && pm)
      hour += 12;
 
    std::stringstream hourss;
    hourss << hour;
    std::string hourstr = hourss.str();
    if (hourstr.size() == 1)
      hourstr = "0" + hourstr;
    std::string minstr = terms2[1];
    if (minstr.size() == 1)
      minstr = "0" + minstr;
    std::string secstr = terms2[2];
    if (secstr.size() == 1)
      secstr = "0" + secstr;
 
    mtdtime = hourstr + ":" + minstr + ":" + secstr;
  }
 
  return mtdtime;
}
 
//----------------------------------------------------------------------------------------------
template <typename T>
void LoadSpiceAscii::addProperty(const API::MatrixWorkspace_sptr &ws, const std::string &pname, T pvalue) {
  ws->mutableRun().addLogData(new PropertyWithValue<T>(pname, pvalue));
}
 
} // namespace Mantid::DataHandling