Skip to content
Snippets Groups Projects
LoadDNSSCDTest.h 114 KiB
Newer Older
    // there are 7 points (the rest is outside of 2theta limits)
    TS_ASSERT_EQUALS(nws->getNPoints(), 7);

    AnalysisDataService::Instance().remove(normWSName);
  }

  void test_TOFWSStructure() {
    std::string outWSName("LoadDNSSCDTest_OutputWS");
    std::string normWSName("LoadDNSSCDTest_OutputWS_norm");

    LoadDNSSCD alg;
    TS_ASSERT_THROWS_NOTHING(alg.initialize());
    TS_ASSERT(alg.isInitialized());
    TS_ASSERT_THROWS_NOTHING(alg.setPropertyValue("Filenames", "dnstof.d_dat"));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("OutputWorkspace", outWSName));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("NormalizationWorkspace", normWSName));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("Normalization", "monitor"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("DeltaEmin", "-2.991993"));
    TS_ASSERT_THROWS_NOTHING(alg.execute(););
    TS_ASSERT(alg.isExecuted());

    // Retrieve the workspace from data service.
    IMDEventWorkspace_sptr iws;
    TS_ASSERT_THROWS_NOTHING(
        iws = AnalysisDataService::Instance().retrieveWS<IMDEventWorkspace>(
            outWSName));
    TS_ASSERT(iws);
    TS_ASSERT_EQUALS(iws->getNumDims(), 4);
    TS_ASSERT_EQUALS(iws->getNPoints(), 1968);
    TS_ASSERT_EQUALS(iws->id(), "MDEventWorkspace<MDEvent,4>");
    // test some metadata
    TS_ASSERT_EQUALS(iws->getNumExperimentInfo(), 1);
    ExperimentInfo_sptr expinfo = iws->getExperimentInfo(0);
    auto &run = expinfo->run();
Marina Ganeva's avatar
Marina Ganeva committed
    TimeSeriesProperty<double> *p = dynamic_cast<TimeSeriesProperty<double> *>(
        run.getProperty("TOF channels"));
    TS_ASSERT_DELTA(p->firstValue(), 100, 1.0e-05);
