Ignore:
Timestamp:
10/04/16 18:20:50 (8 years ago)
Author:
Takeshi Nakazato
Message:

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:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc

    r2407 r3106  
    7474  //logsink_->postLocally( LogMessage("sFreqFrame = "+sFreqFrame,LogOrigin(className_,funcName,WHERE)) ) ;
    7575 
    76   Vector<casa::Double> antpos = table_->getHeader().antennaposition ;
     76  Vector<casacore::Double> antpos = table_->getHeader().antennaposition ;
    7777
    7878  // data selection
    … …  
    172172         
    173173          // fill MOLECULE_ID and add MOLECULES row if necessary
    174           Vector<casa::Double> restFreqs( rf.size() ) ;
     174          Vector<casacore::Double> restFreqs( rf.size() ) ;
    175175          for ( uInt i = 0 ; i < rf.size() ; i++ )
    176             restFreqs[i] = (casa::Double)(rf[i]) ;
     176            restFreqs[i] = (casacore::Double)(rf[i]) ;
    177177          setMolecule( restFreqs ) ;
    178178         
    179179          // time and interval
    180           casa::Double mjd = (casa::Double)(reader_->getTime()) ;
    181           casa::Double interval = (casa::Double)(reader_->getInterval()) ;
     180          casacore::Double mjd = (casacore::Double)(reader_->getTime()) ;
     181          casacore::Double interval = (casacore::Double)(reader_->getInterval()) ;
    182182
    183183          // fill TIME and INTERVAL
    … …  
    185185         
    186186          // fill SRCNAME, SRCTYPE, FIELDNAME, SRCDIRECTION, SRCPROPERMOTION, and SRCVELOCITY
    187           Vector<casa::Double> srcDir( 2 ) ;
    188           srcDir[0] = (casa::Double)(srcDirection[0]) ;
    189           srcDir[1] = (casa::Double)(srcDirection[1]) ;
    190           Vector<casa::Double> srcPM( 2 ) ;
    191           srcPM[0] = (casa::Double)(srcProperMotion[0]) ;
    192           srcPM[1] = (casa::Double)(srcProperMotion[1]) ;
    193           setSource( srcname, srctype, fieldname, srcDir, srcPM, (casa::Double)sysVel ) ;
     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 ) ;
    194194
    195195          // fill FLAGROW
    … …  
    208208                                   windspeed,
    209209                                   windaz ) ;
    210           setWeather2( (casa::Float)temperature,
    211                        (casa::Float)pressure,
    212                        (casa::Float)humidity,
    213                        (casa::Float)windspeed,
    214                        (casa::Float)windaz ) ;
     210          setWeather2( (casacore::Float)temperature,
     211                       (casacore::Float)pressure,
     212                       (casacore::Float)humidity,
     213                       (casacore::Float)windspeed,
     214                       (casacore::Float)windaz ) ;
    215215
    216216          // fill AZIMUTH, ELEVATION, DIRECTION and SCANRATE
    … …  
    223223                                    el,
    224224                                    srate ) ;
    225           Vector<casa::Double> scanRate( 2, 0.0 ) ;
    226           Vector<casa::Double> direction( 2, 0.0 ) ;
     225          Vector<casacore::Double> scanRate( 2, 0.0 ) ;
     226          Vector<casacore::Double> direction( 2, 0.0 ) ;
    227227          if ( srate.size() > 0 ) {
    228             scanRate[0] = (casa::Double)(srate[0]) ;
    229             scanRate[1] = (casa::Double)(srate[1]) ;
     228            scanRate[0] = (casacore::Double)(srate[0]) ;
     229            scanRate[1] = (casacore::Double)(srate[1]) ;
    230230          }
    231231          setScanRate( scanRate ) ;
    232232          if ( dir.size() > 0 ) {
    233             direction[0] = (casa::Double)(dir[0]) ;
    234             direction[1] = (casa::Double)(dir[1]) ;
     233            direction[0] = (casacore::Double)(dir[0]) ;
     234            direction[1] = (casacore::Double)(dir[1]) ;
    235235          }
    236236          else {
    … …  
    238238          }
    239239          //logsink_->postLocally( LogMessage("direction = "+String::toString(direction),LogOrigin(className_,funcName,WHERE)) ) ;
    240           setDirection( direction, (casa::Float)az, (casa::Float)el ) ;
     240          setDirection( direction, (casacore::Float)az, (casacore::Float)el ) ;
    241241
    242242           // REFPIX, REFVAL, INCREMENT
    … …  
    247247          else {
    248248            reader_->getFrequency( refpix, refval, incr, freqref ) ;
    249             refval = (double)toLSRK( casa::Double(refval),
     249            refval = (double)toLSRK( casacore::Double(refval),
    250250                                     String(freqref),
    251251                                     mjd,
    … …  
    258258
    259259          // fill FREQ_ID and add FREQUENCIES row if necessary
    260           setFrequency( (casa::Double)refpix, (casa::Double)refval, (casa::Double)incr ) ;
     260          setFrequency( (casacore::Double)refpix, (casacore::Double)refval, (casacore::Double)incr ) ;
    261261
    262262          // loop on polarization
    … …  
    280280          // OPACITY
    281281          vector<float> tau = reader_->getOpacity() ;
    282           Vector<casa::Float> opacity = toVector( tau, numPol ) ;
     282          Vector<casacore::Float> opacity = toVector( tau, numPol ) ;
    283283
    284284          // SPECTRA, FLAGTRA, TSYS, TCAL
    … …  
    287287          vector< vector<float> > tc ;
    288288          reader_->getTcalAndTsys( tc, ts ) ;
    289           Matrix<casa::Float> spectra = toMatrix( sp, numPol, numChan ) ;
     289          Matrix<casacore::Float> spectra = toMatrix( sp, numPol, numChan ) ;
    290290          Vector<uChar> flagtra( numChan, 0 ) ;
    291           Matrix<casa::Float> tsys = toMatrix( ts, numPol, numChan ) ;
    292           Matrix<casa::Float> tcal = toMatrix( tc, numPol, numChan ) ;
     291          Matrix<casacore::Float> tsys = toMatrix( ts, numPol, numChan ) ;
     292          Matrix<casacore::Float> tcal = toMatrix( tc, numPol, numChan ) ;
    293293//           String caltime = "" ;
    294 //           if ( anyNE( tcal, (casa::Float)1.0 ) )
     294//           if ( anyNE( tcal, (casacore::Float)1.0 ) )
    295295//             caltime = toTcalTime( mjd ) ;
    296296          String caltime = toTcalTime( mjd ) ;
    … …  
    374374    string ref ;
    375375    reader_->getSourceDirection( sdir, ref ) ;
    376     Vector<casa::Double> sourceDir( sdir ) ;
     376    Vector<casacore::Double> sourceDir( sdir ) ;
    377377    hdr.reffreq = toLSRK( hdr.reffreq, hdr.freqref, hdr.utc, hdr.antennaposition, sdir, String(ref) ) ;
    378378    hdr.freqref = "LSRK" ;
    … …  
    410410}
    411411
    412 Matrix<casa::Float> OldASDMFiller::toMatrix( float *sp,
     412Matrix<casacore::Float> OldASDMFiller::toMatrix( float *sp,
    413413                                         unsigned int npol,
    414414                                         unsigned int nchan )
    415415{
    416   Matrix<casa::Float> mSp( npol, nchan ) ;
     416  Matrix<casacore::Float> mSp( npol, nchan ) ;
    417417  if ( npol <= 2 ) {
    418418    // 1 or 2 polarization case
    419419    for ( unsigned int ich = 0 ; ich < nchan ; ich++ ) {
    420420      for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ ) {
    421         mSp(ipol,ich) = (casa::Float)(sp[npol*ich+ipol]) ;
     421        mSp(ipol,ich) = (casacore::Float)(sp[npol*ich+ipol]) ;
    422422      }
    423423    }
    … …  
    426426    // 4 polarization case
    427427    for ( unsigned int ich = 0 ; ich < nchan ; ich++ ) {
    428       mSp(0,ich) = (casa::Float)(sp[4*ich]) ;   // Re(XX)
    429       mSp(1,ich) = (casa::Float)(sp[4*ich+4]) ; // Re(YY)
    430       mSp(2,ich) = (casa::Float)(sp[4*ich+2]) ; // Re(XY)
    431       mSp(3,ich) = (casa::Float)(sp[4*ich+3]) ; // Im(XY)
     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)
    432432    }
    433433  }
    … …  
    435435}
    436436
    437 Matrix<casa::Float> OldASDMFiller::toMatrix( vector< vector<float> > &tsys,
     437Matrix<casacore::Float> OldASDMFiller::toMatrix( vector< vector<float> > &tsys,
    438438                                               unsigned int npol,
    439439                                               unsigned int nchan )
    … …  
    441441  unsigned int numRec = tsys.size() ;
    442442  unsigned int numChan = tsys[0].size() ;
    443   Matrix<casa::Float> ret ;
     443  Matrix<casacore::Float> ret ;
    444444  if ( npol == numRec && nchan == numChan ) {
    445445    ret.resize( npol, nchan ) ;
    446446    for ( unsigned int ip = 0 ; ip < npol ; ip++ )
    447447      for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
    448         ret( ip, ic ) = (casa::Float)(tsys[ip][ic]) ;
     448        ret( ip, ic ) = (casacore::Float)(tsys[ip][ic]) ;
    449449  }
    450450  else if ( npol == numRec && numChan == 1 ) {
    451451    ret.resize( npol, 1 ) ;
    452452    for ( unsigned int ip = 0 ; ip < npol ; ip++ )
    453       ret( ip, 0 ) = (casa::Float)(tsys[0][0]) ;
     453      ret( ip, 0 ) = (casacore::Float)(tsys[0][0]) ;
    454454  }
    455455  else if ( numRec == 1 && nchan == numChan ) {
    … …  
    457457    for ( unsigned int ip = 0 ; ip < npol ; ip++ )
    458458      for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
    459         ret( ip, ic ) = (casa::Float)(tsys[0][ic]) ;
     459        ret( ip, ic ) = (casacore::Float)(tsys[0][ic]) ;
    460460  }
    461461  else if ( numRec == 1 && numChan == 1 ) {
    462462    ret.resize( npol, 1 ) ;
    463463    for ( unsigned int ip = 0 ; ip < npol ; ip++ )
    464       ret( ip, 0 ) = (casa::Float)(tsys[0][0]) ;
     464      ret( ip, 0 ) = (casacore::Float)(tsys[0][0]) ;
    465465  }
    466466  else if ( numRec == 2 && npol == 4 && numChan == nchan ) {
    … …  
    469469    ret.resize( npol, nchan ) ;
    470470    for ( unsigned int ic = 0 ; ic < nchan ; ic++ ) {
    471       casa::Float tsysxy = (casa::Float)(0.5*(tsys[0][ic]+tsys[1][ic])) ;
    472       ret( 0, ic ) = (casa::Float)(tsys[0][ic]) ;
    473       ret( 1, ic ) = (casa::Float)(tsys[1][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]) ;
    474474      ret( 2, ic ) = tsysxy ;
    475475      ret( 3, ic ) = tsysxy ;
    … …  
    480480    //       at the moment Tsys[XY] = 0.5*(Tsys[X]+Tsys[Y])
    481481    ret.resize( npol, 1 ) ;
    482     casa::Float tsysxy = (casa::Float)(0.5*(tsys[0][0]+tsys[1][0])) ;
    483     ret( 0, 0 ) = (casa::Float)(tsys[0][0]) ;
    484     ret( 1, 0 ) = (casa::Float)(tsys[1][0]) ;
     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]) ;
    485485    ret( 2, 0 ) = tsysxy ;
    486486    ret( 3, 0 ) = tsysxy ;
    … …  
    490490    for ( unsigned int ip = 0 ; ip < npol ; ip++ )
    491491      for ( unsigned int ic = 0 ; ic < nchan ; ic++ )
    492         ret( ip, ic ) = (casa::Float)(tsys[0][ic]) ;   
     492        ret( ip, ic ) = (casacore::Float)(tsys[0][ic]) ;   
    493493  }
    494494  return ret ;
    495495}
    496496
    497 Vector<casa::Float> OldASDMFiller::toVector( vector<float> &tau,
     497Vector<casacore::Float> OldASDMFiller::toVector( vector<float> &tau,
    498498                                               unsigned int npol )
    499499{
    500500  String funcName = "toVector" ;
    501501
    502   Vector<casa::Float> ret( npol ) ;
     502  Vector<casacore::Float> ret( npol ) ;
    503503  //logsink_->postLocally( LogMessage("tau0="+String::toString(tau[0]),LogOrigin(className_,funcName,WHERE)) ) ;
    504504  if ( npol == 4 ) {
    505     ret[0] = (casa::Float)tau[0] ;
    506     ret[1] = (casa::Float)tau[1] ;
     505    ret[0] = (casacore::Float)tau[0] ;
     506    ret[1] = (casacore::Float)tau[1] ;
    507507    ret[2] = 0.5 * ( ret[0] + ret[1] ) ;
    508508    ret[3] = ret[2] ;
    … …  
    510510  else if ( npol == tau.size() ) {
    511511    for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ )
    512       ret[ipol] = (casa::Float)tau[ipol] ;
     512      ret[ipol] = (casacore::Float)tau[ipol] ;
    513513  }
    514514  else {
    515515    // I don't know how to handle...
    516516    for ( unsigned int ipol = 0 ; ipol < npol ; ipol++ )
    517       ret[ipol] = (casa::Float)tau[0] ;
     517      ret[ipol] = (casacore::Float)tau[0] ;
    518518  }
    519519  //logsink_->postLocally( LogMessage("tau="+String::toString(ret),LogOrigin(className_,funcName,WHERE)) ) ;
    … …  
    521521}
    522522
    523 String OldASDMFiller::toTcalTime( casa::Double mjd )
     523String OldASDMFiller::toTcalTime( casacore::Double mjd )
    524524{
    525525  return MVTime( mjd ).string( MVTime::YMD ) ;
    526526}
    527527
    528 void OldASDMFiller::toJ2000( Vector<casa::Double> &dir,
     528void OldASDMFiller::toJ2000( Vector<casacore::Double> &dir,
    529529                               double az,
    530530                               double el,
    531                                casa::Double mjd,
    532                                Vector<casa::Double> antpos )
     531                               casacore::Double mjd,
     532                               Vector<casacore::Double> antpos )
    533533{
    534534  String funcName = "toJ2000" ;
    535535
    536   Vector<casa::Double> azel( 2 ) ;
     536  Vector<casacore::Double> azel( 2 ) ;
    537537  azel[0] = az ;
    538538  azel[1] = el ;
    … …  
    553553}
    554554
    555 Vector<casa::Double> OldASDMFiller::toJ2000( Vector<casa::Double> dir,
     555Vector<casacore::Double> OldASDMFiller::toJ2000( Vector<casacore::Double> dir,
    556556                                          String dirref,
    557                                           casa::Double mjd,
    558                                           Vector<casa::Double> antpos )
    559 {
    560   Vector<casa::Double> newd( dir ) ;
     557                                          casacore::Double mjd,
     558                                          Vector<casacore::Double> antpos )
     559{
     560  Vector<casacore::Double> newd( dir ) ;
    561561  if ( dirref != "J2000" ) {
    562562    MEpoch me( Quantity( mjd, "d" ), MEpoch::UTC ) ;
    … …  
    593593}
    594594
    595 casa::Double OldASDMFiller::toLSRK( casa::Double freq,
     595casacore::Double OldASDMFiller::toLSRK( casacore::Double freq,
    596596                                 String freqref,
    597                                  casa::Double utc,
    598                                  Vector<casa::Double> antpos,
    599                                  Vector<casa::Double> dir,
     597                                 casacore::Double utc,
     598                                 Vector<casacore::Double> antpos,
     599                                 Vector<casacore::Double> dir,
    600600                                 String dirref )
    601601{
    … …  
    603603
    604604  //logsink_->postLocally( LogMessage("freqref = "+freqref,LogOrigin(className_,funcName,WHERE)) ) ;
    605   casa::Double newf = freq ;
     605  casacore::Double newf = freq ;
    606606  if ( freqref != "LSRK" ) {
    607     MEpoch me( Quantum<casa::Double>( utc, Unit("d") ), MEpoch::UTC ) ;
    608     Vector< Quantum<casa::Double> > antposQ( 3 ) ;
     607    MEpoch me( Quantum<casacore::Double>( utc, Unit("d") ), MEpoch::UTC ) ;
     608    Vector< Quantum<casacore::Double> > antposQ( 3 ) ;
    609609    for ( int i = 0 ; i < 3 ; i++ )
    610       antposQ[i] = Quantum<casa::Double>( antpos[i], Unit("m") ) ;
     610      antposQ[i] = Quantum<casacore::Double>( antpos[i], Unit("m") ) ;
    611611    MPosition mp( antposQ, MPosition::ITRF ) ;
    612612    MDirection::Types dirtype ;
    … …  
    614614    if ( !b )
    615615      dirtype = MDirection::J2000 ;
    616     MDirection md( Quantum<casa::Double>( dir[0], Unit("rad") ),
    617                    Quantum<casa::Double>( dir[1], Unit("rad") ),
     616    MDirection md( Quantum<casacore::Double>( dir[0], Unit("rad") ),
     617                   Quantum<casacore::Double>( dir[1], Unit("rad") ),
    618618                   dirtype ) ;
    619619    MeasFrame mf( me, mp, md ) ;
    … …  
    624624    MFrequency::Convert tolsr( freqtype,
    625625                               MFrequency::Ref( MFrequency::LSRK, mf ) ) ;
    626     newf = tolsr( Quantum<casa::Double>( freq, Unit("Hz") ) ).get( "Hz" ).getValue() ;
     626    newf = tolsr( Quantum<casacore::Double>( freq, Unit("Hz") ) ).get( "Hz" ).getValue() ;
    627627    //logsink_->postLocally( LogMessage("freq = "+String::toString(freq)+", newf = "+String::toString(newf),LogOrigin(className_,funcName,WHERE)) ) ;
    628628  }
Note: See TracChangeset for help on using the changeset viewer.