source: trunk/src/SDContainer.h@ 73

Last change on this file since 73 was 67, checked in by mmarquar, 20 years ago

Fixed various errors arising when compiled without warning suppression.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.5 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 <aips/aips.h>
37#include <aips/Utilities/String.h>
38#include <aips/Arrays/Array.h>
39#include <aips/Arrays/Vector.h>
40
41template<class T> class Matrix;
42
43namespace atnf_sd {
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 setPointing(const Vector<Double>& point, uInt whichBeam) {;}
118
119 Bool setFrequencyMap(uInt freqslot, uInt whichIF);
120 Bool putFreqMap(const Vector<uInt>& freqs);
121
122 Array<Float> getSpectrum(uInt whichBeam, uInt whichIF) const;
123 Array<uChar> getFlags(uInt whichBeam, uInt whichIF) const;
124 Array<Float> getTsys(uInt whichBeam, uInt whichIF) const;
125
126 const Array<Float>& getSpectrum() const { return spectrum_; }
127 const Array<uChar>& getFlags() const { return flags_; }
128 const Array<Float>& getTsys() const { return tsys_; }
129
130 const Vector<uInt>& getFreqMap() const { return freqidx_; }
131
132 Double timestamp;
133 String sourcename;
134 Double interval;
135 Int scanid;
136
137private:
138 uInt nBeam_,nIF_,nPol_,nChan_;
139
140 // (nBeam,nIF,nPol,nChannel)
141 Array<Float> spectrum_;
142 Array<uChar> flags_;
143 // (nBeam,nIF,nPol,[nChannel]) Tsys is not really a function of
144 // channel, but this makes it easier to work with at the expense of
145 // a little memory
146 Array<Float> tsys_;
147 Array<Float> tcal_;
148
149 //(nBeam) maybe use Measures here...
150 //*** Vector<Vector<Double>> pointing_;
151 //(nIF) indx into "global" frequency table
152 Vector<uInt> freqidx_;
153
154};
155
156} // namespace
157#endif
Note: See TracBrowser for help on using the repository browser.