Marina Ganeva's avatar
Marina Ganeva committed
    p = dynamic_cast<TimeSeriesProperty<double> *>(
        run.getProperty("Time per channel"));
    TS_ASSERT_DELTA(p->firstValue(), 40.1, 1.0e-05);
    // test box controller
    BoxController_sptr bc = iws->getBoxController();
    TS_ASSERT(bc);
    TS_ASSERT_EQUALS(bc->getNumMDBoxes().size(), 6);

    // test dimensions
    std::vector<std::string> v = {"H", "K", "L", "DeltaE"};
    for (auto i = 0; i < 4; i++) {
      auto dim = iws->getDimension(i);
      TS_ASSERT(dim);
      TS_ASSERT_EQUALS(dim->getName(), v[i]);
      TS_ASSERT_EQUALS(dim->getNBins(), 5);
      double d(1.0e-05);
      TS_ASSERT_DELTA(dim->getMinimum(), -2.991993, d);
      if (i < 3) {
Marina Ganeva's avatar
Marina Ganeva committed
        TS_ASSERT_DELTA(dim->getMaximum(), 2.991993, d);
      } else {
        TS_ASSERT_DELTA(dim->getMaximum(), 4.637426, d);
      }
    }
    AnalysisDataService::Instance().remove(outWSName);
  }

  void test_TOFWSData() {
    // test whether the calculation for inelastic data are correct

    std::string outWSName("LoadDNSSCDTest_OutputWS");
    std::string normWSName("LoadDNSSCDTest_OutputWS_norm");

    LoadDNSSCD alg;
    TS_ASSERT_THROWS_NOTHING(alg.initialize());
    TS_ASSERT(alg.isInitialized());
    TS_ASSERT_THROWS_NOTHING(alg.setPropertyValue("Filenames", "dnstof.d_dat"));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("OutputWorkspace", outWSName));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("NormalizationWorkspace", normWSName));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("Normalization", "monitor"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("a", 3.55));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("b", 3.55));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("c", 24.778));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("alpha", 90.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("beta", 90.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("gamma", 120.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("OmegaOffset", 0.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("HKL1", "1,1,0"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("HKL2", "0,0,1"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("TwoThetaLimits", "20.0,55.0"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("DeltaEmin", "-3.0"));
    TS_ASSERT_THROWS_NOTHING(alg.execute(););
    TS_ASSERT(alg.isExecuted());

    // Retrieve the workspace from data service.
    IMDEventWorkspace_sptr iws;
    TS_ASSERT_THROWS_NOTHING(
        iws = AnalysisDataService::Instance().retrieveWS<IMDEventWorkspace>(
            outWSName));
    TS_ASSERT(iws);

    std::vector<API::IMDNode *> boxes(0, NULL);
    iws->getBoxes(boxes, 10000, false);
    TSM_ASSERT_EQUALS("Number of boxes", boxes.size(), 1);
    API::IMDNode *box = boxes[0];
    // there are 7 points (the rest is outside of 2theta limits)
    TS_ASSERT_EQUALS(box->getNPoints(), 574);
    std::vector<coord_t> events;
    size_t ncols;
    box->getEventsData(events, ncols);
    // 8 columns: I, err^2, run_num, det_id, h, k, l, dE
    TS_ASSERT_EQUALS(ncols, 8);
    // 8*574 = 4592
    TS_ASSERT_EQUALS(events.size(), 4592);
    // reference vector, EPP is wrong, since no elastic channel is given
Marina Ganeva's avatar
Marina Ganeva committed
1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772
    const std::vector<coord_t> ref = {0,
                                      0,
                                      0,
                                      3,
                                      -0.0943541043211f,
                                      -0.0943541043211f,
                                      2.51817307323f,
                                      -2.22473993186f, // 0
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.105491719946f,
                                      -0.105491719946f,
                                      2.18460421817f,
                                      -1.53897441518f, // 1
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.115542738924f,
                                      -0.115542738924f,
                                      1.88357866605f,
                                      -0.951095651067f, // 2
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.124658779393f,
                                      -0.124658779393f,
                                      1.61055549086f,
                                      -0.443322919741f, // 3
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.132964505153f,
                                      -0.132964505153f,
                                      1.36180104236f,
                                      -0.00173676689289f, // 4
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.140563360636f,
                                      -0.140563360636f,
                                      1.13421718522f,
                                      0.384685191257f, // 5
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.147541901385f,
                                      -0.147541901385f,
                                      0.925211602133f,
                                      0.724762824962f, // 6
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.153973105606f,
                                      -0.153973105606f,
                                      0.732598613797f,
                                      1.02562126597f, // 7
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.159918935922f,
                                      -0.159918935922f,
                                      0.55452245477f,
                                      1.29306717609f, // 8
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.165432342216f,
                                      -0.165432342216f,
                                      0.389397289126f,
                                      1.53187104043f, // 9
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.170558842805f,
                                      -0.170558842805f,
                                      0.235859854405f,
                                      1.7459813825f, // 10
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.175337784032f,
                                      -0.175337784032f,
                                      0.0927317372921f,
                                      1.93868903792f, // 11
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.179803352064f,
                                      -0.179803352064f,
                                      -0.0410109295185f,
                                      2.11275436309f, // 12
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.183985391967f,
                                      -0.183985391967f,
                                      -0.166261998436f,
                                      2.27050663742f, // 13
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.187910075568f,
                                      -0.187910075568f,
                                      -0.283805309267f,
                                      2.41392239469f, // 14
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.1916004497f,
                                      -0.1916004497f,
                                      -0.394331109003f,
                                      2.54468763779f, // 15
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.1950768891f,
                                      -0.1950768891f,
                                      -0.498449616001f,
                                      2.664247618f, // 16
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.198357472759f,
                                      -0.198357472759f,
                                      -0.596702291619f,
                                      2.77384694056f, // 17
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.20145829841f,
                                      -0.20145829841f,
                                      -0.689571258983f,
                                      2.87456208674f, // 18
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.204393746692f,
                                      -0.204393746692f,
                                      -0.777487214755f,
                                      2.96732794808f, // 19
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.207176704155f,
                                      -0.207176704155f,
                                      -0.86083610789f,
                                      3.05295960034f, // 20
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.209818752378f,
                                      -0.209818752378f,
                                      -0.939964803903f,
                                      3.13217026882f, // 21
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.212330329085f,
                                      -0.212330329085f,
                                      -1.01518590999f,
                                      3.20558622767f, // 22
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.214720865951f,
                                      -0.214720865951f,
                                      -1.08678190253f,
                                      3.27375921711f, // 23
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.216998906964f,
                                      -0.216998906964f,
                                      -1.15500867189f,
                                      3.33717683991f, // 24
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.219172210459f,
                                      -0.219172210459f,
                                      -1.2200985783f,
                                      3.39627130464f, // 25
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.221247837392f,
                                      -0.221247837392f,
                                      -1.28226309565f,
                                      3.45142680925f, // 26
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.223232227977f,
                                      -0.223232227977f,
                                      -1.34169510674f,
                                      3.5029858015f, // 27
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.225131268429f,
                                      -0.225131268429f,
                                      -1.39857090231f,
                                      3.55125430721f, // 28
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.226950349283f,
                                      -0.226950349283f,
                                      -1.45305192753f,
                                      3.59650648186f, // 29
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.228694416494f,
                                      -0.228694416494f,
                                      -1.50528631254f,
                                      3.63898851198f, // 30
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.230368016343f,
                                      -0.230368016343f,
                                      -1.55541021734f,
                                      3.67892197067f, // 31
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.231975335009f,
                                      -0.231975335009f,
                                      -1.60354901701f,
                                      3.71650671254f, // 32
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.233520233533f,
                                      -0.233520233533f,
                                      -1.64981834872f,
                                      3.75192337916f, // 33
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.23500627878f,
                                      -0.23500627878f,
                                      -1.69432503923f,
                                      3.78533557363f, // 34
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.236436770934f,
                                      -0.236436770934f,
                                      -1.73716792823f,
                                      3.81689175332f, // 35
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.237814767963f,
                                      -0.237814767963f,
                                      -1.77843860111f,
                                      3.84672688183f, // 36
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.239143107441f,
                                      -0.239143107441f,
                                      -1.81822204254f,
                                      3.87496387446f, // 37
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.240424426052f,
                                      -0.240424426052f,
                                      -1.85659722056f,
                                      3.90171486617f, // 38
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.24166117706f,
                                      -0.24166117706f,
                                      -1.89363760977f,
                                      3.92708232664f, // 39
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.242855645982f,
                                      -0.242855645982f,
                                      -1.92941166089f,
                                      3.95116004298f, // 40
                                      1,
                                      1,
                                      0,
                                      3,
                                      -0.244009964688f,
                                      -0.244009964688f,
                                      -1.9639832229f,
                                      3.97403398783f, // 41
                                      1,
                                      1,
                                      0,
                                      3,
                                      -0.245126124099f,
                                      -0.245126124099f,
                                      -1.99741192336f,
                                      3.99578308789f, // 42
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.246205985642f,
                                      -0.246205985642f,
                                      -2.0297535116f,
                                      4.01647990566f,
                                      2,
                                      2,
                                      0,
                                      3,
                                      -0.247251291615f,
                                      -0.247251291615f,
                                      -2.06106016902f,
                                      4.03619124547f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.248263674566f,
                                      -0.248263674566f,
                                      -2.09138078999f,
                                      4.05497869322f,
                                      2,
                                      2,
                                      0,
                                      3,
                                      -0.249244665798f,
                                      -0.249244665798f,
                                      -2.12076123667f,
                                      4.07289909803f,
                                      6,
                                      6,
                                      0,
                                      3,
                                      -0.250195703099f,
                                      -0.250195703099f,
                                      -2.14924457046f,
                                      4.09000500276f, // EPP
                                      3,
                                      3,
                                      0,
                                      3,
                                      -0.251118137774f,
                                      -0.251118137774f,
                                      -2.17687126264f,
                                      4.10634502964f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.252013241051f,
                                      -0.252013241051f,
                                      -2.20367938617f,
                                      4.12196422612f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.252882209926f,
                                      -0.252882209926f,
                                      -2.22970479075f,
                                      4.1369043757f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.253726172503f,
                                      -0.253726172503f,
                                      -2.25498126284f,
                                      4.15120427767f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.254546192879f,
                                      -0.254546192879f,
                                      -2.27954067188f,
                                      4.16489999925f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.255343275622f,
                                      -0.255343275622f,
                                      -2.30341310445f,
                                      4.17802510323f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.256118369875f,
                                      -0.256118369875f,
                                      -2.32662698715f,
                                      4.19061085383f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.256872373129f,
                                      -0.256872373129f,
                                      -2.34920919958f,
                                      4.20268640299f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.257606134684f,
                                      -0.257606134684f,
                                      -2.37118517811f,
                                      4.2142789594f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.258320458847f,
                                      -0.258320458847f,
                                      -2.39257901152f,
                                      4.22541394183f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.259016107869f,
                                      -0.259016107869f,
                                      -2.41341352902f,
                                      4.23611511866f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.259693804658f,
                                      -0.259693804658f,
                                      -2.43371038155f,
                                      4.24640473475f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.260354235286f,
                                      -0.260354235286f,
                                      -2.45349011683f,
                                      4.25630362721f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.260998051308f,
                                      -0.260998051308f,
                                      -2.4727722487f,
                                      4.2658313309f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.2616258719f,
                                      -0.2616258719f,
                                      -2.49157532139f,
                                      4.27500617494f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.262238285852f,
                                      -0.262238285852f,
                                      -2.50991696898f,
                                      4.2838453709f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.262835853406f,
                                      -0.262835853406f,
                                      -2.52781397058f,
                                      4.2923650936f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.263419107964f,
                                      -0.263419107964f,
                                      -2.54528230147f,
                                      4.30058055512f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.26398855768f,
                                      -0.26398855768f,
                                      -2.56233718075f,
                                      4.3085060728f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.264544686935f,
                                      -0.264544686935f,
                                      -2.57899311548f,
                                      4.31615513161f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.265087957711f,
                                      -0.265087957711f,
                                      -2.5952639419f,
                                      4.32354044159f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.265618810868f,
                                      -0.265618810868f,
                                      -2.61116286371f,
                                      4.33067399066f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.266137667344f,
                                      -0.266137667344f,
                                      -2.62670248786f,
                                      4.33756709332f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.266644929262f,
                                      -0.266644929262f,
                                      -2.64189485783f,
                                      4.34423043551f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.267140980972f,
                                      -0.267140980972f,
                                      -2.65675148482f,
                                      4.35067411606f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.267626190022f,
                                      -0.267626190022f,
                                      -2.67128337679f,
                                      4.35690768501f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.268100908065f,
                                      -0.268100908065f,
                                      -2.6855010657f,
                                      4.36294017896f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.268565471711f,
                                      -0.268565471711f,
                                      -2.69941463291f,
                                      4.3687801539f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.269020203322f,
                                      -0.269020203322f,
                                      -2.7130337331f,
                                      4.37443571549f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.269465411759f,
                                      -0.269465411759f,
                                      -2.72636761653f,
                                      4.37991454727f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.269901393077f,
                                      -0.269901393077f,
                                      -2.73942515004f,
                                      4.38522393675f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.27032843119f,
                                      -0.27032843119f,
                                      -2.75221483671f,
                                      4.39037079963f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.270746798477f,
                                      -0.270746798477f,
                                      -2.76474483431f,
                                      4.39536170235f,
                                      0,
                                      0,
                                      0,
                                      3,
                                      -0.271156756372f,
                                      -0.271156756372f,
                                      -2.77702297266f,
                                      4.40020288306f};
    double d(1.0e-06);
    sort_Events(events);
Marina Ganeva's avatar
Marina Ganeva committed
    for (auto i = 0; i < 82 * 8; i++) {
      TS_ASSERT_DELTA(events[i], ref[i], d);
    AnalysisDataService::Instance().remove(outWSName);

    // test the normalization workspace as well
    IMDEventWorkspace_sptr nws;
    TS_ASSERT_THROWS_NOTHING(
        nws = AnalysisDataService::Instance().retrieveWS<IMDEventWorkspace>(
            normWSName));
    TS_ASSERT(nws);
    // there are 7 histograms (the rest is outside of 2theta limits)
    TS_ASSERT_EQUALS(nws->getNPoints(), 574);

    AnalysisDataService::Instance().remove(normWSName);
  void test_TOFWSDataRotateEPP() {
    // test whether the calculation for inelastic data are correct

    std::string outWSName("LoadDNSSCDTest_OutputWS");
    std::string normWSName("LoadDNSSCDTest_OutputWS_norm");

    LoadDNSSCD alg;
    TS_ASSERT_THROWS_NOTHING(alg.initialize());
    TS_ASSERT(alg.isInitialized());
    TS_ASSERT_THROWS_NOTHING(alg.setPropertyValue("Filenames", "dnstof.d_dat"));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("OutputWorkspace", outWSName));
    TS_ASSERT_THROWS_NOTHING(
        alg.setPropertyValue("NormalizationWorkspace", normWSName));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("Normalization", "monitor"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("a", 3.55));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("b", 3.55));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("c", 24.778));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("alpha", 90.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("beta", 90.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("gamma", 120.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("OmegaOffset", -43.0));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("HKL1", "1,1,0"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("HKL2", "0,0,1"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("TwoThetaLimits", "25.0,60.0"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("DeltaEmin", "-3.0"));
    TS_ASSERT_THROWS_NOTHING(alg.setProperty("ElasticChannel", "64"));
    TS_ASSERT_THROWS_NOTHING(alg.execute(););
    TS_ASSERT(alg.isExecuted());
    // Retrieve the workspace from data service.
    IMDEventWorkspace_sptr iws;
    TS_ASSERT_THROWS_NOTHING(
        iws = AnalysisDataService::Instance().retrieveWS<IMDEventWorkspace>(
            outWSName));
    TS_ASSERT(iws);

    std::vector<API::IMDNode *> boxes(0, NULL);
    iws->getBoxes(boxes, 10000, false);
    TSM_ASSERT_EQUALS("Number of boxes", boxes.size(), 1);
    API::IMDNode *box = boxes[0];
    // there are 7 points (the rest is outside of 2theta limits)
    TS_ASSERT_EQUALS(box->getNPoints(), 574);
    std::vector<coord_t> events;
    size_t ncols;
    box->getEventsData(events, ncols);
    // 8 columns: I, err^2, run_num, det_id, h, k, l, dE
    TS_ASSERT_EQUALS(ncols, 8);
    // 8*574 = 4592
    TS_ASSERT_EQUALS(events.size(), 4592);
    // reference vector, EPP should be near dE=0
Marina Ganeva's avatar
Marina Ganeva committed
    const std::vector<coord_t> ref = {0,
                                      0,
                                      0,
                                      4,
                                      0.0354153287653f,
                                      0.0354153287653f,
                                      3.30960431031f,
                                      -2.22473993186f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      0.0119271789404f,
                                      0.0119271789404f,
                                      3.14246833446f,
                                      -1.53897441518f,
                                      1,
                                      1,
                                      0,
                                      4,
                                      -0.00926944407241f,
                                      -0.00926944407241f,
                                      2.99163830747f,
                                      -0.951095651067f,
                                      1,
                                      1,
                                      0,
                                      4,
                                      -0.0284942882003f,
                                      -0.0284942882003f,
                                      2.85483898067f,
                                      -0.443322919741f,
                                      1,
                                      1,
                                      0,
                                      4,
                                      -0.0460102572946f,
                                      -0.0460102572946f,
                                      2.73019959403f,
                                      -0.00173676689289f,
                                      3,
                                      3,
                                      0,
                                      4,
                                      -0.0620355056149f,
                                      -0.0620355056149f,
                                      2.61616781476f,
                                      0.384685191257f,
                                      1,
                                      1,
                                      0,
                                      4,
                                      -0.0767525703988f,
                                      -0.0767525703988f,
                                      2.51144475217f,
                                      0.724762824962f,
                                      1,
                                      1,
                                      0,
                                      4,
                                      -0.0903153555918f,
                                      -0.0903153555918f,
                                      2.41493526311f,
                                      1.02562126597f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.102854534355f,
                                      -0.102854534355f,
                                      2.32570950908f,
                                      1.29306717609f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.114481772845f,
                                      -0.114481772845f,
                                      2.24297290079f,
                                      1.53187104043f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.125293064774f,
                                      -0.125293064774f,
                                      2.16604237028f,
                                      1.7459813825f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.135371387759f,
                                      -0.135371387759f,
                                      2.09432746896f,
                                      1.93868903792f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.144788837105f,
                                      -0.144788837105f,
                                      2.02731518411f,
                                      2.11275436309f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.153608353159f,
                                      -0.153608353159f,
                                      1.96455764752f,
                                      2.27050663742f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.161885129764f,
                                      -0.161885129764f,
                                      1.90566211317f,
                                      2.41392239469f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.169667770453f,
                                      -0.169667770453f,
                                      1.85028273012f,
                                      2.54468763779f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.176999243565f,
                                      -0.176999243565f,
                                      1.7981137461f,
                                      2.664247618f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.183917675938f,
                                      -0.183917675938f,
                                      1.74888385976f,
                                      2.77384694056f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.190457016127f,
                                      -0.190457016127f,
                                      1.70235150145f,
                                      2.87456208674f,
                                      0,
                                      0,
                                      0,
                                      4,
                                      -0.196647591505f,