source: trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc@ 3116

Last change on this file since 3116 was 3106, checked in by Takeshi Nakazato, 8 years ago

New Development: No

JIRA Issue: No

Ready for Test: Yes/No

Interface Changes: Yes/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...


Check-in asap modifications from Jim regarding casacore namespace conversion.
File size: 22.5 KB
Line 
1#include<iostream>
2
3#include <STHeader.h>
4
5#include <measures/Measures/MEpoch.h>
6#include <measures/Measures/MPosition.h>
7#include <measures/Measures/MDirection.h>
8#include <measures/Measures/MCDirection.h>
9#include <measures/Measures/MFrequency.h>
10#include <measures/Measures/MCFrequency.h>
11#include <measures/Measures/MeasFrame.h>
12#include <measures/Measures/MeasConvert.h>
13#include <casa/Arrays/Vector.h>
14#include <casa/Arrays/Matrix.h>
15#include <casa/Quanta/MVTime.h>
16#include <casa/Logging/LogMessage.h>
17
18#include "OldASDMFiller.h"
19
20using namespace std ;
21using namespace casa ;
22using namespace asap ;
23
24OldASDMFiller::OldASDMFiller( CountedPtr<Scantable> stable )
25 : FillerBase( stable ),
26 antennaId_( -1 ),
27 antennaName_( "" ),
28 className_("OldASDMFiller")
29{
30 reader_ = new OldASDMReader() ;
31}
32
33OldASDMFiller::~OldASDMFiller()
34{
35 // nothing to do?
36 logsink_ = 0 ;
37}
38
39void OldASDMFiller::setLogger( CountedPtr<LogSinkInterface> &logsink )
40{
41 logsink_ = logsink ;
42 if ( !(reader_.null()) ) {
43 reader_->setLogger( logsink ) ;
44 }
45}
46
47bool OldASDMFiller::open( const string &filename, const Record &rec )
48{
49 String funcName = "open" ;
50 bool status = reader_->open( filename, rec ) ;
51
52 antennaId_ = reader_->getAntennaId() ;
53 antennaName_ = reader_->getAntennaName() ;
54
55 //logsink_->postLocally( LogMessage("antennaId_ = "+String::toString(antennaId_),LogOrigin(className_,funcName,WHERE)) ) ;
56 //logsink_->postLocally( LogMessage("antennaName_ = "+antennaName_,LogOrigin(className_,funcName,WHERE)) ) ;
57
58 return status ;
59}
60
61void OldASDMFiller::fill()
62{
63 String funcName = "fill" ;
64
65 // header
66 fillHeader() ;
67
68 // set Frame for FREQUENCIES table
69 string sFreqFrame = reader_->getFrame() ;
70 //MFrequency::Types freqFrame = toFrameType( sFreqFrame ) ;
71 MFrequency::Types freqFrame = MFrequency::LSRK ;
72 table_->frequencies().setFrame( freqFrame, false ) ;
73 table_->frequencies().setFrame( freqFrame, true ) ;
74 //logsink_->postLocally( LogMessage("sFreqFrame = "+sFreqFrame,LogOrigin(className_,funcName,WHERE)) ) ;
75
76 Vector<casacore::Double> antpos = table_->getHeader().antennaposition ;
77
78 // data selection
79 reader_->select() ;
80
81 // pick up valid configDescriptionId
82 Vector<uInt> configDescIdList = reader_->getConfigDescriptionIdList() ;
83 uInt numConfigDescId = configDescIdList.size() ;
84
85 //logsink_->postLocally( LogMessage("configDescIdList = "+String::toString(configDescIdList),LogOrigin(className_,funcName,WHERE)) ) ;
86
87 // get field list
88 Vector<uInt> fieldIdList = reader_->getFieldIdList() ;
89 uInt numFieldId = fieldIdList.size() ;
90
91 //logsink_->postLocally( LogMessage("fieldIdList = "+String::toString(fieldIdList),LogOrigin(className_,funcName,WHERE)) ) ;
92
93 // BEAMNO is always 0 since ALMA antenna is single beam
94 uInt beamno = 0 ;
95
96 // REFBEAMNO is -1
97 setReferenceBeam() ;
98
99 // fill FOCUS_ID and add FOCUS row if necessary
100 setFocus() ;
101
102 // CYCLENO
103 map< unsigned int, unsigned int > cycleno ;
104 map< unsigned int, unsigned int >::iterator citer ;
105
106 for ( unsigned int ifield = 0 ; ifield < numFieldId ; ifield++ ) {
107 for ( uInt icon = 0 ; icon < numConfigDescId ; icon++ ) {
108 //logsink_->postLocally( LogMessage("start configDescId "+String::toString(configDescIdList[icon])+" fieldId "+String::toString(fieldIdList[ifield]),LogOrigin(className_,funcName,WHERE)) ) ;
109
110 if ( !(reader_->setMainRow( configDescIdList[icon], fieldIdList[ifield] )) ) {
111 //logsink_->postLocally( LogMessage("skip configDescId "+String::toString(configDescIdList[icon])+" fieldId "+String::toString(fieldIdList[ifield]),LogOrigin(className_,funcName,WHERE)) ) ;
112 continue ;
113 }
114
115 // number of rows
116 uInt nrow = reader_->getNumMainRow() ;
117
118 //logsink_->postLocally( LogMessage("There are "+String::toString(nrow)+" rows in Main table corresponding to configDescId "+String::toString(configDescIdList[icon])+" fieldId "+String::toString(fieldIdList[ifield]),LogOrigin(className_,funcName,WHERE)) ) ;
119
120 for ( uInt irow = 0 ; irow < nrow ; irow++ ) {
121
122 // set main row
123 if ( !(reader_->setMainRow( irow )) ) {
124 // skip row since the row doesn't have valid configDescId
125 //logsink_->postLocally( LogMessage("skip "+String::toString(irow)+" since OldASDMReader::setMainrow() returns False",LogOrigin(className_,funcName,WHERE)) ) ;
126 continue ;
127 }
128
129 // scan and subscan
130 unsigned int scanno = reader_->getScanNoOfCurrentRow() ;
131 citer = cycleno.find( scanno ) ;
132 if ( citer == cycleno.end() )
133 cycleno[scanno] = 0 ;
134
135 // set data
136 if ( !(reader_->setData()) ) {
137 // skip row since reader failed to retrieve data
138 //logsink_->postLocally( LogMessage("skip "+String::toString(irow)+" since OldASDMReader::setData() returns False",LogOrigin(className_,funcName,WHERE)) ) ;
139 continue ;
140 }
141
142 unsigned int numData = reader_->getNumData() ;
143 double refpix = 0.0 ;
144 double refval = 0.0 ;
145 double incr = 0.0 ;
146 string freqref = "" ;
147
148 for ( unsigned int idata = 0 ; idata < numData ; idata++ ) {
149 // prepare to extract binary data
150 reader_->prepareData( idata ) ;
151
152 // subscan number
153 unsigned int subscanno = reader_->getSubscanNo() ;
154
155 // IFNO
156 uInt ifno = reader_->getIFNo() ;
157
158 // source spec
159 int srctype = reader_->getSrcType( scanno, subscanno ) ;
160 string srcname ;
161 string fieldname ;
162 vector<double> srcDirection ;
163 vector<double> srcProperMotion ;
164 double sysVel ;
165 vector<double> rf ;
166 reader_->getSourceProperty( srcname,
167 fieldname,
168 srcDirection,
169 srcProperMotion,
170 sysVel,
171 rf ) ;
172
173 // fill MOLECULE_ID and add MOLECULES row if necessary
174 Vector<casacore::Double> restFreqs( rf.size() ) ;
175 for ( uInt i = 0 ; i < rf.size() ; i++ )
176 restFreqs[i] = (casacore::Double)(rf[i]) ;
177 setMolecule( restFreqs ) ;
178
179 // time and interval
180 casacore::Double mjd = (casacore::Double)(reader_->getTime()) ;
181 casacore::Double interval = (casacore::Double)(reader_->getInterval()) ;
182
183 // fill TIME and INTERVAL
184 setTime( mjd, interval ) ;
185
186 // fill SRCNAME, SRCTYPE, FIELDNAME, SRCDIRECTION, SRCPROPERMOTION, and SRCVELOCITY
187 Vector<casacore::Double> srcDir( 2 ) ;
188 srcDir[0] = (casacore::Double)(srcDirection[0]) ;
189 srcDir[1] = (casacore::Double)(srcDirection[1]) ;
190 Vector<casacore::Double> srcPM( 2 ) ;
191 srcPM[0] = (casacore::Double)(srcProperMotion[0]) ;
192 srcPM[1] = (casacore::Double)(srcProperMotion[1]) ;
193 setSource( srcname, srctype, fieldname, srcDir, srcPM, (casacore::Double)sysVel ) ;
194
195 // fill FLAGROW
196 unsigned int flagrow = reader_->getFlagRow() ;
197 setFlagrow( (uInt)flagrow ) ;
198
199 // fill WEATHER_ID and add WEATHER row if necessary
200 float temperature ;
201 float pressure ;
202 float humidity ;
203 float windspeed ;
204 float windaz ;
205 reader_->getWeatherInfo( temperature,
206 pressure,
207 humidity,
208 windspeed,
209 windaz ) ;
210 setWeather2( (casacore::Float)temperature,
211 (casacore::Float)pressure,
212 (casacore::Float)humidity,
213 (casacore::Float)windspeed,
214 (casacore::Float)windaz ) ;
215
216 // fill AZIMUTH, ELEVATION, DIRECTION and SCANRATE
217 vector<double> dir ;
218 double az ;
219 double el ;
220 vector<double> srate ;
221 reader_->getPointingInfo( dir,
222 az,
223 el,
224 srate ) ;
225 Vector<casacore::Double> scanRate( 2, 0.0 ) ;
226 Vector<casacore::Double> direction( 2, 0.0 ) ;
227 if ( srate.size() > 0 ) {
228 scanRate[0] = (casacore::Double)(srate[0]) ;
229 scanRate[1] = (casacore::Double)(srate[1]) ;
230 }
231 setScanRate( scanRate ) ;
232 if ( dir.size() > 0 ) {
233 direction[0] = (casacore::Double)(dir[0]) ;
234 direction[1] = (casacore::Double)(dir[1]) ;
235 }
236 else {
237 toJ2000( direction, az, el, mjd, antpos ) ;
238 }
239 //logsink_->postLocally( LogMessage("direction = "+String::toString(direction),LogOrigin(className_,funcName,WHERE)) ) ;
240 setDirection( direction, (casacore::Float)az, (casacore::Float)el ) ;
241
242 // REFPIX, REFVAL, INCREMENT
243 String ifkey = getIFKey( ifno ) ;
244 if ( ifrec_.isDefined( ifkey ) ) {
245 getFrequencyRec( ifkey, refpix, refval, incr ) ;
246 }
247 else {
248 reader_->getFrequency( refpix, refval, incr, freqref ) ;
249 refval = (double)toLSRK( casacore::Double(refval),
250 String(freqref),
251 mjd,
252 antpos,
253 //direction,
254 srcDir,
255 "J2000" ) ;
256 setFrequencyRec( ifkey, refpix, refval, incr ) ;
257 }
258
259 // fill FREQ_ID and add FREQUENCIES row if necessary
260 setFrequency( (casacore::Double)refpix, (casacore::Double)refval, (casacore::Double)incr ) ;
261
262 // loop on polarization
263 vector<unsigned int> dataShape = reader_->getDataShape() ;
264// ostringstream oss ;
265// for ( unsigned int i = 0 ; i < dataShape.size() ; i++ ) {
266// if ( i == 0 )
267// oss << "dataShape=[" << dataShape[i] << ", " ;
268// else if ( i == dataShape.size()-1 )
269// oss << dataShape[i] << "]" ;
270// else
271// oss << dataShape[i] << ", " ;
272// }
273 //logsink_->postLocally( LogMessage(oss.str(),LogOrigin(className_,funcName,WHERE)) ) ;
274
275 unsigned int numPol = dataShape[0] ;
276 unsigned int numChan = dataShape[1] ;
277
278 //logsink_->postLocally( LogMessage("numPol = "+String::toString(numPol),LogOrigin(className_,funcName,WHERE)) ) ;
279
280 // OPACITY
281 vector<float> tau = reader_->getOpacity() ;
282 Vector<casacore::Float> opacity = toVector( tau, numPol ) ;
283
284 // SPECTRA, FLAGTRA, TSYS, TCAL
285 float *sp = reader_->getSpectrum() ;
286 vector< vector<float> > ts ;
287 vector< vector<float> > tc ;
288 reader_->getTcalAndTsys( tc, ts ) ;
289 Matrix<casacore::Float> spectra = toMatrix( sp, numPol, numChan ) ;
290 Vector<uChar> flagtra( numChan, 0 ) ;
291 Matrix<casacore::Float> tsys = toMatrix( ts, numPol, numChan ) ;
292 Matrix<casacore::Float> tcal = toMatrix( tc, numPol, numChan ) ;
293// String caltime = "" ;
294// if ( anyNE( tcal, (casacore::Float)1.0 ) )
295// caltime = toTcalTime( mjd ) ;
296 String caltime = toTcalTime( mjd ) ;
297
298 for ( unsigned int ipol = 0 ; ipol < numPol ; ipol++ ) {
299
300 // fill SCANNO, CYCLENO, IFNO, POLNO, and BEAMNO
301 setIndex( (uInt)scanno-1, (uInt)cycleno[scanno], ifno, ipol, beamno ) ;
302
303 // fill SPECTRA, FLAGTRA, TSYS
304 setSpectrum( spectra.row(ipol), flagtra, tsys.row(ipol) ) ;
305
306 // fill TCAL_ID and add TCAL row if necessary
307 setTcal2( caltime, tcal.row(ipol) ) ;
308
309 // fill OPACITY
310 setOpacity( opacity[ipol] ) ;
311
312 // commit row
313 commitRow() ;
314 }
315
316 // increment CYCLENO
317 cycleno[scanno]++ ;
318 }
319 }
320 }
321 }
322
323 return ;
324}
325
326void OldASDMFiller::close()
327{
328 reader_->close() ;
329 reader_ = 0 ;
330
331 return ;
332}
333
334void OldASDMFiller::fillHeader()
335{
336 STHeader hdr ;
337
338 reader_->fillHeader( hdr.nchan,
339 hdr.npol,
340 hdr.nif,
341 hdr.nbeam,
342 hdr.observer,
343 hdr.project,
344 hdr.obstype,
345 hdr.antennaname,
346 hdr.antennaposition,
347 hdr.equinox,
348 hdr.freqref,
349 hdr.reffreq,
350 hdr.bandwidth,
351 hdr.utc,
352 hdr.fluxunit,
353 hdr.epoch,
354 hdr.poltype ) ;
355
356 if ( hdr.freqref != "LSRK" ) {
357// String freqref ;
358// if (hdr.freqref == "TOPOCENT") {
359// freqref = "TOPO";
360// } else if (hdr.freqref == "GEOCENTR") {
361// freqref = "GEO";
362// } else if (hdr.freqref == "BARYCENT") {
363// freqref = "BARY";
364// } else if (hdr.freqref == "GALACTOC") {
365// freqref = "GALACTO";
366// } else if (hdr.freqref == "LOCALGRP") {
367// freqref = "LGROUP";
368// } else if (hdr.freqref == "CMBDIPOL") {
369// freqref = "CMB";
370// } else if (hdr.freqref == "SOURCE") {
371// freqref = "REST";
372// }
373 vector<double> sdir ;
374 string ref ;
375 reader_->getSourceDirection( sdir, ref ) ;
376 Vector<casacore::Double> sourceDir( sdir ) ;
377 hdr.reffreq = toLSRK( hdr.reffreq, hdr.freqref, hdr.utc, hdr.antennaposition, sdir, String(ref) ) ;
378 hdr.freqref = "LSRK" ;
379 }
380
381 setHeader( hdr ) ;
382}
383
384String OldASDMFiller::getIFKey( uInt ifno )
385{
386 return "IFNO"+String::toString( ifno ) ;
387}
388
389void OldASDMFiller::getFrequencyRec( String key,
390 double &refpix,
391 double &refval,
392 double &incr )
393{
394 Record frec = ifrec_.asRecord( key ) ;
395 refpix = frec.asdouble( "REFPIX" ) ;
396 refval = frec.asdouble( "REFVAL" ) ;
397 incr = frec.asdouble( "INCREMENT" ) ;
398}
399
400void OldASDMFiller::setFrequencyRec( String key,
401 double refpix,
402 double refval,
403 double incr )
404{
405 Record frec ;
406 frec.define( "REFPIX", refpix ) ;
407 frec.define( "REFVAL", refval ) ;
408 frec.define( "INCREMENT", incr ) ;
409 ifrec_.defineRecord( key, frec ) ;
410}
411
412Matrix<casacore::Float> OldASDMFiller::toMatrix( float *sp,
413 unsigned int npol,
414 unsigned int nchan )
415{
416 Matrix<casacore::Float> mSp( npol, nchan ) ;
417 if ( npol <= 2 ) {
418 // 1 or 2 polarization case
419 for ( unsigned int ich = 0 ; ich < nchan ; ich++ ) {
420 for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ ) {
421 mSp(ipol,ich) = (casacore::Float)(sp[npol*ich+ipol]) ;
422 }
423 }
424 }
425 else {
426 // 4 polarization case
427 for ( unsigned int ich = 0 ; ich < nchan ; ich++ ) {
428 mSp(0,ich) = (casacore::Float)(sp[4*ich]) ; // Re(XX)
429 mSp(1,ich) = (casacore::Float)(sp[4*ich+4]) ; // Re(YY)
430 mSp(2,ich) = (casacore::Float)(sp[4*ich+2]) ; // Re(XY)
431 mSp(3,ich) = (casacore::Float)(sp[4*ich+3]) ; // Im(XY)
432 }
433 }
434 return mSp ;
435}
436
437Matrix<casacore::Float> OldASDMFiller::toMatrix( vector< vector<float> > &tsys,
438 unsigned int npol,
439 unsigned int nchan )
440{
441 unsigned int numRec = tsys.size() ;
442 unsigned int numChan = tsys[0].size() ;
443 Matrix<casacore::Float> ret ;
444 if ( npol == numRec && nchan == numChan ) {
445 ret.resize( npol, nchan ) ;
446 for ( unsigned int ip = 0 ; ip < npol ; ip++ )
447 for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
448 ret( ip, ic ) = (casacore::Float)(tsys[ip][ic]) ;
449 }
450 else if ( npol == numRec && numChan == 1 ) {
451 ret.resize( npol, 1 ) ;
452 for ( unsigned int ip = 0 ; ip < npol ; ip++ )
453 ret( ip, 0 ) = (casacore::Float)(tsys[0][0]) ;
454 }
455 else if ( numRec == 1 && nchan == numChan ) {
456 ret.resize( npol, nchan ) ;
457 for ( unsigned int ip = 0 ; ip < npol ; ip++ )
458 for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
459 ret( ip, ic ) = (casacore::Float)(tsys[0][ic]) ;
460 }
461 else if ( numRec == 1 && numChan == 1 ) {
462 ret.resize( npol, 1 ) ;
463 for ( unsigned int ip = 0 ; ip < npol ; ip++ )
464 ret( ip, 0 ) = (casacore::Float)(tsys[0][0]) ;
465 }
466 else if ( numRec == 2 && npol == 4 && numChan == nchan ) {
467 // TODO: How to determine Tsys for XY?
468 // at the moment Tsys[XY] = 0.5*(Tsys[X]+Tsys[Y])
469 ret.resize( npol, nchan ) ;
470 for ( unsigned int ic = 0 ; ic < nchan ; ic++ ) {
471 casacore::Float tsysxy = (casacore::Float)(0.5*(tsys[0][ic]+tsys[1][ic])) ;
472 ret( 0, ic ) = (casacore::Float)(tsys[0][ic]) ;
473 ret( 1, ic ) = (casacore::Float)(tsys[1][ic]) ;
474 ret( 2, ic ) = tsysxy ;
475 ret( 3, ic ) = tsysxy ;
476 }
477 }
478 else if ( numRec == 2 && npol == 4 && numChan == 1 ) {
479 // TODO: How to determine Tsys for XY?
480 // at the moment Tsys[XY] = 0.5*(Tsys[X]+Tsys[Y])
481 ret.resize( npol, 1 ) ;
482 casacore::Float tsysxy = (casacore::Float)(0.5*(tsys[0][0]+tsys[1][0])) ;
483 ret( 0, 0 ) = (casacore::Float)(tsys[0][0]) ;
484 ret( 1, 0 ) = (casacore::Float)(tsys[1][0]) ;
485 ret( 2, 0 ) = tsysxy ;
486 ret( 3, 0 ) = tsysxy ;
487 }
488 else {
489 // I don't know how to handle ...
490 for ( unsigned int ip = 0 ; ip < npol ; ip++ )
491 for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
492 ret( ip, ic ) = (casacore::Float)(tsys[0][ic]) ;
493 }
494 return ret ;
495}
496
497Vector<casacore::Float> OldASDMFiller::toVector( vector<float> &tau,
498 unsigned int npol )
499{
500 String funcName = "toVector" ;
501
502 Vector<casacore::Float> ret( npol ) ;
503 //logsink_->postLocally( LogMessage("tau0="+String::toString(tau[0]),LogOrigin(className_,funcName,WHERE)) ) ;
504 if ( npol == 4 ) {
505 ret[0] = (casacore::Float)tau[0] ;
506 ret[1] = (casacore::Float)tau[1] ;
507 ret[2] = 0.5 * ( ret[0] + ret[1] ) ;
508 ret[3] = ret[2] ;
509 }
510 else if ( npol == tau.size() ) {
511 for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ )
512 ret[ipol] = (casacore::Float)tau[ipol] ;
513 }
514 else {
515 // I don't know how to handle...
516 for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ )
517 ret[ipol] = (casacore::Float)tau[0] ;
518 }
519 //logsink_->postLocally( LogMessage("tau="+String::toString(ret),LogOrigin(className_,funcName,WHERE)) ) ;
520 return ret ;
521}
522
523String OldASDMFiller::toTcalTime( casacore::Double mjd )
524{
525 return MVTime( mjd ).string( MVTime::YMD ) ;
526}
527
528void OldASDMFiller::toJ2000( Vector<casacore::Double> &dir,
529 double az,
530 double el,
531 casacore::Double mjd,
532 Vector<casacore::Double> antpos )
533{
534 String funcName = "toJ2000" ;
535
536 Vector<casacore::Double> azel( 2 ) ;
537 azel[0] = az ;
538 azel[1] = el ;
539// MEpoch me( Quantity( mjd, "d" ), MEpoch::UTC ) ;
540// Vector<Quantity> qantpos( 3 ) ;
541// qantpos[0] = Quantity( antpos[0], "m" ) ;
542// qantpos[1] = Quantity( antpos[1], "m" ) ;
543// qantpos[2] = Quantity( antpos[2], "m" ) ;
544// MPosition mp( MVPosition( qantpos ),
545// MPosition::ITRF ) ;
546// // mp.print( os_.output() ) ;
547// MeasFrame mf( me, mp ) ;
548// MDirection::Convert toj2000( MDirection::AZELGEO,
549// MDirection::Ref( MDirection::J2000, mf ) ) ;
550// dir = toj2000( azel ).getAngle( "rad" ).getValue() ;
551 dir = toJ2000( azel, "AZELGEO", mjd, antpos ) ;
552 //logsink_->postLocally( LogMessage("dir = "+String::toString(dir),LogOrigin(className_,funcName,WHERE)) ) ;
553}
554
555Vector<casacore::Double> OldASDMFiller::toJ2000( Vector<casacore::Double> dir,
556 String dirref,
557 casacore::Double mjd,
558 Vector<casacore::Double> antpos )
559{
560 Vector<casacore::Double> newd( dir ) ;
561 if ( dirref != "J2000" ) {
562 MEpoch me( Quantity( mjd, "d" ), MEpoch::UTC ) ;
563 Vector<Quantity> qantpos( 3 ) ;
564 qantpos[0] = Quantity( antpos[0], "m" ) ;
565 qantpos[1] = Quantity( antpos[1], "m" ) ;
566 qantpos[2] = Quantity( antpos[2], "m" ) ;
567 MPosition mp( MVPosition( qantpos ),
568 MPosition::ITRF ) ;
569 // mp.print( os_.output() ) ;
570 MeasFrame mf( me, mp ) ;
571 MDirection::Types dirtype ;
572 Bool b = MDirection::getType( dirtype, dirref ) ;
573 if ( b ) {
574 MDirection::Convert toj2000( dirtype,
575 MDirection::Ref( MDirection::J2000, mf ) ) ;
576 newd = toj2000( dir ).getAngle( "rad" ).getValue() ;
577 }
578 }
579 return newd ;
580}
581
582MFrequency::Types OldASDMFiller::toFrameType( string &s )
583{
584 MFrequency::Types ftype = MFrequency::DEFAULT ;
585 if ( s == "LABREST" )
586 ftype = MFrequency::REST ;
587 else {
588 Bool b = MFrequency::getType( ftype, String(s) ) ;
589 if (!b)
590 ftype = MFrequency::DEFAULT ;
591 }
592 return ftype ;
593}
594
595casacore::Double OldASDMFiller::toLSRK( casacore::Double freq,
596 String freqref,
597 casacore::Double utc,
598 Vector<casacore::Double> antpos,
599 Vector<casacore::Double> dir,
600 String dirref )
601{
602 String funcName = "toLSRK" ;
603
604 //logsink_->postLocally( LogMessage("freqref = "+freqref,LogOrigin(className_,funcName,WHERE)) ) ;
605 casacore::Double newf = freq ;
606 if ( freqref != "LSRK" ) {
607 MEpoch me( Quantum<casacore::Double>( utc, Unit("d") ), MEpoch::UTC ) ;
608 Vector< Quantum<casacore::Double> > antposQ( 3 ) ;
609 for ( int i = 0 ; i < 3 ; i++ )
610 antposQ[i] = Quantum<casacore::Double>( antpos[i], Unit("m") ) ;
611 MPosition mp( antposQ, MPosition::ITRF ) ;
612 MDirection::Types dirtype ;
613 Bool b = MDirection::getType( dirtype, dirref ) ;
614 if ( !b )
615 dirtype = MDirection::J2000 ;
616 MDirection md( Quantum<casacore::Double>( dir[0], Unit("rad") ),
617 Quantum<casacore::Double>( dir[1], Unit("rad") ),
618 dirtype ) ;
619 MeasFrame mf( me, mp, md ) ;
620 MFrequency::Types freqtype ;
621 b = MFrequency::getType( freqtype, freqref ) ;
622 if ( !b )
623 freqtype = MFrequency::TOPO ;
624 MFrequency::Convert tolsr( freqtype,
625 MFrequency::Ref( MFrequency::LSRK, mf ) ) ;
626 newf = tolsr( Quantum<casacore::Double>( freq, Unit("Hz") ) ).get( "Hz" ).getValue() ;
627 //logsink_->postLocally( LogMessage("freq = "+String::toString(freq)+", newf = "+String::toString(newf),LogOrigin(className_,funcName,WHERE)) ) ;
628 }
629 return newf ;
630}
Note: See TracBrowser for help on using the repository browser.