Changeset 2301 for trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc

Timestamp:

09/16/11 14:02:21 (13 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-1913

Ready for Test: Yes

Interface Changes: No

What Interface Changed:

Test Programs: test_importasdm_sd

Put in Release Notes: Yes/No?

Module(s): Module Names change impacts.

Description: Describe your changes here...

Refactoring the code.

File:

• trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc (modified) (13 diffs)

trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc

@@ -103 +103 @@
 
   // CYCLENO
-  unsigned int cycleno = 0 ;
+  //unsigned int cycleno = 0 ;
+  map< unsigned int, unsigned int > cycleno ;
+  map< unsigned int, unsigned int >::iterator citer ;
 
   for ( unsigned int ifield = 0 ; ifield < numFieldId ; ifield++ ) {
@@ -125 +127 @@
         if ( !(reader_->setMainRow( irow )) ) {
           // skip row since the row doesn't have valid configDescId
-          //logsink_->postLocally( LogMessage("skip "+String::toString(irow),LogOrigin(className_,funcName,WHERE)) ) ;
+          //logsink_->postLocally( LogMessage("skip "+String::toString(irow)+" since OldASDMReader::setMainrow() returns False",LogOrigin(className_,funcName,WHERE)) ) ;
           continue ;
         }
 
         // scan and subscan
-        unsigned int scanno = reader_->getScanNo() ;
+        //unsigned int scanno = reader_->getScanNo() ;
+        unsigned int scanno = reader_->getScanNoOfCurrentRow() ;
         //uInt subscanno = reader_->getSubscanNo() ;
+        citer = cycleno.find( scanno ) ;
+        if ( citer == cycleno.end() )
+          cycleno[scanno] = 0 ;
 
         // set data
         if ( !(reader_->setData()) ) {
           // skip row since reader failed to retrieve data
-          //logsink_->postLocally( LogMessage("skip "+String::toString(irow),LogOrigin(className_,funcName,WHERE)) ) ;
+          //logsink_->postLocally( LogMessage("skip "+String::toString(irow)+" since OldASDMReader::setData() returns False",LogOrigin(className_,funcName,WHERE)) ) ;
           continue ;
         }
@@ -147 +153 @@
 
         for ( unsigned int idata = 0 ; idata < numData ; idata++ ) {
+          // prepare to extract binary data
+          reader_->prepareData( idata ) ;
 
           // subscan number
-          unsigned int subscanno = reader_->getSubscanNo( idata ) ;
+          //unsigned int subscanno = reader_->getSubscanNo( idata ) ;
+          unsigned int subscanno = reader_->getSubscanNo() ;
 
           // IFNO
-          uInt ifno = reader_->getIFNo( idata ) ;
-
-          // rest frequency
-          vector<double> rf = reader_->getRestFrequency( idata ) ;
+          //uInt ifno = reader_->getIFNo( idata ) ;
+          uInt ifno = reader_->getIFNo() ;
+
+          // source spec
+          int srctype = reader_->getSrcType( scanno, subscanno ) ;
+          string srcname ;
+          string fieldname ;
+          vector<double> srcDirection ;
+          vector<double> srcProperMotion ;
+          double sysVel ;
+          vector<double> rf ;
+          //reader_->getSourceProperty( idata, 
+          //                            srcname, 
+          //                            fieldname,
+          //                            srcDirection,
+          //                            srcProperMotion,
+          //                            sysVel,
+          //                            rf ) ;
+          reader_->getSourceProperty( srcname, 
+                                      fieldname,
+                                      srcDirection,
+                                      srcProperMotion,
+                                      sysVel,
+                                      rf ) ;
           
           // fill MOLECULE_ID and add MOLECULES row if necessary
@@ -164 +193 @@
           
           // time and interval
-          casa::Double mjd = (casa::Double)(reader_->getTime( idata )) ;
-          casa::Double interval = (casa::Double)(reader_->getInterval( idata )) ;
+          //casa::Double mjd = (casa::Double)(reader_->getTime( idata )) ;
+          //casa::Double interval = (casa::Double)(reader_->getInterval( idata )) ;
+          casa::Double mjd = (casa::Double)(reader_->getTime()) ;
+          casa::Double interval = (casa::Double)(reader_->getInterval()) ;
 
           // fill TIME and INTERVAL
           setTime( mjd, interval ) ;
-
-
-          // source spec
-          string srcname = reader_->getSourceName( idata ) ;
-          string fieldname = reader_->getFieldName( idata ) ;
-          int srctype = reader_->getSrcType( scanno, subscanno ) ;
-          vector<double> srcDirection = reader_->getSourceDirection( idata ) ;
-          //vector<double> srcDirection ;
-          //string srcDirRef ;
-          //reader_->getSourceDirection( idata, srcDirection, srcDirRef ) ;
-          vector<double> srcProperMotion = reader_->getSourceProperMotion( idata ) ;
-          double sysVel = reader_->getSysVel( idata ) ;
           
           // fill SRCNAME, SRCTYPE, FIELDNAME, SRCDIRECTION, SRCPROPERMOTION, and SRCVELOCITY
@@ -192 +211 @@
 
           // fill FLAGROW
-          unsigned int flagrow = reader_->getFlagRow( idata ) ;
+          //unsigned int flagrow = reader_->getFlagRow( idata ) ;
+          unsigned int flagrow = reader_->getFlagRow() ;
           setFlagrow( (uInt)flagrow ) ;
 
@@ -201 +221 @@
           float windspeed ;
           float windaz ;
-          reader_->getWeatherInfo( idata,
-                                   temperature, 
+          //reader_->getWeatherInfo( idata,
+          //                         temperature, 
+          //                         pressure,
+          //                         humidity,
+          //                         windspeed,
+          //                         windaz ) ;
+          reader_->getWeatherInfo( temperature, 
                                    pressure,
                                    humidity,
@@ -218 +243 @@
           double el ;
           vector<double> srate ;
-          reader_->getPointingInfo( idata,
-                                    dir,
+          //reader_->getPointingInfo( idata,
+          //                          dir,
+          //                          az,
+          //                          el,
+          //                          srate ) ;
+          reader_->getPointingInfo( dir,
                                     az,
                                     el,
@@ -246 +275 @@
           }
           else {
-            reader_->getFrequency( idata, refpix, refval, incr, freqref ) ;
+            //reader_->getFrequency( idata, refpix, refval, incr, freqref ) ;
+            reader_->getFrequency( refpix, refval, incr, freqref ) ;
             refval = (double)toLSRK( casa::Double(refval), 
                                      String(freqref),
@@ -261 +291 @@
 
           // loop on polarization
-          vector<unsigned int> dataShape = reader_->getDataShape( idata ) ;
+          //vector<unsigned int> dataShape = reader_->getDataShape( idata ) ;
+          vector<unsigned int> dataShape = reader_->getDataShape() ;
 //           ostringstream oss ;
 //           for ( unsigned int i = 0 ; i < dataShape.size() ; i++ ) {
@@ -271 +302 @@
 //               oss << dataShape[i] << ", " ;
 //           }
-//           logsink_->postLocally( LogMessage(oss.str(),LogOrigin(className_,funcName,WHERE)) ) ;
+          //logsink_->postLocally( LogMessage(oss.str(),LogOrigin(className_,funcName,WHERE)) ) ;
                                      
           //int numPol = reader_->getNumPol( idata ) ;
@@ -280 +311 @@
 
           // OPACITY
-          vector<float> tau = reader_->getOpacity( idata ) ;
+          //vector<float> tau = reader_->getOpacity( idata ) ;
+          vector<float> tau = reader_->getOpacity() ;
           Vector<casa::Float> opacity = toVector( tau, numPol ) ;
 
           // SPECTRA, FLAGTRA, TSYS, TCAL
-          float *sp = reader_->getSpectrum( idata ) ;
-          //vector< vector<float> > ts = reader_->getTsys( idata ) ;
-          //vector< vector<float> > tc = reader_->getTcal( idata ) ;
+          //float *sp = reader_->getSpectrum( idata ) ;
+          float *sp = reader_->getSpectrum() ;
           vector< vector<float> > ts ;
           vector< vector<float> > tc ;
-          reader_->getTcalAndTsys( idata, tc, ts ) ;
+          //reader_->getTcalAndTsys( idata, tc, ts ) ;
+          reader_->getTcalAndTsys( tc, ts ) ;
           Matrix<casa::Float> spectra = toMatrix( sp, numPol, numChan ) ;
           Vector<uChar> flagtra( numChan, 0 ) ;
@@ -302 +334 @@
 
             // fill SCANNO, CYCLENO, IFNO, POLNO, and BEAMNO
-            setIndex( (uInt)scanno-1, (uInt)cycleno, ifno, ipol, beamno ) ;
+            setIndex( (uInt)scanno-1, (uInt)cycleno[scanno], ifno, ipol, beamno ) ;
 
             // fill SPECTRA, FLAGTRA, TSYS
@@ -318 +350 @@
 
           // increment CYCLENO
-          cycleno++ ;
+          cycleno[scanno]++ ;
         }
       }