source: trunk/src/SDAttr.cc@ 360

Last change on this file since 360 was 353, checked in by kil064, 20 years ago

new class to hold instrument attributes like beam efficiencies, diameters etc
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.2 KB
RevLine 
[353]1//#---------------------------------------------------------------------------
2//# SDAttr.cc: A collection of attributes for different telescopes
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2004
5//# 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
32#include "SDAttr.h"
33#include <casa/aips.h>
34#include <casa/Arrays/Vector.h>
35#include <casa/Arrays/ArrayMath.h>
36#include <casa/Exceptions.h>
37#include <casa/Quanta/QC.h>
38#include <casa/Quanta/Quantum.h>
39
40#include <measures/Measures/MEpoch.h>
41
42#include <scimath/Mathematics/InterpolateArray1D.h>
43
44using namespace casa;
45using namespace asap;
46
47SDAttr::SDAttr ()
48{
49 init();
50}
51
52SDAttr::SDAttr(const SDAttr& other)
53{
54 init(); // state just private 'static' data
55}
56
57SDAttr& SDAttr::operator=(const SDAttr& other)
58{
59 if (this != &other) {
60 ; // state just private 'static' data
61 }
62 return *this;
63}
64
65SDAttr::~SDAttr()
66{;}
67
68
69Float SDAttr::diameter (Instrument inst) const
70{
71 Float D = 1.0;
72 if (inst==ATPKSMB || inst==ATPKSHOH) {
73 D = 64.0;
74 } else if (inst==ATMOPRA) {
75 D = 22.0;
76 } else if (inst==TIDBINBILLA) {
77 D = 70.0;
78 } else if (inst==CEDUNA) {
79 D = 30.0;
80 } else if (inst==HOBART) {
81 D = 26.0;
82 } else {
83 throw(AipsError("Unknown instrument"));
84 }
85//
86 return D;
87}
88
89Vector<Float> SDAttr::beamEfficiency (Instrument inst, const MEpoch& dateObs, const Vector<Float>& freqs) const
90{
91
92// Look at date where appropriate
93
94 Vector<Float> facs(freqs.nelements(),1.0);
95 if (inst==ATPKSMB) {
96 cerr << "No beam efficiency data for this instrument - assuming unity" << endl;
97 } else if (inst==ATPKSHOH) {
98 cerr << "No beam efficiency data for this instrument - assuming unity" << endl;
99 } else if (inst==ATMOPRA) {
100 facs = interp (freqs/1.0e9f, MopEtaBeamX_, MopEtaBeam2004Y_);
101 } else if (inst==TIDBINBILLA) {
102 cerr << "No beam efficiency data for this instrument - assuming unity" << endl;
103 } else if (inst==CEDUNA) {
104 cerr << "No beam efficiency data for this instrument - assuming unity" << endl;
105 } else if (inst==HOBART) {
106 cerr << "No beam efficiency data for this instrument - assuming unity" << endl;
107 } else {
108 }
109//
110 return facs;
111}
112
113Vector<Float> SDAttr::apertureEfficiency (Instrument inst, const MEpoch& dateObs, const Vector<Float>& freqs) const
114{
115
116// Look at date where appropriate
117
118 Vector<Float> facs(freqs.nelements(),1.0);
119 if (inst==ATPKSMB) {
120 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
121 } else if (inst==ATPKSHOH) {
122 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
123 } else if (inst==ATMOPRA) {
124 facs = interp (freqs/1.0e9f, MopEtaApX_, MopEtaAp2004Y_);
125 } else if (inst==TIDBINBILLA) {
126 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
127 } else if (inst==CEDUNA) {
128 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
129 } else if (inst==HOBART) {
130 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
131 } else {
132 cerr << "No aperture efficiency data for this instrument - assuming unity" << endl;
133 }
134 return facs;
135}
136
137Vector<Float> SDAttr::JyPerK (Instrument inst, const MEpoch& dateObs, const Vector<Float>& freqs) const
138{
139
140// FInd what we need
141
142 Vector<Float> etaAp = apertureEfficiency (inst, dateObs, freqs);
143 Float D = diameter(inst);
144
145// Compute it
146
147 Vector<Float> facs(freqs.nelements(),1.0);
148 for (uInt i=0; i<freqs.nelements(); i++) {
149 facs(i) = SDAttr::findJyPerKFac (etaAp(i), D);
150 }
151//
152 return facs;
153}
154
155
156Float SDAttr::findJyPerKFac (Float etaAp, Float D)
157//
158// Converts K -> Jy
159// D in m
160//
161{
162 Double kb = QC::k.getValue(Unit(String("erg/K")));
163 Float gA = C::pi * D * D / 4.0;
164 return (2.0 * 1.0e19 * kb / etaAp / gA);
165}
166
167
168
169// Private
170
171Vector<Float> SDAttr::interp (const Vector<Float>& xOut, const Vector<Float>& xIn,
172 const Vector<Float>& yIn) const
173{
174 Int method = 1; // Linear
175 Vector<Float> yOut;
176 Vector<Bool> mOut;
177//
178 Vector<Bool> mIn(xIn.nelements(),True);
179//
180 InterpolateArray1D<Float,Float>::interpolate(yOut, mOut, xOut,
181 xIn, yIn, mIn,
182 method, True, True);
183//
184 return yOut;
185}
186
187void SDAttr::init ()
188{
189// Beam efficiencies
190
191 MopEtaBeamX_.resize(3);
192 MopEtaBeamX_(0) = 86.0;
193 MopEtaBeamX_(1) = 100.0;
194 MopEtaBeamX_(2) = 115.0;
195//
196 MopEtaBeam2003Y_.resize(3);
197 MopEtaBeam2003Y_(0) = 0.39;
198 MopEtaBeam2003Y_(1) = 0.37;
199 MopEtaBeam2003Y_(2) = 0.37; // replicated from (1)
200//
201 MopEtaBeam2004Y_.resize(3);
202 MopEtaBeam2004Y_(0) = 0.49;
203 MopEtaBeam2004Y_(1) = 0.44;
204 MopEtaBeam2004Y_(2) = 0.42; // replicated from (1)
205
206
207// Aperture efficiency
208
209 MopEtaApX_.resize(2);
210 MopEtaApX_(0) = 86.0;
211 MopEtaApX_(1) = 115.0;
212//
213 MopEtaAp2004Y_.resize(2);
214 MopEtaAp2004Y_(0) = 0.33;
215 MopEtaAp2004Y_(1) = 0.24;
216}
217
Note: See TracBrowser for help on using the repository browser.