Commit 96413ecf authored by Janik Zikovsky's avatar Janik Zikovsky
Browse files

Refs #1521: Some minor refactoring for clarity.

parent bf76a6e7
......@@ -90,6 +90,9 @@ public:
IndexToIndexMap * getSpectrumToWorkspaceIndexMap() const;
IndexToIndexMap * getWorkspaceIndexToDetectorIDMap() const;
IndexToIndexMap * getDetectorIDToWorkspaceIndexMap() const;
void getIndicesFromSpectra(const std::vector<int>& spectraList, std::vector<int>& indexList) const;
/// Sample accessors
......
......@@ -69,9 +69,6 @@ struct DLLExport WorkspaceHelpers
static bool sharedXData(const MatrixWorkspace_const_sptr WS);
// Divides the data in a workspace by the bin width to make it a distribution (or the reverse)
static void makeDistribution(MatrixWorkspace_sptr workspace, const bool forwards = true);
// Convert a list of spectrum numbers to the corresponding workspace indices
static void getIndicesFromSpectra(const MatrixWorkspace_const_sptr WS, const std::vector<int>& spectraList,
std::vector<int>& indexList);
};
} // namespace API
......
......@@ -236,6 +236,42 @@ IndexToIndexMap * MatrixWorkspace::getWorkspaceIndexToDetectorIDMap() const
}
//---------------------------------------------------------------------------------------
/** Converts a list of spectrum numbers to the corresponding workspace indices.
* Not a very efficient operation, but unfortunately it's sometimes required.
* @param WS The workspace on which to carry out the operation
* @param spectraList The list of spectrum numbers required
* @param indexList Returns a reference to the vector of indices (empty if not a Workspace2D)
*/
void MatrixWorkspace::getIndicesFromSpectra(const std::vector<int>& spectraList, std::vector<int>& indexList) const
{
// Clear the output index list
indexList.clear();
indexList.reserve(this->getNumberHistograms());
// get the spectra axis
SpectraAxis *spectraAxis;
if (this->axes() == 2)
{
spectraAxis = dynamic_cast<SpectraAxis*>(this->getAxis(1));
if (!spectraAxis) return;
}
// Just return an empty list if this isn't a Workspace2D
else return;
std::vector<int>::const_iterator iter = spectraList.begin();
while( iter != spectraList.end() )
{
for (int i = 0; i < this->getNumberHistograms(); ++i)
{
if ( spectraAxis->spectraNo(i) == *iter )
{
indexList.push_back(i);
}
}
++iter;
}
}
//---------------------------------------------------------------------------------------
......
......@@ -384,41 +384,5 @@ void WorkspaceHelpers::makeDistribution(MatrixWorkspace_sptr workspace, const bo
workspace->isDistribution(forwards);
}
/** Converts a list of spectrum numbers to the corresponding workspace indices.
* Not a very efficient operation, but unfortunately it's sometimes required.
* @param WS The workspace on which to carry out the operation
* @param spectraList The list of spectrum numbers required
* @param indexList Returns a reference to the vector of indices (empty if not a Workspace2D)
*/
void WorkspaceHelpers::getIndicesFromSpectra(const MatrixWorkspace_const_sptr WS, const std::vector<int>& spectraList,
std::vector<int>& indexList)
{
// Clear the output index list
indexList.clear();
indexList.reserve(WS->getNumberHistograms());
// get the spectra axis
SpectraAxis *spectraAxis;
if (WS->axes() == 2)
{
spectraAxis = dynamic_cast<SpectraAxis*>(WS->getAxis(1));
if (!spectraAxis) return;
}
// Just return an empty list if this isn't a Workspace2D
else return;
std::vector<int>::const_iterator iter = spectraList.begin();
while( iter != spectraList.end() )
{
for (int i = 0; i < WS->getNumberHistograms(); ++i)
{
if ( spectraAxis->spectraNo(i) == *iter )
{
indexList.push_back(i);
}
}
++iter;
}
}
} // namespace API
} // namespace Mantid
......@@ -79,7 +79,7 @@ void CalculateTransmission::exec()
udets.push_back(getProperty("TransmissionMonitor"));
// Convert UDETs to workspace indices via spectrum numbers
const std::vector<int> sampleSpectra = sampleWS->spectraMap().getSpectra(udets);
WorkspaceHelpers::getIndicesFromSpectra(sampleWS,sampleSpectra,indices);
sampleWS->getIndicesFromSpectra(sampleSpectra,indices);
if (indices.size() < 2)
{
if (indices.size() == 1)
......@@ -88,7 +88,7 @@ void CalculateTransmission::exec()
}
else
{
g_log.debug() << "WorkspaceHelpers::getIndicesFromSpectra() returned empty\n";
g_log.debug() << "sampleWS->getIndicesFromSpectra() returned empty\n";
}
throw std::invalid_argument("Could not find the incident and transmission monitor spectra\n");
}
......@@ -104,7 +104,7 @@ void CalculateTransmission::exec()
MatrixWorkspace_sptr M2_sample = this->extractSpectrum(sampleWS,indices[0]);
MatrixWorkspace_sptr M3_sample = this->extractSpectrum(sampleWS,indices[1]);
const std::vector<int> directSpectra = directWS->spectraMap().getSpectra(udets);
WorkspaceHelpers::getIndicesFromSpectra(sampleWS,directSpectra,indices);
sampleWS->getIndicesFromSpectra(directSpectra,indices);
// Check that given spectra are monitors
if ( !directWS->getDetector(indices.front())->isMonitor() )
{
......
......@@ -95,7 +95,7 @@ void CalculateTransmissionBeamSpreader::exec()
std::vector<int> spectra = sample_scatterWS->spectraMap().getSpectra(udets);
// Get monitors (assume that the detector mapping is the same for all data sets)
WorkspaceHelpers::getIndicesFromSpectra(sample_scatterWS, spectra, indices);
sample_scatterWS->getIndicesFromSpectra(spectra, indices);
if (indices.size() != 1)
{
g_log.error() << "Could not find the incident monitor spectra\n";
......
......@@ -159,7 +159,7 @@ std::vector<int> GetEi::getMonitorSpecIndexs(API::MatrixWorkspace_const_sptr WS,
// get the index number of the histogram for the first monitor
std::vector<int> specNumTemp(&specNum1, &specNum1+1);
WorkspaceHelpers::getIndicesFromSpectra(WS, specNumTemp, specInds);
WS->getIndicesFromSpectra(specNumTemp, specInds);
if ( specInds.size() != 1 )
{// the monitor spectrum isn't present in the workspace, we can't continue from here
g_log.error() << "Couldn't find the first monitor spectrum, number " << specNum1 << std::endl;
......@@ -169,7 +169,7 @@ std::vector<int> GetEi::getMonitorSpecIndexs(API::MatrixWorkspace_const_sptr WS,
// nowe the second monitor
std::vector<int> specIndexTemp;
specNumTemp[0] = specNum2;
WorkspaceHelpers::getIndicesFromSpectra(WS, specNumTemp, specIndexTemp);
WS->getIndicesFromSpectra(specNumTemp, specIndexTemp);
if ( specIndexTemp.size() != 1 )
{// the monitor spectrum isn't present in the workspace, we can't continue from here
g_log.error() << "Couldn't find the second monitor spectrum, number " << specNum2 << std::endl;
......
......@@ -105,7 +105,7 @@ double GetEi2::calculateEi(const double initial_guess)
std::vector<int> mon_indices;
mon_indices.reserve(2);
// get the index number of the histogram for the first monitor
WorkspaceHelpers::getIndicesFromSpectra(m_input_ws, spec_nums, mon_indices);
m_input_ws->getIndicesFromSpectra(spec_nums, mon_indices);
if( mon_indices.size() != 2 )
{
......
......@@ -73,14 +73,14 @@ void GroupDetectors::exec()
// appropriate spectra number and adding the indices they are linked to the list to be processed
if ( ! spectraList.empty() )
{
WorkspaceHelpers::getIndicesFromSpectra(WS,spectraList,indexList);
WS->getIndicesFromSpectra(spectraList,indexList);
}// End dealing with spectraList
else if ( ! detectorList.empty() )
{// Dealing with DetectorList
//convert from detectors to spectra numbers
std::vector<int> mySpectraList = WS->spectraMap().getSpectra(detectorList);
//then from spectra numbers to indices
WorkspaceHelpers::getIndicesFromSpectra(WS,mySpectraList,indexList);
WS->getIndicesFromSpectra(mySpectraList,indexList);
}
if ( indexList.size() == 0 )
......
......@@ -183,8 +183,7 @@ void GroupDetectors2::getGroups(API::MatrixWorkspace_const_sptr workspace,
// only look at these other parameters if the file wasn't set
if ( ! spectraList.empty() )
{
WorkspaceHelpers::getIndicesFromSpectra(workspace,
spectraList, m_GroupSpecInds[0]);
workspace->getIndicesFromSpectra( spectraList, m_GroupSpecInds[0]);
g_log.debug() << "Converted " << spectraList.size() << " spectra numbers into spectra indices to be combined\n";
}
else
......@@ -193,8 +192,7 @@ void GroupDetectors2::getGroups(API::MatrixWorkspace_const_sptr workspace,
{// we are going to group on the basis of detector IDs, convert from detectors to spectra numbers
std::vector<int> mySpectraList = workspace->spectraMap().getSpectra(detectorList);
//then from spectra numbers to indices
WorkspaceHelpers::getIndicesFromSpectra(
workspace, mySpectraList, m_GroupSpecInds[0]);
workspace->getIndicesFromSpectra( mySpectraList, m_GroupSpecInds[0]);
g_log.debug() << "Found " << m_GroupSpecInds[0].size() << " spectra indices from the list of " << detectorList.size() << " detectors\n";
}
else if ( ! indexList.empty() )
......
......@@ -86,11 +86,11 @@ public:
void testSetup()
{
GroupDetectors2 gd;
TS_ASSERT_EQUALS( gd.name(), "GroupDetectors" )
TS_ASSERT_EQUALS( gd.version(), 2 )
TS_ASSERT_EQUALS( gd.category(), "DataHandling\\Detectors" )
TS_ASSERT_THROWS_NOTHING( gd.initialize() )
TS_ASSERT( gd.isInitialized() )
TS_ASSERT_EQUALS( gd.name(), "GroupDetectors" );
TS_ASSERT_EQUALS( gd.version(), 2 );
TS_ASSERT_EQUALS( gd.category(), "DataHandling\\Detectors" );
TS_ASSERT_THROWS_NOTHING( gd.initialize() );
TS_ASSERT( gd.isInitialized() );
gd.setPropertyValue("InputWorkspace", inputWS);
gd.setPropertyValue("OutputWorkspace", outputBase);
......@@ -114,20 +114,20 @@ public:
MatrixWorkspace_sptr outputWS = boost::dynamic_pointer_cast<MatrixWorkspace>(
AnalysisDataService::Instance().retrieve(output));
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 )
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 );
std::vector<double> tens(NBINS+1, 10.0);
std::vector<double> ones(NBINS, 1.0);
TS_ASSERT_EQUALS( outputWS->dataX(0), tens )
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, 1+4) )
TS_ASSERT_EQUALS( outputWS->dataX(0), tens );
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, 1+4) );
for (int i = 0; i < NBINS; ++i)
{
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(double(2)), 0.0001 )
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(double(2)), 0.0001 );
}
boost::shared_ptr<IDetector> det;
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) )
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) )
TS_ASSERT_THROWS_ANYTHING(det = outputWS->getDetector(1))
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) );
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) );
TS_ASSERT_THROWS_ANYTHING(det = outputWS->getDetector(1));
AnalysisDataService::Instance().remove(output);
}
......@@ -142,7 +142,7 @@ public:
// test the algorithm behaves if you give it a non-existent index
grouper3.setPropertyValue("WorkspaceIndexList","4-6");
grouper3.execute();
TS_ASSERT( ! grouper3.isExecuted() )
TS_ASSERT( ! grouper3.isExecuted() );
grouper3.setPropertyValue("WorkspaceIndexList","2-5");
TS_ASSERT_THROWS_NOTHING( grouper3.execute());
......@@ -150,20 +150,20 @@ public:
MatrixWorkspace_sptr outputWS = boost::dynamic_pointer_cast<MatrixWorkspace>(
AnalysisDataService::Instance().retrieve(output));
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 )
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 );
std::vector<double> tens(NBINS+1, 10.0);
std::vector<double> ones(NBINS, 1.0);
TS_ASSERT_EQUALS( outputWS->dataX(0), tens )
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, (3+4+5+6)) )
TS_ASSERT_EQUALS( outputWS->dataX(0), tens );
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, (3+4+5+6)) );
for (int i = 0; i < NBINS; ++i)
{
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(4.0), 0.0001 )
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(4.0), 0.0001 );
}
boost::shared_ptr<IDetector> det;
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) )
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) )
TS_ASSERT_THROWS_ANYTHING(det = outputWS->getDetector(1))
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) );
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) );
TS_ASSERT_THROWS_ANYTHING(det = outputWS->getDetector(1));
AnalysisDataService::Instance().remove(output);
}
......@@ -183,21 +183,21 @@ public:
MatrixWorkspace_const_sptr outputWS = boost::dynamic_pointer_cast<MatrixWorkspace>(
AnalysisDataService::Instance().retrieve(output));
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 )
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), 1 );
std::vector<double> tens(NBINS+1, 10.0);
std::vector<double> ones(NBINS, 1.0);
TS_ASSERT_EQUALS( outputWS->dataX(0), tens )
TS_ASSERT_EQUALS( outputWS->dataX(0), tens );
TS_ASSERT_EQUALS( outputWS->dataY(0),
std::vector<double>(NBINS, (3+1)+(1+1)+(4+1)+(0+1)+(2+1)+(5+1)) )
std::vector<double>(NBINS, (3+1)+(1+1)+(4+1)+(0+1)+(2+1)+(5+1)) );
for (int i = 0; i < NBINS; ++i)
{// assume that we have grouped all the spectra in the input workspace
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(double(NHIST)), 0.0001 )
TS_ASSERT_DELTA( outputWS->dataE(0)[i], std::sqrt(double(NHIST)), 0.0001 );
}
boost::shared_ptr<IDetector> det;
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) )
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) )
TS_ASSERT_THROWS_ANYTHING( det = outputWS->getDetector(1) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) );
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) );
TS_ASSERT_THROWS_ANYTHING( det = outputWS->getDetector(1) );
AnalysisDataService::Instance().remove(output);
}
......@@ -220,49 +220,49 @@ public:
MatrixWorkspace_const_sptr outputWS = boost::dynamic_pointer_cast<MatrixWorkspace>(
AnalysisDataService::Instance().retrieve(output));
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), NHIST-1 )
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), NHIST-1 );
std::vector<double> tens(NBINS+1, 10.0);
std::vector<double> ones(NBINS, 1.0);
// check the two grouped spectra
TS_ASSERT_EQUALS( outputWS->dataX(0), tens )
TS_ASSERT_EQUALS( outputWS->dataX(0), tens );
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, 1+3) );
for (int i = 0; i < NBINS; ++i)
{
TS_ASSERT_DELTA(outputWS->dataE(0)[i], std::sqrt(static_cast<double>(2)), 1e-6)
TS_ASSERT_DELTA(outputWS->dataE(0)[i], std::sqrt(static_cast<double>(2)), 1e-6);
}
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(0), 1 )
TS_ASSERT_EQUALS( outputWS->dataX(1), tens )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(0), 1 );
TS_ASSERT_EQUALS( outputWS->dataX(1), tens );
TS_ASSERT_EQUALS( outputWS->dataY(1), std::vector<double>(NBINS, 4 ) );
TS_ASSERT_EQUALS( outputWS->dataE(1), ones )
TS_ASSERT_EQUALS( outputWS->dataE(1), ones );
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(1), 4 )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(1), 4 );
//check the unmoved spectra
TS_ASSERT_EQUALS( outputWS->dataX(2), tens )
TS_ASSERT_EQUALS( outputWS->dataY(2), std::vector<double>(NBINS, 2) )
TS_ASSERT_EQUALS( outputWS->dataE(2), ones )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(2), 2 )
TS_ASSERT_EQUALS( outputWS->dataX(3), tens )
TS_ASSERT_EQUALS( outputWS->dataY(3), std::vector<double>(NBINS, 5) )
TS_ASSERT_EQUALS( outputWS->dataE(3), ones )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(3), 5 )
TS_ASSERT_EQUALS( outputWS->dataY(4), std::vector<double>(NBINS, 6) )
TS_ASSERT_EQUALS( outputWS->dataE(4), ones )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(4), 6 )
TS_ASSERT_EQUALS( outputWS->dataX(2), tens );
TS_ASSERT_EQUALS( outputWS->dataY(2), std::vector<double>(NBINS, 2) );
TS_ASSERT_EQUALS( outputWS->dataE(2), ones );
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(2), 2 );
TS_ASSERT_EQUALS( outputWS->dataX(3), tens );
TS_ASSERT_EQUALS( outputWS->dataY(3), std::vector<double>(NBINS, 5) );
TS_ASSERT_EQUALS( outputWS->dataE(3), ones );
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(3), 5 );
TS_ASSERT_EQUALS( outputWS->dataY(4), std::vector<double>(NBINS, 6) );
TS_ASSERT_EQUALS( outputWS->dataE(4), ones );
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(4), 6 );
// the first two spectra should have a group of detectors the other spectra a single detector
boost::shared_ptr<IDetector> det;
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) )
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(1) )
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(2) )
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(3) )
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(4) )
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) )
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(0) );
TS_ASSERT( boost::dynamic_pointer_cast<DetectorGroup>(det) );
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(1) );
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) );
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(2) );
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) );
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(3) );
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) );
TS_ASSERT_THROWS_NOTHING( det = outputWS->getDetector(4) );
TS_ASSERT( boost::dynamic_pointer_cast<Detector>(det) );
AnalysisDataService::Instance().remove(output);
remove(inputFile.c_str());
......@@ -286,31 +286,31 @@ public:
MatrixWorkspace_const_sptr outputWS = boost::dynamic_pointer_cast<MatrixWorkspace>(
AnalysisDataService::Instance().retrieve(output));
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), NHIST-3 )
TS_ASSERT_EQUALS( outputWS->getNumberHistograms(), NHIST-3 );
std::vector<double> tens(NBINS+1, 10.0);
std::vector<double> ones(NBINS, 1.0);
// check the first grouped spectrum
TS_ASSERT_EQUALS( outputWS->dataX(0), tens )
TS_ASSERT_EQUALS( outputWS->dataX(0), tens );
TS_ASSERT_EQUALS( outputWS->dataY(0), std::vector<double>(NBINS, 1+2+3) );
for (int i = 0; i < NBINS; ++i)
{
TS_ASSERT_DELTA(outputWS->dataE(0)[i], std::sqrt(static_cast<double>(3)), 1e-6)
TS_ASSERT_DELTA(outputWS->dataE(0)[i], std::sqrt(static_cast<double>(3)), 1e-6);
}
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(0), 1 )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(0), 1 );
// check the second grouped spectrum
TS_ASSERT_EQUALS( outputWS->dataX(1), tens )
TS_ASSERT_EQUALS( outputWS->dataX(1), tens );
TS_ASSERT_EQUALS( outputWS->dataY(1), std::vector<double>(NBINS, 4 ) );
TS_ASSERT_EQUALS( outputWS->dataE(1), ones )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(1), 4 )
TS_ASSERT_EQUALS( outputWS->dataE(1), ones );
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(1), 4 );
// check the third grouped spectrum
TS_ASSERT_EQUALS( outputWS->dataX(2), tens )
TS_ASSERT_EQUALS( outputWS->dataX(2), tens );
TS_ASSERT_EQUALS( outputWS->dataY(2), std::vector<double>(NBINS, 5+6 ) );
for (int i = 0; i < NBINS; ++i)
{
TS_ASSERT_DELTA(
outputWS->dataE(2)[i], std::sqrt(static_cast<double>(2)), 1e-6)
outputWS->dataE(2)[i], std::sqrt(static_cast<double>(2)), 1e-6);
}
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(2), 5 )
TS_ASSERT_EQUALS( outputWS->getAxis(1)->spectraNo(2), 5 );
AnalysisDataService::Instance().remove(output);
remove(inputFile.c_str());
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment