GenerateDenseArray.cc

Go to the documentation of this file.

#include <sstream>
#include <fstream>
#include <math.h>

#include <utl/Reader.h>

using namespace std;
using namespace utl;

//
// Simple stand-alone dense array generator.
// ----------------------------------------
//   T. Paul
//   13 Jan 2006
//
// This program writes out an XML file
// defining a so-called "dense" array.  In this
// case, the dense array has a ring of stations
// placed at some distance from the shower core.
// This distance and the number of stations
// is defined in the conf.xml file.
// The format used is that expected by the
// SDetenseStationListXMLManager.
//

int
main()
{
  Reader theReader("conf.xml", Reader::eSCHEMA);

  double radius;
  theReader.GetTopBranch().GetChild("radius").GetData(radius);
  int numStations;
  theReader.GetTopBranch().GetChild("numberStations").GetData(numStations);

  ofstream outp("SDenseStationRingAt1000.xml", ios::out);

  // Leading XML boilerplate
  //
  outp <<
    "<?xml version=\"1.0\" encoding=\"iso-8859-1\"?>\n"
    "<SDenseStationList\n"
    "  xmlns:xsi=\"http://www.w3.org/2001/XMLSchema-instance\"\n"
    "  xsi:noNamespaceSchemaLocation=\""
    SCHEMALOCATION
    "/SDenseStationList.xsd\">\n\n";

  outp <<
    "<!--\n"
    "  deltaPhi = polar angle measured in the shower frame\n"
    "  r        = distance from the shower core in plane perpindicular to shower axis\n"
    "  z        = height above the core position\n"
    "  -->\n\n"
    "<UTMInfo>\n"
    "  <ellipsoid> WGS84 </ellipsoid>\n"
    "  <zone>      19    </zone>\n"
    "  <band>      H     </band>\n"
    "</UTMInfo>\n\n";

  // Create list of "dense" stations situated on a ring with
  // radius specified in conf file
  //
  const int firstStation = 90000;   // dense station id's have to be > 90000
  for (int i = 0; i < numStations; ++i) {
    const double deltaPhi = i*2*M_PI/numStations;

    outp <<
      "<denseStation id=\"" << firstStation+i << "\">\n"
      "  <deltaPhi unit=\"radian\"> " << deltaPhi << " </deltaPhi>\n"
      "  <r unit=\"meter\"> " << radius << "</r>\n"
      "</denseStation>\n";
  }

  // Trailing XML boilerplate
  //
  outp << "</SDenseStationList>\n";

  outp.close();

  return EXIT_SUCCESS;
}