source: trunk/external-alma/atnf/PKSIO/NRODataset.h@ 2770

Last change on this file since 2770 was 2766, checked in by Takeshi Nakazato, 12 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...

Use CountedPtr for NRODataRecord object instead of raw C pointer.


File size: 15.9 KB
Line 
1//#---------------------------------------------------------------------------
2//# NRODataset.h: Base class for NRO dataset.
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2000-2006
5//# Associated Universities, Inc. Washington DC, USA.
6//#
7//# This library is free software; you can redistribute it and/or modify it
8//# under the terms of the GNU Library General Public License as published by
9//# the Free Software Foundation; either version 2 of the License, or (at your
10//# option) any later version.
11//#
12//# This library is distributed in the hope that it will be useful, but WITHOUT
13//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14//# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public
15//# License for more details.
16//#
17//# You should have received a copy of the GNU Library General Public License
18//# along with this library; if not, write to the Free Software Foundation,
19//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
20//#
21//# Correspondence concerning AIPS++ should be addressed as follows:
22//# Internet email: aips2-request@nrao.edu.
23//# Postal address: AIPS++ Project Office
24//# National Radio Astronomy Observatory
25//# 520 Edgemont Road
26//# Charlottesville, VA 22903-2475 USA
27//#
28//# $Id$
29//#---------------------------------------------------------------------------
30//# Original: 2009/02/27, Takeshi Nakazato, NAOJ
31//#---------------------------------------------------------------------------
32
33#ifndef NRO_DATASET_H
34#define NRO_DATASET_H
35
36#include <string>
37#include <stdio.h>
38#include <vector>
39#include <iostream>
40
41//#include <casa/aips.h>
42#include <casa/Logging/LogIO.h>
43#include <casa/Containers/Record.h>
44#include <casa/Utilities/CountedPtr.h>
45#include <atnf/PKSIO/NRODataRecord.h>
46#include <casa/namespace.h>
47#include <casa/iomanip.h>
48
49#define SCAN_HEADER_SIZE 424
50
51// <summary>
52// Base class for NRO accessor classes.
53// </summary>
54//
55// <prerequisite>
56// <li> <linkto class=NROReader>NROReader</linkto>
57// <li> <linkto class=NRODataRecord>NRODataRecord</linkto>
58// </prerequisite>
59//
60// <reviewed reviewer="" date="" tests="" demos="">
61// </reviewed>
62//
63// <etymology>
64// This class is a base class for classes that actually access data from NRO telescopes.
65// Concrete classes are defiened for each data type (OTF format or NRO FITS) and/or
66// telescopes (45m or ASTE).
67// The class have two filler method: fillHeader and fillRecord. The former reads header
68// information from the data. Since header data depends on the telescope and its configuration,
69// it is an abstract in this class and is defined in each concrete class.
70// On the other hand, the later reads each scan record (set of meta data
71// and spectral data). The method uses <linkto class=NRODataRecord>NRODataRecord</linkto>
72// to access scan record. It is implemented here since contents of scan record is
73// quite similar for several types of data.
74// </etymology>
75//
76// <synopsis>
77// Abstract class that is designed as a base class for all accessor classes.
78// </synopsis>
79//
80
81class NRODataset
82{
83 public:
84 // Constructor
85 NRODataset( std::string name ) ;
86
87 // Destructor
88 virtual ~NRODataset() ;
89
90 // Data initialization
91 virtual void initialize() ;
92
93 // open file
94 virtual int open() ;
95
96 // close file
97 virtual void close() ;
98
99 // Fill data header from file
100 virtual int fillHeader() = 0 ;
101
102 // Fill data record
103 virtual int fillRecord( int i ) ;
104
105 // simple getter
106 const std::string getLOFIL() const { return LOFIL ; } ;
107 const std::string getVER() const { return VER ; } ;
108 const std::string getGROUP() const { return GROUP ; } ;
109 const std::string getPROJ() const { return PROJ ; } ;
110 const std::string getSCHED() const { return SCHED ; } ;
111 const std::string getOBSVR() const { return OBSVR ; } ;
112 const std::string getLOSTM() const { return LOSTM ; } ;
113 const std::string getLOETM() const { return LOETM ; } ;
114 const int getARYNM() const { return ARYNM ; } ;
115 const int getNSCAN() const { return NSCAN ; } ;
116 const std::string getTITLE() const { return TITLE ; } ;
117 const std::string getOBJ() const { return OBJ ; } ;
118 const std::string getEPOCH() const { return EPOCH ; } ;
119 const double getRA0() const { return RA0 ; } ;
120 const double getDEC0() const { return DEC0 ; } ;
121 const double getGLNG0() const { return GLNG0 ; } ;
122 const double getGLAT0() const { return GLAT0 ; } ;
123 const int getNCALB() const { return NCALB ; } ;
124 const int getSCNCD() const { return SCNCD ; } ;
125 const std::string getSCMOD() const { return SCMOD ; } ;
126 const double getURVEL() const { return URVEL ; } ;
127 const std::string getVREF() const { return VREF ; } ;
128 const std::string getVDEF() const { return VDEF ; } ;
129 const std::string getSWMOD() const { return SWMOD ; } ;
130 const double getFRQSW() const { return FRQSW ; } ;
131 const double getDBEAM() const { return DBEAM ; } ;
132 const double getMLTOF() const { return MLTOF ; } ;
133 const double getCMTQ() const { return CMTQ ; } ;
134 const double getCMTE() const { return CMTE ; } ;
135 const double getCMTSOM() const { return CMTSOM ; } ;
136 const double getCMTNODE() const { return CMTNODE ; } ;
137 const double getCMTI() const { return CMTI ; } ;
138 const std::string getCMTTM() const { return CMTTM ; } ;
139 const double getSBDX() const { return SBDX ; } ;
140 const double getSBDY() const { return SBDY ; } ;
141 const double getSBDZ1() const { return SBDZ1 ; } ;
142 const double getSBDZ2() const { return SBDZ2 ; } ;
143 const double getDAZP() const { return DAZP ; } ;
144 const double getDELP() const { return DELP ; } ;
145 const int getCHBIND() const { return CHBIND ; } ;
146 const int getNUMCH() const { return NUMCH ; } ;
147 const int getCHMIN() const { return CHMIN ; } ;
148 const int getCHMAX() const { return CHMAX ; } ;
149 const double getALCTM() const { return ALCTM ; } ;
150 const double getIPTIM() const { return IPTIM ; } ;
151 const double getPA() const { return PA ; } ;
152 const int getSCNLEN() const { return SCNLEN ; } ;
153 const int getSBIND() const { return SBIND ; } ;
154 const int getIBIT() const { return IBIT ; } ;
155 const std::string getSITE() const { return SITE ; } ;
156 const std::vector<std::string> getRX() const { return RX ; } ;
157 const std::vector<double> getHPBW() const { return HPBW ; } ;
158 const std::vector<double> getEFFA() const { return EFFA ; } ;
159 const std::vector<double> getEFFB() const { return EFFB ; } ;
160 const std::vector<double> getEFFL() const { return EFFL ; } ;
161 const std::vector<double> getEFSS() const { return EFSS ; } ;
162 const std::vector<double> getGAIN() const { return GAIN ; } ;
163 const std::vector<std::string> getHORN() const { return HORN ; } ;
164 const std::vector<std::string> getPOLTP() const { return POLTP ; } ;
165 const std::vector<double> getPOLDR() const { return POLDR ; } ;
166 const std::vector<double> getPOLAN() const { return POLAN ; } ;
167 const std::vector<double> getDFRQ() const { return DFRQ ; } ;
168 const std::vector<std::string> getSIDBD() const { return SIDBD ; } ;
169 const std::vector<int> getREFN() const { return REFN ; } ;
170 const std::vector<int> getIPINT() const { return IPINT ; } ;
171 const std::vector<int> getMULTN() const { return MULTN ; } ;
172 const std::vector<double> getMLTSCF() const { return MLTSCF ; } ;
173 const std::vector<std::string> getLAGWIND() const { return LAGWIND ; } ;
174 const std::vector<double> getBEBW() const { return BEBW ; } ;
175 const std::vector<double> getBERES() const { return BERES ; } ;
176 const std::vector<double> getCHWID() const { return CHWID ; } ;
177 const std::vector<int> getARRY() const { return ARRY ; } ;
178 const std::vector<int> getNFCAL() const { return NFCAL ; } ;
179 const std::vector<double> getF0CAL() const { return F0CAL ; } ;
180 const std::vector< std::vector<double> > getFQCAL() const { return FQCAL ; } ;
181 const std::vector< std::vector<double> > getCHCAL() const { return CHCAL ; } ;
182 const std::vector< std::vector<double> > getCWCAL() const { return CWCAL ; } ;
183 const std::string getCDMY1() const { return CDMY1 ; } ;
184 const std::vector<double> getDSBFC() const { return DSBFC ;} ;
185 const int getDataSize() const { return datasize_ ; } ;
186 const int getRowNum() const { return rowNum_ ; } ;
187
188 // get various parameters
189 NRODataRecord *getRecord( int i ) ;
190 virtual std::vector< std::vector<double> > getSpectrum() ;
191 virtual std::vector<double> getSpectrum( int i ) ;
192 virtual int getIndex( int irow ) ;
193 virtual int getPolarizationNum() ;
194 virtual std::vector<double> getStartIntTime() ;
195 virtual double getStartIntTime( int i ) ;
196 virtual double getScanTime( int i ) ;
197 virtual double getMJD( const char *time ) ;
198 virtual std::vector<bool> getIFs() ;
199 virtual std::vector<double> getFrequencies( int i ) ;
200 virtual uInt getArrayId( std::string type ) ;
201 virtual uInt getPolNo( int irow ) ;
202
203 protected:
204 // fill header information
205 virtual int fillHeader( int sameEndian ) = 0 ;
206
207 // Endian conversion for int variable
208 void convertEndian( int &value ) ;
209
210 // Endian convertion for float variable
211 void convertEndian( float &value ) ;
212
213 // Endian conversion for double variable
214 void convertEndian( double &value ) ;
215
216 // Endian conversion for NRODataRecord
217 void convertEndian( NRODataRecord &r ) ;
218
219 // Read char data
220 int readHeader( char *v, int size ) ;
221
222 // Read int data
223 int readHeader( int &v, int b ) ;
224
225 // Read float data
226 int readHeader( float &v, int b ) ;
227
228 // Read double data
229 int readHeader( double &v, int b ) ;
230
231 // Release DataRecord
232 void releaseRecord() ;
233
234 // show primary information
235 void show() ;
236
237 // convert frequency frame
238 virtual double toLSR( double v, double t, double x, double y ) ;
239
240 // POLNO from RX
241 //uInt polNoFromRX( const char *rx ) ;
242 uInt polNoFromRX( const std::string &rx ) ;
243
244 // Type of file record
245 std::string LOFIL ;
246
247 // Version
248 std::string VER ;
249
250 // Group name
251 std::string GROUP ;
252
253 // Project name
254 std::string PROJ ;
255
256 // Name of observation scheduling file
257 std::string SCHED ;
258
259 // Name of observer
260 std::string OBSVR ;
261
262 // Observation start time with format of "YYYYMMDDHHMMSS" (UTC)
263 std::string LOSTM ;
264
265 // observation end time with format of "YYYYMMDDHHMMSS" (UTC)
266 std::string LOETM ;
267
268 // Number of arrays (beams and IFs)
269 int ARYNM ;
270
271 // Number of scans
272 int NSCAN ;
273
274 // Title of observation
275 std::string TITLE ;
276
277 // Name of target object
278 std::string OBJ ;
279
280 // Equinox (B1950 or J2000)
281 std::string EPOCH ;
282
283 // Right ascension [rad]
284 double RA0 ;
285
286 // Declination [rad]
287 double DEC0 ;
288
289 // Galactic longitude [rad]
290 double GLNG0 ;
291
292 // Galactic latitude [rad]
293 double GLAT0 ;
294
295 // Calibration interval
296 int NCALB ;
297
298 // Scan coordinate (0: RADEC 1: LB 2: AZEL)
299 int SCNCD ;
300
301 // Scan sequence pattern
302 std::string SCMOD ;
303
304 // User-defined recessional velocity [m/s]
305 double URVEL ;
306
307 // Reference frame for recessional velocity (LSR or HEL or GAL)
308 std::string VREF ;
309
310 // Definition of recessional velocity (RAD or OPT)
311 std::string VDEF ;
312
313 // Switching mode (POS or BEAM or FREQ)
314 std::string SWMOD ;
315
316 // Switching frequency [Hz]
317 double FRQSW ;
318
319 // Off-beam angle of beam switching [rad]
320 double DBEAM ;
321
322 // Initial inclination angle of multi-beam array
323 double MLTOF ;
324
325 // Comet: Perihelion distance
326 double CMTQ ;
327
328 // Comet: Eccentricity
329 double CMTE ;
330
331 // Comet: Argument of perihelion
332 double CMTSOM ;
333
334 // Comet: Longitude of the ascending node
335 double CMTNODE ;
336
337 // Comet: Orbital inclination angle
338 double CMTI ;
339
340 // Comet: Time of the perihelion passage
341 std::string CMTTM ;
342
343 // Correction for position of subreflector DX [mm]
344 double SBDX ;
345
346 // Correction for position of subreflector DY [mm]
347 double SBDY ;
348
349 // Correction for position of subreflector DZ1 [mm]
350 double SBDZ1 ;
351
352 // Correction for position of subreflector DZ2 [mm]
353 double SBDZ2 ;
354
355 // Correction for pointing on azimuth [rad]
356 double DAZP ;
357
358 // Correction for pointing on elevation [rad]
359 double DELP ;
360
361 // Number of channel binding
362 int CHBIND ;
363
364 // Number of channel after binding
365 int NUMCH ;
366
367 // Channel range (minimum)
368 int CHMIN ;
369
370 // Channel range (maximum)
371 int CHMAX ;
372
373 // ALC time constant
374 double ALCTM ;
375
376 // Interval to get data from spectrometer
377 double IPTIM ;
378
379 // Position angle of the map
380 double PA ;
381
382 // Length of scan record [bytes]
383 int SCNLEN ;
384
385 // Range of space binding
386 int SBIND ;
387
388 // Quantization bit number (fixed to 12)
389 int IBIT ;
390
391 // Site (antenna) name (45m or ASTE)
392 std::string SITE ;
393
394 // Dummy data
395 std::string CDMY1 ;
396
397 // Type of detector frontend
398 std::vector<std::string> RX ;
399
400 // HPBW [rad]
401 std::vector<double> HPBW ;
402
403 // Aperture efficiencies
404 std::vector<double> EFFA ;
405
406 // Beam efficiencies
407 std::vector<double> EFFB ;
408
409 // Antenna efficiencies
410 std::vector<double> EFFL ;
411
412 // FSS efficiencies
413 std::vector<double> EFSS ;
414
415 // Antenna gain
416 std::vector<double> GAIN ;
417
418 // Type of polarization at feed horn (R or L or H or V)
419 std::vector<std::string> HORN ;
420
421 // Type of polarization (CIRC or LINR)
422 std::vector<std::string> POLTP ;
423
424 // Rotation direction of circular polarization
425 std::vector<double> POLDR ;
426
427 // Polarization angle of linear polarization
428 std::vector<double> POLAN ;
429
430 // Switching frequency of frequcency switching [Hz]
431 std::vector<double> DFRQ ;
432
433 // Type of sideband (LSB or USB or DSB)
434 std::vector<std::string> SIDBD ;
435
436 // Identifier of reference synthesizer
437 std::vector<int> REFN ;
438
439 // Temperature of calibrator
440 std::vector<int> IPINT ;
441
442 // Beam id of the multi-beam detector
443 std::vector<int> MULTN ;
444
445 // Scaling factor of the multi-beam detector
446 std::vector<double> MLTSCF ;
447
448 // Type of LAG window (NONE or HANN or HAMM or BLCK)
449 std::vector<std::string> LAGWIND ;
450
451 // Bandwidth at backend
452 std::vector<double> BEBW ;
453
454 // Spectral resolution at backend
455 std::vector<double> BERES ;
456
457 // Channel width at backend
458 std::vector<double> CHWID ;
459
460 // Array usage (1: used 0: not used)
461 std::vector<int> ARRY ;
462
463 // Frequency calibration: Number of measurement (max 10)
464 std::vector<int> NFCAL ;
465
466 // Frequency calibration: Central frequency [Hz]
467 std::vector<double> F0CAL ;
468
469 // Frequency calibration: Measured central frequency [Hz]
470 std::vector< std::vector<double> > FQCAL ;
471
472 // Frequency calibration: Measured channel number
473 std::vector< std::vector<double> > CHCAL ;
474
475 // Frequency calibration: Measured channel width [Hz]
476 std::vector< std::vector<double> > CWCAL ;
477
478 // DSB scaling factor
479 std::vector<double> DSBFC ;
480
481 // number of scan
482 int scanNum_ ;
483
484 // number of row
485 int rowNum_ ;
486
487 // length of scan (byte)
488 int scanLen_ ;
489
490 // length of spectral data (byte)
491 int dataLen_ ;
492
493 // Data size of the header [bytes]
494 int datasize_ ;
495
496 // maximum channel number
497 int chmax_ ;
498
499 // Current data id
500 int dataid_ ;
501
502 // Data record
503 CountedPtr<NRODataRecord> record_ ;
504
505 // input filename
506 std::string filename_ ;
507
508 // file pointer
509 FILE *fp_ ;
510
511 // OS endian
512 int same_ ;
513
514 // Logger
515 //LogIO os ;
516
517 // reference frequency for each array
518 std::vector<double> refFreq_ ;
519
520 // record to store REFPIX, REFVAL, INCREMENT pair for each array
521 Record frec_ ;
522} ;
523
524// debug message output
525template<class T> inline void nro_debug_output( char *name, int len, std::vector<T> &val )
526{
527 for ( int i = 0 ; i < len ; i++ ) {
528 if ( i == 0 ) {
529 cout << setw(8) << left << name ;
530 }
531 else if ( ( i % 5 ) == 0 ) {
532 cout << endl << " " ;
533 }
534 cout << "\'" << val[i] << "\' " ;
535 }
536 cout << endl ;
537}
538
539template<class T> inline void nro_debug_output( char *name, int len1, int len2, std::vector< std::vector<T> > &val )
540{
541 for ( int i = 0 ; i < len1 ; i++ ) {
542 for ( int j = 0 ; j < len2 ; j++ ) {
543 if ( j == 0 ) {
544 if ( i < 10 )
545 cout << name << "0" << i << " " ;
546 else
547 cout << name << i << " " ;
548 }
549 else if ( ( j % 5 ) == 0 ) {
550 cout << endl << " " ;
551 }
552 cout << "\'" << val[i][j] << "\' " ;
553 }
554 cout << endl ;
555 }
556}
557
558
559#endif /* NRO_HEADER_H */
Note: See TracBrowser for help on using the repository browser.