| [2] | 1 | //#---------------------------------------------------------------------------
|
|---|
| 2 | //# SDMemTable.h: A MemoryTable container for single dish integrations
|
|---|
| 3 | //#---------------------------------------------------------------------------
|
|---|
| 4 | //# Copyright (C) 2004
|
|---|
| 5 | //# Malte Marquarding, ATNF
|
|---|
| 6 | //#
|
|---|
| 7 | //# This program is free software; you can redistribute it and/or modify it
|
|---|
| 8 | //# under the terms of the GNU General Public License as published by the Free
|
|---|
| 9 | //# Software Foundation; either version 2 of the License, or (at your option)
|
|---|
| 10 | //# any later version.
|
|---|
| 11 | //#
|
|---|
| 12 | //# This program is distributed in the hope that it will be useful, but
|
|---|
| 13 | //# WITHOUT ANY WARRANTY; without even the implied warranty of
|
|---|
| 14 | //# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
|
|---|
| 15 | //# Public License for more details.
|
|---|
| 16 | //#
|
|---|
| 17 | //# You should have received a copy of the GNU General Public License along
|
|---|
| 18 | //# with this program; if not, write to the Free Software Foundation, Inc.,
|
|---|
| 19 | //# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
|
|---|
| 20 | //#
|
|---|
| 21 | //# Correspondence concerning this software should be addressed as follows:
|
|---|
| 22 | //# Internet email: Malte.Marquarding@csiro.au
|
|---|
| 23 | //# Postal address: Malte Marquarding,
|
|---|
| 24 | //# Australia Telescope National Facility,
|
|---|
| 25 | //# P.O. Box 76,
|
|---|
| 26 | //# Epping, NSW, 2121,
|
|---|
| 27 | //# AUSTRALIA
|
|---|
| 28 | //#
|
|---|
| 29 | //# $Id:
|
|---|
| 30 | //#---------------------------------------------------------------------------
|
|---|
| 31 | #ifndef _SDMEMTABLE_H_
|
|---|
| 32 | #define _SDMEMTABLE_H_
|
|---|
| 33 |
|
|---|
| 34 | // STL
|
|---|
| 35 | #include <string>
|
|---|
| 36 | #include <vector>
|
|---|
| 37 | // AIPS++
|
|---|
| 38 | #include <aips/aips.h>
|
|---|
| 39 | #include <aips/Utilities/String.h>
|
|---|
| 40 | #include <aips/Tables/Table.h>
|
|---|
| 41 | #include <aips/Arrays/MaskedArray.h>
|
|---|
| 42 |
|
|---|
| [39] | 43 | #include <trial/Coordinates/SpectralCoordinate.h>
|
|---|
| 44 |
|
|---|
| [2] | 45 | namespace atnf_sd {
|
|---|
| 46 |
|
|---|
| 47 | class SDContainer;
|
|---|
| 48 | class SDHeader;
|
|---|
| 49 | class SDFrequencyTable;
|
|---|
| 50 |
|
|---|
| 51 |
|
|---|
| 52 | class SDMemTable {
|
|---|
| 53 | public:
|
|---|
| [19] | 54 | // create a new (empty) SDMemTable
|
|---|
| [18] | 55 | SDMemTable();
|
|---|
| [19] | 56 | // create a SDMemTable from a )aips++) table on disk
|
|---|
| [18] | 57 | SDMemTable(const std::string& name);
|
|---|
| [19] | 58 |
|
|---|
| 59 | // Copy Construct a SDMemTable, if clear==True only header and
|
|---|
| 60 | // skeleton are copied, otherwise the whole table is copied.
|
|---|
| [16] | 61 | SDMemTable(const SDMemTable& other, Bool clear=False);
|
|---|
| [2] | 62 |
|
|---|
| [19] | 63 | // Copy Construct a SDMemTable, give a scanid constraint
|
|---|
| 64 | // see also getScan()
|
|---|
| [2] | 65 | SDMemTable(const Table& tab, Int scanID);
|
|---|
| [19] | 66 |
|
|---|
| [2] | 67 | virtual ~SDMemTable();
|
|---|
| [19] | 68 |
|
|---|
| 69 | // put data from meta conatiner into the table
|
|---|
| [39] | 70 | bool putSDContainer(const SDContainer& sdc);
|
|---|
| 71 | bool putSDHeader(const SDHeader& sdh);
|
|---|
| 72 | bool putSDFreqTable(const SDFrequencyTable& sdft);
|
|---|
| [21] | 73 |
|
|---|
| 74 | //get the dat wrapped up in a meta container
|
|---|
| [39] | 75 | SDContainer getSDContainer(uInt whichRow=0) const;
|
|---|
| 76 | SDHeader getSDHeader() const;
|
|---|
| 77 | SDFrequencyTable getSDFreqTable() const;
|
|---|
| [19] | 78 | // get spectrum,mask and tsys for the given row, at the selected
|
|---|
| [21] | 79 | // cursor - all as stl vectors
|
|---|
| 80 | virtual std::vector<float> getSpectrum(Int whichRow);
|
|---|
| [16] | 81 | virtual std::vector<bool> getMask(Int whichRow) const;
|
|---|
| [21] | 82 |
|
|---|
| [19] | 83 | virtual Float getTsys(Int whichRow) const;
|
|---|
| [21] | 84 | // get all as aips++ Vectors
|
|---|
| 85 | virtual void getSpectrum(Vector<Float>& spectrum, Int whichRow=0);
|
|---|
| 86 | virtual void getMask(Vector<Bool>& mask,Int whichRow=0) const;
|
|---|
| [2] | 87 |
|
|---|
| [21] | 88 | // get info for current row
|
|---|
| [2] | 89 | virtual Double getTime(Int whichRow) const ;
|
|---|
| [19] | 90 | virtual std::string getSourceName(Int whichRow) const;
|
|---|
| [2] | 91 |
|
|---|
| [19] | 92 | // set the current value
|
|---|
| [2] | 93 | virtual bool setIF(Int whichIF=0);
|
|---|
| 94 | virtual bool setBeam(Int whichBeam=0);
|
|---|
| 95 | virtual bool setPol(Int whichPol=0);
|
|---|
| [18] | 96 | //sets the user mask
|
|---|
| 97 | virtual bool setMask(const std::vector<int>& whichChans);
|
|---|
| 98 |
|
|---|
| [19] | 99 | // return the currently selected values
|
|---|
| [2] | 100 | virtual Int getIF() { return IFSel_; }
|
|---|
| 101 | virtual Int getBeam() { return beamSel_; }
|
|---|
| 102 | virtual Int getPol() { return polSel_; }
|
|---|
| [18] | 103 |
|
|---|
| [19] | 104 | // print a summary to stdout
|
|---|
| [2] | 105 | virtual void summary() const;
|
|---|
| 106 |
|
|---|
| [19] | 107 | // write to disk as aips++ table
|
|---|
| [2] | 108 | void makePersistent(const std::string& filename);
|
|---|
| [19] | 109 |
|
|---|
| 110 | // get a new SDMemTable containg all rows with the same give SCANID
|
|---|
| [2] | 111 | SDMemTable getScan(Int scanID);
|
|---|
| [18] | 112 |
|
|---|
| [21] | 113 | const TableRecord& getHeader() const {return table_.keywordSet();}
|
|---|
| [19] | 114 | // get a handle to the "raw" aips++ table
|
|---|
| [2] | 115 | const Table& table() { return table_; }
|
|---|
| 116 |
|
|---|
| [19] | 117 | // return the number of values
|
|---|
| [18] | 118 | Int nBeam() const;
|
|---|
| 119 | Int nIF() const;
|
|---|
| 120 | Int nPol() const;
|
|---|
| 121 | Int nChan() const;
|
|---|
| 122 |
|
|---|
| [22] | 123 | // return the number of rows (integrations) in the table
|
|---|
| 124 | Int nRows() const { return table_.nrow(); }
|
|---|
| 125 |
|
|---|
| [19] | 126 | // return a row as a Masked array, internally converting uChar flags
|
|---|
| 127 | // to bool mask
|
|---|
| 128 | MaskedArray<Float> rowAsMaskedArray(uInt whichRow,
|
|---|
| 129 | Bool useSelection = False);
|
|---|
| 130 |
|
|---|
| [39] | 131 | SpectralCoordinate getCoordinate(uInt whichIdx) const;
|
|---|
| 132 | std::vector<double> getAbscissa(int whichRow,
|
|---|
| 133 | const std::string& whichUnit="GHz",
|
|---|
| 134 | double restfreq=0.0);
|
|---|
| [2] | 135 | private:
|
|---|
| [22] | 136 | // set up table structure
|
|---|
| [2] | 137 | void setup();
|
|---|
| [19] | 138 | // the current cursor into the array
|
|---|
| [2] | 139 | Int IFSel_,beamSel_,polSel_;
|
|---|
| [19] | 140 | std::vector<bool> chanMask_;
|
|---|
| [2] | 141 | String name_;
|
|---|
| [19] | 142 | // the unerlying memory table
|
|---|
| [2] | 143 | Table table_;
|
|---|
| 144 | };
|
|---|
| 145 |
|
|---|
| 146 | }// namespace
|
|---|
| 147 | #endif
|
|---|
