Refs #22420 integration working for satellites

Framework/DataObjects/src/PeaksWorkspace.cpp

@@ -527,6 +527,7 @@ Peak *PeaksWorkspace::createPeakHKL(const V3D &HKL) const {
 
   // We need to set HKL separately to keep things consistent.
   peak->setHKL(HKL[0], HKL[1], HKL[2]);
+  peak->setIntHKL(peak->getHKL());
 
   // Set the goniometer
   peak->setGoniometerMatrix(goniometer.getR());

Framework/MDAlgorithms/CMakeLists.txt

@@ -318,6 +318,7 @@ set ( TEST_FILES
 	Integrate3DEventsTest.h
 	IntegrateEllipsoidsTest.h
 	IntegrateEllipsoidsTwoStepTest.h
+	IntegrateEllipsoidsWithSatellitesTest.h
 	IntegrateFluxTest.h
 	IntegrateMDHistoWorkspaceTest.h
 	IntegratePeaksCWSDTest.h

Framework/MDAlgorithms/inc/MantidMDAlgorithms/IntegrateEllipsoidsWithSatellites.h

+// 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
+// SPDX - License - Identifier: GPL - 3.0 +
 #ifndef MANTID_MDALGORITHMS_INTEGRATE_ELLIPSOIDS_with_Satellites_H_
 #define MANTID_MDALGORITHMS_INTEGRATE_ELLIPSOIDS_with_Satellites_H_
 

Framework/MDAlgorithms/src/Integrate3DEvents.cpp

@@ -752,7 +752,7 @@ int64_t Integrate3DEvents::getHklKey(int h, int k, int l) {
   int64_t key(0);
 
   if (h != 0 || k != 0 || l != 0)
-    key = 1000000000000 * h + 1000000 * k + l;
+    key = 1000000000000 * h + 100000000 * k + 10000 * l;
 
   return key;
 }
@@ -801,10 +801,6 @@ int64_t Integrate3DEvents::getHklMnpKey2(V3D const &hkl) {
   V3D modvec2 = V3D(m_ModHKL[0][1], m_ModHKL[1][1], m_ModHKL[2][1]);
   V3D modvec3 = V3D(m_ModHKL[0][2], m_ModHKL[1][2], m_ModHKL[2][2]);
   if (Geometry::IndexingUtils::ValidIndex(hkl, m_radius))
-
-  //            if (fabs(hkl[0]-boost::math::iround<double>(hkl[0]))<m_radius &&
-  //            fabs(hkl[1]-boost::math::iround<double>(hkl[1]))<m_radius &&
-  //            fabs(hkl[2]-boost::math::iround<double>(hkl[2]))<m_radius)
   {
     int h = boost::math::iround<double>(hkl[0]);
     int k = boost::math::iround<double>(hkl[1]);
@@ -821,12 +817,6 @@ int64_t Integrate3DEvents::getHklMnpKey2(V3D const &hkl) {
         hkl1[1] -= order * modvec1[1];
         hkl1[2] -= order * modvec1[2];
         if (Geometry::IndexingUtils::ValidIndex(hkl1, s_radius))
-        //                        if
-        //                        (fabs(hkl1[0]-boost::math::iround<double>(hkl1[0]))<m_radius
-        //                        &&
-        //                        fabs(hkl1[1]-boost::math::iround<double>(hkl1[1]))<m_radius
-        //                        &&
-        //                        fabs(hkl1[2]-boost::math::iround<double>(hkl1[2]))<m_radius)
         {
           int h = boost::math::iround<double>(hkl1[0]);
           int k = boost::math::iround<double>(hkl1[1]);

Framework/MDAlgorithms/src/IntegrateEllipsoidsWithSatellites.cpp

+// 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
+// SPDX - License - Identifier: GPL - 3.0 +
 #include "MantidMDAlgorithms/IntegrateEllipsoidsWithSatellites.h"
 #include "MantidAPI/AnalysisDataService.h"
 #include "MantidAPI/FileProperty.h"
@@ -281,7 +287,7 @@ void IntegrateEllipsoidsWithSatellites::init() {
                   "major radius for second pass when SpecifySize is false.");
 
   declareProperty(
-      "IntegrateIfOnEdge", false,
+      "IntegrateIfOnEdge", true,
       "Set to false to not integrate if peak radius is off edge of detector."
       "Background will be scaled if background radius is off edge.");
 
@@ -387,17 +393,11 @@ void IntegrateEllipsoidsWithSatellites::exec() {
     if (mnp[2] != 0 && ModDim == 2)
       ModDim = 3;
 
-    if (Geometry::IndexingUtils::ValidIndex(hkl, 1.0) ||
-        Geometry::IndexingUtils::ValidIndex(
-            mnp, 1.0)) // use tolerance == 1 to
+    if (Geometry::IndexingUtils::ValidIndex(hkl, 1.0)) // use tolerance == 1 to
                        // just check for (0,0,0,0,0,0)
     {
       peak_q_list.emplace_back(peaks[i].getQLabFrame());
       qList.emplace_back(1., V3D(peaks[i].getQLabFrame()));
-      //                    V3D
-      //                    miller_ind(static_cast<double>(boost::math::iround<double>(hkl[0])),
-      //                                   static_cast<double>(boost::math::iround<double>(hkl[1])),
-      //                                   static_cast<double>(boost::math::iround<double>(hkl[2])));
       hkl_vectors.push_back(hkl);
       mnp_vectors.push_back(mnp);
       indexed_count++;
@@ -416,14 +416,11 @@ void IntegrateEllipsoidsWithSatellites::exec() {
   Geometry::IndexingUtils::Optimize_6dUB(UB, modUB, hkl_vectors, mnp_vectors,
                                          ModDim, peak_q_list);
 
-  OrientedLattice lattice;
-  //            lattice = peak_ws->mutableSample().getOrientedLattice();
-  //            UB = lattice.getUB();
+  OrientedLattice lattice = peak_ws->mutableSample().getOrientedLattice();
   lattice.setUB(UB);
   lattice.setModUB(modUB);
   modHKL = lattice.getModHKL();
 
-  lattice = peak_ws->mutableSample().getOrientedLattice();
   int maxOrder = lattice.getMaxOrder();
   bool CT = lattice.getCrossTerm();
 

Framework/MDAlgorithms/test/IntegrateEllipsoidsWithSatellitesTest.h

+// 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
+// SPDX - License - Identifier: GPL - 3.0 +
+#include "MantidAPI/AlgorithmManager.h"
+#include "MantidAPI/FrameworkManager.h"
+#include "MantidAPI/NumericAxis.h"
+#include "MantidAPI/Sample.h"
+#include "MantidDataObjects/EventWorkspace.h"
+#include "MantidDataObjects/PeakShapeEllipsoid.h"
+#include "MantidDataObjects/PeaksWorkspace.h"
+#include "MantidDataObjects/WorkspaceSingleValue.h"
+#include "MantidGeometry/Crystal/OrientedLattice.h"
+#include "MantidMDAlgorithms/IntegrateEllipsoidsWithSatellites.h"
+#include "MantidTestHelpers/ComponentCreationHelper.h"
+#include <boost/make_shared.hpp>
+#include <boost/tuple/tuple.hpp>
+#include <cxxtest/TestSuite.h>
+
+using namespace Mantid;
+using namespace Mantid::MDAlgorithms;
+using namespace Mantid::Kernel;
+using namespace Mantid::Geometry;
+using namespace Mantid::DataObjects;
+using Mantid::Types::Event::TofEvent;
+
+namespace {
+// Add A Fake 'Peak' to both the event data and to the peaks workspace
+void addFakeEllipsoid(const V3D &peakHKL, const int &totalNPixels,
+                      const int &nEvents, const double tofGap,
+                      EventWorkspace_sptr &eventWS,
+                      PeaksWorkspace_sptr &peaksWS) {
+  // Create the peak and add it to the peaks ws
+  Peak *peak = peaksWS->createPeakHKL(peakHKL);
+  peaksWS->addPeak(*peak);
+  const int detectorId = peak->getDetectorID();
+  const double tofExact = peak->getTOF();
+  delete peak;
+
+  EventList &el = eventWS->getSpectrum(detectorId - totalNPixels);
+
+  // Add more events to the event list corresponding to the peak centre
+  double start = tofExact - (double(nEvents) / 2 * tofGap);
+  for (int i = 0; i < nEvents; ++i) {
+    const double tof = start + (i * tofGap);
+    el.addEventQuickly(TofEvent(tof));
+  }
+}
+
+// Create diffraction data for test schenarios
+boost::tuple<EventWorkspace_sptr, PeaksWorkspace_sptr>
+createDiffractionData(const int nPixels = 100, const int nEventsPerPeak = 20,
+                      const double tofGapBetweenEvents = 10) {
+  Mantid::Geometry::Instrument_sptr inst =
+      ComponentCreationHelper::createTestInstrumentRectangular(
+          1 /*num_banks*/, nPixels /*pixels in each direction yields n by n*/,
+          0.01, 1.0);
+
+  // Create a peaks workspace
+  auto peaksWS = boost::make_shared<PeaksWorkspace>();
+  // Set the instrument to be the fake rectangular bank above.
+  peaksWS->setInstrument(inst);
+  // Set the oriented lattice for a cubic crystal
+  OrientedLattice ol(6, 6, 6, 90, 90, 90);
+  ol.setUFromVectors(V3D(6, 0, 0), V3D(0, 6, 0));
+  peaksWS->mutableSample().setOrientedLattice(&ol);
+
+  // Make an event workspace and add fake peak data
+  auto eventWS = boost::make_shared<EventWorkspace>();
+  eventWS->setInstrument(inst);
+  eventWS->initialize(nPixels * nPixels /*n spectra*/, 3 /* x-size */,
+                      3 /* y-size */);
+  eventWS->getAxis(0)->setUnit("TOF");
+  // Give the spectra-detector mapping for all event lists
+  const int nPixelsTotal = nPixels * nPixels;
+  for (int i = 0; i < nPixelsTotal; ++i) {
+    EventList &el = eventWS->getSpectrum(i);
+    el.setDetectorID(i + nPixelsTotal);
+  }
+
+  // Add some peaks which should correspond to real reflections (could
+  // calculate these). Same function also adds a fake ellipsoid
+  addFakeEllipsoid(V3D(1, -5, -3), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+  addFakeEllipsoid(V3D(1, -4, -4), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+  addFakeEllipsoid(V3D(1, -3, -5), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+  addFakeEllipsoid(V3D(1, -4, -2), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+  addFakeEllipsoid(V3D(1, -4, 0), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+  addFakeEllipsoid(V3D(2, -3, -4), nPixelsTotal, nEventsPerPeak,
+                   tofGapBetweenEvents, eventWS, peaksWS);
+
+  // Return test data.
+  return boost::tuple<EventWorkspace_sptr, PeaksWorkspace_sptr>(eventWS,
+                                                                peaksWS);
+}
+} // namespace
+
+class IntegrateEllipsoidsWithSatellitesTest : public CxxTest::TestSuite {
+
+private:
+  Mantid::DataObjects::EventWorkspace_sptr m_eventWS;
+  Mantid::DataObjects::PeaksWorkspace_sptr m_peaksWS;
+  Mantid::API::MatrixWorkspace_sptr m_histoWS;
+
+  // Check that n-peaks from the workspace are integrated as we expect
+  void do_test_n_peaks(PeaksWorkspace_sptr &integratedPeaksWS,
+                       const int nPeaks) {
+    auto instrument = integratedPeaksWS->getInstrument();
+    const V3D samplePos = instrument->getComponentByName("sample")->getPos();
+    const V3D sourcePos = instrument->getComponentByName("source")->getPos();
+    V3D beamDir = samplePos - sourcePos;
+    beamDir.normalize();
+
+    // Just test the first few peaks
+    for (int i = 0; i < nPeaks; ++i) {
+      const Peak &peak = integratedPeaksWS->getPeak(i);
+      const PeakShape &peakShape = peak.getPeakShape();
+
+      TSM_ASSERT_RELATION("Peak should be integrated", std::greater<double>,
+                          peak.getIntensity(), 0);
+
+      std::stringstream stream;
+      stream << "Wrong shape name for peak " << i;
+      TSM_ASSERT_EQUALS(stream.str(), PeakShapeEllipsoid::ellipsoidShapeName(),
+                        peakShape.shapeName());
+
+      // Calculate the q direction based on geometry
+      const V3D detPos = peak.getDetectorPosition();
+      V3D detDir = detPos - samplePos;
+      detDir.normalize();
+      V3D qDir = detDir - beamDir;
+      qDir.normalize();
+
+      // Get the q-direction off the ellipsoid
+      PeakShapeEllipsoid const *const ellipsoid =
+          dynamic_cast<const PeakShapeEllipsoid *>(&peakShape);
+      auto dirs = ellipsoid->directions();
+
+      /* We have set the fake ellipsoids up to be lines along a single detectors
+       * TOF (see setup).
+       * We therefore expect the principle axis of the ellipsoid to be the same
+       * as the q-dir!
+       */
+      TS_ASSERT_EQUALS(qDir, dirs[0]);
+    }
+  }
+
+public:
+  static void destroySuite(IntegrateEllipsoidsWithSatellitesTest *suite) { delete suite; }
+
+  // This pair of boilerplate methods prevent the suite being created statically
+  // This means the constructor isn't called when running other tests
+  static IntegrateEllipsoidsWithSatellitesTest *createSuite() {
+    return new IntegrateEllipsoidsWithSatellitesTest();
+  }
+
+  IntegrateEllipsoidsWithSatellitesTest() {
+
+    // Need to get and run algorithms from elsewhere in the framework.
+    Mantid::API::FrameworkManager::Instance();
+
+    auto data = createDiffractionData();
+
+    m_eventWS = data.get<0>();
+    m_peaksWS = data.get<1>();
+
+    /*
+     Simply rebin the event workspace to a histo workspace to create the input
+     we need.
+    */
+    auto rebinAlg =
+        Mantid::API::AlgorithmManager::Instance().createUnmanaged("Rebin");
+    rebinAlg->setChild(true);
+    rebinAlg->initialize();
+    rebinAlg->setProperty("InputWorkspace", m_eventWS);
+    auto params = std::vector<double>();
+    params.push_back(950);
+    params.push_back(10);
+    params.push_back(2500);
+    rebinAlg->setProperty("Params", params);
+    rebinAlg->setProperty("PreserveEvents", false); // Make a histo workspace
+    rebinAlg->setPropertyValue("OutputWorkspace", "dummy");
+    rebinAlg->execute();
+
+    m_histoWS = rebinAlg->getProperty("OutputWorkspace");
+  }
+
+  void test_init() {
+    Mantid::MDAlgorithms::IntegrateEllipsoidsWithSatellites alg;
+    TS_ASSERT_THROWS_NOTHING(alg.initialize());
+  }
+
+  void test_ws_has_instrument() {
+    auto inputWorkspaceNoInstrument = boost::make_shared<EventWorkspace>();
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    TS_ASSERT_THROWS(
+        alg.setProperty("InputWorkspace", inputWorkspaceNoInstrument),
+        std::invalid_argument &);
+  }
+
+  void test_execution_events() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_eventWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+
+    do_test_n_peaks(integratedPeaksWS, 3 /*check first 3 peaks*/);
+  }
+
+  void test_execution_histograms() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_histoWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    alg.setProperty("RegionRadius", 0.35);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundInnerSize", 0.23));
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+
+    do_test_n_peaks(integratedPeaksWS, 3 /*check first 3 peaks*/);
+
+    const auto &peak1 = integratedPeaksWS->getPeak(0);
+    const auto &peak2 = integratedPeaksWS->getPeak(1);
+    const auto &peak3 = integratedPeaksWS->getPeak(2);
+
+    TS_ASSERT_DELTA(peak1.getIntensity(), 1., 1e-6);
+    TS_ASSERT_DELTA(peak2.getIntensity(), 1., 1e-6);
+    TS_ASSERT_DELTA(peak3.getIntensity(), 1., 1e-6);
+  }
+
+  void test_execution_histograms_distribution_data() {
+    using namespace Mantid::API;
+    const auto &algManager = AlgorithmManager::Instance();
+
+    auto cloneWorkspace = algManager.createUnmanaged("CloneWorkspace");
+    cloneWorkspace->setChild(true);
+    cloneWorkspace->initialize();
+    cloneWorkspace->setProperty("InputWorkspace", m_histoWS);
+    cloneWorkspace->setPropertyValue("OutputWorkspace", "dist_workspace");
+    cloneWorkspace->execute();
+    Workspace_sptr temp = cloneWorkspace->getProperty("OutputWorkspace");
+    auto distWS = boost::dynamic_pointer_cast<MatrixWorkspace>(temp);
+
+    auto convertToDist = algManager.createUnmanaged("ConvertToDistribution");
+    convertToDist->setChild(true);
+    convertToDist->initialize();
+    convertToDist->setProperty("Workspace", distWS);
+    convertToDist->execute();
+    distWS = convertToDist->getProperty("Workspace");
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", distWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    alg.setProperty("RegionRadius", 0.35);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundInnerSize", 0.23));
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+
+    do_test_n_peaks(integratedPeaksWS, 3 /*check first 3 peaks*/);
+
+    const auto &peak1 = integratedPeaksWS->getPeak(0);
+    const auto &peak2 = integratedPeaksWS->getPeak(1);
+    const auto &peak3 = integratedPeaksWS->getPeak(2);
+
+    const double binWidth{10.};
+    TS_ASSERT_DELTA(peak1.getIntensity(), 1. / binWidth, 1e-6);
+    TS_ASSERT_DELTA(peak2.getIntensity(), 1. / binWidth, 1e-6);
+    TS_ASSERT_DELTA(peak3.getIntensity(), 1. / binWidth, 1e-6);
+  }
+
+  void test_execution_events_adaptive() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_eventWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("SpecifySize", true));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("PeakSize", 0.20));
+    alg.setProperty("RegionRadius", 0.35);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundInnerSize", 0.23));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundOuterSize", 0.26));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("AdaptiveQMultiplier", 0.01));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("AdaptiveQBackground", true));
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+
+    TSM_ASSERT_DELTA("Wrong intensity for peak 0",
+                     integratedPeaksWS->getPeak(0).getIntensity(), 2, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 1",
+                     integratedPeaksWS->getPeak(1).getIntensity(), 0.8, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 2",
+                     integratedPeaksWS->getPeak(2).getIntensity(), 2, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 3",
+                     integratedPeaksWS->getPeak(3).getIntensity(), 7, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 4",
+                     integratedPeaksWS->getPeak(4).getIntensity(), 0, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 5",
+                     integratedPeaksWS->getPeak(5).getIntensity(), 5.83, 0.01);
+  }
+
+  void test_execution_histograms_adaptive() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_histoWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("SpecifySize", true));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("PeakSize", 0.20));
+    alg.setProperty("RegionRadius", 0.35);
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundInnerSize", 0.23));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("BackgroundOuterSize", 0.26));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("AdaptiveQMultiplier", 0.01));
+    TS_ASSERT_THROWS_NOTHING(alg.setProperty("AdaptiveQBackground", true));
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+    TSM_ASSERT_DELTA("Wrong intensity for peak 0",
+                     integratedPeaksWS->getPeak(0).getIntensity(), 3, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 1",
+                     integratedPeaksWS->getPeak(1).getIntensity(), 3, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 2",
+                     integratedPeaksWS->getPeak(2).getIntensity(), 3, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 3",
+                     integratedPeaksWS->getPeak(3).getIntensity(), 10, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 4",
+                     integratedPeaksWS->getPeak(4).getIntensity(), 0, 0.01);
+    TSM_ASSERT_DELTA("Wrong intensity for peak 5",
+                     integratedPeaksWS->getPeak(5).getIntensity(), 9.94, 0.01);
+  }
+};
+
+class IntegrateEllipsoidsWithSatellitesTestPerformance : public CxxTest::TestSuite {
+
+private:
+  Mantid::API::MatrixWorkspace_sptr m_eventWS;
+  Mantid::DataObjects::PeaksWorkspace_sptr m_peaksWS;
+  Mantid::API::MatrixWorkspace_sptr m_histoWS;
+
+public:
+  static void destroySuite(IntegrateEllipsoidsWithSatellitesTestPerformance *suite) {
+    delete suite;
+  }
+
+  // This pair of boilerplate methods prevent the suite being created statically
+  // This means the constructor isn't called when running other tests
+  static IntegrateEllipsoidsWithSatellitesTestPerformance *createSuite() {
+    return new IntegrateEllipsoidsWithSatellitesTestPerformance();
+  }
+
+  IntegrateEllipsoidsWithSatellitesTestPerformance() {
+    // Need to get and run algorithms from elsewhere in the framework.
+    Mantid::API::FrameworkManager::Instance();
+
+    auto data = createDiffractionData(200 /*sqrt total pixels*/,
+                                      60 /*events per peak*/, 2 /*tof gap*/);
+
+    m_eventWS = data.get<0>();
+    m_peaksWS = data.get<1>();
+
+    /*
+     Simply rebin the event workspace to a histo workspace to create the input
+     we need.
+    */
+    auto rebinAlg =
+        Mantid::API::AlgorithmManager::Instance().createUnmanaged("Rebin");
+    rebinAlg->setChild(true);
+    rebinAlg->initialize();
+    rebinAlg->setProperty("InputWorkspace", m_eventWS);
+    auto params = std::vector<double>();
+    params.push_back(950);
+    params.push_back(5);
+    params.push_back(2500);
+    rebinAlg->setProperty("Params", params);
+    rebinAlg->setProperty("PreserveEvents", false); // Make a histo workspace
+    rebinAlg->setPropertyValue("OutputWorkspace", "dummy");
+    rebinAlg->execute();
+
+    m_histoWS = rebinAlg->getProperty("OutputWorkspace");
+  }
+
+  void test_execution_events() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_eventWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+  }
+
+  void test_execution_histograms() {
+
+    IntegrateEllipsoidsWithSatellites alg;
+    alg.setChild(true);
+    alg.setRethrows(true);
+    alg.initialize();
+    alg.setProperty("InputWorkspace", m_histoWS);
+    alg.setProperty("PeaksWorkspace", m_peaksWS);
+    alg.setPropertyValue("OutputWorkspace", "dummy");
+    alg.execute();
+    PeaksWorkspace_sptr integratedPeaksWS = alg.getProperty("OutputWorkspace");
+
+    TSM_ASSERT_EQUALS("Wrong number of peaks in output workspace",
+                      integratedPeaksWS->getNumberPeaks(),
+                      m_peaksWS->getNumberPeaks());
+  }
+};

docs/source/algorithms/IntegrateEllipsoidsWithSatellites-v1.rst

+.. algorithm::
+
+.. summary::
+
+.. relatedalgorithms::