[1974] | 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> |
---|
[1987] | 18 | |
---|
[1974] | 19 | // AIPS++ |
---|
| 20 | #include <casa/aips.h> |
---|
| 21 | #include <casa/Utilities/CountedPtr.h> |
---|
| 22 | #include <casa/Arrays/Vector.h> |
---|
[1987] | 23 | #include <casa/Arrays/Matrix.h> |
---|
[1993] | 24 | #include <casa/Arrays/Cube.h> |
---|
[1974] | 25 | #include <casa/Logging/LogIO.h> |
---|
[2291] | 26 | #include <casa/Containers/RecordField.h> |
---|
[1974] | 27 | #include <casa/Containers/Record.h> |
---|
[2002] | 28 | #include <casa/Containers/Block.h> |
---|
[2291] | 29 | #include <casa/Quanta/MVTime.h> |
---|
[1974] | 30 | |
---|
| 31 | #include <ms/MeasurementSets/MeasurementSet.h> |
---|
[1987] | 32 | #include <ms/MeasurementSets/MSPointing.h> |
---|
[1974] | 33 | |
---|
[2291] | 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" |
---|
[1974] | 44 | #include "Scantable.h" |
---|
[2291] | 45 | #include "MathUtils.h" |
---|
[1974] | 46 | |
---|
[2291] | 47 | using namespace casa; |
---|
| 48 | |
---|
[1974] | 49 | namespace asap |
---|
| 50 | { |
---|
| 51 | |
---|
[2291] | 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 | |
---|
[1974] | 123 | class MSFiller |
---|
| 124 | { |
---|
| 125 | public: |
---|
[2291] | 126 | explicit MSFiller(CountedPtr<Scantable> stable) ; |
---|
[1974] | 127 | virtual ~MSFiller() ; |
---|
| 128 | |
---|
[2291] | 129 | virtual bool open(const std::string& filename, const Record& rec) ; |
---|
[1974] | 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() ; |
---|
[2291] | 149 | void fillTcal() ; |
---|
[1974] | 150 | |
---|
[1990] | 151 | // create key for TCAL table |
---|
[2291] | 152 | String keyTcal( Int feedid, Int spwid, String stime ) ; |
---|
[2003] | 153 | |
---|
[2217] | 154 | // get frequency frame |
---|
| 155 | std::string getFrame() ; |
---|
[2237] | 156 | |
---|
| 157 | // initialize header |
---|
| 158 | void initHeader( STHeader &header ) ; |
---|
| 159 | |
---|
[2291] | 160 | CountedPtr<Scantable> table_ ; |
---|
| 161 | MeasurementSet mstable_ ; |
---|
| 162 | String tablename_ ; |
---|
| 163 | Int antenna_ ; |
---|
| 164 | String antennaStr_ ; |
---|
| 165 | Bool getPt_ ; |
---|
[2237] | 166 | |
---|
[2291] | 167 | Bool isFloatData_ ; |
---|
| 168 | Bool isData_ ; |
---|
[2237] | 169 | |
---|
[2291] | 170 | Bool isDoppler_ ; |
---|
| 171 | Bool isFlagCmd_ ; |
---|
| 172 | Bool isFreqOffset_ ; |
---|
| 173 | Bool isHistory_ ; |
---|
| 174 | Bool isProcessor_ ; |
---|
| 175 | Bool isSysCal_ ; |
---|
| 176 | Bool isWeather_ ; |
---|
[2237] | 177 | |
---|
[2291] | 178 | String colTsys_ ; |
---|
| 179 | String colTcal_ ; |
---|
[1974] | 180 | |
---|
[2291] | 181 | LogIO os_ ; |
---|
[1974] | 182 | |
---|
[2291] | 183 | Vector<Double> mwTime_ ; |
---|
| 184 | Vector<Double> mwInterval_ ; |
---|
| 185 | Vector<uInt> mwIndex_ ; |
---|
[1974] | 186 | |
---|
[1987] | 187 | // Record for TCAL_ID |
---|
| 188 | // "FIELD0": "SPW0": Vector<uInt> |
---|
| 189 | // "SPW1": Vector<uInt> |
---|
| 190 | // ... |
---|
[2291] | 191 | Record tcalrec_ ; |
---|
| 192 | //map< String,Vector<uInt> > tcalrec_ ; |
---|
[1974] | 193 | }; |
---|
| 194 | |
---|
| 195 | |
---|
| 196 | }; |
---|
| 197 | #endif |
---|