void LoadLLB::loadDataIntoTheWorkSpace(NeXus::NXEntry &entry) { // read in the data NXData dataGroup = entry.openNXData("nxdata"); NXFloat data = dataGroup.openFloatData(); data.load(); // EPP int calculatedDetectorElasticPeakPosition = getDetectorElasticPeakPosition(data); std::vector<double> timeBinning = getTimeBinning(calculatedDetectorElasticPeakPosition, m_channelWidth); // Assign time bin to first X entry m_localWorkspace->dataX(0).assign(timeBinning.begin(), timeBinning.end()); Progress progress(this, 0, 1, m_numberOfTubes * m_numberOfPixelsPerTube); size_t spec = 0; for (size_t i = 0; i < m_numberOfTubes; ++i) { for (size_t j = 0; j < m_numberOfPixelsPerTube; ++j) { if (spec > 0) { // just copy the time binning axis to every spectra m_localWorkspace->dataX(spec) = m_localWorkspace->readX(0); } // Assign Y float *data_p = &data(static_cast<int>(i), static_cast<int>(j)); m_localWorkspace->dataY(spec).assign(data_p, data_p + m_numberOfChannels); // Assign Error MantidVec &E = m_localWorkspace->dataE(spec); std::transform(data_p, data_p + m_numberOfChannels, E.begin(), LoadLLB::calculateError); ++spec; progress.report(); } } g_log.debug() << "Data loading inti WS done....\n"; }