Skip to content

Commit

Permalink
Add tests for the SofQW spectrum mapping.
Browse files Browse the repository at this point in the history 
Refs #12584
  • Loading branch information
@martyngigg
martyngigg committed Jun 25, 2015
1 parent 1f9d645 commit 137727e
Show file tree
Hide file tree
Showing 2 changed files with 122 additions and 81 deletions.
103 changes: 62 additions & 41 deletions Code/Mantid/Framework/Algorithms/test/SofQWNormalisedPolygonTest.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,64 +5,85 @@
#include "MantidAlgorithms/SofQWNormalisedPolygon.h"

#include "SofQWTest.h"
#include <boost/assign.hpp>

using namespace Mantid::Algorithms;
using namespace Mantid::API;

class SofQWNormalisedPolygonTest : public CxxTest::TestSuite
{
class SofQWNormalisedPolygonTest : public CxxTest::TestSuite {
public:

void test_Init()
{
void test_Init() {
SofQWNormalisedPolygon alg;
TS_ASSERT_THROWS_NOTHING( alg.initialize() )
TS_ASSERT( alg.isInitialized() )
TS_ASSERT_THROWS_NOTHING(alg.initialize())
TS_ASSERT(alg.isInitialized())
}

void test_Aliased_To_SofQW3()
{
void test_Aliased_To_SofQW3() {
SofQWNormalisedPolygon alg;
TS_ASSERT_EQUALS("SofQW3", alg.alias())
}

void test_exec()
{
auto result = SofQWTest::runSQW<Mantid::Algorithms::SofQWNormalisedPolygon>();
void test_exec() {
auto result =
SofQWTest::runSQW<Mantid::Algorithms::SofQWNormalisedPolygon>();

TS_ASSERT_EQUALS( result->getAxis(0)->length(), 1904 );
TS_ASSERT_EQUALS( result->getAxis(0)->unit()->unitID(), "DeltaE" );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(0), -0.5590, 0.0001 );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(999), -0.0971, 0.0001 );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(1900), 0.5728, 0.0001 );
TS_ASSERT_EQUALS(result->getAxis(0)->length(), 1904);
TS_ASSERT_EQUALS(result->getAxis(0)->unit()->unitID(), "DeltaE");
TS_ASSERT_DELTA((*(result->getAxis(0)))(0), -0.5590, 0.0001);
TS_ASSERT_DELTA((*(result->getAxis(0)))(999), -0.0971, 0.0001);
TS_ASSERT_DELTA((*(result->getAxis(0)))(1900), 0.5728, 0.0001);

TS_ASSERT_EQUALS( result->getAxis(1)->length(), 7 );
TS_ASSERT_EQUALS( result->getAxis(1)->unit()->unitID(), "MomentumTransfer" );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(0), 0.5 );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(3), 1.25 );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(6), 2.0 );
TS_ASSERT_EQUALS(result->getAxis(1)->length(), 7);
TS_ASSERT_EQUALS(result->getAxis(1)->unit()->unitID(), "MomentumTransfer");
TS_ASSERT_EQUALS((*(result->getAxis(1)))(0), 0.5);
TS_ASSERT_EQUALS((*(result->getAxis(1)))(3), 1.25);
TS_ASSERT_EQUALS((*(result->getAxis(1)))(6), 2.0);

const double delta(1e-08);
TS_ASSERT_DELTA( result->readY(0)[1160], 22.8567683273, delta);
TS_ASSERT_DELTA( result->readE(0)[1160], 0.2568965638, delta);

TS_ASSERT_DELTA( result->readY(1)[1145], 7.5942160104, delta);
TS_ASSERT_DELTA( result->readE(1)[1145], 0.1490079010 ,delta);

TS_ASSERT_DELTA( result->readY(2)[1200], 2.0249626546, delta);
TS_ASSERT_DELTA( result->readE(2)[1200], 0.0752776593, delta);

TS_ASSERT_DELTA( result->readY(3)[99], 0.0419939169, delta);
TS_ASSERT_DELTA( result->readE(3)[99], 0.0175106375, delta);

TS_ASSERT_DELTA( result->readY(4)[1654], 0.0167189448, delta);
TS_ASSERT_DELTA( result->readE(4)[1654], 0.0056801131, delta);

TS_ASSERT_DELTA( result->readY(5)[1025], 0.0808168496, delta);
TS_ASSERT_DELTA( result->readE(5)[1025], 0.0161117732, delta);

