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 | |
---|
47 | using namespace casa; |
---|
48 | |
---|
49 | namespace asap |
---|
50 | { |
---|
51 | |
---|
52 | class MSFillerUtils { |
---|
53 | protected: |
---|
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 | |
---|
123 | class MSFiller |
---|
124 | { |
---|
125 | public: |
---|
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 | |
---|
133 | protected: |
---|
134 | |
---|
135 | |
---|
136 | private: |
---|
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 |
---|