source: trunk/src/SDMemTableWrapper.h@ 525

Last change on this file since 525 was 505, checked in by kil064, 20 years ago

remove pa offset from getSTokesSpectrum

  • 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, bool linPol=false) const {
89 return table_->getStokesSpectrum(whichRow, linPol);
90 }
91
92 std::vector<float> stokesToPolSpectrum (int whichRow, bool linear,
93 int polIdx) const {
94 return table_->stokesToPolSpectrum(whichRow, linear, polIdx);
95 }
96
97// std::string getPolarizationLabel (bool linear, bool stokes, bool linPol, int polIdx) const {
98// Boost fails with 4 arguments.
99 std::string getPolarizationLabel (bool linear, bool stokes, bool linPol) const {
100 int polIdx = -1;
101 return table_->getPolarizationLabel (linear, stokes, linPol, polIdx);
102 }
103//
104 std::vector<double> getAbcissa(int whichRow=0) const {
105 return table_->getAbcissa(whichRow);
106 }
107 std::string getAbcissaString(int whichRow=0) const {
108 return table_->getAbcissaString(whichRow);
109 }
110
111 std::vector<float> getTsys() {
112 int nRow = table_->nRow();
113 std::vector<float> result(nRow);
114 for (uint i=0; i<nRow; i++) {
115 result[i] = table_->getTsys(i);
116 }
117 return result;
118 }
119
120 std::string getTime(int whichRow=0) {return table_->getTime(whichRow);}
121
122 std::string getFluxUnit() const {return table_->getFluxUnit();}
123 void setFluxUnit(const std::string& unit) {table_->setFluxUnit(unit);}
124
125 void setInstrument(const std::string& name) {table_->setInstrument(name);}
126
127 std::vector<bool> getMask(int whichRow=0) const {
128 return table_->getMask(whichRow);
129 }
130 bool setMask(const std::vector<int> mvals) const {
131 return table_->setMask(mvals);
132 }
133
134 void flag(int whichRow=-1) {
135 table_->flag(whichRow);
136 }
137 std::string getSourceName(int whichRow=0) {
138 return table_->getSourceName(whichRow);
139 }
140
141 void setSpectrum(std::vector<float> spectrum, int whichRow=0) {
142 table_->setSpectrum(spectrum, whichRow);
143 }
144
145 bool setIF(int whichIF=0) {return table_->setIF(whichIF);}
146 bool setBeam(int whichBeam=0) {return table_->setBeam(whichBeam);}
147 bool setPol(int whichPol=0) {return table_->setPol(whichPol);}
148
149 int getIF() {return table_->getIF();}
150 int getBeam() {return table_->getBeam();}
151 int getPol() {return table_->getPol();}
152
153 int nIF() {return table_->nIF();}
154 int nBeam() {return table_->nBeam();}
155 int nPol() {return table_->nPol();}
156 int nChan() {return table_->nChan();}
157 int nScan() {return table_->nScan();}
158 int nRow() {return table_->nRow();}
159 int nStokes() {return table_->nStokes();}
160
161 //sets the mask
162 bool setChannels(const std::vector<int>& whichChans) {
163 return setChannels(whichChans);
164 }
165 void makePersistent(const std::string& fname) {
166 table_->makePersistent(fname);
167 }
168
169 bool setRestFreqs(const std::vector<double>& restfreqs,
170 const std::string& unit,
171 const std::vector<std::string>& lines,
172 const std::string& source, int whichIF) {
173 casa::Vector<casa::Double> restFreqs2(restfreqs);
174 return table_->setRestFreqs(restFreqs2, casa::String(unit),
175 lines, casa::String(source),
176 casa::Int(whichIF));
177 }
178
179 void spectralLines() const {table_->spectralLines();}
180
181 std::vector<double> getRestFreqs() {
182 return table_->getRestFreqs();
183 }
184
185 void setCoordInfo(std::vector<string> theinfo) {
186 table_->setCoordInfo(theinfo);
187 }
188 std::vector<string> getCoordInfo() const {
189 return table_->getCoordInfo();
190 }
191
192 casa::CountedPtr<SDMemTable> getCP() const {return table_;}
193 SDMemTable* getPtr() {return &(*table_);}
194
195 std::string summary(bool verbose=false) const {
196 return table_->summary(verbose);
197 }
198
199 std::vector<std::string> getHistory() {
200 return table_->getHistory();
201 }
202 void addHistory(const std::string& hist) {
203 table_->addHistory(hist);
204 }
205
206 void addFit(int whichRow, const std::vector<double>& p,
207 const std::vector<bool>& m, const std::vector<string>& f,
208 const std::vector<int>& c) {
209
210 casa::Vector<casa::Double> p2(p);
211 casa::Vector<casa::Bool> m2(m);
212 casa::Vector<casa::String> f2 = mathutil::toVectorString(f);
213 casa::Vector<casa::Int> c2(c);
214 table_->addFit(casa::uInt(whichRow), p2,m2,f2,c2);
215 }
216 SDFitTable getSDFitTable(int whichRow) {
217 return table_->getSDFitTable(casa::uInt(whichRow));
218 }
219
220
221private:
222 casa::CountedPtr<SDMemTable> table_;
223};
224
225} // namespace
226#endif
227
Note: See TracBrowser for help on using the repository browser.