Changeset 2407 for trunk/external-alma/oldasdm2ASAP
- Timestamp:
- 02/14/12 13:08:36 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/external-alma/oldasdm2ASAP/OldASDMFiller.cc
r2301 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 //uInt subscanno = reader_->getSubscanNo() ;137 131 citer = cycleno.find( scanno ) ; 138 132 if ( citer == cycleno.end() ) … … 157 151 158 152 // subscan number 159 //unsigned int subscanno = reader_->getSubscanNo( idata ) ;160 153 unsigned int subscanno = reader_->getSubscanNo() ; 161 154 162 155 // IFNO 163 //uInt ifno = reader_->getIFNo( idata ) ;164 156 uInt ifno = reader_->getIFNo() ; 165 157 … … 172 164 double sysVel ; 173 165 vector<double> rf ; 174 //reader_->getSourceProperty( idata,175 // srcname,176 // fieldname,177 // srcDirection,178 // srcProperMotion,179 // sysVel,180 // rf ) ;181 166 reader_->getSourceProperty( srcname, 182 167 fieldname, … … 193 178 194 179 // time and interval 195 //casa::Double mjd = (casa::Double)(reader_->getTime( idata )) ;196 //casa::Double interval = (casa::Double)(reader_->getInterval( idata )) ;197 180 casa::Double mjd = (casa::Double)(reader_->getTime()) ; 198 181 casa::Double interval = (casa::Double)(reader_->getInterval()) ; … … 211 194 212 195 // fill FLAGROW 213 //unsigned int flagrow = reader_->getFlagRow( idata ) ;214 196 unsigned int flagrow = reader_->getFlagRow() ; 215 197 setFlagrow( (uInt)flagrow ) ; … … 221 203 float windspeed ; 222 204 float windaz ; 223 //reader_->getWeatherInfo( idata,224 // temperature,225 // pressure,226 // humidity,227 // windspeed,228 // windaz ) ;229 205 reader_->getWeatherInfo( temperature, 230 206 pressure, … … 243 219 double el ; 244 220 vector<double> srate ; 245 //reader_->getPointingInfo( idata,246 // dir,247 // az,248 // el,249 // srate ) ;250 221 reader_->getPointingInfo( dir, 251 222 az, … … 275 246 } 276 247 else { 277 //reader_->getFrequency( idata, refpix, refval, incr, freqref ) ;278 248 reader_->getFrequency( refpix, refval, incr, freqref ) ; 279 249 refval = (double)toLSRK( casa::Double(refval), … … 291 261 292 262 // loop on polarization 293 //vector<unsigned int> dataShape = reader_->getDataShape( idata ) ;294 263 vector<unsigned int> dataShape = reader_->getDataShape() ; 295 264 // ostringstream oss ; … … 304 273 //logsink_->postLocally( LogMessage(oss.str(),LogOrigin(className_,funcName,WHERE)) ) ; 305 274 306 //int numPol = reader_->getNumPol( idata ) ;307 275 unsigned int numPol = dataShape[0] ; 308 276 unsigned int numChan = dataShape[1] ; … … 311 279 312 280 // OPACITY 313 //vector<float> tau = reader_->getOpacity( idata ) ;314 281 vector<float> tau = reader_->getOpacity() ; 315 282 Vector<casa::Float> opacity = toVector( tau, numPol ) ; 316 283 317 284 // SPECTRA, FLAGTRA, TSYS, TCAL 318 //float *sp = reader_->getSpectrum( idata ) ;319 285 float *sp = reader_->getSpectrum() ; 320 286 vector< vector<float> > ts ; 321 287 vector< vector<float> > tc ; 322 //reader_->getTcalAndTsys( idata, tc, ts ) ;323 288 reader_->getTcalAndTsys( tc, ts ) ; 324 289 Matrix<casa::Float> spectra = toMatrix( sp, numPol, numChan ) ;
Note:
See TracChangeset
for help on using the changeset viewer.