testAnalyticWindows.cc

Go to the documentation of this file.

#include <iostream>

#include <tst/Verify.h>

#include <revt/Channel.h>
#include <revt/Station.h>
#include <utl/HannWindow.h>
#include <utl/AugerException.h>
#include <utl/AugerUnits.h>
#include <utl/PhysicalConstants.h>
#include <cmath>

#include <cppunit/extensions/HelperMacros.h>

using namespace std;
using namespace utl;
using namespace tst;
using namespace revt;


class testAnalyticWindows : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testAnalyticWindows);
  CPPUNIT_TEST(testHannWindow);
  CPPUNIT_TEST_SUITE_END();

private:

public:
  void
  setUp()
  {
  }

  void
  tearDown()
  {
  }

  void
  testHannWindow()
  {
    long traceLength = 1000;
    double windowWidth = 0.1;

    HannWindow ourWindow(windowWidth, traceLength);

    // check edges of the window
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(long(windowWidth*traceLength)+1),1.0));
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(traceLength-long(windowWidth*traceLength+1)),1.0));
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(0),0.0));
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(traceLength-1),0.0));

    double testLocation = 0.3; // relative location inside each side of the Hann window to test for

    // check weights inside the window
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(testLocation*long(windowWidth*traceLength)),0.5*(1-std::cos(2*kPi*testLocation*long(windowWidth*traceLength)/(long(2*windowWidth*traceLength)-1)))));
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(traceLength-1-(testLocation*long(windowWidth*traceLength))),0.5*(1-std::cos(2*kPi*testLocation*long(windowWidth*traceLength)/(long(2*windowWidth*traceLength)-1)))));
  
    // redefine trace length of Hann window
    traceLength = 2500;
    ourWindow.SetTraceLength(traceLength);
    
    // check redefinition of trace length
    CPPUNIT_ASSERT(Verify<Equal>(ourWindow.GetTraceLength(),traceLength));
        
    testLocation = 0.4; // relative location inside each side of the Hann window to test for
  
    // check weights inside the window with the new trace length
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(testLocation*long(windowWidth*traceLength)),0.5*(1-std::cos(2*kPi*testLocation*long(windowWidth*traceLength)/(long(2*windowWidth*traceLength)-1)))));
    CPPUNIT_ASSERT(Verify<CloseTo>(ourWindow.GetWeightAtBin(traceLength-1-(testLocation*long(windowWidth*traceLength))),0.5*(1-std::cos(2*kPi*testLocation*long(windowWidth*traceLength)/(long(2*windowWidth*traceLength)-1)))));
  }

};

CPPUNIT_TEST_SUITE_REGISTRATION(testAnalyticWindows);

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: