Changeset 2407 for trunk/external-alma/asdm2ASAP
- Timestamp:
- 02/14/12 13:08:36 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/external-alma/asdm2ASAP/ASDMFiller.cc
r2355 r2407 76 76 Vector<casa::Double> antpos = table_->getHeader().antennaposition ; 77 77 78 //STHeader hdr = table_->getHeader() ;79 80 78 // data selection 81 79 reader_->select() ; … … 103 101 104 102 // CYCLENO 105 //unsigned int cycleno = 0 ;106 103 map< unsigned int, unsigned int > cycleno ; 107 104 map< unsigned int, unsigned int >::iterator citer ; … … 111 108 //logsink_->postLocally( LogMessage("start configDescId "+String::toString(configDescIdList[icon])+" fieldId "+String::toString(fieldIdList[ifield]),LogOrigin(className_,funcName,WHERE)) ) ; 112 109 113 //Bool status = reader_->setMainRow( configDescIdList[icon], fieldIdList[ifield] ) ;114 110 if ( !(reader_->setMainRow( configDescIdList[icon], fieldIdList[ifield] )) ) { 115 111 //logsink_->postLocally( LogMessage("skip configDescId "+String::toString(configDescIdList[icon])+" fieldId "+String::toString(fieldIdList[ifield]),LogOrigin(className_,funcName,WHERE)) ) ; … … 132 128 133 129 // scan and subscan 134 //unsigned int scanno = reader_->getScanNo() ;135 130 unsigned int scanno = reader_->getScanNoOfCurrentRow() ; 136 131 //logsink_->postLocally( LogMessage("scanno = "+String::toString(scanno),LogOrigin(className_,funcName,WHERE)) ) ; 137 //uInt subscanno = reader_->getSubscanNo() ;138 132 citer = cycleno.find( scanno ) ; 139 133 if ( citer == cycleno.end() ) … … 160 154 161 155 // subscan number 162 //unsigned int subscanno = reader_->getSubscanNo( idata ) ;163 156 unsigned int subscanno = reader_->getSubscanNo() ; 164 157 //logsink_->postLocally( LogMessage("subscanno = "+String::toString(subscanno),LogOrigin(className_,funcName,WHERE)) ) ; 165 158 166 159 // IFNO 167 //uInt ifno = reader_->getIFNo( idata ) ;168 160 uInt ifno = reader_->getIFNo() ; 169 161 //logsink_->postLocally( LogMessage("ifno = "+String::toString(ifno),LogOrigin(className_,funcName,WHERE)) ) ; … … 177 169 double sysVel ; 178 170 vector<double> rf ; 179 //reader_->getSourceProperty( idata,180 // srcname,181 // fieldname,182 // srcDirection,183 // srcProperMotion,184 // sysVel,185 // rf ) ;186 171 reader_->getSourceProperty( srcname, 187 172 fieldname, … … 199 184 200 185 // time and interval 201 //casa::Double mjd = (casa::Double)(reader_->getTime( idata )) ;202 //casa::Double interval = (casa::Double)(reader_->getInterval( idata )) ;203 186 casa::Double mjd = (casa::Double)(reader_->getTime()) ; 204 187 casa::Double interval = (casa::Double)(reader_->getInterval()) ; … … 218 201 219 202 // fill FLAGROW 220 //unsigned int flagrow = reader_->getFlagRow( idata ) ;221 203 unsigned int flagrow = reader_->getFlagRow() ; 222 204 setFlagrow( (uInt)flagrow ) ; … … 228 210 float windspeed ; 229 211 float windaz ; 230 //reader_->getWeatherInfo( idata,231 // temperature,232 // pressure,233 // humidity,234 // windspeed,235 // windaz ) ;236 212 reader_->getWeatherInfo( temperature, 237 213 pressure, … … 250 226 double el ; 251 227 vector<double> srate ; 252 //reader_->getPointingInfo( idata,253 // dir,254 // az,255 // el,256 // srate ) ;257 228 reader_->getPointingInfo( dir, 258 229 az, … … 282 253 } 283 254 else { 284 //reader_->getFrequency( idata, refpix, refval, incr, freqref ) ;285 255 reader_->getFrequency( refpix, refval, incr, freqref ) ; 286 256 refval = (double)toLSRK( casa::Double(refval), … … 298 268 299 269 // loop on polarization 300 //vector<unsigned int> dataShape = reader_->getDataShape( idata ) ;301 270 vector<unsigned int> dataShape = reader_->getDataShape() ; 302 271 // ostringstream oss ; … … 311 280 // logsink_->postLocally( LogMessage(oss.str(),LogOrigin(className_,funcName,WHERE)) ) ; 312 281 313 //int numPol = reader_->getNumPol( idata ) ;314 282 unsigned int numPol = dataShape[0] ; 315 283 unsigned int numChan = dataShape[1] ; … … 318 286 319 287 // OPACITY 320 //vector<float> tau = reader_->getOpacity( idata ) ;321 288 vector<float> tau = reader_->getOpacity() ; 322 289 Vector<casa::Float> opacity = toVector( tau, numPol ) ; 323 290 324 291 // SPECTRA, FLAGTRA, TSYS, TCAL 325 //float *sp = reader_->getSpectrum( idata ) ;326 //logsink_->postLocally( LogMessage("getting spectra...",LogOrigin(className_,funcName,WHERE)) ) ;327 292 float *sp = reader_->getSpectrum() ; 328 293 //logsink_->postLocally( LogMessage("sp[0] = "+String::toString(sp[0]),LogOrigin(className_,funcName,WHERE)) ) ; 329 294 vector< vector<float> > ts ; 330 295 vector< vector<float> > tc ; 331 //reader_->getTcalAndTsys( idata, tc, ts ) ;332 296 reader_->getTcalAndTsys( tc, ts ) ; 333 297 Matrix<casa::Float> spectra = toMatrix( sp, numPol, numChan ) ;
Note:
See TracChangeset
for help on using the changeset viewer.