TS_ASSERT_DELTA(result->readY(0)[1160], 22.8567683273, delta);
TS_ASSERT_DELTA(result->readE(0)[1160], 0.2568965638, delta);

TS_ASSERT_DELTA(result->readY(1)[1145], 7.5942160104, delta);
TS_ASSERT_DELTA(result->readE(1)[1145], 0.1490079010, delta);

TS_ASSERT_DELTA(result->readY(2)[1200], 2.0249626546, delta);
TS_ASSERT_DELTA(result->readE(2)[1200], 0.0752776593, delta);

TS_ASSERT_DELTA(result->readY(3)[99], 0.0419939169, delta);
TS_ASSERT_DELTA(result->readE(3)[99], 0.0175106375, delta);

TS_ASSERT_DELTA(result->readY(4)[1654], 0.0167189448, delta);
TS_ASSERT_DELTA(result->readE(4)[1654], 0.0056801131, delta);

TS_ASSERT_DELTA(result->readY(5)[1025], 0.0808168496, delta);
TS_ASSERT_DELTA(result->readE(5)[1025], 0.0161117732, delta);

// Spectra-detector mapping
const size_t nspectra(6);
using namespace boost::assign;
typedef std::set<int> IDSet;
std::vector<IDSet> expectedIDs(nspectra);
IDSet s1 = list_of(3);
expectedIDs[0] = s1;
IDSet s2 = list_of(13);
expectedIDs[1] = s2;
IDSet s3 = list_of(13)(23);
expectedIDs[2] = s3;
IDSet s4 = list_of(23)(33);
expectedIDs[3] = s4;
IDSet s5 = list_of(33)(43);
expectedIDs[4] = s5;
IDSet s6 = list_of(43);
expectedIDs[5] = s6;

for (size_t i = 0; i < nspectra; ++i) {
const auto *spectrum = result->getSpectrum(i);
Mantid::specid_t specNoActual = spectrum->getSpectrumNo();
Mantid::specid_t specNoExpected = static_cast<Mantid::specid_t>(i + 1);
TS_ASSERT_EQUALS(specNoExpected, specNoActual);
TS_ASSERT_EQUALS(expectedIDs[i], spectrum->getDetectorIDs());
}
}
};


#endif /* MANTID_ALGORITHMS_SOFQW2TEST_H_ */
100 changes: 60 additions & 40 deletions Code/Mantid/Framework/Algorithms/test/SofQWPolygonTest.h
View file
Original file line number Diff line number Diff line change
Expand Up @@ -5,64 +5,84 @@
#include "MantidAlgorithms/SofQWPolygon.h"

#include "SofQWTest.h"
#include <boost/assign.hpp>

using namespace Mantid::Algorithms;
using namespace Mantid::API;

