source: branches/Release12/src/SDMemTableWrapper.h@ 2861

Last change on this file since 2861 was 800, checked in by mar637, 19 years ago

Request: added channel based flagging

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