testSDenseDetectorAzimuth.cc

Go to the documentation of this file.

#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/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 <tst/Verify.h>

#include <cppunit/extensions/HelperMacros.h>

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

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


// This file contains a few additional tests of the dense array
// geometry.  These tests are run separately from the other geometry
// tests as the SDenseArrayXMLManager accesses
// shower info through RunController::GetCurrentEvent(), and the
// RunController is a singleton.  There are probably other ways to
// deal with this, but this was convenient and quick to implement.

class SDenseDetectorTestAzimuth : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(SDenseDetectorTestAzimuth);
  CPPUNIT_TEST(testBasicAzimuth);
  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
  testBasicAzimuth()
  {
    Detector& theDetector = Detector::GetInstance();

    // get reference coordinate system of detector (usually PampaAmarilla)
    const utl::CoordinateSystemPtr referenceCS = theDetector.GetReferenceCoordinateSystem();

    // 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(90*degree);
    simShower.SetGroundParticleCoordinateSystemZenith(0);
    const double x = 8*km;
    const double y = -12*km;
    const double z = 0;
    const utl::Point core(x, y, z, referenceCS);
    simShower.MakeGeometry(core);

    /*
      RU Do 27. Jul 12:22:46 CEST 2017

      interesting observation:

      in the new ShowerSimData logic, where the shower azimuth is obtained from the utl::Vector of the axis,
      for the special case of zenith==0 the azimuth is not well defined and not stored in the axis.
      It is always zero (and in fact irrelevant)

      Thus, in this test the symmetry around z is perfect and the azimuth doesn't matter. I modify the
      test in a way to only check for the radial distance, not the exact location of stations...
     */

    const Station& stat1 = theDetector.GetSDetector().GetStation(99998);
    const Point& pos1 = stat1.GetPosition();
    const double distance1 = (pos1 - core).GetMag();
    CPPUNIT_ASSERT(Verify<CloseTo>(distance1, 1000*m, 1*m));

    //CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetX(referenceCS), 0.*m));
    //CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetY(referenceCS), 1000.*m));
    //CPPUNIT_ASSERT(Verify<CloseTo>(pos1.GetZ(referenceCS), 1400.*m));

    const Station& stat2 = theDetector.GetSDetector().GetStation(99997);
    const Point& pos2 = stat2.GetPosition();
    const double distance2 = (pos2 - core).GetMag();
    CPPUNIT_ASSERT(Verify<CloseTo>(distance2, 1000*m, 1*m));

    //CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetX(referenceCS), -1000.*m));
    //CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetY(referenceCS), 0.));
    //CPPUNIT_ASSERT(Verify<CloseTo>(pos2.GetZ(referenceCS), 1400.*m));
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(SDenseDetectorTestAzimuth);