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