source: branches/Release-1-fixes/src/SDMemTableWrapper.h@ 1435

Last change on this file since 1435 was 539, checked in by mar637, 20 years ago

Reworked stokestopolspectrum to have polIdx=-1 to enable polSel_.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 7.2 KB
Line 
1//#---------------------------------------------------------------------------
2//# SDMemTableWrapper.h: Wrapper classes to use CountedPtr
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2004
5//# 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 SDMEMTABLEWRAPPER_H
32#define SDMEMTABLEWRAPPER_H
33
34#include <vector>
35#include <string>
36#include <casa/Arrays/Vector.h>
37
38#include "MathUtils.h"
39#include "SDFitTable.h"
40#include "SDMemTable.h"
41
42namespace asap {
43
44class SDMemTableWrapper {
45
46public:
47 SDMemTableWrapper(const std::string& name) :
48 table_(new SDMemTable(name)) {;}
49 SDMemTableWrapper() :
50 table_(new SDMemTable()) {;}
51
52 SDMemTableWrapper(casa::CountedPtr<SDMemTable> cp) : table_(cp) {;}
53 //SDMemTableWrapper(SDMemTable* sdmt) : table_(sdmt) {;}
54 SDMemTableWrapper(const SDMemTableWrapper& mt) :
55 table_(mt.getCP()) {;}
56
57 SDMemTableWrapper(const SDMemTableWrapper& mt, const std::string& expr) :
58 table_(new SDMemTable(mt.getCP()->table(), expr)) {;}
59
60 SDMemTableWrapper copy() {
61 //CountedPtr<SDMemTable> cp = new SDMemTable(*this, False);
62 return SDMemTableWrapper(new SDMemTable(*(this->getCP()), casa::False));
63 }
64
65 SDMemTableWrapper getScan(std::vector<int> scan) {
66 casa::String cond("SELECT FROM $1 WHERE SCANID IN ");
67 casa::Vector<casa::Int> v(scan);
68 casa::ostringstream oss;
69 oss << v;
70 cond += casa::String(oss);
71 return SDMemTableWrapper(*this, cond);
72 }
73
74 SDMemTableWrapper getSource(const std::string& source) {
75 casa::String cond("SELECT * from $1 WHERE SRCNAME == pattern('");
76 cond += source;cond += "')";
77 return SDMemTableWrapper(*this, cond);
78 }
79
80 SDMemTableWrapper getSQL(const std::string& sqlexpr) {
81 return SDMemTableWrapper(*this, sqlexpr);
82 }
83
84 std::vector<float> getSpectrum(int whichRow=0) const {
85 return table_->getSpectrum(whichRow);
86 }
87
88 std::vector<float> getStokesSpectrum(int whichRow=0,
89 bool linPol=false) const {
90 return table_->getStokesSpectrum(whichRow, linPol);
91 }
92
93 std::vector<float> stokesToPolSpectrum(int whichRow=0, bool linear=false,
94 int polIdx=-1) const {
95 return table_->stokesToPolSpectrum(whichRow, linear, polIdx);
96 }
97
98 // std::string getPolarizationLabel(bool linear, bool stokes, bool linPol, int polIdx) const {
99 // Boost fails with 4 arguments.
100 std::string getPolarizationLabel(bool linear, bool stokes,
101 bool linPol) const {
102 int polIdx = -1;
103 return table_->getPolarizationLabel(linear, stokes, linPol, polIdx);
104 }
105
106 std::vector<double> getAbcissa(int whichRow=0) const {
107 return table_->getAbcissa(whichRow);
108 }
109 std::string getAbcissaString(int whichRow=0) const {
110 return table_->getAbcissaString(whichRow);
111 }
112
113 std::vector<float> getTsys() {
114 int nRow = table_->nRow();
115 std::vector<float> result(nRow);
116 for (uint i=0; i<nRow; i++) {
117 result[i] = table_->getTsys(i);
118 }
119 return result;
120 }
121
122 std::string getTime(int whichRow=0) {return table_->getTime(whichRow);}
123
124 std::string getFluxUnit() const {return table_->getFluxUnit();}
125 void setFluxUnit(const std::string& unit) {table_->setFluxUnit(unit);}
126
127 void setInstrument(const std::string& name) {table_->setInstrument(name);}
128
129 std::vector<bool> getMask(int whichRow=0) const {
130 return table_->getMask(whichRow);
131 }
132 bool setMask(const std::vector<int> mvals) const {
133 return table_->setMask(mvals);
134 }
135
136 void flag(int whichRow=-1) {
137 table_->flag(whichRow);
138 }
139 std::string getSourceName(int whichRow=0) {
140 return table_->getSourceName(whichRow);
141 }
142
143 void setSpectrum(std::vector<float> spectrum, int whichRow=0) {
144 table_->setSpectrum(spectrum, whichRow);
145 }
146
147 bool setIF(int whichIF=0) {return table_->setIF(whichIF);}
148 bool setBeam(int whichBeam=0) {return table_->setBeam(whichBeam);}
149 bool setPol(int whichPol=0) {return table_->setPol(whichPol);}
150
151 int getIF() {return table_->getIF();}
152 int getBeam() {return table_->getBeam();}
153 int getPol() {return table_->getPol();}
154
155 int nIF() {return table_->nIF();}
156 int nBeam() {return table_->nBeam();}
157 int nPol() {return table_->nPol();}
158 int nChan() {return table_->nChan();}
159 int nScan() {return table_->nScan();}
160 int nRow() {return table_->nRow();}
161 int nStokes() {return table_->nStokes();}
162
163 //sets the mask
164 bool setChannels(const std::vector<int>& whichChans) {
165 return setChannels(whichChans);
166 }
167 void makePersistent(const std::string& fname) {
168 table_->makePersistent(fname);
169 }
170
171 bool setRestFreqs(const std::vector<double>& restfreqs,
172 const std::string& unit,
173 const std::vector<std::string>& lines,
174 const std::string& source, int whichIF) {
175 casa::Vector<casa::Double> restFreqs2(restfreqs);
176 return table_->setRestFreqs(restFreqs2, casa::String(unit),
177 lines, casa::String(source),
178 casa::Int(whichIF));
179 }
180
181 void spectralLines() const {table_->spectralLines();}
182
183 std::vector<double> getRestFreqs() {
184 return table_->getRestFreqs();
185 }
186
187 void setCoordInfo(std::vector<string> theinfo) {
188 table_->setCoordInfo(theinfo);
189 }
190 std::vector<string> getCoordInfo() const {
191 return table_->getCoordInfo();
192 }
193
194 casa::CountedPtr<SDMemTable> getCP() const {return table_;}
195 SDMemTable* getPtr() {return &(*table_);}
196
197 std::string summary(bool verbose=false) const {
198 return table_->summary(verbose);
199 }
200
201 std::vector<std::string> getHistory() {
202 return table_->getHistory();
203 }
204 void addHistory(const std::string& hist) {
205 table_->addHistory(hist);
206 }
207
208 void addFit(int whichRow, const std::vector<double>& p,
209 const std::vector<bool>& m, const std::vector<string>& f,
210 const std::vector<int>& c) {
211
212 casa::Vector<casa::Double> p2(p);
213 casa::Vector<casa::Bool> m2(m);
214 casa::Vector<casa::String> f2 = mathutil::toVectorString(f);
215 casa::Vector<casa::Int> c2(c);
216 table_->addFit(casa::uInt(whichRow), p2,m2,f2,c2);
217 }
218 SDFitTable getSDFitTable(int whichRow) {
219 return table_->getSDFitTable(casa::uInt(whichRow));
220 }
221
222
223private:
224 casa::CountedPtr<SDMemTable> table_;
225};
226
227} // namespace
228#endif
229
Note: See TracBrowser for help on using the repository browser.