DivideMDTest.h

#ifndef MANTID_MDALGORITHMS_DIVIDEMDTEST_H_
#define MANTID_MDALGORITHMS_DIVIDEMDTEST_H_

#include <cxxtest/TestSuite.h>
#include "MantidKernel/Timer.h"
#include "MantidKernel/System.h"
#include <iostream>
#include <iomanip>

#include "MantidMDAlgorithms/DivideMD.h"
#include "MantidTestHelpers/BinaryOperationMDTestHelper.h"
#include "MantidDataObjects/MDHistoWorkspace.h"

using namespace Mantid;
using namespace Mantid::MDAlgorithms;
using namespace Mantid::API;
using Mantid::DataObjects::MDHistoWorkspace_sptr;

/** Note: More detailed tests for the underlying
 * operations are in BinaryOperationMDTest and
 * MDHistoWorkspaceTest.
 *
 */
class DivideMDTest : public CxxTest::TestSuite
{
public:

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

  void test_histo_histo()
  {
    MDHistoWorkspace_sptr out;
    out = BinaryOperationMDTestHelper::doTest("DivideMD", "histo_A", "histo_B", "out");
    TS_ASSERT_DELTA( out->getSignalAt(0), 2./3., 1e-5);
  }

  void test_histo_scalar()
  {
    MDHistoWorkspace_sptr out;
    out = BinaryOperationMDTestHelper::doTest("DivideMD", "histo_A", "scalar", "out");
    TS_ASSERT_DELTA( out->getSignalAt(0), 2./3., 1e-5);
    BinaryOperationMDTestHelper::doTest("DivideMD", "scalar", "histo_A", "out", false /*fails*/);
  }

  void test_event_event_or_histo_fails()
  {
    BinaryOperationMDTestHelper::doTest("DivideMD", "event_A", "histo_A", "out", false /*fails*/);
    BinaryOperationMDTestHelper::doTest("DivideMD", "histo_A", "event_A", "out", false /*fails*/);
    BinaryOperationMDTestHelper::doTest("DivideMD", "event_A", "event_A", "out", false /*fails*/);
    BinaryOperationMDTestHelper::doTest("DivideMD", "scalar", "event_A", "out", false /*fails*/);
  }

  /** Get a MDEventWorkspace and check that all events have the given signal/error */
  void checkMDEWSignal(std::string wsName, signal_t expectedSignal, signal_t expectedError)
  {
    IMDEventWorkspace_sptr ws = AnalysisDataService::Instance().retrieveWS<IMDEventWorkspace>(wsName);
    TS_ASSERT(ws); if (!ws) return;
    IMDIterator* it = ws->createIterator(NULL);
    do
    {
      TS_ASSERT_EQUALS( it->getNumEvents(), 1);
      TS_ASSERT_DELTA( it->getInnerSignal(0), expectedSignal, 1e-5);
      TS_ASSERT_DELTA( it->getInnerError(0), expectedError, 1e-5);
    }
    while (it->next());
  }

  /// Divide events by a scalar
  void test_event_scalar()
  {
    BinaryOperationMDTestHelper::doTest("DivideMD", "event_A", "scalar", "out");
    checkMDEWSignal("out", 1./3., sqrt(12./81.));
  }


};


#endif /* MANTID_MDALGORITHMS_DIVIDEMDTEST_H_ */