source: trunk/src/SDContainer.h@ 120

Last change on this file since 120 was 104, checked in by mar637, 20 years ago

added azimuth,elevation,parangle,refbeam,fieldname,tcal,tcaltime

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.8 KB
Line 
1//#---------------------------------------------------------------------------
2//# SDContainer.h: A container class for single dish integrations
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2004
5//# Malte Marquarding, ATNF
6//#
7//# This program is free software; you can redistribute it and/or modify it
8//# under the terms of the GNU General Public License as published by the Free
9//# Software Foundation; either version 2 of the License, or (at your option)
10//# any later version.
11//#
12//# This program is distributed in the hope that it will be useful, but
13//# WITHOUT ANY WARRANTY; without even the implied warranty of
14//# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
15//# Public License for more details.
16//#
17//# You should have received a copy of the GNU General Public License along
18//# with this program; if not, write to the Free Software Foundation, Inc.,
19//# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
20//#
21//# Correspondence concerning this software should be addressed as follows:
22//# Internet email: Malte.Marquarding@csiro.au
23//# Postal address: Malte Marquarding,
24//# Australia Telescope National Facility,
25//# P.O. Box 76,
26//# Epping, NSW, 2121,
27//# AUSTRALIA
28//#
29//# $Id:
30//#---------------------------------------------------------------------------
31#ifndef _SDCONTAINER_H
32#define _SDCONTAINER_H
33
34#include <vector>
35
36#include <casa/aips.h>
37#include <casa/BasicSL/String.h>
38#include <casa/Arrays/Array.h>
39#include <casa/Arrays/Vector.h>
40
41template<class T> class Matrix;
42
43namespace asap {
44
45
46struct SDHeader {
47 Int nchan;
48 Int npol;
49 Int nif;
50 Int nbeam;
51 String observer;
52 String project;
53 String obstype;
54 String antennaname;
55 Vector<Double> antennaposition;
56 Float equinox;
57 String freqref;
58 Double reffreq;
59 Double bandwidth;
60 Double utc;
61 void print() const ;
62};
63
64class SDFrequencyTable {
65
66public:
67
68 SDFrequencyTable() : nFreq_(0) {;}
69 virtual ~SDFrequencyTable() {;}
70
71 Int length() const { return nFreq_;};// # of stored Frequencies
72
73 Double referencePixel(uInt which) const { return refPix_[which];}
74 Double referenceValue(uInt which) const { return refVal_[which];}
75 Double increment(uInt which) const { return increment_[which];}
76 Float equinox() const { return equinox_; }
77 String refFrame() const { return refFrame_; }
78
79 // returns the index into the table
80 // this creates a new one or returns an existing one
81 Int addFrequency(Int refPix, Double refVal, Double inc);
82 void setEquinox(Float eq) { equinox_ = eq; }
83 void setRefFrame(const String& reff) { refFrame_ = reff; }
84
85private:
86 Int nFreq_;
87 Vector<Double> refPix_;
88 Vector<Double> refVal_;
89 Vector<Double> increment_;
90 Float equinox_;
91 String refFrame_;
92};
93
94
95class SDContainer {
96
97public:
98 SDContainer(uInt nBeam, uInt nIF, uInt nPol, uInt nChan);
99 SDContainer(IPosition shp);
100
101 virtual ~SDContainer();
102
103 Bool resize(IPosition shp);
104
105 Bool setSpectrum(const Matrix<Float>& spec,
106 uInt whichBeam, uInt whichIF);
107 Bool putSpectrum(const Array<Float>& spec);
108
109 Bool setFlags(const Matrix<uChar>& flgs,
110 uInt whichBeam, uInt whichIF);
111 Bool putFlags(const Array<uChar>& spec);
112
113 Bool setTsys(const Vector<Float>& ts,
114 uInt whichBeam, uInt whichIF);
115 Bool putTsys(const Array<Float>& spec);
116
117 Bool setDirection(const Vector<Double>& point, uInt whichBeam);
118 Bool putDirection(const Array<Double>& dir);
119
120 Bool setFrequencyMap(uInt freqslot, uInt whichIF);
121 Bool putFreqMap(const Vector<uInt>& freqs);
122
123 Array<Float> getSpectrum(uInt whichBeam, uInt whichIF) const;
124 Array<uChar> getFlags(uInt whichBeam, uInt whichIF) const;
125 Array<Float> getTsys(uInt whichBeam, uInt whichIF) const;
126 Array<Double> getDirection(uInt whichBeam) const;
127
128 const Array<Float>& getSpectrum() const { return spectrum_; }
129 const Array<uChar>& getFlags() const { return flags_; }
130 const Array<Float>& getTsys() const { return tsys_; }
131 const Array<Double>& getDirection() const { return direction_; }
132
133 const Vector<uInt>& getFreqMap() const { return freqidx_; }
134
135 Double timestamp;
136 //Double bandwidth;
137 String sourcename;
138 String fieldname;
139 Double interval;
140 Int scanid;
141 Vector<Float> tcal;
142 String tcaltime;
143 Float azimuth;
144 Float elevation;
145 Float parangle;
146 Int refbeam;
147
148private:
149 uInt nBeam_,nIF_,nPol_,nChan_;
150
151 // (nBeam,nIF,nPol,nChannel)
152 Array<Float> spectrum_;
153 Array<uChar> flags_;
154 // (nBeam,nIF,nPol,[nChannel]) Tsys is not really a function of
155 // channel, but this makes it easier to work with at the expense of
156 // a little memory
157 Array<Float> tsys_;
158 Array<Float> tcal_;
159
160 //(nIF) indx into "global" frequency table
161 Vector<uInt> freqidx_;
162 //(nBeam,2) maybe use Measures here...
163 Array<Double> direction_;
164
165};
166
167} // namespace
168#endif
Note: See TracBrowser for help on using the repository browser.