- Timestamp:
- 02/14/12 13:08:36 (13 years ago)
- Location:
- trunk/external-alma
- Files:
-
- 2 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 ) ; -
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.