casa6

/*

 * Scantable2MSReader.cc

 *

 *  Created on: Jan 5, 2016

 *      Author: nakazato

 */

​

#include <singledishfiller/Filler/Scantable2MSReader.h>

​

#include <iostream>

#include <string>

​

#include <casacore/casa/OS/File.h>

#include <casacore/casa/Containers/Record.h>

#include <casacore/casa/Utilities/Regex.h>

#include <casacore/tables/Tables/TableRecord.h>

#include <casacore/casa/Arrays/Vector.h>

#include <casacore/casa/Arrays/ArrayMath.h>

#include <casacore/casa/IO/ArrayIO.h>

#include <casacore/tables/Tables/Table.h>

#include <casacore/measures/Measures/Stokes.h>

​

using namespace casacore;

​

namespace {

Double queryAntennaDiameter(String const &name) {

  String capitalized = name;

  capitalized.upcase();

  Double diameter = 0.0;

  if (capitalized.matches(Regex(".*(DV|DA|PM)[0-9]+$"))) {

    diameter = 12.0;

  } else if (capitalized.matches(Regex(".*CM[0-9]+$"))) {

    diameter = 7.0;

  } else if (capitalized.contains("GBT")) {

    diameter = 104.9;

  } else if (capitalized.contains("MOPRA")) {

    diameter = 22.0;

  } else if (capitalized.contains("PKS") || capitalized.contains("PARKS")) {

    diameter = 64.0;

  } else if (capitalized.contains("TIDBINBILLA")) {

    diameter = 70.0;

  } else if (capitalized.contains("CEDUNA")) {

    diameter = 30.0;

  } else if (capitalized.contains("HOBART")) {

    diameter = 26.0;

  } else if (capitalized.contains("APEX")) {

    diameter = 12.0;

  } else if (capitalized.contains("ASTE")) {

    diameter = 10.0;

  } else if (capitalized.contains("NRO")) {

    diameter = 45.0;

  }

​

  return diameter;

}

​

template<class T, class U>

U getMapValue(std::map<T, U> const mymap, T const key, U const default_value) {

  auto iter = mymap.find(key);

  if (iter != mymap.end()) {

    return iter->second;

  } else {

    return default_value;