//#--------------------------------------------------------------------------- //# SDMath.h: A collection of single dish mathematical operations //#--------------------------------------------------------------------------- //# Copyright (C) 2004 //# 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: //#--------------------------------------------------------------------------- #ifndef SDMATH_H #define SDMATH_H #include #include #include #include #include "SDDefs.h" class casa::Table; namespace asap { class SDMemTable; class SDMath { public: // Default constructor SDMath(); // Copy Constructor (copy semantics) SDMath (const SDMath& other); // Assignment (copy semantics) SDMath &operator=(const SDMath& other); // Destructor ~SDMath(); // Binary Table operators. op=ADD, SUB, MUL, DIV, QUOTIENT casa::CountedPtr binaryOperate (const casa::CountedPtr& left, const casa::CountedPtr& right, const casa::String& op, casa::Bool preserve) const; // Average in time casa::CountedPtr average(const casa::Block >& in, const casa::Vector& mask, casa::Bool scanAverage, const casa::String& weightStr, casa::Bool alignVelocity=casa::False) const; // Statistics. If row<0, all rows are done otherwise, just the // specified row. std::vector statistic(const casa::CountedPtr& in, const casa::Vector& mask, const casa::String& which, casa::Int row) const; // Bin up spectra SDMemTable* bin(const SDMemTable& in, casa::Int width) const; // Smooth SDMemTable* smooth (const SDMemTable& in, const casa::String& kernel, casa::Float width, casa::Bool doAll) const; // Flux conversion between Jansky and Kelvin SDMemTable* convertFlux (const SDMemTable& in, casa::Float area, casa::Float eta, casa::Bool doAll) const; // Gain-elevation correction SDMemTable* gainElevation (const SDMemTable& in, const casa::Vector& coeffs, const casa::String& fileName, const casa::String& method, casa::Bool doAll) const; // Opacity correction SDMemTable* opacity (const SDMemTable& in, casa::Float tau, casa::Bool doAll) const; // Simple unary mathematical operations. what=0 (mul) or 1 (add) SDMemTable* unaryOperate(const SDMemTable& in, casa::Float offset, casa::Bool doAll, casa::uInt what) const; // Average polarizations SDMemTable* averagePol(const SDMemTable& in, const casa::Vector& mask) const; private: // Weighting type for time averaging enum WeightType {NONE,VAR,TSYS}; // Function to use accumulate data during time averaging void accumulate (casa::Double& timeSum, casa::Double& intSum, casa::Int& nAccum, casa::MaskedArray& sum, casa::Array& sumSq, casa::Array& nPts, casa::Array& tSysSum, const casa::Array& tSys, const casa::Array& nInc, const casa::Vector& mask, casa::Double time, casa::Double interval, const casa::Block >& in, casa::uInt iTab, casa::uInt iRow, casa::uInt axis, casa::uInt nAxesSub, casa::Bool useMask, WeightType wtType) const; // Function to fill Scan Container when averaging in time void fillSDC (SDContainer& sc, const casa::Array& mask, const casa::Array& data, const casa::Array& tSys, casa::Int scanID, casa::Double timeStamp, casa::Double interval, const casa::String& sourceName, const casa::Vector& freqID) const; // Put the data and mask into the SDContainer void putDataInSDC (SDContainer& sc, const casa::Array& data, const casa::Array& mask) const; // Function to normalize data when averaging in time void normalize (casa::MaskedArray& data, const casa::Array& sumSq, const casa::Array& nPts, WeightType wtType, casa::Int axis, casa::Int nAxes) const; // Function to get the current cursor location void getCursorLocation (casa::IPosition& start, casa::IPosition& end, const SDMemTable& in) const; // Convert weight string to enum value void convertWeightString (WeightType& wt, const casa::String& weightStr) const; // Convert interpolation type string void convertInterpString(casa::Int& type, const casa::String& interp) const; // Correct data from an ascii Table void correctFromAsciiTable(SDMemTable* pTabOut, const SDMemTable& in, const casa::String& fileName, const casa::String& col0, const casa::String& col1, const casa::String& methodStr, casa::Bool doAll, const casa::Vector& xOut) const; // Correct data from a Table void correctFromTable(SDMemTable* pTabOut, const SDMemTable& in, const casa::Table& tTable, const casa::String& col0, const casa::String& col1, const casa::String& methodStr, casa::Bool doAll, const casa::Vector& xOut) const; // Correct data from a Vector void correctFromVector (SDMemTable* pTabOut, const SDMemTable& in, casa::Bool doAll, const casa::Vector& factor) const; // Read ascii file into a Table casa::Table readAsciiFile (const casa::String& fileName) const; }; } // namespace #endif