class SofQWPolygonTest : public CxxTest::TestSuite
{
class SofQWPolygonTest : public CxxTest::TestSuite {
public:

void test_Init()
{
void test_Init() {
SofQWPolygon alg;
TS_ASSERT_THROWS_NOTHING( alg.initialize() )
TS_ASSERT( alg.isInitialized() )
TS_ASSERT_THROWS_NOTHING(alg.initialize())
TS_ASSERT(alg.isInitialized())
}

void test_Aliased_To_SofQW2()
{
void test_Aliased_To_SofQW2() {
SofQWPolygon alg;
TS_ASSERT_EQUALS("SofQW2", alg.alias())
}

void test_exec()
{
void test_exec() {
auto result = SofQWTest::runSQW<Mantid::Algorithms::SofQWPolygon>();

TS_ASSERT_EQUALS( result->getAxis(0)->length(), 1904 );
TS_ASSERT_EQUALS( result->getAxis(0)->unit()->unitID(), "DeltaE" );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(0), -0.5590, 0.0001 );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(999), -0.0971, 0.0001 );
TS_ASSERT_DELTA( (*(result->getAxis(0)))(1900), 0.5728, 0.0001 );
TS_ASSERT_EQUALS(result->getAxis(0)->length(), 1904);
TS_ASSERT_EQUALS(result->getAxis(0)->unit()->unitID(), "DeltaE");
TS_ASSERT_DELTA((*(result->getAxis(0)))(0), -0.5590, 0.0001);
TS_ASSERT_DELTA((*(result->getAxis(0)))(999), -0.0971, 0.0001);
TS_ASSERT_DELTA((*(result->getAxis(0)))(1900), 0.5728, 0.0001);

TS_ASSERT_EQUALS( result->getAxis(1)->length(), 7 );
TS_ASSERT_EQUALS( result->getAxis(1)->unit()->unitID(), "MomentumTransfer" );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(0), 0.5 );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(3), 1.25 );
TS_ASSERT_EQUALS( (*(result->getAxis(1)))(6), 2.0 );
TS_ASSERT_EQUALS(result->getAxis(1)->length(), 7);
TS_ASSERT_EQUALS(result->getAxis(1)->unit()->unitID(), "MomentumTransfer");
TS_ASSERT_EQUALS((*(result->getAxis(1)))(0), 0.5);
TS_ASSERT_EQUALS((*(result->getAxis(1)))(3), 1.25);
TS_ASSERT_EQUALS((*(result->getAxis(1)))(6), 2.0);

const double delta(1e-08);
TS_ASSERT_DELTA( result->readY(0)[1160], 17.5583314826, delta);
TS_ASSERT_DELTA( result->readE(0)[1160], 0.197345265992, delta);

TS_ASSERT_DELTA( result->readY(1)[1145], 4.61301046588, delta);
TS_ASSERT_DELTA( result->readE(1)[1145], 0.0721823446635,delta);

TS_ASSERT_DELTA( result->readY(2)[1200], 1.33394133548, delta);
TS_ASSERT_DELTA( result->readE(2)[1200], 0.0419839252961, delta);

TS_ASSERT_DELTA( result->readY(3)[99], 0.0446085388561, delta);
TS_ASSERT_DELTA( result->readE(3)[99], 0.0185049423467, delta);

TS_ASSERT_DELTA( result->readY(4)[1654], 0.0171136490957 , delta);
TS_ASSERT_DELTA( result->readE(4)[1654], 0.005007299861, delta);

TS_ASSERT_DELTA( result->readY(5)[1025], 0.0516113202152, delta);
TS_ASSERT_DELTA( result->readE(5)[1025], 0.0102893133461, delta);

TS_ASSERT_DELTA(result->readY(0)[1160], 17.5583314826, delta);
TS_ASSERT_DELTA(result->readE(0)[1160], 0.197345265992, delta);

TS_ASSERT_DELTA(result->readY(1)[1145], 4.61301046588, delta);
TS_ASSERT_DELTA(result->readE(1)[1145], 0.0721823446635, delta);

TS_ASSERT_DELTA(result->readY(2)[1200], 1.33394133548, delta);
TS_ASSERT_DELTA(result->readE(2)[1200], 0.0419839252961, delta);

TS_ASSERT_DELTA(result->readY(3)[99], 0.0446085388561, delta);
TS_ASSERT_DELTA(result->readE(3)[99], 0.0185049423467, delta);

TS_ASSERT_DELTA(result->readY(4)[1654], 0.0171136490957, delta);
TS_ASSERT_DELTA(result->readE(4)[1654], 0.005007299861, delta);

TS_ASSERT_DELTA(result->readY(5)[1025], 0.0516113202152, delta);
TS_ASSERT_DELTA(result->readE(5)[1025], 0.0102893133461, delta);

// Spectra-detector mapping
const size_t nspectra(6);
using namespace boost::assign;
typedef std::set<int> IDSet;
std::vector<IDSet> expectedIDs(nspectra);
IDSet s1 = list_of(3);
expectedIDs[0] = s1;
IDSet s2 = list_of(13);
expectedIDs[1] = s2;
IDSet s3 = list_of(23);
expectedIDs[2] = s3;
IDSet s4 = list_of(23)(33);
expectedIDs[3] = s4;
IDSet s5 = list_of(33)(43);
expectedIDs[4] = s5;
IDSet s6 = list_of(43);
expectedIDs[5] = s6;

for (size_t i = 0; i < nspectra; ++i) {
const auto *spectrum = result->getSpectrum(i);
Mantid::specid_t specNoActual = spectrum->getSpectrumNo();
Mantid::specid_t specNoExpected = static_cast<Mantid::specid_t>(i + 1);
TS_ASSERT_EQUALS(specNoExpected, specNoActual);
TS_ASSERT_EQUALS(expectedIDs[i], spectrum->getDetectorIDs());
}
}
};


#endif /* MANTID_ALGORITHMS_SofQWPolygonTEST_H_ */

0 comments on commit 137727e

Please sign in to comment.