source: trunk/external-alma/asdm2ASAP/ASDMFiller.h@ 2815

Last change on this file since 2815 was 2754, checked in by Takeshi Nakazato, 12 years ago

New Development: No

JIRA Issue: Yes CSV-2532 (may be related to CSV-1908 and CSV-2161)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed: In tool level, added parameter 'freq_tolsr' to

scantable constructor and function sd.splitant.

Test Programs: test_sdsave, test_importasdm_sd

Put in Release Notes: Yes

Module(s): Module Names change impacts.

Description: Describe your changes here...

In importing MS to Scantable, frequency frame information is
imported as is by default, i.e., base frame in Scantable is
TOPO for ALMA data, which is forcibly converted to LSRK with
wrong time and direction reference.

Some functions have a boolean parameter 'freq_tolsr' that controls
the above behavior. If freq_tolsr is False (default), frequency
is imported as is, while frequency is converted to LSRK (wrongly)
when it is True.


File size: 3.1 KB
Line 
1#ifndef ASAP_ASDM_FILLER_H
2#define ASAP_ASDM_FILLER_H
3
4#include <string>
5
6#include <casa/Logging/LogSinkInterface.h>
7
8#include <FillerBase.h>
9#include "ASDMReader.h"
10
11class ASDMFiller : public asap::FillerBase
12{
13public:
14 // constructor and destructor
15 ASDMFiller( casa::CountedPtr<asap::Scantable> stable ) ;
16 ~ASDMFiller() ;
17
18 // open data
19 bool open( const std::string &filename, const casa::Record &rec ) ;
20
21 // fill data
22 void fill() ;
23
24 // close data
25 void close() ;
26
27 // get reader object
28 casa::CountedPtr<ASDMReader> getReader() { return reader_ ; } ;
29
30 // set logger
31 void setLogger( casa::CountedPtr<casa::LogSinkInterface> &logsink ) ;
32
33private:
34 // fill header
35 void fillHeader() ;
36
37 // get IF key
38 casa::String getIFKey( casa::uInt ifno ) ;
39
40 // get FREQUENCIES attributes from ifrec_
41 void getFrequencyRec( casa::String key,
42 double &refpix,
43 double &refval,
44 double &incr ) ;
45
46 // set FREQUENCIES attributes to ifrec_
47 void setFrequencyRec( casa::String key,
48 double refpix,
49 double refval,
50 double incr ) ;
51
52 // reshape float array spectra to Matrix<Float>
53 casa::Matrix<casa::Float> toMatrix( float *sp,
54 unsigned int npol,
55 unsigned int nchan ) ;
56
57 // reshape 2d vector Tsys to Matrix<Float>
58 casa::Matrix<casa::Float> toMatrix( std::vector< std::vector<float> > &tsys,
59 unsigned int npol,
60 unsigned int nchan ) ;
61
62 // reshape vector<float> to Vector<Float> with appropriate length
63 casa::Vector<casa::Float> toVector( std::vector<float> &tau,
64 unsigned int npol ) ;
65
66 // create TCAL time string from MJD
67 casa::String toTcalTime( casa::Double mjd ) ;
68
69 // AZEL to J2000
70 void toJ2000( casa::Vector<casa::Double> &dir,
71 double az,
72 double el,
73 casa::Double mjd,
74 casa::Vector<casa::Double> antpos ) ;
75
76 // to J2000
77 casa::Vector<casa::Double> toJ2000( casa::Vector<casa::Double> dir,
78 casa::String dirref,
79 casa::Double mjd,
80 casa::Vector<casa::Double> antpos ) ;
81
82 // get frequency frame enum value from string
83 casa::MFrequency::Types toFrameType( std::string &s ) ;
84
85 // to LSRK
86 // utc must be UTC time in "d" (day)
87 // antpos must be ITRF value in "m"
88 casa::Double toLSRK( casa::Double freq,
89 casa::String freqref,
90 casa::Double utc,
91 casa::Vector<casa::Double> antpos,
92 casa::Vector<casa::Double> dir,
93 casa::String dirref ) ;
94
95 casa::CountedPtr<ASDMReader> reader_ ;
96 casa::Int antennaId_ ;
97 casa::String antennaName_ ;
98
99 casa::Record ifrec_ ;
100
101 casa::CountedPtr<casa::LogSinkInterface> logsink_ ;
102
103 casa::String className_ ;
104 casa::Bool freqToLsr_ ;
105
106} ;
107#endif // ASAP_ASDM_FILLER_H
Note: See TracBrowser for help on using the repository browser.