testRSQLManager.cc

Go to the documentation of this file.

// Test giving out infos on Stations

#include <iostream>
#include <sstream>
#include <string>
#include <vector>
#include <det/Detector.h>
#include <rdet/RDetector.h>
#include <rdet/Station.h>
#include <rdet/Channel.h>
#include <sdet/SDetector.h>
#include <sdet/Station.h>
#include <fwk/CentralConfig.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <cppunit/extensions/HelperMacros.h>
#include <tst/Verify.h>
#include <utl/Point.h>
#include <utl/UTMPoint.h>
#include <utl/ReferenceEllipsoid.h>
#include <fwk/CoordinateSystemRegistry.h>
#include <fwk/LocalCoordinateSystem.h>
#include <TROOT.h>
#include <TH2F.h>
#include <TFile.h>
#include <TCanvas.h>
#include <TPaveText.h>
#include <TAttText.h>

using namespace det;
using namespace std;
using namespace utl;
using namespace tst;
using namespace rdet;


#define ASSERT_CLOSE(x, y, eps) CPPUNIT_ASSERT(Verify<CloseTo>(x, y, eps))
#define ASSERT_EQUAL(x, y) CPPUNIT_ASSERT(Verify<Equal>(x, y))


class RSQLManagerTest : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(RSQLManagerTest);
  CPPUNIT_TEST(testStationListManager);
  CPPUNIT_TEST(testBeaconManager);
  CPPUNIT_TEST_SUITE_END();

public:
  void setUp()
  {
    ErrorLogger::GetInstance().SetVerbosity(Verbosity::eVerbose);

    // CentralConfig must be created first.
    //
    /*CentralConfig* fCentralConfig =*/
    fwk::CentralConfig::GetInstance (BOOTSTRAPFILE);

    // Create a detector
    //
    /*Detector& theDetector =*/
    Detector::GetInstance();
  }

  void testStationListManager()
  {
    Detector& detector = Detector::GetInstance();
    TimeStamp detTime = UTCDateTime(2010, 10, 1, 16, 45, 12).GetTimeStamp();
    detector.Update(detTime);

    // try to get the name of the RDetector.Station id=1
    const RDetector& RDetector = detector.GetRDetector();
    const Station& RStation = RDetector.GetStation(101);
    cout << "Name of rdet::Station with Id 101: " << RStation.GetName() << endl;

    utl::Point pos1 = RDetector.GetStation(101).GetPosition();
    utl::Point pos2 = RDetector.GetStation(110).GetPosition();
    utl::Point pos3 = RDetector.GetStation(112).GetPosition();

    const UTMPoint utmpos1(pos1, ReferenceEllipsoid::GetWGS84());
    const UTMPoint utmpos2(pos2, ReferenceEllipsoid::GetWGS84());
    const UTMPoint utmpos3(pos3, ReferenceEllipsoid::GetWGS84());

    ASSERT_CLOSE(utmpos1.GetEasting(), 450885.827, 1e-8);
    ASSERT_CLOSE(utmpos1.GetNorthing(), 6114721.732, 1e-8);
    ASSERT_CLOSE(utmpos1.GetHeight(), 1566.993, 1e-6);

    ASSERT_CLOSE(utmpos2.GetEasting(), 450755.970, 1e-8);
    ASSERT_CLOSE(utmpos2.GetNorthing(), 6114509.317, 1e-8);
    ASSERT_CLOSE(utmpos2.GetHeight(), 1567.973, 1e-6);

    ASSERT_CLOSE(utmpos3.GetEasting(), 451251.865, 1e-8);
    ASSERT_CLOSE(utmpos3.GetNorthing(), 6114503.188, 1e-8);
    ASSERT_CLOSE(utmpos3.GetHeight(), 1562.512, 1e-6);

    ASSERT_CLOSE(RDetector.GetStation(101).GetChannel(1).GetSamplingFrequency(), 0.2, 1e-4);
    ASSERT_CLOSE(RDetector.GetStation(110).GetChannel(1).GetSamplingFrequency(), 0.2, 1e-4);

    cout << "update detector" << endl;
    detTime = UTCDateTime(2011, 5, 1, 16, 45, 12).GetTimeStamp();
    detector.Update(detTime);
    pos1 = RDetector.GetStation(101).GetPosition();

    cout << "update detector" << endl;
    detTime = UTCDateTime(2012, 5, 1, 16, 45, 12).GetTimeStamp();
    detector.Update(detTime);

    pos1 = RDetector.GetStation(101).GetPosition();
    pos2 = RDetector.GetStation(110).GetPosition();
    pos3 = RDetector.GetStation(112).GetPosition();

    const UTMPoint utmpos4(pos1, ReferenceEllipsoid::GetWGS84());
    const UTMPoint utmpos5(pos2, ReferenceEllipsoid::GetWGS84());
    const UTMPoint utmpos6(pos3, ReferenceEllipsoid::GetWGS84());

    ASSERT_CLOSE(utmpos4.GetEasting(), 450885.827, 1e-8);
    ASSERT_CLOSE(utmpos4.GetNorthing(), 6114721.732, 1e-8);
    ASSERT_CLOSE(utmpos4.GetHeight(), 1566.993, 1e-6);

    ASSERT_CLOSE(utmpos5.GetEasting(), 450755.970, 1e-8);
    ASSERT_CLOSE(utmpos5.GetNorthing(), 6114509.317, 1e-8);
    ASSERT_CLOSE(utmpos5.GetHeight(), 1567.973, 1e-6);

    ASSERT_CLOSE(utmpos6.GetEasting(), 451251.865, 1e-8);
    ASSERT_CLOSE(utmpos6.GetNorthing(), 6114503.188, 1e-8);
    ASSERT_CLOSE(utmpos6.GetHeight(), 1562.512, 1e-6);

    ASSERT_CLOSE(RDetector.GetStation(101).GetChannel(1).GetSamplingFrequency(), 0.18, 1e-4);
    ASSERT_CLOSE(RDetector.GetStation(110).GetChannel(1).GetSamplingFrequency(), 0.18, 1e-4);
  }

  void testBeaconManager()
  {
    Detector& detector = Detector::GetInstance();
    const TimeStamp detTime = UTCDateTime(2013, 11, 1, 16, 45, 12).GetTimeStamp();
    detector.Update(detTime);
    const RDetector& rDetector = detector.GetRDetector();
    const std::vector<double>& beaconFrequencies = rDetector.GetBeaconFrequencies();
    std::cout << "RDetector::GetBeaconFrequencies() returns frequencies (in GHz):" << std::endl;
    for (const auto& f : beaconFrequencies)
      std::cout << f << std::endl;
    // test Rdetector::GetBeaconReferencePhase
    const double beaconRefPhase = rDetector.GetBeaconReferencePhase(101, 0.061523);
    ASSERT_CLOSE(beaconRefPhase, 0.261232, 1e-5);
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(RSQLManagerTest);