testSDenseDetectorBasic.cc

Go to the documentation of this file.

#include <map>
#include <sstream>
#include <iostream>
#include <string>
#include <vector>

#include <det/Detector.h>
#include <sdet/SDetector.h>
#include <fwk/CentralConfig.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/TabulatedFunction.h>
#include <utl/Reader.h>
#include <utl/CoordinateSystem.h>
#include <utl/Point.h>
#include <utl/TimeStamp.h>
#include <utl/UTCDateTime.h>
#include <utl/AugerException.h>
#include <utl/ErrorLogger.h>

#include <fwk/CoordinateSystemRegistry.h>
#include <fwk/RunController.h>

#include <evt/Event.h>
#include <evt/ShowerSimData.h>
#include <evt/DefaultShowerGeometryProducer.h>

#include <sevt/SEvent.h>
#include <sevt/Station.h>

#include <cppunit/extensions/HelperMacros.h>
#include <tst/Verify.h>

using namespace det;
using namespace std;
using namespace fwk;
using namespace utl;
using namespace tst;
using namespace sdet;
using namespace evt;


class SDenseDetectorTestBasic : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(SDenseDetectorTestBasic);
  CPPUNIT_TEST_EXCEPTION(testNoSimShower, MissingEventDataException);
  CPPUNIT_TEST_EXCEPTION(testNoCorePosition, MissingEventDataException);
  CPPUNIT_TEST(testBasic1);
  CPPUNIT_TEST_SUITE_END();

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

    // CentralConfig must be created first.
    //
    CentralConfig::GetInstance(BOOTSTRAPFILE);

    // Create a detector
    //
    Detector::GetInstance().Update(UTCDateTime(2005,1,1).GetTimeStamp());
  }

  void tearDown() { }

  void
  testBasic1()
  {
    const Detector& theDetector = Detector::GetInstance();

    /*const SDetector& sDet = theDetector.GetSDetector();
    cout << "stations:";
    for (const auto& id : sDet.GetFullStationList())
      cout << ' ' << id;
    cout << endl;*/
    const CoordinateSystemPtr pampaCS = theDetector.GetReferenceCoordinateSystem();
    const CoordinateSystemPtr malargueCS = theDetector.GetSiteCoordinateSystem();

    // Set a sim shower data so that dense station position can
    // be computed relative to it
    //
    Event& evt = RunController::GetInstance().GetCurrentEvent();
    if (!evt.HasSimShower())
      evt.MakeSimShower(DefaultShowerGeometryProducer());

    ShowerSimData& simShower = evt.GetSimShower();
    simShower.SetGroundParticleCoordinateSystemAzimuth(0);
    simShower.SetGroundParticleCoordinateSystemZenith(0);
    const double x = 123*m;
    const double y = 456*m;
    const double z = -13*m;
    //const PointUTM coreUTM(
    simShower.MakeGeometry(Point(x, y, z, pampaCS));
    const double tolerance = 1*m;

    // Check that non-dense stations are not screwed up by presence
    // of the dense manager.
    const Station& nonDense = theDetector.GetSDetector().GetStation(57);
    CPPUNIT_ASSERT(Verify<CloseTo>(nonDense.GetPosition().GetX(malargueCS), 14563.7, 1*mm));
    CPPUNIT_ASSERT(Verify<Equal>(nonDense.GetName(), string("Skeeter")));
    CPPUNIT_ASSERT(Verify<Equal>(nonDense.IsDense(), false));
    CPPUNIT_ASSERT(Verify<Equal>(nonDense.IsInGrid(), true));

    // Check basics of dense stations
    const Station& stat1 = theDetector.GetSDetector().GetStation(99999); // a dense-array station
    CPPUNIT_ASSERT(Verify<Equal>(stat1.GetId(), 99999));
    CPPUNIT_ASSERT(Verify<Equal>(stat1.GetName(), string("DenseStation99999")));
    CPPUNIT_ASSERT(Verify<Equal>(stat1.IsDense(), true));
    CPPUNIT_ASSERT(Verify<Equal>(stat1.IsInGrid(), false)); // dense stations are off grid

    const Point pos1 = stat1.GetPosition();
    CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetX(pampaCS), x, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetY(pampaCS), y, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetZ(pampaCS), z, tolerance));

    const Station& stat2 = theDetector.GetSDetector().GetStation(99998);
    const Point pos2 = stat2.GetPosition();
    CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetX(pampaCS), x + 1000*m, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetY(pampaCS), y, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetZ(pampaCS), z, tolerance));

    const Station& stat3 = theDetector.GetSDetector().GetStation(99996);
    const Point pos3 = stat3.GetPosition();
    CPPUNIT_ASSERT(Verify<CloseTo>(pos3.GetX(pampaCS), x, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos3.GetY(pampaCS), y + 1000*m, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos3.GetZ(pampaCS), z, tolerance));

    const Station& stat4 = theDetector.GetSDetector().GetStation(99990);
    const Point pos4 = stat4.GetPosition();
    CPPUNIT_ASSERT(Verify<CloseTo>(pos4.GetX(pampaCS), x + 1000*m, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos4.GetY(pampaCS), y, tolerance));
    CPPUNIT_ASSERT(Verify<CloseTo>(pos4.GetZ(pampaCS), z, tolerance));
  }

  void
  testNoSimShower()
  {
    Detector& theDetector = Detector::GetInstance();

    // should throw exception, since no SimShower exists.
    Point pos = theDetector.GetSDetector().GetStation(99999).GetPosition();
  }

  void
  testNoCorePosition()
  {
    Event& evt = RunController::GetInstance().GetCurrentEvent();
    if (!evt.HasSimShower())
      evt.MakeSimShower(DefaultShowerGeometryProducer());

    Detector& theDetector = Detector::GetInstance();

    // should throw exception, since shower core position is not set
    Point pos = theDetector.GetSDetector().GetStation(99999).GetPosition();
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(SDenseDetectorTestBasic);