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