source: trunk/src/MSFiller.h@ 2166

Last change on this file since 2166 was 2021, checked in by Takeshi Nakazato, 14 years ago

New Development: No

JIRA Issue: Yes CAS-2718

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...

Bug fix related to ALMA data.

  1. bug fix when NUM_LINES is 0 in SOURCE subtable
  2. support variable number of polarization data (e.g. including WVR data)
  3. bug fix on WEATHER_ID numbering in MSFiller


File size: 3.4 KB
RevLine 
[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
19// Boost
20#include <boost/pool/object_pool.hpp>
21
[1974]22// AIPS++
23#include <casa/aips.h>
24#include <casa/Utilities/CountedPtr.h>
25#include <casa/Arrays/Vector.h>
[1987]26#include <casa/Arrays/Matrix.h>
[1993]27#include <casa/Arrays/Cube.h>
[1974]28#include <casa/Logging/LogIO.h>
29
30#include <casa/Containers/Record.h>
[2002]31#include <casa/Containers/Block.h>
[1974]32
[1987]33//#include <tables/Tables/TableColumn.h>
34
[1974]35#include <ms/MeasurementSets/MeasurementSet.h>
[1987]36#include <ms/MeasurementSets/MSPointing.h>
[1974]37
38#include "Scantable.h"
39
40namespace asap
41{
42
43class MSFiller
44{
45public:
46 explicit MSFiller(casa::CountedPtr<Scantable> stable) ;
47 virtual ~MSFiller() ;
48
49 virtual bool open(const std::string& filename, const casa::Record& rec) ;
50 virtual void fill() ;
51 virtual void close() ;
52
53protected:
54
55
56private:
57
58 MSFiller();
59 MSFiller(const MSFiller&);
60 MSFiller& operator=(const MSFiller&);
61
62 // fill subtables
63 //void fillFrequencies() ;
64 //void fillMolecules() ;
65 void fillWeather() ;
66 void fillFocus() ;
67 //void fillHistory() ;
68 //void fillFit() ;
[2002]69 void fillTcal( boost::object_pool<casa::ROTableColumn> *poolr ) ;
[1974]70
71 // get SRCTYPE from STATE_ID
[1987]72 casa::Int getSrcType( casa::Int stateId, boost::object_pool<casa::ROTableColumn> *pool ) ;
[1974]73
74 // get POLNO from CORR_TYPE
[2002]75 casa::Block<casa::uInt> getPolNo( casa::Int corrType ) ;
[1974]76
77 // get poltype from CORR_TYPE
78 casa::String getPolType( casa::Int corrType ) ;
79
80 // get WEATHER_ID
81 casa::uInt getWeatherId( casa::uInt idx, casa::Double wtime ) ;
82
83 // get time stamp in SYSCAL table
84 // assume that tab is selected by ANTENNA_ID, FEED_ID, SPECTRAL_WINDOW_ID
85 // and sorted by TIME
[2017]86 void getSysCalTime( casa::Vector<casa::MEpoch> &scTimeIn, casa::Vector<casa::Double> &scInterval, casa::Block<casa::MEpoch> &tcol, casa::Block<casa::Int> &tidx ) ;
[1974]87
88 // get TCAL_ID
[2004]89 casa::Block<casa::uInt> getTcalId( casa::Int feedId, casa::Int spwId, casa::MEpoch &t ) ;
[1974]90
91 // get direction for DIRECTION, AZIMUTH, and ELEVATION columns
[1987]92 casa::uInt getDirection( casa::uInt idx, casa::Vector<casa::Double> &dir, casa::Vector<casa::Double> &srate, casa::String &ref, casa::MSPointing &tab, casa::Double t ) ;
[1990]93
94 // create key for TCAL table
95 casa::String keyTcal( casa::Int feedid, casa::Int spwid, casa::String stime ) ;
[2003]96
97 // binary search
98 casa::uInt binarySearch( casa::Vector<casa::MEpoch> &timeList, casa::Double target ) ;
[2017]99
100 // tool for HPC
101 double gettimeofday_sec() ;
[1974]102
103 casa::CountedPtr<Scantable> table_ ;
104 casa::MeasurementSet mstable_ ;
[1987]105 casa::String tablename_ ;
[1974]106 casa::Int antenna_ ;
[2021]107 casa::String antennaStr_ ;
[1974]108 casa::Bool getPt_ ;
109
110 casa::Bool isFloatData_ ;
111 casa::Bool isData_ ;
112
113 casa::Bool isDoppler_ ;
114 casa::Bool isFlagCmd_ ;
115 casa::Bool isFreqOffset_ ;
116 casa::Bool isHistory_ ;
117 casa::Bool isProcessor_ ;
118 casa::Bool isSysCal_ ;
119 casa::Bool isWeather_ ;
120
[1993]121 casa::String colTsys_ ;
122 casa::String colTcal_ ;
123
[1974]124 casa::LogIO os_ ;
125
126 casa::Vector<casa::Double> mwTime_ ;
127 casa::Vector<casa::Double> mwInterval_ ;
128
[1987]129 // Record for TCAL_ID
130 // "FIELD0": "SPW0": Vector<uInt>
131 // "SPW1": Vector<uInt>
132 // ...
[1974]133 casa::Record tcalrec_ ;
[1987]134
135 //casa::ROTableColumn *scCol_ ;
[1974]136};
137
138
139};
140#endif
Note: See TracBrowser for help on using the repository browser.