source: trunk/src/SDWriter.cc@ 33

//#---------------------------------------------------------------------------
//# SDWriter.cc: ASAP class to write out single dish spectra.
//#---------------------------------------------------------------------------
//# Copyright (C) 2004
//# Mark Calabretta, ATNF
//#
//# This program is free software; you can redistribute it and/or modify it
//# under the terms of the GNU General Public License as published by the Free
//# Software Foundation; either version 2 of the License, or (at your option)
//# any later version.
//#
//# This program is distributed in the hope that it will be useful, but
//# WITHOUT ANY WARRANTY; without even the implied warranty of
//# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General
//# Public License for more details.
//#
//# You should have received a copy of the GNU General Public License along
//# with this program; if not, write to the Free Software Foundation, Inc.,
//# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning this software should be addressed as follows:
//#        Internet email: Malte.Marquarding@csiro.au
//#        Postal address: Malte Marquarding,
//#                        Australia Telescope National Facility,
//#                        P.O. Box 76,
//#                        Epping, NSW, 2121,
//#                        AUSTRALIA
//#
//# $Id: SDWriter.cc 28 2004-07-07 01:24:23Z mcalabre $
//#---------------------------------------------------------------------------

#include <string>

#include <aips/aips.h>
#include <aips/Arrays.h>
#include <aips/Mathematics/Complex.h>
#include <aips/Utilities/CountedPtr.h>

#include <atnf/PKSIO/PKSMS2writer.h>
#include <atnf/PKSIO/PKSSDwriter.h>

#include <SDMemTable.h>
#include <SDWriter.h>


using namespace atnf_sd;

//--------------------------------------------------------- SDWriter::SDWriter

// Default constructor.

SDWriter::SDWriter(
        const string &format = "SDFITS")
{
  cFormat = format;

  if (cFormat == "MS2") {
    cWriter = new PKSMS2writer();
  } else {
    cWriter = new PKSSDwriter();
  }
}

//-------------------------------------------------------- SDWriter::~SDWriter

// Destructor.

SDWriter::~SDWriter()
{
  // Nothing.
}

//-------------------------------------------------------- SDWriter::setFormat

// Reset the output format.

Int SDWriter::setFormat(
        const string &format = "SDFITS")
{
  if (format != cFormat) {
    delete cWriter;

    cFormat = format;
    if (cFormat == "MS2") {
      cWriter = new PKSMS2writer();
    } else {
      cWriter = new PKSSDwriter();
    }
  }

  return 0;
}

//------------------------------------------------------------ SDWriter::write

// Write an SDMemTable to file in the desired format, closing the file when
// finished.

Int SDWriter::write(
        const CountedPtr<SDMemTable> table,
        const string &filename)
{
  // Extract the header from the table.
  SDHeader hdr = table->getSDHeader();
  Int nPol  = hdr.npol;
  Int nChan = hdr.nchan;

  // Create the output file and write static data.
  Int status;
  if (status = cWriter->create(filename, hdr.observer, hdr.project,
                               hdr.antennaname, hdr.antennaposition,
                               hdr.obstype, hdr.equinox, hdr.freqref,
                               nChan, nPol, False, False)) {
    cerr << "Failed to create output file." << endl;
    return 1;
  }

  Int scanNo = -1;
  Int cycleNo;
  Double mjd0 = 0.0;

  Int count = 0;
  for (Int iRow = 0; iRow < table->nRows(); iRow++) {
    // Extract the next integration from the table.
    SDContainer sd = table->getSDContainer(iRow);
    if (sd.scanid != scanNo) {
      scanNo = sd.scanid;
      mjd0 = sd.timestamp;
      cycleNo = 1;
    } else if (fabs(sd.timestamp-mjd0) > sd.interval) {
      cycleNo++;
    }

    // Write it out beam by beam.
    for (Int iBeam = 0; iBeam < hdr.nbeam; iBeam++) {

      // Write it out IF by IF.
      for (Int iIF = 0; iIF < hdr.nif; iIF++) {
        // None of these are stored in SDMemTable by SDReader.
        String          fieldName = "";
        Vector<Double>  srcDir(2, 0.0);
        Vector<Double>  srcPM(2, 0.0);
        Double          srcVel = 0.0;
        Double          freqInc = 0.0;
        Vector<Float>   tcal(2, 0.0f);
        String          tcalTime = "";
        Float           azimuth = 0.0f;
        Float           elevation = 0.0f;
        Float           parAngle = 0.0f;
        Float           focusAxi = 0.0f;
        Float           focusTan = 0.0f;
        Float           focusRot = 0.0f;
        Float           temperature = 0.0f;
        Float           pressure = 0.0f;
        Float           humidity = 0.0f;
        Float           windSpeed = 0.0f;
        Float           windAz = 0.0f;
        Int             refBeam = 0;
        Vector<Double>  direction(2, 0.0);
        Vector<Double>  scanRate(2, 0.0);
        Vector<Float>   sigma(nPol, 0.0f);
        Vector<Float>   calFctr(nPol, 0.0f);
        Matrix<Float>   baseLin(nPol,2, 0.0f);
        Matrix<Float>   baseSub(nPol,9, 0.0f);
        Complex         xCalFctr;
        Vector<Complex> xPol;

        if (status = cWriter->write(sd.scanid, cycleNo, sd.timestamp,
                                    sd.interval, fieldName, sd.sourcename,
                                    srcDir, srcPM, srcVel, iIF+1, hdr.reffreq,
                                    hdr.bandwidth, freqInc, tcal, tcalTime,
                                    azimuth, elevation, parAngle,
                                    focusAxi, focusTan, focusRot, temperature,
                                    pressure, humidity, windSpeed, windAz,
                                    refBeam, iBeam+1, direction, scanRate,
                                    sd.getTsys(iBeam, iIF),
                                    sigma, calFctr, baseLin, baseSub,
                                    sd.getSpectrum(iBeam, iIF),
                                    sd.getFlags(iBeam, iIF),
                                    xCalFctr, xPol)) {
          cerr << "Error writing output file." << endl;
          return 1;
        }

        count++;
      }
    }
  }

  cout << "SDWriter: wrote " << count << " rows to " << filename << endl;
  cWriter->close();

  return 0;
}