source: branches/casa-prerelease/pre-asap/src/MSFiller.h@ 2308

Last change on this file since 2308 was 2293, checked in by Takeshi Nakazato, 13 years ago

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

Removed unused methods in MSWriter.
Cleaned up comments.


File size: 5.1 KB
Line 
1//
2// C++ Interface: MSFiller
3//
4// Description:
5//
6// This class is specific filler for MS format
7//
8// Takeshi Nakazato <takeshi.nakazato@nao.ac.jp>, (C) 2010
9//
10// Copyright: See COPYING file that comes with this distribution
11//
12//
13#ifndef ASAPMSFILLER_H
14#define ASAPMSFILLER_H
15
16// STL
17#include <string>
18
19// AIPS++
20#include <casa/aips.h>
21#include <casa/Utilities/CountedPtr.h>
22#include <casa/Arrays/Vector.h>
23#include <casa/Arrays/Matrix.h>
24#include <casa/Arrays/Cube.h>
25#include <casa/Logging/LogIO.h>
26#include <casa/Containers/RecordField.h>
27#include <casa/Containers/Record.h>
28#include <casa/Containers/Block.h>
29#include <casa/Quanta/MVTime.h>
30
31#include <ms/MeasurementSets/MeasurementSet.h>
32#include <ms/MeasurementSets/MSPointing.h>
33
34#include <tables/Tables/ScalarColumn.h>
35#include <tables/Tables/ArrayColumn.h>
36#include <tables/Tables/TableRow.h>
37
38#include <measures/TableMeasures/ScalarMeasColumn.h>
39#include <measures/TableMeasures/ArrayMeasColumn.h>
40#include <measures/TableMeasures/ScalarQuantColumn.h>
41#include <measures/TableMeasures/ArrayQuantColumn.h>
42
43#include "TableTraverse.h"
44#include "Scantable.h"
45#include "MathUtils.h"
46
47using namespace casa;
48
49namespace asap
50{
51
52class MSFillerUtils {
53protected:
54 template<class T> void getScalar( const String &name,
55 const uInt &idx,
56 const Table &tab,
57 T &val )
58 {
59 ROScalarColumn<T> col( tab, name ) ;
60 val = col( idx ) ;
61 }
62 template<class T> void getArray( const String &name,
63 const uInt &idx,
64 const Table &tab,
65 Array<T> &val )
66 {
67 ROArrayColumn<T> col( tab, name ) ;
68 val = col( idx ) ;
69 }
70 template<class T> void getScalarMeas( const String &name,
71 const uInt &idx,
72 const Table &tab,
73 T &val )
74 {
75 ROScalarMeasColumn<T> measCol( tab, name ) ;
76 val = measCol( idx ) ;
77 }
78 template<class T> void getArrayMeas( const String &name,
79 const uInt &idx,
80 const Table &tab,
81 Array<T> &val )
82 {
83 ROArrayMeasColumn<T> measCol( tab, name ) ;
84 val = measCol( idx ) ;
85 }
86 template<class T> void getScalarQuant( const String &name,
87 const uInt &idx,
88 const Table &tab,
89 Quantum<T> &val )
90 {
91 ROScalarQuantColumn<T> quantCol( tab, name ) ;
92 val = quantCol( idx ) ;
93 }
94 template<class T> void getArrayQuant( const String &name,
95 const uInt &idx,
96 const Table &tab,
97 Array< Quantum<T> > &val )
98 {
99 ROArrayQuantColumn<T> quantCol( tab, name ) ;
100 val = quantCol( idx ) ;
101 }
102 template<class T> T interp( Double x0, Double x1, Double x, T y0, T y1 )
103 {
104 Double dx0 = x - x0 ;
105 Double dx1 = x1 - x ;
106 return ( y0 * dx1 + y1 * dx0 ) / ( x1 - x0 ) ;
107 }
108 String keyTcal( const Int &feedid, const Int &spwid, const Double &time )
109 {
110 String stime = MVTime( Quantity(time,Unit("s")) ).string( MVTime::YMD ) ;
111 String sfeed = "FEED" + String::toString( feedid ) ;
112 String sspw = "SPW" + String::toString( spwid ) ;
113 return sfeed+":"+sspw+":"+stime ;
114 }
115 String keyTcal( const Int &feedid, const Int &spwid, const String &stime )
116 {
117 String sfeed = "FEED" + String::toString( feedid ) ;
118 String sspw = "SPW" + String::toString( spwid ) ;
119 return sfeed+":"+sspw+":"+stime ;
120 }
121};
122
123class MSFiller
124{
125public:
126 explicit MSFiller(CountedPtr<Scantable> stable) ;
127 virtual ~MSFiller() ;
128
129 virtual bool open(const std::string& filename, const Record& rec) ;
130 virtual void fill() ;
131 virtual void close() ;
132
133protected:
134
135
136private:
137
138 MSFiller();
139 MSFiller(const MSFiller&);
140 MSFiller& operator=(const MSFiller&);
141
142 // fill subtables
143 //void fillFrequencies() ;
144 //void fillMolecules() ;
145 void fillWeather() ;
146 void fillFocus() ;
147 //void fillHistory() ;
148 //void fillFit() ;
149 void fillTcal() ;
150
151 // create key for TCAL table
152 String keyTcal( Int feedid, Int spwid, String stime ) ;
153
154 // get frequency frame
155 std::string getFrame() ;
156
157 // initialize header
158 void initHeader( STHeader &header ) ;
159
160 CountedPtr<Scantable> table_ ;
161 MeasurementSet mstable_ ;
162 String tablename_ ;
163 Int antenna_ ;
164 String antennaStr_ ;
165 Bool getPt_ ;
166
167 Bool isFloatData_ ;
168 Bool isData_ ;
169
170 Bool isDoppler_ ;
171 Bool isFlagCmd_ ;
172 Bool isFreqOffset_ ;
173 Bool isHistory_ ;
174 Bool isProcessor_ ;
175 Bool isSysCal_ ;
176 Bool isWeather_ ;
177
178 String colTsys_ ;
179 String colTcal_ ;
180
181 LogIO os_ ;
182
183 Vector<Double> mwTime_ ;
184 Vector<Double> mwInterval_ ;
185 Vector<uInt> mwIndex_ ;
186
187 // Record for TCAL_ID
188 // "FIELD0": "SPW0": Vector<uInt>
189 // "SPW1": Vector<uInt>
190 // ...
191 Record tcalrec_ ;
192 //map< String,Vector<uInt> > tcalrec_ ;
193};
194
195
196};
197#endif
Note: See TracBrowser for help on using the repository browser.